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.

Mathematical Analysis of Some Reaction Networks Inducing Biological Growth/Decay Functions.

A Scalable Computational Framework for Establishing Long-Term Behavior of Stochastic Reaction Networks

A statistical interpretation of biologically inspired growth models

The two-step exponential decay reaction network: analysis of the solutions and relation to epidemiological SIR models with logistic and Gompertz type infection contact patterns

Markovian Dynamics on Complex Reaction Networks

Analytical Equilibrium Solutions of Biochemical Systems with Synthesis and Degradation

Product-Form Poisson-Like Distributions and Complex Balanced Reaction Systems

Multiscale Analysis of Reaction Networks

Biomass transfer on autocatalytic reaction network: a delay differential equation formulation

Metabolic Network Modelling: Including Stochastic Effects.

Mathematical Analysis of Chemical Reaction Systems

Reaction-diffusion kinetics in growing domains

Growth of microbial populations. Mathematical modeling, laboratory exercises, and model-based data analysis

Stochastic reaction networks with input processes: Analysis and applications to reporter gene systems

Qualitative analysis on a reaction-diffusion model arising from population dynamics

Global Dynamics of a Mathematical Model of Competition in the Chemostat: General Response Functions and Differential Death Rates

Long-term analysis of phenotypically structured models

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

An Analytical Method For Modeling First-Order Decay Networks

Characterizing the Conditions for Indefinite Growth in Open Chemical Reaction Networks

The Steady State Bifurcation for General Network-Organized Reaction-Diffusion Systems and Its Application in a Metapopulation Epidemic Model.