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.

Differentiable Programming of Reaction-Diffusion Patterns

Differentiable Programming of Chemical Reaction Networks

Programmable reaction-diffusion fronts

Bifurcation Analysis of Reaction Diffusion Systems on Arbitrary Surfaces

Bespoke Turing Systems

A Selection Criterion for Patterns in Reaction-Diffusion Systems

Network Controllability of Turing Reaction and Diffusion Model

Pattern control via Diffussion interaction

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

Isolating Patterns in Open Reaction-Diffusion Systems

Implementation of Lenia as a Reaction-Diffusion System

Solution landscape of reaction-diffusion systems reveals a nonlinear mechanism and spatial robustness of pattern formation

Adaptive POD-DEIM correction for Turing pattern approximation in reaction-diffusion PDE systems

Differentiable Simulations for Enhanced Sampling of Rare Events

DiffPD: Differentiable Projective Dynamics

Model-Based Reinforcement Learning Control of Reaction-Diffusion Problems

Spatial pattern formation in reaction–diffusion models: a computational approach

Parametric pattern selection in a reaction-diffusion model

Designing reaction-cross-diffusion systems with Turing and wave instabilities

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

A first-order computational algorithm for reaction-diffusion type equations via primal-dual hybrid gradient method