Introduction
- Aims of the modeling part
- First approach proposed
- Hill functions
- first model + score function
- bibliography
- findparam
- experiments
- Second approach
- bibliography
- equations
- results
- experiments
- Continue the previous model
- Synchronyzation
- Estimation of the FIFO processes
- Stochastic modeling (Gilespie)
- Test of robustness
- Enhancing the system
- Better FIFO behaviour
- Other interactions to increase the robustness
An Oscillatory Biological Model
Introduction
The goal here is to present the differential equations we used for our system modelization. At each step, we shall describe why we chose this precise model, its drawbacks and possible improvments, the parameters involved and enventually a biologically coherent value.
The key problem with a differential system consists in the fact that adding a new equation gives a more detailed idea of the overall process, but one looses precision by doing so since new parameters appear.
In that respect (in all cases but one) we chose not to model the mRNA steps, that is translation and transcription. We then assumed that we could act as if a protein would directly beget another one, without loosing too much precision.
This first approach only refers to a single cell. We shall examine later on what happens if we put more cells together.
Bibliography
We whall refer to those three articles :
[1] Shiraz Kalir, Uri Alon. Using quantitative blueprint to reprogram the dynamics of the flagella network. Cell, June 11, 2004, Vol.117, 713-720.
[2]Jordi Garcia-Ojalvo, Michael B. Elowitz, Steven H. Strogratz. Modeling a synthetic multicellular clock : repressilators coupled by quorum sensing. PNAS, July 27, 204, Vol. 101, no. 30.
[3]Nitzan Rosenfeld, Uri Alon. Response delays and the structure of transcription networks. JMB, 2003, 329, 645-654.
[4]Nitzan Rosenfeld, Michael B. Elowitz, Uri Alon. Negative autoregulation speeds the response times of transcription networks. JMB, 2003, 323, 785-793.
Equations
- flhDC ---> fliL ---> Fluorescent Protein 1 (FP1) _____ (1)
- flhDC ---> flgA ---> Fluorescent Protein 2 (FP2) _____ (2)
- flhDC ---> flhB ---> Fluorescent Protein 3 (FP3) _____ (3)
- flhDC ---> flhB ---> lasI (4)
- fliA ---> fliL ---> Fluorescent Protein 1 (FP1) (5)
- fliA ---> flgA ---> Fluorescent Protein 2 (FP2) (6)
- fliA ---> flhB ---> Fluorescent Protein 3 (FP3) (7)
- fliA ---> flhB ---> lasI (8)
For all these equations, we found in [1] that in that precise case, the promoter activity the seven class 2 operons, among which fLiL, flgA, flhB, may be mathematically described in that way :
where [X] denotes the effective protein-level activity at time.
For each operon, Shiraz Kalir and Uri Alon came up with numerical values of β and β', available in [1].
Furthermore, the protein-level activity can be presented (for a more detailed presentation, see[4]) as
Thus :
Even though we considered a single cell, we decided to model both HSL inside and outside the cell. In a first approach, we assumed that HSL could be modelized in the same fashion as AHL. The process was well detailed in [2].
- lasI ---> HSLext
- lasI ---> HSLint
Parameters summary
Graph screenshots
Roadmap
If you want to have a look at our roadmap : Roadmap
First Mathematical Approach
Bibliography
In order to choose a proper modeling approach for our system, we have decided to list all the chemical reactions we will take into account. Afterwards, we will find the needed parameters reading articles or devising the required experiments.
An overview of the work that has to be done can be found here : Bibliography
Estimation of parameters
Estimation of the parameters
First Approach
First Approach
More precise Bio-Mathematical Description
Bio-Mathematical Description
|