Team:Tsinghua-A/Modeling/Sensitivity

From 2012.igem.org

(Difference between revisions)
 
(13 intermediate revisions not shown)
Line 22: Line 22:
<div id="topTitle">
<div id="topTitle">
-
<h2 id="titleText" style="color:rgb(137,202,154);">Tsinghua-A::Modeling</h2>
+
<h2 id="titleText" style="color:rgb(137,202,154);">Tsinghua-A::Modeling::<span style="color: rgb(12,117,34);">Sensitivity</span></h2>
</div>
</div>
<div id="topWrapper">
<div id="topWrapper">
Line 28: Line 28:
</div>
</div>
        <div class="normalTextWrapper">
        <div class="normalTextWrapper">
-
     <h2  style="margin-top:135px; font-size:38px;">Sensitivity analysis</h2>
+
     <h2  style="margin-top:135px; ">Sensitivity analysis</h2>
-
                 <p>To get a better result, which means more genes in state ‘B’ finally, we analyze on some key parameters to find how to improve our system.</p>
+
                 <p>To get a better result, which means more genes in state ‘B’ int the end, we analyze on some key parameters to find how to improve our system.</p>
                 <h2>Analysis of  K3-</h2>
                 <h2>Analysis of  K3-</h2>
                 <p><img src="https://static.igem.org/mediawiki/2012/1/16/THU-AMS1.png"/>is the rate constant that AraC-arabinose complexes separate from araI site. If we change it from the original 1 to 0.01, we can get the following result.</br>
                 <p><img src="https://static.igem.org/mediawiki/2012/1/16/THU-AMS1.png"/>is the rate constant that AraC-arabinose complexes separate from araI site. If we change it from the original 1 to 0.01, we can get the following result.</br>
Line 39: Line 39:
<img src="https://static.igem.org/mediawiki/2012/f/fe/THU-AMS4.png"/>
<img src="https://static.igem.org/mediawiki/2012/f/fe/THU-AMS4.png"/>
<img src="https://static.igem.org/mediawiki/2012/9/9d/THU-AMS5.png"/>
<img src="https://static.igem.org/mediawiki/2012/9/9d/THU-AMS5.png"/>
-
The result shows that when  <img src="https://static.igem.org/mediawiki/2012/1/16/THU-AMS1.png"/> is to small which means that the rate that AraC-arabinose complexes binds to araI site (  <img src="https://static.igem.org/mediawiki/2012/a/a5/THU-AMS0.png"/>) is much bigger than the rate it separate from aral site (<img src="https://static.igem.org/mediawiki/2012/1/16/THU-AMS1.png"/>), there will be fluctuations before Cre degrade to zero.
+
</br>
 +
The result shows that when  <img src="https://static.igem.org/mediawiki/2012/1/16/THU-AMS1.png"/> is very small which means that the rate that AraC-arabinose complexes binds to araI site (  <img src="https://static.igem.org/mediawiki/2012/a/a5/THU-AMS0.png"/>) is much bigger than that it separates from aral site (<img src="https://static.igem.org/mediawiki/2012/1/16/THU-AMS1.png"/>), there will be fluctuations before Cre degrade to zero.
</p>
</p>
<h2>Analysis of the rate of flip</h2>
<h2>Analysis of the rate of flip</h2>
Line 49: Line 50:
</p>
</p>
<img src="https://static.igem.org/mediawiki/2012/0/0c/THU-AMS9.png"/>
<img src="https://static.igem.org/mediawiki/2012/0/0c/THU-AMS9.png"/>
-
<h2>Analysis of  d</h2><h2 style="font-size:20px;">pcre</h2>
+
<h2>Analysis of  d(pcre)</h2>
-
<p> <img src="https://static.igem.org/mediawiki/2012/f/fb/THU-AMS10.png"/> is the degradation constant of Cre protein. If we change<img src="https://static.igem.org/mediawiki/2012/f/fb/THU-AMS10.png"/>  from 0.14/min to 1.4/min, we can get the following result.</p>
+
<p> <img src="https://static.igem.org/mediawiki/2012/f/fb/THU-AMS10.png"/> is the degradation constant of Cre protein. If we change<img src="https://static.igem.org/mediawiki/2012/f/fb/THU-AMS10.png"/>  from 0.14/min to 1.4/min, we can get the following result.</br>
 +
The concentration of Cre :
 +
</p>
 +
<img src="https://static.igem.org/mediawiki/2012/e/e5/THU-AMS11.png"/>
 +
<p>The concentration of Cre-Loxp :</p>
 +
<img src="https://static.igem.org/mediawiki/2012/4/44/THU-AMS12.png"/>
 +
<p>The number of the genes in different states (green line for original state and red line for state after inversion):
 +
</br>
 +
1 trajectory:
 +
</p>
 +
<img src="https://static.igem.org/mediawiki/2012/6/64/THU-AMS13.png"/>
 +
<p>10 trajectories:</p>
 +
<img src="https://static.igem.org/mediawiki/2012/8/8d/THU-AMS14.png"/>
 +
<p>100 trajectories:</p>
 +
<img src="https://static.igem.org/mediawiki/2012/c/c7/THU-AMS15.png"/>
 +
<p>1000 trajectories:</p>
 +
<img src="https://static.igem.org/mediawiki/2012/a/a1/THU-AMS16.png"/>
 +
<p>From the above results , we can see that if the rate of the degradation of Cre protein increases so that the Cre declines to zero faster and the time when Cre exists decreases , the number of times of flip will be smaller and the inversion may even happen only once when the degradation constant of Cre protein is proper. We can also see that the number of genes in state ‘B’ which we want is bigger for the reason that when lots of genes flip to the state ‘B’ from state ’A’ and ready to flip back , Cre and Cre-Loxp have declined to zero , so these genes keep the state ‘B’ finally.
 +
</br>
 +
To understand the affection of the degradation constant of Cre protein more clearly, we then simulate the process when the degradation constant of Cre protein changes from 0.1 to 1.5. And we can get the following result.
 +
</p>
 +
<img src="https://static.igem.org/mediawiki/2012/7/7f/THU-AMS17.png"/>
 +
<p>To see the result more clearly, we magnify the above figure and just see a part of it. </p>
 +
<img src="https://static.igem.org/mediawiki/2012/e/e8/THU-AMS18.png"/>
 +
<p>The result shows that there is more genes in state ‘B’ which we want when the degradation constant of Cre protein is smaller. So we can expect that if we introduce a feed forward to the system which can accelerate the degradation of Cre protein after a certain time we add arabinose, we can also get more genes in state ‘B’ finally. Based on this idea, we designed a system with feed forward.
 +
 
 +
</p>
        </div>
        </div>
<a href="https://2012.igem.org/Team:Tsinghua-A/Modeling" style="margin-left:40px;font-size:20px;">Return</a>
<a href="https://2012.igem.org/Team:Tsinghua-A/Modeling" style="margin-left:40px;font-size:20px;">Return</a>

Latest revision as of 17:24, 26 September 2012

Tsinghua-A::Modeling::Sensitivity

Sensitivity analysis

To get a better result, which means more genes in state ‘B’ int the end, we analyze on some key parameters to find how to improve our system.

Analysis of K3-

is the rate constant that AraC-arabinose complexes separate from araI site. If we change it from the original 1 to 0.01, we can get the following result.
The concentration of Cre :
The concentration of Cre-Loxp :
The number of genes in different state:

The result shows that when is very small which means that the rate that AraC-arabinose complexes binds to araI site ( ) is much bigger than that it separates from aral site (), there will be fluctuations before Cre degrade to zero.

Analysis of the rate of flip

We use to denote the rate that genes in state ‘C’ change to state ‘B’, as well as the rate that genes in state ‘D’ change to state ‘A’. Then we change from 100 to 1, we can get the following result.

The result, especially the bottom of the number of genes in state ‘A’, shows clearly that genes in state ‘A’ first change to state ‘C’ before they become state ‘B’.
The following figure can give us more details about the affection of .

Analysis of d(pcre)

is the degradation constant of Cre protein. If we change from 0.14/min to 1.4/min, we can get the following result.
The concentration of Cre :

The concentration of Cre-Loxp :

The number of the genes in different states (green line for original state and red line for state after inversion):
1 trajectory:

10 trajectories:

100 trajectories:

1000 trajectories:

From the above results , we can see that if the rate of the degradation of Cre protein increases so that the Cre declines to zero faster and the time when Cre exists decreases , the number of times of flip will be smaller and the inversion may even happen only once when the degradation constant of Cre protein is proper. We can also see that the number of genes in state ‘B’ which we want is bigger for the reason that when lots of genes flip to the state ‘B’ from state ’A’ and ready to flip back , Cre and Cre-Loxp have declined to zero , so these genes keep the state ‘B’ finally.
To understand the affection of the degradation constant of Cre protein more clearly, we then simulate the process when the degradation constant of Cre protein changes from 0.1 to 1.5. And we can get the following result.

To see the result more clearly, we magnify the above figure and just see a part of it.

The result shows that there is more genes in state ‘B’ which we want when the degradation constant of Cre protein is smaller. So we can expect that if we introduce a feed forward to the system which can accelerate the degradation of Cre protein after a certain time we add arabinose, we can also get more genes in state ‘B’ finally. Based on this idea, we designed a system with feed forward.

Return