Paris/Continuous modelb
From 2007.igem.org
We developed simple models to evaluate the relative benefits of two different designs in term of robustness and optimization capabilities.
Contents |
Introduction
This model aims at describing the dynamic evolution of populations of germen and soma type bacteria. It is based on a set of differential equations describing DAP synthesis, DAP transport, differentiation of germen bacteria into soma and bacteria death. This approach differs form the precedent one by the level of description of the model and the numerical analysis done on the model.
By conducting robustness and optimization analysis on two different systems, one with a constant rate of differentiation against one with a rate of differentiation driven by the concentration of DAP, we evaluate benefits and drawbacks of both systems.
We first display the consistency of the biological system by providing a set of kinetic parameters such that the numerical simulation validate a given minimal behavior. Then we analyze the system robustness with regard to its kinetic parameters and finally we try to optimize the system output by adjusting some biologically relevant parameters.
Problem description
System
The system studied is made of two populations of bacteria G and S. Bacteria G can either divide or differentiate into S. Dap is only synthesized by bacteria S and diffuses in the environment and in G. G is dependent of Dap for division. Differentiation of G is controlled by Cre. Cre synthesis in G is either constant or dependant of Dap.
Expected behavior
We consider initial conditions in which there is some bacteria G and no soma type bacteria, so that the system first has to go through an intialisation phase where a part of the initial population of G differentiate into S. Adding Dap in the environment can only have a benefic effect on the population so we choose to start this study with an initial condition lacking Dap. The requested behavior for this system is that it grows fast enough and has a robust initialization phase.
Modeling purpose
- Is minimal behavior robust ?
- How can we tune the system to improve its output ?
- Which system should we prefer ?
Modelling
Differential equations
To model the dependence of division and differentiation of G we use Michaelis-Menten kinetics (terms λdiv and λdiff in equations 1 and 2). Equation (3) represents the synthesis of Cre and its degradation. We consider two models : the synthesis rate of Cre is either constant or inhibited by DAP (red term).
Term λdiv in equations (3,4) accounts for the dilution of Cre and DAP when the population of G grows.Equations involving DAP (4,5,6) represent DAP synthesis in somabacteria, DAP degradation in bacteria and in the environment, DAP transport and the release of DAP in the environment from dying bacteria. Terms λdiff in equations (4,5) are corrective terms to take into account the mean variation of concentration of DAP in G and S when some bacteria G differentiate into S.
