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.

RTF: An R package for modelling time course data

Dynamic modelling of signalling pathways when ODEs are not feasible

Impulse Model-Based Differential Expression Analysis of Time Course Sequencing Data

Data2Dynamics: a Modeling Environment Tailored to Parameter Estimation in Dynamical Systems.

Inferring Population Dynamics from Single-Cell RNA-sequencing Time Series Data.

TIDAL: Tool to Implement Developmental Analysis of Longitudinal data

Beyond pseudotime: Following T-cell maturation in single-cell RNAseq time series

Diffusion model analysis with MATLAB: A DMAT primer

Parameter Estimation for Dynamical Systems with Discrete Events and Logical Operations

Automated analysis of activity, sleep, and rhythmic behaviour in various animal species with the Rtivity software

Alignment of time-resolved data from high throughput experiments

Cost-effectiveness analysis of treatments for chronic disease: using R to incorporate time dependency of treatment response

pulseR: Versatile computational analysis of RNA turnover from metabolic labeling experiments

PyDTS: A Python Package for Discrete-Time Survival (Regularized) Regression with Competing Risks

Classification using functional data analysis for temporal gene expression data

Deriving quantitative dynamics information for proteins and RNAs using ROTDIF with a graphical user interface

Bayesian modelling of time series data (BayModTS)-a FAIR workflow to process sparse and highly variable data

Inferring the perturbation time from biological time course data

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

Fitting the ratcliff diffusion model to experimental data

pulseTD: RNA life cycle dynamics analysis based on pulse model of 4sU-seq time course sequencing data