The two-regime method for optimizing stochastic reaction-diffusion simulations.

Author: 

Flegg, M
Chapman, S
Erban, R

Publication Date: 

7 May 2012

Journal: 

J R Soc Interface

Last Updated: 

2019-08-18T13:43:03.78+01:00

Issue: 

70

Volume: 

9

DOI: 

10.1098/rsif.2011.0574

page: 

859-868

abstract: 

Spatial organization and noise play an important role in molecular systems biology. In recent years, a number of software packages have been developed for stochastic spatio-temporal simulation, ranging from detailed molecular-based approaches to less detailed compartment-based simulations. Compartment-based approaches yield quick and accurate mesoscopic results, but lack the level of detail that is characteristic of the computationally intensive molecular-based models. Often microscopic detail is only required in a small region (e.g. close to the cell membrane). Currently, the best way to achieve microscopic detail is to use a resource-intensive simulation over the whole domain. We develop the two-regime method (TRM) in which a molecular-based algorithm is used where desired and a compartment-based approach is used elsewhere. We present easy-to-implement coupling conditions which ensure that the TRM results have the same accuracy as a detailed molecular-based model in the whole simulation domain. Therefore, the TRM combines strengths of previously developed stochastic reaction-diffusion software to efficiently explore the behaviour of biological models. Illustrative examples and the mathematical justification of the TRM are also presented.

Symplectic id: 

193095

Download URL: 

Submitted to ORA: 

Not Submitted

Publication Type: 

Journal Article