Team:Paris/Modeling

From 2008.igem.org

(Difference between revisions)
(An Oscillatory Biological Model)
 
(186 intermediate revisions not shown)
Line 1: Line 1:
-
{{Paris/Menu}}
+
{{:Team:Paris/MenuBackup}}
-
{|cellspacing="5" cellpadding="10" style="background:#649CD7; width: 965px;"
+
-
|-valign="top"
+
-
|style="background:#ffffff"|
+
-
== Introduction ==
 
-
* Aims of the modeling part
+
= Our train of thoughts... =
-
* First approach proposed
+
We hereby propose different and complementary approaches to model the biological system. We found interesting to explain two of the paths that we chose to follow in order to understand and predict our system. It is important to note that both models aim at different goals in the process of understanding our system.
-
** Hill functions
+
Furthermore, we wished to describe our thought process, the way these models interact, their respective roles. 
-
** first model + score function
+
An overall description of the way we model our biological system can be found below :
-
** bibliography
+
<center>[[Team:Paris/Modeling/History|Read more !]]</center>
-
** findparam
+
-
**experiments
+
-
*Second approach
+
-
**bibliography
+
-
**equations
+
-
**results
+
-
**experiments
+
-
* Continue the previous model
+
-
**Synchronyzation
+
-
**Estimation of the FIFO processes
+
-
**Stochastic modeling (Gilespie)
+
-
*Test of robustness
+
-
**repressilator
+
-
**comparison
+
-
*Enhancing the system
+
-
**Better FIFO behaviour
+
-
**Other interactions to increase the robustness
+
-
===Roadmap===
+
= BOB (Based On Bibliography) Approach =
 +
[[Image:BOB.jpg|250px|thumb]]
-
If you want to have a look at our modeling notebook: [[Team:Paris/Modeling/Roadmap|Notebook]]
+
Our first approach is quite rough and simple but effective. The goal here was to guess the behavior of our Bacteri'OClock, considering the overall system. Since it is a preliminary approach, we could not yet fill the model with data from the wet lab. This is why our work is mainly based on a bibliographic work, which allows us to use parameters and data from scientific articles.
-
==An Oscillatory Biological Model==
+
The key points of this approach:
-
[[Image:Globalmodel.jpg|frame|<center><b>Overall Model</b></center>|center]]
+
* Simplicity for itself is not that important. In fact, what we were looking for was understandability at first.
 +
* We used linear equations as much as possible: wherever it had been proved in a paper than an interaction could be efficiently modeled with a elementary expression, we kept it.
 +
* Too many parameters in a model mean less relevancy. In addition, the more parameters you have, the hardest it is to tune the system in order to have the behavior you are looking for.
-
===Introduction===
 
-
<br>The goal here is to present the differential equations we used in the 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.
+
<center>[[Team:Paris/Modeling/BOB|Read more]]</center>
-
<br>The key problem with a dynamic 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.
+
-
<br>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 act upon the follower, without loosing too much precision.
+
-
<br>This first approach only refers to a single cell. We shall examine later on what happens if we put more cells together.
+
-
===Bibliography===
+
= APE (APE Parameters Estimation) Approach=
-
Musch of our inspiration comes from four articles to which we shall refer in the next subsections :
+
[[Image:APE.jpg|250px|thumb]]
-
<br>
+
The second approach was motivated by our will to characterize our system in the most precise way. What is at stake here is to determine the "real parameters" that govern the dynamics of our system.
-
* [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===
+
* Each step is taken into account at a fundamental kinetic processes level or at a more global level by a function describing the complexation, which is a simple way to take into account multiple interactions and more especially cooperative binding.
-
* <b>Steps</b>
+
-
**  flhDC &rarr; fliA &nbsp;&nbsp;&nbsp; (1)
+
-
**  flhDC &rarr; fliL &rarr; Fluorescent Protein 1 (FP1)  &nbsp;&nbsp;&nbsp;    (2)
+
-
**  flhDC &rarr; flgA &rarr; Fluorescent Protein 2 (FP2)  &nbsp;&nbsp;&nbsp;    (3)
+
-
**  flhDC &rarr; flhB &rarr; Fluorescent Protein 3 (FP3)  &nbsp;&nbsp;&nbsp;    (4)
+
-
**  flhDC &rarr; flhB &rarr; lasI                        &nbsp;&nbsp;&nbsp;    (5)
+
-
**  fliA &rarr; fliL &rarr; Fluorescent Protein 1 (FP1)  &nbsp;&nbsp;&nbsp;    (6)
+
<center> >> [[Team:Paris/Modeling/hill_approach|Explanations and description]] </center>
-
**  fliA &rarr; flgA &rarr; Fluorescent Protein 2 (FP2)  &nbsp;&nbsp;&nbsp;    (7)
+
-
**  fliA &rarr; flhB &rarr; Fluorescent Protein 3 (FP3)  &nbsp;&nbsp;&nbsp;    (8)
+
-
**  fliA &rarr; flhB &rarr; lasI                          &nbsp;&nbsp;&nbsp;    (9)
+
-
<br>
+
* Specific experiments focused on finding relevant parameters have been designed and planned.
-
* <b>Model</b>
+
<center> >> [[Team:Paris/Modeling/estimation|Estimation]] </center>
-
The transcription hierarchy of the flagella genes have been studied in [1]. It resulted that the promoter activity of the seven class 2 operons, among which fLiL, flgA, flhB, may be mathematically described in that way :
+
-
<br>
+
-
[[Image:eqpromact.jpg|center]]
+
-
<br>
+
-
<br> where [X] denotes the effective protein-level activity.
+
-
<br>
+
-
<br>
+
-
For each operon, Shiraz Kalir and Uri Alon came up with numerical values of β and β', available in [1]. 
+
-
<br>
+
-
<br> Furthermore, the protein-level activity can be presented (for a more detailed presentation, see[4]) as
+
-
<br> [[Image:equation1.jpg|center]]
+
-
<br>Thus :
+
-
[[Image:fliA.jpg|center]]
+
-
[[Image:FP1.jpg|center]]
+
-
[[Image:FP2.jpg|center]]
+
-
[[Image:FP3.jpg|center]]
+
-
[[Image:lasI.jpg|center]]
+
-
* <b>Discution and Asumptions</b>
+
= Old but still usefull pages =
-
** As presented in the first equation below, there is a retroaction from fliA over fliA. In a first approach, we decided not to take into account this interaction, that is setting β'<sub>FliA</sub> to zero. In a second approach, we will consider this retroaction, and find a way to describe it mathematically.
+
-
** Another step consists in choosing the proper units, so as to have a homogenous set of parameters. In fact we need temporal parameters, but those given in the article are measured with the derivative of GFP over OD. Since we will work with a constant population density, we should only need to add a multiplicative factor, that we shall determine thanks to biological processes.
+
-
<br>
+
*[[Team:Paris/Modeling/Bibliography|Bibliographic References]]
-
 
+
*[[Team:Paris/Modeling/linear_approach|Preliminary approach]]
-
----
+
*[[Team:Paris/Modeling/Roadmap|Roadmap]]
-
* <b>Steps</b>
+
-
** lasI &rarr; HSL<sub>ext</sub> &nbsp;&nbsp;&nbsp; (10)
+
-
** lasI &rarr; HSL<sub>int</sub> &nbsp;&nbsp;&nbsp; (11)
+
-
 
+
-
* <b>Model</b>
+
-
[[Image:HSLint.jpg|center]]
+
-
[[Image:HSLext.jpg|center]]
+
-
 
+
-
* <b>Discution and Asumptions</b>
+
-
** Even though we have considered a single cell, we have decided to model both HSL inside and outside the cell. In a first approach, we have assumed that HSL could be modelized in the same fashion as AHL. The process was well detailed in [2] and we adapted their model to our present situation. Furthermore, a biologically coherent set of parameters is also given in this article. We decided then to use them as well for our simulation.
+
-
** As explained before, in this model, we only considered a single cell. The method to model the interactions with other cells is also described.
+
-
 
+
-
<br>
+
-
----
+
-
* <b>Steps</b>
+
-
** HSL<sub>int</sub> &rarr; tetR mRNA &nbsp;&nbsp;&nbsp; (12)
+
-
** tetR mRNA &rarr; tet R &nbsp;&nbsp;&nbsp; (13)
+
-
 
+
-
* <b>Model</b>
+
-
 
+
-
[[Image:tetR mRNA.jpg|center]]
+
-
[[Image:tetR.jpg|center]]
+
-
 
+
-
* <b>Discution and Asumptions</b>
+
-
** Here is the only case where we have not yet found out how we could manage to skip the mRNA step. However, this is acceptable since every parameter should be available. Then, (12) represents the binding and transcription steps ; (13) represents the translation step.
+
-
** The effect of HSL as well as the link between the mRNA activity and protein activity is well studied in [2]. We were the able to adapt the introduced situation to our system.
+
-
** The ρ parameter corresponds to maximal contribution to tetR transcription in presence of saturating amounts of HSL. Yet, to begin with, we have used the parameter given in the text. We shall then find a more accurate parameter corresponding to our system. The HSL activation is hereby chosen to follow Michalelis Menten kinetics since we assumed that there would be no cooperativity.
+
-
** Finally, we assumed that lasR would exist permanently in the medium in sufficient amounts, so that we needn't take it into account mathematically speaking.
+
-
 
+
-
 
+
-
<br>
+
-
----
+
-
* <b>Steps</b>
+
-
** tetR &rarr; flhDC &nbsp;&nbsp;&nbsp; (14)
+
-
 
+
-
* <b>Equations</b>
+
-
[[Image:flhDC.jpg|center]]
+
-
 
+
-
* <b>Discution and Asumptions</b>
+
-
** To determine this part of the model, we used [4] where the influence of tetR is described.
+
-
** Finally, to close the loop, we saw the necessity to use a little refinment. the influence of tetR over pTet (pTet being the promoter of flhDC) is a repression. To express the fact that a certain concentration is needed for the expression of flhDC, we used a step function. This choice was influenced by [3]. To fully define this function, two parameters are needed, so as to describe the threshold when the step has to occur, and the initial value of this step. We are facing two solutions. The more accurate procedure would consist in getting those values from the wet-lab. However, we may also look it up in the literature. Anyhow, as presented in the conclusive part, there is a really wide range of parameters that enable the system to oscillate.
+
-
** The key problem was that those two parameters were to be determined by the overall evolution of tetR, and the program cannot wait for tetR to be entirely determined to evaluate flhDC, since we have a loop. We then decided that these parameters could be obtanied by lab experiments.
+
-
 
+
-
<br>
+
-
 
+
-
===Parameters summary===
+
-
 
+
-
explication pourquoi on a pris des taux de degradation toujours égaux a 1
+
-
 
+
-
{|
+
-
|- style="background: #649CD7;"
+
-
! colspan="7" style="background: #649CD7;" | Parameter Table
+
-
|- style="background: #649CD7; text-align: center;"
+
-
| Parameter
+
-
| Meaning
+
-
| Value
+
-
| Unit
+
-
| Source
+
-
| Accuracy
+
-
| Comments
+
-
|- style="background: #dddddd;" 
+
-
| style="background: #D4E2EF;" |γ<sub>FliA</sub>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β<sub>FliA</sub>
+
-
| FlhDC activation coefficient
+
-
| 50
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β'<sub>FliA</sub>
+
-
| FliA activation coefficient
+
-
| 0
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |γ<sub>FP1</sub>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β<sub>FP1</sub>
+
-
| FlhDC activation coefficient
+
-
| 1200
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β'<sub>FP1</sub>
+
-
| FliA activation coefficient
+
-
| 250
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |γ<sub>FP2</sub>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β<sub>FP2</sub>
+
-
| FlhDC activation coefficient
+
-
| 150
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β'<sub>FP2</sub>
+
-
| FliA activation coefficient
+
-
| 300
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |γ<sub>FP3</sub>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β<sub>FP3</sub>
+
-
| FlhDC activation coefficient
+
-
| 100
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β'<sub>FP3</sub>
+
-
| FliA activation coefficient
+
-
| 350
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |γ<sub>lasI</sub>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β<sub>lasI</sub>
+
-
| FlhDC activation coefficient
+
-
| 100
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β'<sub>lasI</sub>
+
-
| FliA activation coefficient
+
-
| 350
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |γ<sub>HSL<sub>int</sub></sub> <br> <br>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |β<sub>HSL<sub>int</sub></sub> <br> <br>
+
-
| Production rate
+
-
| 0,01
+
-
|
+
-
| [2]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |δ<sub>int</sub>
+
-
| Diffusion rate in the cell
+
-
| 2
+
-
|
+
-
| [2]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |γ<sub>HSL<sub>ext</sub></sub> <br> <br>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |δ<sub>ext</sub>
+
-
| Diffusion rate out off the cell
+
-
| 2
+
-
|
+
-
| [2]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |γ<sub>tetR mRNA</sub>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |ρ
+
-
| Maximal contribution to tetR in the presence of saturating amounts of HSL
+
-
| 20
+
-
|
+
-
| [2]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |θ
+
-
| Ratio between mRNA and protein lifetimes
+
-
| 1
+
-
|
+
-
| [2]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |γ<sub>flhDC</sub>
+
-
| Degradation rate
+
-
| 1
+
-
|
+
-
| [1]
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |τ
+
-
| Threshold for tetR repression
+
-
| (20000)
+
-
|
+
-
| lab
+
-
|
+
-
|
+
-
|- style="background: #dddddd;"
+
-
| style="background: #D4E2EF;" |∆
+
-
| Promoter activity without the repressor
+
-
|(100)
+
-
|
+
-
| lab
+
-
|
+
-
|
+
-
|}
+
-
<br>
+
-
<br>
+
-
 
+
-
===Graph screenshots===
+
-
 
+
-
[[Image:OscillationsBio.jpg|center]]
+
-
 
+
-
===Corresponding codes===
+
-
 
+
-
*To change the parameters, you have to use this file, and refer to the previous array : [[Team:Paris/Modeling/modele|Model]]
+
-
*To launch the numerical simulation and plot the curves : [[Team:Paris/Modeling/plot|Plot the curves]]
+
-
*To compute the numerical simulations (does not provide any visual output, but necessary for the previous script) : [[Team:Paris/Modeling/simul|Numerical simulations]]
+
-
|}<br style="clear:both" />
+
-
 
+
-
==First Mathematical Approach==
+
-
 
+
-
===Introduction===
+
-
 
+
-
As a first approach, we had decided to take into account the binding to the promoters steps. Moreover, the transcription rates were expected to be Hill functions. Obvisouly, this modeling requires a huge number of parameters. To obtain them, we had planed to devise specific experiments (described below).
+
-
 
+
-
Nonetheless, after reading some more articles, we have decided to change several asumptions of the modeling choice. Therefore, we have devised a perhaps more biologically relevant framework (see above).
+
-
 
+
-
This part describes in detail the first approach and the codes that have been produced.
+
-
 
+
-
===First Approach===
+
-
 
+
-
As a first approximation, we have proposed a set of 5 ordinary differential equations, without taking into account the transcription step. Besides, we have had not introduced yet a synchronizaton module. Therefore, the repression of FlhDC is directly modeled by the presence of the 'Z3' gene (that is the last that is turned on).
+
-
 
+
-
In this framework, we have found parameters that have provided oscillations as well as a function that automatically detects whether the output of the ode system is oscillating. This has allowed to screen a little the parameters used, in order to evaluate the robustness of the system.
+
-
 
+
-
The methods employed are described there : [[Team:Paris/Modeling/first approach| First Approach]].
+
-
 
+
-
===More precise Bio-Mathematical Description===
+
-
 
+
-
After trying to obtain oscillations from a simple model, we have tried to described more precisely the studied system. Therefore, we have obtained the following formalism : [[Team:Paris/Modeling/description|Bio-Mathematical Description]].
+
-
 
+
-
===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 : [[Team:Paris/Modeling/Bibliography|Parameters we have to use]].
+
-
 
+
-
===Estimation of parameters===
+
-
 
+
-
Then, we will need many parameters to fully desribe the system according to the asumptions of the previous section. A natural way to have access to their value, after searching in the litterature, is to devise specific experiments. As a consequence of the characterization of the promoters activity, some Hill functions could be obtained.
+
-
 
+
-
Thus, we have described the experimental approach required : [[Team:Paris/Modeling/estimation|Estimation of the parameters]].
+
-
 
+
-
Nonetheless, as mentioned above, we have changed the way to model the biological reactions. As a result have stopped investigating in this way to focus on the [[Team:Paris/Modeling#An_Oscillatory_Biological_Model|An Oscillatory Biological Model]].
+
-
 
+
-
 
+
-
|}<br style="clear:both" />
+

Latest revision as of 04:46, 30 October 2008

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


Contents

Our train of thoughts...

We hereby propose different and complementary approaches to model the biological system. We found interesting to explain two of the paths that we chose to follow in order to understand and predict our system. It is important to note that both models aim at different goals in the process of understanding our system. Furthermore, we wished to describe our thought process, the way these models interact, their respective roles. An overall description of the way we model our biological system can be found below :

Read more !

BOB (Based On Bibliography) Approach

BOB.jpg

Our first approach is quite rough and simple but effective. The goal here was to guess the behavior of our Bacteri'OClock, considering the overall system. Since it is a preliminary approach, we could not yet fill the model with data from the wet lab. This is why our work is mainly based on a bibliographic work, which allows us to use parameters and data from scientific articles.

The key points of this approach:

  • Simplicity for itself is not that important. In fact, what we were looking for was understandability at first.
  • We used linear equations as much as possible: wherever it had been proved in a paper than an interaction could be efficiently modeled with a elementary expression, we kept it.
  • Too many parameters in a model mean less relevancy. In addition, the more parameters you have, the hardest it is to tune the system in order to have the behavior you are looking for.


Read more

APE (APE Parameters Estimation) Approach

APE.jpg

The second approach was motivated by our will to characterize our system in the most precise way. What is at stake here is to determine the "real parameters" that govern the dynamics of our system.

  • Each step is taken into account at a fundamental kinetic processes level or at a more global level by a function describing the complexation, which is a simple way to take into account multiple interactions and more especially cooperative binding.
>> Explanations and description
  • Specific experiments focused on finding relevant parameters have been designed and planned.
>> Estimation

Old but still usefull pages