Bounds on stationary moments in stochastic chemical kinetics
نویسندگان
چکیده
In the stochastic formulation of chemical kinetics, the stationary moments of the population count of species can be described via a set of linear equations. However, except for some specific cases such as systems with linear reaction propensities, the moment equations are underdetermined as a lower order moment might depend upon a higher order moment. Here, we propose a method to find lower, and upper bounds on stationary moments of molecular counts in a chemical reaction system. The method exploits the fact that statistical moments of any positive-valued random variable must satisfy some constraints. Such constraints can be expressed as nonlinear inequalities on moments in terms of their lower order moments, and solving them in conjugation with the stationary moment equations results in bounds on the moments. Using two examples of biochemical systems, we illustrate that not only one obtains upper and lower bounds on a given stationary moment, but these bounds also improve as one uses more moment equations and utilizes the inequalities for the corresponding higher order moments. Our results provide avenues for development of moment approximations that provide explicit bounds on moment dynamics for systems whose dynamics are otherwise intractable.
منابع مشابه
Exact lower and upper bounds on stationary moments in stochastic biochemical systems.
In the stochastic description of biochemical reaction systems, the time evolution of statistical moments for species population counts is described by a linear dynamical system. However, except for some ideal cases (such as zero- and first-order reaction kinetics), the moment dynamics is underdetermined as lower-order moments depend upon higher-order moments. Here, we propose a novel method to ...
متن کاملDynamic Bounds on Stochastic Chemical Kinetic Systems Using Semidefinite Programming
The method of moments has been proposed as a potential means to reduce the dimensionality of the chemical master equation (CME) appearing in stochastic chemical kinetics. However, attempts to apply the method of moments to the CME usually result in the so-called closure problem. Several authors have proposed moment closure schemes, which allow them to obtain approximations of quantities of inte...
متن کاملRigorous bounds on the stationary distributions of the chemical master equation via mathematical programming
The stochastic dynamics of networks of biochemical reactions in living cells are typically modelled using chemical master equations (CMEs). The stationary distributions of CMEs are seldom solvable analytically, and few methods exist that yield numerical estimates with computable error bounds. Here, we present two such methods based on mathematical programming techniques. First, we use semidefin...
متن کاملMoment Growth Bounds on Continuous Time Markov Processes on Non-negative Integer Lattices
We consider time homogeneous Markov processes in continuous time with state space Z+ and provide two sufficient conditions and one necessary condition for the existence of moments E(‖X(t)‖) of all orders r ∈ N for all t ≥ 0. The sufficient conditions also guarantee an exponential in time growth bound for the moments. The class of processes studied have finitely many state independent jumpsize v...
متن کاملA general moment expansion method for stochastic kinetic models.
Moment approximation methods are gaining increasing attention for their use in the approximation of the stochastic kinetics of chemical reaction systems. In this paper we derive a general moment expansion method for any type of propensities and which allows expansion up to any number of moments. For some chemical reaction systems, more than two moments are necessary to describe the dynamic prop...
متن کامل