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

From 2014.igem.org




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
{{:Team:TU_Delft-Leiden/Templates/End}}

Retrieved from "http://2014.igem.org/Team:TU_Delft-Leiden/Modeling/CodeRepository/StandardActivation"