Team:Valencia/Modeling

From 2011.igem.org

(Difference between revisions)
 
(27 intermediate revisions not shown)
Line 1: Line 1:
 +
{{Template:Valencia-0-3}}
<html>
<html>
-
<body style="background-color:PowderBlue;">
 
<body>
<body>
-
</html>
+
<style type="text/css">
 +
.center {
-
{| style="color:#1b2c8a;background-color:powderblue;" cellpadding="3" cellspacing="1" border="1" bordercolor="#fff" width="85%" align="center"
+
width:980px;
-
!align="center"|[[Team:Valencia|Home]]
+
height:auto;
-
!align="center"|[[Team:Valencia/Team|Team]]
+
margin-left:px;
-
!align="center"|[https://igem.org/Team.cgi?year=2010&team_name=Valencia Official Team Profile]
+
        margin-right:px;
-
!align="center"|[[Team:Valencia/Project|Project]]
+
        margin-top:20px;
-
!align="center"|[[Team:Valencia/Parts|Parts Submitted to the Registry]]
+
margin-bottom:0px;
-
!align="center"|[[Team:Valencia/Modeling|Modeling]]
+
}
-
!align="center"|[[Team:Valencia/Notebook|Notebook]]
+
 
-
!align="center"|[[Team:Valencia/Safety|Safety]]
+
.col_center{
-
!align="center"|[[Team:Valencia/Attributions|Attributions]]
+
width:970px;
-
|}
+
        height:auto;
 +
float:left;
 +
background-color:transparent;
 +
 +
}
 +
 
 +
.col_center_top{
 +
width:908px;
 +
background-color:#0087ba;
 +
height:20px;
 +
font-family:Verdana, Arial, Helvetica, sans-serif;
 +
font-size:20px;
 +
color:#fff;
 +
padding:20px;
 +
        float:center;
 +
        margin-right:15px;
 +
 
 +
}
 +
 
 +
.col_center_center{
 +
width:800px;
 +
        text-align:justify;
 +
        height: auto;
 +
font-family:Verdana, Arial, Helvetica, sans-serif;
 +
font-size:14px;
 +
color:black;
 +
padding:20px;
 +
}
 +
 
 +
.col_center_bottom{
 +
width:900px;
 +
        text-align:center;
 +
        height: auto;
 +
font-family:Verdana, Arial, Helvetica, sans-serif;
 +
font-size:14px;
 +
color:#fff;
 +
padding:20px;
 +
}
 +
 
 +
 
 +
 
 +
 
 +
 
 +
 
 +
 
 +
html, body {
 +
margin: 0;
 +
padding: 0;
 +
}
 +
 
 +
body {
 +
background:darkblue url(https://static.igem.org/mediawiki/2011/6/65/Valencia_Aguamitadinferior1500x800.jpg) repeat-x top fixed;
 +
 
 +
}
 +
 
 +
#globalWrapper { background-color: transparent; border: none; margin: 0; padding: 0;}
 +
#content{
 +
                background:transparent url("https://static.igem.org/mediawiki/2011/f/fd/Valencia_Transpa30a.png") center fixed repeat-y;
 +
                padding: 0px;
 +
        margin:0 auto;
 +
                border: none; 
 +
           
 +
 +
            }
 +
 
 +
</style>
 +
<script type="text/javascript"
 +
src="http://ajax.googleapis.com/ajax/libs/jquery/1.3/jquery.min.js"></script>
 +
<script type="text/javascript">
 +
$(document).ready(function() {
 +
$("ul#topnav li").hover(function() {
 +
$(this).css({ 'background' : 'transparent url(https://static.igem.org/mediawiki/2010/0/03/Valencia_navbar_active.gif) repeat-x'});  
 +
$(this).find("span").show();
 +
} , function() { //on hover out...
 +
$(this).css({ 'background' : 'none'});
 +
$(this).find("span").hide();
 +
});
 +
 +
});
 +
</script>
 +
 
 +
<style type="text/css">
 +
 
 +
#headerlinks{
 +
      margin-left: -12px;
 +
}
 +
 
 +
#nav, #nav ul {
 +
padding: 0;
 +
list-style: none;
 +
}
 +
 
 +
.abstractext{
 +
float:right;
 +
width="400px"
 +
}
 +
 
 +
#nav a {
 +
display: block;
 +
        width: 149px;
 +
        height: 40px;
 +
        font-family:Verdana, Arial, Helvetica, sans-serif;
 +
}
 +
 
 +
#nav li {
 +
float: left;
 +
width: 149px;
 +
margin-right: 41px;
 +
}
 +
 
 +
#nav li ul {
 +
position: absolute;
 +
width: 130px;
 +
left: -999em;
 +
 +
}
 +
 
 +
#nav li:hover ul {
 +
left: auto;
 +
      text-decoration: none;
 +
      color: white !important;
 +
 
 +
}
 +
 
 +
#nav li:hover ul, #nav li.sfhover ul {
 +
left: auto;
 +
      text-decoration: none;
 +
      color: white !important;
 +
}
 +
 
 +
sfHover = function() {
 +
var sfEls = document.getElementById("nav").getElementsByTagName("LI");
 +
for (var i=0; i<sfEls.length; i++) {
 +
sfEls[i].onmouseover=function() {
 +
this.className+=" sfhover";
 +
}
 +
sfEls[i].onmouseout=function() {
 +
this.className=this.className.replace(new RegExp(" sfhover\\b"), "");
 +
}
 +
}
 +
}
 +
if (window.attachEvent) window.attachEvent("onload", sfHover);
 +
 +
body {
 +
font: 10px normal Verdana, Arial, Helvetica, sans-serif;
 +
 
 +
padding: 0;
 +
}
 +
h1 {
 +
font: 3em normal Verdana, 'Times New Roman', Times, serif;
 +
color: #333;
 +
margin: 0;
 +
padding: 5px 0;
 +
}
 +
h3 {
 +
font: 17px normal Verdana, 'Times New Roman', Times, serif;
 +
color: #0000;
 +
margin: 0;
 +
padding: 5px 0;
 +
margin: 10px;
 +
        text-align:justify;
 +
        margin-top: 8px;
 +
        margin-bottom: 20px;   
 +
        margin-left: 20px;
 +
        margin-right: 20px;
 +
}
 +
h1 small{
 +
font: 0.2em normal Verdana, Arial, Helvetica, sans-serif;
 +
text-transform:uppercase;
 +
letter-spacing: 1.5em;
 +
display: block;
 +
color: #fff;
 +
}
 +
.container {width: 970px; margin: 0 auto;}
 +
 
 +
.curved {
 +
        -moz-border-radius: 10px;
 +
        -webkit-border-radius: 10px;
 +
        behavior:url(border-radius.htc);
 +
}
 +
   
 +
 
 +
 
 +
 
 +
 
 +
ul#topnav {
 +
margin: 0; padding: 0;
 +
float: left;
 +
width: 970px;
 +
list-style: none;
 +
position: relative;
 +
font-size: 18.8px;
 +
background: url(topnav_stretch.gif) repeat-x;
 +
}
 +
ul#topnav li {
 +
float: left;
 +
margin: 0; padding: 0;
 +
border-right: 1px solid #555;
 +
}
 +
ul#topnav li a {
 +
padding: 10px 15px;
 +
display: block;
 +
color: darkblue;
 +
text-decoration: none;
 +
}
 +
li span {
 +
float: left;
 +
padding: 15px 0;
 +
position: absolute;
 +
left: 0; top:40px;
 +
display: none;
 +
width: 965px;
 +
background: lightblue;
 +
color: white;
 +
-moz-border-radius-bottomright: 5px;
 +
-khtml-border-radius-bottomright: 5px;
 +
-webkit-border-bottom-right-radius: 5px;
 +
-moz-border-radius-bottomleft: 5px;
 +
-khtml-border-radius-bottomleft: 5px;
 +
-webkit-border-bottom-left-radius: 5px;
 +
}
 +
ul#topnav li:hover span { display: block; }
 +
ul#topnav li span a { display: inline; }
 +
ul#topnav li span a:hover {text-decoration: underline;}
 +
</style>
 +
<script type="text/javascript"
 +
src="http://ajax.googleapis.com/ajax/libs/jquery/1.3/jquery.min.js"></script>
 +
<script type="text/javascript">
 +
$(document).ready(function() {
 +
 +
$("ul#topnav li").hover(function() { //Hover over event on list item
 +
$(this).css({ 'background' : '#1376c9 url(topnav_active.gif) repeat-x'}); //Add background color + image on hovered list item
 +
$(this).find("span").show(); //Show the subnav
 +
} , function() { //on hover out...
 +
$(this).css({ 'background' : 'none'}); //Ditch the background
 +
$(this).find("span").hide(); //Hide the subnav
 +
});
 +
 +
});
 +
</script>
 +
 
 +
 
 +
 
 +
 
 +
 
 +
<div class=center>
 +
 
 +
<div class=col_center>
 +
<div class=col_center_top><b>Modeling</b></div><!-- fin clase col_left_top-->
 +
<div class=col_center_center>
 +
 
 +
<h2>pH variations with light</h2>
 +
<p>We have successfully simulated the behaviour of pH changes in response to irradiation on cyanobacteria <i>Synechocystis</i> sp. PCC 6803, as you can see by this graph:
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/f/fd/VLC_Resultat.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/f/fd/VLC_Resultat.png" width="800" align="center" border="0" /></a></span></div></center></p>
 +
 
 +
<b>Why is this important?</b>
 +
<p>A model is not good or wrong, is more or less useful. Our model apart from making us think about the causes and relationships among the effects involved in this mechanism, has help us understand some results of the <a href="https://2011.igem.org/Team:Valencia/Project2" TARGET="_blank" title="pH variations">pH variations</a> in our pH-stat.</p>
 +
 
 +
<b>How have done this?</b>
 +
<p>First of all, we have to think on the physics of the behaviour we wanted to describe. We have light arriving upon a cyanobacteria and pH changing more or less proportionately to the time of irradiation.</p>
 +
<p>Thus, we have a light input and a pH variation output.</p>
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/9/90/VLC_Un_escalo.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/9/90/VLC_Un_escalo.png" width="800" align="center" border="0" /></a></span></div></center>
 +
<p>Now we need a transfer function that giving a light input delivers and output such as the one observed. To that purpose we can observe that the response (pH variation) behaves like a first order system with a lagging time.</p><br>
 +
<p>We work on the frequency plane (using Laplace transforms). Thus, transfer function of a first order system with a lagging time is such as:</p><br>
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/a/a6/VLC_Eq1.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/a/a6/VLC_Eq1.png" width="400" align="center" border="0" /></a></span></div></center>
 +
<p>Using the transfer function of the system we can get the response such as:</p><br>
 +
y(t) = pH - pH<sub>0</sub><br>
 +
y(s)= L [y(t)] where L[] function is Laplace transform<br>
 +
y(s) = G(s)·u(s)<br>
 +
<p>With the experimental data we have obtained the parameters values:</p><br>
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/9/93/VLC_Param.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/9/93/VLC_Param.png" width="400" align="" border="0" /></a></span></div></center><br>
 +
And with them we can adjust almost perfectly the experimental upwards dynamics:<br>
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/6/6c/VLC_Realisim.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/6/6c/VLC_Realisim.png" width="500" align="" border="0" /></a></span></div></center><br>
 +
We assume that the dynamics downwards should have the same underlying physics, so parameters should remain the same, we obtain a theoretical simulation and tallies real data:<br>
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/4/41/VLC_Baixada.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/4/41/VLC_Baixada.png" width="800" align="" border="0" /></a></span></div></center>
 +
As we can see, they are pretty much equivalent.<br>
 +
Anyway, let's consider now more than one peak. If we have different pulses of light, it's straight ofrward to think that we'll have different peaks of pH variations, but, if we look closer to the time series, we observe that gain response is weaker as time goes by and the basal value follows a slow, but constant upward tendency.<br>
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/4/47/VLC_Llarga.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/4/47/VLC_Llarga.png" width="700" align="" border="0" /></a></span></div></center><br>
 +
In order to model this behaviour we'll use two functions:<br>
 +
On one hand a buffering function that will multiply the response function
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/f/fe/VLC_At.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/f/fe/VLC_At.png" width="150" align="" border="0" /></a></span></div></center><br>
 +
 
 +
On the other hand a function that models appropriately the increase of a basal point is the following one and should be added to the response function.
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/4/48/VLC_Bt.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/4/48/VLC_Bt.png" width="200" align="" border="0" /></a></span></div></center><br>
 +
 
 +
As we have done before, we can calculate the following parameters:
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/3/31/VLC_Tresparam.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/3/31/VLC_Tresparam.png" width="300" align="" border="0" /></a></span></div></center><br>
 +
 
 +
With them, we can get the final equations:<br>
 +
y(s) = G(s)·u(s)<br>
 +
y(s)= L<sup>-1</sup>[y(s)]   (reversed Laplace transform)<br>
 +
Y(t) = y(t)·A(t)+B(t)<br>
 +
pH=Y(t) + pH<sub>0</sub><br><br>
 +
 
 +
Function y(t) is not easy to describe as it is a stepwise function depending of the time of the simulation, it is like:
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/1/13/VLC_Yt.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/1/13/VLC_Yt.png" width="875" align="" border="0" /></a></span></div></center>
 +
 
 +
Now, just run the set of equations and plot the results accordingly and you have the graph:
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/f/fd/VLC_Resultat.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/f/fd/VLC_Resultat.png" width="800" align="center" border="0" /></a></span></div></center>
 +
</br></br>
 +
<h2>Bacteriocin efficiency on killing bacteria</h2>
 +
<p>It is really important to understand the bacterias resistance against antimicrobial peptides for many bioremediation and health applications. With the following modeling, we want to know which is the response from our pathogens towards different bacteriocins concentrations and deduce which is the approximated amount of peptides needed to kill around the whole culture, which is 0.7µM of bacteriocin in our simulation.</p>
 +
<center><div><span><a href="https://static.igem.org/mediawiki/2011/3/34/Pepti.png" class="image" ><img alt="" src="https://static.igem.org/mediawiki/2011/3/34/Pepti.png" width="800" align="center" border="0" /></a></span></div></center>
 +
</br></br>
 +
The model above consists of a logistic regression with following expression:
 +
 
 +
<center><IMG SRC="https://static.igem.org/mediawiki/2011/1/14/Valeqmodel1.JPG" ALT="Overall Setup"></center>
 +
where we simulate the % of targeted killed bacteria vs the concentration of peptides.
 +
The fitting of the experimental values has an AIC of 49,27.<br/><br/>
 +
<center><IMG SRC="https://static.igem.org/mediawiki/2011/c/cd/Valeqmodel2.JPG" ALT="Overall Setup"></center>
 +
</br></br>
 +
Own experimental values where not validated, so similar results have been taken for the model from an external source: </br></br>
 +
"Clavanin permeabilizes target membranes via two distinct ph-dependent mechanisms", Ellen J. M. van Kan.

Latest revision as of 03:47, 22 September 2011



Modeling

pH variations with light

We have successfully simulated the behaviour of pH changes in response to irradiation on cyanobacteria Synechocystis sp. PCC 6803, as you can see by this graph:

Why is this important?

A model is not good or wrong, is more or less useful. Our model apart from making us think about the causes and relationships among the effects involved in this mechanism, has help us understand some results of the pH variations in our pH-stat.

How have done this?

First of all, we have to think on the physics of the behaviour we wanted to describe. We have light arriving upon a cyanobacteria and pH changing more or less proportionately to the time of irradiation.

Thus, we have a light input and a pH variation output.

Now we need a transfer function that giving a light input delivers and output such as the one observed. To that purpose we can observe that the response (pH variation) behaves like a first order system with a lagging time.


We work on the frequency plane (using Laplace transforms). Thus, transfer function of a first order system with a lagging time is such as:


Using the transfer function of the system we can get the response such as:


y(t) = pH - pH0
y(s)= L [y(t)] where L[] function is Laplace transform
y(s) = G(s)·u(s)

With the experimental data we have obtained the parameters values:



And with them we can adjust almost perfectly the experimental upwards dynamics:

We assume that the dynamics downwards should have the same underlying physics, so parameters should remain the same, we obtain a theoretical simulation and tallies real data:
As we can see, they are pretty much equivalent.
Anyway, let's consider now more than one peak. If we have different pulses of light, it's straight ofrward to think that we'll have different peaks of pH variations, but, if we look closer to the time series, we observe that gain response is weaker as time goes by and the basal value follows a slow, but constant upward tendency.

In order to model this behaviour we'll use two functions:
On one hand a buffering function that will multiply the response function

On the other hand a function that models appropriately the increase of a basal point is the following one and should be added to the response function.

As we have done before, we can calculate the following parameters:

With them, we can get the final equations:
y(s) = G(s)·u(s)
y(s)= L-1[y(s)] (reversed Laplace transform)
Y(t) = y(t)·A(t)+B(t)
pH=Y(t) + pH0

Function y(t) is not easy to describe as it is a stepwise function depending of the time of the simulation, it is like:
Now, just run the set of equations and plot the results accordingly and you have the graph:


Bacteriocin efficiency on killing bacteria

It is really important to understand the bacterias resistance against antimicrobial peptides for many bioremediation and health applications. With the following modeling, we want to know which is the response from our pathogens towards different bacteriocins concentrations and deduce which is the approximated amount of peptides needed to kill around the whole culture, which is 0.7µM of bacteriocin in our simulation.



The model above consists of a logistic regression with following expression:
Overall Setup
where we simulate the % of targeted killed bacteria vs the concentration of peptides. The fitting of the experimental values has an AIC of 49,27.

Overall Setup


Own experimental values where not validated, so similar results have been taken for the model from an external source:

"Clavanin permeabilizes target membranes via two distinct ph-dependent mechanisms", Ellen J. M. van Kan.