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

From 2014.igem.org

(Difference between revisions)
DucoBlanken (Talk | contribs)
(Created page with " <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> <html><head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <!-- This HTML was a...")
Newer edit →

Revision as of 13:30, 11 October 2014


<!DOCTYPE html

 PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">

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