An algebraic method to calculate parameter regions for constrained steady-state distribution in stochastic reaction networks

Abstract

Steady state is an essential concept in reaction networks. Its stability reflects fundamental characteristics of several biological phenomena such as cellular signal transduction and gene expression. Because biochemical reactions occur at the cellular level, they are affected by unavoidable fluctuations. Although several methods have been proposed to detect and analyze the stability of steady states for deterministic models, these methods cannot be applied to stochastic reaction networks. In this paper, we propose an algorithm based on algebraic computations to calculate parameter regions for constrained steady-state distribution of stochastic reaction networks, in which the means and variances satisfy some given inequality constraints. To evaluate our proposed method, we perform computer simulations for three typical chemical reactions and demonstrate that the results obtained with our method are consistent with the simulation results.

arXiv comment17 pages, 4 figures

References in corpus (9)

• Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
• A scalable computational framework for establishing long-term behavior of stochastic reaction networks
• Transitions Induced by the Discreteness of Molecules in a Small Autocatalytic System
• Comparison of different moment-closure approximations for stochastic chemical kinetics
• A finite state projection algorithm for the stationary solution of the chemical master equation
• A general moment expansion method for stochastic kinetic models
• The linear-noise approximation and the chemical master equation exactly agree up to second-order moments for a class of chemical systems
• Validity conditions for moment closure approximations in stochastic chemical kinetics
• Optimal implementations for reliable circadian clocks