IGEM:IMPERIAL/2006/project/Oscillator/Theoretical Analyses
From OpenWetWare
< IGEM:IMPERIAL | 2006 | project | Oscillator
Contents |
Analysis of the Model of the Molecular Predation Oscillator
- Introduction
- Our Approach
- Model Simplication
- Our Results
- Conclusion
- Appendix
Introduction
- The present part of the i-coli Reporter deals with the analysis of a complex 3 dimensional dynamical system that model the dynamic of the Molecular Predation System
- Since the analysis was long and complex, we have split its results in different sections that are accessible with the blue tabs above.
- These different sections present the logical progression in our analysis
- We suggest you browse them in the order shown above
- We hope you enjoy your reading
Presentation of our Dynamical System
- The molecular predator-prey system has been modelled by the following 3D Dynamical System
- The molecular predator-prey system has been modelled by the following 3D Dynamical System
- As is often the case in Mathematics, we renamed the parameters in order to make the system more symmetric and easier to interprete. We formalised the dynamical system as follows
How to read the model (the basics)
- The variables U,V and W respectively stand for the concentrations of AHL, aiiA and LuxR.
- The time derivatives of U,V and W stand for their growth rates.
- Their bounded growth terms is due to gene expression.
- Degradation of the prey (AHL) is partly due to an enzymatic reaction
- Finally The model comprises dissipative terms of the model (eU,d1V,d2W) due to the washout in the chemostat.
- Different types of parameters
- - ao,bo and co are assumed constant (that is without our control)
- - The washout parameters (d1,d2 and e) are assumed within our control. However, experimental contraints prevent them from getting very small or very large
- - a, b and c are within our control. Out of design (2 cell system) they vary accordingly to the concentrations of cells.
General Remarks on the Model of Interest
- In the analysis part we only deal with the model presented above.
- We do not consider its possible (and natural) extensions with an exponent in the model and with leakage terms.
- From an experimental point of view, restricting our analysis makes sense. We have built test constructs to measure the growth of AHL, aiiA and LuxR as well as the degradation of AHL by aiiA. Our experimental results generally showed that our model was valid.
- However, we are very aware of the pitfalls of limiting our analysis, but due to the time constrain, we will not tackle them this year.