Stochastic Simulation of Biomolecular Networks in Dynamic Environments

Abstract

Simulation of biomolecular networks is now indispensable for studying biological systems, from small reaction networks to large ensembles of cells. Here we present a novel approach for stochastic simulation of networks embedded in the dynamic environment of the cell and its surroundings. We thus sample trajectories of the stochastic process described by the chemical master equation with time-varying propensities. A comparative analysis shows that existing approaches can either fail dramatically, or else can impose impractical computational burdens due to numerical integration of reaction propensities, especially when cell ensembles are studied. Here we introduce the Extrande method which, given a simulated time course of dynamic network inputs, provides a conditionally exact and several orders-of-magnitude faster simulation solution. The new approach makes it feasible to demonstrate, using decision-making by a large population of quorum sensing bacteria, that robustness to fluctuations from upstream signaling places strong constraints on the design of networks determining cell fate. Our approach has the potential to significantly advance both understanding of molecular systems biology and design of synthetic circuits.

arXiv commentDownload source code for SI information

References in corpus (6)

• A modified Next Reaction Method for simulating chemical systems with time dependent propensities and delays
• Colored extrinsic fluctuations and stochastic gene expression
• Pooling single-cell recordings: Scalable inference through heterogeneous kinetics
• Mitochondrial Variability as a Source of Extrinsic Cellular Noise
• Uncoupled Analysis of Stochastic Reaction Networks in Fluctuating Environments
• Molecular finite-size effects in stochastic models of equilibrium chemical systems

Cited by in corpus (11)

• Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
• Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art
• Macroscopic phase resetting-curves determine oscillatory coherence and signal transfer in inter-coupled neural circuits
• Hybrid framework for the simulation of stochastic chemical kinetics
• Model reduction for stochastic chemical systems with abundant species
• Extending the linear-noise approximation to biochemical systems influenced by intrinsic noise and slow lognormally distributed extrinsic noise
• BioSimulator.jl: Stochastic simulation in Julia
• The exit time finite state projection scheme: bounding exit distributions and occupation measures of continuous-time Markov chains
• Uniformisation techniques for stochastic simulation of chemical reaction networks
• Low variance couplings for stochastic models of intracellular processes with time-dependent rate functions
• Simulation-based inference for stochastic nonlinear mixed-effects models with applications in systems biology