Team:CIDEB-UANL Mexico/math overview

From 2014hs.igem.org

(Difference between revisions)
 
(25 intermediate revisions not shown)
Line 342: Line 342:
<div class="container-text">
<div class="container-text">
-
  <p><b><h2>When Biology meets Math</h2>
+
  <p><b><h2>When Biology meets Math</h2></b></p>
-
<p>&nbsp;</p>
+
 
-
</b>For  iGEM projects, the use of mathematical models is necessary to predict the  behavior of a biological machine, representing the quantitative relations  between two or more variables involved in the function and expression of a gene  or a set of genes in organisms like <i>E. coli</i>.
+
<p>In an iGEM project, mathematical models are necessary to predict the  behavior of a biological machine, representing the quantitative relations  between two or more variables involved in the function and expression of a gene  or a set of genes in organisms like <i>E. coli</i>.</p>
  <p>Our  team decided to use a deterministic model to simulate and represent the  function of the four modules, assuming that the variables (mRNA and protein  concentrations) adopt a continuous behavior and obey kinetic rules that can be  represented with constant values.</p>
  <p>Our  team decided to use a deterministic model to simulate and represent the  function of the four modules, assuming that the variables (mRNA and protein  concentrations) adopt a continuous behavior and obey kinetic rules that can be  represented with constant values.</p>
-
<p>This  type of mathematical model is used to include variables that considerate both,  the gene expression and the physiological cycles (chemical process, transport  of proteins, etc.). By the use of traditional differential equations we were  able to construct the description and analyze the behavior of mRNAs and protein  production. We were focused in determine through a mathematical model the ideal  behavior of the four modules independently.
+
<p>This  type of mathematical model is used to include variables that considerate both,  the gene expression and the physiological cycles (chemical process, transport  of proteins, etc.). By using traditional differential equations, the description was able to be constructed and it also permitted the analysis of the behavior of mRNAs and protein  production. The focus was to determine through a mathematical model the ideal  behavior of the four modules independently.
<p><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/project_capture">Capture</a></p><p><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/project_aroma">A</b>roma</a></p><p><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/project_resistance">Resistance </a></p><p><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/project_union">Union</a></p>
<p><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/project_capture">Capture</a></p><p><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/project_aroma">A</b>roma</a></p><p><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/project_resistance">Resistance </a></p><p><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/project_union">Union</a></p>
<p><b><h2>Deterministic Modelling: Equations and Parameters</h2></b></p>
<p><b><h2>Deterministic Modelling: Equations and Parameters</h2></b></p>
-
<p>As it was previously stated, it was necessary to work with different equations, focused in the production and degradation rate of proteins, according to the length of the genes in each module. This was performed in order to obtain both, the concentration rates of mRNA and protein based on system (gene) length and protein length respectively.</p>
+
<p>As it was previously stated, it was necessary to work with different equations focused in the production and degradation rate of mRNAs and proteins, according to the length of the genes in each module. This was performed in order to obtain both, the concentration rates of mRNA and protein based on system (gene) length and protein length respectively.</p>
<br><p><b>Equations</b></p>
<br><p><b>Equations</b></p>
Line 376: Line 376:
\end{equation}<br>
\end{equation}<br>
-
<br><p>Where: <b><i>“&#945;<sub>2 </sub>  [mRNA]”</i></b> means the translation rate of a protein based on the amount of mRNA available; <b><i>“d<sub>2</sub>[P]</i></b> represents the degradation rate of that protein; and <b><i>“f<sub>post</sub>"</i></b> the posttranslational variables which affect the production of the final (functional) protein.<p>
+
<br><p>Where: <b><i>“&#945;<sub>2 </sub>  [mRNA]”</i></b> means the translation rate of a protein based on the amount of mRNA available; <b><i>“d<sub>2</sub>[P]</i></b> represents the degradation rate of that protein; and <b><i>“f<sub>post</sub>"</i></b> the post-translational variables which affect the production of the final (functional) protein.<p>
<br><p><b>Parameters</b></p>
<br><p><b>Parameters</b></p>
-
<p>To determine both the transcription <b>(1)</b> and translation <b>(2)</b> rates, the parameters from Singapore 2008 iGEM team, at wich<i> E. coli</i> carries out transcription and translation were used ; as well as the degradation rates from Beijing PKU 2009. The parameters <i></i> from Singapore 2008 iGEM were used, assuming a transcription speed of <i>70nt/s</i> and a translation speed of <i>40aa/s</i>. The speeds were multiplied  by 60, because  minutes were used as units in the simulations of the modules. We assumed that <i>E. coli</i> division time was 30min based on PKU Beijing 2009 iGEM team.</p>
+
<p>To determine both the transcription (1) and translation (2) rates, the parameters from <a href="https://2008.igem.org/Team:NTU-Singapore">Singapore 2008</a> team at which <i> E. coli</i> carries out transcription and translation were used ; as well as the degradation rates from <a href="https://2009.igem.org/Team:PKU_Beijing">PKU Beijing  2009</a> team. The parameters <i></i> from <a href="https://2008.igem.org/Team:NTU-Singapore">Singapore 2008</a> were used, assuming a transcription speed of <i>70nt/s</i> and a translation speed of <i>40aa/s</i>. The speeds were multiplied  by 60 because  minutes were used as units in the simulations of the modules. It was determined that <i>E. coli's</i> division time was 30min, based on <a href="https://2009.igem.org/Team:PKU_Beijing">PKU Beijing 2009</a> team.</p>
<br>
<br>
Line 397: Line 397:
\end{equation}<br>
\end{equation}<br>
\begin{equation}
\begin{equation}
-
\large    d_{1} =  \frac{1}{half-life} + \frac{1}{division time} min
+
\large    d_{1} =  \frac{1}{half-life(min)} + \frac{1}{division time (min)}
\end{equation}<br>
\end{equation}<br>
\begin{equation}
\begin{equation}
-
\large    d_{2} =  \frac{1}{half-life} + \frac{1}{division time} min
+
\large    d_{2} =  \frac{1}{half-life(min)} + \frac{1}{division time (min)}
\end{equation}<br>
\end{equation}<br>
\begin{equation}
\begin{equation}
-
\large    d_{1} =  \frac{1}{half-life} + \frac{1}{30} min
+
\large    d_{1} =  \frac{1}{half-life(min)} + \frac{1}{30min}  
\end{equation}<br>
\end{equation}<br>
\begin{equation}
\begin{equation}
-
\large    d_{2} =  \frac{1}{half-life} + \frac{1}{30} min
+
\large    d_{2} =  \frac{1}{half-life(min)} + \frac{1}{30min}
\end{equation}<br>
\end{equation}<br>
Line 415: Line 415:
<font size="2">
<font size="2">
-
<p>● David Dibden, J. G. (2005). In vivo cycling of the Escherichia coli transcription factor FNR between active and inactive states. <i>Microbiology</i>, 4063-4070.</p>
+
<p>● BERGANT, Martina  N. M. (2010). <i>Modification of Human Papillomavirus Minor Capsid Protein L2 by Sumoylation Journal of Virology</i>, 11585-11589.</p>
-
<p>● Douglas Selinger, R. M. (2003). Global RNA Half-Life Analysis in Escherichia coli Reveals Positional Patterns of Transcript Degradation. <i>Genome Research</i>, 216-223.</p>
+
<p>● DIBDEN, David, J. G. (2005). In vivo cycling of the Escherichia coli transcription factor FNR between active and inactive states. <i>Microbiology</i>, 4063-4070.</p>
-
<p>● iGEM CIDEB UANL. (2012). <i>Modelling: Equations.</i> Retrieved from: <a href="https://2013hs.igem.org/Team:CIDEB-UANL_Mexico/Math-Equations"></a> </p>
+
<p>● HARE, James  K. T. (1991). Mechanisms of plasma membrane protein degradation: Recycling proteins are degraded more rapidly than those confined to the cell surface. <i>PNAS</i>, 5902-5906.</p>
-
<p>● James Hare, K. T. (1991). Mechanisms of plasma membrane protein degradation: Recycling proteins are degraded more rapidly than those confined to the cell surface. <i>PNAS</i>, 5902-5906.</p>
+
<p>● iGEM CIDEB UANL. (2012). <i>Modelling: Equations.</i>Retrieved from: <a href="https://2012hs.igem.org/Team:CIDEB-UANL_Mexico/Math/Overview">https://2012hs.igem.org/Team:CIDEB-UANL_Mexico/Math/Overview</a>.</p>
-
<p>● Martina Bergant, N. M. (2010). Modification of Human Papillomavirus Minor Capsid Protein L2 by Sumoylation <i>Journal of Virology</i>, 11585-11589.</p>
+
<p>● NTU Singapore. (2008). <i>Modelling: Parameters</i>. Retrieved from: <a href="https://2008.igem.org/Team:NTU-Singapore/Modelling/Parameter">https://2008.igem.org/Team:NTU-Singapore/Modelling/Parameter</a>.</p>
-
<p>● NTU Singapore. (2008). <i>Modelling: Parameters</i> Retrieved from: <a href="https://2008.igem.org/Team:NTU-Singapore/Modelling/Parameter">https://2008.igem.org/Team:NTU-Singapore/Modelling/Parameter</p>
+
<p>● PKU Beijing. (2009).<i>Modelling: Parameters</i>. Retrieved from:   <a href="https://2009.igem.org/Team:PKU_Beijing/Modeling/Parameters">https://2009.igem.org/Team:PKU_Beijing/Modeling/Parameters</a>.</p>
-
<p>● PKU Beijing. (2009).<i>Modelling: Parameters.</i> Retrieved from: https://2009.igem.org/Team:PKU_Beijing/Modeling/Parameters</p>
+
<p>● SELINGER, Douglas  R. M. (2003). Global RNA Half-Life Analysis in Escherichia coli Reveals Positional Patterns of Transcript Degradation. <i>Genome Research</i>, 216-223.</p>
-
<p>● Vorackova Irena, S. S. (2011). Purification of proteins containing zinc finger domains using Immobilized Metal Ion Affinity Chromatography. <i>Protein Expression and Purification</i>, 88-95</p>
+
<p>● VORACKOVA Irena, S. S. (2011). Purification of proteins containing zinc finger domains using Immobilized Metal Ion Affinity Chromatography. <i>Protein Expression and Purification</i>, 88-95.</p>
 +
 
<div style="text-align: right;"><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/math_overview#"><font color="blue">Return to the Top</font></a></p></div>
<div style="text-align: right;"><a href="https://2014hs.igem.org/Team:CIDEB-UANL_Mexico/math_overview#"><font color="blue">Return to the Top</font></a></p></div>
 +
</div>
</div>

Latest revision as of 01:43, 21 June 2014

iGEM CIDEB 2014 - Project

Overview

When Biology meets Math

In an iGEM project, mathematical models are necessary to predict the behavior of a biological machine, representing the quantitative relations between two or more variables involved in the function and expression of a gene or a set of genes in organisms like E. coli.

Our team decided to use a deterministic model to simulate and represent the function of the four modules, assuming that the variables (mRNA and protein concentrations) adopt a continuous behavior and obey kinetic rules that can be represented with constant values.

This type of mathematical model is used to include variables that considerate both, the gene expression and the physiological cycles (chemical process, transport of proteins, etc.). By using traditional differential equations, the description was able to be constructed and it also permitted the analysis of the behavior of mRNAs and protein production. The focus was to determine through a mathematical model the ideal behavior of the four modules independently.

Capture

Aroma

Resistance

Union

Deterministic Modelling: Equations and Parameters

As it was previously stated, it was necessary to work with different equations focused in the production and degradation rate of mRNAs and proteins, according to the length of the genes in each module. This was performed in order to obtain both, the concentration rates of mRNA and protein based on system (gene) length and protein length respectively.


Equations

mRNA

Generally, to describe the amount of mRNA produced over t time, the equation implemented is shown below:


\begin{equation} \large \frac{d[mRNA]}{dt} = \alpha_{1} \cdot f_{y} - d_{1}[mRNA] \end{equation}

Where: “α1 means the transcription rate of a given gene; “fy represents a regulatory function (if there is) that can activate or inhibit the system; and “d1 [mRNA]” the degradation rate of the mRNA produced.

Protein

The same happens with the protein production, but differs in the formula because it also involves a post-translational variable:



\begin{equation} \large \frac{d[P]}{dt} = \alpha_{2} \cdot[mRNA] - d_{2}[P] - f_{post} \end{equation}

Where: “α2 [mRNA]” means the translation rate of a protein based on the amount of mRNA available; “d2[P]” represents the degradation rate of that protein; and “fpost" the post-translational variables which affect the production of the final (functional) protein.


Parameters

To determine both the transcription (1) and translation (2) rates, the parameters from Singapore 2008 team at which E. coli carries out transcription and translation were used ; as well as the degradation rates from PKU Beijing 2009 team. The parameters from Singapore 2008 were used, assuming a transcription speed of 70nt/s and a translation speed of 40aa/s. The speeds were multiplied by 60 because minutes were used as units in the simulations of the modules. It was determined that E. coli's division time was 30min, based on PKU Beijing 2009 team.



\begin{equation} \large \alpha_{1} = \frac{transcription speed}{gene length (nt)} \end{equation}
\begin{equation} \large \alpha_{2} = \frac{translation speed}{protein length (aa)} \end{equation}
\begin{equation} \large \alpha_{1} = \frac{70 \frac {nt}{s} \cdot(60)}{gene length (nt)} \end{equation}
\begin{equation} \large \alpha_{2} = \frac{40 \frac {aa}{s} \cdot(60)}{protein length (aa)} \end{equation}
\begin{equation} \large d_{1} = \frac{1}{half-life(min)} + \frac{1}{division time (min)} \end{equation}
\begin{equation} \large d_{2} = \frac{1}{half-life(min)} + \frac{1}{division time (min)} \end{equation}
\begin{equation} \large d_{1} = \frac{1}{half-life(min)} + \frac{1}{30min} \end{equation}
\begin{equation} \large d_{2} = \frac{1}{half-life(min)} + \frac{1}{30min} \end{equation}

Bibliography/References

● BERGANT, Martina N. M. (2010). Modification of Human Papillomavirus Minor Capsid Protein L2 by Sumoylation Journal of Virology, 11585-11589.

● DIBDEN, David, J. G. (2005). In vivo cycling of the Escherichia coli transcription factor FNR between active and inactive states. Microbiology, 4063-4070.

● HARE, James K. T. (1991). Mechanisms of plasma membrane protein degradation: Recycling proteins are degraded more rapidly than those confined to the cell surface. PNAS, 5902-5906.

● iGEM CIDEB UANL. (2012). Modelling: Equations.Retrieved from: https://2012hs.igem.org/Team:CIDEB-UANL_Mexico/Math/Overview.

● NTU Singapore. (2008). Modelling: Parameters. Retrieved from: https://2008.igem.org/Team:NTU-Singapore/Modelling/Parameter.

● PKU Beijing. (2009).Modelling: Parameters. Retrieved from: https://2009.igem.org/Team:PKU_Beijing/Modeling/Parameters.

● SELINGER, Douglas R. M. (2003). Global RNA Half-Life Analysis in Escherichia coli Reveals Positional Patterns of Transcript Degradation. Genome Research, 216-223.

● VORACKOVA Irena, S. S. (2011). Purification of proteins containing zinc finger domains using Immobilized Metal Ion Affinity Chromatography. Protein Expression and Purification, 88-95.

iGEM CIDEB 2014 - Footer