Team:Grenoble/Modeling/Amplification/ODE

From 2012.igem.org

(Difference between revisions)
Line 98: Line 98:
<div class="encadre" id="red">  
<div class="encadre" id="red">  
<center><img src="https://static.igem.org/mediawiki/2012/8/85/Eq13_grenoble.png" alt="" /></center>
<center><img src="https://static.igem.org/mediawiki/2012/8/85/Eq13_grenoble.png" alt="" /></center>
 +
</div>
 +
<br>
 +
<br>
 +
(CRP-cAMP) is the transcription factor of the gene arac. When it appears in the network, it activates the production of the protein Arac. This is modeled by a Hill function. In addition, there is some outflow linked to the promoter pAraBAD, which is the promoter regulating arac, thus there is a basal production of Arac. We take into account this basal production, because we need to know if because of them our system will always be turned on, thus useless. Arac is also naturally degraded by the bacterium. Thus, we get as the equation of evolution of Arac
 +
concentration:
 +
<br>
 +
<br>
 +
<div class="encadre" id="red">
 +
<center><img src="https://static.igem.org/mediawiki/2012/8/84/Eq14_grenoble.png" alt="" /></center>
 +
</div>
 +
<br>
 +
<br>
 +
Then, the protein Arac complexes with arabinose to create Arac active, written Arac*. It is modeled by the following chemical equation:
 +
</br>
 +
</br>
 +
<center><img src="https://static.igem.org/mediawiki/2012/a/a8/Complexation_arac_arabinose_grenoble.png" alt="" /></center>
 +
</br>
 +
</br>
 +
We get the evolution of Arac* r<SUB>Arac*<SUB>:
 +
</br>
 +
</br>
 +
<center><img src="https://static.igem.org/mediawiki/2012/8/88/Eq15_grenoble.png" alt="" /></center>
 +
</br>
 +
</br>
 +
With the qssa, we get:
 +
</br>
 +
</br>
 +
<center><img src="https://static.igem.org/mediawiki/2012/e/e2/Eq16_grenoble.png" alt="" /></center>
 +
</br>
 +
</br>
 +
In addition we have the conservation equation of Arac:
 +
</br>
 +
</br>
 +
<center><img src="https://static.igem.org/mediawiki/2012/d/d1/Eq17_grenoble.png" alt="" /></center>
 +
</br>
 +
</br>
 +
We assume that we have <img src="https://static.igem.org/mediawiki/2012/8/82/Eq18_grenoble.png" alt="" />
 +
<br>
 +
<br>
 +
<div class="encadre" id="red">
 +
<center><img src="https://static.igem.org/mediawiki/2012/3/32/Eq19_grenoble.png" alt="" /></center>
</div>
</div>
<br>
<br>

Revision as of 16:19, 20 September 2012

iGEM Grenoble 2012

iGEM Grenoble iGEM
iGEM 2012
main page
Project

Preliminary

We will use the quasi steady state approximation (QSSA) then. The idea is that there are quick reactions, such as enzymatic ones, complexations, etc… And there are slow reactions such as protein production. We assume that the evolution speed of an element that is created only by quick reaction is null.



Indeed, we have these types of evolution for the biological elements. The ones involved only in quick reactions are most of the time in a steady state, and there jump from one steady state to an other has an infinite speed, which doesn’t interest us.

Goal

In this part, we want to answer to three questions:

      What is the sensitivity of our system?
      What is the time response?
      What steady states will our system always reach?

The system

Here is the schema of the real system, in orange are the reactions which didn’t appear in the simplified system of the overview:


cAMP is the quorum sensing molecule. When we put some cAMP out of the system, it enters into the system. Then, it complexes with CRP to create (CRP-cAMP), which is the transcription factor of the gene arac. When some Arac is created, it will complex with arabinose to create Arac*. Arac* is the active form of Arac. Arac* with (CRP-cAMP) are the transcription factors of the gene cya. Then when some protein of adenylate cyclase is produced, it will catalyze the production of cAMP.

We first have the complexation of cAMP with CRP to get (CRP-CAMP). It is modeled by this biological reaction:



Thus, we get the evolution speed of (CRP-CAMP), r(CRP-cAMP):



Then, we have the conservation equations:



Thus, we get:



Eventually, we get:



where
We have

Eventually, we have to choose between the two solutions



These two solutions are positive, because:



To chose which solution is the good one, we know that:





Eventually, we get:



(CRP-cAMP) is the transcription factor of the gene arac. When it appears in the network, it activates the production of the protein Arac. This is modeled by a Hill function. In addition, there is some outflow linked to the promoter pAraBAD, which is the promoter regulating arac, thus there is a basal production of Arac. We take into account this basal production, because we need to know if because of them our system will always be turned on, thus useless. Arac is also naturally degraded by the bacterium. Thus, we get as the equation of evolution of Arac concentration:



Then, the protein Arac complexes with arabinose to create Arac active, written Arac*. It is modeled by the following chemical equation:



We get the evolution of Arac* rArac*:



With the qssa, we get:



In addition we have the conservation equation of Arac:



We assume that we have



Retrieved from "http://2012.igem.org/Team:Grenoble/Modeling/Amplification/ODE"