Team:Paris/Analysis/Construction2

From 2008.igem.org

(Difference between revisions)
(Kinetics)
(Introduction)
 
(209 intermediate revisions not shown)
Line 1: Line 1:
{{Paris/Menu}}
{{Paris/Menu}}
-
{{Paris/Header|Model Construction}}
+
{{Paris/Header|System Improvements Description}}
-
== Description ==
+
{{Paris/Section_contents_analysis}}
-
<!--why use a chemostat-->
+
-
* We chose to use a chemostat. We want to impose to our model the fact that the rate of production has to be proportional to the existing population and to the amount of available resources.
+
-
<!--which model of chemostat we choose-->
+
== Introduction ==
-
* We assume a logistic model for the population kinetics in the chemostat (reference).
+
-
<!--Autoinducer -->
+
In this section, we investigate possible improvements of the [[Team:Paris/Network_analysis_and_design/Core_system/Model_construction|core system]]. Our objective is twofold : the system has to provide sustained oscillations and these oscillations should be synchronized amongst a population of cells. To this aim we explore designs inspired by quorum sensing and model in a chemostat cell growth and species diffusion outside cells. We consider two models relying on different designs principles.
-
* To archieve synchronization we use QS. Explain the principle?
+
-
== Kinetics ==
+
{{Paris/Toggle|To join our quest for alternatives click here!|Team:Paris/DescriptionMoreModelConstruction2}}
-
* Following a logistic model of the population in the chamostat, the concentration of cells in the medium can be expressed in terms of a production (positive) term and degradation (negative) terms, as:
+
-
<br>
+
-
[[Image:population_evolution_full_2.jpg|center]]
+
-
<br>
+
-
On the one hand, the concentration of cells (c) increase exponentialy with a velocity given by the growth rate (α<sub>cell</sub>). After a certain time, the population reaches a maximum concentration (c<sub>max</sub>).
+
-
On the other hand, c is reduced due to the dilution cause by the renewal of the chemostat (D<sub>renewal</sub>) and the cell death (d).
+
 +
== Modeling Alternatives ==
-
*To model the quorum sensing dynamics, we use:
+
The proposed systems are:
 +
{|align="center" cellpadding="0" cellspacing="20"
 +
|-style="background: #D4E2EF; text-align: center;"
 +
|'''HSL mediated coupled oscillators'''
 +
|'''HSL mediated simple oscillator'''
 +
|-style="background: #ffffff; text-align: center;"
 +
|[[Image:Bimo.png|400px]]
 +
|[[Image:Unimo.png|400px]]
 +
|}
 +
*In the first system, we use a modular design. Since the core system is a 'poor' oscillator we hope that by coupling it with a 'good' oscillator we could obtain sustained oscillations in the whole system. Thus, we consider that the core system is one of the modules of the system and that the other module is a two gene oscillator system presented in [[Team:Paris/Bibliography|[2]]] that accounts for quorum sensing. We call this alternative the 'HSL mediated coupled oscillators'.
-
Production of HSL:
+
*In the second system, namely the 'HSL mediated simple oscillator', we rewire the architecture of the core system to introduce delay via HSL export in the environment in a single circuit.
-
eqHSL 1
+
Both the coupled and simple oscillators describe events that happen not only at the cellular level (as in the core system) but also at the population level due to interactions needed between a cell and its environment.
-
Transport of HSL:
+
In the following sections, we first describe the [[Team:Paris/Network_analysis_and_design/System_improvements/Description#Common_Description|common]] part among the two proposed models and then focus our attention to the  [[Team:Paris/Network_analysis_and_design/System_improvements/Description#Alternatives_Description| alternatives description]] where the characteristics that are specific to each of the modeling alternatives are presented.
 +
<!-- MORE DETAILS ABOUT THE CHEMOST
 +
* A chemostat is generally used to keep bacteria volume constant in the medium. The constant conditions provided by the chemostat help us to control bacteria growth rate. 
-
eqHSL 2
+
* We assume a logistic model to determine bacteria growth in the chemostat, in agreement with standar procedures (reference). This hypothesis implies that bacteria growth rate has to be proportional to the existing bacteria population size and to the amount of available resources in the medium.
 +
-->
 +
<!-- MORE DETAILS ABOUT QS
 +
* As discussed [[Team:Paris/Network_analysis_and_design/System_improvements|previously]], we choose to use quorum sensing as a way to improve the oscillating behaviour of our [[Team:Paris/Network_analysis_and_design/Core_system|core system]] and, at the same time, as a way to archeive population synchronization. When a cultive of bacteria is synchronized, it means that every single cell express in average the same genes in unison. As a result, a maximum level of fluorescense is obtained.
 +
-->
-
Activation of ''envZ'' depends of the concentration of HSL:
+
== Common Description ==
 +
{|align="center"
 +
|-style="background: #ffffff; text-align: left;"
 +
|Logistic growth in chemostat:
 +
|-style="background: #ffffff; text-align: center;"
 +
|[[Image:DescriptionCommonDynamicsPart1.png|center|600px]]
 +
|-style="background: #ffffff; text-align: left;"
 +
|Quorum sensing by HSL diffusion:
 +
|-style="background: #ffffff; text-align: center;"
 +
|[[Image:DescriptionCommonDynamicsPart2.png|center|600px]]
 +
|-style="background: #ffffff; text-align: left;"
 +
|As in [[Team:Paris/Network_analysis_and_design/Core_system|core system]] flhDC inhibited by EnvZ and fliA activated by FlhDC and FliA:
 +
|-style="background: #ffffff; text-align: center;"
 +
|[[Image:DescriptionCommonDynamicsPart3.png|center|600px]]
 +
|}
 +
{{Paris/Toggle|read more on common dynamics...|Team:Paris/DescriptionDetailsS4-S2Part1}}
-
eqHSL 3
+
== Alternatives Description ==
 +
{|align="center" cellpadding="2" cellspacing="10"
 +
|-style="background: #D4E2EF; text-align: center;"
 +
|'''HSL mediated coupled oscillators'''
 +
|'''HSL mediated simple oscillator'''
 +
|-style="background: #ffffff; text-align: left;"
 +
|Core system coupled with an oscillator:
 +
|Modified core system that accounts for quorum sensing:
 +
|-style="background: #ffffff; text-align: center;"
 +
|[[Image:SumaryBiMo.png|450px]]
 +
|[[Image:SumaryUniMo.png|350px]]
 +
|}
 +
{{Paris/Toggle|read more on alternatives description...|Team:Paris/DescriptionDetailsS4-S2Part2}}
-
<!-- equations in code:
+
== Kinetic parameter values==
-
d[c]/dt=alpha_cell*[c]*(c_max-[c])/c_max-(D_renewal+death)*[c]
+
-
d[LasI]/dt=beta_FlhDC*theta_FlhDC^n_FlhDC/(theta_FlhDC^n_FlhDC+[TetR]^n_FlhDC)-gamma*[LasI]
+
-
d[TetR]/dt=beta_TetR*[HSL]^n_TetR/(theta_TetR^n_TetR+[HSL]^n_TetR)-gamma*[TetR]
+
-
d[HSL::ext]/dt=eta*v_cell*[c]*([HSL]-[HSL::ext])-(D_renewal+gamma_HSLext)*[HSL::ext]
+
-
d[HSL]/dt=beta_HSL*[LasI]+eta*[HSL::ext]-(gamma+gamma_HSLint+eta)*[HSL]
+
-
-->
+
-
== Parameters Search==
+
Remarkably, almost all parameter values are available from experimental measurements in [[Team:Paris/Bibliography|[1]]] and from the work of Garcia-Ojalvo in [[Team:Paris/Bibliography|[2]]]. Additionally minimizing the number of parameters is possible by rescaling as done for the core system. Relevant values for the sole two missing parameters are found by exploiting a modularity assumption as described in the  [[Team:Paris/Network_analysis_and_design/System_improvements/Analysis|next section]]. The following table summarize our findings:
-
manly from literature but also from  [[Team:Paris/Network_analysis_and_design/System_improvements/Description/QS_Module|S0 analysis]].
 
<center>
<center>
Line 54: Line 80:
! colspan="8" | &nbsp;
! colspan="8" | &nbsp;
|- style="background: #649CD7;" |
|- style="background: #649CD7;" |
-
|style="background: #ffffff; text-align:center;" |'''Chemostat'''
+
|style="background: #ffffff; text-align:center;" width=20% |'''Chemostat'''
|Parameter
|Parameter
|Meaning
|Meaning
Line 62: Line 88:
|Source  
|Source  
|- style="background: #dddddd;"|     
|- style="background: #dddddd;"|     
-
| ! rowspan="5" style="background: #ffffff;" | figure / equations of Chemostat
+
| ! rowspan="5" style="background: #ffffff;" |  
|style="background: #D4E2EF;" | α<sub>cell</sub>
|style="background: #D4E2EF;" | α<sub>cell</sub>
|style="background: #dddddd;"| Growth rate
|style="background: #dddddd;"| Growth rate
Line 68: Line 94:
| style="background: #dddddd;"|1
| style="background: #dddddd;"|1
|style="background: #dddddd;"|min<sup>-1</sup>
|style="background: #dddddd;"|min<sup>-1</sup>
-
|style="background: #dddddd;"| [[Team:Paris/Modeling/BOB#Remarks|wet-lab]]  
+
|style="background: #dddddd;"| [[Team:Paris/Remarks|wet-lab]]  
|- style="background: #dddddd;"|     
|- style="background: #dddddd;"|     
| style="background: #D4E2EF;"| c<sub>max</sub>
| style="background: #D4E2EF;"| c<sub>max</sub>
Line 75: Line 101:
| style="background: #dddddd;"| 0.1
| style="background: #dddddd;"| 0.1
| style="background: #dddddd;"| µm<sup>3</sup>
| style="background: #dddddd;"| µm<sup>3</sup>
-
| style="background: #dddddd;"| [[Team:Paris/Modeling/BOB#Bibliography|[3]]]
+
| style="background: #dddddd;"| [[Team:Paris/Bibliography|[3]]]
|- style="background: #dddddd;"|     
|- style="background: #dddddd;"|     
| style="background: #D4E2EF;"| D<sub>renewal</sub>
| style="background: #D4E2EF;"| D<sub>renewal</sub>
Line 82: Line 108:
| style="background: #dddddd;"| 0.1
| style="background: #dddddd;"| 0.1
| style="background: #dddddd;"| min<sup>-1</sup>
| style="background: #dddddd;"| min<sup>-1</sup>
-
| style="background: #dddddd;"| [[Team:Paris/Modeling/BOB#Remarks|wet-lab]] ([[Team:Paris/Modeling/BOB#Bibliography|[3]]])
+
| style="background: #dddddd;"| [[Team:Paris/Remarks|wet-lab]] ([[Team:Paris/Bibliography|[3]]])
|- style="background: #dddddd;"|     
|- style="background: #dddddd;"|     
| style="background: #D4E2EF;"| d
| style="background: #D4E2EF;"| d
Line 89: Line 115:
| style="background: #dddddd;"| 0.5
| style="background: #dddddd;"| 0.5
| style="background: #dddddd;"| min<sup>-1</sup>
| style="background: #dddddd;"| min<sup>-1</sup>
-
| style="background: #dddddd;"|[[Team:Paris/Modeling/BOB#Remarks|wet-lab]]
+
| style="background: #dddddd;"|[[Team:Paris/Remarks|wet-lab]]
|- style="background: #ffffff;"
|- style="background: #ffffff;"
! colspan="8" | &nbsp;
! colspan="8" | &nbsp;
Line 101: Line 127:
|Source  
|Source  
|- style="background: #dddddd;"|     
|- style="background: #dddddd;"|     
-
| ! rowspan="6" style="background: #ffffff;" | figure / equations of Quorum Sensing
+
| ! rowspan="6" style="background: #ffffff;" |  
| style="background: #D4E2EF;" |γ<sub>HSL</sub>  
| style="background: #D4E2EF;" |γ<sub>HSL</sub>  
| style="background: #dddddd;"|Degradation rate
| style="background: #dddddd;"|Degradation rate
Line 107: Line 133:
| style="background: #dddddd;"|0.2690
| style="background: #dddddd;"|0.2690
| style="background: #dddddd;"|min<sup>-1</sup>
| style="background: #dddddd;"|min<sup>-1</sup>
-
| style="background: #dddddd;"|[[Team:Paris/Modeling/BOB#Remarks|wet-lab]]  
+
| style="background: #dddddd;"|[[Team:Paris/Remarks|see remarks]]  
|- style="background: #dddddd;"|  
|- style="background: #dddddd;"|  
| style="background: #D4E2EF;"|γ<sub>HSL<sub>ext</sub></sub> <br> <br>
| style="background: #D4E2EF;"|γ<sub>HSL<sub>ext</sub></sub> <br> <br>
Line 114: Line 140:
| style="background: #dddddd;"|0.5380
| style="background: #dddddd;"|0.5380
|style="background: #dddddd;"| min<sup>-1</sup>
|style="background: #dddddd;"| min<sup>-1</sup>
-
| style="background: #dddddd;"|[[Team:Paris/Modeling/BOB#Bibliography|[6]]]
+
| style="background: #dddddd;"|[[Team:Paris/Bibliography|[2]]]
|- style="background: #dddddd;"|  
|- style="background: #dddddd;"|  
| style="background: #D4E2EF;" |β<sub>HSL</sub>  
| style="background: #D4E2EF;" |β<sub>HSL</sub>  
Line 121: Line 147:
| style="background: #dddddd;"|16
| style="background: #dddddd;"|16
| style="background: #dddddd;"|min<sup>-1</sup>
| style="background: #dddddd;"|min<sup>-1</sup>
-
| style="background: #dddddd;"|[[Team:Paris/Modeling/BOB#Bibliography|∅]]
+
| style="background: #dddddd;"|[[Team:Paris/Bibliography|∅]]
|- style="background: #dddddd;"|  
|- style="background: #dddddd;"|  
| style="background: #D4E2EF;" |η
| style="background: #D4E2EF;" |η
Line 128: Line 154:
| style="background: #dddddd;"|505
| style="background: #dddddd;"|505
| style="background: #dddddd;"|min<sup>-1</sup>
| style="background: #dddddd;"|min<sup>-1</sup>
-
| style="background: #dddddd;"|[[Team:Paris/Modeling/BOB#Bibliography|[2]]]
+
| style="background: #dddddd;"|[[Team:Paris/Bibliography|[2]]]
|- style="background: #dddddd;"|  
|- style="background: #dddddd;"|  
| style="background: #D4E2EF;" |n<sub>HSL</sub>
| style="background: #D4E2EF;" |n<sub>HSL</sub>
Line 135: Line 161:
| style="background: #dddddd;"|4
| style="background: #dddddd;"|4
| style="background: #dddddd;"|&nbsp;
| style="background: #dddddd;"|&nbsp;
-
| style="background: #dddddd;"|[[Team:Paris/Modeling/BOB#Bibliography|[3]]]
+
| style="background: #dddddd;"|[[Team:Paris/Bibliography|[3]]]
|- style="background: #dddddd;"|  
|- style="background: #dddddd;"|  
| style="background: #D4E2EF;" |θ<sub>HSL</sub>
| style="background: #D4E2EF;" |θ<sub>HSL</sub>
Line 141: Line 167:
| style="background: #dddddd;"|&nbsp;
| style="background: #dddddd;"|&nbsp;
| style="background: #dddddd;"|0.5
| style="background: #dddddd;"|0.5
-
| style="background: #dddddd;"|[[Team:Paris/Modeling/BOB#Remarks|c.u]]
+
| style="background: #dddddd;"|[[Team:Paris/Remarks|c.u.]]
-
| style="background: #dddddd;"|[[Team:Paris/Modeling/BOB#Bibliography|[3]]]
+
| style="background: #dddddd;"|[[Team:Paris/Bibliography|[3]]]
|}</center>
|}</center>
 +
 +
 +
{{Paris/Navig|Team:Paris/Analysis}}

Latest revision as of 03:18, 30 October 2008

System Improvements Description


Contents

Introduction

In this section, we investigate possible improvements of the core system. Our objective is twofold : the system has to provide sustained oscillations and these oscillations should be synchronized amongst a population of cells. To this aim we explore designs inspired by quorum sensing and model in a chemostat cell growth and species diffusion outside cells. We consider two models relying on different designs principles.

↓ To join our quest for alternatives click here! ↑


Same environment
Chemostat.png
Using a chemostat
Same starting point, different synchronization mechanisms
HSL mediated simple oscillator (learn more...)
HSL mediated coupled oscillators (learn more...)

Modeling Alternatives

The proposed systems are:

HSL mediated coupled oscillators HSL mediated simple oscillator
Bimo.png Unimo.png
  • In the first system, we use a modular design. Since the core system is a 'poor' oscillator we hope that by coupling it with a 'good' oscillator we could obtain sustained oscillations in the whole system. Thus, we consider that the core system is one of the modules of the system and that the other module is a two gene oscillator system presented in [2] that accounts for quorum sensing. We call this alternative the 'HSL mediated coupled oscillators'.
  • In the second system, namely the 'HSL mediated simple oscillator', we rewire the architecture of the core system to introduce delay via HSL export in the environment in a single circuit.

Both the coupled and simple oscillators describe events that happen not only at the cellular level (as in the core system) but also at the population level due to interactions needed between a cell and its environment.

In the following sections, we first describe the common part among the two proposed models and then focus our attention to the alternatives description where the characteristics that are specific to each of the modeling alternatives are presented.

Common Description

Logistic growth in chemostat:
DescriptionCommonDynamicsPart1.png
Quorum sensing by HSL diffusion:
DescriptionCommonDynamicsPart2.png
As in core system flhDC inhibited by EnvZ and fliA activated by FlhDC and FliA:
DescriptionCommonDynamicsPart3.png
↓ read more on common dynamics... ↑


Common Dynamics: Chemostat
The variation of cells' concentration in the chemostat over time can be expressed in terms of a production (positive) term and degradation (negative) terms:
Final chemostat.png
    For the production term, we use a logistic equation to model cell growth, according to standard assumptions. The behaviour obtained is the following one: at low population density, the concentration of cells in the chemostat (c) increase exponentially with a growth rate αcell and at high population density, the population reaches a maximum concentration, cmax.
    For the degradation term, we consider that c decrease proportionally to both a dilution phenomena cause by the renewal of the medium in the chemostat (Drenewal) and cell death (d).
Common Dynamics: Quorum Sensing
In order to model the quorum sensing dynamics, we consider that:
     1) Inside a cell, the HSL concentration increases proportionally to the concentration of LasI and decreases according to both a degradation term (proportional to the internal HSL concentration) and a transport term (proportional to the difference between the internal and external concentration of HSL). Thus, the equation for the internal HSL concentration is:
Final HSL.png
     2) Outside the cells, HSL is accumulated with the same transport term that we use in the previous equation. The degradation of HSL in the external medium and the dilution controlled via the chemostat accounts for HSL external decrease. So the external HSL concentration is given by:
HSLext.png
that is equivalent to:
HSLext2.png
where    Hsl average.png    and     Cell number volume.png
Common Network Dynamics: FlhDC and Flia
FlhDC and Flia are regulated in the same way in both systems. FlhDC is produced under the influence of EnvZ via an inhibition. Flia is regulated for its self and FlhDC.
FlhDC-FliaEq.png

Alternatives Description

HSL mediated coupled oscillators HSL mediated simple oscillator
Core system coupled with an oscillator: Modified core system that accounts for quorum sensing:
SumaryBiMo.png SumaryUniMo.png
↓ read more on alternatives description... ↑


HSL mediated coupled oscillators HSL mediated simple oscillator
     a) The expression of lasI is under the control of the same promotor that used for FlhDC.      a) The expression of lasI is regulated by FlhDC and Flia (as in the core system)
LasIEqInBIMOdularSys.png LasIEqInUNIModularSys.png
     b) The expression of envZ depends on both the activation from FlhDC and Flia (as in the core system) and the concentration of HSL present in the cell.      b) The expression of envZ depends only on the concentration of HSL present in the cell.
EnvZInBIMOdularSys.png EnvZInUNIModularSys.png

Kinetic parameter values

Remarkably, almost all parameter values are available from experimental measurements in [1] and from the work of Garcia-Ojalvo in [2]. Additionally minimizing the number of parameters is possible by rescaling as done for the core system. Relevant values for the sole two missing parameters are found by exploiting a modularity assumption as described in the next section. The following table summarize our findings:


Parameters
 
Chemostat Parameter Meaning Original Value Normalized Value Unit Source
αcell Growth rate 0.0198 1 min-1 wet-lab
cmax Carrying capacity for cell growth 0.1 0.1 µm3 [3]
Drenewal Dilution rate 0.00198 0.1 min-1 wet-lab ([3])
d Death rate 0.0099 0.5 min-1 wet-lab
 
Quorum Sensing Parameter Meaning Original Value Normalized Value Unit Source
γHSL Degradation rate 0.0053 0.2690 min-1 see remarks
γHSLext

Degradation rate 0.0106 0.5380 min-1 [2]
βHSL Production rate 0.3168 16 min-1
η Diffusion rate 10 505 min-1 [2]
nHSL Hill coefficient   4   [3]
θHSL Hill characteristic concentration for the second operator   0.5 c.u. [3]