Phase Reduction of Stochastic Biochemical Oscillators

A common method for analyzing the effects of molecular noise in chemical reaction networks is to approximate the underlying chemical master equation by a Fokker-Planck equation, and to study the statistics of the associated chemical Langevin equation. This so-called diffusion approximation involves performing a perturbation expansion with respect to a small dimensionless parameter $\epsilon =\Omega^{-1}$, where $\Omega$ characterizes the system size. For example, $\Omega$ could be the mean number of proteins produced by a gene regulatory network. In the deterministic limit $\Omega\rightarrow \infty$, the chemical reaction network evolves according to a system of ordinary differential equations based on classical mass action kinetics. In this paper we develop a phase reduction method for chemical reaction networks that support a stable limit cycle in the deterministic limit. We present a variational principle for the phase reduction, yielding an exact analytic expression for the resulting phase dynamics. We demonstrate that this decomposition is accurate over timescales that are exponential in the system size $\Omega$. This contrasts with the phase equation obtained under the diffusion approximation, which is only accurate up to times $O(\Omega)$. In particular, we show that for a constant $C$, the probability that the system leaves an $O(\zeta)$ neighborhood of the limit cycle before time $T$ scales as $T\exp\big(-C\Omega b\zeta^2 \big)$, where $b$ is the rate of attraction to the limit cycle. We illustrate our analysis using the example of a chemical Brusselator.

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

[2]  J. Elf,et al.  Fast evaluation of fluctuations in biochemical networks with the linear noise approximation. , 2003, Genome research.

[3]  Hiroya Nakao,et al.  Phase reduction approach to synchronisation of nonlinear oscillators , 2016, 1704.03293.

[4]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[5]  A. Winfree The geometry of biological time , 1991 .

[6]  Thomas G. Kurtz,et al.  Stochastic Analysis of Biochemical Systems , 2015 .

[7]  Stephen Coombes,et al.  Mathematical Frameworks for Oscillatory Network Dynamics in Neuroscience , 2015, The Journal of Mathematical Neuroscience.

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

[9]  Arkady Pikovsky,et al.  Synchronization and desynchronization of self-sustained oscillators by common noise. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Paul C. Bressloff,et al.  A Variational Method for Analyzing Stochastic Limit Cycle Oscillators , 2017, SIAM J. Appl. Dyn. Syst..

[11]  L. Tsimring Noise in biology , 2014, Reports on progress in physics. Physical Society.

[12]  Paul C. Bressloff,et al.  Stochastic switching in biology: from genotype to phenotype , 2017 .

[13]  Ruth E Baker,et al.  Turing's model for biological pattern formation and the robustness problem , 2012, Interface Focus.

[14]  Alper Demir,et al.  Phase computations and phase models for discrete molecular oscillators , 2012, EURASIP J. Bioinform. Syst. Biol..

[15]  T. Kurtz Limit theorems for sequences of jump Markov processes approximating ordinary differential processes , 1971, Journal of Applied Probability.

[16]  L. Popovic,et al.  Large deviations for multi-scale jump-diffusion processes , 2015, 1503.05990.

[17]  Johan Paulsson,et al.  Models of stochastic gene expression , 2005 .

[18]  David F. Anderson,et al.  Product-Form Stationary Distributions for Deficiency Zero Chemical Reaction Networks , 2008, Bulletin of mathematical biology.

[19]  Peter Hänggi,et al.  Bistable systems: Master equation versus Fokker-Planck modeling , 1984 .

[20]  T. Elston,et al.  Stochasticity in gene expression: from theories to phenotypes , 2005, Nature Reviews Genetics.

[21]  Yoshiki Kuramoto,et al.  Chemical Oscillations, Waves, and Turbulence , 1984, Springer Series in Synergetics.

[22]  Tobias Galla,et al.  Limit cycles, complex Floquet multipliers, and intrinsic noise. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[23]  A. Oudenaarden,et al.  Nature, Nurture, or Chance: Stochastic Gene Expression and Its Consequences , 2008, Cell.

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

[25]  C. Gardiner Handbook of Stochastic Methods , 1983 .

[26]  P. Protter Stochastic integration and differential equations , 1990 .

[27]  T. Kurtz Representations of Markov Processes as Multiparameter Time Changes , 1980 .

[28]  Michele Bonnin,et al.  Amplitude and phase dynamics of noisy oscillators , 2015, Int. J. Circuit Theory Appl..

[29]  Amir Dembo,et al.  Large deviations theory for Markov jump models of chemical reaction networks , 2017, The Annals of Applied Probability.

[30]  P. Gaspard,et al.  Biochemical Clocks and Molecular Noise: Theoretical Study of Robustness Factors , 2002 .

[31]  D. Gillespie Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .

[32]  H Koeppl,et al.  Deterministic characterization of phase noise in biomolecular oscillators , 2011, Physical biology.

[33]  Paul C. Bressloff,et al.  Stochastic Processes in Cell Biology , 2014, Interdisciplinary Applied Mathematics.

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

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

[36]  Andreas Hellander,et al.  Perspective: Stochastic algorithms for chemical kinetics. , 2013, The Journal of chemical physics.

[37]  P. Bressloff,et al.  A variational method for analyzing limit cycle oscillations in stochastic hybrid systems. , 2018, Chaos.

[38]  Richard P Boland,et al.  How limit cycles and quasi-cycles are related in systems with intrinsic noise , 2008, 0805.1607.

[39]  L. Reichl,et al.  An Introduction to Nonlinear Chemical Dynamics: Oscillations, Waves, Patterns, and Chaos By Irving R. Epstein (Brandeis University) and John A. Pojman (University of S. Mississippi). Oxford University Press: New York. 1998. 408 pp. $75.00. ISBN 0-19-509670-3. , 2000 .

[40]  David A. Rand,et al.  Long-time analytic approximation of large stochastic oscillators: Simulation, analysis and inference , 2016, bioRxiv.