Team:Valencia/Modeling
From 2011.igem.org
Line 267: | Line 267: | ||
<h2>pH variations with light</h2> | <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> | <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> | ||
- | 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.<br> | + | <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> |
- | We work on the frequency plane (using Laplace transforms). Thus, transfer function of a first order system with a lagging time is such as:<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> | <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> | ||
- | Using the transfer function of the system we can get the response such as:<br> | + | <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(t) = pH - pH<sub>0</sub><br> | ||
y(s)= L [y(t)] where L[] function is Laplace transform<br> | y(s)= L [y(t)] where L[] function is Laplace transform<br> | ||
y(s) = G(s)·u(s)<br> | y(s) = G(s)·u(s)<br> | ||
- | With the experimental data we have obtained the parameters values:<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> | + | <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: | + | 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> | + | <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: | + | 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> | <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. | + | 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. | + | 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: | + | In order to model this behaviour we'll use two functions:<br> |
On one hand a buffering function that will multiply the response function | 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. | 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: | 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: | + | With them, we can get the final equations:<br> |
y(s) = G(s)·u(s)<br> | y(s) = G(s)·u(s)<br> | ||
- | y(s)= L<sup>-1<sup>[y(s)] (reversed Laplace transform)<br> | + | y(s)= L<sup>-1</sup>[y(s)] (reversed Laplace transform)<br> |
Y(t) = y(t)·A(t)+B(t)<br> | Y(t) = y(t)·A(t)+B(t)<br> | ||
- | pH=Y(t) + pH<sub>0</sub><br> | + | pH=Y(t) + pH<sub>0</sub><br><br> |
- | Function y(t) is not easy | + | 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> | |
<h2>Bacteriocin efficiency on killing bacteria</h2> | <h2>Bacteriocin efficiency on killing bacteria</h2> |
Revision as of 02:00, 22 September 2011
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:
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:
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: