A data-integrated method for analyzing stochastic biochemical networks

Abstract

Variability and fluctuations among genetically identical cells under uniform experimental conditions stem from the stochastic nature of biochemical reactions. Understanding network function for endogenous biological systems or designing robust synthetic genetic circuits requires accounting for and analyzing this variability. Stochasticity in biological networks is usually represented using a continuous-time discrete-state Markov formalism, where the chemical master equation (CME) and its kinetic Monte Carlo equivalent, the stochastic simulation algorithm (SSA), are used. These two representations are computationally intractable for many realistic biological problems. Fitting parameters in the context of these stochastic models is particularly challenging and has not been accomplished for any but very simple systems. In this work, we propose that moment equations derived from the CME, when treated appropriately in terms of higher order moment contributions, represent a computationally efficient framework for estimating the kinetic rate constants of stochastic network models and subsequent analysis of their dynamics. To do so, we present a practical data-derived moment closure method for these equations. In contrast to previous work, this method does not rely on any assumptions about the shape of the stochastic distributions or a functional relationship among their moments. We use this method to analyze a stochastic model of a biological oscillator and demonstrate its accuracy through excellent agreement with CME/SSA calculations. By coupling this moment-closure method with a parameter search procedure, we further demonstrate how a model's kinetic parameters can be iteratively determined in order to fit measured distribution data.

Extrinsic Noise Suppression in Micro RNA mediated Incoherent Feedforward Loops

“…It turns out that the nonlinear product term ms in (43) makes computations non-trivial due to unclosed moment dynamics -the time evolution of the lower-order moments always depends on high-order moments. Moving forward, one general approach is to exploit closure schemes that approximate higher-order moments as nonlinear functions of lower-order moments [21]- [35]. Here we use an alternative approach based on linearizing the nonlinearity…”

A general moment expansion method for stochastic kinetic models

“…Several other approaches have been previously evaluated. Chevalier and El-Samad, 38 for example, approximated higher order moments using experimental data. Azunre et al 39 showed that for very small molecule numbers, using only two moments can lead to unstable results.…”

Bootstrapping least-squares estimates in biochemical reaction networks

“…RT-PCR or flow cytometry assays [4, 21]. Although this particular problem has been discussed extensively in the computational biology and statistical literature ([9, 12] and references therein), its asymptotic analysis was often dismissed due to a small number of longitudinal data points available from typical laboratory experiments.…”