Team:Carnegie Mellon/SensorModel

From 2014.igem.org

(Difference between revisions)
 
(98 intermediate revisions not shown)
Line 53: Line 53:
                     <li><a href="https://2014.igem.org/Team:Carnegie_Mellon/Our_Sensor">Our Sensor</a></li>
                     <li><a href="https://2014.igem.org/Team:Carnegie_Mellon/Our_Sensor">Our Sensor</a></li>
                     <li><a href="https://2014.igem.org/Team:Carnegie_Mellon/Achievements">Achievements</a></li>
                     <li><a href="https://2014.igem.org/Team:Carnegie_Mellon/Achievements">Achievements</a></li>
-
                    <li><a href="https://2014.igem.org/Team:Carnegie_Mellon/Future_Plans">Future Plans</a></li>
 
                   </ul>
                   </ul>
                 </li>                <li class="dropdown">
                 </li>                <li class="dropdown">
Line 65: Line 64:
                 </li>
                 </li>
<li class="dropdown">
<li class="dropdown">
-
                   <a href="#" class="dropdown-toggle" data-toggle="dropdown">Outreach <b class="caret"></b></a>
+
                   <a href="#" class="dropdown-toggle" data-toggle="dropdown">Policies & Practices <b class="caret"></b></a>
                   <ul class="dropdown-menu">
                   <ul class="dropdown-menu">
                     <li><a href="https://2014.igem.org/Team:Carnegie_Mellon/Our_Projects">Our Projects</a></li>
                     <li><a href="https://2014.igem.org/Team:Carnegie_Mellon/Our_Projects">Our Projects</a></li>
Line 98: Line 97:
<h1><center>Purpose</center> </h1>
<h1><center>Purpose</center> </h1>
-
<p> <center> The wet-lab model was written in BioNetGen Language, a rule-based modeling language that is useful for generating differential equations from a description of how various biological components interact with one another. The wet-lab model was constructed from both data found in the literature and data from the wet-lab itself. This allows the model to run various simulations of the experiment under different conditions in order to find the most optimal conditions for running the experiment. Also if the wet-lab run of the experiment is not producing any measurable results, the model can be used to identify exactly which component of the wet-lab experiment could be the problem. <center></p>
+
<p align="justify"> The bacterial cell model was written in the BioNetGen Language, a rule-based modeling language that is useful
 +
 
 +
for generating differential equations from a description of how various biological components interact with one  
 +
 
 +
another. The model was constructed from both data found in the literature and experimental data from the lab.  
 +
 
 +
This allows the model to run various simulations of the experiment under different conditions in order to find  
 +
 
 +
the optimal conditions for running the experiment. Also if the experiments are not producing any measurable  
 +
 
 +
results, the model can be used to identify which component of the experiment could be the problem. The model was run in Rule-Bender, an environment which is dedicated to running, analyzing, visualizing, and debugging BioNetGen Language models.</p>
<hr>
<hr>
<h2> <center>The Outline</center></h2>
<h2> <center>The Outline</center></h2>
-
<p align="left"><center> The model captures a total of 17 different interactions.</br></br>  
+
<p align="center"> <img src="https://static.igem.org/mediawiki/2014/8/80/Contact_Map.png"> </p>
-
1. The rate at which mRNA T7-Intein is transcribed.</br>
+
 
 +
<p align="left"> Rule-Bender provides users with a visual contact map of how each molecular component interfaces with one another. The above figure is a contact map for the T7-Intein complex, and was generated by Rule-Bender prior to the simulation of the model.The T7-Intein complex is composed of the N terminal T7 RNAP - N terminal intein – Estrogen Ligand Binding Domain – C terminal
 +
intein – C terminal T7 RNAP.</br></br>
 +
The model captures a total of 17 different interactions:</br></br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 +
1. The rate at which mRNA T7-Intein and YFP is transcribed.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
2. The rate at which mRNA T7-Intein is degraded.</br>
2. The rate at which mRNA T7-Intein is degraded.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
3. The rate at which T7-Intein is translated.</br>
3. The rate at which T7-Intein is translated.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
4. The rate at which T7-Intein is degraded.</br>
4. The rate at which T7-Intein is degraded.</br>
-
5. The rate at which T7 polymerase is degraded.</br>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 +
5. The rate at which T7 RNAP polymerase is degraded.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
6. The rate at which estrogen enters the cell.</br>
6. The rate at which estrogen enters the cell.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
7. The rate at which estrogen leaves the cell.</br>
7. The rate at which estrogen leaves the cell.</br>
-
8. The rate at which estrogen binds to the intein.</br>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-
9. The rate at which estrogen disassociates from the intein.</br>
+
8. The rate at which estrogen binds to the intein-LBD.</br>
-
10. The rate at which the intein is spliced out and T7 is formed.</br>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-
11. The rate at which T7 binds to the promoter of plasmid 2.</br>
+
9. The rate at which estrogen disassociates from the intein-LBD.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 +
10. The rate at which the intein is spliced out and T7 RNAP is formed.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 +
11. The rate at which T7 RNAP binds to the T7 promoter of plasmid 2.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
12. The rate at which mRNA RFP is transcribed.</br>
12. The rate at which mRNA RFP is transcribed.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
13. The rate at which mRNA RFP is degraded.</br>
13. The rate at which mRNA RFP is degraded.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
14. The rate at which RFP is translated.</br>
14. The rate at which RFP is translated.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
15. The rate at which RFP is degraded.</br>
15. The rate at which RFP is degraded.</br>
 +
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
16. The rate at which YFP is translated.</br>
16. The rate at which YFP is translated.</br>
-
17. The rate at which YFP is degraded.</br></br></center></p>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-
<p> <center>To run the model simply open the latest version of Rulebender, go to the simulation tab, provide the file path for the file you wish to run, and hit run.</center>
+
17. The rate at which YFP is degraded.
 +
</br></br></p>
 +
<p> To run the model simply open the latest version of Rulebender, go to the simulation tab, provide the file path for the file you wish to run, and hit run.
</p>
</p>
<hr>
<hr>
<h2> <center>The Parameters</center></h2>
<h2> <center>The Parameters</center></h2>
-
<table border="1" cellspacing="0" cellpadding="0">
+
<table border="1" cellspacing="0" cellpadding="10" align="center">
     <tbody>
     <tbody>
         <tr>
         <tr>
-
             <td width="180" valign="top">
+
             <td width="250" valign="top">
                 <p align="center">
                 <p align="center">
                     <strong>Parameter</strong>
                     <strong>Parameter</strong>
                 </p>
                 </p>
             </td>
             </td>
-
             <td width="180" valign="top">
+
             <td width="250" valign="top">
                 <p align="center">
                 <p align="center">
                     <strong>Value</strong>
                     <strong>Value</strong>
                 </p>
                 </p>
             </td>
             </td>
-
             <td width="180" valign="top">
+
             <td width="250" valign="top">
                 <p align="center">
                 <p align="center">
                     <strong>Reasoning</strong>
                     <strong>Reasoning</strong>
Line 155: Line 186:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     The experiment can be carried out under various temperatures. Temperature affects the rate of intein splicing.
+
                     The experiment can be carried out under various temperatures. Temperature affects the rate of intein splicing and RFP levels of
 +
                    fluorescence.
                 </p>
                 </p>
             </td>
             </td>
Line 179: Line 211:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     rfp_copies
+
                     RFP<sub>copies</sub>
                 </p>
                 </p>
             </td>
             </td>
Line 189: Line 221:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Number of rfp containing plasmids per cell.
+
                     Number of RFP containing plasmids per cell (pSB3K3).
                 </p>
                 </p>
             </td>
             </td>
Line 196: Line 228:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     int_copies
+
                     Int<sub>copies</sub>
                 </p>
                 </p>
             </td>
             </td>
Line 206: Line 238:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Number of intein containing plasmids per cell.
+
                     Number of intein containing plasmids per cell (pSB1C3).
                 </p>
                 </p>
             </td>
             </td>
Line 213: Line 245:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     estrogen_out
+
                     Estrogen<sub>out</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     ppt * 5.6e-02 nM
+
                     ppt * 5.6*10<sup>-2</sup> nM
                 </p>
                 </p>
             </td>
             </td>
Line 230: Line 262:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     estrogen_in
+
                     Estrogen<sub>in</sub>
                 </p>
                 </p>
             </td>
             </td>
Line 247: Line 279:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     t7_intein
+
                     T7-Intein
                 </p>
                 </p>
             </td>
             </td>
Line 264: Line 296:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     t7_estrogen_u
+
                     T7-Estrogen(U)
                 </p>
                 </p>
             </td>
             </td>
Line 281: Line 313:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     t7_estrogen_s
+
                     T7-Estrogen(S)
                 </p>
                 </p>
             </td>
             </td>
Line 298: Line 330:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     t7_nc
+
                     T7-NC
                 </p>
                 </p>
             </td>
             </td>
Line 315: Line 347:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     plas_int
+
                     Plas<sub>Int</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     int_copies * 1.0e-03 nM
+
                     Int<sub>copies</sub> * 1.0*10<sup>-3</sup> nM
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Concentration of T7-Intein plasmids.
+
                     Concentration of T7-Intein plasmids (pSB1C3).
                 </p>
                 </p>
             </td>
             </td>
Line 332: Line 364:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     plas_rfp
+
                     Plas<sub>RFP</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     rfp_copies * 1.0e-03 nM
+
                     RFP<sub>copies</sub> * 1.0*10<sup>-3</sup> nM
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Concentration of T7-RFP plasmids.
+
                     Concentration of T7-RFP plasmids (pSB3K3).
                 </p>
                 </p>
             </td>
             </td>
Line 349: Line 381:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     t7_rfp
+
                     T7-RFP
                 </p>
                 </p>
             </td>
             </td>
Line 359: Line 391:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Assume initial concentration of T7 bound to plasmid 2 is 0 nM.
+
                     Assume initial concentration of T7 RNAP bound to plasmid 2 is 0 nM.
                 </p>
                 </p>
             </td>
             </td>
Line 366: Line 398:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     mrna_rfp
+
                     mRNA<sub>RFP</sub>
                 </p>
                 </p>
             </td>
             </td>
Line 383: Line 415:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     prot_rfp
+
                     Prot<sub>RFP</sub>
                 </p>
                 </p>
             </td>
             </td>
Line 400: Line 432:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     mrna_int
+
                     mRNA<sub>Int</sub>
                 </p>
                 </p>
             </td>
             </td>
Line 417: Line 449:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     prot_yfp
+
                     Prot<sub>YFP</sub>
                 </p>
                 </p>
             </td>
             </td>
Line 434: Line 466:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_e
+
                     k<sub>e</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     1.7e-02 s<sup>-1</sup>
+
                     1.7*10<sup>-2</sup> s<sup>-1</sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
Rate at which estrogen diffuses through membrane. The diffusion coefficient of a steroid hormone in an aqueous phase = 10<sup>-13</sup> m<sup>2</sup>/s and the surface area of a bacteria is 6 * 10<sup>-12 </sup>m<sup>2</sup>. Thus (10<sup>-13</sup> m<sup>2</sup>/s) / (6 * 10                   <sup>-12 </sup>m<sup>2</sup>) = 1.7e-02 s<sup>-1</sup>
+
Rate at which estrogen diffuses through membrane. The diffusion coefficient of a steroid hormone in an aqueous phase = 10<sup>-13</sup> m                   <sup>2</sup>/s [1] and the surface area of a bacteria is 6 * 10<sup>-12 </sup>m<sup>2 </sup>[2]. Thus (10<sup>-13</sup> m<sup>2</sup>/s) /
 +
                    (6 * 10<sup>-12 </sup>m<sup>2</sup>) = 1.7*10<sup>-2</sup> s<sup>-1</sup>
                 </p>
                 </p>
             </td>
             </td>
Line 451: Line 484:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_e_t7_on
+
                     k<sub>e-T7(on)</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     1.3e-03 nM<sup>-1</sup>s<sup>-1</sup>
+
                     1.3*10<sup>-3</sup> nM<sup>-1</sup>s<sup>-1</sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which estrogen binds to the intein. The intein is the human estrogen receptor, so the value is from the literature.
+
                     Rate at which estrogen binds to the intein. The intein contains the ligand binding domain of the human estrogen receptor, so the value is
 +
                    from the literature [3].
                 </p>
                 </p>
             </td>
             </td>
Line 468: Line 502:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_e_t7_off
+
                     k<sub>e-T7(off)</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     1.2e-03 s<sup>-1</sup>
+
                     1.2*10<sup>-3</sup> s<sup>-1</sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which estrogen dissociates from the intein. Based on literature value of human estrogen receptor.
+
                     Rate at which estrogen dissociates from the intein. Based on literature value of human estrogen receptor [3].
                 </p>
                 </p>
             </td>
             </td>
Line 485: Line 519:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_splice
+
                     k<sub>splice</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     7.1e-04 s<sup>-1</sup>
+
                     7.1*10<sup>-4</sup> s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which intein splices out. Based on literature values of other inteins’ splicing kinetics.
+
                     Rate at which intein splices out. Based on literature values of other inteins’ splicing kinetics [4].
                 </p>
                 </p>
             </td>
             </td>
Line 502: Line 536:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_t7_rfp_on
+
                     k<sub>T7-RFP(on)</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     3.3e-01 nM<sup>-1</sup>s<sup>-1</sup>
+
                     3.3*10<sup>-1</sup> nM<sup>-1</sup>s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which T7 binds to the RFP plasmid. Value obtained from literature.
+
                     Rate at which T7 RNAP binds to the T7 promoter of the RFP plasmid. Value obtained from literature [5].
                 </p>
                 </p>
             </td>
             </td>
Line 519: Line 553:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_t7_rfp_off
+
                     k<sub>T7-RFP(off)</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     1.0e-01 s<sup>-1</sup>
+
                     1.0*10<sup>-1</sup> s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which T7 dissociates from the RFP plasmid. Value obtained from literature.
+
                     Rate at which T7 RNAP dissociates from the T7 promoter of the RFP plasmid. Value obtained from literature [5].
                 </p>
                 </p>
             </td>
             </td>
Line 536: Line 570:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_rna_rfp
+
                     k<sub>RNA<sub>RFP</sub></sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     8.8e-01 nM s<sup>-1</sup>
+
                     8.8*10<sup>-1</sup> nM s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which mRNA is synthesized from RFP plasmid. From iGEM team PKU ’09.
+
                     Rate at which mRNA is synthesized from RFP plasmid. From iGEM team PKU ’09 [6].
                 </p>
                 </p>
             </td>
             </td>
Line 553: Line 587:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_rna_deg_rfp
+
                     k<sub>RNA<sub>RFP</sub>(deg)</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     4.3e-03 s<sup>-1</sup>
+
                     4.3*10<sup>-3</sup> s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which mRNA RFP degrades. From iGEM team PKU ’09.
+
                     Rate at which mRNA RFP degrades. From iGEM team PKU ’09 [6].
                 </p>
                 </p>
             </td>
             </td>
Line 570: Line 604:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_prot_rfp
+
                     k<sub>Prot<sub>RFP</sub></sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     9.0e03 s<sup>-1</sup>
+
                     9.0*10<sup>-3</sup> s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which RFP protein is made. From iGEM team PKU ’09.
+
                     Rate at which RFP protein is made. From iGEM team PKU ’09 [6].
                 </p>
                 </p>
             </td>
             </td>
Line 587: Line 621:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_prot_deg_rfp
+
                     k<sub>Prot<sub>RFP</sub>(deg)</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     8.3e-04 s<sup>-1</sup>
+
                     8.3*10<sup>-4</sup> s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which RFP protein degrades. From iGEM team PKU ’09.
+
                     Rate at which RFP protein degrades. From iGEM team PKU ’09 [6].
                 </p>
                 </p>
             </td>
             </td>
Line 604: Line 638:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_rna_int
+
                     k<sub>RNA<sub>Int</sub></sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     2.6e-02 nM s<sup>-1</sup>
+
                     2.6*10<sup>-2</sup> nM s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which mRNA is synthesized from T7-Intein plasmid. From iGEM team PKU ’09.
+
                     Rate at which mRNA is synthesized from T7-Intein plasmid. From iGEM team PKU ’09 [6].
                 </p>
                 </p>
             </td>
             </td>
Line 621: Line 655:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_rna_deg_int
+
                     k<sub>RNA<sub>Int</sub>(deg)</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     4.3e-03 s<sup>-1</sup>
+
                     4.3*10<sup>-3</sup> s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which mRNA of T7-Intein degrades. From iGEM team PKU ’09.
+
                     Rate at which mRNA of T7-Intein degrades. From iGEM team PKU ’09 [6].
                 </p>
                 </p>
             </td>
             </td>
Line 638: Line 672:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_prot_int
+
                     k<sub>Prot<sub>Int</sub></sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     2.2e-03 s<sup>-1</sup>
+
                     2.2*10<sup>-3</sup> s<sup>-1</sup><sup></sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which T7-Intein is made from mRNA. From iGEM team PKU ’09.
+
                     Rate at which T7-Intein is made from mRNA. From iGEM team PKU ’09 [6].
                 </p>
                 </p>
             </td>
             </td>
Line 655: Line 689:
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     k_prot_deg_int
+
                     k<sub>Prot<sub>Int</sub>(deg)</sub>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     9.7e-04
+
                     9.7*10<sup>-4</sup> s<sup>-1</sup>
                 </p>
                 </p>
             </td>
             </td>
             <td width="180" valign="top">
             <td width="180" valign="top">
                 <p align="center">
                 <p align="center">
-
                     Rate at which T7 polymerase is degraded. From iGEM team PKU ’09.
+
                     Rate at which T7 polymerase is degraded. From iGEM team PKU ’09 [6].
                 </p>
                 </p>
             </td>
             </td>
Line 675: Line 709:
<p> Insight:
<p> Insight:
-
In order to accurately describe the behavior of certain interactions in the cell, it is necessary to use more than just simple rate constants. New functions which modify the rate constant based on cellular conditions must be introduced.</br>
+
In order to accurately describe the behavior of certain interactions in the cell, it is necessary to use more than just simple rate constants. New functions which modify the rate constant based on cellular conditions must be introduced.</br></br>
The rate at which mRNA RFP is synthesized is the product of the rate constant and the number of T7 polymerases bound to the RFP plasmid. Thus:</br></br>
The rate at which mRNA RFP is synthesized is the product of the rate constant and the number of T7 polymerases bound to the RFP plasmid. Thus:</br></br>
-
mRNA_RFP_synth() = k_rna_rfp * T7_RFP</br></br>
+
<center>mRNA<sub>RFP</sub> Synthesis = k<sub>RNA<sub>RFP</sub></sub> * T7_RFP</center></br></br>
The rate at which RFP protein is produced increases proportionally with the concentration of mRNA RFP present and decreases proportionally with the square root of the concentration of RFP currently present in the cell. Essentially the cell will make less RFP if there is already a lot of RFP currently present. Thus:</br></br>
The rate at which RFP protein is produced increases proportionally with the concentration of mRNA RFP present and decreases proportionally with the square root of the concentration of RFP currently present in the cell. Essentially the cell will make less RFP if there is already a lot of RFP currently present. Thus:</br></br>
-
protein_RFP_synth() = k_prot_rfp * mRNA_RFP * sqrt(1/(1 + prot_RFP))</br></br>
+
<center>Protein<sub>RFP</sub> Synthesis = k<sub>Prot<sub>RFP</sub></sub> * mRNA<sub>RFP</sub> * (1 + Prot<sub>RFP</sub>)
-
 
+
<sup><font size="2">-&frac12;</font></sup></center>
 +
</br></br></br>
The rate at which estrogen enters and exits the cell also depends on the concentration gradient of estrogen. If the concentration of estrogen outside the cell is much greater than the concentration of estrogen inside the cell, then estrogen will enter the cell at a much faster rate than if the concentration gradient is not as large. Thus:</br></br>
The rate at which estrogen enters and exits the cell also depends on the concentration gradient of estrogen. If the concentration of estrogen outside the cell is much greater than the concentration of estrogen inside the cell, then estrogen will enter the cell at a much faster rate than if the concentration gradient is not as large. Thus:</br></br>
-
e_in() = k_e * (E_out)/(E_out + E_in)</br>
+
<center>E<sub>in</sub> = k<sub>e</sub> * (Estrogen<sub>out</sub>) / (Estrogen<sub>out</sub> + Estrogen<sub>in</sub>)</br>
-
e_out() = k_e * (E_in)/(E_out + E_in)</br></br>
+
E<sub>out</sub> = k<sub>e</sub> * (Estrogen<sub>in</sub>) / (Estrogen<sub>out</sub> + Estrogen<sub>in</sub>)</center></br></br>
The rate at which the T7-Intein complex is produced increases proportionally with the concentration of mRNA T7-Intein present and decreases proportionally with the square root of the concentration of T7-Intein complex currently present in the cell. Essentially the cell will make less T7-Intein complex if there is already a lot of T7-Intein currently present. Thus:</br></br>
The rate at which the T7-Intein complex is produced increases proportionally with the concentration of mRNA T7-Intein present and decreases proportionally with the square root of the concentration of T7-Intein complex currently present in the cell. Essentially the cell will make less T7-Intein complex if there is already a lot of T7-Intein currently present. Thus:</br></br>
-
protein_INT_synth() = k_prot_int * mRNA_INT * sqrt(1/(1 + T7_unbound))</br></br>
+
<center>Protein<sub>Int</sub> Synthesis = k<sub>Prot<sub>Int</sub></sub> * mRNA<sub>Int</sub> * (1 + T7<sub>unbound</sub>)
 +
<sup><font size="2">-&frac12;</font></sup></center></br></br>
The rate at which YFP is produced increases proportionally with the concentration of mRNA T7-Intein present (as the YFP sequence is located on the same mRNA) and decreases proportionally with the square root of the concentration of YFP currently present in the cell. Essentially the cell will make less YFP if there is already a lot of YFP currently present. Thus:</br></br>
The rate at which YFP is produced increases proportionally with the concentration of mRNA T7-Intein present (as the YFP sequence is located on the same mRNA) and decreases proportionally with the square root of the concentration of YFP currently present in the cell. Essentially the cell will make less YFP if there is already a lot of YFP currently present. Thus:</br></br>
-
protein_INT_synth() = k_prot_int * mRNA_INT * sqrt(1/(1 + T7_unbound))</br></br>
+
 
 +
<center>Protein<sub>YFP</sub> Synthesis = k<sub>Prot<sub>Int</sub></sub> * mRNA<sub>Int</sub> * (1 + Prot<sub>YFP</sub>)
 +
<sup><font size="2">-&frac12;</font></sup></center></br></br>
The rate at which the intein splices out is dependent on the temperature of environment surrounding the cell. Since the cell is an E. coli cell, the optimal temperature for the surrounding environment is 37 degrees Celsius. This was confirmed in the wet-lab. Thus:</br></br>
The rate at which the intein splices out is dependent on the temperature of environment surrounding the cell. Since the cell is an E. coli cell, the optimal temperature for the surrounding environment is 37 degrees Celsius. This was confirmed in the wet-lab. Thus:</br></br>
-
int_Splice() = k_splice * (1/(1 + sqrt(abs(37 – Temp)))</br>
+
<center>Int Splice() = k<sub>splice</sub> * (1 + |37 - Temp| / &deg;C)
 +
<sup><font size="2">-&frac12;</font></sup></center>
 +
 
</p>
</p>
-
<p> The units of concentration in all of the graphs are in nM. The units of time in all of the graphs are in seconds. According to the literature enhanced GFP can be detected fairly reasonably at concentrations greater than 1 µm. Thus we can assume that regular RFP will be detected at concentrations of around 100 µm. All of the simulations were carried out at 37 degrees Celsius for the optimal results. The sensor should detect estrogen within a couple of hours, thus the simulation ends in a couple of hours.</p></br>
 
 +
<hr>
 +
<h2><p align="center">Results</p></h2>
 +
<p>The units of concentration in all of the graphs are in nM. The units of time in all of the graphs are in seconds. According to the literature enhanced GFP can be detected fairly reasonably at concentrations greater than 1 µm [7]. Thus we can assume that regular RFP will be detected at concentrations of around 100 µm. All of the simulations were carried out at 37 degrees Celsius for the optimal results. The sensor should detect estrogen within a couple of hours, thus the simulation ends in a couple of hours.</br></br>
 +
<center><i><b><font size="4">Scenario 1: Low amount of estrogen in the environment (< 100 ppt)</font></b></i></br></br>
-
<p>protein_INT_synth() = k_prot_int * mRNA_INT * sqrt(1/(1 + T7_unbound))</br></br>
+
<img src="https://static.igem.org/mediawiki/2014/e/e8/Scen_1_Wet_lab.PNG"></center></br></br>
-
The rate at which YFP is produced increases proportionally with the concentration of mRNA T7-Intein present (as the YFP sequence is located on the same mRNA) and decreases proportionally with the square root of the concentration of YFP currently present in the cell. Essentially the cell will make less YFP if there is already a lot of YFP currently present. Thus:</br></br>
+
The sensor is not sensitive enough to detect such low concentrations of estrogen in the environment, as the concentration of RFP produced never exceeds 100 µm. Also eventually the rate at which RFP degrades eclipses the rate at which RFP is produced thus making the sensor not very valuable at detecting incredibly low estrogenic concentrations.</br></br>
-
protein_INT_synth() = k_prot_int * mRNA_INT * sqrt(1/(1 + T7_unbound))</br></br>
+
<center><i><b><font size="4">Scenario 2: Ideal scenario. 1000 ppt of estrogen in the environment and all of the splicing kinetics work at a reasonably fast rate. Essentially the same scenario as if the T7 polymerase was constitutively expressed.</font></b></i></br></br>
-
The rate at which the intein splices out is dependent on the temperature of environment surrounding the cell. Since the cell is an E. coli cell, the optimal temperature for the surrounding environment is 37 degrees Celsius. This was confirmed in the wet-lab. Thus:</br></br>
+
<img src="https://static.igem.org/mediawiki/2014/5/55/Scen_2_Wet_Lab.PNG"></center></br></br>
 +
 
 +
This scenario is based on the wet-lab data we received for the constitutively expressed T7 polymerase. RFP is produced at detectable amounts and the sensor does its job of detecting estrogen in water.</br></br>
 +
 
 +
<center><i><b><font size="4">Scenario 3: The intein is not splicing properly. This is representative of what is occurring in the wet-lab at the moment when T7 is not constitutively expressed and its formation solely depends on the splicing of the intein.</font></b></i></br></br>
 +
 
 +
<img src="https://static.igem.org/mediawiki/2014/e/ea/Scen_3_Wet_Lab.PNG"></center></br></br>
 +
 
 +
Currently in the lab the sensor does not produce any measurable amount of RFP in the presence of estrogen. This graph represents what could possibly be occurring. Using the model, the bottleneck was identified to be the rate at which the intein splices. Thus the model saves the wet-lab biologists’ time by showing exactly which part of the system needs to be improved.
-
int_Splice() = k_splice * (1/(1 + sqrt(abs(37 – Temp)))</br>
 
</p>
</p>
-
<hr>
+
<h2><center>Code</center></h2>
-
<h2><p>Results</p></h2>
+
A commented, working version of the code can be found at <a href="https://static.igem.org/mediawiki/2014/9/91/WetLab_Model.pdf">
-
<p>The units of concentration in all of the graphs are in nM. The units of time in all of the graphs are in seconds. According to the literature enhanced GFP can be detected fairly reasonably at concentrations greater than 1 µm. Thus we can assume that regular RFP will be detected at concentrations of around 100 µm. All of the simulations were carried out at 37 degrees Celsius for the optimal results. The sensor should detect estrogen within a couple of hours, thus the simulation ends in a couple of hours.</br></br>INSERT GRAPHS</p>
+
E. coli Cell Model</a>
 +
 
 +
<h2><b><center>References</center></b></h2>
 +
<p align="left">1. E. Oren et al. Free Diffusion of Steroid Hormones Across Biomembranes: A Simplex Search with Implicit Solvent Model Calculations. Biophys J. Vol. 87, pp. 768–779. August 2004.</br>
 +
 
 +
2. http://bionumbers.hms.harvard.edu/search.aspx?log=y&task=searchbytrmorg&trm=surface+area+of+cell&org=&rpp=100</br>
 +
 
 +
3. R. L. Rich et al. Kinetic analysis of estrogen receptor/ligand interactions. PNAS. Vol. 99. No. 13. pp. 8562–8567, June 2002.</br>
 +
 
 +
4. S. Brenzel, T. Kurpiers, H. D. Mootz. Engineering artificially split inteins for applications in protein chemistry: biochemical characterization of the split Ssp DnaB intein and comparison to the split Sce VMA intein. Biochemistry. Vol 45(6), pp. 1571-1578, February 2006.</br>
 +
 
 +
5. http://www.jbc.org/content/271/48/30451/T1.expansion.html</br>
 +
 
 +
6. iGEM ’09 PKU-Beijing. https://2009.igem.org/Team:PKU_Beijing/Modeling/Parameters</br>
 +
 
 +
7. A. Furtado, and R. Henry. Measurement of green fluorescent protein concentration in single cells by image analysis. Vol. 310, pp. 84–92, May 2002.</br>
 +
 
 +
8. J. R. Faeder et al. Rule-based modeling of biochemical systems with BioNetGen. Methods Mol Biol. Vol. 500, pp. 113-167, 2009.</br>
 +
 
 +
9. RuleBender. http://visualizlab.org/rulebender</br>
 +
 
 +
10. BioNetGen. http://bionetgen.org/index.php/Main_Page</br>

Latest revision as of 02:45, 18 October 2014

Carousel Template · Bootstrap

Purpose

The bacterial cell model was written in the BioNetGen Language, a rule-based modeling language that is useful for generating differential equations from a description of how various biological components interact with one another. The model was constructed from both data found in the literature and experimental data from the lab. This allows the model to run various simulations of the experiment under different conditions in order to find the optimal conditions for running the experiment. Also if the experiments are not producing any measurable results, the model can be used to identify which component of the experiment could be the problem. The model was run in Rule-Bender, an environment which is dedicated to running, analyzing, visualizing, and debugging BioNetGen Language models.

The Outline

Rule-Bender provides users with a visual contact map of how each molecular component interfaces with one another. The above figure is a contact map for the T7-Intein complex, and was generated by Rule-Bender prior to the simulation of the model.The T7-Intein complex is composed of the N terminal T7 RNAP - N terminal intein – Estrogen Ligand Binding Domain – C terminal intein – C terminal T7 RNAP.

The model captures a total of 17 different interactions:

                       1. The rate at which mRNA T7-Intein and YFP is transcribed.
                       2. The rate at which mRNA T7-Intein is degraded.
                       3. The rate at which T7-Intein is translated.
                       4. The rate at which T7-Intein is degraded.
                       5. The rate at which T7 RNAP polymerase is degraded.
                       6. The rate at which estrogen enters the cell.
                       7. The rate at which estrogen leaves the cell.
                       8. The rate at which estrogen binds to the intein-LBD.
                       9. The rate at which estrogen disassociates from the intein-LBD.
                       10. The rate at which the intein is spliced out and T7 RNAP is formed.
                       11. The rate at which T7 RNAP binds to the T7 promoter of plasmid 2.
                       12. The rate at which mRNA RFP is transcribed.
                       13. The rate at which mRNA RFP is degraded.
                       14. The rate at which RFP is translated.
                       15. The rate at which RFP is degraded.
                       16. The rate at which YFP is translated.
                       17. The rate at which YFP is degraded.

To run the model simply open the latest version of Rulebender, go to the simulation tab, provide the file path for the file you wish to run, and hit run.

The Parameters

Parameter

Value

Reasoning

Temp

Variable (default value of 37 Celsius)

The experiment can be carried out under various temperatures. Temperature affects the rate of intein splicing and RFP levels of fluorescence.

ppt

Variable (default value of 1000)

Parts per trillion of estrogen in the water you are testing.

RFPcopies

15 copies/cell

Number of RFP containing plasmids per cell (pSB3K3).

Intcopies

2984 copies/cell

Number of intein containing plasmids per cell (pSB1C3).

Estrogenout

ppt * 5.6*10-2 nM

Concentration of estrogen outside the cell.

Estrogenin

0 nM

Concentration of estrogen inside the cell. Since e. coli is a prokaryote assume it is 0 nM.

T7-Intein

0 nM

Assume initial concentration of T7-Intein complex is 0 nM.

T7-Estrogen(U)

0 nM

Assume initial concentration of T7-Intein complex bound to estrogen is 0 nM.

T7-Estrogen(S)

0 nM

Assume initial concentration of spliced out Intein-Estrogen complex is 0 nM.

T7-NC

0 nM

Assume initial concentration of T7 polymerase is 0 nM.

PlasInt

Intcopies * 1.0*10-3 nM

Concentration of T7-Intein plasmids (pSB1C3).

PlasRFP

RFPcopies * 1.0*10-3 nM

Concentration of T7-RFP plasmids (pSB3K3).

T7-RFP

0 nM

Assume initial concentration of T7 RNAP bound to plasmid 2 is 0 nM.

mRNARFP

0 nM

Assume initial concentration of mRNA RFP is 0 nM.

ProtRFP

0 nM

Assume initial concentration of protein RFP is 0 nM.

mRNAInt

0 nM

Assume initial concentration of mRNA intein is 0 nM.

ProtYFP

0 nM

Assume initial concentration of protein YFP is 0 nM.

ke

1.7*10-2 s-1

Rate at which estrogen diffuses through membrane. The diffusion coefficient of a steroid hormone in an aqueous phase = 10-13 m 2/s [1] and the surface area of a bacteria is 6 * 10-12 m2 [2]. Thus (10-13 m2/s) / (6 * 10-12 m2) = 1.7*10-2 s-1

ke-T7(on)

1.3*10-3 nM-1s-1

Rate at which estrogen binds to the intein. The intein contains the ligand binding domain of the human estrogen receptor, so the value is from the literature [3].

ke-T7(off)

1.2*10-3 s-1

Rate at which estrogen dissociates from the intein. Based on literature value of human estrogen receptor [3].

ksplice

7.1*10-4 s-1

Rate at which intein splices out. Based on literature values of other inteins’ splicing kinetics [4].

kT7-RFP(on)

3.3*10-1 nM-1s-1

Rate at which T7 RNAP binds to the T7 promoter of the RFP plasmid. Value obtained from literature [5].

kT7-RFP(off)

1.0*10-1 s-1

Rate at which T7 RNAP dissociates from the T7 promoter of the RFP plasmid. Value obtained from literature [5].

kRNARFP

8.8*10-1 nM s-1

Rate at which mRNA is synthesized from RFP plasmid. From iGEM team PKU ’09 [6].

kRNARFP(deg)

4.3*10-3 s-1

Rate at which mRNA RFP degrades. From iGEM team PKU ’09 [6].

kProtRFP

9.0*10-3 s-1

Rate at which RFP protein is made. From iGEM team PKU ’09 [6].

kProtRFP(deg)

8.3*10-4 s-1

Rate at which RFP protein degrades. From iGEM team PKU ’09 [6].

kRNAInt

2.6*10-2 nM s-1

Rate at which mRNA is synthesized from T7-Intein plasmid. From iGEM team PKU ’09 [6].

kRNAInt(deg)

4.3*10-3 s-1

Rate at which mRNA of T7-Intein degrades. From iGEM team PKU ’09 [6].

kProtInt

2.2*10-3 s-1

Rate at which T7-Intein is made from mRNA. From iGEM team PKU ’09 [6].

kProtInt(deg)

9.7*10-4 s-1

Rate at which T7 polymerase is degraded. From iGEM team PKU ’09 [6].

Insights

Insight: In order to accurately describe the behavior of certain interactions in the cell, it is necessary to use more than just simple rate constants. New functions which modify the rate constant based on cellular conditions must be introduced.

The rate at which mRNA RFP is synthesized is the product of the rate constant and the number of T7 polymerases bound to the RFP plasmid. Thus:

mRNARFP Synthesis = kRNARFP * T7_RFP


The rate at which RFP protein is produced increases proportionally with the concentration of mRNA RFP present and decreases proportionally with the square root of the concentration of RFP currently present in the cell. Essentially the cell will make less RFP if there is already a lot of RFP currently present. Thus:

ProteinRFP Synthesis = kProtRFP * mRNARFP * (1 + ProtRFP) -½



The rate at which estrogen enters and exits the cell also depends on the concentration gradient of estrogen. If the concentration of estrogen outside the cell is much greater than the concentration of estrogen inside the cell, then estrogen will enter the cell at a much faster rate than if the concentration gradient is not as large. Thus:

Ein = ke * (Estrogenout) / (Estrogenout + Estrogenin)
Eout = ke * (Estrogenin) / (Estrogenout + Estrogenin)


The rate at which the T7-Intein complex is produced increases proportionally with the concentration of mRNA T7-Intein present and decreases proportionally with the square root of the concentration of T7-Intein complex currently present in the cell. Essentially the cell will make less T7-Intein complex if there is already a lot of T7-Intein currently present. Thus:

ProteinInt Synthesis = kProtInt * mRNAInt * (1 + T7unbound) -½


The rate at which YFP is produced increases proportionally with the concentration of mRNA T7-Intein present (as the YFP sequence is located on the same mRNA) and decreases proportionally with the square root of the concentration of YFP currently present in the cell. Essentially the cell will make less YFP if there is already a lot of YFP currently present. Thus:

ProteinYFP Synthesis = kProtInt * mRNAInt * (1 + ProtYFP) -½


The rate at which the intein splices out is dependent on the temperature of environment surrounding the cell. Since the cell is an E. coli cell, the optimal temperature for the surrounding environment is 37 degrees Celsius. This was confirmed in the wet-lab. Thus:

Int Splice() = ksplice * (1 + |37 - Temp| / °C) -½

Results

The units of concentration in all of the graphs are in nM. The units of time in all of the graphs are in seconds. According to the literature enhanced GFP can be detected fairly reasonably at concentrations greater than 1 µm [7]. Thus we can assume that regular RFP will be detected at concentrations of around 100 µm. All of the simulations were carried out at 37 degrees Celsius for the optimal results. The sensor should detect estrogen within a couple of hours, thus the simulation ends in a couple of hours.

Scenario 1: Low amount of estrogen in the environment (< 100 ppt)



The sensor is not sensitive enough to detect such low concentrations of estrogen in the environment, as the concentration of RFP produced never exceeds 100 µm. Also eventually the rate at which RFP degrades eclipses the rate at which RFP is produced thus making the sensor not very valuable at detecting incredibly low estrogenic concentrations.

Scenario 2: Ideal scenario. 1000 ppt of estrogen in the environment and all of the splicing kinetics work at a reasonably fast rate. Essentially the same scenario as if the T7 polymerase was constitutively expressed.



This scenario is based on the wet-lab data we received for the constitutively expressed T7 polymerase. RFP is produced at detectable amounts and the sensor does its job of detecting estrogen in water.

Scenario 3: The intein is not splicing properly. This is representative of what is occurring in the wet-lab at the moment when T7 is not constitutively expressed and its formation solely depends on the splicing of the intein.



Currently in the lab the sensor does not produce any measurable amount of RFP in the presence of estrogen. This graph represents what could possibly be occurring. Using the model, the bottleneck was identified to be the rate at which the intein splices. Thus the model saves the wet-lab biologists’ time by showing exactly which part of the system needs to be improved.

Code

A commented, working version of the code can be found at E. coli Cell Model

References

1. E. Oren et al. Free Diffusion of Steroid Hormones Across Biomembranes: A Simplex Search with Implicit Solvent Model Calculations. Biophys J. Vol. 87, pp. 768–779. August 2004.
2. http://bionumbers.hms.harvard.edu/search.aspx?log=y&task=searchbytrmorg&trm=surface+area+of+cell&org=&rpp=100
3. R. L. Rich et al. Kinetic analysis of estrogen receptor/ligand interactions. PNAS. Vol. 99. No. 13. pp. 8562–8567, June 2002.
4. S. Brenzel, T. Kurpiers, H. D. Mootz. Engineering artificially split inteins for applications in protein chemistry: biochemical characterization of the split Ssp DnaB intein and comparison to the split Sce VMA intein. Biochemistry. Vol 45(6), pp. 1571-1578, February 2006.
5. http://www.jbc.org/content/271/48/30451/T1.expansion.html
6. iGEM ’09 PKU-Beijing. https://2009.igem.org/Team:PKU_Beijing/Modeling/Parameters
7. A. Furtado, and R. Henry. Measurement of green fluorescent protein concentration in single cells by image analysis. Vol. 310, pp. 84–92, May 2002.
8. J. R. Faeder et al. Rule-based modeling of biochemical systems with BioNetGen. Methods Mol Biol. Vol. 500, pp. 113-167, 2009.
9. RuleBender. http://visualizlab.org/rulebender
10. BioNetGen. http://bionetgen.org/index.php/Main_Page

Retrieved from "http://2014.igem.org/Team:Carnegie_Mellon/SensorModel"