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.

A new class of highly efficient exact stochastic simulation algorithms for chemical reaction networks

Exact rule-based stochastic simulations for systems with unlimited number of molecular species

Coarse-Grained Stochastic Particle-based Reaction-Diffusion Simulation Algorithm

Constant-complexity Stochastic Simulation Algorithm with Optimal Binning

Direct ab initio MD simulation of silver ion diffusion in chalcogenide glasses

Exact Parallelization of the Stochastic Simulation Algorithm for Scalable Simulation of Large Biochemical Networks

Nested stochastic simulation algorithms for chemical kinetic systems with multiple time scales

The regime-conversion method: a hybrid technique for simulating well-mixed chemical reaction networks

A hybrid tau-leap for simulating chemical kinetics with applications to parameter estimation

Optimisation of Simulations of Stochastic Processes by Removal of Opposing Reactions

Robust Stochastic Chemical Reaction Networks and Bounded Tau-Leaping

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

Differentiable Simulations for Enhanced Sampling of Rare Events

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

Reaction Sampling and Reactivity Prediction Using the Stochastic Surface Walking Method

A general method for numerically simulating the stochastic time evolution of coupled chemical reactions

A low-rank complexity reduction algorithm for the high-dimensional kinetic chemical master equation

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

Reaction Prediction Via Atomistic Simulation: from Quantum Mechanics to Machine Learning

Automated Importance Sampling via Optimal Control for Stochastic Reaction Networks: A Markovian Projection-based Approach