Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art

Abstract

Stochasticity is a key characteristic of intracellular processes such as gene regulation and chemical signalling. Therefore, characterising stochastic effects in biochemical systems is essential to understand the complex dynamics of living things. Mathematical idealisations of biochemically reacting systems must be able to capture stochastic phenomena. While robust theory exists to describe such stochastic models, the computational challenges in exploring these models can be a significant burden in practice since realistic models are analytically intractable. Determining the expected behaviour and variability of a stochastic biochemical reaction network requires many probabilistic simulations of its evolution. Using a biochemical reaction network model to assist in the interpretation of time course data from a biological experiment is an even greater challenge due to the intractability of the likelihood function for determining observation probabilities. These computational challenges have been subjects of active research for over four decades. In this review, we present an accessible discussion of the major historical developments and state-of-the-art computational techniques relevant to simulation and inference problems for stochastic biochemical reaction network models. Detailed algorithms for particularly important methods are described and complemented with MATLAB implementations. As a result, this review provides a practical and accessible introduction to computational methods for stochastic models within the life sciences community.

References in corpus (18)

• Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems
• The pseudo-marginal approach for efficient Monte Carlo computations
• Adaptive approximate Bayesian computation
• A modified Next Reaction Method for simulating chemical systems with time dependent propensities and delays
• Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
• Stochastic modelling of reaction-diffusion processes: algorithms for bimolecular reactions
• Approximate Bayesian computation (ABC) gives exact results under the assumption of model error
• Multilevel Sequential Monte Carlo Samplers
• Stochastic Simulation of Biomolecular Networks in Dynamic Environments
• Error analysis of tau-leap simulation methods
• Optimal quantification of contact inhibition in cell populations
• An adaptive multi-level simulation algorithm for stochastic biological systems
• Multilevel rejection sampling for approximate Bayesian computation
• The breakdown of the reaction-diffusion master equation with non-elementary rates
• Extending the multi-level method for the simulation of stochastic biological systems
• Multilevel Markov Chain Monte Carlo Method for High-Contrast Single-Phase Flow Problems
• Models, measurement and inference in epithelial tissue dynamics
• Approximate maximum likelihood estimation using data-cloning ABC

Cited by in corpus (9)

• A practical guide to pseudo-marginal methods for computational inference in systems biology
• Rapid Bayesian inference for expensive stochastic models
• Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes
• A new approach to simulating stochastic delayed systems
• Vector operations for accelerating expensive Bayesian computations -- a tutorial guide
• Identifiability of Stochastically Modelled Reaction Networks
• Asymptotic expansion approximation for spatial structure arising from directionally biased movement
• Inferring birth versus death dynamics for ecological interactions in stochastic heterogeneous populations
• Stochastic Model of siRNA Endosomal Escape Mediated by Fusogenic Peptides in OVCAR-3