Team:Paris/Modeling

From 2008.igem.org

Revision as of 15:22, 6 August 2008 by YannLC (Talk | contribs)

Contents

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
    • repressilator
    • comparison
  • Enhancing the system
    • Better FIFO behaviour
    • Other interactions to increase the robustness

Roadmap

If you want to have a look at our modeling notebook: Notebook

An Oscillatory Biological Model

Overall Model

Introduction


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.
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.
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

Musch of our inspiration comes from four articles to which we shall refer in the next subsections :

  • [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

  • Steps
    • flhDC → fliA     (1)
    • flhDC → fliL → Fluorescent Protein 1 (FP1)     (2)
    • flhDC → flgA → Fluorescent Protein 2 (FP2)     (3)
    • flhDC → flhB → Fluorescent Protein 3 (FP3)     (4)
    • flhDC → flhB → lasI     (5)
    • fliA → fliL → Fluorescent Protein 1 (FP1)     (6)
    • fliA → flgA → Fluorescent Protein 2 (FP2)     (7)
    • fliA → flhB → Fluorescent Protein 3 (FP3)     (8)
    • fliA → flhB → lasI     (9)


  • Model

The transcription hierarchy of the flagella genes have been studied in [1]. It resuted that the promoter activity of the seven class 2 operons, among which fLiL, flgA, flhB, may be mathematically described in that way :

Eqpromact.jpg



where [X] denotes the effective protein-level activity.

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


Equation1.jpg


Thus :

FliA.jpg
FP1.jpg
FP2.jpg
FP3.jpg
LasI.jpg
  • Discution and Asumptions
    • 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 β'FliA 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.



  • Steps
    • lasI → HSLext     (10)
    • lasI → HSLint     (11)
  • Model
HSLint.jpg
HSLext.jpg
  • Discution and Asumptions
    • 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] 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.



  • Steps
    • HSLint → tetR mRNA     (12)
    • tetR mRNA → tet R     (13)
  • Model
TetR mRNA.jpg
TetR.jpg
  • Discution and Asumptions
    • Here is the only case where we have not yet found out how we could manage to skip the mRNA step. Then, (12) represents the binding and transcription steps ; (13) represents the translation step.




Finally, the close the loop, we saw the necessity to use a little refinment. tetR influence over pTet, being the promoter of flhDC, is a repression. To express the non-instantaneous character of this interaction, we used a step function. To fully define this function, two parameters are needed, so as to describe the time 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 for it in the literature. Anyhow, as presented in the conclusive part, there is a really wide range of parameters that enable the system to oscillate.


FlhDC.jpg

Parameters summary

Graph screenshots

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.

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