Team:Grenoble-EMSE-LSU/Project/Modelling

From 2013.igem.org

(Difference between revisions)
 
(48 intermediate revisions not shown)
Line 27: Line 27:
<li id="titre">
<li id="titre">
<h1>Modelling</h1>
<h1>Modelling</h1>
-
<p>Modeling took a large place in the project; it was not only used for the characterization of KillerRed and the Voigt plasmids, it was needed for the control of the bacteria’s population. With our device, we cannot control a population with a simple closed-loop transfer function. First, because it doesn’t measure the amount of cells (the OD 600) but the fluorescence of the medium. This fluorescence gives clues about the activity of the cells and so the look of its evolution permits to find their number, but there is no simple relation between them. Second, there is a big lag between an action and its effect: one can count one or two hours between the illumination and the deceleration of fluorescence. In those conditions, a simple closed-loop transfer function is bound to be unstable, a model predictive control is needed to stabilize the population.</p>
+
<p>Modelling took a large place in the project. We used modelling for the characterization of KillerRed and the Voigt plasmids and also needed it for the control of the bacterial population. With our device, we cannot control a population of living cells with a simple closed-loop transfer function. The first reason is that optical measurements (OD at 600 nm or fluorescence) originate from all cells, whether they are alive or not. As a consequence, it is not easy to reconstruct the size of a population of living cells from a fluorescent intensity or an $OD_{600}$ reading. Second, there is a large delay between an action and its effect: there are about one or two hours between the onset of illumination and a decrease in fluorescence growth rate. In those conditions, a simple closed-loop transfer function is predictably unstable, and a predictive model is needed to stabilize the population of living cells.</p>
-
</li>
+
                                </li>
-
+
-
<li>
+
-
<h1>Building the Model</h1>
+
-
<h2>Initial Model</h2>
+
-
<h3>The equation</h3>
 
-
<p>    Our system is made of Bacteria and ‘KillerRed’ proteins. A bacteria divides and produces KillerRed, and a KillerRed reacts to light: it degrades (photobleaching) and produces Radical Oxygen Species or ROS (phototoxicity). These reactions are shared by all fluorescent proteins, but the shape of KillerRed makes the degradation quicker and permits ROS to wander in the bacteria and damage its vital functions (proteins, cell wall or worse: DNA)</p>
 
-
<br>
 
-
<p>$C$ the amount of living bacteria </p>
 
-
<p>$K$ the amount of KillerRed inside the living bacteria</p>
 
-
<p>$I$ the quantity of received light</p>
 
 +
                                <li>
 +
<a href="/Team:Grenoble-EMSE-LSU/Project/Modelling/Building" title="BuildingLink"><h3>Building the Model</h3></a>
 +
<a href="/Team:Grenoble-EMSE-LSU/Project/Modelling/Building" title="BuildingLink2">
 +
<img src="https://static.igem.org/mediawiki/2013/5/5b/Gre_Mod_Lego.png" style="float:left"></a>
 +
<p style="height:100px;float:none;padding-left:125px;padding-right:150px"> Come, sit, and listen to the story of the model's construction! You will hear about the journey of this model, from its genesis to its completion ! Going through the reasons that drove it to consider the time of <a href="https://2013.igem.org/Team:Grenoble-EMSE-LSU/Project/Modelling/Building#MatTime">maturation of KillerRed</a> or the <a href="https://2013.igem.org/Team:Grenoble-EMSE-LSU/Project/Modelling/Building#AccDam">resilience of bacteria</a>.</p>
-
<center style="font-size:150%;">
 
-
 
-
$
 
-
\left\{
 
-
  \begin{array}{l l}
 
-
    \frac{dC}{dt}=rC-kIK \\
 
-
    \frac{dK}{dt}=aC-bIK-kI\frac{K^2}{C} \\
 
-
  \end{array}
 
-
\right.
 
-
$
 
-
</center>
 
-
 
-
<p> $rC$ is the term of growth.</p>
 
-
<p> $kIK=kI\frac{K}{C}C$ the amount of bacteria killed by KillerRed and light.</p>
 
-
<p> $aC$ the production of KillerRed.</p>
 
-
<p> $bIK$ the amount of KillerRed photobleached.</p>
 
-
<p> $kIK\frac{K}{C}$ the amount of KillerRed in the bacteria killed in the last step of time.</p>
 
<br>
<br>
-
<p>    Unfortunately, $C$ and $K$ are not measurable variables. The only thing we can quickly and easily measure are the optical density (OD) bound to the amount of bacteria dead AND alive, and the global fluoresence bound to the amount of KillerRed in the bacteria dead AND alive. To be able to compare our model to experimental results, we need two other variables : </p>
 
<br>
<br>
-
<p>$D$ the amount of dead bacteria </p>
 
-
<p>$K_d$ the amount of KillerRed inside the dead bacteria</p>
 
-
<center style="font-size:150%;">
+
<a href="/Team:Grenoble-EMSE-LSU/Project/Modelling/Parameters" title="ParamLink"><h3>Finding Parameters</h3></a>
-
 
+
<a href="/Team:Grenoble-EMSE-LSU/Project/Modelling/Parameters" title="ParamLink2">
-
$
+
<img src="https://static.igem.org/mediawiki/2013/d/d4/Gre_Mod_Tron.png" style="float:left"></a>
-
\left\{
+
<p style="height:100px;float:none;padding-left:125px;padding-right:150px"> Here, you traveler, you will read about the way the parameters were chosen to best fit the experiments, and thus make the model properly predict the evolution of the bacterial concentration. You will see how a <a href="https://2013.igem.org/Team:Grenoble-EMSE-LSU/Project/Modelling/Parameters#AlgoGen">genetic algorithm</a> can be used to understand a genetic network. And of course you will appreciate the <a href="https://2013.igem.org/Team:Grenoble-EMSE-LSU/Project/Modelling/Parameters#Results">final results</a>.</p>
-
  \begin{array}{l l}
+
-
    \frac{dD}{dt}=kIK \\
+
-
    \frac{dK_d}{dt}=kI\frac{K^2}{C}-bIK_d\\
+
-
  \end{array}
+
-
\right.
+
-
$
+
-
</center>
+
<br>
<br>
-
<br>
 
-
<p> The simplest units possible were used, directly taken from our measure instruments : </p>
 
-
<p> $C$ and $D$ are in 'OD'.</p>
 
-
<p> $K$ and $K_d$ are in 'units of fluorescence' : UF. The auto-fluorescence of bacteria are considered as negligeable compared to the fluorescence of KillerRed. </p>
 
<br>
<br>
-
 
-
<h3>Analytic resolution</h3>
 
-
 
-
<p> This simple model can be partially solved, for $C(t)$ or $I(t)$ constant for example : </p>
 
-
<br>
 
-
<p>If we say : $\forall t, C(t)=C_0$, we have :</p>
 
-
<p>
 
-
$
 
-
\left\{
 
-
  \begin{array}{l l}
 
-
    rC_0=kIK \\
 
-
    \frac{dK}{dt}=aC_0-bIK-kI\frac{K^2}{C_0}\\
 
-
  \end{array}
 
-
\right.
 
-
$
 
-
</p>
 
-
<br>
 
-
<p> and so : $\frac{dK}{dt}=\left(a-\frac{br}{k}\right)C_0-rK$</p>
 
-
<p> which gives : $K=\left(\frac{a}{r}-\frac{b}{k}\right)C_0+Be^{-rt}$</p>
 
-
<br>
 
-
<p> then $I(t)=\frac{rC_0}{kK(t)}$ should give a constant concentration of living cell.</p>
 
-
<br>
 
-
<br>
 
-
<p> But if we say : $\forall t, I(t)=I_0$, we need another variable to solve easily our equation : </p>
 
-
<p> We define : $Y=\frac{K}{C}$ the amount of KillerRed per bacteria.</p>
 
-
<p> $\frac{dY}{dt}=\frac{d}{dt}\left(\frac{K}{C}\right)=a-(bI_0+r)Y$</p>
 
-
<p> which gives : $Y=\frac{a}{bI_0+r}+Be{-(bI_0+r)t} $</p>
 
-
<br>
 
-
<p> $Y$ is then stabilizing, let's look $C$ :</p>
 
-
<p>$\frac{dC}{dt}=C(r-kI_0Y)$</p>
 
-
<p>$\frac{dC}{dt}=C\left(r-\frac{kI_0a}{bI_0+r}+Be^{-(bI_0+r)t}\right)$</p>
 
-
<p> and so, for $I_0=\frac{r^2}{ak-rb}$, we have : $\lim_{t\to\infty}\frac{dC}{dt}(t)=0$</p>
 
-
<br>
 
-
<p> The resolution of this equation have shown the possibility to stabilize the system.</p> 
 
-
 
-
 
-
<h3>Comparison with experiments</h3>
 
-
 
-
<p> This first model is noteworthy for solutions and conditions on a stabilized state appear. But it is not complete enough to explain the experiments : </p>
 
-
 
-
 
-
 
 +
<a href="/Team:Grenoble-EMSE-LSU/Project/Modelling/Density" title="ParamLink"><h3>Density Control</h3></a>
 +
<a href="/Team:Grenoble-EMSE-LSU/Project/Modelling/Density" title="ControlLink">
 +
<img src="https://static.igem.org/mediawiki/2013/d/d7/Gre_Mod_Manette.png" style="float:left"></a>
 +
<p style="height:100px;float:none;padding-left:125px;padding-right:150px""> At last, you will discover <a href="https://2013.igem.org/Team:Grenoble-EMSE-LSU/Project/Modelling/Density#RbL">the power and the impact</a> that light has on our system. Then you will be introduced to a <a href="https://2013.igem.org/Team:Grenoble-EMSE-LSU/Project/Modelling/Density#MPC">crafty way</a> to master it and control a bacterial population with its help. You will be shown <a href="https://2013.igem.org/Team:Grenoble-EMSE-LSU/Project/Modelling/Density#MCBP">how it is supposed to be done</a>.</p>
<br>
<br>
<br>
<br>
-
<p> Whereas we observe a lag between the light switch and the decrease of fluorescence, the model predicts an immediate and clear reaction. A new way to explain and describe the lag between the stimulus (the light) and the reaction (the decrease of fluorescence and the sabilization of DO) has to be found. Of course this way shall have a biological explanation. </p>
 
-
      </li>
 
-
 
-
<li>
 
-
<h2>Maturation Time</h2>
 
-
<h3>The maturation of fluorescent proteins</h3>
 
-
<p>After the traduction of mRNA into a polypeptide, then its correct foldind, a fluorescent protein still has one step to do before ot works properly : the maturation. The protein is matured after an oxydation reaction in which it gets its chromophore. The time of a maturation is significant for our experiment, for a GFP, it is 30 minutes. </p>
 
-
<h3>Second model </h3>
 
-
<p> We consider maturation to be a simple chemical reaction, and the concentration on immature protein $[K_im]$ is governed by the equation : $[K_im]'=-m[K_im]$. It is a realistic hypothesis since maturation is an oxydation. A new variable is needed : </p>
 
-
 
-
<br>
 
-
<p>$K_m$ the amount of mature KillerRed inside the living bacteria</p>
 
-
<p>$K_i$ the amount of immature KillerRed inside the living bacteria</p>
 
-
 
-
<center style="font-size:150%;">
 
-
$
 
-
\left\{
 
-
  \begin{array}{l l}
 
-
    \frac{dC}{dt}=rC-kIK_m \\
 
-
    \frac{dK_i}{dt}=aC-kI\frac{K_i^2}{C}-mK_i\\
 
-
    \frac{dK_m}{dt}=-kI\frac{K_m^2}{C}-bIK_m+mK_i\\
 
-
  \end{array}
 
-
\right.
 
-
$
 
-
</center>
 
-
 
-
 
-
<p> $mK_i$ is the term traducing the maturation of KillerRed </p>
 
-
 
-
<br>
 
-
<p> Here too, other equations are required to find the measureable variables : </p>
 
-
 
-
<center style="font-size:150%;">
 
-
$
 
-
\left\{
 
-
  \begin{array}{l l}
 
-
    \frac{dD}{dt}=kIK_m \\
 
-
    \frac{dK_{di}}{dt}=kI\frac{K_i^2}{C}-mK_{di}\\
 
-
    \frac{dK_{dm}}{dt}=kI\frac{K_m^2}{C}-bIK_{dm}+mK_{di}\\
 
-
  \end{array}
 
-
\right.
 
-
$
 
-
</center>
 
-
 
-
<br>
 
-
<h3> Numeric results</h3>
 
-
<p> This model draws curves that properly follow the tendancy of the experiments : the lag of the reaction, the peak of fluorescence and then the brutal decrease are there. </p>
 
-
<br>
 
-
<br>
 
-
<p>Nonetheless, it is impossible to get a good fit between the prediction of the model and the experiment. The integration of maturation doesn't explain why the production of KillerRed is so hight even two hours after the start of illumination and the decrease of fluorescence is so brutal four hours after the illumination.</p>
 
-
 
</li>
</li>
-
 
-
<li>
 
-
<h2>Damage Accumulation</h2>
 
-
<p> Until now, in the equations, KillerRed and light produced ROS, those reacted in the instant, killed or not the bacteria, and disappeared. Now, we consider that bacteria are unable to instantly repear all the damages caused by ROS, and with damages mounting up, they are more and more fragile and close to death. Considering this accumulation allows us to shift the effect of illumination, and so to have a model more accurate.</p>
 
-
<br>
 
-
<p> Writen in its numeric form, the evolution of C was : </p>
 
-
<center style="font-size:150%;">
 
-
$
 
-
C(t+1)-C(t)=rC(t)-kI(t)K(t)
 
-
$
 
-
</center>
 
-
<p> It is now : </p>
 
-
<center style="font-size:150%;">
 
-
$
 
-
\left\{
 
-
  \begin{array}{l l}
 
-
    C(t+1)-C(t)=rC(t)-\mbox{tox}(t) \\
 
-
    \mbox{tox}(t+1)=l.\mbox{tox}(t)+k'I(t)K(t)\\
 
-
  \end{array}
 
-
\right.
 
-
$
 
-
</center>
 
-
<center>
 
-
with $l\in[0,1]$
 
-
</center>
 
-
<p>
 
-
The variable 'tox' is representative of the amount of damages inflicted to bacteria, and so their probability to die. With each step of time (for us, a minute), bacteria cures a part of their injuries ($l<1$) and suffers other damages ($k'I(t)K(t)$). At first, few bacteria die : $\mbox{tox}(t)\cong k'I(t)K(t)$, then, 'tox' increases until it reach $\mbox{tox}(t)\cong\frac{k'}{1-l}I(t)K(t)$
 
-
</p>
 
-
<br>
 
-
<p> This new equation can also be writen in an analytic form :</p>
 
-
<center style="font-size:150%;">
 
-
$
 
-
\frac{dC}{dt}(t)=rC(t)-\int_{u=0}^t k'I(u)K(u)e^{\ln(l)(t-u)}du
 
-
$
 
-
</center>
 
-
<p> But it is unsolvable and brings only complexity.</p>
 
-
 
-
<br>
 
-
 
-
<p> With this model, we can explain properly the look of our curves : </p>
 
-
 
-
 
-
</li>
 
-
<li id="next"><a href="/Team:Grenoble-EMSE-LSU/Project/Modelling/Density">Next Page</a></li>
+
<li id="next"><a href="/Team:Grenoble-EMSE-LSU/Project/Modelling/Building">Next Page</a></li>
</ul>
</ul>
</div>
</div>
</html>
</html>

Latest revision as of 01:55, 5 October 2013

Grenoble-EMSE-LSU, iGEM


Grenoble-EMSE-LSU, iGEM

Retrieved from "http://2013.igem.org/Team:Grenoble-EMSE-LSU/Project/Modelling"