Development of a Two-Dimension Manifold to Represent High Dimension Mathematical Models of the Intracellular Mammalian Circadian Clock

Abstract

A new focus for mathematical models of the circadian pacemaker involves the encapsulation within the models of detailed biological processes responsible for generating those circadian rhythms. Representing greater biological detail requires more mathematical equations, which pose a greater challenge for the analysis of such systems. Development of a method that retains the predominant dynamics while still providing biologically detailed information is advantageous. Two high-dimension mathematical models of intracellular mammalian circadian pacemakers, Leloup-Goldbeter and Forger-Peskin, with 19 and 73 differential equations, respectively, have been published. The authors projected each of these high-dimension models onto their respective manifold using proper orthogonal functions (POFs) obtained from the empirical decomposition of the model’s phase space to obtain a 2-dimension model. The resulting 2-dimension model, represented by 2 differential equations, predicts most of the salient characteristics of a biological clock including ~24-h oscillations, entrainment to an LD cycle, phase response curves, and the amplitude recovery dynamics that emerge following amplitude suppression. The manifold representation simplifies the mathematical analysis, since only 2 variables need to be observed and analyzed to understand the behavior of the biological clock. This reduced model derived from a model based on biological variables can be used for the development and analysis of mathematical models of the coupled mammalian oscillators to understand the dynamics of the integrated circadian pacemaker.

Keywords

circadian rhythm mammalian clock models model reduction proper orthogonal functions eigenvalue decomposition manifold models

References

Andronov AA , Vitt AA , and Khaikin SE (1966) Theory of Oscillators. New York: Dover.

Antle MC , Foley DK , Foley NC , and Silver R (2003) Gates and oscillators: A network model of the brain clock. J Biol Rhythms 18:339-350.

Balsers MA , Wriggers W , Oono Y , and Schulten K (1996) Principal component analysis and long time protein dynamics. J Phys Chem 100:2567-2572.

Broomhead DS and Kirby MA (2000) New approach to dimensionality reduction: Theory and algorithms. SIAM J Appl Math 60:2114-2142.

Chen J and Kang S-M (2000) Model reduction of nonlinear MEMS devices through arclength based Karhunen-Loe’ve decomposition. Proc IEEE Int Symposium on Circuits and Systems 2;457-460.

Forger DB , Dean DA , Gurdziel K , Leloup JC , Lee C , Von Gall C , Etchegaray JP , Kronauer RE , Goldbeter A , Peskin CS , et al. (2003) Development and validation of computational models for mammalian circadian oscillators. OMICS 7:387-400.

Forger DB and Kronauer RE (2002) Reconciling mathematical models of biological clocks by averaging on approximate manifolds. SIAM J Appl Math 62:1281-1296.

Forger DB and Peskin CS (2003) A detailed predictive model of the mammalian circadian clock. Proc Natl Acad Sci U S A 100:14806-14811.

Forger DB and Peskin CS (2004) Model based conjectures on mammalian clock controversies. J Theor Biol 230:533-539.

10.

Goldbeter A (1995) A model for circadian oscillations in the Drosophila period protein (PER). Proc Biol Sci 22:319-324.

11.

Goldbeter A (2002) Computational approaches to cellular rhythms. Nature 420:238-245.

12.

Guckenheimer J and Holmes P (1983) Nonlinear Oscillations, Dynamical Systems, and Bifurcation of Vector Fields. New York: Springer-Verlag.

13.

Jewett ME , Kronauer RE , and Czeisler CA (1991) Light-induced suppression of endogenous circadian amplitude in humans. Nature 350:59-62.

14.

Karhunen K (1946) Zur Spektraltheorie stochastischer prozesse. Annales Academiae Scientarum Fennicae 37.

15.

Kowalski ME and Jin J-M (2002) Karhunen-Loe’ve based model order reduction of nonlinear systems. Proc IEEE Antennas and Propagation Soc Int Symposium 1:552-555.

16.

Kronauer RE , Forger DB , and Jewett ME (1999) Quantifying human circadian pacemaker response to brief, extended, and repeated light stimuli over the photopic range. J Biol Rhythms 14:500-515.

17.

Kunz H and Achermann P (2003) Simulation of circadian rhythm generation in the suprachiasmatic nucleus with locally coupled self-sustained oscillators. J Theor Biol 224:63-78.

18.

Kwasniok F (1996) The reduction of complex dynamical systems using principal interaction patterns. Physica D 92:28-60.

19.

Lall S , Marsden JE , and Glavaski S (2002) A subspace approach to balanced truncation for model reduction of nonlinear control systems. Int J on Robust and Nonlinear Control 12;519-535.

20.

Leloup J-C and Goldbeter A (1998) A model for circadian rhythms in Drosophila incorporating the formation of a complex between the PER and TIM proteins. J Biol Rhythms 13:70-87.

21.

Leloup J-C and Goldbeter A (2003) Toward a detailed computational model for the mammalian circadian clock. Proc Natl Acad Sci U S A 100:7051-7056.

22.

Leloup J-C and Goldbeter A (2004) Modeling the mammalian circadian clock: sensitivity analysis and multiplicity of oscillatory mechanisms. J Theor Biol 230:541-562.

23.

Okamura H , Miyake S , Sumi Y , Yamaguchi S , Yasui A , Muijtjens M , Hoeijmakers JH , and van der Horst GT (1999) Photic induction of mPer1and mPer2 in cry-deficient mice lacking a biological clock. Science 286:2531-2534.

24.

Press WH , Teukolsky SA , Vetterling WT , and Flannery BP (2002) Numerical Recipes in C: The Art of Scientific Computing. Cambridge, UK: Cambridge University Press.

25.

Reppert SM and Weaver DR (2002) Coordination of circadian timing in mammals. Nature 418:935-941.

26.

Shearman LP , Sriram S , Weaver DR , Maywood ES , Chaves I , Zheng B , Kume K , Lee CC , van der Horst GT , Hastings MH , et al. (2000) Interacting molecular loops in the mammalian circadian clock. Science 288:1013-1019.

27.

Sirovich L and Kirby M (1989) An eigenfunction approach to large-scale transitional structures in jet flow. Phys Fluids A2:127-136.

28.

Smolen P , Baxter DA , and Byrne JH (2001) Modeling circadian oscillations with interlocking positive and negative feedback loops. J Neurosci 21:6644-6656.

29.

Strogatz SH (1994) Nonlinear Dynamics and Chaos. Reading, MA: Addison-Wesley.

30.

Tyson JJ , Hong CI , Thron CD , and Novak B (1999) A simple model of circadian rhythms based on dimerization and proteolysis of PER and TIM proteins. Biophys J 77:2411-2417.

31.

Ueda HR , Hagiwara M , and Kitano H (2001) Robust oscillations within the interlocked feedback model of Drosophila circadian rhythm. J Theor Biol 210:401-406.

32.

Winfree AT (1975) Resetting biological clocks. Phys Today 28:34-39.

33.

Winfree AT (2000) The Geometry of Biological Time, 2nd ed. New York: Springer-Verlag.

34.

Zhang Y , Henson MA , and Kevrekidis YG (2003) Nonlinear model reduction for dynamical analysis of cell population models. Chem Eng Sci 58:429-445.