Model reduction for stochastic CaMKII reaction kinetics in synapses by graph-constrained correlation dynamics.

A stochastic reaction network model of Ca(2+) dynamics in synapses (Pepke et al PLoS Comput. Biol. 6 e1000675) is expressed and simulated using rule-based reaction modeling notation in dynamical grammars and in MCell. The model tracks the response of calmodulin and CaMKII to calcium influx in synapses. Data from numerically intensive simulations is used to train a reduced model that, out of sample, correctly predicts the evolution of interaction parameters characterizing the instantaneous probability distribution over molecular states in the much larger fine-scale models. The novel model reduction method, 'graph-constrained correlation dynamics', requires a graph of plausible state variables and interactions as input. It parametrically optimizes a set of constant coefficients appearing in differential equations governing the time-varying interaction parameters that determine all correlations between variables in the reduced model at any time slice.

[1]  C. Rao,et al.  Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm , 2003 .

[2]  S. Orszag,et al.  Moment expansions in spatial ecological models and moment closure through Gaussian approximation , 2000, Bulletin of mathematical biology.

[3]  Prodromos Daoutidis,et al.  Model Reduction of Multiscale Chemical Langevin Equations: A Numerical Case Study , 2009, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[4]  Stefan Mihalas,et al.  A Dynamic Model of Interactions of Ca2+, Calmodulin, and Catalytic Subunits of Ca2+/Calmodulin-Dependent Protein Kinase II , 2010, PLoS Comput. Biol..

[5]  Ioannis G Kevrekidis,et al.  Equation-free implementation of statistical moment closures. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  Trevor Hastie,et al.  Regularization Paths for Generalized Linear Models via Coordinate Descent. , 2010, Journal of statistical software.

[7]  Eric Mjolsness,et al.  Time-ordered product expansions for computational stochastic system biology , 2012, Physical biology.

[8]  H. Risken.  Fokker-Planck Equation , 1984 .

[9]  Arjan van der Schaft,et al.  A model reduction method for biochemical reaction networks , 2014, BMC Systems Biology.

[10]  Scott B. Baden,et al.  Fast Monte Carlo Simulation Methods for Biological Reaction-Diffusion Systems in Solution and on Surfaces , 2008, SIAM J. Sci. Comput..

[11]  Chang Hyeong Lee,et al.  A moment closure method for stochastic reaction networks. , 2009, The Journal of chemical physics.

[12]  Daniel T. Gillespie,et al.  The multivariate Langevin and Fokker–Planck equations , 1996 .

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

[14]  Padhraic Smyth,et al.  Belief networks, hidden Markov models, and Markov random fields: A unifying view , 1997, Pattern Recognit. Lett..

[15]  Geoffrey E. Hinton,et al.  A Learning Algorithm for Boltzmann Machines , 1985, Cogn. Sci..

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

[17]  A. Singer,et al.  Maximum entropy formulation of the Kirkwood superposition approximation. , 2004, The Journal of chemical physics.

[18]  Mary B. Kennedy,et al.  Integration of biochemical signalling in spines , 2005, Nature Reviews Neuroscience.

[19]  Ilya Nemenman,et al.  Adiabatic coarse-graining and simulations of stochastic biochemical networks , 2009, Proceedings of the National Academy of Sciences.

[20]  T. Kurtz,et al.  Separation of time-scales and model reduction for stochastic reaction networks. , 2010, 1011.1672.

[21]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[22]  J. Hespanha.  Moment closure for biochemical networks , 2008, 2008 3rd International Symposium on Communications, Control and Signal Processing.

[23]  D. Gillespie.  The chemical Langevin equation , 2000 .

[24]  Cosimo Laneve,et al.  Formal molecular biology , 2004, Theor. Comput. Sci..

[25]  Heinz Koeppl,et al.  Lumpability abstractions of rule-based systems , 2010, Theor. Comput. Sci..

[26]  Ulf Dieckmann,et al.  A multiscale maximum entropy moment closure for locally regulated space–time point process models of population dynamics , 2011, Journal of mathematical biology.

[27]  Eric Mjolsness,et al.  Stochastic process semantics for dynamical grammars , 2006, Annals of Mathematics and Artificial Intelligence.

[28]  M. Khammash,et al.  The finite state projection algorithm for the solution of the chemical master equation. , 2006, The Journal of chemical physics.

[29]  Gábor Szederkényi,et al.  Model reduction in bio-chemical reaction networks with Michaelis-Menten kinetics , 2013, 2013 European Control Conference (ECC).

[30]  Ruth E. Baker,et al.  Modelling collective cell behaviour , 2014 .

[31]  Terrence J. Sejnowski,et al.  Synthesis of models for excitable membranes, synaptic transmission and neuromodulation using a common kinetic formalism , 1994, Journal of Computational Neuroscience.

[32]  Guido Sanguinetti,et al.  Validity conditions for moment closure approximations in stochastic chemical kinetics. , 2014, The Journal of chemical physics.

[33]  Dirk Lebiedz,et al.  Automatic Complexity Analysis and Model Reduction of Nonlinear Biochemical Systems , 2008, CMSB.

[34]  Michael Hucka,et al.  A Correction to the Review Titled "Rules for Modeling Signal-Transduction Systems" by W. S. Hlavacek et al. , 2006, Science's STKE.

[35]  Irad Yavneh,et al.  Why Multigrid Methods Are So Efficient , 2006, Computing in Science & Engineering.

[36]  J.P. Hespanha,et al.  Lognormal Moment Closures for Biochemical Reactions , 2006, Proceedings of the 45th IEEE Conference on Decision and Control.

[37]  Eric Mjolsness,et al.  Dependency diagrams and graph-constrained correlation dynamics: new systems for probabilistic graphical modeling , 2012 .