Team:Tsinghua-A/Modeling/P2A

From 2011.igem.org

(Difference between revisions)
 
(10 intermediate revisions not shown)
Line 100: Line 100:
h2{
h2{
line-height:2;
line-height:2;
-
font-weight:bold;
 
font-family:"Times New Roman",Georgia,Serif;
font-family:"Times New Roman",Georgia,Serif;
-
color:#0431B4;
+
color:#555555;
-
font-size:175%;
+
font-weight:bold;
 +
font-size:130%;
}
}
Line 130: Line 130:
#menu a {
#menu a {
display: block;
display: block;
-
width: 155px;
+
width: 133px;
height: 33px;
height: 33px;
padding-top: 16px;
padding-top: 16px;
Line 168: Line 168:
   <div id="header">
   <div id="header">
 +
<div style="position:fixed;right:0px;bottom: 0px;display:inline"><a href="#back-top">
 +
<img src="https://static.igem.org/mediawiki/2011/5/52/Back-top.png" width="60px" height="60px"></a></div>
 +
<div style="position:fixed;right:80px;bottom: 0px;display:inline"><a href="https://2011.igem.org/Team:Tsinghua-A/Homepage">
 +
<img src="https://static.igem.org/mediawiki/2011/8/82/Back-home.png" width="60px" height="60px"></a></div>
-
<p align="left"><img src="https://static.igem.org/mediawiki/2011/f/f9/THU_banner_mod.jpg" alt="" width="960"/></p>
+
<p align="CENTER" id="back-top" style="text-indent:0em"><img src="https://static.igem.org/mediawiki/2011/f/f9/THU_banner_mod.jpg" alt="" width="960"/></p>
</div>
</div>
Line 175: Line 179:
<div id="menu">
<div id="menu">
<ul>
<ul>
-
<li><a href="https://2011.igem.org/Team:Tsinghua-A/Team">Team</a></li>
+
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Homepage">Home</a></li>
-
<li><a href="https://2011.igem.org/Team:Tsinghua-A/Project">Project</a></li>
+
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Team">Team</a></li>
-
<li><a href="https://2011.igem.org/Team:Tsinghua-A/safety">Safety</a></li>
+
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Project">Project</a></li>
<li class="current_page_item"><a href="https://2011.igem.org/Team:Tsinghua-A/Modeling">Modeling</a></li>
<li class="current_page_item"><a href="https://2011.igem.org/Team:Tsinghua-A/Modeling">Modeling</a></li>
-
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Parts">Parts</a></li>
+
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Wetlab">Wetlab</a></li>
 +
<li><a href="https://2011.igem.org/Team:Tsinghua-A/Safety">Safety</a></li>
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Notebook">Notebook</a></li>
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Notebook">Notebook</a></li>
</ul>
</ul>
Line 191: Line 196:
</HEAD>
</HEAD>
<BODY LANG="zh-CN">
<BODY LANG="zh-CN">
 +
<div style="position:relative;left:60px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling"><img src="https://static.igem.org/mediawiki/2011/2/2f/ThuA_A1.png" width="110px" height="105px"></A></div>
<div style="position:relative;left:60px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling"><img src="https://static.igem.org/mediawiki/2011/2/2f/ThuA_A1.png" width="110px" height="105px"></A></div>
Line 200: Line 206:
<div style="position:relative;left:140px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P4"><img src="https://static.igem.org/mediawiki/2011/6/6a/ThuA_E1.png" width="110px" height="120px"></A></div>
<div style="position:relative;left:140px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P4"><img src="https://static.igem.org/mediawiki/2011/6/6a/ThuA_E1.png" width="110px" height="120px"></A></div>
-
 
<div class="slider">
<div class="slider">
-
<p><IMG SRC="https://static.igem.org/mediawiki/2011/1/13/ModelBack.jpg" ALIGN=bottom WIDTH=710 HEIGHT=175 BORDER=0 ISMAP></div>
+
<p><IMG SRC="https://static.igem.org/mediawiki/2011/1/13/ModelBack.jpg" ALIGN=bottom WIDTH=710 HEIGHT=250 BORDER=0 ISMAP></div>
-
 
+
<br>
<table id="toc" class="toc">
<table id="toc" class="toc">
Line 236: Line 241:
<P><B>●Protein which is translated from mRNA is proportional to corresponding mRNA in a previous time.</B></P>
<P><B>●Protein which is translated from mRNA is proportional to corresponding mRNA in a previous time.</B></P>
<p>Consider gene lasR and luxR which are all constant genes, since AHL has a relatively low concentration compared to other substances, we assume that protein concentrations of PlasR and PluxR are irrelevant to AHLs. Therefore, concentrations of these proteins will reach to constant quantities. Moreover, certain mRNAs will also be constant. That is:</p>
<p>Consider gene lasR and luxR which are all constant genes, since AHL has a relatively low concentration compared to other substances, we assume that protein concentrations of PlasR and PluxR are irrelevant to AHLs. Therefore, concentrations of these proteins will reach to constant quantities. Moreover, certain mRNAs will also be constant. That is:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/5/58/Part2-1.png" width = "311px" height="63px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/5/58/Part2-1.png" width = "311px" height="63px" />
-
<p>We have:</p>
+
</p><p>We have:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/6/67/Part2-2.png" width = "314px" height="166px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/6/67/Part2-2.png" width = "314px" height="166px" /></p>
<p>Then focus on compound LA1 and LA2. The binding process is far quicker than other reactions such as translation reactions, thus we assume that the compounds are Quasi-equilibrium:</p>
<p>Then focus on compound LA1 and LA2. The binding process is far quicker than other reactions such as translation reactions, thus we assume that the compounds are Quasi-equilibrium:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e4/Part2-3.png" width = "159px" height="60px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e4/Part2-3.png" width = "159px" height="60px" />
-
<p>We have:</p>
+
</p><p>We have:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/2/2c/Part2-4.png" width = "304px" height="121px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/2/2c/Part2-4.png" width = "304px" height="121px" /></p>
<p>Take the expressions to equation (6) (15), the functions are transformed to: </p>
<p>Take the expressions to equation (6) (15), the functions are transformed to: </p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/3/34/Part2-5.png" width = "384px" height="63px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/3/34/Part2-5.png" width = "384px" height="63px" />
-
<p>The transformed functions have nothing to do with feedback factors compared to original ones. Knowing that such feedback is essential for our system, we add additional feedback factors k<sub>a</sub>,k<sub>b</sub> manipulatively.</p>
+
</p><p>The transformed functions have nothing to do with feedback factors compared to original ones. Knowing that such feedback is essential for our system, we add additional feedback factors k<sub>a</sub>,k<sub>b</sub> manipulatively.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/4/40/Part2-6.png" width = "313px" height="120px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/4/40/Part2-6.png" width = "313px" height="120px" /></p>
<p>Then we want to come up a clear expression of A2<sub>c1</sub>,A1<sub>c2</sub>.They are generated by protein PluxI and PlasI, whose translation rate is much slower than the transcription rate. So we assume that certain mRNA is Quasi-equilibrium:</p>
<p>Then we want to come up a clear expression of A2<sub>c1</sub>,A1<sub>c2</sub>.They are generated by protein PluxI and PlasI, whose translation rate is much slower than the transcription rate. So we assume that certain mRNA is Quasi-equilibrium:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/1/11/Part2-7.png" width = "179px" height="58px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/1/11/Part2-7.png" width = "179px" height="58px" />
-
<p>We have:</p>
+
</p><p>We have:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/3/38/Part2-8.png" width = "521px" height="196px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/3/38/Part2-8.png" width = "521px" height="196px" /></p>
<p>By ignoring the basal expression:</p>
<p>By ignoring the basal expression:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/b/bc/Part2-9.png" width = "320px" height="195px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/b/bc/Part2-9.png" width = "320px" height="195px" /></p>
<p>Knowing that protein PtetR is controlled by LA2 through Hill Function, we assume that concentration of mRNA is directly related to LA2:</p>
<p>Knowing that protein PtetR is controlled by LA2 through Hill Function, we assume that concentration of mRNA is directly related to LA2:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e0/Part2-10.png" width = "302px" height="81px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e0/Part2-10.png" width = "302px" height="81px" /></p>
<p>Parameter KM4 is relevant to both KM1 and KM3, however accurate function depicting the relationship is unknown. We estimate the quantity of KM4 by multiplying KM1 and KM3. Noting LA1 and LA2 have already been expressed, we have:</p>
<p>Parameter KM4 is relevant to both KM1 and KM3, however accurate function depicting the relationship is unknown. We estimate the quantity of KM4 by multiplying KM1 and KM3. Noting LA1 and LA2 have already been expressed, we have:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/6/65/Part2-11.png" width = "351px" height="226px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/6/65/Part2-11.png" width = "351px" height="226px" /></p>
<p>Then we have trouble in expressing protein PluxI and PlasI. Since concentration of protein is the integrals of its mRNA, we assume that it is proportional to concentration of mRNA in a previous time. Thus we have a DDE function:</p>
<p>Then we have trouble in expressing protein PluxI and PlasI. Since concentration of protein is the integrals of its mRNA, we assume that it is proportional to concentration of mRNA in a previous time. Thus we have a DDE function:</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/9/9d/Part2-12.png" width = "342px" height="166px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/9/9d/Part2-12.png" width = "342px" height="166px" /></p>
<p>Equations concerning environmental AHLs remain unchanged.</p>
<p>Equations concerning environmental AHLs remain unchanged.</p>
<p>The original ODEs can be transformed to a much simple DDEs by above deductions.</p>
<p>The original ODEs can be transformed to a much simple DDEs by above deductions.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/b/b3/Part2-13.png" width = "722px" height="393px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/b/b3/Part2-13.png" width = "722px" height="393px" /></p>
<br><br>
<br><br>
<h1 id="Parameters">Parameters</h1><hr width="100%" size=2 color=gray>
<h1 id="Parameters">Parameters</h1><hr width="100%" size=2 color=gray>
<p>Some of the parameters are derived from original ones, some of them are created to describe the new equations, and others are set for further testing.</p>
<p>Some of the parameters are derived from original ones, some of them are created to describe the new equations, and others are set for further testing.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e6/Part2-14.png" width = "743px" height="689px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e6/Part2-14.png" width = "743px" height="689px" /></p>
<p align="center" class="cite">Table 2 Parameters of DDEs</p>
<p align="center" class="cite">Table 2 Parameters of DDEs</p>
Line 273: Line 278:
<h2 id="General results">General results</h2>
<h2 id="General results">General results</h2>
<p>We coded the system in MATLAB. The result shows that every signal AHL is oscillating.</p>
<p>We coded the system in MATLAB. The result shows that every signal AHL is oscillating.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/d/db/Part2-15.png" width = "900px" height="675px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/d/db/Part2-15.png" width = "800px" height="600px" /></p>
-
<p class="cite">Figure 3 Signal AHL Oscillations</p>
+
<p class="cite" align="CENTER" style="text-indent:0em">Figure 3 Signal AHL Oscillations</p>
<h2 id="Stability analysis">Stability analysis</h2>
<h2 id="Stability analysis">Stability analysis</h2>
Line 281: Line 286:
<p> k<sub>a</sub>=0,k<sub>b</sub>=0,μ=10min<sup>-1</sup>,γ=2.3min<sup>-1</sup>,β<sub>1</sub>=30nm/min,k<sub>m1</sub>=1,n<sub>1</sub>=n<sub>2</sub>=2,p=0.2,arab=20,t<sub>1</sub>=t<sub>2</sub>=0,τ<sub>1</sub>=6,τ<sub>2</sub>=10. </p>
<p> k<sub>a</sub>=0,k<sub>b</sub>=0,μ=10min<sup>-1</sup>,γ=2.3min<sup>-1</sup>,β<sub>1</sub>=30nm/min,k<sub>m1</sub>=1,n<sub>1</sub>=n<sub>2</sub>=2,p=0.2,arab=20,t<sub>1</sub>=t<sub>2</sub>=0,τ<sub>1</sub>=6,τ<sub>2</sub>=10. </p>
<p>We changed the binary parameter (β<sub>2</sub>,k<sub>m2</sub>) and got the result below.</p>
<p>We changed the binary parameter (β<sub>2</sub>,k<sub>m2</sub>) and got the result below.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/d/d2/Part2-16.png" width = "900px" height="450px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/d/d2/Part2-16.png" width = "900px" height="450px" /></p>
<p align="center" class="cite">Figure 4 C12 and C6 when (β<sub>2</sub>,k<sub>m2</sub>)=(30,0.1)</p>
<p align="center" class="cite">Figure 4 C12 and C6 when (β<sub>2</sub>,k<sub>m2</sub>)=(30,0.1)</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e1/Part2-17.png" width = "900px" height="450px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e1/Part2-17.png" width = "900px" height="450px" /></p>
<p align="center" class="cite">Figure 5 C12 and C6 when (β<sub>2</sub>,k<sub>m2</sub>)=(120,0.1)</p>
<p align="center" class="cite">Figure 5 C12 and C6 when (β<sub>2</sub>,k<sub>m2</sub>)=(120,0.1)</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/0/0d/Part2-18.png" width = "900px" height="450px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/0/0d/Part2-18.png" width = "900px" height="450px" /></p>
<p align="center" class="cite">Figure 6 C12 and C6 when (β<sub>2</sub>,k<sub>m2</sub>)=(60,0.3)</p>
<p align="center" class="cite">Figure 6 C12 and C6 when (β<sub>2</sub>,k<sub>m2</sub>)=(60,0.3)</p>
<p>After simulating the system at different parameters, we recorded several critical points for oscillation and made a table as follows.</p>
<p>After simulating the system at different parameters, we recorded several critical points for oscillation and made a table as follows.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/c/c9/Part2-19.png" width = "712px" height="101px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/c/c9/Part2-19.png" width = "712px" height="101px" /></p>
<p align="center" class="cite">Table 3 Critical points (β<sub>2</sub>,k<sub>m2</sub>) for oscillation</p>
<p align="center" class="cite">Table 3 Critical points (β<sub>2</sub>,k<sub>m2</sub>) for oscillation</p>
<p>Depicting those critical points on an axis, we immediately got the bifurcation line of parameters (β<sub>2</sub>,k<sub>m2</sub>), which indicates the parameters’ value range when our system can oscillate stably, being marked in ‘bistable’.</p>
<p>Depicting those critical points on an axis, we immediately got the bifurcation line of parameters (β<sub>2</sub>,k<sub>m2</sub>), which indicates the parameters’ value range when our system can oscillate stably, being marked in ‘bistable’.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/9/9b/Part2-20.png" width = "900px" height="675px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/9/9b/Part2-20.png" width = "800px" height="600px" /></p>
<p align="center" class="cite">Figure 7 Bifurcation Analysis on (β<sub>2</sub>,k<sub>m2</sub>)</p>
<p align="center" class="cite">Figure 7 Bifurcation Analysis on (β<sub>2</sub>,k<sub>m2</sub>)</p>
<h2 id="Proportion of cell volume">Proportion of cell volume</h2>
<h2 id="Proportion of cell volume">Proportion of cell volume</h2>
<p>  In actual vivo experiment, volumes of two separate cells might not be the same. We simulate the system by changing the proportion of two types of cells. What we know from the simulation is that cell proportions only affect the amplitude of signal molecules’ oscillation, but no influence to the period or stability of oscillation. The following graphics are drawn under cell proportion 1 and 0.7.</p>
<p>  In actual vivo experiment, volumes of two separate cells might not be the same. We simulate the system by changing the proportion of two types of cells. What we know from the simulation is that cell proportions only affect the amplitude of signal molecules’ oscillation, but no influence to the period or stability of oscillation. The following graphics are drawn under cell proportion 1 and 0.7.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/8/8a/Part2-21.png" width = "900px" height="335px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/8/8a/Part2-21.png" width = "900px" height="335px" /></p>
<p align="center" class="cite">Figure 8 oscillations under different cell proportions</p>
<p align="center" class="cite">Figure 8 oscillations under different cell proportions</p>
<p>We can find that under different cell proportions, our system can oscillate at different amplitude, and varying in cell proportions may also lead to the change of oscillation period.</p>
<p>We can find that under different cell proportions, our system can oscillate at different amplitude, and varying in cell proportions may also lead to the change of oscillation period.</p>
Line 303: Line 308:
<p>As figure 2 indicates, we expected to adjust oscillation period by adding signal molecule aTc into our system. A small molecule as aTc is, it can easily bind to protein TetR and quickly depletes TetR in cell 2, which results in reduction of TetR net production rate, and indirectly, the protein’s inhibition on promoter5 is crippled. Thus, it would take a longer time for our system to reach each threshold, which is equivalent to prolonging the time delay τ<sub>2</sub> in our simplified model. So we can deduce that changing the amount of aTc added into the system in precise model is equivalent to varying the time delay τ<sub>2</sub> in simplified model.</p>
<p>As figure 2 indicates, we expected to adjust oscillation period by adding signal molecule aTc into our system. A small molecule as aTc is, it can easily bind to protein TetR and quickly depletes TetR in cell 2, which results in reduction of TetR net production rate, and indirectly, the protein’s inhibition on promoter5 is crippled. Thus, it would take a longer time for our system to reach each threshold, which is equivalent to prolonging the time delay τ<sub>2</sub> in our simplified model. So we can deduce that changing the amount of aTc added into the system in precise model is equivalent to varying the time delay τ<sub>2</sub> in simplified model.</p>
<p>Simulation results under distinct τ<sub>2</sub> are presented as follows.</p>
<p>Simulation results under distinct τ<sub>2</sub> are presented as follows.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/6/66/Part2-22.png" width = "900px" height="335px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/6/66/Part2-22.png" width = "900px" height="335px" /></p>
<p align="center" class="cite">Figure 9 Oscillation cycle’s regulation</p>
<p align="center" class="cite">Figure 9 Oscillation cycle’s regulation</p>
<p>The result is exactly what we expected, which clearly demonstrates our system can truly be controlled by adding in external signal molecules.</p>
<p>The result is exactly what we expected, which clearly demonstrates our system can truly be controlled by adding in external signal molecules.</p>
Line 309: Line 314:
<h2 id="Phase adjustment">Phase adjustment</h2>
<h2 id="Phase adjustment">Phase adjustment</h2>
<p>In cell 2’s gene circuit, we designed a promoter induced by arabinose, marked by promoter 6 (see in figure 2). Promoter 6 is in a suppressed state until being induced by adding arabinose, and after the inhibition is relieved, signal molecule 3O12HSL will be generated extra. We can also analyze the differential equations describing dA1C2/dt, when adding arabinose during period from t<sub>1</sub> to t<sub>2</sub>, dA1C2/dt contains an extra item arab?(t>t<sub>1</sub>)*(t<t<sub>2</sub>). The parameter arab can reflect the rate of adding arabinose nonlinearly. Here we set arab=20,t<sub>1</sub>=80,t<sub>2</sub>=120 and simulation result is presented as follow.</p>
<p>In cell 2’s gene circuit, we designed a promoter induced by arabinose, marked by promoter 6 (see in figure 2). Promoter 6 is in a suppressed state until being induced by adding arabinose, and after the inhibition is relieved, signal molecule 3O12HSL will be generated extra. We can also analyze the differential equations describing dA1C2/dt, when adding arabinose during period from t<sub>1</sub> to t<sub>2</sub>, dA1C2/dt contains an extra item arab?(t>t<sub>1</sub>)*(t<t<sub>2</sub>). The parameter arab can reflect the rate of adding arabinose nonlinearly. Here we set arab=20,t<sub>1</sub>=80,t<sub>2</sub>=120 and simulation result is presented as follow.</p>
-
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e5/Part2-23.png" width = "900px" height="675px" />
+
<p align="CENTER" style="text-indent:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/e/e5/Part2-23.png" width = "800px" height="600px" /></p>
<p align="center" class="cite">Figure 10 Oscillation phase’s regulation</p>
<p align="center" class="cite">Figure 10 Oscillation phase’s regulation</p>

Latest revision as of 03:52, 28 October 2011

Republic by Free CSS Templates



Contents


Simplification


Although ODEs provide a thorough, precise description of the whole system, they contain too many equations and parameters which would act as a barrier for simulation and further analysis. A simplification of complicated ODEs is necessary. We simplify every single ODE according to certain appropriate assumptions. Finally, we came up with a set of DDE equations.

Assumptions we have made to deduct the original equations involve:

●Relatively faster reactions such as transcription reactions and binding reactions will reach to Quasi-equilibrium.

●Basal expression of protein is so meager that they can be ignored in modeling.

●Two series of Hill Kinetics equations can be estimated through a single Hill Kinetics equation.

●Protein which is translated from mRNA is proportional to corresponding mRNA in a previous time.

Consider gene lasR and luxR which are all constant genes, since AHL has a relatively low concentration compared to other substances, we assume that protein concentrations of PlasR and PluxR are irrelevant to AHLs. Therefore, concentrations of these proteins will reach to constant quantities. Moreover, certain mRNAs will also be constant. That is:

We have:

Then focus on compound LA1 and LA2. The binding process is far quicker than other reactions such as translation reactions, thus we assume that the compounds are Quasi-equilibrium:

We have:

Take the expressions to equation (6) (15), the functions are transformed to:

The transformed functions have nothing to do with feedback factors compared to original ones. Knowing that such feedback is essential for our system, we add additional feedback factors ka,kb manipulatively.

Then we want to come up a clear expression of A2c1,A1c2.They are generated by protein PluxI and PlasI, whose translation rate is much slower than the transcription rate. So we assume that certain mRNA is Quasi-equilibrium:

We have:

By ignoring the basal expression:

Knowing that protein PtetR is controlled by LA2 through Hill Function, we assume that concentration of mRNA is directly related to LA2:

Parameter KM4 is relevant to both KM1 and KM3, however accurate function depicting the relationship is unknown. We estimate the quantity of KM4 by multiplying KM1 and KM3. Noting LA1 and LA2 have already been expressed, we have:

Then we have trouble in expressing protein PluxI and PlasI. Since concentration of protein is the integrals of its mRNA, we assume that it is proportional to concentration of mRNA in a previous time. Thus we have a DDE function:

Equations concerning environmental AHLs remain unchanged.

The original ODEs can be transformed to a much simple DDEs by above deductions.



Parameters


Some of the parameters are derived from original ones, some of them are created to describe the new equations, and others are set for further testing.

Table 2 Parameters of DDEs

Results


General results

We coded the system in MATLAB. The result shows that every signal AHL is oscillating.

Figure 3 Signal AHL Oscillations

Stability analysis

In a nonlinear system, proper value ranges of the parameters are vital for producing a periodic oscillating solution. It’s not difficult to find the nonlinear part in our model, thus, we sought to analyze the sensitivity of Hill parameters and thoroughly reveal the internal relationship between Hill parameters and the system’s robustness. For simplicity, here we denoteKM1/ρ1 , KM4/ρ2 and kTL2,kTL4 by km1,km2and β1, β2 respectively. Without loss of generality, we determined to search the affection to the system’s robustness caused by the fluctuation of binary parameter (km22). After observing distinct oscillation mode (we choose two types of concentration of AHLs in environment as our observed object) and its phase trajectory, we depicted the bifurcation of our system.

Set basic parameters as follows:

ka=0,kb=0,μ=10min-1,γ=2.3min-11=30nm/min,km1=1,n1=n2=2,p=0.2,arab=20,t1=t2=0,τ1=6,τ2=10.

We changed the binary parameter (β2,km2) and got the result below.

Figure 4 C12 and C6 when (β2,km2)=(30,0.1)

Figure 5 C12 and C6 when (β2,km2)=(120,0.1)

Figure 6 C12 and C6 when (β2,km2)=(60,0.3)

After simulating the system at different parameters, we recorded several critical points for oscillation and made a table as follows.

Table 3 Critical points (β2,km2) for oscillation

Depicting those critical points on an axis, we immediately got the bifurcation line of parameters (β2,km2), which indicates the parameters’ value range when our system can oscillate stably, being marked in ‘bistable’.

Figure 7 Bifurcation Analysis on (β2,km2)

Proportion of cell volume

In actual vivo experiment, volumes of two separate cells might not be the same. We simulate the system by changing the proportion of two types of cells. What we know from the simulation is that cell proportions only affect the amplitude of signal molecules’ oscillation, but no influence to the period or stability of oscillation. The following graphics are drawn under cell proportion 1 and 0.7.

Figure 8 oscillations under different cell proportions

We can find that under different cell proportions, our system can oscillate at different amplitude, and varying in cell proportions may also lead to the change of oscillation period.

Period adjustment

As figure 2 indicates, we expected to adjust oscillation period by adding signal molecule aTc into our system. A small molecule as aTc is, it can easily bind to protein TetR and quickly depletes TetR in cell 2, which results in reduction of TetR net production rate, and indirectly, the protein’s inhibition on promoter5 is crippled. Thus, it would take a longer time for our system to reach each threshold, which is equivalent to prolonging the time delay τ2 in our simplified model. So we can deduce that changing the amount of aTc added into the system in precise model is equivalent to varying the time delay τ2 in simplified model.

Simulation results under distinct τ2 are presented as follows.

Figure 9 Oscillation cycle’s regulation

The result is exactly what we expected, which clearly demonstrates our system can truly be controlled by adding in external signal molecules.

Phase adjustment

In cell 2’s gene circuit, we designed a promoter induced by arabinose, marked by promoter 6 (see in figure 2). Promoter 6 is in a suppressed state until being induced by adding arabinose, and after the inhibition is relieved, signal molecule 3O12HSL will be generated extra. We can also analyze the differential equations describing dA1C2/dt, when adding arabinose during period from t1 to t2, dA1C2/dt contains an extra item arab?(t>t1)*(t2). The parameter arab can reflect the rate of adding arabinose nonlinearly. Here we set arab=20,t1=80,t2=120 and simulation result is presented as follow.

Figure 10 Oscillation phase’s regulation