- Research article
- Open Access
Global parameter search reveals design principles of the mammalian circadian clock
BMC Systems Biologyvolume 2, Article number: 22 (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
We carried out a global parameter screen of a previously developed simplified model of the SCN  in order to study in depth the design principles of the mammalian circadian clock. In this model (Figure 1A), within each clock cell, a clock gene mRNA (X) produces a clock protein (Y) which goes onto activate a transcriptional repressor (Z). Z represses the transcription of the clock gene X, closing a negative feedback loop. It is assumed that clock cells synthesize a neuropeptide denoted by V and that production is induced by the activation of the clock gene (X). The neurotransmitter release is assumed to be fast in comparison to the timescale of the oscillations (~24 h). The neurotransmitter levels thus become homogeneous, resulting in an average neurotransmitter level, or a mean field F. The mean field level of the neuropeptide then acts as a synchronising factor by activating X transcription both in an autocrine and paracrine manner (see Figure 1B for model equations, and Additional file 3 for detailed description). In reality, the SCN is much more complex than this simplified model. The clock circuit within each cell contains multiple feedback loops, and the strength of the neurotransmitter feedback appears to vary from neuron to neuron. For example, only a subsection of the actual SCN neurons express the VPAC2 receptor for the neurotransmitter VIP . However, by analysing in depth an idealized model such as this, it is often possible to gain important insights about putative mechanisms .
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.
Figure 2A shows a simulation using our optimal parameter set of X, Y and Z levels, as well as neurotransmitter levels, for a single neuron under constant dark conditions. X levels peak 3–4 h before Y, and there is another 2–3 h delay between Y levels and Z similar to experimental findings. The neurotransmitter has a similar profile as X, but is slightly delayed. Figure 2B shows X levels for 5 randomly selected oscillators from a set of 200 coupled oscillators. Each individual oscillator has its parameters rescaled by a scaling factor to represent the fluctuations in period and amplitude seen experimentally in isolated neurons . The oscillators are extremely well synchronised, with the mean field of the transmitter oscillating with a period of 24 h (Figure 2C). The oscillators remain synchronised even when the individual parameters are rescaled to represent 5 fold larger fluctuations in period (Figure 2D). A simulation of 10000 coupled oscillators similar to the total number of neurons in the SCN was also extremely well synchronised (data not shown), but for this study we concentrate on a more numerically tractable number of 200 oscillators for in depth analysis.
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
We went on to investigate the entrainment properties of our optimal parameter set, to see whether we were able to simulate the fast re-entrainment of the SCN seen experimentally. Figure 3A shows that the mean field of the neurotransmitter is stably entrained in a 12:12 LD cycle with light intensity of 1.0, showing no signs of quasi-periodicity as previously reported . We then tested the response of our network with the optimal parameter set to 'jet lag' type experiments, where the LD cycle is advanced or delayed by 6 h, and the resulting phase difference to an unperturbed system is measured. Excitingly, using a light intensity of 1 for the LD cycles which gives only a maximal phase shift of 2.5 h if applied for 1 h (Figure 3B), the mean field oscillations show very fast re-entrainment similar to that seen experimentally: After 1–2 days, for both a 6 h advance and delay, (Figure 3C), there is virtually no difference in phase between the perturbed and unperturbed systems, as seen in experiment [20, 21]. We see a slight 'overshoot' effect after a 6 h advance in the light cycle, where the advance in phase of the clock becomes larger than 6 h, then retards to the appropriate steady-state advance over the following days, which is qualitatively similar to the overshoot seen experimentally . The overshoot size was dependent on the light level used in our simulation. Using a lower light level, that results in a maximum phase shift of 0.5 hours if applied for one hour, the system takes around 6 days to re-entrain (Additional file 1: Supplementary Figure S4), with no overshoot effect seen. The actual SCN shows a more pronounced and prolonged desynchronization after phase advances than delays compared to our simulations [21, 22], and also shows heterogeneity in the phase resetting properties of individual neurons . It is not possible to simulate this complexity using the simplified SCN model we study. However, these results do show that the experimentally seen extremely fast phase resetting of the SCN after a shift in LD cycle is possible when the neuropeptide positive feedback loop drives the oscillation.
Light activated pacemaker cells can synchronize the SCN
It has been shown that only a subset of the SCN neurons receives direct light input from the retina. We tested whether the results we observe still hold true under the more physiological assumption that only a subset of the simulated neurons are light activated . Figure 4A shows the results of a 'jet lag' type experiment where only 25% of the cells are light activated with a light intensity of 1.0. Compared to a simulation where all the neurons are light activated (Figure 3), re-entrainment takes longer, about 4 days. However, if we double the strength of the light activation into the SCN, (Figure 4B), we regain fast re-entrainment of 1–2 days. This fits with experimental data showing that PER1 is strongly light activated in light responsive cells [37, 38].
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.
Dunlap JC: Molecular bases for circadian clocks. Cell. 1999, 96 (2): 271-290. 10.1016/S0092-8674(00)80566-8
Hastings MH, Herzog ED: Clock genes, oscillators, and cellular networks in the suprachiasmatic nuclei. J Biol Rhythms. 2004, 19 (5): 400-413. 10.1177/0748730404268786
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.022
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-7
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.2036281100
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.040824
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.1132112100
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.001
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.
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.047
Strogatz SH, Stewart I: Coupled oscillators and biological synchronization. Scientific American. 1993, 269 (6): 102-109.
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.
Somers D, Kopell N: Rapid Synchronization through Fast Threshold Modulation. Biological Cybernetics. 1993, 68 (5): 393-407. 10.1007/BF00198772
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-3
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-3
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/074873099129001028
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.058388
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.0030068
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.094086
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.
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.2005
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.
Best JD, Maywood ES, Smith KL, Hastings MH: Rapid resetting of the mammalian circadian clock. J Neurosci. 1999, 19 (2): 828-835.
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-5
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.
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.
Piggins HD, Antle MC, Rusak B: Neuropeptides phase shift the mammalian circadian pacemaker. J Neurosci. 1995, 15 (8): 5612-5622.
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/msb4100018
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.038
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.
Bornholdt S: Systems biology. Less is more in modeling large genetic networks. Science (New York, NY. 2005, 310 (5747): 449-451.
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.023
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/417405a
Goodwin BC: Oscillatory behavior in enzymatic control processes. Adv Enzyme Regul. 1965, 3: 425-438. 10.1016/0065-2571(65)90067-1
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-6
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.051
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.
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-8
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.0408465102
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.397006
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.0604511103
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.003
Quintero JE, Kuhlman SJ, McMahon DG: The biological clock nucleus: a multiphasic oscillator network regulated by light. J Neurosci. 2003, 23 (22): 8070-8076.
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.
JCWL, PW and HH conceived and designed the study. JCWL and PW performed the numerical experiments. JCWL, PW, HH, and AK analyzed the data. All authors contributed to writing the paper, and read and approved the final manuscript.