Team:UANL Mty-Mexico/Modeling

From 2013.igem.org

(Difference between revisions)
Line 14: Line 14:
     <script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
     <script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
-
<p><br><h2><a name="Biosensor"></a>Biosensor</h2><br></p>
 
-
 
-
<p>Our biosensor model is a continuation of the previously exposed transport and accumulation model with the addition of the production of luciferase as a function of ArsR|As concentration. In this way, we can relate the concentration of extracellular As to the production of luciferase. </p>
 
-
 
-
<p><br><h3><a name="ODEBIOSENSOR"></a>ODEs</h3></p><hr align="center" width="33%"/></br>
 
-
 
-
<p>The ODEs that describe the kinetics of the biosensor modules include all the assumptions presented in the <a href="https://2012.igem.org/Team:UANL_Mty-Mexico/Modeling/transport_and_accumulation#core_model">core model</a> and the equations presented for the <a href="https://2012.igem.org/Team:UANL_Mty-Mexico/Modeling/transport_and_accumulation#ODEANT"> accumulation and transport</a> module, along with their corresponding <a href="https://2012.igem.org/Team:UANL_Mty-Mexico/Modeling/transport_and_accumulation#parametersANT"> parameters</a>.
 
-
</p>
 
-
 
-
<p>Taking into account the equations from the accumulation and transport module, we introduce a generator for Luciferase:
 
-
</p>
 
-
<br><p><b><i>Luciferase mRNA</i></b></p>
 
\begin{equation}
\begin{equation}
\large \frac{d[mRNA_{Luc}]}{dt} = \alpha _{mLuc}\cdot (pro_{ars})\cdot(\frac{k_{D1}^{h_{1}}}{k_{D1}^{h_{1}}+[ArsR]^{h_{1}}})- \delta _{mRNA_{Luc}}[mRNA_{Luc}]
\large \frac{d[mRNA_{Luc}]}{dt} = \alpha _{mLuc}\cdot (pro_{ars})\cdot(\frac{k_{D1}^{h_{1}}}{k_{D1}^{h_{1}}+[ArsR]^{h_{1}}})- \delta _{mRNA_{Luc}}[mRNA_{Luc}]
\end{equation} </br></p>
\end{equation} </br></p>
-
<br><p><b><i>Luciferase protein
 
</p></b></i>
</p></b></i>
\begin{equation}
\begin{equation}
\large \frac{d[Luc]}{dt} = \alpha _{pLuc}\cdot[mRNA_{Luc}]- \delta _{Luc}[Luc]
\large \frac{d[Luc]}{dt} = \alpha _{pLuc}\cdot[mRNA_{Luc}]- \delta _{Luc}[Luc]
\end{equation} </p>
\end{equation} </p>
-
 
-
<p><br><h3><a name="Parbiosensor"></a>Parameters</h3><hr align="center" width="33%"/></p><hr align="center" width="33%"/></br>
 
-
 
-
<p>The parameters for the ODEs describing the production of luciferase mRNA and protein are shown in the next table.</p>
 
-
 
-
 
-
<br><p><b>Parameter table</b></p>
 
-
<br>
 
-
<center>
 
-
<table border="1" cellspacing="1" cellpadding="1">
 
-
  <tbody>
 
-
    <tr  style="background-color:#99CCFF">
 
-
      <td valign="top">
 
-
        <p>Parameter</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Description</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Value</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>References</p>
 
-
      </td>
 
-
    </tr>
 
-
 
 
-
  <tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>α<sub>mArsR</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Maximal transcription rate of ArsR</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>3.74 nM/min</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
  <tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>α<sub>mFluc</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Maximal transcription rate of Fluc</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>1.40 nM/min</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
  <tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>α<sub>pArsR</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Maximal translation rate of ArsR</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>9.74 nM/min</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
  <tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>α<sub>pFluc</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Maximal translation rate of Fluc</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>2.06 nM/min</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
  <tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>δ<sub>mRNAArsR</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Degradation rate of ArsR mRNA</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>2.16x10<sup>-1</sup> min<sup>-1</sup></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions, <a href="http://genome.cshlp.org/content/13/2/216.long">Selinger, <i>et al.</i> (2003)</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
  <tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>δ<sub>mRNAFluc</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Degradation rate of Fluc mRNA</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>3.37x10<sup>-2</sup> min<sup>-1</sup></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions, <a href="http://www.sciencemag.org/content/suppl/2011/03/16/science.1198817.DC1/Suter.SOM.pdf">Suter, <i>et al.</i> (2011), <a href="http://jgv.sgmjournals.org/content/88/1/286.full">Karetnivok and Lehto (2006).</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
  <tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>δ<sub>ArsR</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Degradation rate of ArsR</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>2.31x10<sup>-2</sup> min<sup>-1</sup></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
  <tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>δ<sub>Fluc</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Degradation rate of Fluc</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>2.60x10<sup>-2</sup> min<sup>-1</sup></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions, <a href="http://www.pnas.org/content/early/2009/02/10/0813345106.full.pdf">Auld, <i>et al. </i> (2008)</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
   
 
-
<tr  border="1" style="background-color:#CCCCCC">
 
-
      <td valign="top">
 
-
        <p>pro<sub>ars</sub></p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Concentration of ars promoter</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>aprox. 1*plasmid copy (nM)</p>
 
-
      </td>
 
-
      <td valign="top">
 
-
        <p>Assumptions</p>
 
-
      </td>
 
-
    </tr>
 
-
 
-
 
-
  </tbody>
 
-
</table>
 
-
</center>
 
-
 
-
 
-
<p><br><h3><a name="Parbiosensor"></a>Simulations</h3><hr align="center" width="33%"/></p><hr align="center" width="33%"/></br>
 
-
 
-
<p>For the constructed ODEs, we built a Simulink model. With the extracellular As set to zero and the initial conditions for all other variables set to zero, as well, we obtained the values for the variables ArsR and Fluc (both in protein and mRNA) at steady state and used them as initial conditions for further simulations:</p><p>
 
-
mRNAs: ArsR = 2.1058 nM; Fluc = 5 nM</p><p>
 
-
Proteins: ArsR = 887.7 nM; Fluc = 398.6 nM</p><p>
 
-
We also set all the numerical integrators to have lower saturation limits equal to zero, to be in concordance with the mass conservation law. Here are the results for simulations at extracellular As set to 0, 0.1, 1, 5 and 10 μM (0, 100, 1000, 5000, 10000 nM). </p>
 
-
 
-
<table class="image" align="center">
 
-
<caption align="bottom"><b>Figure 1.</b> Time vs Luciferase production (nM)</caption>
 
-
<tr><td><img src="https://static.igem.org/mediawiki/2012/3/3a/Arsenic_vs_flucUANL2012.jpg" style="width:800px;"></td></tr>
 
-
</table>
 
-
<br>
 
-
 
-
 
-
<br>
 
-
<p>The last simulation shows the dynamics of the Luc production at different As concentrations; now we want to know the dynamic range of our biosensor. In order to do this, we performed simulations at varying extracellular As concentrations (ranging 0, 0.1, 1, 5, 10, 50 and 100 μM) and recorded the level of Luc at 150 min. Next, we plotted the pairs of extracellular As concentration and Luc levels at 150 min; this pairs were then fitted to a rational curve using the Matlab fitting tool, so we ended with a formula for the theoretical relation between extracellular arsenic and Luc levels at 150 min
 
-
</p>
 
-
 
-
<table class="image" align="center">
 
-
<caption align="bottom"><b>Figure 2.</b> Extracellular arsenic vs Luc(nM)</caption>
 
-
<tr><td><img src="https://static.igem.org/mediawiki/2012/thumb/5/54/As_vs_flucUANL2012.jpg/800px-As_vs_flucUANL2012.jpg" style="width:800px;"></td></tr>
 
-
</table>
 
-
<br>
 
-
 
-
<br><p><b><i>Function relating extracellular arsenic and Luc
 
-
</p></b></i>
 
-
\begin{equation}
 
-
\large [Luc] = \frac{2816[As]^{2} + 316.6[As] - 607}{[As]^{2} + 10.24[As] - 3.521}
 
-
\end{equation} </p>
 
-
 
-
<p>The R<sup>2</sup> for this equation as a regresion of the actual data was 0.998.
 
-
 
-
<p><br><h2><a name="Simulink_biosensor"></a>Simulink model: biosensor</h2><br></p>
 
-
 
-
See our simulink model <a href="https://dl.dropbox.com/u/105944108/Biosensor_Rubengood.mdl">here</a>.
 
-
 
-
</div>
 
-
<div class="sidebar2">
 
-
    <ul class="nav">
 
-
<a href="https://2012.igem.org/Team:UANL_Mty-Mexico/Modeling"><li>Modeling</a></li>
 
-
<a href="https://2012.igem.org/Team:UANL_Mty-Mexico/Modeling/transport_and_accumulation"><li>Transport and accumulation</a></li>
 
-
<a href="#Biosensor"><li><b>>>Biosensor</b></a></li>
 
-
<a href="#ODEBIOSENSOR"><li><i>ODEs</i></a></li>
 
-
<a href="#Parbiosensor"><li><i>Parameters</i></a></li>
 
-
<a href="#Parbiosensor"><li><i>Simulations</i></a></li>
 
-
<a href="#Simulink_biosensor"><li><i>Simulink model: biosensor</i></a></li>
 
-
<a href="https://2012.igem.org/Team:UANL_Mty-Mexico/Modeling/Silica_binding"><li>Silica binding</a></li>
 
-
 
-
  </ul>
 
-
 
-
 
-
 
-
 
-
</div>
 
-
</html>
 
-
{{:Team:UANL_Mty-Mexico/Templates:Footer}}
 

Revision as of 19:41, 1 September 2013

UANL México

\begin{equation} \large \frac{d[mRNA_{Luc}]}{dt} = \alpha _{mLuc}\cdot (pro_{ars})\cdot(\frac{k_{D1}^{h_{1}}}{k_{D1}^{h_{1}}+[ArsR]^{h_{1}}})- \delta _{mRNA_{Luc}}[mRNA_{Luc}] \end{equation}

\begin{equation} \large \frac{d[Luc]}{dt} = \alpha _{pLuc}\cdot[mRNA_{Luc}]- \delta _{Luc}[Luc] \end{equation}