Team:TU Delft-Leiden/Modeling/CodeRepository/StandardActivation

Contents

function RR = StandardActivation(beta,x)

R = zeros(size(x));
RR =zeros(size(x));

Parameters

d_R = 4E-6; %protein degradation rate

% beta
% beta(1) = K_d (Microscopic dissociation constant and also threshold value)
% beta(2) = k1 (transcription rate active promoter)
% beta(3) = k2 (transcription rate repressed promoter)

Determine reporter protein (R) equilibrium concentration

    if x >= 0

        R = (1/d_R).*((beta(2).*x + beta(1).*beta(3))./(beta(1) + x)); %Concentration of R in equilibrium
        R0 = (beta(3)/d_R); %Concentration of R in equilibrium without DNT
        RR = R./(R0);

    else
        RR = 0;
    end

end

Published with MATLAB® R2014a