Team:UNIPV-Pavia/Modeling

From 2008.igem.org

(Difference between revisions)
Line 139: Line 139:
According to the equations described so far, every protein and mRNA is a state variable for the model. Mux and Demux models are described in the next paragraphs.
According to the equations described so far, every protein and mRNA is a state variable for the model. Mux and Demux models are described in the next paragraphs.
<br>
<br>
-
We considered full dimensional models, in which transcription and translation processes are separated. In many contexts, anyway, we can write only one differential equation to describe a "regulated protein production", in which trascription process is hidden. We preferred to consider the transcription dynamic because our aim in this section is to provide a model for future parameter identification and so
+
We considered full ODE models, in which transcription and translation processes are separated. In many contexts, anyway, we can write only one differential equation to describe a "regulated protein production", in which trascription process is hidden. In this section we preferred to be general and we considered the transcription dynamic because our aim is to provide a mathematical model for future parameter identification. We also decided to neglect AHL synthesis dynamic, considering luxI protein a "direct" activator for luxR.
-
 
+
-
dimensionless equations, in which parameters have been calibrated qualitatively, also taking into account values found in literature, to study system trends and to build up a scaffold for future parameter identification.
+
<br>
<br>
In the following paragraphs, the complete ODE systems and the Simulink block diagrams will be reported.
In the following paragraphs, the complete ODE systems and the Simulink block diagrams will be reported.
Line 163: Line 161:
== '''Simulations''' ==
== '''Simulations''' ==
 +
The written models have been simulated in Matlab and Simulink, in order to perform a preliminary heuristic study of the system behavior and to search for critical parameters. For simulations, we used dimensionless ODEs to avoid numeric errors due to small parameter values, so results are reported in arbitrary units. We defined ODE systems in Simulink (block diagrams showed above) and then solved them using ode15s, a multistep ODE solver for stiff problems.
 +
<br>
 +
Parameters have been calibrated considering the values reported in literature and in some other iGEM wiki pages. In some cases values were missing, so we fixed these missing values at 1*10^om, where "om" is the median order of magnitude of the parameters having the same meaning (e.g. protein degradation rates).
 +
In other cases we found different values for the same parameter, so we computed the mean value.
 +
<br>
 +
An important feature that will have to be studied is systems' basic activity. In other words, can genetic Mux discriminate the output in response to a low input signal and the output in response to a high input signal? Moreover, does Demux have a non negligible baseline activity that turns on the output channel that should be unactive? One of the parameters that plays a very important role in basic activities is the leakage factor ("a"). Its future estimation for all
=== '''Mux'''===
=== '''Mux'''===
=== '''Demux'''===
=== '''Demux'''===

Revision as of 11:02, 27 October 2008


Home.jpg Home Unipv logo.jpg The Team And.jpg The Project Safety.jpg Biological Safety Dna.png Parts Submitted to the Registry
Laptop.jpg Dry Lab Pipette.jpg Wet Lab Math.gif Modeling Note.jpg Protocols Notebook.gif Activity Notebook


Contents

Mathematical modeling page

In this section we explain two dynamic models that can be used to describe the gene networks in our project. After a brief overview about the motivation of a mathematical model, we will illustrate the general formulas we used, we will show the complete ODE models for Mux and Demux gene networks, and then we will report results of some simulations performed with Matlab and Simulink.

Why writing a mathematical model?

The purposes of writing mathematical models for gene networks can be:

  • Prediction: a good and well identificated model can be used in simulations to predict real system behavior. In particular we could be interested in system output in response to never seen inputs. In this way, the system can be tested 'in silico', without performing real experiments 'in vitro' or 'in vivo'.
  • Parameter identification: we already wrote that it is very important to estimate all the parameters involved in the model, in order to perform realistic simulations. Another goal that can be reached with parameter identification is 'network summarization', in fact estimated parameters can be used as 'behavior indexes' for the network (or a part of it). These indexes can be very useful for synthetic biologists to choose and compare BioBrick standard parts for genetic circuits design, just like electronic engineers choose, for example, a Zener diode, knowing its Zener voltage.



Equations for gene networks

In this paragraph, mathematical modeling for gene and protein interactions will be described.

Binding of a ligand to a molecule: Hill equation

The fraction of a molecule saturated by a ligand ( = the probability that the molecule is bound to a ligand) can be expressed as a function of ligand concentration using Hill equation:

Pv hill.jpg
  • Prob: the probability that the molecule is bound to the ligand
  • [T]: concentration of the ligand
  • K50: dissociation constant. It is the concentration producing a probability of 0.5
  • n: Hill coefficient. It describes binding cooperativity


Regulated transcription

Transcription of a mRNA molecule can be regulated by transcription factors in active form. These factors can be activators or inhibitors: they respectively increase and decrease the probability that RNA polymerase binds the promoter. This probability must be function of active transcription factor concentration and can be modeled using Hill equation. If factor in active form activates transcription, we can write:

Pv acttranscr.jpg
  • Prob: probability that the gene is transcripted
  • [T]: concentration of the transcription activator in active form
  • K50: dissociation constant. It is the concentration producing a probability of 0.5
  • n: Hill coefficient (n>0)

while, if factor in active form inhibits transcription, we are interested in unbound promoter and so we can write:

Pv inhtranscr.jpg
  • Prob: probability that the gene is transcripted
  • [T]: concentration of the transcription inhibitor in active form
  • K50: dissociation constant. It is the concentration producing a probability of 0.5
  • n: Hill coefficient (n>0)

The equations show that:

  • In activation formula, if [T]=0 the trascription probability is Prob=0, while the maximum probability, Prob=1, is reached asymptotically for [T]->Inf
  • In inhibition formula, if [T]=0 the trascription probability is Prob=1, while the minimum probability, Prob=0, is reached asymptotically for [T]->Inf

These formulas can be included in a differential equation that can describe the dynamic behavior of a mRNA molecule as a function of the factor that regulates its transcription.

ODE for activated transcription

Pv ode acttranscr.jpg
  • Vmax: maximum transcription rate
  • [m]: concentration of mRNA molecule
  • [T]: concentration of the transcription activator in active form
  • K50: dissociation constant
  • n: Hill coefficient (n>0)
  • delta: mRNA degradation constant
  • a: leakage factor that can modellize promoter basic activity. It is a percentage of Vmax


ODE for inhibited transcription

Pv ode inhtranscr.jpg
  • Vmax: maximum transcription rate
  • [m]: concentration of mRNA molecule
  • [T]: concentration of the transcription inhibitor in active form
  • K50: dissociation constant
  • n: Hill coefficient (n>0)
  • delta: mRNA degradation rate
  • a: leakage factor that can modellize promoter basic activity. It is a percentage of Vmax


ODE for constitutive transcription

Pv odecontranscr.jpg
  • c: constant transcription rate
  • [m]: concentration of mRNA molecule
  • delta: mRNA degradation rate


Transcription factors in active form

Formulas described so far are dependent on the concentration of transcription factor in active form. But how can we calculate the amount of transcription factor in active form? Assuming that the transcription factor can be activated or inhibited by a inducer factor (that can be an exogenous input), we can use Hill equation again to calculate the concentration of bound and unbound transcription factor as a function of the inducer factor:

Pv inducers.jpg
  • [I]: concentration of the inducer
  • [Tt]: total concentration of the transcription factor
  • [Tu]: concentration of unbound transcription factor
  • [Tb]: concentration of bound transcription factor
  • K50: dissociation constant
  • n: Hill coefficient (n>0)

If active form is bound transcripion factor, we are interested in Tb, while if active form is unbound transcription factor, we are interested in Tu.

ODE for protein production

We have to modellize translation of mRNA molecules. We chose to describe it as a differential equation in which protein production and protein degradation are linear:

Pv protprod.jpg
  • [P]: concentration of protein
  • [m]: concentration of mRNA molecule
  • alpha: translation rate
  • beta protein degradation rate


Our model

We considered Tetracycline, IPTG, Arabinose and GFP respectively as CH0, CH1, SEL and OUT signals for Mux, while we considered IPTG, Tetracycline, RFP and GFP respectively as IN, SEL, OUT0 and OUT1 signals for Demux. For all these signals, logic 1 corresponds to the presence of the molecule.
According to the equations described so far, every protein and mRNA is a state variable for the model. Mux and Demux models are described in the next paragraphs.
We considered full ODE models, in which transcription and translation processes are separated. In many contexts, anyway, we can write only one differential equation to describe a "regulated protein production", in which trascription process is hidden. In this section we preferred to be general and we considered the transcription dynamic because our aim is to provide a mathematical model for future parameter identification. We also decided to neglect AHL synthesis dynamic, considering luxI protein a "direct" activator for luxR.
In the following paragraphs, the complete ODE systems and the Simulink block diagrams will be reported.

Mux

200%
200%

Demux

200%
200%


Simulations

The written models have been simulated in Matlab and Simulink, in order to perform a preliminary heuristic study of the system behavior and to search for critical parameters. For simulations, we used dimensionless ODEs to avoid numeric errors due to small parameter values, so results are reported in arbitrary units. We defined ODE systems in Simulink (block diagrams showed above) and then solved them using ode15s, a multistep ODE solver for stiff problems.
Parameters have been calibrated considering the values reported in literature and in some other iGEM wiki pages. In some cases values were missing, so we fixed these missing values at 1*10^om, where "om" is the median order of magnitude of the parameters having the same meaning (e.g. protein degradation rates). In other cases we found different values for the same parameter, so we computed the mean value.
An important feature that will have to be studied is systems' basic activity. In other words, can genetic Mux discriminate the output in response to a low input signal and the output in response to a high input signal? Moreover, does Demux have a non negligible baseline activity that turns on the output channel that should be unactive? One of the parameters that plays a very important role in basic activities is the leakage factor ("a"). Its future estimation for all

Mux

Demux


References

[1] Modelli multiscala per studi di proteomica e genomica funzionale, S. Cavalcanti, S. Furini, E. Giordano. "Genomica e proteomica computazionale", GNB Book, Publisher: Pàtron, 2007

[2] ETH Zurich iGEM 2007 Wiki - Modeling basics page