Development of a two-dimension manifold to represent high dimension mathematical models of the intracellular Mammalian circadian clock

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 approximately 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.