Global parameter search reveals design principles of the mammalian circadian clock

Background: 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. Results: 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. Conclusion: 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

Network Features of the Mammalian Circadian Clock

The mammalian circadian clock is a cell-autonomous system that drives oscillations in behavior and physiology in anticipation of daily environmental change. To assess the robustness of a human molecular clock, we systematically depleted known clock components and observed that circadian oscillations are maintained over a wide range of disruptions. We developed a novel strategy termed Gene Dosage Network Analysis (GDNA) in which small interfering RNA (siRNA)-induced dose-dependent changes in gene expression were used to build gene association networks consistent with known biochemical constraints. The use of multiple doses powered the analysis to uncover several novel network features of the circadian clock, including proportional responses and signal propagation through interacting genetic modules. We also observed several examples where a gene is up-regulated following knockdown of its paralog, suggesting the clock network utilizes active compensatory mechanisms rather than simple redundancy to confer robustness and maintain function. We propose that these network features act in concert as a genetic buffering system to maintain clock function in the face of genetic and environmental perturbation

Model-based Predictions of the Influence of Circadian Clock Genes Knock-Outs on the Cell Cycle

The existence of links between the cell division cycle and the circadian clock has been recently discovered. In this research report, we perform a small in silico analysis of how mutations on the clock genes composing the mammalian circadian clock affect the phenotype of the cell cycle. For this purpose, we use a coupled model of the mammalian circadian clock and the mammalian cell cycle where the latter one is entrained by the former one via a common protein kinase WEE1. Our in silico experiments exploit the modeling environment BIOCHAM

Mammalian Molecular Clocks

As a consequence of the Earth's rotation, almost all organisms experience day and night cycles within a 24-hr period. To adapt and synchronize biological rhythms to external daily cycles, organisms have evolved an internal time-keeping system. In mammals, the master circadian pacemaker residing in the suprachiasmatic nucleus (SCN) of the anterior hypothalamus generates circadian rhythmicity and orchestrates numerous subsidiary local clocks in other regions of the brain and peripheral tissues. Regardless of their locations, these circadian clocks are cell-autonomous and self-sustainable, implicating rhythmic oscillations in a variety of biochemical and metabolic processes. A group of core clock genes provides interlocking molecular feedback loops that drive the circadian rhythm even at the single-cell level. In addition to the core transcription/translation feedback loops, post-translational modifications also contribute to the fine regulation of molecular circadian clocks. In this article, we briefly review the molecular mechanisms and post-translational modifications of mammalian circadian clock regulation. We also discuss the organization of and communication between central and peripheral circadian oscillators of the mammalian circadian clock

Dynamics at the serine loop underlie differential affinity of cryptochromes for CLOCK:BMAL1 to control circadian timing.

Mammalian circadian rhythms are generated by a transcription-based feedback loop in which CLOCK:BMAL1 drives transcription of its repressors (PER1/2, CRY1/2), which ultimately interact with CLOCK:BMAL1 to close the feedback loop with ~24 hr periodicity. Here we pinpoint a key difference between CRY1 and CRY2 that underlies their differential strengths as transcriptional repressors. Both cryptochromes bind the BMAL1 transactivation domain similarly to sequester it from coactivators and repress CLOCK:BMAL1 activity. However, we find that CRY1 is recruited with much higher affinity to the PAS domain core of CLOCK:BMAL1, allowing it to serve as a stronger repressor that lengthens circadian period. We discovered a dynamic serine-rich loop adjacent to the secondary pocket in the photolyase homology region (PHR) domain that regulates differential binding of cryptochromes to the PAS domain core of CLOCK:BMAL1. Notably, binding of the co-repressor PER2 remodels the serine loop of CRY2, making it more CRY1-like and enhancing its affinity for CLOCK:BMAL1

A mammalian circadian clock model incorporating daytime expression elements

Models of the mammalian clock have traditionally been based around two feedback loops-the self-repression of Per/Cry by interfering with activation by BMAL/CLOCK, and the repression of Bmal/Clock by the REV-ERB proteins. Recent experimental evidence suggests that the D-box, a transcription factor binding site associated with daytime expression, plays a larger role in clock function than has previously been understood. We present a simplified clock model that highlights the role of the D-box and illustrate an approach for finding maximum-entropy ensembles of model parameters, given experimentally imposed constraints. Parameter variability can be mitigated using prior probability distributions derived from genome-wide studies of cellular kinetics. Our model reproduces predictions concerning the dual regulation of Cry1 by the D-box and Rev-ErbA/ROR response element (RRE) promoter elements and allows for ensemble-based predictions of phase response curves (PRCs). Nonphotic signals such as Neuropeptide Y (NPY) may act by promoting Cry1 expression, whereas photic signals likely act by stimulating expression from the E/E' box. Ensemble generation with parameter probability restraints reveals more about a model's behavior than a single optimal parameter set

Feedback Loops of the Mammalian Circadian Clock Constitute Repressilator

Mammals evolved an endogenous timing system to coordinate their physiology and behaviour to the 24h period of the solar day. While it is well accepted that circadian rhythms are generated by intracellular transcriptional feedback loops, it is still debated which network motifs are necessary and sufficient for generating self-sustained oscillations. Here, we systematically explore a data-based circadian oscillator model with multiple negative and positive feedback loops and identify a series of three subsequent inhibitions known as “repressilator” as a core element of the mammalian circadian oscillator. The central role of the repressilator motif is consistent with time-resolved ChIP- seq experiments of circadian clock transcription factors and loss of rhythmicity in core clock gene knockouts

A Multi-Layered Study on Harmonic Oscillations in Mammalian Genomics and Proteomics

Cellular, organ, and whole animal physiology show temporal variation predominantly featuring 24-h (circadian) periodicity. Time-course mRNA gene expression profiling in mouse liver showed two subsets of genes oscillating at the second (12-h) and third (8-h) harmonic of the prime (24-h) frequency. The aim of our study was to identify specific genomic, proteomic, and functional properties of ultradian and circadian subsets. We found hallmarks of the three oscillating gene subsets, including different (i) functional annotation, (ii) proteomic and electrochemical features, and (iii) transcription factor binding motifs in upstream regions of 8-h and 12-h oscillating genes that seemingly allow the link of the ultradian gene sets to a known circadian network. Our multifaceted bioinformatics analysis of circadian and ultradian genes suggests that the different rhythmicity of gene expression impacts physiological outcomes and may be related to transcriptional, translational and post-translational dynamics, as well as to phylogenetic and evolutionary components