Abstract:BackgroundIt has been established that stochastic effects play an important role in spatio-temporal biochemical networks. A popular method of representing such stochastic systems is the Reaction Diffusion Master Equation (RDME). However, simulating sample paths from the RDME can be computationally expensive, particularly at large populations. Here we investigate an uncommon, but much faster alternative: the Spatial Chemical Langevin Equation (SCLE).MethodsWe investigate moment equations and correlation functions analytically, then we compare sample paths and moments of the SCLE to the RDME and associated deterministic solutions. Sample paths are generated computationally by the Next Subvolume method (RDME) and the Euler-Maruyama method (SCLE), while a deterministic solution is obtained with an Euler method. We consider the Gray-Scott model, a well-known pattern generating system, and a predator–prey system with spatially inhomogeneous parameters as sample applications.ResultsFor linear reaction networks, it is well known that the first order moments of all three approaches match, that the RDME and SCLE match to the second moment, and that all approaches diverge at third order moments. For non-linear reaction networks, differential equations governing moments do not form a closed system, but a general moment equation can be compared term wise. All approaches match at the leading order, and the RDME and SCLE match at the second leading order. As expected, the SCLE captures many dynamics of the RDME where deterministic methods fail to represent them. However, areas of the parameter space in the Gray-Scott model exist where either the SCLE and RDME give qualitatively different predictions, or the RDME predicts patterns, while the SCLE does not.ConclusionsThe SCLE provides a fast alternative to existing methods for simulation of spatial stochastic biochemical networks, capturing many aspects of dynamics represented by the RDME. This becomes very useful in search of quantitative parameters yielding desired qualitative solutions. However, there exist parameter sets where both the qualitative and quantitative behaviour of the SCLE can differ when compared to the RDME, so care should be taken in its use for applications demanding greater accuracy.

On the Reaction Diffusion Master Equation in the Microscopic Limit

Reaction rates for mesoscopic reaction-diffusion kinetics

A Convergent Reaction-Diffusion Master Equation

Reaction-diffusion kinetics on lattice at the microscopic scale

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

An Unstructured Mesh Reaction-Drift-Diffusion Master Equation with Reversible Reactions

A route to the hydrodynamic limit of a reaction-diffusion master equation using gradient structures

Macroscopic limit for stochastic chemical reactions involving diffusion and spatial heterogeneity

Multi-grid reaction-diffusion master equation: applications to morphogen gradient modelling

A Comparison of Bimolecular Reaction Models for Stochastic Reaction Diffusion Systems

Surface reaction-diffusion kinetics on lattice at the microscopic scale

The Spatial Chemical Langevin Equation and Reaction Diffusion Master Equations: moments and qualitative solutions

Grand canonical diffusion-influenced reactions: a stochastic theory with applications to multiscale reaction-diffusion simulations

Reaction-controlled diffusion: Monte Carlo simulations

A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter

A microscopic approach to nonlinear Reaction-Diffusion: the case of morphogen gradient formation

Reaction-diffusion and reaction-subdiffusion equations on arbitrarily evolving domains

Mean field limits of particle-based stochastic reaction-drift-diffusion models

Chemical diffusion master equation: formulations of reaction--diffusion processes on the molecular level

Diffusion approximations to the chemical master equation only have a consistent stochastic thermodynamics at chemical equilibrium

Effects of different discretisations of the Laplacian upon stochastic simulations of reaction-diffusion systems on both static and growing domains