Team/CINVESTAV-IPN-UNAM MX/Differential.htm

From 2012.igem.org

(Difference between revisions)
(Created page with "Differential")
 
(6 intermediate revisions not shown)
Line 1: Line 1:
-
Differential
+
<html xmlns="http://www.w3.org/1999/xhtml">
 +
<head>
 +
<meta http-equiv="content-type" content="text/html; charset=utf-8" />
 +
<title>Rho</title>
 +
 
 +
<style type="text/css">
 +
 
 +
#top-section {
 +
position: relative;
 +
height: 20px;
 +
width: 975px;
 +
margin: 0 auto 0 auto;
 +
border-left: hidden;
 +
border-right: hidden;
 +
border-bottom: hidden;
 +
}
 +
 
 +
* {
 +
margin: 0;
 +
padding: 0;
 +
}
 +
 
 +
body {
 +
background: #FFFFFF url(https://static.igem.org/mediawiki/2012/f/f3/Fondo_bod.jpg) repeat-x;
 +
text-align: justify;
 +
font-family: Georgia, "Times New Roman", Times, serif;
 +
font-size: 12px;
 +
color: #000000;
 +
}
 +
 
 +
h1, h2, h3 {
 +
color: #333333;
 +
}
 +
 
 +
h1, h2 {
 +
margin-bottom: 20px;
 +
font-weight: normal;
 +
}
 +
 
 +
h1 {
 +
font-size: 197%;
 +
}
 +
 
 +
h2 {
 +
font-size: 152%;
 +
}
 +
 
 +
h3 {
 +
font-size: 100%;
 +
}
 +
 
 +
p, ul, ol {
 +
margin-bottom: 20px;
 +
line-height: 1.8em;
 +
}
 +
 
 +
ul, ol {
 +
margin-left: 3em;
 +
}
 +
 
 +
blockquote {
 +
padding-left: 3em;
 +
background: url(https://static.igem.org/mediawiki/2012/1/11/Img08.gif) no-repeat;
 +
font-style: italic;
 +
}
 +
 
 +
a:link {
 +
color: #1953A7;
 +
}
 +
 
 +
a:active {
 +
color: #ED5C0C;
 +
}
 +
 
 +
a:visited {
 +
color: #666666;
 +
}
 +
 
 +
a:hover {
 +
text-decoration: none;
 +
color: #067600;
 +
}
 +
 
 +
img {
 +
border: none;
 +
}
 +
 
 +
img.left {
 +
float: left;
 +
margin: 0 20px 10px 0;
 +
}
 +
 
 +
/* Menu */
 +
 
 +
#menu {
 +
width: 770px;
 +
height: 65px;
 +
margin: 0 auto;
 +
border-right: 1px solid #EDEBD5;
 +
}
 +
 
 +
#menu ul {
 +
margin: 0;
 +
padding: 0;
 +
list-style: none;
 +
}
 +
 
 +
#menu li {
 +
display: inline;
 +
font-size: 12px;
 +
color:#4E4E4E;
 +
}
 +
 
 +
#menu a {
 +
display: block;
 +
float: left;
 +
width: 95px;
 +
height: 45px;
 +
padding: 10px 0 0 0;
 +
border-left: 1px solid #EDEBD5;
 +
text-transform: lowercase;
 +
text-decoration: none;
 +
text-align: center;
 +
font-size: 135%;
 +
color: #BABABA;
 +
}
 +
 
 +
#menu a:hover {
 +
background: #F7F7F4;
 +
}
 +
 
 +
#menu .active a {
 +
background: #5dbb10;
 +
color: #FFFFFF;
 +
}
 +
 
 +
/* Header */
 +
 
 +
#header {
 +
width: 770px;
 +
height: 170px;
 +
margin: 0 auto;
 +
}
 +
 
 +
#header h1 {
 +
float: left;
 +
width: 220px;
 +
height: 100px;
 +
margin: 0;
 +
padding: 70px 0 0 0;
 +
background: url(https://static.igem.org/mediawiki/2012/0/0a/Img02.jpg) no-repeat;
 +
text-align: center;
 +
font-size: 36px;
 +
}
 +
 
 +
#header h2 {
 +
float: right;
 +
width: 530px;
 +
height: 86px;
 +
margin: 0;
 +
background: url(https://static.igem.org/mediawiki/2012/4/4a/Img03.jpg) no-repeat;
 +
padding: 84px 20px 0 0;
 +
text-align: right;
 +
}
 +
 
 +
#header a {
 +
text-transform: lowercase;
 +
text-decoration: none;
 +
font-style: italic;
 +
font-weight: normal;
 +
color: #FFFFFF;
 +
}
 +
 
 +
/* Page */
 +
 
 +
#page {
 +
width: 770px;
 +
margin: 0 auto;
 +
padding: 20px 0;
 +
}
 +
#context{
 +
float: left;
 +
width: 570px;
 +
padding-top:5px;
 +
}       
 +
/* menudecontenido */
 +
 
 +
#sidebar {
 +
float: right;
 +
width: 160px;
 +
padding: 15px;
 +
background: #EDEBD5;
 +
}
 +
 
 +
#sidebar h2 {
 +
margin: 0;
 +
padding: 5px 15px;
 +
color: #FFFFFF;
 +
}
 +
 
 +
#sidebar h3 {
 +
font-size: 60%;
 +
color: #FFFFFF;
 +
}
 +
 
 +
#sidebar ul {
 +
margin: 0;
 +
padding: 15px;
 +
list-style: none;
 +
}
 +
 
 +
#sidebar li {
 +
              font-size: 12px;
 +
  text-align:left;
 +
}
 +
 
 +
#sidebar a {
 +
color: #FFFFFF;
 +
}
 +
 
 +
.orangebox {
 +
background: url(https://static.igem.org/mediawiki/2012/b/b0/Img05.gif) repeat-y;
 +
}
 +
 
 +
.orangebox h2 {
 +
background: url(https://static.igem.org/mediawiki/2012/c/c2/Img06.gif) repeat-y;
 +
}
 +
 
 +
/* Footer */
 +
 
 +
#piep {
 +
clear: both;
 +
height: 150px;
 +
padding: 25px 0;
 +
background: #EDEBD5 url(https://static.igem.org/mediawiki/2012/a/a5/Img04.gif) repeat-x;
 +
}
 +
 
 +
#sponsors{
 +
 +
float: left;
 +
width: 85px;
 +
padding-top:5px;
 +
background: #EDEBD5;
 +
 +
      }
 +
#text2
 +
{
 +
color: #333333;
 +
font-size: 152%;
 +
margin-bottom: 20px;
 +
font-weight: normal;
 +
}
 +
#refe
 +
{
 +
color: #000000;
 +
font-size:9;
 +
margin-bottom: 20px;
 +
font-weight: normal;
 +
}
 +
.Estilo2 {color: #FDFDFD; font-weight: bold; }
 +
.Estilo3 {
 +
color: #FFFFFF;
 +
font-weight: bold;
 +
}
 +
 
 +
.firstHeading {
 +
margin-bottom: .1em;
 +
visibility:hidden;
 +
}
 +
 
 +
img[src*="/wiki/skins/common/images/wiki.png"]{
 +
    width:0;
 +
    height:0px;
 +
    padding-right:975px;
 +
    background:transparent url(https://2011.igem.org/wiki/skins/common/images/wiki.png) top left no-repeat;
 +
}
 +
 
 +
</style>
 +
 
 +
<meta name="keywords" content="" />
 +
<meta name="description" content="" />
 +
 
 +
</head>
 +
<body>
 +
<div id="menu">
 +
<ul>
 +
<li><a href="home.htm">Home</a></li>
 +
<li><a href="rhodofactory.htm">Rhodofactory</a></li>
 +
<li><a href="promoterselection.htm">Promoter selection </a></li>
 +
<li><a href="results.htm">results </a></li>
 +
<li><a href="perspectives.htm">Perspectives</a></li>
 +
<li class="active"><a href="Modelling.htm">Modeling</a></li>
 +
<li><a href="aboutus.htm">About us</a> </li>
 +
<li><a href="outreach.htm">outreach</a></li>
 +
</ul>
 +
</div>
 +
<!-- end #menu -->
 +
<div id="header">
 +
<h1><a href="#"></a></h1>
 +
<h2><span></span></h2>
 +
</div>
 +
<!-- end #header -->
 +
<div id="page">
 +
<div id="context">
 +
 +
<h1><em>Differential equations modeling the <br /><br />interactions between AppA and PpsR.!</em></h1>
 +
    <p>The proteins in our regulatory system constitute a network of intermingling elements
 +
whose behavior can be narrow into seven reactions. Primarily in this part of the modeling
 +
we focused on the PpsR and AppA anti-repressor/repressor system and chose to model it
 +
using a system of ordinary differential equations.The system’s repressive activity varies
 +
with the redox (O2) state of the cell and Light intensity and thus we wondered how exactly
 +
these two variables affect the overall repressive level of PpsR over the expression of PS
 +
genes.
 +
</p>
 +
<p id="text2">Model Construction.</p>
 +
 
 +
<p>To do this we based ourselves on a previous work that had attempted to study this same
 +
regulatory system. [Pandey et al, 2011]. Therefore the mathematical model we developed is
 +
based on the ODEs and kinetic parameters outlined there.<br><br>
 +
 
 +
 
 +
The following are its assumptions and basis:
 +
AppA inhibits the DNA-binding activity of oxidized PpsR by two mechanisms:<br>
 +
<br>
 +
 
 +
1. By reducing a disulfide bond in PpsR.<br>
 +
 
 +
2. By a blue-light-dependent sequestration of PpsR proteins into transcriptionally inactive
 +
complexes.<br>
 +
 
 +
At first stage, the reduced form of AppA (A-) reduces a disulfide bond in oxidized PpsR
 +
(P4+), which occurs independently of the light conditions.The molecular mechanism of
 +
this two-electron transfer is not yet clear.<br><br>
 +
 
 +
 
 +
Some experiments have shown that both PpsR and AppA have two redox-active Thiol
 +
groups that can form intramolecular disulfide bonds with a similar midpoint potential,
 +
according to this equation:</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/4/4b/Dife01.png" width="301" height="76"><br>
 +
 
 +
</div>
 +
<p><br>
 +
Then at the second level of regulation, the reduced form of AppA can form a complex with
 +
reduced PpsR. During the complex formation, one AppA molecule is associated with two
 +
PpsR monomers corresponding to half of a PpsR molecule, which exists as a stable
 +
tetramer in solution.<br>
 +
<br>
 +
 
 +
However there is a point were complex formation is inhibited by high intensities of blue-
 +
light irradiation (LI ¼ 900mmol/m2s). A subsequent study found that AppA responds to
 +
blue light over several orders of magnitude down to 0.2 mmol/m2s.</p>
 +
<p>Other experiments indicate that light absorption induces a structural change in the BLUF
 +
domain of AppA, which results in interactions with its C-terminal part, thereby causing the
 +
dissociation of PpsR.</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/3/33/Dife02.png"><br>
 +
 
 +
</div>
 +
<p><br>
 +
To implement the redox-sensing capabilities of AppA, we use the model proposed by Han
 +
et al, 2007 according to which AppA utilizes heme as a cofactor, bound to its C-terminal
 +
domain, to sense the cytosolic redox conditions, according to this equation:</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/6/6f/Dife03.png">
 +
</div>
 +
<p><br>
 +
  If the electron transfer from AppA to PpsR in Eq. 1 was indeed effectively irreversible:
 +
(kpr <sup>-</sup> &lt;&lt; kpr <sup>+</sup>), as suggested by the experiments of Masuda and Bauer, 2002.PpsR
 +
would have to be reoxidized through an AppA-independent mechanism. To account for
 +
this possibility, the assumption is that PpsR is reoxidized proportional to the oxygen
 +
concentration as:</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/5/5c/Dife04.png">
 +
</div>
 +
<p><br>
 +
Finally we assumed mass-action kinetics for the reactions above and the following are the
 +
set of ordinary differential equations established:</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/f/fc/Dife05.png" width="565" height="313"></div>
 +
<p><br>
 +
Were the total amounts of PpsR and AppA molecules are conserved according to:</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/d/d7/Dife06.png"></div>
 +
  <P><br>
 +
    Later on we reduced this set of ODEs into seven basic reactions, to simplify the code
 +
written in Mathematica.</P>
 +
<p>This is the code written in Mathematica.</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/e/ef/Dife07.png">
 +
</div>
 +
<p><br>
 +
In addition to contacting the responsible authors and discussing this model, we chose to
 +
improve it based on recent publications and preliminary experimental results i.e. we
 +
introduced a more subtle way light and oxygen affects the protein concentration dynamic,
 +
which is exactly what we were interested in understanding better.</p>
 +
<p id="text2">Differential equations Results:</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/d/d9/Grafidif04.jpg" width="400" height="291">
 +
<img src="https://static.igem.org/mediawiki/2012/d/dc/Grafidif02.jpg" width="414" height="321">
 +
<img src="https://static.igem.org/mediawiki/2012/6/63/Grafidif03.jpg" width="414" height="321">
 +
<img src="https://static.igem.org/mediawiki/2012/d/db/Grafidif01.jpg" width="420" height="307">
 +
<p>Figure1. Steady state behavior of reduced PpsR (P4+), oxidized PpsR (P4-) and the
 +
AppA-PpsR complex as a function of oxygen concentration and light irradiance.</p>
 +
</div>
 +
<p>If we assumed that the reduction of PpsR by AppA in Equation 1.</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/4/4b/Dife01.png" width="301" height="76"><br></div>
 +
<p><br>
 +
Is effectively irreversible when <em><strong>Keq>>1</strong></em></p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/f/f1/Dife09.png" width="100" height="66"></div>
 +
<p><br>
 +
  We obtain a fifth order polynomial equation that admits at most five real roots,
 +
corresponding to five possible stationary states. (Graphics above).</p>
 +
<p>For them to be biologically meaningful we require that they fall within the interval [0,1].
 +
Using our differential equations system we calculated protein concentrations at stable
 +
steady states and assumed this to be a natural “resting”, or homeostatic, state given a set
 +
of parameters, oxygen and light intensities. As such, this gave us a steady state total
 +
repressive force, in essence, as a function of oxygen tension and light intensity. This is
 +
because PpsR/AppA’s total repressive strength is a function of the concentrations of both
 +
its oxidative and reduced state.<br>
 +
<br>
 +
 
 +
Since this was all done in Mathematica, we were also able to generate an interactive
 +
graphical representation of the repressive strength over a range of conditions.</p>
 +
<p>Download the code:<a href="https://static.igem.org/mediawiki/2012/8/8b/Interactive_Graphical_AP.zip" target="_parent">Interactive Graphical AP</a></p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2012/6/64/Dife10.png" width="561" height="282">
 +
  <p>This is how it will look like when you run the program.</p></div>
 +
<p><br>
 +
  With this model, we can make predictions in sillico as to how the organism will react to a
 +
given environmental surrounding and then go back to the lab and validate that prediction.
 +
Also, since we have are able to manipulate the values of the parameters in an easy and
 +
intuitive way, we can explore the parameter space and visually observe the change in the
 +
concentration curves, which can be of interest to, not only us, but other people that might
 +
not be interested in all the math behind interphase.<br>
 +
<br>
 +
 
 +
 
 +
This parameter tweaking allows us to compare experimental data against sillico
 +
predictions; the only difference is that if we focus on the parameters, this can actually give
 +
us hints respect to the still-blurred nature of Protein-Protein and Protein-Environment
 +
mechanisms that exist in this extremely versatile cell. </p>
 +
 
 +
</div>
 +
<div id="sidebar">
 +
  <div id="updates" class="orangebox">
 +
 +
    <h2>Modeling</h2>
 +
<ul>
 +
<li><a href="Tools.htm" target="_parent">Tools we used</a></li>
 +
<li><a href="Differential.htm">Differential equations modeling</a></li>
 +
<li><a href="Rule-Based.htm">Rule-Based model</a></li>
 +
<li><a href="Cellular.htm">Cellular Automata</a></li>
 +
</ul>
 +
</div>
 +
  </div>  
 +
  <div style="clear: both;">&nbsp;</div>
 +
</div>
 +
<!-- end #page -->
 +
<div id="piep">
 +
<p align="center"> <strong>Rhodofactory 2012</strong> </p>
 +
<div id="sponsors">
 +
  <div align="center"><img src="https://static.igem.org/mediawiki/2012/8/8a/Icytdf.png" alt="icytdf" width="90" height="82" /></div>
 +
</div>
 +
<div id="sponsors">
 +
  <div align="center"><img src="https://static.igem.org/mediawiki/2012/6/67/Osli.png" alt="osli" width="88" height="46" /></div>
 +
</div>
 +
 +
<div id="sponsors">
 +
  <div align="center"><img src="https://static.igem.org/mediawiki/2012/f/fb/Bio.png" alt="bio" width="80" height="97" /></div>
 +
</div>
 +
<div id="sponsors"><img src="https://static.igem.org/mediawiki/2012/d/d6/Fermentas.png" alt="fermentAS" width="82" height="45" /></div>
 +
<div id="sponsors">
 +
  <div align="center"><img src="https://static.igem.org/mediawiki/2012/a/a3/Cinestav.png" alt="cinestav" width="74" height="91" /></div>
 +
</div>
 +
<div id="sponsors"><img src="https://static.igem.org/mediawiki/2012/9/9b/Genscript.png" alt="genscript" width="83" height="45" /></div>
 +
<div id="sponsors">
 +
  <div align="center"><img src="https://static.igem.org/mediawiki/2012/1/16/Unam.png" alt="unam" width="83" height="93" /></div>
 +
</div>
 +
<div id="sponsors"><img src="https://static.igem.org/mediawiki/2012/c/c9/Gto.png" alt="gto" width="84" height="46" /></div>
 +
    <div id="sponsors">
 +
      <div align="center"><img src="https://static.igem.org/mediawiki/2012/f/fa/Qimica.png" alt="quimica" width="80" height="72" /></div>
 +
    </div>
 +
<div id="sponsors"><img src="https://static.igem.org/mediawiki/2012/7/73/Valaner.png" alt="valaner" width="84" height="48" /></div>
 +
<div id="sponsors">
 +
  <div align="center"><img src="https://static.igem.org/mediawiki/2012/7/79/Ipn.png" alt="ipn" width="70" height="96" /></div>
 +
</div>
 +
 +
</div>
 +
<!-- end #footer -->
 +
</body>
 +
</html>

Latest revision as of 02:05, 27 October 2012

Rho

Differential equations modeling the

interactions between AppA and PpsR.!

The proteins in our regulatory system constitute a network of intermingling elements whose behavior can be narrow into seven reactions. Primarily in this part of the modeling we focused on the PpsR and AppA anti-repressor/repressor system and chose to model it using a system of ordinary differential equations.The system’s repressive activity varies with the redox (O2) state of the cell and Light intensity and thus we wondered how exactly these two variables affect the overall repressive level of PpsR over the expression of PS genes.

Model Construction.

To do this we based ourselves on a previous work that had attempted to study this same regulatory system. [Pandey et al, 2011]. Therefore the mathematical model we developed is based on the ODEs and kinetic parameters outlined there.

The following are its assumptions and basis: AppA inhibits the DNA-binding activity of oxidized PpsR by two mechanisms:

1. By reducing a disulfide bond in PpsR.
2. By a blue-light-dependent sequestration of PpsR proteins into transcriptionally inactive complexes.
At first stage, the reduced form of AppA (A-) reduces a disulfide bond in oxidized PpsR (P4+), which occurs independently of the light conditions.The molecular mechanism of this two-electron transfer is not yet clear.

Some experiments have shown that both PpsR and AppA have two redox-active Thiol groups that can form intramolecular disulfide bonds with a similar midpoint potential, according to this equation:



Then at the second level of regulation, the reduced form of AppA can form a complex with reduced PpsR. During the complex formation, one AppA molecule is associated with two PpsR monomers corresponding to half of a PpsR molecule, which exists as a stable tetramer in solution.

However there is a point were complex formation is inhibited by high intensities of blue- light irradiation (LI ¼ 900mmol/m2s). A subsequent study found that AppA responds to blue light over several orders of magnitude down to 0.2 mmol/m2s.

Other experiments indicate that light absorption induces a structural change in the BLUF domain of AppA, which results in interactions with its C-terminal part, thereby causing the dissociation of PpsR.



To implement the redox-sensing capabilities of AppA, we use the model proposed by Han et al, 2007 according to which AppA utilizes heme as a cofactor, bound to its C-terminal domain, to sense the cytosolic redox conditions, according to this equation:


If the electron transfer from AppA to PpsR in Eq. 1 was indeed effectively irreversible: (kpr - << kpr +), as suggested by the experiments of Masuda and Bauer, 2002.PpsR would have to be reoxidized through an AppA-independent mechanism. To account for this possibility, the assumption is that PpsR is reoxidized proportional to the oxygen concentration as:


Finally we assumed mass-action kinetics for the reactions above and the following are the set of ordinary differential equations established:


Were the total amounts of PpsR and AppA molecules are conserved according to:


Later on we reduced this set of ODEs into seven basic reactions, to simplify the code written in Mathematica.

This is the code written in Mathematica.


In addition to contacting the responsible authors and discussing this model, we chose to improve it based on recent publications and preliminary experimental results i.e. we introduced a more subtle way light and oxygen affects the protein concentration dynamic, which is exactly what we were interested in understanding better.

Differential equations Results:

Figure1. Steady state behavior of reduced PpsR (P4+), oxidized PpsR (P4-) and the AppA-PpsR complex as a function of oxygen concentration and light irradiance.

If we assumed that the reduction of PpsR by AppA in Equation 1.



Is effectively irreversible when Keq>>1


We obtain a fifth order polynomial equation that admits at most five real roots, corresponding to five possible stationary states. (Graphics above).

For them to be biologically meaningful we require that they fall within the interval [0,1]. Using our differential equations system we calculated protein concentrations at stable steady states and assumed this to be a natural “resting”, or homeostatic, state given a set of parameters, oxygen and light intensities. As such, this gave us a steady state total repressive force, in essence, as a function of oxygen tension and light intensity. This is because PpsR/AppA’s total repressive strength is a function of the concentrations of both its oxidative and reduced state.

Since this was all done in Mathematica, we were also able to generate an interactive graphical representation of the repressive strength over a range of conditions.

Download the code:Interactive Graphical AP

This is how it will look like when you run the program.


With this model, we can make predictions in sillico as to how the organism will react to a given environmental surrounding and then go back to the lab and validate that prediction. Also, since we have are able to manipulate the values of the parameters in an easy and intuitive way, we can explore the parameter space and visually observe the change in the concentration curves, which can be of interest to, not only us, but other people that might not be interested in all the math behind interphase.

This parameter tweaking allows us to compare experimental data against sillico predictions; the only difference is that if we focus on the parameters, this can actually give us hints respect to the still-blurred nature of Protein-Protein and Protein-Environment mechanisms that exist in this extremely versatile cell.

 

Rhodofactory 2012

icytdf
osli
bio
fermentAS
cinestav
genscript
unam
gto
quimica
valaner
ipn