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.

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

Modeling of chemical reaction systems with detailed balance using gradient structures

A Convergent Reaction-Diffusion Master Equation

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

On the Reaction Diffusion Master Equation in the Microscopic Limit

Reaction rates for mesoscopic reaction-diffusion kinetics

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

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

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

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

Drift in Diffusion Gradients.

On microscopic origins of generalized gradient structures

Gradient flows on metric graphs with reservoirs: Microscopic derivation and multiscale limits

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

Multiscale Mechanisms of Reaction-Diffusion Process in Electrode Systems: A Classical Density Functional Study

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

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

Chemotaxis and Reactions in Anomalous Diffusion Dynamics

A new flow-kinematics-based model for time-dependent effective dispersion in mixing-limited reactions

Gamma-convergence of a gradient-flow structure to a non-gradient-flow structure

A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter