Oscillator model reduction preserving the phase response: application to the circadian clock.

Mathematical model reduction is a long-standing technique used both to gain insight into model subprocesses and to reduce the computational costs of simulation and analysis. A reduced model must retain essential features of the full model, which, traditionally, have been the trajectories of certain state variables. For biological clocks, timing, or phase, characteristics must be preserved. A key performance criterion for a clock is the ability to adjust its phase correctly in response to external signals. We present a novel model reduction technique that removes components from a single-oscillator clock model and discover that four feedback loops are redundant with respect to its phase response behavior. Using a coupled multioscillator model of a circadian clock, we demonstrate that by preserving the phase response behavior of a single oscillator, we preserve timing behavior at the multioscillator level.

[1]  David E. Goldberg,et al.  Genetic and evolutionary algorithms come of age , 1994, CACM.

[2]  Eric Shea-Brown,et al.  On the Phase Reduction and Response Dynamics of Neural Oscillator Populations , 2004, Neural Computation.

[3]  Rudiyanto Gunawan,et al.  Isochron-based phase response analysis of circadian rhythms. , 2006, Biophysical journal.

[4]  Ulrich Maas,et al.  Simplifying chemical kinetics: Intrinsic low-dimensional manifolds in composition space , 1992 .

[5]  Premananda Indic,et al.  Development of a Two-Dimension Manifold to Represent High Dimension Mathematical Models of the Intracellular Mammalian Circadian Clock , 2006, Journal of biological rhythms.

[6]  Francis J Doyle,et al.  A novel computational model of the circadian clock in Arabidopsis that incorporates PRR7 and PRR9 , 2006, Molecular systems biology.

[7]  Erik D Herzog,et al.  Vasoactive intestinal polypeptide mediates circadian rhythmicity and synchrony in mammalian clock neurons , 2005, Nature Neuroscience.

[8]  M. Kramer,et al.  Sensitivity analysis of oscillatory systems , 1984 .

[9]  J. Tyson,et al.  A simple model of circadian rhythms based on dimerization and proteolysis of PER and TIM. , 1999, Biophysical journal.

[10]  B. Rusak,et al.  Neuropeptides phase shift the mammalian circadian pacemaker , 1995, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[11]  A. Goldbeter,et al.  Toward a detailed computational model for the mammalian circadian clock , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Daniel B Forger,et al.  Model based conjectures on mammalian clock controversies. , 2004, Journal of theoretical biology.

[13]  John Durkin,et al.  Tools and applications , 2002 .

[14]  Goldberg,et al.  Genetic algorithms , 1993, Robust Control Systems with Genetic Algorithms.

[15]  David Paydarfar,et al.  Starting, stopping, and resetting biological oscillators: in search of optimum perturbations. , 2004, Journal of theoretical biology.

[16]  M. Mavrovouniotis,et al.  Simplification of Mathematical Models of Chemical Reaction Systems. , 1998, Chemical reviews.

[17]  Daniel B. Forger,et al.  A detailed predictive model of the mammalian circadian clock , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[18]  T. Turányi Sensitivity analysis of complex kinetic systems. Tools and applications , 1990 .

[19]  Francis J. Doyle,et al.  ROBUST TIMEKEEPING IN CIRCADIAN NETWORKS: FROM GENES TO CELLS , 2007 .

[20]  Achim Kramer,et al.  Synchronization-Induced Rhythmicity of Circadian Oscillators in the Suprachiasmatic Nucleus , 2007, PLoS Comput. Biol..

[21]  Ioannis P. Androulakis,et al.  Kinetic mechanism reduction based on an integer programming approach , 2000 .

[22]  Erik D Herzog,et al.  GABA and Gi/o differentially control circadian rhythms and synchrony in clock neurons , 2006, Proceedings of the National Academy of Sciences.

[23]  Marianthi G. Ierapetritou,et al.  Design of flexible reduced kinetic mechanisms , 2001 .

[24]  Lawrence. Davis,et al.  Handbook Of Genetic Algorithms , 1990 .

[25]  Zbigniew Michalewicz,et al.  Genetic Algorithms + Data Structures = Evolution Programs , 1996, Springer Berlin Heidelberg.

[26]  Richard E. Kronauer,et al.  Reconciling Mathematical Models of Biological Clocks by Averaging on Approximate Manifolds , 2002, SIAM J. Appl. Math..

[27]  S. Lam,et al.  The CSP method for simplifying kinetics , 1994 .

[28]  S. Bernard,et al.  Spontaneous synchronization of coupled circadian oscillators. , 2005, Biophysical journal.

[29]  G. Briggs,et al.  A Note on the Kinetics of Enzyme Action. , 1925, The Biochemical journal.

[30]  Ignacio E. Grossmann,et al.  Part II. Future perspective on optimization , 2004, Comput. Chem. Eng..

[31]  Vasilios Manousiouthakis,et al.  KINETIC MODEL REDUCTION USING GENETIC ALGORITHMS , 1998 .

[32]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[33]  S. Strogatz From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators , 2000 .

[34]  Marc R. Roussel,et al.  Geometry of the steady-state approximation: Perturbation and accelerated convergence methods , 1990 .

[35]  Yoshiyuki Sakaki,et al.  Temporal Precision in the Mammalian Circadian System: A Reliable Clock from Less Reliable Neurons , 2004, Journal of biological rhythms.

[36]  Simon J. Fraser,et al.  The steady state and equilibrium approximations: A geometrical picture , 1988 .

[37]  Paul E. Brown,et al.  Extension of a genetic network model by iterative experimentation and mathematical analysis , 2005, Molecular systems biology.

[38]  Habib N. Najm,et al.  Model Reduction and Physical Understanding of Slowly Oscillating Processes: The Circadian Cycle , 2006, Multiscale Model. Simul..

[39]  L. Petzold,et al.  Model reduction for chemical kinetics: an optimization approach , 1999 .

[40]  Anthony Hall,et al.  Disruption of Hepatic Leptin Signaling Protects Mice From Age- and Diet-Related Glucose Intolerance , 2010, Diabetes.

[41]  Ignacio E. Grossmann,et al.  Retrospective on optimization , 2004, Comput. Chem. Eng..

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

[43]  E D Gilles,et al.  Reduction of mathematical models of signal transduction networks: simulation-based approach applied to EGF receptor signalling. , 2004, Systems biology.

[44]  M S Turner,et al.  Modelling genetic networks with noisy and varied experimental data: the circadian clock in Arabidopsis thaliana. , 2005, Journal of theoretical biology.

[45]  Eve Marder,et al.  Reduction of conductance-based neuron models , 1992, Biological Cybernetics.

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

[47]  M. A. Henson,et al.  A molecular model for intercellular synchronization in the mammalian circadian clock. , 2007, Biophysical journal.

[48]  Hans-Paul Schwefel,et al.  Evolution and Optimum Seeking: The Sixth Generation , 1993 .

[49]  Rudiyanto Gunawan,et al.  Sensitivity Measures for Oscillating Systems: Application to Mammalian Circadian Gene Network , 2008, IEEE Transactions on Automatic Control.