Goal

Stochastic modeling is a technique of predicting outcomes that takes into account a certain degree of randomness or unpredictability. In a stochastic modeling, a small amount of randomness is added at each time step of the simulation.
It is the probabilistic counterpart to a deterministic process.

Why

The deterministic model, who studies a population of bacteria, considers continuous concentrations of molecules. However, in a single bacteria, the quantity of the different proteins is of the order of 100, and the concentrations take thus discrete values. These values depend on events (production, degradation) which are hard to predict, and must therefore be approached in terms of probability of occurence (or preopensity).

How

A stochastic model incorporates random variations to predict future conditions and to see what they might be like.
To introduce that randomness we use a new function : propensities.

As an example, to illustrate this we consider 4 possible reactions. Each reaction has a probability to happen in the next time step.

We randomly chose the next reaction regarding the propensities.

When we run a stochastic simulation once, we get a trajectory. This graph represents the random evolution of the variables. Because of the randomness, if we run the script another time we will get a different trajectory. That is why to be able to interpret the results we have to run the cripts hundreds or thousands of times.

Instead of describing a process which can only evolve in one way, in a stochastic or random process there is some indeterminacy : even if the initial condition is known, there are several directions in which the process may evolve.

To model that randomness we use a Gillepsie algorithm or Stoachastic Simulation Algorithm (SSA).

Traditional continuous and deterministic biochemical rate equations do not accurately predict cellular reactions since they rely on random reactions that require the interactions of millions of molecules. In contrast, the Gillespie algorithm allows a discrete and stochastic simulation of a system with few reactants because every molecule is explicitly simulated. When simulated, a Gillespie realization represents a random walk of the entire system.

We will use the following notations :

Each reaction Rj is characterized mathematically by two quantities :

Once we have presented the theory behind the stochastic approach, let us have a look at the algorithm.

The algorithm comprises 5 steps.

Conerning the parameters, we use the same as in the ODE modeling but we have to make a conversion to remove the dependence in volume and mol. Thus, in most cases, we need to multiply those parameters by the Avogadro number and by the cell's volume.

You can find our script for the stochastic modeling here.

Team:Grenoble/Modeling/Amplification/Stochastic/what

From 2012.igem.org

Goal

Why

How