A scalable computational framework for establishing long-term behavior of stochastic reaction networks

Abstract

Reaction networks are systems in which the populations of a finite number of species evolve through predefined interactions. Such networks are found as modeling tools in many biological disciplines such as biochemistry, ecology, epidemiology, immunology, systems biology and synthetic biology. It is now well-established that, for small population sizes, stochastic models for biochemical reaction networks are necessary to capture randomness in the interactions. The tools for analyzing such models, however, still lag far behind their deterministic counterparts. In this paper, we bridge this gap by developing a constructive framework for examining the long-term behavior and stability properties of the reaction dynamics in a stochastic setting. In particular, we address the problems of determining ergodicity of the reaction dynamics, which is analogous to having a globally attracting fixed point for deterministic dynamics. We also examine when the statistical moments of the underlying process remain bounded with time and when they converge to their steady state values. The framework we develop relies on a blend of ideas from probability theory, linear algebra and optimization theory. We demonstrate that the stability properties of a wide class of biological networks can be assessed from our sufficient theoretical conditions that can be recast as efficient and scalable linear programs, well-known for their tractability. It is notably shown that the computational complexity is often linear in the number of species. We illustrate the validity, the efficiency and the wide applicability of our results on several reaction networks arising in biochemistry, systems biology, epidemiology and ecology. The biological implications of the results as well as an example of a non-ergodic biological network are also discussed.

31 pages, 9 figures

References in corpus (4)

• A general moment expansion method for stochastic kinetic models
• Exact protein distributions for stochastic models of gene expression using partitioning of Poisson processes
• Dynamical Properties of Discrete Reaction Networks
• Determining the long-term behavior of cell populations: A new procedure for detecting ergodicity in large stochastic reaction networks

Cited by in corpus (32)

• Antithetic Integral Feedback ensures robust perfect adaptation in noisy biomolecular networks
• A finite state projection algorithm for the stationary solution of the chemical master equation
• State and parameter estimation from exact partial state observation in stochastic reaction networks
• Stochastic filtering for multiscale stochastic reaction networks based on hybrid approximations
• Replica-Mean-Field Limits of Fragmentation-Interaction-Aggregation Processes
• Phenomenology and dynamics of competitive ecosystems beyond the niche-neutral regimes
• Multilevel Hybrid Split-Step Implicit Tau-Leap
• Bounding Transient Moments of Stochastic Chemical Reactions
• Strong convergence for split-step methods in stochastic jump kinetics
• A finite state projection method for steady-state sensitivity analysis of stochastic reaction networks
• Stationary distributions of continuous-time Markov chains: a review of theory and truncation-based approximations
• Robust and structural ergodicity analysis of stochastic biomolecular networks involving synthetic antithetic integral controllers
• Learning-Based Importance Sampling via Stochastic Optimal Control for Stochastic Reaction Networks
• Robust ergodicity and tracking in antithetic integral control of stochastic biochemical reaction networks
• Classification and threshold dynamics of stochastic reaction networks
• Determining the long-term behavior of cell populations: A new procedure for detecting ergodicity in large stochastic reaction networks
• The asymptotic tails of limit distributions of continuous time Markov chains
• Long-time asymptotics of stochastic reaction systems
• Moment growth bounds on continuous time Markov processes on non-negative integer lattices
• Stochastic filters based on hybrid approximations of multiscale stochastic reaction networks
• Importance sampling for a robust and efficient multilevel Monte Carlo estimator for stochastic reaction networks
• Robust control of biochemical reaction networks via stochastic morphing
• Mathematical Models of Gene Expression
• Some network conditions for positive recurrence of stochastically modeled reaction networks
• An algebraic method to calculate parameter regions for constrained steady-state distribution in stochastic reaction networks
• Stability and performance analysis of linear positive systems with delays using input-output methods
• Steady State Sensitivity Analysis of Continuous Time Markov Chains
• Variance of finite difference methods for reaction networks with non-Lipschitz rate functions
• Optimal and $H_{\infty }$ Control of Stochastic Reaction Networks
• Parameter Invariance Analysis of Moment Equations Using Dulmage-Mendelsohn Decomposition
• Hybrid $L_{\infty }\times {\ell }_{\infty }$-Performance Analysis and Control of Linear Time-Varying Impulsive and Switched Positive Systems
• Optimization-based Lyapunov function construction for continuous-time Markov chains with affine transition rates