Difference between revisions of "Team:Tokyo Tech/Modeling Details"

Line 6: Line 6:
 
/********************** All style of this page ************************/
 
/********************** All style of this page ************************/
 
body{
 
body{
    background: #052d5d;
 
 
     background-image: url("https://static.igem.org/mediawiki/2016/8/83/Tokyo_Tech_Background.png");
 
     background-image: url("https://static.igem.org/mediawiki/2016/8/83/Tokyo_Tech_Background.png");
background-size: 75%;
 
// background-size: 60%;
 
 
}
 
}
 
#globalWrapper, #content{
 
#globalWrapper, #content{
Line 15: Line 12:
 
}
 
}
 
#content{
 
#content{
//    background: #052d5d;
 
 
//    background-image: url("https://static.igem.org/mediawiki/2016/8/83/Tokyo_Tech_Background.png");
 
//    background-image: url("https://static.igem.org/mediawiki/2016/8/83/Tokyo_Tech_Background.png");
    margin: 0 auto;
 
    border: 0px solid #052d5d;
 
 
}
 
}
 
/**********************************************************************/
 
/**********************************************************************/
Line 26: Line 20:
 
     min-width: 880px;
 
     min-width: 880px;
 
     max-width: 1100px;
 
     max-width: 1100px;
//    background: #052d5d;
 
//    border:1px solid #052d5d;
 
 
}
 
}
 
.container{
 
.container{
 
     margin: 5px auto;
 
     margin: 5px auto;
 
     background: #ffffff;
 
     background: #ffffff;
 +
border: 1px solid #583612;
 
}
 
}
 
.container_top h1{
 
.container_top h1{
Line 42: Line 35:
 
}
 
}
 
.container_header span{
 
.container_header span{
     border-bottom: 2px solid #052d5d;
+
     border-bottom: 2px solid #583612;
 
}
 
}
 
.container_contents{
 
.container_contents{
Line 57: Line 50:
 
     font-family: 'Century', serif;
 
     font-family: 'Century', serif;
 
}
 
}
.on{
+
.non_dotted_list{
display: block;
+
list-style: none;
-webkit-transiton: all 0.6s ease;
+
transition: all 0.6s ease;
+
}
+
.off{
+
display: none;
+
-webkit-transiton: all 0.6s ease;
+
transition: all 0.6s ease;
+
 
}
 
}
 
--></style>
 
--></style>
Line 72: Line 58:
 
<div id="main_contents">
 
<div id="main_contents">
 
<div id="page_header" class="container container_top">
 
<div id="page_header" class="container container_top">
<h1 align="center">Detail Descriptions</h1>
+
<h1 align="center">Model Development</h1>
 +
<div id="page_header_contents" class="container_contents">
 +
<p class="normal_text">To simulate the cell-cell communication system, we developed an ordinary differential equation model.
 +
The following sentences describe how the equations were developed.
 +
And in this page we expound not only on the model with the Maz system, which we selected as the best TA system for our project, but also on the one with the Yaf system, which we chose as an alternative.</p>
 +
</div><!-- /page_header_contents -->
 
</div><!-- /page_header -->
 
</div><!-- /page_header -->
<div id="development" class="container">
+
<div id="modeling_maz_system" class="container">
<div id="development_header" class="container_header">
+
<div id="modeling_maz_header" class="container_header">
<h2><span>Model Development</span></h2>
+
<h2><span>1. Maz System</span></h1>
</div><!-- /contents_header -->
+
</div><!-- /modeling_maz_header -->
<div id="development_contents" class="container_contents">
+
<div id="modeling_maz_contents" class="container_contents">
<p class="normal_text">
+
<ul id="modeling_maz_list" class="non_dotted_list">
 
+
<li><h2>1.1. Cell Population</h2>
</p>
+
<p class="normal_text">The equations above describe how cells grow in the culture.
</div><!-- /contents_contents -->
+
Equations (1), (2) and (3) describe the populations of Snow White, the Queen and the Prince. (3) is described by the logistic growth equation, but (1) and (2) are represented by the growth inhibition by MazF dimers.
</div><!-- /contents -->
+
This factor is designed so that its value is small when the concentration of MazF dimers is low, and its value converges to 1 when the concentration of MazF dimers is high.</p>
 
+
</li><!-- /1.1. Cell Population -->
 +
<li><h2>1.2. Maz System</h2>
 +
<ul id="modeling_maz_system" class="non_dotted_list">
 +
<li><h3>1.2.1. Expression of Maz System</h3>
 +
<p class="normal_text">After translation, MazE and MazF each form  an stable dimer which can be activated to exert its function.</p>
 +
<p class="normal_text">Two MazE dimers sandwich the MazF dimer, forming MazF2-MazE2-MazF2 heterohexamers and suppressing the toxicity of the MazF dimers.</p>
 +
<p class="normal_text">The mRNAs of Snow White and the Queen decrease by their original degradation and by the cleavage at ACA sequences by MazF dimers.</p>
 +
<p class="normal_text">Applying mass action kinetic laws, we obtain the following set of differential equations.</p>
 +
<p class="normal_text">Equations (1) and (8) describe the concentration of mRNAs under the AHL inducing promoters.
 +
Thus, they comprise terms of production by leaky expressions of promoters, terms of production by Hill function dependent on the concentration of C12/C4, terms of original degradation and terms of degradation from cleavage at ACA sequences by MazF dimers.
 +
Since Equations (2), (3), (5), (6), (7), (9), (10), (12), (13) and (14) describe the concentrations of complexes, mainly they comprise terms of production and terms of binding and dissociation.</p>
 +
</li><!-- /1.2.1. Expression of Maz System -->
 +
<li><h3>1.2.2. Cleavage by MazF dimers</h3>
 +
<p class="normal_text">MazF dimers recognize and cleave ACAs in mRNAs, thus acting as Toxin.</p>
 +
<p class="normal_text">We estimated the rate of recognitions of ACA sequences by MazF dimers at 1-(1-f)^n, where the number of ACA sequences in mRNA.</p>
 +
<p class="normal_text">Then, we expressed the rate of degradation by MazF dimers in F(1-(1-f)^{f_{mRNA}}) and obtain the following set of differential equations.</p>
 +
<p class="normal_text">The equations above comprise terms of production, terms of original degradation and terms of degradation from cleavage at ACA sequences by MazF dimers.
 +
</p>
 +
</li><!-- /1.2.2. Cleavage by MazF dimers -->
 +
</ul><!-- /modeling_maz_system -->
 +
</li><!-- /1.2. Maz System -->
 +
<li><h2>1.3. Signal Molecules</h2>
 +
<p class="normal_text">Snow White expresses RhlI under Plux induced by C12, the Queen expresses LasI under Prhl induced by C4 and the Prince expresses AmiE under Plux induced by C12.</p>
 +
<p class="normal_text">The mRNAs of Snow White and the Queen decrease from original degradation and the cleavage at ACA sequences by MazF dimers.
 +
On the other hand, those of the Prince don’t have any MazF gene so they decrease from only original degradation.</p>
 +
<p class="normal_text">After translation, C12AHL and C4 are enzymatically synthesized by LasI and RhlI from some substrates respectively.
 +
For simplicity, we assumed that the amount of substrates is sufficient so that the C12AHL / C4 synthesis rate per cell is estimated to be proportional to the LasI and RhlI concentrations.</p>
 +
<p class="normal_text">C4 decreases from original degradation meanwhile C12AHL decreases from both original degradation and degradation by AmiE, which Prince products.</p>
 +
<p class="normal_text">Applying mass action kinetic laws, we obtain the following set of differential equations.</p>
 +
<p class="normal_text">Equations (1), (4) and (7) describe the concentrations of mRNAs under the AHL inducing promoters.
 +
Thus, they comprise terms of production by leaky expressions of promoters, terms of production by Hill function dependent on the concentration of C12/C4, terms of original degradation and terms of degradation from cleavage at ACA sequences by MazF dimers.</p>
 +
<p class="normal_text">The other ODEs describe how the concentrations of materials change in individuals, on the other hand (3), (6) describe the concentrations of C4 C12AHL in the whole culture medium.</p>
 +
</li>
 +
</ul><!-- /modeling_maz_list -->
 +
</div><!-- /modeling_maz_contents -->
 +
</div><!-- /modeling_maz_system -->
 +
<div id="modeling_yaf_system" class="container container_bottom">
 +
<div id="modeling_yaf_header" class="container_header">
 +
<h2><span>2. Yaf System (Alternative Design)</span></h2>
 +
</div><!-- /modeling_yaf_header -->
 +
<div id="modeling_yaf_contents" class="container_contents">
 +
<p class="normal_text">We also designed alternative design with Yaf system as our plan B.</p>
 +
<p class="normal_text">After translation, YafO and YafN exert their function alone.</p>
 +
<p class="normal_text">YafO forms an heterodimer with YafO and suppress the toxicity of YafO.</p>
 +
<p class="nprmal_text">Applying mass action kinetic laws, we obtain the following set of differential equations.</p>
 +
<p class="normal_text">Equations (1) and (6) describe the concentration of mRNAs under the AHL inducing promoters.
 +
Thus, they comprise terms of production by leaky expressions of promoters, terms of production by Hill function dependent on the concentration of C12/C4, terms of original degradation and terms of degradation by YafO.
 +
Since Equations (5) and (10) describe the concentrations of complexes, mainly they comprise terms of production and terms of binding dissociation.</p>
 +
</div><!-- /modeling_yaf_contents -->
 +
</div><!-- /modeling_yaf_system -->
 +
</div><!-- /main_contents -->
 
</body>
 
</body>
 
</html>
 
</html>
{{:Team:Tokyo_Tech/tpl/head|Parameters and Tools}}
 
{{:Team:Tokyo _Tecj/tpl/fortables}}
 
 
 
<html><article></html>
 
==Parameters==
 
No model is complete without parameters. Our exhaustive list of parameters are summarised in the table below.
 
 
{|class="wikitable sortable"
 
!Parameter!!Value!!Description!!Reference
 
|-
 
|g||0.0123||Growth rate of each cells||[ここにリンク先 Fitted to experimental data]
 
|-
 
|P<sub>max</sub>||3.6641||Carrying capacity||[ここにリンク先 Fitted to experimental data]
 
|-
 
|E<sub>DiMazF</sub>||0.101715||Effect of MazF dimer on growth rate of each cells||[ここにリンク先 Fitted to experimental data]
 
|-
 
|k||5 min<sup>-1</sup>||Transcriptional rate of mRNA under P<sub>tet</sub>||Assumed
 
|-
 
|leak<sub>P<sub>lux</sub></sub>||0.113 nMmin<sup>-1</sup>||Leakage of P<sub>lux</sub>||[ここにリンク先 Fitted to experimental data]
 
|-
 
|leak<sub>P<sub>rhl</sub></sub>||3.5 nMmin<sup>-1</sup>||Leakage of P<sub>rhl</sub>||[ここにリンク先 Fitted to experimental data]
 
|-
 
|κ<sub>P<sub>lux</sub></sub>||10.3 nM<sup>-1</sup>min<sup>-1</sup>||Maximum transcriptional rate of mRNA under P<sub>lux</sub>||[ここにリンク先 Fitted to experimental data]
 
|-
 
|κ<sub>P<sub>rhl</sub></sub>||9.73 nM<sup>-1</sup>min<sup>-1</sup>||Maximum transcriptional rate of mRNA under P<sub>rhl</sub>||[ここにリンク先 Fitted to experimental data]
 
|-
 
|n<sub>Lux</sub>||2.5||Hill coefficient for P<sub>lux</sub>||[ここにリンク先 Fitted to experimental data]
 
|-
 
|n<sub>Rhl</sub>||2.5||Hill coefficient for P<sub>rhl</sub>||[ここにリンク先 Fitted to experimental data]
 
|-
 
|K<sub>mLux</sub>||78.32 nM||Lumped parameter for the Lux system||Estimated from Literature
 
|-
 
|K<sub>mRhl</sub>||1969 nM||Lumped parameter for the Rhl system||Estimated from Literature
 
|-
 
|F<sub>DiMazF</sub>||5 nM<sup>-1</sup>min<sup>-1</sup>||Cutting rate at ACA sequences on mRNA by MazF dimer||Assumed
 
|-
 
|f||0.409213||The probability of distinction of ACA sequences on mRNA by MazF dimer||Estimated
 
|-
 
|f<sub>mRNA<sub>RFP</sub></sub>||10||The number of ACA arrays in mRNA<sub>RFP</sub>||Extraction of data
 
|-
 
|f<sub>mRNA<sub>GFP</sub></sub>||23||The number of ACA arrays in mRNA<sub>GFP</sub>||Extraction of data
 
|-
 
|f<sub>mRNA<sub>RhlI</sub></sub>||1||The number of ACA arrays in mRNA<sub>RhlI</sub>||Extraction of data
 
|-
 
|f<sub>mRNA<sub>LasI</sub></sub>||10||The number of ACA arrays in mRNA<sub>LasI</sub>||Extraction of data
 
|-
 
|f<sub>mRNA<sub>MazF</sub></sub>||2||The number of ACA arrays in mRNA<sub>MazF</sub>||Extraction of data
 
|-
 
|f<sub>mRNA<sub>MazE</sub></sub>||10||The number of ACA arrays in mRNA<sub>MazE</sub>||Extraction of data
 
|-
 
|α||5 min<sup>-1</sup>||Translation rate of protein||Assumed
 
|-
 
|k<sub>DiMazF</sub>||5.226124 nM<sup>-1</sup>min<sup>-1</sup>||Formation rate of MazF dimer||Literature
 
|-
 
|k<sub>-DiMazF</sub>||3.151962 min<sup>-1</sup>||Dissociation rate of MazF dimer||Literature
 
|-
 
|k<sub>DiMazE</sub>||6.759816 nM<sup>-1</sup>min<sup>-1</sup>||Formation rate of MazE dimer||Literature
 
|-
 
|k<sub>-DiMazE</sub>||3.808504 min<sup>-1</sup>||Dissociation rate of MazE dimer||Literature
 
|-
 
|k<sub>Hexa</sub>||3.170544 nM<sup>-1</sup>min<sup>-1</sup>||Formation rate of Maz hexamer||Literature
 
|-
 
|k<sub>-Hexa</sub>||8.344358 min<sup>-1</sup>||Dissociation rate of Maz hexamer||Literature
 
|-
 
|p<sub>Las</sub>||1 min<sup>-1</sup>||Production rate of 3OC12AHL by LasI||Assumed
 
|-
 
|p<sub>Rhl</sub>||1 min<sup>-1</sup>||Production rate of C4AHL by RhlI||Assumed
 
|-
 
|D||1 nM<sup>-1</sup>min<sup>-1</sup>||Decomposition rate of 3OC12AHL by AmiE||[Assumed
 
|-
 
|d||0.2773 min<sup>-1</sup>||Degradation rate of mRNA||Literature
 
|-
 
|d<sub>RFP</sub>||0.017 min<sup>-1</sup>||Degradation rate of RFP||Literature
 
|-
 
|d<sub>GFP</sub>||0.017 min<sup>-1</sup>||Degradation rate of GFP||Literature
 
|-
 
|d<sub>RhlI</sub>||0.0167 min<sup>-1</sup>||Degradation rate of RhlI||Literature
 
|-
 
|d<sub>LasI</sub>||0.0167 min<sup>-1</sup>||Degradation rate of LasI||Literature
 
|-
 
|d<sub>MazF</sub>||0.173286 min<sup>-1</sup>||Degradation rate of RhlI||[ここにリンク先 Fitted to experimental data]
 
|-
 
|d<sub>DiMazF</sub>||0.700037 min<sup>-1</sup>||Degradation rate of MazF dimer||[ここにリンク先 Fitted to experimental data]
 
|-
 
|d<sub>MazE</sub>||0.340082 min<sup>-1</sup>||Degradation rate of MazE||[ここにリンク先 Fitted to experimental data]
 
|-
 
|d<sub>DiMazE</sub>||0.166615 min<sup>-1</sup>||Degradation rate of MazE dimer||[ここにリンク先 Fitted to experimental data]
 
|-
 
|d<sub>Hexa</sub>||0.155088 min<sup>-1</sup>||Degradation rate of Maz hexamer||[ここにリンク先 Fitted to experimental data]
 
|-
 
|d<sub>C4AHL</sub>||0.000222 min<sup>-1</sup>||Degradation rate of C4AHL||Literature
 
|}
 
|d<sub>C12AHL</sub>||0.004 min<sup>-1</sup>||Degradation rate of 3OC12AHL||Literature
 
|}
 
|d<sub>AmiE</sub>||0.01 min<sup>-1</sup>||Degradation rate of mRNA||Assumed
 
|}
 
<html></article></html>
 
 
<html><article></html>
 

Revision as of 00:32, 13 October 2016

1. Maz System

  • 1.1. Cell Population

    The equations above describe how cells grow in the culture. Equations (1), (2) and (3) describe the populations of Snow White, the Queen and the Prince. (3) is described by the logistic growth equation, but (1) and (2) are represented by the growth inhibition by MazF dimers. This factor is designed so that its value is small when the concentration of MazF dimers is low, and its value converges to 1 when the concentration of MazF dimers is high.

  • 1.2. Maz System

    • 1.2.1. Expression of Maz System

      After translation, MazE and MazF each form an stable dimer which can be activated to exert its function.

      Two MazE dimers sandwich the MazF dimer, forming MazF2-MazE2-MazF2 heterohexamers and suppressing the toxicity of the MazF dimers.

      The mRNAs of Snow White and the Queen decrease by their original degradation and by the cleavage at ACA sequences by MazF dimers.

      Applying mass action kinetic laws, we obtain the following set of differential equations.

      Equations (1) and (8) describe the concentration of mRNAs under the AHL inducing promoters. Thus, they comprise terms of production by leaky expressions of promoters, terms of production by Hill function dependent on the concentration of C12/C4, terms of original degradation and terms of degradation from cleavage at ACA sequences by MazF dimers. Since Equations (2), (3), (5), (6), (7), (9), (10), (12), (13) and (14) describe the concentrations of complexes, mainly they comprise terms of production and terms of binding and dissociation.

    • 1.2.2. Cleavage by MazF dimers

      MazF dimers recognize and cleave ACAs in mRNAs, thus acting as Toxin.

      We estimated the rate of recognitions of ACA sequences by MazF dimers at 1-(1-f)^n, where the number of ACA sequences in mRNA.

      Then, we expressed the rate of degradation by MazF dimers in F(1-(1-f)^{f_{mRNA}}) and obtain the following set of differential equations.

      The equations above comprise terms of production, terms of original degradation and terms of degradation from cleavage at ACA sequences by MazF dimers.

  • 1.3. Signal Molecules

    Snow White expresses RhlI under Plux induced by C12, the Queen expresses LasI under Prhl induced by C4 and the Prince expresses AmiE under Plux induced by C12.

    The mRNAs of Snow White and the Queen decrease from original degradation and the cleavage at ACA sequences by MazF dimers. On the other hand, those of the Prince don’t have any MazF gene so they decrease from only original degradation.

    After translation, C12AHL and C4 are enzymatically synthesized by LasI and RhlI from some substrates respectively. For simplicity, we assumed that the amount of substrates is sufficient so that the C12AHL / C4 synthesis rate per cell is estimated to be proportional to the LasI and RhlI concentrations.

    C4 decreases from original degradation meanwhile C12AHL decreases from both original degradation and degradation by AmiE, which Prince products.

    Applying mass action kinetic laws, we obtain the following set of differential equations.

    Equations (1), (4) and (7) describe the concentrations of mRNAs under the AHL inducing promoters. Thus, they comprise terms of production by leaky expressions of promoters, terms of production by Hill function dependent on the concentration of C12/C4, terms of original degradation and terms of degradation from cleavage at ACA sequences by MazF dimers.

    The other ODEs describe how the concentrations of materials change in individuals, on the other hand (3), (6) describe the concentrations of C4 C12AHL in the whole culture medium.

2. Yaf System (Alternative Design)

We also designed alternative design with Yaf system as our plan B.

After translation, YafO and YafN exert their function alone.

YafO forms an heterodimer with YafO and suppress the toxicity of YafO.

Applying mass action kinetic laws, we obtain the following set of differential equations.

Equations (1) and (6) describe the concentration of mRNAs under the AHL inducing promoters. Thus, they comprise terms of production by leaky expressions of promoters, terms of production by Hill function dependent on the concentration of C12/C4, terms of original degradation and terms of degradation by YafO. Since Equations (5) and (10) describe the concentrations of complexes, mainly they comprise terms of production and terms of binding dissociation.