Abstract:Background: Experiments in silico using stochastic reaction-diffusion models have emerged as an important tool in molecular systems biology. Designing computational software for such applications poses several challenges. Firstly, realistic lattice-based modeling for biological applications requires a consistent way of handling complex geometries, including curved inner- and outer boundaries. Secondly, spatiotemporal stochastic simulations are computationally expensive due to the fast time scales of individual reaction- and diffusion events when compared to the biological phenomena of actual interest. We therefore argue that simulation software needs to be both computationally efficient, employing sophisticated algorithms, yet in the same time flexible in order to meet present and future needs of increasingly complex biological modeling. Results: We have developed URDME, a flexible software framework for general stochastic reaction-transport modeling and simulation. URDME uses Unstructured triangular and tetrahedral meshes to resolve general geometries, and relies on the Reaction-Diffusion Master Equation formalism to model the processes under study. An interface to a mature geometry and mesh handling external software (Comsol Multiphysics) provides for a stable and interactive environment for model construction. The core simulation routines are logically separated from the model building interface and written in a low-level language for computational efficiency. The connection to the geometry handling software is realized via a Matlab interface which facilitates script computing, data management, and post-processing. For practitioners, the software therefore behaves much as an interactive Matlab toolbox. At the same time, it is possible to modify and extend URDME with newly developed simulation routines. Since the overall design effectively hides the complexity of managing the geometry and meshes, this means that newly developed methods may be tested in a realistic setting already at an early stage of development. Conclusions: In this paper we demonstrate, in a series of examples with high relevance to the molecular systems biology community, that the proposed software framework is a useful tool for both practitioners and developers of spatial stochastic simulation algorithms. Through the combined efforts of algorithm development and improved modeling accuracy, increasingly complex biological models become feasible to study through computational methods. URDME is freely available at http://www.urdme.org.

Larkin — A Software Package for the Numerical Simulation of LARge Systems Arising in Chemical Reaction KINetics

Diffusion Limited Kinetics In Reactive Systems

CERENA: ChEmical REaction Network Analyzer—A Toolbox for the Simulation and Analysis of Stochastic Chemical Kinetics

ReactionKinetics - A Mathematica package with applications II. Computational problems when building a reaction kinetics package

Predicting long timescale kinetics under variable experimental conditions with Kinetica.jl

ReactionMechanismSimulator.jl: A Modern Approach to Chemical Kinetic Mechanism Simulation and Analysis

Dynamic Adaptive Chemistry with Operator Splitting Schemes for Reactive Flow Simulations.

Predicting Long-Time-Scale Kinetics under Variable Experimental Conditions with Kinetica.jl

Efficient time stepping for reactive turbulent simulations with stiff chemistry

An overview of computational methods for chemical equilibrium and kinetic calculations for geochemical and reactive transport modeling

IMPLICIT-EXPLICIT RUNGE-KUTTA METHOD FOR COMBUSTION SIMULATION

Fourth-Order Runge–Kutta Schemes for Fluid Mechanics Applications

URDME: a modular framework for stochastic simulation of reaction-transport processes in complex geometries

Lattice Boltzmann Simulation of Diffusion-Convection Systems with Surface Chemical Reaction

A Gas-Kinetic Scheme for Multimaterial Flows and Its Application in Chemical Reactions

Catalytic flow with a coupled Finite Difference -- Lattice Boltzmann scheme

A Machine Learning Based Approach to Reaction Rate Estimation

Analytical methods for fitting integrated rate equations. A discontinuous assay

Dynamic Adaptive Chemistry for Turbulent Flame Simulations

ChemDyME: Kinetically Steered, Automated Mechanism Generation through Combined Molecular Dynamics and Master Equation Calculations

Arrhenius.jl: A Differentiable Combustion SimulationPackage