Paris/Continuous model
From 2007.igem.org
Line 28: | Line 28: | ||
the cellular growth, so that at this time scale we can take d' = 0 and write d = s α5/ | the cellular growth, so that at this time scale we can take d' = 0 and write d = s α5/ | ||
α6 . This gives the two-variable system : | α6 . This gives the two-variable system : | ||
+ | |||
+ | <center>[[Image:eq2.jpg]]</center> | ||
+ | |||
+ | Redefinition of parameters k → kα6 /α5 and α3 → α2 + α3 leads to the simpler writing : | ||
+ | |||
+ | <center>[[Image:eq3.jpg]]</center> | ||
+ | |||
+ | Let us do some rewriting : | ||
+ | |||
+ | <center>[[Image:eq4.jpg]]</center> | ||
+ | |||
+ | and by redefining the time and most of the parameters we get : | ||
+ | |||
+ | <center>[[Image:eq5.jpg]]</center> | ||
+ | |||
+ | and by redefining the time and most of the parameters we get : | ||
+ | |||
+ | |||
+ | <center>[[Image:eq6.jpg]]</center> | ||
+ | |||
+ | The fixed points are (g = 0, s = 0) and the solution of : | ||
+ | |||
+ | |||
+ | <center>[[Image:eq7.jpg]]</center> | ||
+ | |||
+ | |||
+ | <center>[[Image:eq8.jpg]]</center> | ||
+ | |||
+ | |||
+ | <center>[[Image:eq9.jpg]]</center> | ||
+ | |||
+ | |||
+ | <center>[[Image:eq10.jpg]]</center> | ||
+ | |||
+ | |||
+ | <center>[[Image:eq11.jpg]]</center> | ||
+ | |||
+ | |||
+ | <center>[[Image:eq12.jpg]]</center> |
Revision as of 20:54, 23 October 2007
Well mixed population model
We present here a theoretical approach based on population dynamics. We consider here the case of a well mixed, homogenous, culture of the SMB organism, i.e. there is no space in this analysis and we follow only the variation of the different cell lines concentrations in the culture volume.
Derivation of the model
Let the variables g, s and d describe respectively the concentrations of germinal and so- matic cells and the conentration of DAP, then we can write for our system :
![Eq1.jpg](/wiki/images/5/5f/Eq1.jpg)
Equation (1a) describes the growth of the germinal cell population in presence of sufficient DAP (interaction represented by the Michaelis-Menten function) ; the term proportional to α2 is the differentiation into somatic cells by recombination of the CRE/LOX box, and the last term proportional to α3 is the germinal cells’ death. Equation (1b) is the variation of the somatic cells’ population, with the term proportional to α4 for somatic cells’ death. The last line describes DAP production by the somatic cells, and includes a degradation term. In absence of any quantitative details on assimilation of DAP by germinal cells and response to DAP levels, n and k are are to be considered as arbirary phenomenological parameters. We take however in the following n = 1 neglecting potential saturation related non-linearities for high DAP concentrations. The value of k corresponds to the DAP concentration for half-maximal growth rate, and could set experimenntally. We simplify the previous system by assuming that the evolution of d is rapid compared to the cellular growth, so that at this time scale we can take d' = 0 and write d = s α5/ α6 . This gives the two-variable system :
![Eq2.jpg](/wiki/images/4/4a/Eq2.jpg)
Redefinition of parameters k → kα6 /α5 and α3 → α2 + α3 leads to the simpler writing :
![Eq3.jpg](/wiki/images/0/01/Eq3.jpg)
Let us do some rewriting :
![Eq4.jpg](/wiki/images/0/07/Eq4.jpg)
and by redefining the time and most of the parameters we get :
![Eq5.jpg](/wiki/images/c/ca/Eq5.jpg)
and by redefining the time and most of the parameters we get :
![Eq6.jpg](/wiki/images/9/9b/Eq6.jpg)
The fixed points are (g = 0, s = 0) and the solution of :
![Eq7.jpg](/wiki/images/3/31/Eq7.jpg)
![Eq8.jpg](/wiki/images/c/cd/Eq8.jpg)
![Eq9.jpg](/wiki/images/b/b9/Eq9.jpg)
![Eq10.jpg](/wiki/images/d/d1/Eq10.jpg)