Stochastic parameter search for events

BackgroundWith recent increase in affordability and accessibility of high-performance computing (HPC), the use of large stochastic models has become increasingly popular for its ability to accurately mimic the behavior of the represented biochemical system. One important application of such models is to predict parameter configurations that yield an event of scientific significance. Due to the high computational requirements of Monte Carlo simulations and dimensionality of parameter space, brute force search is computationally infeasible for most large models.ResultsWe have developed a novel parameter estimation algorithm—Stochastic Parameter Search for Events (SParSE)—that automatically computes parameter configurations for propagating the system to produce an event of interest at a user-specified success rate and error tolerance. Our method is highly automated and parallelizable. In addition, computational complexity does not scale linearly with the number of unknown parameters; all reaction rate parameters are updated concurrently at the end of each iteration in SParSE. We apply SParSE to three systems of increasing complexity: birth-death, reversible isomerization, and Susceptible-Infectious-Recovered-Susceptible (SIRS) disease transmission. Our results demonstrate that SParSE substantially accelerates computation of the parametric solution hyperplane compared to uniform random search. We also show that the novel heuristic for handling over-perturbing parameter sets enables SParSE to compute biasing parameters for a class of rare events that is not amenable to current algorithms that are based on importance sampling.ConclusionsSParSE provides a novel, efficient, event-oriented parameter estimation method for computing parametric configurations that can be readily applied to any stochastic systems obeying chemical master equation (CME). Its usability and utility do not diminish with large systems as the algorithmic complexity for a given system is independent of the number of unknown reaction rate parameters.

[1]  J. Griffith.  Mathematics of cellular control processes. II. Positive feedback to one gene. , 1968, Journal of theoretical biology.

[2]  Michael Baym,et al.  The Separatrix Algorithm for Synthesis and Analysis of Stochastic Simulations with Applications in Disease Modeling , 2014, PloS one.

[3]  Concetto Spampinato,et al.  Parallel stochastic systems biology in the cloud , 2013, Briefings Bioinform..

[4]  Jose M. G. Vilar,et al.  Modeling network dynamics: the lac operon, a case study , 2004 .

[5]  Fabian J Theis,et al.  Method of conditional moments (MCM) for the Chemical Master Equation , 2013, Journal of Mathematical Biology.

[6]  Radek Erban,et al.  Fat versus Thin Threading Approach on GPUs: Application to Stochastic Simulation of Chemical Reactions , 2012, IEEE Transactions on Parallel and Distributed Systems.

[7]  Lorenzo Dematté,et al.  GPU computing for systems biology , 2010, Briefings Bioinform..

[8]  Linda R Petzold,et al.  Refining the weighted stochastic simulation algorithm. , 2009, The Journal of chemical physics.

[9]  Bernie J Daigle,et al.  Automated estimation of rare event probabilities in biochemical systems. , 2011, The Journal of chemical physics.

[10]  Linda R. Petzold,et al.  Accelerated maximum likelihood parameter estimation for stochastic biochemical systems , 2012, BMC Bioinformatics.

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

[12]  M. Mackey,et al.  Feedback regulation in the lactose operon: a mathematical modeling study and comparison with experimental data. , 2003, Biophysical journal.

[13]  D. Sherrington.  Stochastic Processes in Physics and Chemistry , 1983 .

[14]  Roshan M. D'Souza,et al.  Accelerating the Gillespie Exact Stochastic Simulation Algorithm Using Hybrid Parallel Execution on Graphics Processing Units , 2012, PloS one.

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

[16]  W. Ebeling.  Stochastic Processes in Physics and Chemistry , 1995 .

[17]  N. G. Van Kampen,et al.  Chapter III – STOCHASTIC PROCESSES , 2007 .

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

[19]  Chandra Krintz,et al.  Neptune: a domain specific language for deploying hpc software on cloud platforms , 2011, ScienceCloud '11.

[20]  Reuven Y. Rubinstein,et al.  Optimization of computer simulation models with rare events , 1997 .

[21]  Hong Li,et al.  Efficient Parallelization of the Stochastic Simulation Algorithm for Chemically Reacting Systems On the Graphics Processing Unit , 2010, Int. J. High Perform. Comput. Appl..

[22]  András Horváth,et al.  Parameter Estimation of Kinetic Rates in Stochastic Reaction Networks by the EM Method , 2008, 2008 International Conference on BioMedical Engineering and Informatics.