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.

Stochastic theory of large-scale enzyme-reaction networks: Finite copy number corrections to rate equation models

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

Investigating the robustness of the classical enzyme kinetic equations in small intracellular compartments

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Stochastic quasi-steady state approximations for asymptotic solutions of the chemical master equation

Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion

Multiscale Analysis of Reaction Networks

A Stein's Method Approach to the Linear Noise Approximation for Stationary Distributions of Chemical Reaction Networks

Variance-corrected Michaelis-Menten equation predicts transient rates of single-enzyme reactions and response times in bacterial gene-regulation

Sensitivity of chemical reaction networks: a structural approach 3. Regular multimolecular systems

Approximation and inference methods for stochastic biochemical kinetics - a tutorial review

Stochastic simulation of biochemical systems with randomly fluctuating rate constants

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

Model reduction for stochastic chemical systems with abundant species

Michaelis-Menten dynamics in protein subnetworks

Coarse-graining stochastic biochemical networks: quasi-stationary approximation and fast simulations using a stochastic path integral technique

Noise and critical phenomena in biochemical signaling cycles at small molecule numbers

Reduction for stochastic biochemical reaction networks with multiscale conservations

Fast reactions with non-interacting species in stochastic reaction networks

Quasiequilibrium approximation of fast reaction kinetics in stochastic biochemical systems.

A Generic Rate Equation for modeling Enzymatic Reactions under Living Conditions