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.

What problem does this paper attempt to address?

This paper attempts to solve the problem of stochastic simulation of biomolecular networks in dynamic environments. Specifically, it focuses on how to effectively simulate the behavior of these networks in the case of continuous changes in cells and their surrounding environments. Traditional stochastic simulation algorithms (such as SSA, that is, Stochastic Simulation Algorithm) have limitations when dealing with such dynamic inputs, either producing significant errors or bringing a huge computational burden. Therefore, this paper proposes a new method - Extrande (Extra Reaction Algorithm for Networks in Dynamic Environments), for condition - exact and computationally efficient stochastic simulation of biomolecular networks in dynamic environments. The core of the Extrande method lies in the introduction of an additional "virtual" reaction channel, whose occurrence does not change the number of molecules in the system. The propensity of this additional channel is designed to fluctuate over time, so that when it is added to the sum of the propensities of all other reaction channels, it becomes constant between events and is equal to the upper bound of the sum of the propensities in the original system. This method takes advantage of the exogenous characteristics of dynamic inputs, that is, assuming that the influence of the network on external inputs is negligible, so that future input trajectories can be used to find an effective upper bound. In this way, the Extrande method can achieve accurate stochastic simulation of downstream reaction networks without performing analytical or numerical integration. This not only improves the accuracy of the simulation but also greatly reduces the computational time, especially when simulating large cell populations. For example, the paper shows the distribution of 1,000 bacterial cells differentiating into two phenotypes under stress conditions simulated by the Extrande method. This process would require more than 14 years of computational time if using traditional integration methods, while it only takes 2.3 months using Extrande. In addition, the paper also explores the influence of upstream signal fluctuations on the design of cell - fate - determination networks. Through a comparative study of wild - type and synthetic decision - making networks (SynDM), it is found that the wild - type network exhibits higher robustness in the face of upstream signal fluctuations and downstream process heterogeneity. These findings not only enhance our understanding of molecular systems biology but also provide valuable insights for the design of synthetic circuits.

Stochastic Simulation of Biomolecular Networks in Dynamic Environments

Computationally Efficient Modelling of Stochastic Spatio-Temporal Dynamics in Biomolecular Networks

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

Stochastic Differential Equations for Practical Simulation of Gene Circuits

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

Quantifying the roles of space and stochasticity in computer simulations for cell biology and cellular biochemistry

The Roles of Space and Stochasticity in Computational Simulations of Cellular Biochemistry: Quantitative Analysis and Qualitative Insights

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

Stochastic single-particle based simulations of cellular signaling embedded into computational models of cellular morphology

Single molecule simulations in complex geometries with embedded dynamic one-dimensional structures

Modelling non-Markovian dynamics in biochemical reactions

Parallel Replica Dynamics Method for Bistable Stochastic Reaction Networks: Simulation and Sensitivity Analysis.

Computational Framework for In-Silico Study of Virtual Cell Biology Via Process Simulation and Multiscale Modeling

Global Optimization Approach for Parameter Estimation in Stochastic Dynamic Models of Biosystems

A Stochastic Framework for the Design of Transient and Steady State Behavior of Biochemical Reaction Networks

Identifying Latent Dynamic Components in Biological Systems

Coarse-graining stochastic biochemical networks: quasi-stationary approximation and fast simulations using a stochastic path integral technique

Direct ab initio MD simulation of silver ion diffusion in chalcogenide glasses

Exact lattice-based stochastic cell culture simulation algorithms incorporating spontaneous and contact-dependent reactions

Green's Function Reaction Dynamics: a new approach to simulate biochemical networks at the particle level and in time and space

Exact Parallelization of the Stochastic Simulation Algorithm for Scalable Simulation of Large Biochemical Networks