Abstract:Biochemical reactions often occur at low copy numbers but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. While particle-based models guarantee the level of detail necessary to accurately describe the microscopic dynamics at very low copy numbers, the algorithms used to simulate them typically imply trade-offs between computational efficiency and biochemical accuracy. eGFRD (enhanced Green's Function Reaction Dynamics) is an exact algorithm that evades such trade-offs by partitioning the N-particle system into M ≤ N analytically tractable one- and two-particle systems; the analytical solutions (Green's functions) then are used to implement an event-driven particle-based scheme that allows particles to make large jumps in time and space while retaining access to their state variables at arbitrary simulation times. Here we present "eGFRD2," a new eGFRD version that implements the principle of eGFRD in all dimensions, thus enabling efficient particle-based simulation of biochemical reaction-diffusion processes in the 3D cytoplasm, on 2D planes representing membranes, and on 1D elongated cylinders representative of, e.g., cytoskeletal tracks or DNA; in 1D, it also incorporates convective motion used to model active transport. We find that, for low particle densities, eGFRD2 is up to 6 orders of magnitude faster than conventional Brownian dynamics. We exemplify the capabilities of eGFRD2 by simulating an idealized model of Pom1 gradient formation, which involves 3D diffusion, active transport on microtubules, and autophosphorylation on the membrane, confirming recent experimental and theoretical results on this system to hold under genuinely stochastic conditions.

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

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

pSpatiocyte: A Parallel Stochastic Method for Particle Reaction-Diffusion Systems

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

Grid-free stochastic simulations of reaction-diffusion processes at cell-cell contacts

Constant-complexity Stochastic Simulation Algorithm with Optimal Binning

Green's Function Reaction Dynamics: a new approach to simulate biochemical networks at the particle level and in time and space

eGFRD in all dimensions

A practical guide to stochastic simulations of reaction-diffusion processes

Fine-Grained Parallel And Distributed Spatial Stochastic Simulation Of Biological Reactions

Smoldyn on Graphics Processing Units: Massively Parallel Brownian Dynamics Simulations

Hybrid methods in reaction-diffusion equations

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

Realistic boundary conditions for stochastic simulations of reaction-diffusion processes

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

A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter

Optimisation of Simulations of Stochastic Processes by Removal of Opposing Reactions

First-Passage Kinetic Monte Carlo method

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

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

Single molecule simulations in complex geometries with embedded dynamic one-dimensional structures