Team:uOttawa/modeling

From 2013.igem.org

Modeling the Type-I Incoherent Feedforward Loop

The Feedforward loop can be broken down into 3 main steps

    1. Toxin will stimulate production of X (rtTA).

    2. X (rtTA) stimulates production of both Z (sfGFP) and Y (LacI).

    3. Sufficient repression from Y terminates Z (sfGFP) production.

You can take a better look at the system by clicking on the button below!





Equations involved in modeling the biological system


The mathematical model was developed using basic rate formulas of activation and repression as shown in the system of equations above. We then expanded them to include the different parameters required for our model.

These additions are as follows

IPTG concentration to calculate the amount of Y (LacI) repressing Z (sfGFP).

ATC concentration to calculate the amount of X (rtTA) available to bind Y promoter and Z promoter.

Transfer rates of proteins between the cytosol and the nuclei to calculate transfer rates by size of each protein.


Acknowledgements

The differential equations model for the system and implementation in MATLAB were done by Paul Reginato from the University of Waterloo iGEM team. This was only possible as a result of the great collaboration between the University of Waterloo iGEM team and the uOttawa iGEM teams learn more. Parameterization and analysis of the model to show pulse and fold change detection were done by Nicolas Tremblay from the University of Ottawa iGEM team with extensive guidance of Daniel Charlebois advisor for the uOttawa iGEM team.

Equations

The four rate equations for the production of mRNA over time of GEV (1.1), X (rtTA) (2.1), Y (LacI) (3.1) and Z (sfGFP) (4.1) in the nucleus are shown below.




The first sets of equations were used to describe the amount of mRNA transcribed. These equations also took into consideration that the mRNA would be transported out of the nuclei. In the model, we assumed that all of the mRNA molecules would be exported to the cytosol. The differential equations representing the variation in mRNA in the cytosol are represented below in equations 1.2, 2.2, 3.2 and 4.2.




Computing the variation and amounts of mRNA in the cytosol over time then allowed us to compute the variation of protein translated. The equations of translation shown below in equations 1.3, 2.3, 3.3, 4.3 represent the variation of protein in the cytosol over time. The translational rates were calculated based on the average translation rate of a ribosome (9.5 base pair per second 12). Here, we assume that the quantity of ribosome is not a limiting factor and that there are no translational regulations involved in the model. The degradation rates are taken from the average degradation rate of protein in yeast14.





Once the amount of protein calculated, we assumed that they would diffuse at an average rate through the nuclear pore to then affect the transcription of the other genes in the system. The calculations leading to the transfer rates are described in the parameter section under transf. Equations 1.4, 2.4 and 3.4 represent the variation of proteins present in the nucleus over time.





To have better control on the feed-forward loop system, LacI was regulated with IPTG. Changing the amount of ATC regulated the amount of active rtTA. And the Amount of GEV available to activate production of rtTA is controlled with ß-estradiol. The following equations used to compute the amount of each compound are shown below. These numbers are imputed in the mRNA equations and vary accordingly with the amount of each transcription factors expressed in the ODE system.



Parameters

Some parameters were found in the literature. Others were derived from other known constants as explained below.

Transcription rate

The transcription rate was calculated using the known rate of 2kb per minute6 which translates to 33.3 nucleotides per second. This rate was then divided by the amount of nucleotides in one mRNA strand. The values are in the parameter table.

Translation rate

The translation rate was calculated using the known rate of 9.5 amino acid per second12. This rate was then divided by the amount of amino acids per protein. The values are in the parameter table.

Transfer rate of protein from the cytosol to the nucleus

The rate of transfer from the cytosol to the nucleus was extrapolated from the formula rate = 0.01*(5000-size of protein in Dalton)-50. This is the multiplied by 119 or the average number of nuclear pores in yeast7. The first equation was determined by assuming that transfer rate of protein 5000 Da takes 0 seconds (its size being insignificant ) and a 17000 Da protein takes 120 seconds 11.



 

GEV

rtTA

LacI

sfGFP

IPTG

ß-estradiol

ATC

Translation rates (protein/min)

1.25

0.901

0.959

2.395

 

 

 

n (Hill coefficient)

1

2

2

-

2

1

1

K

10^-6 13

10^-6 13

10^-10 1

-

6.0*10^-6 3

10^-94

1.56*10^-913

Transcription rate (mRNA/min)

1.465

1.055

1.122

2.801

 

 

 

Transf

0.00249

0.00168

0.00181

 

 

 

 

Leakage

1.5 molecule/hour for one promoter

References

Retrieved from "http://2013.igem.org/Team:uOttawa/modeling"