Abstract:Biochemical networks play a crucial role in biological systems, implementing a broad range of vital functions. They normally operate at low copy numbers and in spatial settings, but this is often ignored and well-stirred conditions are assumed. Yet, it is increasingly becoming clear that even microscopic spatial inhomogeneities oftentimes can induce significant differences on the macroscopic level. Since experimental observation of single-molecule behavior is extremely challenging, theoretical modeling of biochemical reactions on the single-particle level is an important tool for understanding spatial effects in biochemical systems. While purely analytical models quickly become intractable here, spatial-stochastic simulations can capture a wide range of biochemical processes with the necessary levels of detail. Here we discuss different techniques for spatial-stochastic simulation of reaction-diffusion systems, and explain important precautions required to make them biochemically accurate and efficient. We illustrate non-negligible accuracy issues arising even in the most simple approaches to biochemical simulation, and present methods to deal with them. We first explain how Brownian Dynamics, a widely used particle-based diffusion simulation technique with fixed propagation time, can be adapted to incorporate chemical reactions, and portray a range of schemes that elaborate on this idea. We then introduce event-driven spatial-stochastic simulation methods, in which system updates are performed asynchronously with situation-dependent, varying time steps; here we particularly focus on eGFRD, a computationally efficient particle-based algorithm that makes use of analytical functions to accurately sample interparticle reactions and diffusive motion with large jumps in time and space. We end by briefly presenting recent developments in the field of spatial-stochastic biochemical simulation.

Macroscopic limit for stochastic chemical reactions involving diffusion and spatial heterogeneity

Diffusion-limited reactions in dynamic heterogeneous media

Computing macroscopic reaction rates in reaction-diffusion systems using Monte Carlo simulations

On the Reaction Diffusion Master Equation in the Microscopic Limit

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

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

Drift in Diffusion Gradients.

A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter

Homogenization results for a coupled system of reaction-diffusion equations

Revisit of macroscopic dynamics for some non-equilibrium chemical reactions from a Hamiltonian viewpoint

Exponential time decay of solutions to reaction-cross-diffusion systems of Maxwell-Stefan type

Fast Diffusion Limit for Reaction-Diffusion Systems with Stochastic Neumann Boundary Conditions

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

Macroscopic dynamics of biological cells interacting via chemotaxis and direct contact

A Comparison of Bimolecular Reaction Models for Stochastic Reaction Diffusion Systems

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

Fractional Chemotaxis Diffusion Equations

Macroscopic Stochastic Thermodynamics

Modeling of chemical reaction systems with detailed balance using gradient structures

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

Chemical Continuous Time Random Walks under Anomalous Diffusion