Team:Paris/Modeling/hill approach

From 2008.igem.org

(Difference between revisions)
(Bio-Chemical General Assumptions)
 
(14 intermediate revisions not shown)
Line 1: Line 1:
-
{{Paris/Menu}}
+
{{:Team:Paris/MenuBackup}}
= Model of the ''APE modelisation'' =
= Model of the ''APE modelisation'' =
==What kind of Mathematical Simulation ?==
==What kind of Mathematical Simulation ?==
-
We decided to use mostly Ordinary Differential Equation approach, at least for the study of the Oscillations and of the FIFO. For the Synchronisation module, we will probably use Probabilistic Differential Equations, in order to introduce the differences between the cells.
+
One of the strength of the synthetic biology is that precise knowledge and caracterisation of certains interactions allow very good predictions and simulations. So, our second model intends to get the best precision in the modelisation, consistent with the simpliest (but still logical) hypothesis possible.
-
We assume in this section that the concentrations of the products of complexation reactions are at steady-states (which is reached ''immediatly''), described as ''hill function'' (see [[Team:Paris/Modeling/Programs|precision on hill functions]]).
+
To determine interactions like Michaëlis-Menten's or Hill's, we start from the basical chemicals equations and try to caracterise their consequences on the behaviour of the system with few parameters. For instance, each ''complexation reactions'' will be caracterised at their steady-state, for all sets of initial concentrations (see [[Team:Paris/Modeling/Programs|complexations]]).
 +
 
 +
We decided to use mostly Ordinary Differential Equation approach, at least for the study of the Oscillations and of the FIFO. For the Synchronisation module, we will probably use Probabilistic Differential Equations, in order to introduce the differences between the cells.
==Bio-Chemical General Assumptions==
==Bio-Chemical General Assumptions==
-
We know that the following equations do not describe properly what ''really'' happens in the cells. For exemple, we know that the transcription factor flhD-flhC is actually an ''hexamere'' FlhD<sub>4</sub>C<sub>2</sub>. But, as we will surely not get access to the ''dissociation constant'' of the ''hexamerisation'', we just treat it as an ''abstract'' inducer protein "FlhDC", with an order in its ''Hill function'' probably between 3 and 6 (but perhaps completly different; the estimation of the error by the [[Team:Paris/Modeling/Programs|'findparam']] program will tell us if we are right to do so).
+
We know that the following equations do not describe properly what ''really'' happens in the cells. For exemple, we know that the transcription factor flhD-flhC is actually an ''hexamere'' FlhD<sub>4</sub>C<sub>2</sub>. But, as we will surely not get access to the ''dissociation constant'' of the ''hexamerisation'', we just treat it as an ''abstract'' inducer protein "FlhDC", with an order (''n'') in its ''complexation caracterisation'' probably around 2*4 = 8 (but perhaps completly different ! ; the estimation of the error by the [[Team:Paris/Modeling/Programs#Finding_Parameters|'findparam']] program will tell us if we are right to do so).
For the moment, at each part of our modelisation, we reduce the expression of a gene at its '''transcription'''. The '''translation''' process is not taken into acount (see however [[Team:Paris/Modeling/estimation#RBS_Issue|considerations on RBS]]).
For the moment, at each part of our modelisation, we reduce the expression of a gene at its '''transcription'''. The '''translation''' process is not taken into acount (see however [[Team:Paris/Modeling/estimation#RBS_Issue|considerations on RBS]]).
 +
Besause of that, our caracterisations doesn't allow us to know the ''real concentrations'' of the proteins we produce, but only their "real times effects" on the promoters they influence.
 +
 +
Every complexations we deal with are the supposed to reach ''immediatly'' their '''steady-states'''. So, the only phenomenon we are observing along time is the protein production ; that's why this model can't really simulate the possible delay a complexation could introduce, and which could be important for the stiffness of our oscillations...
As we need to keep our cells in the exponential phase of growth (and since we can't use the already mentioned system of Ron Weiss, see [[Team:Paris/Project|description of the project]]), our system works in a '''chemostat'''. We will also be able to estimate the Cell Density, and we will have to take into acount the ''renewal phenomenon''. Under this conditions, we assume the following constants to be true :
As we need to keep our cells in the exponential phase of growth (and since we can't use the already mentioned system of Ron Weiss, see [[Team:Paris/Project|description of the project]]), our system works in a '''chemostat'''. We will also be able to estimate the Cell Density, and we will have to take into acount the ''renewal phenomenon''. Under this conditions, we assume the following constants to be true :
* Cell Division : every 35 minutes
* Cell Division : every 35 minutes
-
** Dilution Rate : 0.0198 min<sup>-1</sup>
+
&rarr; Dilution Rate : 0.0198 min<sup>-1</sup>
-
* Cell Density : XXX cell.L<sup>-1</sup>
+
* Cell Density : cell.L<sup>-1</sup>
-
** Average Intracellular Volume : XXX L
+
&rarr; Average Intracellular Volume : L <br>
-
** Average Extracellular Volume (in the chemostat) : XXX L
+
&rarr; Average Extracellular Volume (in the chemostat) : L
-
* Renewal Rate : XXX L.min<sup>-1</sup>
+
* Renewal Rate : L.min<sup>-1</sup>
-
** Dilution Rate : XXX min<sup>-1</sup>
+
(&rarr; Dilution Rate : min<sup>-1</sup>)
-
To see more details about the modelisation and the values of the involved constants, see [[Team:Paris/Modeling/Bibliography|the bibliography]] and the [[Team:Paris/Modeling/estimation|estimation section]].
+
To see more details about the values of the involved constants, see [[Team:Paris/Modeling/Bibliography|the bibliography]] and the [[Team:Paris/Modeling/estimation|estimation section]].
==Incrementally detailed Parts of our Project==
==Incrementally detailed Parts of our Project==

Latest revision as of 07:06, 30 October 2008

You are currently on the Wiki Museum
Go back to the normal Wiki

Contents

Model of the APE modelisation

What kind of Mathematical Simulation ?

One of the strength of the synthetic biology is that precise knowledge and caracterisation of certains interactions allow very good predictions and simulations. So, our second model intends to get the best precision in the modelisation, consistent with the simpliest (but still logical) hypothesis possible.

To determine interactions like Michaëlis-Menten's or Hill's, we start from the basical chemicals equations and try to caracterise their consequences on the behaviour of the system with few parameters. For instance, each complexation reactions will be caracterised at their steady-state, for all sets of initial concentrations (see complexations).

We decided to use mostly Ordinary Differential Equation approach, at least for the study of the Oscillations and of the FIFO. For the Synchronisation module, we will probably use Probabilistic Differential Equations, in order to introduce the differences between the cells.

Bio-Chemical General Assumptions

We know that the following equations do not describe properly what really happens in the cells. For exemple, we know that the transcription factor flhD-flhC is actually an hexamere FlhD4C2. But, as we will surely not get access to the dissociation constant of the hexamerisation, we just treat it as an abstract inducer protein "FlhDC", with an order (n) in its complexation caracterisation probably around 2*4 = 8 (but perhaps completly different ! ; the estimation of the error by the 'findparam' program will tell us if we are right to do so).

For the moment, at each part of our modelisation, we reduce the expression of a gene at its transcription. The translation process is not taken into acount (see however considerations on RBS). Besause of that, our caracterisations doesn't allow us to know the real concentrations of the proteins we produce, but only their "real times effects" on the promoters they influence.

Every complexations we deal with are the supposed to reach immediatly their steady-states. So, the only phenomenon we are observing along time is the protein production ; that's why this model can't really simulate the possible delay a complexation could introduce, and which could be important for the stiffness of our oscillations...

As we need to keep our cells in the exponential phase of growth (and since we can't use the already mentioned system of Ron Weiss, see description of the project), our system works in a chemostat. We will also be able to estimate the Cell Density, and we will have to take into acount the renewal phenomenon. Under this conditions, we assume the following constants to be true :

  • Cell Division : every 35 minutes

→ Dilution Rate : 0.0198 min-1

  • Cell Density : cell.L-1

→ Average Intracellular Volume : L
→ Average Extracellular Volume (in the chemostat) : L

  • Renewal Rate : L.min-1

(→ Dilution Rate : min-1)

To see more details about the values of the involved constants, see the bibliography and the estimation section.

Incrementally detailed Parts of our Project