Team:CIDEB-UANL Mexico/math overview

From 2014hs.igem.org

(Difference between revisions)
 
(104 intermediate revisions not shown)
Line 324: Line 324:
</style>
</style>
<body>
<body>
 +
<!-- MathJax (LaTeX for the web) -->
 +
    <script type="text/x-mathjax-config">
 +
        MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}});
 +
        MathJax.Hub.Config({
 +
            TeX: {
 +
                equationNumbers: {  autoNumber: "AMS"  }
 +
            }
 +
        });
 +
    </script>
 +
    <script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
<div class="wrapper">
<div class="wrapper">
<div id="welcome" class="container">
<div id="welcome" class="container">
Line 332: Line 342:
<div class="container-text">
<div class="container-text">
-
<p style="font-size20px"><b>When Biology Meets Math</b></p>
+
<p><b><h2>When Biology meets Math</h2></b></p>
-
<p>For iGEM projects, the use of mathematical models is performed to describe the behavior of biological systems, representing the quantitative relations between two or more variables involved in the function of gene(s) in organisms as <i>E. coli</i>.</p>
+
<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 functioning 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 (Loera, 2012).</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 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>Then, we use the parameters for degradation rates for proteins and mRNAs from Beijing PKU 2009 iGEM team:</p>
+
<p><b><h2>Deterministic Modelling: Equations and Parameters</h2></b></p>
-
<p>As the protein was the fusion of two we need to search for each half-life. The half-life of membrane proteins range between 2 to 20 hours in E. coli (Hare, 1991), and as AIDA-I is a membrane protein its half-life must be between that range since it is not determined the specific half-life of AIDA.  To find the half-life of L2 we assumed it was 7.8 hours (Bergant, 2010). Bergant’s team made test with a homologous protein but found in the minor capsid of the Human Papillomavirus (HPV). Although the function of the L2 strand in HPV is viral, and in E. coli is ribosomal, both share similar structures and sequences. Once we have decided to use the half-life from the homologous L2 we determined to use it as the half-life for the fusion protein because it was between the range of AIDA-I, and also because it was the lower half-life assuming as E. coli start the L2 degradation, it would degrade the whole protein.</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>
-
<p>This type of mathematical model is used to involve variables that abstract both the gene expression and 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 proteins production.  We were focused in determine through a mathematical model of the four modules (one per each):<b> Capture, Aroma, Resistance</b> and <b>Union</b>.</p>
+
-
 
+
-
<p style="font-size20px"><b>Deterministic Modelling: Equations and Parameters</b></p>
+
-
 
+
-
<p>As it was established previously, is necessary work with different equations focusing in the production rate and degradation rate according to the length of genes in each module of the project. This is performed in order to obtain both the concentration rates of mRNA and protein based on system (gene) length and protein length respectively.</p>
+
-
 
+
-
<p><b>Equations</b></p>
+
<p>● <b>mRNA</b><p>
<p>● <b>mRNA</b><p>
Line 355: Line 360:
<p>Generally, to describe the amount of mRNA produced over <i>t</i> time, the equation implemented is shown below:</p>
<p>Generally, to describe the amount of mRNA produced over <i>t</i> time, the equation implemented is shown below:</p>
-
<center><p><img width=200 height=100 src="http://farm6.staticflickr.com/5581/14421576775_81d23038df_b.jpg"
+
<br>
-
align=center hspace=12 alt="IMG_0317"></p></center>
+
\begin{equation}
 +
\large \frac{d[mRNA]}{dt} = \alpha_{1} \cdot f_{y} - d_{1}[mRNA]
 +
\end{equation}<br>
-
<p>Where: <i>a1</i>” means the transcription rate of a givengene; <i>fy</i>” represents a regulatory function (if there is) that can activate or inhibit the system; and <i>d1[mRNA]</i>the degradation rate of the mRNA produced.</p>
+
<br><p>Where: <b><i>“&#945;<sub>1</sub>”</i></b> means the transcription rate of a given gene; <b><i>“f<sub>y</sub>”</i></b> represents a regulatory function (if there is) that can activate or inhibit the system; and <b><i>“d<sub>1 </sub>  [mRNA]</i></b> the degradation rate of the mRNA produced.</p>
<p>● <b>Protein</b><p>
<p>● <b>Protein</b><p>
-
<p>The same happens with the protein production, but differs in the formula because it also involves a posttranslational variable:</p>
+
<p>The same happens with the protein production, but differs in the formula because it also involves a post-translational variable:</p><br>
-
<p>Where: “<i>a2[mRNA]</i>” means the translation rate of a protein based on the amount of mRNA available; “<i>d2[P]</i>” represents the degradation rate of that protein; and “<i>f(post)</i>" the posttranslational variables which affect the production of the final (functional) protein.<p>
+
<br>
 +
\begin{equation}
 +
\large \frac{d[P]}{dt} = \alpha_{2} \cdot[mRNA] - d_{2}[P] - f_{post}
 +
\end{equation}<br>
-
<p><b>Parameters</b></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>
-
<p>To determine both the transcription (1) and translation (2) rates, we used the parameters from Singapore 2008 iGEM team, as well as the degradation rates from Beijing PKU 2009. We used the parameters for the speeds at which <i>E. coli</i> carries out transcription and translation from Singapore 2008 iGEM, assuming a transcription speed of <i>70nt/s</i> and a translation speed of <i>40aa/s</i>. We multiply the speeds by 60 because we use minutes 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>
+
<br><p><b>Parameters</b></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>
-
<p><b>References</b></p>
 
-
<p>● Douglas Selinger, R. M. (2003). Global RNA Half-Life Analysis in Escherichia coli Reveals Positional Patterns of Transcript Degradation. Genome Research, 216-223.</p>
+
<br>
-
<p>● iGEM CIDEB UANL. (2012). Parameters. Retrieved on June 1st, 2014. From: https://2013hs.igem.org/Team:CIDEB-UANL_Mexico/Math-Equations.</p>
+
\begin{equation}
-
<p>● Loera, Miguel. (2012). “Representación matemática de circuitos genéticos”.</p>
+
\large    \alpha_{1} =  \frac{transcription speed}{gene length (nt)}
-
<p>● PKU Beijing. (2009). Parameters. Retrieved on June 1st, 2014. From: https://2009.igem.org/Team:PKU_Beijing/Modeling/Parameters.</p>
+
\end{equation}<br>
-
<p>● Singapore. (2008). Parameters. Retrieved on June 1st, 2014. From:https://2008.igem.org/Team:NTU-Singapore/Modelling/Parameter.</p>
+
\begin{equation}
 +
\large    \alpha_{2} =  \frac{translation speed}{protein length (aa)}
 +
\end{equation}<br>
 +
\begin{equation}
 +
\large    \alpha_{1} =  \frac{70 \frac {nt}{s} \cdot(60)}{gene length (nt)}
 +
\end{equation}<br>
 +
\begin{equation}
 +
\large    \alpha_{2} =  \frac{40 \frac {aa}{s} \cdot(60)}{protein length (aa)}
 +
\end{equation}<br>
 +
\begin{equation}
 +
\large    d_{1} =  \frac{1}{half-life(min)} + \frac{1}{division time (min)}
 +
\end{equation}<br>
 +
\begin{equation}
 +
\large    d_{2} =  \frac{1}{half-life(min)} + \frac{1}{division time (min)}
 +
\end{equation}<br>
 +
\begin{equation}
 +
\large    d_{1} =  \frac{1}{half-life(min)} + \frac{1}{30min}
 +
\end{equation}<br>
 +
\begin{equation}
 +
\large    d_{2} =  \frac{1}{half-life(min)} + \frac{1}{30min}
 +
\end{equation}<br>
 +
 
 +
 
 +
<p><b>
 +
<h2>Bibliography/References</h2>
 +
</b></p>
 +
 
 +
<font size="2">
 +
<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>● 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>● 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>● 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>● 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>● 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>● 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>
<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