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.

Reaction Kinetics in a Tight Spot

Reaction kinetics in intracellular environments with macromolecular crowding: simulations and rate laws

Strong defocusing of molecular reaction times results from an interplay of geometry and reaction control

Multiscale Analysis of Reaction Networks

Universal kinetics of imperfect reactions in confinement

Fast reactions with non-interacting species in stochastic reaction networks

Effective Kinetics of Chemical Reaction Networks

Diffusion Limited Kinetics In Reactive Systems

Anomalous Kinetics of a Multi-Species Reaction-Diffusion System: Effect of Random Velocity Fluctuations

Diffusion-controlled reactions with non-Markovian binding/unbinding kinetics

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

Nonequilibrium reaction rate theory: Formulation and implementation within the hierarchical equations of motion approach

Reaction kinetics in open reactors and serial transfers between closed reactors

Reaction-diffusion dynamics: confrontation between theory and experiment in a microfluidic reactor

Enzyme kinetics simulation at the scale of individual particles

Smoluchowski Reaction Kinetics for Reactions of Any Order

Learning heterogenous reaction rates from stochastic simulations

Universal Slow Dynamics of Chemical Reaction Networks

On The Effect Of Size On The Kinetics Of Reactions In Solutions

Imperfect Narrow Escape problem