Team:Trieste/project/modeling
From 2012.igem.org
(modeling update) |
(add modeling) |
||
Line 2: | Line 2: | ||
<html> | <html> | ||
<div id="body"> | <div id="body"> | ||
- | + | <div id="container"> <!-- start container --> | |
- | + | </html>{{Team:Trieste/menu}}<html> | |
<div id="content"> <!-- start content --> | <div id="content"> <!-- start content --> | ||
- | + | <h1 id="h1_lf" class="main_tit"><div>Modeling</div></h1> | |
- | + | <h1 id="h1_rt" class="main_tit"><div>More</div></h1> | |
- | + | <div id="box_main"> <!-- start box_main --> | |
- | + | <div class="box_contenuti"> | |
- | + | <script type="text/javascript" | |
- | + | src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"> | |
+ | </script> | ||
+ | <h2>Model description</h2> | ||
+ | <p> | ||
We chose to use a deterministic, empirical model because it gives simple, | We chose to use a deterministic, empirical model because it gives simple, | ||
- | clean solutions. | + | clean solutions. |
- | + | </p> | |
- | + | <p> | |
The model is composed of 6 differential equations, each one describing | The model is composed of 6 differential equations, each one describing | ||
the concentration of a substance inside or outside the bacteria. | the concentration of a substance inside or outside the bacteria. | ||
- | + | </p> | |
- | + | <p> | |
All the $\delta$ parameters represent the variable's decay over time. | All the $\delta$ parameters represent the variable's decay over time. | ||
+ | </p> | ||
+ | <p> | ||
The first equation, describing the bacteria concentration in the media (bacteria growth) is based | The first equation, describing the bacteria concentration in the media (bacteria growth) is based | ||
on the logistics function: | on the logistics function: | ||
- | + | \begin{equation} | |
- | + | \dot{x} = x (1-x) | |
- | + | \end{equation} | |
We added a term to model bacteria death due to toxin, so the complete equation | We added a term to model bacteria death due to toxin, so the complete equation | ||
is: | is: | ||
- | + | \begin{equation} | |
- | + | \dot{b} = k_1 b (N - b) - k_2 b (l_i + l_e) | |
- | + | \end{equation} | |
There are three parameters that control this model's behaviour. | There are three parameters that control this model's behaviour. | ||
The $N$ parameter specifies the maximum concentration of bacteria. | The $N$ parameter specifies the maximum concentration of bacteria. | ||
The $k_1$ parameter controls the growth speed and the $k_2$ parameter sets how | The $k_1$ parameter controls the growth speed and the $k_2$ parameter sets how | ||
powerful the toxin is. | powerful the toxin is. | ||
- | + | </p> | |
- | + | ||
+ | <p> | ||
The next equation models the CymR concentration inside the bacteria. | The next equation models the CymR concentration inside the bacteria. | ||
- | + | \begin{equation} | |
- | + | \dot{p} = k_3 - \delta_1 p - \alpha p c | |
- | + | \end{equation} | |
The $k_3$ parameter determines the rate of synthesis of CymR and $\alpha$ the rate of cumate ligation. | The $k_3$ parameter determines the rate of synthesis of CymR and $\alpha$ the rate of cumate ligation. | ||
- | + | </p> | |
- | + | <p> | |
- | The cumate concentration function is in fact a parameter, because it is | + | The p-cumate concentration function is in fact a parameter, because it is |
- | directly controllable by adding cumate into the bacteria solution. | + | directly controllable by adding p-cumate into the bacteria solution. |
- | + | \begin{equation} | |
- | + | \dot{c} = k_4 - \delta_2 c | |
- | + | \end{equation} | |
- | The $k_4$ parameter models | + | The $k_4$ parameter models p-cumate increase (adding p-cumate at a steady rate). |
- | + | </p> | |
- | + | <p> | |
We use two equations to better model the toxin concentration, taking into | We use two equations to better model the toxin concentration, taking into | ||
account the concentration of toxin both inside and outside the bacteria. The equation | account the concentration of toxin both inside and outside the bacteria. The equation | ||
Line 60: | Line 66: | ||
concentration depends on bacteria death: when the bacteria dies, all the toxin | concentration depends on bacteria death: when the bacteria dies, all the toxin | ||
it has produced is released outside. | it has produced is released outside. | ||
- | + | \begin{align} | |
- | + | \dot{l_i} &= \frac{A}{k_5 + p} -\delta_3 l_i \\ | |
- | + | \dot{l_e} &= v (l_i + l_e) b l_i - \delta_4 l_e | |
- | + | \end{align} | |
The $A$ parameter sets the maximum toxin concentration inside the bacteria, while | The $A$ parameter sets the maximum toxin concentration inside the bacteria, while | ||
the $v$ parameter takes into account the toxin dilution in the bacteria | the $v$ parameter takes into account the toxin dilution in the bacteria | ||
solution. | solution. | ||
- | + | </p> | |
- | + | <p> | |
Last but not least, the antibody density, which we suppose always increasing in | Last but not least, the antibody density, which we suppose always increasing in | ||
time. | time. | ||
- | + | \begin{equation} | |
- | + | \dot{a} = k_6 - \delta_5 a | |
- | + | \end{equation} | |
- | + | </p> | |
- | + | </div> | |
- | + | </div> <!-- end box_main --> | |
- | + | <div id="box_right"> <!-- start box_right --> | |
- | + | <ul id="sub_menu"> | |
- | + | <li><a href="#">Sotto menu 1</a></li> | |
- | + | <li class="select"><a href="#">Sotto menu 2</a></li> | |
- | + | <li><a href="#">Sotto menu 3</a></li> | |
- | + | <li><a href="#">Sotto menu 4</a></li> | |
- | + | <li><a href="#">Sotto menu 5</a></li> | |
- | + | </ul> | |
- | + | <img src="https://static.igem.org/mediawiki/2012/b/b0/Team_trieste.jpg" alt="Team iGEM 2012" id="igem_team" /> | |
- | + | <div class="box_contacts"> | |
- | + | <h2>Contact us</h2> | |
- | + | <p>For other information, write to:</p> | |
- | + | <a href="mailto:igem2012@gmail.com" class="btn">igem2012@gmail.com</a> | |
- | + | <div class="social"> | |
- | + | Follow us also: | |
- | + | <a href="https://www.facebook.com/IGEMUNITS?ref=ts" target="_blank"><img src="https://static.igem.org/mediawiki/2012/f/f4/Ico_fb.png" alt="Facebook" class="fb" /></a> | |
- | + | <a href="https://twitter.com/igemunits" target="_blank"><img src="https://static.igem.org/mediawiki/2012/2/21/Ico_twitter.png" alt="twitter" class="tw" /></a> | |
- | + | </div> | |
- | + | </div> | |
- | + | </div> <!-- end box_right --> | |
- | + | </html>{{Team:Trieste/footer}}<html> | |
- | + | <!-- end content is included in template --> | |
- | + | </div> <!-- end container --> | |
</div> | </div> | ||
</html> | </html> |
Revision as of 20:16, 25 September 2012
Modeling
More
Model description
We chose to use a deterministic, empirical model because it gives simple, clean solutions.
The model is composed of 6 differential equations, each one describing the concentration of a substance inside or outside the bacteria.
All the $\delta$ parameters represent the variable's decay over time.
The first equation, describing the bacteria concentration in the media (bacteria growth) is based on the logistics function: \begin{equation} \dot{x} = x (1-x) \end{equation} We added a term to model bacteria death due to toxin, so the complete equation is: \begin{equation} \dot{b} = k_1 b (N - b) - k_2 b (l_i + l_e) \end{equation} There are three parameters that control this model's behaviour. The $N$ parameter specifies the maximum concentration of bacteria. The $k_1$ parameter controls the growth speed and the $k_2$ parameter sets how powerful the toxin is.
The next equation models the CymR concentration inside the bacteria. \begin{equation} \dot{p} = k_3 - \delta_1 p - \alpha p c \end{equation} The $k_3$ parameter determines the rate of synthesis of CymR and $\alpha$ the rate of cumate ligation.
The p-cumate concentration function is in fact a parameter, because it is directly controllable by adding p-cumate into the bacteria solution. \begin{equation} \dot{c} = k_4 - \delta_2 c \end{equation} The $k_4$ parameter models p-cumate increase (adding p-cumate at a steady rate).
We use two equations to better model the toxin concentration, taking into account the concentration of toxin both inside and outside the bacteria. The equation for the inside concentration is a Hill equation with a decay term. We used it because the concentration increases until the saturation is reached, and from then on the bacteria stops producing the toxin. The point of saturation is never reached, though, because the bacteria dies due to the toxin itself. The outside equation is the most complex one; the outside toxin concentration depends on bacteria death: when the bacteria dies, all the toxin it has produced is released outside. \begin{align} \dot{l_i} &= \frac{A}{k_5 + p} -\delta_3 l_i \\ \dot{l_e} &= v (l_i + l_e) b l_i - \delta_4 l_e \end{align} The $A$ parameter sets the maximum toxin concentration inside the bacteria, while the $v$ parameter takes into account the toxin dilution in the bacteria solution.
Last but not least, the antibody density, which we suppose always increasing in time. \begin{equation} \dot{a} = k_6 - \delta_5 a \end{equation}
![Team iGEM 2012](https://static.igem.org/mediawiki/2012/b/b0/Team_trieste.jpg)
![iGEM 2012](https://static.igem.org/mediawiki/2012/f/f0/Trieste-ombra_left.png)
![iGEM 2012](https://static.igem.org/mediawiki/2012/c/cb/Ombra_right.png)
![iGEM 2012](https://static.igem.org/mediawiki/2012/3/3a/Trieste-ombra_vert_left.png)
![iGEM 2012](https://static.igem.org/mediawiki/2012/3/3a/Trieste-ombra_vert_left.png)
![iGEM 2012](https://static.igem.org/mediawiki/2012/3/3a/Trieste-ombra_vert_left.png)
![iGEM 2012](https://static.igem.org/mediawiki/2012/b/b3/Trieste-ombra_vert_right.png)