Quick search Find article
Quick search
Find article

Stochastic modelling of reaction–diffusion processes: algorithms for bimolecular reactions

Radek Erban and S Jonathan Chapman

Show affiliations


Several stochastic simulation algorithms (SSAs) have recently been proposed for modelling reaction–diffusion processes in cellular and molecular biology. In this paper, two commonly used SSAs are studied. The first SSA is an on-lattice model described by the reaction–diffusion master equation. The second SSA is an off-lattice model based on the simulation of Brownian motion of individual molecules and their reactive collisions. In both cases, it is shown that the commonly used implementation of bimolecular reactions (i.e. the reactions of the form A + BC or A + AC) might lead to incorrect results. Improvements of both SSAs are suggested which overcome the difficulties highlighted. In particular, a formula is presented for the smallest possible compartment size (lattice spacing) which can be correctly implemented in the first model. This implementation uses a new formula for the rate of bimolecular reactions per compartment (lattice site).


PACS

87.15.R- Reactions and kinetics

87.17.-d Cell processes

87.15.A- Theory, modeling, and computer simulation

02.50.Ey Stochastic processes

Subjects

Computational physics

Biological physics

Dates

Issue 4 (December 2009)

Received 6 March 2009, accepted for publication 3 July 2009

Published 21 August 2009

 
Image from Stochastic modelling of reaction-diffusion processes: algorithms for bimolecular reactions


View by subject




Export








Please login to access our web services, or create an account if you don't yet have one.

You must have cookies enabled in your web browser to be able to login.

Username
Password

Forgotten your password? Get a new one here.