Efficient Multiscale Simulations of Circadian Rhythms Using Automated Phase Macomodelling Techniques

Circadian rhythm mechanisms involve multi-scale interactions between endogenous DNA-transcription oscillators. We present the application of efficient, numerically wellconditioned algorithms for abstracting (potentially large) systems of differential equation models of circadian oscillators into compact, accurate phase-only macromodels. We apply and validate our auto-extracted phase macromodelling technique on mammalian and Drosophila circadian systems, obtaining speedups of 9 − 13× over conventional timecourse simulation, with insignificant loss of accuracy, for single oscillators being synchronized by day/night light variations. Further, we apply the macromodels to simulate a system of 400 coupled circadian oscillators, achieving speedups of 240× and accurately reproducing synchronization and locking phenomena amongst the oscillators. We also present the use of parameterized phase macromodels for these circadian systems, and elucidate insights into circadian timing effects directly provided by our auto-extracted macromodels.

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

[2]  Jaijeet S. Roychowdhury,et al.  An efficient and robust technique for tracking amplitude and frequency envelopes in oscillators , 2005, ICCAD-2005. IEEE/ACM International Conference on Computer-Aided Design, 2005..

[3]  Jaijeet S. Roychowdhury,et al.  A robust envelope following method applicable to both non-autonomous and oscillatory circuits , 2006, 2006 43rd ACM/IEEE Design Automation Conference.

[4]  A. Goldbeter,et al.  Robustness of circadian rhythms with respect to molecular noise , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[5]  S. Usui,et al.  Range of entrainment of rat circadian rhythms to sinusoidal light-intensity cycles. , 2000, American journal of physiology. Regulatory, integrative and comparative physiology.

[6]  A. Demir,et al.  Phase noise in oscillators: a unifying theory and numerical methods for characterization , 2000 .

[7]  Jaijeet S. Roychowdhury,et al.  Fast, accurate prediction of PLL jitter induced by power grid noise , 2004, Proceedings of the IEEE 2004 Custom Integrated Circuits Conference (IEEE Cat. No.04CH37571).

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

[9]  Jaijeet S. Roychowdhury,et al.  Fast simulation of large networks of nanotechnological and biochemical oscillators for investigating self-organization phenomena , 2006, Asia and South Pacific Conference on Design Automation, 2006..

[10]  Jaijeet S. Roychowdhury,et al.  Macromodelling oscillators using Krylov-subspace methods , 2006, Asia and South Pacific Conference on Design Automation, 2006..

[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]  Alper Demir,et al.  A reliable and efficient procedure for oscillator PPV computation, with phase noise macromodeling applications , 2003, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[13]  A. Winfree.  Biological rhythms and the behavior of populations of coupled oscillators. , 1967, Journal of theoretical biology.

[14]  J. Roychowdhury,et al.  Capturing oscillator injection locking via nonlinear phase-domain macromodels , 2004, IEEE Transactions on Microwave Theory and Techniques.

[15]  A. Goldbeter,et al.  Modeling the mammalian circadian clock: sensitivity analysis and multiplicity of oscillatory mechanisms. , 2004, Journal of Theoretical Biology.

[16]  Zhichun Wang,et al.  PV-PPV: Parameter Variability Aware, Automatically Extracted, Nonlinear Time-Shifted Oscillator Macromodels , 2007, 2007 44th ACM/IEEE Design Automation Conference.