Comparison of Parameter Estimation Methods in Stochastic Chemical Kinetic Models: Examples in Systems Biology.

Stochastic chemical kinetics has become a staple for mechanistically modeling various phenomena in systems biology. These models, even more so than their deterministic counterparts, pose a challenging problem in the estimation of kinetic parameters from experimental data. As a result of the inherent randomness involved in stochastic chemical kinetic models, the estimation methods tend to be statistical in nature. Three classes of estimation methods are implemented and compared in this paper. The first is the exact method, which uses the continuous-time Markov chain representation of stochastic chemical kinetics and is tractable only for a very restricted class of problems. The next class of methods is based on Markov chain Monte Carlo (MCMC) techniques. The third method, termed conditional density importance sampling (CDIS), is a new method introduced in this paper. The use of these methods is demonstrated on two examples taken from systems biology, one of which is a new model of single-cell viral infection. The applicability, strengths and weaknesses of the three classes of estimation methods are discussed. Using simulated data for the two examples, some guidelines are provided on experimental design to obtain more information from a limited number of measurements.

[1]  T. Kalbfleisch,et al.  A stochastic model of gene transcription: an application to L1 retrotransposition events. , 2006, Journal of theoretical biology.

[2]  Suresh Kumar Poovathingal,et al.  Global parameter estimation methods for stochastic biochemical systems , 2010, BMC Bioinformatics.

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

[4]  Asger Hobolth,et al.  SIMULATION FROM ENDPOINT-CONDITIONED, CONTINUOUS-TIME MARKOV CHAINS ON A FINITE STATE SPACE, WITH APPLICATIONS TO MOLECULAR EVOLUTION. , 2009, The annals of applied statistics.

[5]  J. Rawlings,et al.  Approximate simulation of coupled fast and slow reactions for stochastic chemical kinetics , 2002 .

[6]  Eric L Haseltine,et al.  On the origins of approximations for stochastic chemical kinetics. , 2005, The Journal of chemical physics.

[7]  Darren J Wilkinson,et al.  Bayesian parameter inference for stochastic biochemical network models using particle Markov chain Monte Carlo , 2011, Interface Focus.

[8]  W. Huisinga,et al.  Solving the chemical master equation for monomolecular reaction systems analytically , 2006, Journal of mathematical biology.

[9]  Luonan Chen,et al.  Optimization meets systems biology , 2010, BMC Systems Biology.

[10]  Jared E. Toettcher,et al.  Stochastic Gene Expression in a Lentiviral Positive-Feedback Loop: HIV-1 Tat Fluctuations Drive Phenotypic Diversity , 2005, Cell.

[11]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[12]  Eric L Haseltine,et al.  Two classes of quasi-steady-state model reductions for stochastic kinetics. , 2007, The Journal of chemical physics.

[13]  Larry Wasserman,et al.  All of Statistics , 2004 .

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

[15]  Eric L Haseltine,et al.  The stochastic quasi-steady-state assumption: reducing the model but not the noise. , 2011, The Journal of chemical physics.

[16]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[17]  Darren J. Wilkinson,et al.  Bayesian Emulation and Calibration of a Stochastic Computer Model of Mitochondrial DNA Deletions in Substantia Nigra Neurons , 2009 .

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

[19]  Darren J. Wilkinson,et al.  Bayesian Sequential Inference for Stochastic Kinetic Biochemical Network Models , 2006, J. Comput. Biol..

[20]  John Yin,et al.  Kinetics of virus production from single cells. , 2012, Virology.

[21]  James B. Rawlings,et al.  Stochastic Kinetic Modeling of Vesicular Stomatitis Virus Intracellular Growth , 2009, Bulletin of mathematical biology.

[22]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[23]  Jack Y. Yang,et al.  BindN+ for accurate prediction of DNA and RNA-binding residues from protein sequence features , 2010, BMC Systems Biology.

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

[25]  Ankur Gupta.  Parameter estimation in deterministic and stochastic models of biological systems , 2013 .

[26]  A. Arkin,et al.  Stochastic mechanisms in gene expression. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[27]  E. Cox,et al.  Real-Time Kinetics of Gene Activity in Individual Bacteria , 2005, Cell.

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

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

[30]  A. Arkin,et al.  Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. , 1998, Genetics.

[31]  Qiang Zheng,et al.  Comparison of deterministic and stochastic kinetics for nonlinear systems , 1991 .

[32]  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.

[33]  L. T. Fan,et al.  Stochastic modeling of carbon oxidation , 1999 .

[34]  Linda R Petzold,et al.  The slow-scale stochastic simulation algorithm. , 2005, The Journal of chemical physics.

[35]  Junbin Gao,et al.  Simulated maximum likelihood method for estimating kinetic rates in gene expression , 2007, Bioinform..

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

[37]  Boseung Choi,et al.  Inference for discretely observed stochastic kinetic networks with applications to epidemic modeling. , 2012, Biostatistics.