Global parameter search reveals design principles of the mammalian circadian clock
© Locke et al; licensee BioMed Central Ltd. 2008
Received: 30 May 2007
Accepted: 29 February 2008
Published: 29 February 2008
Virtually all living organisms have evolved a circadian (~24 hour) clock that controls physiological and behavioural processes with exquisite precision throughout the day/night cycle. The suprachiasmatic nucleus (SCN), which generates these ~24 h rhythms in mammals, consists of several thousand neurons. Each neuron contains a gene-regulatory network generating molecular oscillations, and the individual neuron oscillations are synchronised by intercellular coupling, presumably via neurotransmitters. Although this basic mechanism is currently accepted and has been recapitulated in mathematical models, several fundamental questions about the design principles of the SCN remain little understood. For example, a remarkable property of the SCN is that the phase of the SCN rhythm resets rapidly after a 'jet lag' type experiment, i.e. when the light/dark (LD) cycle is abruptly advanced or delayed by several hours.
Here, we describe an extensive parameter optimization of a previously constructed simplified model of the SCN in order to further understand its design principles. By examining the top 50 solutions from the parameter optimization, we show that the neurotransmitters' role in generating the molecular circadian rhythms is extremely important. In addition, we show that when a neurotransmitter drives the rhythm of a system of coupled damped oscillators, it exhibits very robust synchronization and is much more easily entrained to light/dark cycles. We were also able to recreate in our simulations the fast rhythm resetting seen after a 'jet lag' type experiment.
Our work shows that a careful exploration of parameter space for even an extremely simplified model of the mammalian clock can reveal unexpected behaviours and non-trivial predictions. Our results suggest that the neurotransmitter feedback loop plays a crucial role in the robustness and phase resetting properties of the mammalian clock, even at the single neuron level.
Almost all organisms have evolved to co-ordinate their activities with the many changes in the environment caused by the earth's rotation . A circadian clock generates biological rhythms with a period of approximately 24 hours (h), and is found in most eukaryotes and photosynthetic bacteria. In mammals, a central pacemaker exists in the suprachiasmatic nucleus (SCN) of the hypothalamus, which generates and communicates a circadian rhythm to other parts of the brain, and to peripheral tissues. The clock is entrained to the day/night cycle, with the SCN receiving light information from the retina through the retinohypothalamic tract, but is also self-sustained under constant conditions .
The SCN consists of two nuclei located above the optic chiasm, each containing around 10000 neurons. Dissociated neurons are able to produce circadian oscillations under constant conditions, although less robustly, with a period of between 20 and 28 h [3, 4]. The oscillation in individual neurons is generated by a network of interlocking feedback loops involving clock genes and proteins, similar to that seen in other organisms. Key components in this machinery are the clock proteins PERIOD1 (PER1), PER2 and CRYPTOCHROME1 (CRY1), CRY2, and the transcription factors CLOCK and BMAL1. These feedback loops have been extensively studied  and modelled [5–7]. A coupling mechanism must exist between the individual neurons in order to generate a robust 24 h collective self-sustained rhythm under constant conditions of the SCN tissue. Several candidate coupling factors have been discussed . Recently, vasoactive intestinal polypeptide (VIP) has been implicated as one of the possible coupling factors, since mutations in genes coding for VIP as well as for its receptor VPAC2 disrupt normal synchronisation of the clock neurons within the SCN, and maintenance of molecular timekeeping within individual SCN neurons . It has also recently been shown that SCN intercellular coupling can make the system more robust against genetic perturbations .
Populations of individual oscillators with variable periods have been shown to be synchronised by a coupling factor in a variety of contexts , including pacemaker cells in the heart, synchronization of firefly flashes, and circadian rhythms in Drosophila. The properties of coupled oscillators have been researched in great detail [12, 13], and have long been studied in relation to circadian rhythms . Several theoretical studies have also been made specifically examining the coupling mechanism of the SCN [15–19]. Gonze et al. previously studied theoretically the synchronising properties of the SCN, and showed that a global coupling mechanism relying on a mean field of neurotransmitter is sufficient to synchronize a population of 10,000 cells . These results have been generalized to more complex single cell models and a variety of coupling schemes . A recent theoretical study has also examined in detail the role of VIP in synchronising the SCN .
A fundamental property of the SCN is that the phase of the SCN rhythm resets rapidly after a 'jet lag' type experiment, where the light/dark (LD) cycle is abruptly advanced or delayed by several hours. It has been shown experimentally that the levels of the clock protein PER1 in the SCN are almost completely re-entrained within one day of a 6 h delay or advance in LD cycle, being completely reset in around 6 days [20–22]. It remains an open question how ~10000 neurons are able to respond so quickly to a shift in the entraining signal, even though phase response curve experiments have shown that the SCN rhythm appears to be phase shifted only by a maximum of a few hours by light and neurotransmitter pulses. A Phase Response Curve (PRC) is a plot of phase-shifts as a function of circadian phase of a stimulus such as light pulses, temperature pulses, or pulses of drugs or chemicals. The SCN shows a weak (type 1) PRC to pulses of light and a variety of neurotransmitters including VIP and gastrin-releasing peptide (GRP) [23–27].
The phenomenon of fast phase resetting, and in fact, the balance between the two primary synchronisers of the SCN, light, and intrinsic neurotransmitters, is not well understood. In the Gonze et al. model, entrainment to a 12:12 h LD cycle resulted in quasiperiodic behaviour, with small variations in the amplitude from cycle to cycle. A shift in LD cycles caused very long transients, taking around twenty days for the phase to reset to a 6 h shift in LD cycle, rather than the 1–6 days seen in experiment (see Additional file 1: Supplementary Figure S1, Additional file 2: Supplementary figure legends). However, it remains unclear whether this behaviour is a general feature of the network equations (described in Additional file 3, which contains the supplementary information), or if this is due to the parameters chosen.
In this work, we use a global parameter search scheme [28, 29] to investigate the resetting properties of a simplified model of the SCN . We investigate the top 50 parameter sets found using our parameter optimization scheme, and show that the neurotransmitter feedback loop can be crucial for oscillations, even at the single cell level. We also show that if the neurotransmitter drives a system of damped, rather than self-sustained oscillators, it is possible to achieve fast phase resetting in simulated 'jet lag' type experiments even with this simple model of the SCN. Finally, we go on to show that fast resetting properties are still possible even in the more experimentally realistic case of only a sub fraction of the oscillators being light-activated.
Results and Discussion
Global parameter screen reveals that the neurotransmitter feedback loop can be essential for single neuron oscillations in the SCN
In order to estimate the 16 unknown parameters in this model, we defined a cost function (CF) that gave a low score for solutions that had a limit cycle oscillation of the mean field in constant conditions with a period of approximately 24 h, and were well synchronised (see Additional file 3, and  for details of the optimisation scheme). We calculated the CF for a million candidate parameter sets chosen quasi-randomly, and then we picked the parameter set giving the lowest CF score as our optimal solution for further detailed analysis. We also used the top 50 parameter sets to further explore the generalised properties of synchronisation.
Simulations using our optimal parameter set initially appear similar to simulations using handpicked parameters . However, in contrast to that work, a reduction of neurotransmitter strength leads to a very rapid reduction in the amplitude of the oscillation, even for a single neuron simulation. A reduction in neurotransmitter strength could be caused experimentally by several factors, such as pharmacological or genetic perturbations of the VIP/VPAC2 signalling pathway. If the coupling constant is reduced to 80% of the original value a single neuron generates a self-sustaining oscillation, but with a much smaller amplitude (Figure 2E), and complete removal of the neurotransmitter loop results in arrhythmia (Figure 2F). This means the neurotransmitter positive feedback loop substantially contributes to even single cell oscillations. Isolated neurons should oscillate, albeit out of phase with each other, unless the autocrine feedback loop is disrupted. Other modelling work has shown similar results concerning the possibility of the neurotransmitter helping to drive neuronal oscillations [17–19], and recent experimental data also offer support for this finding: in mice with genetic lesions of the VPAC2 receptor oscillations of individual neurons are compromised . A subset of neurons is unable to oscillate, whilst other neurons remain rhythmic, although quite desynchronized from each other. In future experimental work it would be interesting to see whether knocking out other candidate coupling mechanisms, such as Prokineticin 2 , in tandem could increase the number of neurons not capable of self-sustained oscillations.
We investigated further the properties of the top fifty parameter sets found through our optimization. First we examined the oscillatory state of the solutions in order to see whether the oscillations were driven by the neurotransmitter, as in the optimal parameter set. (This was carried out by calculating the eigenvalues of the Jacobian, evaluated at the unique fixed point given by the single-cell core oscillator, see Additional file 3 for details). When the neurotransmitter loop was included (it is conceivable that there is an autocrine feedback of the neurotransmitter released by a given cell), all parameter sets exhibited an unstable fixpoint and sustained oscillations. However, without the neurotransmitter loop but with a fixed neurotransmitter mean field level, corresponding to the average level for the four-dimensional single-cell core oscillator, 44 parameter sets showed a stable fixed point with imaginary eigenvalues, corresponding to damped (not self-sustained) oscillations (See Additional file 3: Supplementary Table S1). It thus appears that damped single-cell oscillators could be advantageous for efficient synchronization. It would be interesting to conceive experimental strategies to test this prediction. Further, the present study suggests quite general qualitative properties of coupled damped oscillators, and sets the stage for a more rigorous analytical approach to our numerical findings.
We then went on to examine the distributions of the optimized parameters and their correlations, in order to gain some insight into the design of the 50 networks. The parameters of these are quite variable (Additional file 1: Supplementary Figure S2). Of particular note is the fact that synchronised oscillations were possible for Hill coefficients varying from 0.7 – 9.8. Previous models have required high Hill coefficients, representing highly cooperative binding of proteins, in order to generate oscillations, although there is no experimental evidence for this. In fact, for a single three-variable genetic circuit of Goodwin-type with linear decay terms , the Hill coefficient must be greater than 8 in order to generate self-sustained oscillations .
We found no significant correlations between a single variable and the CF, which is not unexpected due to the nonlinearity and complexity of the system. However, two weak correlations deserve mentioning: First, a low value of the Michaelis constant for the mRNA decay seems to be slightly advantageous for effective synchronization, (K2 and the CF had a correlation coefficient of 0.16). Second, a low Vmax for Y protein degradation also seems to be slightly favourable for the synchronization (v4) and the CF had a correlation coefficient of 0.19). Correlations were stronger between the parameters themselves (Additional file 3: Supplementary Table S2), most notably between the Hill coefficient and v6, the Vmax value for Z degradation (correlation coefficient -0.52). In other words, a capability for a high degradation of the transcriptional inhibitor relaxes the demand for a switch-like transcriptional inhibition of X. Furthermore, although it seems that it is an advantage for effective synchronization if the single cell core oscillator exhibits damped oscillations (with a constant mean-field level), it seems not to matter how damped it is. More technically, this is because the real part of the largest eigenvalue of the system, using the fixed average mean field, was uncorrelated to the cost function (correlation coefficient 0.05), which suggests that the absolute value of the real part of this eigenvalue is not important for a robust multicellular oscillator, as long as it is negative. Last, we calculated the half-lives for the clock mRNA X, clock protein Y, and transcriptional inhibitor Z, respectively (see Additional file 3: Supplementary Table S3). We found a significant positive correlation between the mRNA half-life and the value of the cost function (see Additional file 1: Supplementary Figure S3). Thus, a prediction from our model is that a not too long half-life for a clock mRNA molecule, perhaps less than an hour, makes efficient synchronization easier.
Fast re-entrainment is possible with optimal parameter set
Light activated pacemaker cells can synchronize the SCN
We have used a global parameter scan of a previously developed simplified model of the mammalian circadian clock (Figure 1)  to reveal several important design principles. First, we have shown that synchronisation is possible and perhaps more efficient if the individual oscillators are damped (i.e. are not self-sustained) without an oscillatory neuropeptide concentration (Figure 2) – 88% of the best parameter combinations corresponded to damped individual oscillators. This is supported by recent experimental results , which showed that without efficient VIP/VPAC2 signalling, a subset of cells are unable to oscillate in a self-sustained manner. We argue that this result, showing that damped core oscillators are preferred for efficient synchronization, has particular weight, since it in no way is due to model assumptions, but rather emerges from the evaluation of a million randomly selected parameter combinations. Second, we have observed that the experimentally seen extremely fast phase resetting of the SCN after a shift in LD cycle can occur in our model when the neuropeptide positive feedback loop drives the oscillation (Figure 3). Future work should test whether fast resetting of the SCN really requires the neuropeptide feedback loop to drive the oscillation. Finally, we see that fast resetting is still possible under the more physiological assumption that only a subset of the simulated neurons are light activated (Figure 4).
Several mathematical models have been constructed that have revealed important aspects of the function of the feedback loops of the mammalian circadian clock [5–7, 39]. Recently modelling has also been used to directly understand as well as propose experiments [40, 41]. Our work shows that unexpected behaviours can be revealed through a thorough scan of parameter space, even with an extremely simplified model of the mammalian clock. The fact that fast phase resetting of the mammalian clock appears possible when the individual oscillators are not self-sustained has profound implications for our understanding of the experimentally observed result that sub-populations of 'uncoupled' neurons have been seen to be arrhythmic . Simulations of the coupled system of damped oscillators showed robust oscillations even with large variations in the periods of the individual oscillators (Figure 2D), which fits with recent experimental evidence that coupling of the SCN neurons can make the system more robust to genetic perturbations . It also appears that having weakly rhythmic individual oscillators that are strongly coupled can make the system more responsive to external cues such as light (Figure 3).
The SCN structure is extremely complex . There is a great deal of spatial heterogeneity, and individual neurons have been shown to differ in their neuropeptide expression, light responsiveness, phase , and free running period. This complexity cannot be reproduced with the extremely simplified model of the SCN used in our study. In future it would be interesting to study the heterogeneity of the SCN, and to examine a more realistic model of its function. However, as we have seen, simplified models are extremely useful for uncovering general mechanisms of genetic networks. More detailed models of the SCN require an increase in the number of unknown parameters [18, 19], and this makes our global parameter search method less comprehensive.
We found the optimized parameters for the core negative feedback loop to be extremely varied (Additional file 1: Supplementary Figure S2). A single-cell model with a Hill coefficient of repression as low as 0.7 was still able to function as the building block of a robustly oscillating cellular network when driven by the neuropeptide mean-field. This is very different behaviour from previous circadian clock models, where a high value of the Hill coefficient was crucial to obtain limit cycle oscillations . Our work shows that future experiments and mathematical models must take into account the neuropeptide feedback loop in order to understand the function of the mammalian clock, even at the single neuron level.
We would like to thank S. Bernard and S. Legewie for useful discussions. This work was supported by the European Commission (BioSim Network, contract No. LSHB-CT-2004.005137 and the Deutsche Forschungsgemeinschaft (SFB 618). Research in AKs lab is also supported by the 6th EU framework program EUCLOCK. Computer facilities were provided by the Centre for Scientific Computing at the University of Warwick.
- Dunlap JC: Molecular bases for circadian clocks. Cell. 1999, 96 (2): 271-290. 10.1016/S0092-8674(00)80566-8View ArticlePubMedGoogle Scholar
- Hastings MH, Herzog ED: Clock genes, oscillators, and cellular networks in the suprachiasmatic nuclei. J Biol Rhythms. 2004, 19 (5): 400-413. 10.1177/0748730404268786View ArticlePubMedGoogle Scholar
- Honma S, Nakamura W, Shirakawa T, Honma K: Diversity in the circadian periods of single neurons of the rat suprachiasmatic nucleus depends on nuclear structure and intrinsic period. Neurosci Lett. 2004, 358 (3): 173-176. 10.1016/j.neulet.2004.01.022View ArticlePubMedGoogle Scholar
- Welsh DK, Logothetis DE, Meister M, Reppert SM: Individual neurons dissociated from rat suprachiasmatic nucleus express independently phased circadian firing rhythms. Neuron. 1995, 14 (4): 697-706. 10.1016/0896-6273(95)90214-7View ArticlePubMedGoogle Scholar
- Forger DB, Peskin CS: A detailed predictive model of the mammalian circadian clock. Proc Natl Acad Sci U S A. 2003, 100 (25): 14806-14811. 10.1073/pnas.2036281100PubMed CentralView ArticlePubMedGoogle Scholar
- Becker-Weimann S, Wolf J, Herzel H, Kramer A: Modeling feedback loops of the Mammalian circadian oscillator. Biophysical journal. 2004, 87 (5): 3023-3034. 10.1529/biophysj.104.040824PubMed CentralView ArticlePubMedGoogle Scholar
- Leloup JC, Goldbeter A: Toward a detailed computational model for the mammalian circadian clock. Proc Natl Acad Sci U S A. 2003, 100 (12): 7051-7056. 10.1073/pnas.1132112100PubMed CentralView ArticlePubMedGoogle Scholar
- Aton SJ, Herzog ED: Come together, right...now: synchronization of rhythms in a mammalian circadian clock. Neuron. 2005, 48 (4): 531-534. 10.1016/j.neuron.2005.11.001PubMed CentralView ArticlePubMedGoogle Scholar
- Aton SJ, Colwell CS, Harmar AJ, Waschek J, Herzog ED: Vasoactive intestinal polypeptide mediates circadian rhythmicity and synchrony in mammalian clock neurons. Nat Neurosci. 2005, 8 (4): 476-483.PubMed CentralPubMedGoogle Scholar
- Liu AC, Welsh DK, Ko CH, Tran HG, Zhang EE, Priest AA, Buhr ED, Singer O, Meeker K, Verma IM, Doyle FJ, Takahashi JS, Kay SA: Intercellular coupling confers robustness against mutations in the SCN circadian clock network. Cell. 2007, 129 (3): 605-616. 10.1016/j.cell.2007.02.047PubMed CentralView ArticlePubMedGoogle Scholar
- Strogatz SH, Stewart I: Coupled oscillators and biological synchronization. Scientific American. 1993, 269 (6): 102-109.View ArticlePubMedGoogle Scholar
- Strogatz SH: From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators. Physica D. 2000, 143 (1-4): 1-20. 10.1016/S0167-2789(00)00094-4.View ArticleGoogle Scholar
- Somers D, Kopell N: Rapid Synchronization through Fast Threshold Modulation. Biological Cybernetics. 1993, 68 (5): 393-407. 10.1007/BF00198772View ArticlePubMedGoogle Scholar
- Winfree AT: Biological rhythms and the behavior of populations of coupled oscillators. Journal of theoretical biology. 1967, 16 (1): 15-42. 10.1016/0022-5193(67)90051-3View ArticlePubMedGoogle Scholar
- Kunz H, Achermann P: Simulation of circadian rhythm generation in the suprachiasmatic nucleus with locally coupled self-sustained oscillators. Journal of theoretical biology. 2003, 224 (1): 63-78. 10.1016/S0022-5193(03)00141-3View ArticlePubMedGoogle Scholar
- Achermann P, Kunz H: Modeling circadian rhythm generation in the suprachiasmatic nucleus with locally coupled self-sustained oscillators: phase shifts and phase response curves. Journal of biological rhythms. 1999, 14 (6): 460-468. 10.1177/074873099129001028View ArticlePubMedGoogle Scholar
- Gonze D, Bernard S, Waltermann C, Kramer A, Herzel H: Spontaneous synchronization of coupled circadian oscillators. Biophysical journal. 2005, 89 (1): 120-129. 10.1529/biophysj.104.058388PubMed CentralView ArticlePubMedGoogle Scholar
- Bernard S, Gonze D, Cajavec B, Herzel H, Kramer A: Synchronization-Induced Rhythmicity of Circadian Oscillators in the Suprachiasmatic Nucleus. PLoS Comput Biol. 2007, 3 (4): e68- 10.1371/journal.pcbi.0030068PubMed CentralView ArticlePubMedGoogle Scholar
- To TL, Henson MA, Herzog ED, Doyle FJ: A molecular model for intercellular synchronization in the mammalian circadian clock. Biophysical journal. 2007, 92 (11): 3792-3803. 10.1529/biophysj.106.094086PubMed CentralView ArticlePubMedGoogle Scholar
- Yamazaki S, Numano R, Abe M, Hida A, Takahashi R, Ueda M, Block GD, Sakaki Y, Menaker M, Tei H: Resetting central and peripheral circadian oscillators in transgenic rats. Science (New York, NY. 2000, 288 (5466): 682-685.View ArticleGoogle Scholar
- Nakamura W, Yamazaki S, Takasu NN, Mishima K, Block GD: Differential response of Period 1 expression within the suprachiasmatic nucleus. J Neurosci. 2005, 25 (23): 5481-5487. 10.1523/JNEUROSCI.0889-05.2005View ArticlePubMedGoogle Scholar
- Nagano M, Adachi A, Nakahama K, Nakamura T, Tamada M, Meyer-Bernstein E, Sehgal A, Shigeyoshi Y: An abrupt shift in the day/night cycle causes desynchrony in the mammalian circadian center. J Neurosci. 2003, 23 (14): 6141-6151.PubMedGoogle Scholar
- Best JD, Maywood ES, Smith KL, Hastings MH: Rapid resetting of the mammalian circadian clock. J Neurosci. 1999, 19 (2): 828-835.PubMedGoogle Scholar
- Wisor JP, Takahashi JS: Regulation of the vgf gene in the golden hamster suprachiasmatic nucleus by light and by the circadian clock. J Comp Neurol. 1997, 378 (2): 229-238. 10.1002/(SICI)1096-9861(19970210)378:2<229::AID-CNE6>3.0.CO;2-5View ArticlePubMedGoogle Scholar
- Prosser RA, Gillette MU: The mammalian circadian clock in the suprachiasmatic nuclei is reset in vitro by cAMP. J Neurosci. 1989, 9 (3): 1073-1081.PubMedGoogle Scholar
- Gribkoff VK, Pieschl RL, Wisialowski TA, van den Pol AN, Yocca FD: Phase shifting of circadian rhythms and depression of neuronal activity in the rat suprachiasmatic nucleus by neuropeptide Y: mediation by different receptor subtypes. J Neurosci. 1998, 18 (8): 3014-3022.PubMedGoogle Scholar
- Piggins HD, Antle MC, Rusak B: Neuropeptides phase shift the mammalian circadian pacemaker. J Neurosci. 1995, 15 (8): 5612-5622.PubMedGoogle Scholar
- Locke JCW, Southern MM, Kozma-Bognar L, Hibberd V, Brown PE, Turner MS, Millar AJ: Extension of a genetic network model by iterative experimentation and mathematical analysis. Mol Syst Biol. 2005, 1: 2005 0013- 10.1038/msb4100018PubMed CentralView ArticlePubMedGoogle Scholar
- Locke JCW, Millar AJ, Turner MS: Modelling genetic networks with noisy and varied experimental data: the circadian clock in Arabidopsis thaliana. Journal of theoretical biology. 2005, 234 (3): 383-393. 10.1016/j.jtbi.2004.11.038View ArticlePubMedGoogle Scholar
- King VM, Chahad-Ehlers S, Shen S, Harmar AJ, Maywood ES, Hastings MH: A hVIPR transgene as a novel tool for the analysis of circadian function in the mouse suprachiasmatic nucleus. The European journal of neuroscience. 2003, 17 (11): 822-832. 10.1046/j.1460-9568.2003.02487.x.View ArticleGoogle Scholar
- Bornholdt S: Systems biology. Less is more in modeling large genetic networks. Science (New York, NY. 2005, 310 (5747): 449-451.View ArticleGoogle Scholar
- Maywood ES, Reddy AB, Wong GK, O'Neill JS, O'Brien JA, McMahon DG, Harmar AJ, Okamura H, Hastings MH: Synchronization and maintenance of timekeeping in suprachiasmatic circadian clock cells by neuropeptidergic signaling. Curr Biol. 2006, 16 (6): 599-605. 10.1016/j.cub.2006.02.023View ArticlePubMedGoogle Scholar
- Cheng MY, Bullock CM, Li C, Lee AG, Bermak JC, Belluzzi J, Weaver DR, Leslie FM, Zhou QY: Prokineticin 2 transmits the behavioural circadian rhythm of the suprachiasmatic nucleus. Nature. 2002, 417 (6887): 405-410. 10.1038/417405aView ArticlePubMedGoogle Scholar
- Goodwin BC: Oscillatory behavior in enzymatic control processes. Adv Enzyme Regul. 1965, 3: 425-438. 10.1016/0065-2571(65)90067-1View ArticlePubMedGoogle Scholar
- Griffith JS: Mathematics of cellular control processes. I. Negative feedback to one gene. Journal of theoretical biology. 1968, 20 (2): 202-208. 10.1016/0022-5193(68)90189-6View ArticlePubMedGoogle Scholar
- Albus H, Vansteensel MJ, Michel S, Block GD, Meijer JH: A GABAergic mechanism is necessary for coupling dissociable ventral and dorsal regional oscillators within the circadian clock. Curr Biol. 2005, 15 (10): 886-893. 10.1016/j.cub.2005.03.051View ArticlePubMedGoogle Scholar
- Kuhlman SJ, Silver R, Le Sauter J, Bult-Ito A, McMahon DG: Phase resetting light pulses induce Per1 and persistent spike activity in a subpopulation of biological clock neurons. J Neurosci. 2003, 23 (4): 1441-1450.PubMed CentralPubMedGoogle Scholar
- Shigeyoshi Y, Taguchi K, Yamamoto S, Takekida S, Yan L, Tei H, Moriya T, Shibata S, Loros JJ, Dunlap JC, Okamura H: Light-induced resetting of a mammalian circadian clock is associated with rapid induction of the mPer1 transcript. Cell. 1997, 91 (7): 1043-1053. 10.1016/S0092-8674(00)80494-8View ArticlePubMedGoogle Scholar
- Forger DB, Peskin CS: Stochastic simulation of the mammalian circadian clock. Proc Natl Acad Sci U S A. 2005, 102 (2): 321-324. 10.1073/pnas.0408465102PubMed CentralView ArticlePubMedGoogle Scholar
- Vanselow K, Vanselow JT, Westermark PO, Reischl S, Maier B, Korte T, Herrmann A, Herzel H, Schlosser A, Kramer A: Differential effects of PER2 phosphorylation: molecular basis for the human familial advanced sleep phase syndrome (FASPS). Genes Dev. 2006, 20 (19): 2660-2672. 10.1101/gad.397006PubMed CentralView ArticlePubMedGoogle Scholar
- Gallego M, Eide EJ, Woolf MF, Virshup DM, Forger DB: An opposite role for tau in circadian rhythms revealed by mathematical modeling. Proc Natl Acad Sci U S A. 2006, 103 (28): 10618-10623. 10.1073/pnas.0604511103PubMed CentralView ArticlePubMedGoogle Scholar
- Antle MC, Silver R: Orchestrating time: arrangements of the brain circadian clock. Trends in neurosciences. 2005, 28 (3): 145-151. 10.1016/j.tins.2005.01.003View ArticlePubMedGoogle Scholar
- Quintero JE, Kuhlman SJ, McMahon DG: The biological clock nucleus: a multiphasic oscillator network regulated by light. J Neurosci. 2003, 23 (22): 8070-8076.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.