Team:CIDEB-UANL Mexico/Math-Improvement

From 2013hs.igem.org

(Difference between revisions)
 
(27 intermediate revisions not shown)
Line 54: Line 54:
</table>
</table>
-
<div style="background-color: #FFFFFF;">
+
 
-
<table width="100%" id="texto">
+
<table style="background-color: #FFFFFF;" width="100%" id="texto">
<tr>
<tr>
<td>
<td>
-
<p><b>Probability Equation</b></p>
+
<p><a href="https://2013hs.igem.org/Team:CIDEB-UANL_Mexico/Math-Improvement#Equation"><font color="blue">Probability Equation</font></a> - <a href="https://2013hs.igem.org/Team:CIDEB-UANL_Mexico/Math-Improvement#Production"><font color="blue">Concentration of Vip in <i>E.Coli</i></font></a><br>
-
<p>We tried to modify our model with a probability equation from an article of the National Institute for Mathematical and Biological Synthesis, Knoxville, Tennessee, United States of America published in the internet page PLOS ONE with the name <a href=" http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0011308"><font color="blue"> Is Thermosensing Property of RNA Thermometers Unique? </font></a>This work looks for the probability of finding a RNA at an openness way using its temperature, intercept and slope parameters time. This would describe the chance of finding our riboswitch k115017 open or closed. Making a smooth change in the fold of the thermosensor and the increase decrease slopes in our simulations.</p>
+
<a name="Equation"></a><b>Probability Equation</b> - <a href="https://2013hs.igem.org/Team:CIDEB-UANL_Mexico/Math-Improvement#"><font color="blue">Return</font></a></p></p>
 +
<p align="justify">We tried to modify our model with a probability equation from an article of the National Institute for Mathematical and Biological Synthesis, Knoxville, Tennessee, United States of America published on the internet page PLOS ONE with the name <a href=" http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0011308"><font color="blue"> Is Thermosensing Property of RNA Thermometers Unique? </font></a>This work looks for the probability of finding a RNA at an openness way using its temperature, intercept and slope parameters time. This would describe the chance of finding our riboswitch k115017 open or closed. Making a smooth change in the fold of the thermosensor and the increase, decrease slopes in our simulations.</p>
<br>
<br>
Line 66: Line 67:
\end{equation} </br>
\end{equation} </br>
-
<p>"P<sub>i</sub>(T) is the probability of finding the window at position “i” in a gene, open at temperature  (C), a<sub>i</sub> and b<sub>i</sub> are the intercept and slope parameters of how the log-odds of finding an open window at position “i”, log(p<sub>i</sub>(T)/1- p<sub>i</sub>(T)), changes with temperature. The ratio –a<sub>i</sub>/b<sub>i</sub>  indicates the temperature at which the probability of openness of a window is 0.5."
+
<p align="justify">"P<sub>i</sub>(T) is the probability of finding the window at position “i” in a gene, open at temperature  (C), a<sub>i</sub> and b<sub>i</sub> are the intercept and slope parameters of how the log-odds of finding an open window at position “i”, log(p<sub>i</sub>(T)/1- p<sub>i</sub>(T)), changes with temperature. The ratio –a<sub>i</sub>/b<sub>i</sub>  indicates the temperature at which the probability of openness of a window is 0.5."
</p>
</p>
-
<p>Our team used the last data of –a<sub>i</sub>/b<sub>i</sub> = T when openness is 0.5, to get “a<sub>i</sub>” expressed in b ( a<sub>i</sub>= -T*b<sub>i</sub>) and have our own goods ; to find b the authors of the articles,  compared across the 100 genes,  found that   values were not significantly greater than zero for 24 genes at p value is 0.05, so b value is 0.157. This implies that a small fraction of genes did not show a significant change in openness of its RBS with temperature over the range of temperatures considered.  Then we place our control temperature in the isolated equation and get a= -32*0.157 to get the “a” value that was negative 5.024.</p>
+
<p align="justify">Our team used the last data of –a<sub>i</sub>/b<sub>i</sub> = T when openness is 0.5, to get “a<sub>i</sub>” expressed in b ( a<sub>i</sub>= -T*b<sub>i</sub>) and have our own goods ; to find b the authors of the articles,  compared across 100 genes,  found that values were not significantly greater than zero for 24 genes at p value is 0.05, so b value is 0.157. This implies that a small fraction of genes did not show a significant change in openness of its RBS with temperature over the range of temperatures considered.  Then we place our control temperature in the isolated equation and get a= -32*0.157 to get the “a” value that was negative 5.024.</p>
-
<p>But there were some problems using the p<sub>i</sub>(T) formula in the matlab simulations like malfunction of the equation at the moment of seen graphs or the incorrect insertion of parameters in the simulation, so finally the team decided to leave it apart of the graphics. </p>
+
<p align="justify">But there were some problems using the p<sub>i</sub>(T) formula in the matlab simulations like malfunctioning of the equation when we tried seeing graphs, or the incorrect insertion of parameters in the simulation, so finally the team decided to leave it separate from the graphics. </p>
-
<p>We mention this special formula because it was the solution of the fast growth and decrease of the concentrations of the RNAs and proteins in E.coli. We run out of time and it was almost imposible to experiment with it and get better results, but heres the small advance we made with all this.  </p>
+
</td></tr></table>
-
+
-
<img "center"= src="https://static.igem.org/mediawiki/2013hs/9/95/Bi.png">
+
-
<p>The graph shows the probability of the increase in openness of an RNA thermometer. b was compared across the 100 genes, they found that values were not significantly greater than zero for 24 genes at p value= 0.05 and remember that ai and bi are the intercept and slope parameters of how the log-odds of finding an open window at position “i”, log(pi(T)/1- pi(T)), changes with temperature. So it can be conclude that a small fraction of genes did not show a significant change in openness of its RBS considering the range of temperature and this indicates that RNA thermometers do not have a great change from non-thermometers in increasing the openness of RBS with temperature. </p>
+
<table style="background-color: #FFFFFF;" width="100%" id="texto">
 +
<tr>
 +
<td>
 +
<p>We mention this special formula because it was the solution to the fast growth and decrease of the concentrations of the RNAs and proteins in <i>E.Coli</i>. We ran out of time and it was almost impossible to experiment with it and get better results, here's an advance of our job.   </p>
-
<h3><b>Concentration of Vip in E.Coli</h3></b>
+
<p align="justify">The graph shows the probability of the increase in openness of an RNA thermometer. b was compared across the 100 genes, they found that values were not significantly greater than zero for 24 genes at p value= 0.05 and remember that ai and bi are the intercept and slope parameters of how the log-odds of finding an open window at position “i”, log(pi(T)/1- pi(T)), changes with temperature. So it can be conclude that a small fraction of genes did not show a significant change in openness of its RBS considering the range of temperature and this indicates that RNA thermometers do not have a great change from non-thermometers in increasing the openness of RBS with temperature.  </p></td>
 +
<td style="padding:12px;"><center><img src="https://static.igem.org/mediawiki/2013hs/9/95/Bi.png"></center></td></tr></table>
 +
 
 +
<table style="background-color: #FFFFFF;" width="100%" id="texto">
 +
<tr>
 +
<td>
 +
<p align="justify"><b>Concentration of Vip in <i>E.Coli</i></b>
 +
 
 +
<br>In our variable of a stochastic model that we try to introduce it says that the overproduction of Vip3Ca3 will kill the host cell, in this case E.Coli, and in In the following article it supports our hypothesis about this situation, they investigated “the over production of seven membrane proteins in a E. Coli bacteriophage. In all seven cases, when the expression of the target membrane protein was induced, most of the host cell died, this make us think or suppose that the E. coli of our system can be affected by a similar case”. <a href=" http://www.ncbi.nlm.nih.gov/pubmed/8757792"><font color="blue">http://www.ncbi.nlm.nih.gov/pubmed/8757792</font></a>
 +
 
 +
<p align="justify"> In another article it states that “The difficulties in recombinant overexpression of membrane proteins arise partly from the physical space limitations of the cell membrane, but overexpression can also be toxic to the host cell”, also we can inferred that the overproduction can cause a toxicity that can affect after the host cell die and affect the system or the worm white.</p>
 +
<a href="http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2271156/#!po=16.6667"><font color="blue"> http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2271156/#!po=16.6667</font></a>
 +
 
 +
<p align="justify">Also we suggest that after the host cell die, we can´t know if the worm white dies will be for the production of Vip3Ca3 or the toxicity that it release after its death.</p>
 +
<p align="justify">We have to recall that the model it has to be in a control environment, in other words, in a lab with all the control of the variables in order that the entire model has a real function. The improvement it reaches a more reality to have other variables that are more specifically to an open environment.</p>
-
<p>Another variable that we have to consider is the saturation of Vip3Ca3 in E.coli in the environment that we want to put it on, in this case, is the gel that serves as a way of transportation. At the beginning, we assume that our variables are going to behave continuously and without noise, so it is a deterministic model but we consider that the overproduction in this medium is a variable of an stochastic model because it is a factor that affect our system and the production, in other words, is a noise that we have consider if we want to check the way that the system is going to work.</p>
 
-
<p>If the production of Vip3Ca3 oversaturates in the bacteria E.coli, then it is going to break down, the medium where our bacteria is will disappear and the production of Vip3Ca3 will stop. Also, another situation that we need to observe... if the white worm dies, the reasons could be the production of Vip3Ca3 or the aftermath of the rests of E.coli and other factors of the system. The implement of new stochastic equations to this model should answer these variables and make a more realistic modeling of our equation, but there isn’t a perfect model, you can only approach to it.</p>
 
-
</td>
 
-
<td style="background-color: #FFFFFF;">
 
-
<img src="https://static.igem.org/mediawiki/2013hs/d/d3/D9b5b18226b32476e3e16622aae4da9e.jpg" height="300px" />
 
</td>
</td>
</tr>
</tr>
</table>
</table>
-
</div>
 

Latest revision as of 02:29, 22 June 2013

Math Model
Improvement

Probability Equation - Concentration of Vip in E.Coli
Probability Equation - Return

We tried to modify our model with a probability equation from an article of the National Institute for Mathematical and Biological Synthesis, Knoxville, Tennessee, United States of America published on the internet page PLOS ONE with the name Is Thermosensing Property of RNA Thermometers Unique? This work looks for the probability of finding a RNA at an openness way using its temperature, intercept and slope parameters time. This would describe the chance of finding our riboswitch k115017 open or closed. Making a smooth change in the fold of the thermosensor and the increase, decrease slopes in our simulations.


\begin{equation} \large p_{i}\left ( T \right )= \frac{e^{a_{i}+b_{i}T}}{1+e^{a_{i}+b_{i}T}} \end{equation}

"Pi(T) is the probability of finding the window at position “i” in a gene, open at temperature (C), ai and bi are the intercept and slope parameters of how the log-odds of finding an open window at position “i”, log(pi(T)/1- pi(T)), changes with temperature. The ratio –ai/bi indicates the temperature at which the probability of openness of a window is 0.5."

Our team used the last data of –ai/bi = T when openness is 0.5, to get “ai” expressed in b ( ai= -T*bi) and have our own goods ; to find b the authors of the articles, compared across 100 genes, found that values were not significantly greater than zero for 24 genes at p value is 0.05, so b value is 0.157. This implies that a small fraction of genes did not show a significant change in openness of its RBS with temperature over the range of temperatures considered. Then we place our control temperature in the isolated equation and get a= -32*0.157 to get the “a” value that was negative 5.024.

But there were some problems using the pi(T) formula in the matlab simulations like malfunctioning of the equation when we tried seeing graphs, or the incorrect insertion of parameters in the simulation, so finally the team decided to leave it separate from the graphics.

We mention this special formula because it was the solution to the fast growth and decrease of the concentrations of the RNAs and proteins in E.Coli. We ran out of time and it was almost impossible to experiment with it and get better results, here's an advance of our job.

The graph shows the probability of the increase in openness of an RNA thermometer. b was compared across the 100 genes, they found that values were not significantly greater than zero for 24 genes at p value= 0.05 and remember that ai and bi are the intercept and slope parameters of how the log-odds of finding an open window at position “i”, log(pi(T)/1- pi(T)), changes with temperature. So it can be conclude that a small fraction of genes did not show a significant change in openness of its RBS considering the range of temperature and this indicates that RNA thermometers do not have a great change from non-thermometers in increasing the openness of RBS with temperature.

Concentration of Vip in E.Coli
In our variable of a stochastic model that we try to introduce it says that the overproduction of Vip3Ca3 will kill the host cell, in this case E.Coli, and in In the following article it supports our hypothesis about this situation, they investigated “the over production of seven membrane proteins in a E. Coli bacteriophage. In all seven cases, when the expression of the target membrane protein was induced, most of the host cell died, this make us think or suppose that the E. coli of our system can be affected by a similar case”. http://www.ncbi.nlm.nih.gov/pubmed/8757792

In another article it states that “The difficulties in recombinant overexpression of membrane proteins arise partly from the physical space limitations of the cell membrane, but overexpression can also be toxic to the host cell”, also we can inferred that the overproduction can cause a toxicity that can affect after the host cell die and affect the system or the worm white.

http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2271156/#!po=16.6667

Also we suggest that after the host cell die, we can´t know if the worm white dies will be for the production of Vip3Ca3 or the toxicity that it release after its death.

We have to recall that the model it has to be in a control environment, in other words, in a lab with all the control of the variables in order that the entire model has a real function. The improvement it reaches a more reality to have other variables that are more specifically to an open environment.

cideb
cideb
Contact us! Follow us on twitter and facebook or send us a mail.
CIDEB UANL Team. Centro de Investigación y Desarrollo de Educación Bilingüe
facebooktwitterenvelope