Team:UT-Tokyo/Modeling

From 2013.igem.org

(Difference between revisions)
(Prototype team page)
 
(41 intermediate revisions not shown)
Line 1: Line 1:
-
<!-- *** What falls between these lines is the Alert Box!  You can remove it from your pages once you have read and understood the alert *** -->
+
{{:Team:UT-Tokyo/Template/Style}}
-
 
+
<html>
<html>
-
<div id="box" style="width: 700px; margin-left: 137px; padding: 5px; border: 3px solid #000; background-color: #fe2b33;">
+
<link
-
<div id="template" style="text-align: center; font-weight: bold; font-size: large; color: #f6f6f6; padding: 5px;">
+
href='http://fonts.googleapis.com/css?family=Source+Sans+Pro:200,300,300italic,400,600,700'
-
This is a template page. READ THESE INSTRUCTIONS.
+
rel='stylesheet' type='text/css'>
-
</div>
+
 
-
<div id="instructions" style="text-align: center; font-weight: normal; font-size: small; color: #f6f6f6; padding: 5px;">
+
<div id="global_wrapper">
-
You are provided with this team page template with which to start the iGEM season.  You may choose to personalize it to fit your team but keep the same "look." Or you may choose to take your team wiki to a different level and design your own wiki.  You can find some examples <a href="https://2008.igem.org/Help:Template/Examples">HERE</a>.
+
<div id="header">
-
</div>
+
<div class="logo">
-
<div id="warning" style="text-align: center; font-weight: bold; font-size: small; color: #f6f6f6; padding: 5px;">
+
<img
-
You <strong>MUST</strong>  have all of the pages listed in the menu below with the names specified.  PLEASE keep all of your pages within your teams namespace. 
+
src="https://static.igem.org/mediawiki/2013/9/93/Clock_ecoli_paper5.png"
-
</div>
+
height="120px">
-
</div>
+
</div>
 +
<div class="wrapper">
 +
<div class="box">We conducted the following modeling!!</div>
 +
</div>
</html>
</html>
 +
{{:Team:UT-Tokyo/Template/Header1|subpagename=MODELING}}
 +
<html>
 +
<div id="contents">
-
<!-- *** End of the alert box *** -->
+
<p class="ini">Our concept of Multicellular Analog Clock is based on qualitative assumption such as how negative feedback loop behaves, how AHL diffuses, and so on. To ascertain our system can function as an analogue clock, namely, to confirm the feasibility of our cell-cell communication included gene circuit, and to deepen understanding of behavior of the system, we conducted the following simulation.</p>
 +
<h1>DDE model derivation</h1>
-
{| style="color:#1b2c8a;background-color:#0c6;" cellpadding="3" cellspacing="1" border="1" bordercolor="#fff" width="62%" align="center"
 
-
!align="center"|[[Team:UT-Tokyo|Home]]
 
-
!align="center"|[[Team:UT-Tokyo/Team|Team]]
 
-
!align="center"|[https://igem.org/Team.cgi?year=2013&team_name=UT-Tokyo Official Team Profile]
 
-
!align="center"|[[Team:UT-Tokyo/Project|Project]]
 
-
!align="center"|[[Team:UT-Tokyo/Parts|Parts Submitted to the Registry]]
 
-
!align="center"|[[Team:UT-Tokyo/Modeling|Modeling]]
 
-
!align="center"|[[Team:UT-Tokyo/Notebook|Notebook]]
 
-
!align="center"|[[Team:UT-Tokyo/Safety|Safety]]
 
-
!align="center"|[[Team:UT-Tokyo/Attributions|Attributions]]
 
-
|}
 
 +
<p class="ini">Using DDE (delay-differential equation) analysis, we can conduct gene network simulation which accounts for delayed production of proteins. It is true that ODE (ordinary differential equation) analysis enables us to simulate rough behavior of complex gene networks, but if you want to simulate temporally more detailed behavior, ODE analysis is insufficient. Since when cells produce proteins including transcription factors, they must pass through processes of transcription, translation, maturation, transcription rate of genes reflects production of transcription factor after these processes' time delay. For this reason, we derived delay-differential equation in reference to [1].</p>
 +
<p>For the reason stated above, we developed the model based on partial delay-differential equations with Hill functions that captured the activation and repression of protein synthesis, and Michaelis-Menten kinetics that captured enzyme catalysis. Finally we derived the following set of delay-differential equation model for concentrations of AiiA (A), TetR (T), LuxI (I), internal AHL (H<sub>i</sub>), and external AHL (H<sub>e</sub>)</p>
-
If you choose to include a '''Modeling''' page, please write about your modeling adventures here.  This is not necessary but it may be a nice list to include.
+
<img src="https://static.igem.org/mediawiki/2013/f/f1/PDDE.png"  class="figure">
 +
<p><b>Figure 1</b> delay-differential equations</p>
 +
 
 +
 
 +
 
 +
<p>We did not include an equation for LuxR assuming that it is constitutively produced at a constant level, which is adopted based on some references in which model development is performed on the same assumption ([1] and [2]). For model simplicity, concentrations of four proteins are approximated as having the same behavior. In equation (1), [H<sub>i</sub>]<sub>τ</sub> and [T]<sub>τ</sub> mean corresponding proteins' concentrations before τ minutes. The parameters are described in the following.
 +
</p>
 +
 
 +
<a href="https://static.igem.org/mediawiki/2013/d/da/Parameter.png"><img src="https://static.igem.org/mediawiki/2013/d/da/Parameter.png" class="figure"></a>
 +
<p><b>Figure 2</b> Parameters</p>
 +
 
 +
<h1>Parameter sensitivity</h1>
 +
<p class="ini">To figure out the sensitive parameters in our system, we varied the parameters and investigated how the behavior was affected. This parameter sensitivity analysis let us know that LuxI degradation rate (γ<sub>I</sub>) is sensitive, and the range of the parameter which realizes sustained moving of clock hand is limited. Therefore, it is thought that when our system doesn't experimentally function as we expect, we may improve it by varying this parameter in diverse ways.</p>
 +
 
 +
<img src="https://static.igem.org/mediawiki/2013/4/46/Fig1_1_1_1_1_1_1_1_1_1_0.5_1_1_1_1_1_1.png"  class="figure">
 +
<p><b>Figure 3</b> The simulation result with γ<sub>I</sub> 0.5 times of standard value</p>
 +
 
 +
<img src="https://static.igem.org/mediawiki/2013/4/44/Fig1_1_1_1_1_1_1_1_1_1_1_1_1_1_1_1_1.png"  class="figure">
 +
<p><b>Figure 4</b> The simulation result with γ<sub>I</sub> standard value</p>
 +
 
 +
<img src="https://static.igem.org/mediawiki/2013/5/53/Fig1_1_1_1_1_1_1_1_1_1_1.4_1_1_1_1_1_1.png"  class="figure">
 +
<p><b>Figure 5</b> The simulation result with γ<sub>I</sub> 1.4 times of standard value</p>
 +
 
 +
<h1>Reference</h1>
 +
 
 +
<p class="ini">[1] Danino, T., Mondragon-Palomino, O., Tsimring, L., & Hasty, J. (2010). A synchronized quorum of genetic clocks. Nature, 463(7279), 326-330.</p>
 +
 
 +
<p class="ini">[2] Tian, T., & Burrage, K. (2006). Stochastic models for regulatory networks of the genetic toggle switch. Proceedings of the National Academy of Sciences, 103(22), 8372-8377.</p>
 +
 
 +
</div>
 +
</div>
 +
 
 +
<div id="footer"></div>
 +
 
 +
</div>
 +
</body>
 +
 
 +
</html>

Latest revision as of 11:54, 14 October 2013

           MODELING
       

Our concept of Multicellular Analog Clock is based on qualitative assumption such as how negative feedback loop behaves, how AHL diffuses, and so on. To ascertain our system can function as an analogue clock, namely, to confirm the feasibility of our cell-cell communication included gene circuit, and to deepen understanding of behavior of the system, we conducted the following simulation.

DDE model derivation

Using DDE (delay-differential equation) analysis, we can conduct gene network simulation which accounts for delayed production of proteins. It is true that ODE (ordinary differential equation) analysis enables us to simulate rough behavior of complex gene networks, but if you want to simulate temporally more detailed behavior, ODE analysis is insufficient. Since when cells produce proteins including transcription factors, they must pass through processes of transcription, translation, maturation, transcription rate of genes reflects production of transcription factor after these processes' time delay. For this reason, we derived delay-differential equation in reference to [1].

For the reason stated above, we developed the model based on partial delay-differential equations with Hill functions that captured the activation and repression of protein synthesis, and Michaelis-Menten kinetics that captured enzyme catalysis. Finally we derived the following set of delay-differential equation model for concentrations of AiiA (A), TetR (T), LuxI (I), internal AHL (Hi), and external AHL (He)

Figure 1 delay-differential equations

We did not include an equation for LuxR assuming that it is constitutively produced at a constant level, which is adopted based on some references in which model development is performed on the same assumption ([1] and [2]). For model simplicity, concentrations of four proteins are approximated as having the same behavior. In equation (1), [Hi]τ and [T]τ mean corresponding proteins' concentrations before τ minutes. The parameters are described in the following.

Figure 2 Parameters

Parameter sensitivity

To figure out the sensitive parameters in our system, we varied the parameters and investigated how the behavior was affected. This parameter sensitivity analysis let us know that LuxI degradation rate (γI) is sensitive, and the range of the parameter which realizes sustained moving of clock hand is limited. Therefore, it is thought that when our system doesn't experimentally function as we expect, we may improve it by varying this parameter in diverse ways.

Figure 3 The simulation result with γI 0.5 times of standard value

Figure 4 The simulation result with γI standard value

Figure 5 The simulation result with γI 1.4 times of standard value

Reference

[1] Danino, T., Mondragon-Palomino, O., Tsimring, L., & Hasty, J. (2010). A synchronized quorum of genetic clocks. Nature, 463(7279), 326-330.

[2] Tian, T., & Burrage, K. (2006). Stochastic models for regulatory networks of the genetic toggle switch. Proceedings of the National Academy of Sciences, 103(22), 8372-8377.