"
Team:NTNU Trondheim/Modeling
From 2011.igem.org
(→Steady State) |
(→Tau - leap Algorithm) |
||
| Line 80: | Line 80: | ||
This approach was used to write our one computer program for simulating any process in the cell using the tau - leap algorithm (see result for further information (coming soon) ). | This approach was used to write our one computer program for simulating any process in the cell using the tau - leap algorithm (see result for further information (coming soon) ). | ||
| + | |||
| + | ==Results== | ||
== Model Validation == | == Model Validation == | ||
Revision as of 19:55, 18 September 2011
Modeling
To describe and understand the biological reactions and processes as the bacteria turns red under stress, we developed multiple mathematical and statistical models. As the most basic model we used a system of ordinary differential equations (ODE). This is a fully deterministic model, describing the change in concentration for all molecules involved. The process for a bacteria to go from normal to glowing red is involved, and one can think of this as a stochastic process, that at each step the process can either succeed or fail with a given probability. This gives rise to a Bayesian model.
Combining these two models will result in a system of stochastic differential equations (SDE), which will be solved using numerical algorithms. As a last model we will explore the relationship between variants of stress and fluorescence intensity using (non-)conventional regression.
Model introduction
At the heart of the modeling lies biological consistency and data integration. The modeling will be focused on interpretation simplicity and data consistency. That is to develop models that can be easily interpreted by biologist and mathematicians, but also describe what is observed at the laboratory.
The two main ways to model biological systems are deterministic and stochastic. In this project we will attempt to approach the problems in both ways. Using a deterministic model, with fixed parameters, and a stochastic model to integrate data more dynamically.
Process description
As described in the introduction, ppGpp will repress the production of LacI and LacI represses the production of mCherry. These are the dominating processes. In short, when ppGpp are not present there will be little mCherry and when ppGpp are present the level of mCherry will be substantially higher.
In addition to these processes there are additional processes which might be of importance. First of all ppGpp affects the RNAP (RNA polymerase) and can therefore affect the production of mCherry as well as the production of LacI. As described in the introduction LacI is expected to be heavily downregulated by ppGpp. The production of mCherry can also in some cases lead to stress and therefore more ppGpp which in turn leads to more mCherry, in other words a positive feedback loop. These effects are assumed to be small, but might still affect the outcome.
The Models
Four basic models where constructed, which will be described below.
Systems of ODE
Models based on Ordinary Differential Equations (ODEs) are one of the most used methods of describing genetic circuits, while the different processes and reactions taking place in the cell are described by a set of coupled differential equations. This method might give both qualitative and quantitative information about the system and can therefore be very useful. It is however dependent on accurate kinetic reaction parameters and in some cases one also has to take into account the stochastic nature of genetic circuits. The equations are then solved either in Dizzy or Matlab.
Basic Model
Based on the model above the most important processes are:
The first process describes ppGpp attaching to the RNA polymerase, RNApA denotes the active part which is not repressed by ppGpp. RNApR is the number of repressed RNA polymerase molecules. This process can also be reversed. The transcription of mRNA is described in one step. D01 is the promoter determining the production of LacI mRNA and M1 is LacI mRNA. The LacI mRNA in turn leads to production of LacI transcription factor denoted as TF. Both mRNA and LacI is degraded as well. LacI will then inhibit mCherry:
Where D12 denotes promotors inhibited by LacI transcription factor. ξ is a measure of how strongly ppGpp affects mCherry directly, if 1 it is unaffected if 0 it is affected as strongly as LacI.
Steady State
Since it is difficult to find accurate parameters for all the processes involved and since it is often only the concentrations in steady state (long after the stress was first induced) that can be measured, all quantities involved were assumed to be constant. This assumption simplifies the system greatly and reduced the problem to these two equations:
Where Ci=αMiαpiRNApT/dMidpi, RNApT is the total number of RNA-polymerase molecules in the cell and Ki=ki/qi. To check that the equations for steady state are correct a comparison was made between the numerical solution in Dizzy and the analytical solution for different levels of ppGpp. The comparison shown in the figure below shows that there is excellent agreement.
Stochastic Differential equations
In the cell, all movement of the molecules is random, this gives rise to the stochasticity observed in gene expression experiment. In this section we outline two methods where the randomness is accounted for; the Gillespie algorithm and the approximation τ - leap algorithm.
When the number of participating molecules are low (which do happen in the cell), then stochasticity really matters. If there are many molecules, then the behavior of the reactions goes "smooth" and looks like an ODE. However when to number of reactants gets small, the number of reactions in a small time frame varies, this gives rise to the irregularity seen in the time series for the number of created products in the cell.
Gillespie
The Gillespie algorithm is a model that simulates the number of reactions and the time between the reactions exactly (under some assumptions). The most important assumptions is that the concentration of the reactants is distributed uniformly in the cell, the other important assumption is that the time between the reactions occur is markovian, that is it has no memory of how long the previous time step was.
An improved version of the original algorithm, the Gibson-Bruck algorithm, was used with Dizzy.
Tau - leap Algorithm
Instead of calculating the time between each reaction, we can fix a time frame and estimate how many reactions occur in that time frame.
It can be shown that setting the production and destruction rate from the ODE in as a rate parameter in an Poisson distribution, will give the same result. This is a nice way of reconfiguring the tau - leap such that it gives a clear and precise relation to the ODE models, and mathematical consistency.
This approach was used to write our one computer program for simulating any process in the cell using the tau - leap algorithm (see result for further information (coming soon) ).
Results
Model Validation
To see if the results agree with the theoretical model is vital in all parts of science. In this section we will give a short reasoning of why our models seems to fit the observed data. Here we will consider our model in two distinct parts, one deterministic model and one stochastic model, both of equally importance.
Stochastic -- Flowcytometry
Deterministic -- mass spectrometry
References
