On Differences between Deterministic and Stochastic Models of Chemical Reactions: Schlögl Solved with ZI-Closure

Deterministic and stochastic models of chemical reaction kinetics can give starkly different results when the deterministic model exhibits more than one stable solution. For example, in the stochastic Schlögl model, the bimodal stationary probability distribution collapses to a unimodal distribution when the system size increases, even for kinetic constant values that result in two distinct stable solutions in the deterministic Schlögl model. Using zero-information (ZI) closure scheme, an algorithm for solving chemical master equations, we compute stationary probability distributions for varying system sizes of the Schlögl model. With ZI-closure, system sizes can be studied that have been previously unattainable by stochastic simulation algorithms. We observe and quantify paradoxical discrepancies between stochastic and deterministic models and explain this behavior by postulating that the entropy of non-equilibrium steady states (NESS) is maximum.

[1]  C. Gardiner,et al.  Stochastic models of firstorder nonequilibrium phase transitions in chemical reactions , 1975 .

[2]  Philipp Thomas,et al.  How accurate are the nonlinear chemical Fokker-Planck and chemical Langevin equations? , 2011, The Journal of chemical physics.

[3]  Elijah Roberts,et al.  Approximation and inference methods for stochastic biochemical kinetics—a tutorial review , 2017 .

[4]  Hong Qian,et al.  Stochastic dynamics and non-equilibrium thermodynamics of a bistable chemical system: the Schlögl model revisited , 2009, Journal of The Royal Society Interface.

[5]  Ramon Grima,et al.  A study of the accuracy of moment-closure approximations for stochastic chemical kinetics. , 2012, The Journal of chemical physics.

[6]  E. Jaynes.  The Minimum Entropy Production Principle , 1980 .

[7]  Yiannis N. Kaznessis,et al.  A closure scheme for chemical master equations , 2013, Proceedings of the National Academy of Sciences.

[8]  Hong Qian,et al.  Non-equilibrium phase transition in mesoscopic biochemical systems: from stochastic to nonlinear dynamics and beyond , 2011, Journal of The Royal Society Interface.

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

[10]  Eberhard O. Voit,et al.  150 Years of the Mass Action Law , 2015, PLoS Comput. Biol..

[11]  Hong Qian,et al.  Physical origins of entropy production, free energy dissipation, and their mathematical representations. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  H. Weiss,et al.  CO Oxidation on a CeOx/Pt(111) Inverse Model Catalyst Surface: Catalytic Promotion and Tuning of Kinetic Phase Diagrams , 2008 .

[13]  Jayajit Das,et al.  Purely stochastic binary decisions in cell signaling models without underlying deterministic bistabilities , 2007, Proceedings of the National Academy of Sciences.

[14]  Eduardo Sontag,et al.  Untangling the wires: A strategy to trace functional interactions in signaling and gene networks , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Y. Kaznessis,et al.  Efficient Moment Matrix Generation for Arbitrary Chemical Networks. , 2012, Chemical engineering science.

[16]  Metastable states in a nonlinear stochastic model , 1984 .

[17]  T. Kurtz.  Limit theorems and diffusion approximations for density dependent Markov chains , 1976 .

[18]  H. Sakaguchi.  Creation and Reproduction of Model Cells with Semipermeable Membrane(Cross-disciplinary physics and related areas of science and technology) , 2008, 0810.4745.

[19]  I Prigogine,et al.  The thermodynamic stability theory of non-equilibrium States. , 1974, Proceedings of the National Academy of Sciences of the United States of America.

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

[21]  I W Richardson.  On the principle of minimum entropy production. , 1969, Biophysical journal.

[22]  P. H. Constantino,et al.  Modeling stochasticity in biochemical reaction networks , 2016 .

[23]  W. Ebeling,et al.  Stochastic dynamics of a bistable reaction system , 1979 .

[24]  P. H. Constantino,et al.  Maximum Entropy Prediction of Non-Equilibrium Stationary Distributions for Stochastic Reaction Networks with Oscillatory Dynamics. , 2017, Chemical engineering science.

[25]  Vassilios Sotiropoulos,et al.  Analytical Derivation of Moment Equations in Stochastic Chemical Kinetics. , 2011, Chemical engineering science.

[26]  F. Schlögl,et al.  On thermodynamics near a steady state , 1971 .

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

[28]  P. Zhu.  Associated relaxation time and intensity correlation function of a bistable system driven by cross-correlation additive and multiplicative coloured noise sources , 2007 .

[29]  Hong Qian,et al.  Thermodynamic limit of a nonequilibrium steady state: Maxwell-type construction for a bistable biochemical system. , 2009, Physical review letters.

[30]  P. Gaspard.  Hamiltonian dynamics, nanosystems, and nonequilibrium statistical mechanics , 2006, cond-mat/0603382.

[31]  Farren J. Isaacs,et al.  Computational studies of gene regulatory networks: in numero molecular biology , 2001, Nature Reviews Genetics.

[32]  K. Christmann,et al.  Macroscopic and mesoscopic characterization of a bistable reaction system: CO oxidation on Pt(111) surface , 1999 .

[33]  T. Kurtz.  Strong approximation theorems for density dependent Markov chains , 1978 .

[34]  D. Gillespie.  Markov Processes: An Introduction for Physical Scientists , 1991 .

[35]  Yu-Chiun Wang,et al.  Spatial bistability of Dpp–receptor interactions during Drosophila dorsal–ventral patterning , 2005, Nature.

[36]  Grégoire Nicolis,et al.  EFFECT OF FLUCTUATIONS ON BIFURCATION PHENOMENA , 1979 .

[37]  A. Goldbeter,et al.  Bistability without Hysteresis in Chemical Reaction Systems: A Theoretical Analysis of Irreversible Transitions between Multiple Steady States , 1997 .

[38]  J. Collins,et al.  Construction of a genetic toggle switch in Escherichia coli , 2000, Nature.

[39]  H. Qian,et al.  A Quasistationary Analysis of a Stochastic Chemical Reaction: Keizer’s Paradox , 2007, Bulletin of mathematical biology.

[40]  Y. Kaznessis,et al.  On a theory of stability for nonlinear stochastic chemical reaction networks. , 2015, The Journal of chemical physics.

[41]  Ilya Prigogine,et al.  Modération et transformation irréversible des systèmes ouverts , 1945 .

[42]  A. Duncan,et al.  Noise-induced multistability in chemical systems: Discrete versus continuum modeling. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[43]  J. N. Kapur.  Maximum-entropy models in science and engineering , 1992 .

[44]  Thomas Wilhelm,et al.  The smallest chemical reaction system with bistability , 2009, BMC Systems Biology.

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

[46]  Mark Ptashne,et al.  A Genetic Switch, Phage Lambda Revisited , 2004 .

[47]  Grégoire Nicolis,et al.  Asymptotic properties of Markovian master equations , 1981 .

[48]  Philipp Thomas,et al.  Distribution Approximations for the Chemical Master Equation: Comparison of the Method of Moments and the System Size Expansion , 2015, 1509.09104.

[49]  John Ross,et al.  Thermodynamics and Fluctuations Far From Equilibrium , 2008, Entropy.

[50]  R. Fox,et al.  Qualms regarding the range of validity of the glansdorff-prigogine criterion for stability of non-equilibrium States. , 1974, Proceedings of the National Academy of Sciences of the United States of America.