Accounting for extrinsic variability in the estimation of stochastic rate constants

Single-cell recordings of transcriptional and post-transcriptional processes reveal the inherent stochasticity of cellular events. However, to a large extent, the observed variability in isogenic cell populations is due to extrinsic factors, such as difference in expression capacity, cell volume and cell cycle stage—to name a few. Thus, such experimental data represents a convolution of effects from stochastic kinetics and extrinsic noise sources. Recent parameter inference schemes for single-cell data just account for variability because of molecular noise. Here, we present a Bayesian inference scheme that deconvolutes the two sources of variability and enables us to obtain optimal estimates of stochastic rate constants of low copy-number events and extract statistical information about cell-to-cell variability. In contrast to previous attempts, we model extrinsic noise by a variability in the abundance of mass-conserved species, rather than a variability in kinetic parameters. We apply the scheme to a simple model of the osmostress-induced transcriptional activation in budding yeast.

[1]  Darren J. Wilkinson Stochastic Modelling for Systems Biology , 2006 .

[2]  E. Klavins,et al.  Model reduction of stochastic processes using Wasserstein pseudometrics , 2008, 2008 American Control Conference.

[3]  S. Kay Fundamentals of statistical signal processing: estimation theory , 1993 .

[4]  J Timmer,et al.  Parameter estimation in stochastic biochemical reactions. , 2006, Systems biology.

[5]  Darren J. Wilkinson,et al.  Bayesian inference for nonlinear multivariate diffusion models observed with error , 2008, Comput. Stat. Data Anal..

[6]  David F. Anderson,et al.  Continuous Time Markov Chain Models for Chemical Reaction Networks , 2011 .

[7]  Toshio Yanagida,et al.  Single molecule dynamics in life science , 2008 .

[8]  Xiaohui Xie,et al.  Parameter inference for discretely observed stochastic kinetic models using stochastic gradient descent , 2010, BMC Systems Biology.

[9]  D. A. Mcquarrie Stochastic approach to chemical kinetics , 1967, Journal of Applied Probability.

[10]  C. Pesce,et al.  Regulated cell-to-cell variation in a cell-fate decision system , 2005, Nature.

[11]  Christopher M. Bishop,et al.  Pattern Recognition and Machine Learning (Information Science and Statistics) , 2006 .

[12]  P. Swain,et al.  Stochastic Gene Expression in a Single Cell , 2002, Science.

[13]  S. Chib,et al.  Understanding the Metropolis-Hastings Algorithm , 1995 .

[14]  Variable-at-a-time Implementations of Metropolis-Hastings , 2009 .

[15]  Matthew J. Beal Variational algorithms for approximate Bayesian inference , 2003 .

[16]  Stefan Hohmann,et al.  Yeast osmoregulation. , 2007, Methods in enzymology.

[17]  Steven Kay,et al.  Fundamentals Of Statistical Signal Processing , 2001 .

[18]  Heinz Koeppl,et al.  Probability metrics to calibrate stochastic chemical kinetics , 2010, Proceedings of 2010 IEEE International Symposium on Circuits and Systems.

[19]  Frank Allgöwer,et al.  A maximum likelihood estimator for parameter distributions in heterogeneous cell populations , 2010, ICCS.

[20]  C. Rao,et al.  Control, exploitation and tolerance of intracellular noise , 2002, Nature.

[21]  Lani F. Wu,et al.  Cellular Heterogeneity: Do Differences Make a Difference? , 2010, Cell.

[22]  Nasser M. Nasrabadi,et al.  Pattern Recognition and Machine Learning , 2006, Technometrics.

[23]  J. Raser,et al.  Control of Stochasticity in Eukaryotic Gene Expression , 2004, Science.

[24]  David Welch,et al.  Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems , 2009, Journal of The Royal Society Interface.

[25]  David F Anderson,et al.  A modified next reaction method for simulating chemical systems with time dependent propensities and delays. , 2007, The Journal of chemical physics.

[26]  Kunihiko Kaneko,et al.  Ubiquity of log-normal distributions in intra-cellular reaction dynamics , 2005, Biophysics.

[27]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[28]  N. Friedman,et al.  Structure and function of a transcriptional network activated by the MAPK Hog1 , 2008, Nature Genetics.

[29]  Manfred Opper,et al.  Approximate Inference for Stochastic Reaction processes , 2010, Learning and Inference in Computational Systems Biology.

[30]  Darren J. Wilkinson,et al.  Bayesian inference for a discretely observed stochastic kinetic model , 2008, Stat. Comput..

[31]  Mads Kærn,et al.  Noise in eukaryotic gene expression , 2003, Nature.

[32]  David F. Anderson,et al.  Error analysis of tau-leap simulation methods , 2009, 0909.4790.

[33]  Jim Freeman Probability Metrics and the Stability of Stochastic Models , 1991 .

[34]  Eulàlia de Nadal,et al.  Multilayered control of gene expression by stress‐activated protein kinases , 2010, The EMBO journal.

[35]  Edda Klipp,et al.  Biophysical properties of Saccharomyces cerevisiae and their relationship with HOG pathway activation , 2010, European Biophysics Journal.

[36]  Fabian Rudolf,et al.  Transient Activation of the HOG MAPK Pathway Regulates Bimodal Gene Expression , 2011, Science.

[37]  Pablo A. Iglesias,et al.  MAPK-mediated bimodal gene expression and adaptive gradient sensing in yeast , 2007, Nature.

[38]  D. Gillespie A rigorous derivation of the chemical master equation , 1992 .