The Goodwin model: simulating the effect of light pulses on the circadian sporulation rhythm of Neurospora crassa

The seasons within: a theoretical perspective on photoperiodic entrainment and encoding

Circadian clocks are internal timing devices that have evolved as an adaption to the omnipresent natural 24 h rhythmicity of daylight intensity. Properties of the circadian system are photoperiod dependent. The phase of entrainment varies systematically with season. Plastic photoperiod-dependent re-arrangements in the mammalian circadian core pacemaker yield an internal representation of season. Output pathways of the circadian clock regulate photoperiodic responses such as flowering time in plants or hibernation in mammals. Here, we review the concepts of seasonal entrainment and photoperiodic encoding. We introduce conceptual phase oscillator models as their high level of abstraction, but, yet, intuitive interpretation of underlying parameters allows for a straightforward analysis of principles that determine entrainment characteristics. Results from this class of models are related and discussed in the context of more complex conceptual amplitude-phase oscillators as well as contextual molecular models that take into account organism, tissue, and cell-type-specific details.

Are circadian amplitudes and periods correlated? A new twist in the story

Three parameters are important to characterize a circadian and in general any biological clock: period, phase and amplitude. While circadian periods have been shown to correlate with entrainment phases, and clock amplitude influences the phase response of an oscillator to pulse-like zeitgeber signals, the co-modulations of amplitude and periods, which we term twist, have not been studied in detail. In this paper we define two concepts: parametric twist refers to amplitude-period correlations arising in ensembles of self-sustained, limit cycle clocks in the absence of external inputs, and phase space twist refers to the co-modulation of an individual clock's amplitude and period in response to external zeitgebers. Our findings show that twist influences the interaction of oscillators with the environment, facilitating entrainment, speeding upfastening recovery to pulse-like perturbations or modifying the response of an individual clock to coupling. This theoretical framework might be applied to understand the emerging properties of other oscillating systems.

Exploiting Information and Control Theory for Directing Gene Expression in Cell Populations

Microbial populations can adapt to adverse environmental conditions either by appropriately sensing and responding to the changes in their surroundings or by stochastically switching to an alternative phenotypic state. Recent data point out that these two strategies can be exhibited by the same cellular system, depending on the amplitude/frequency of the environmental perturbations and on the architecture of the genetic circuits involved in the adaptation process. Accordingly, several mitigation strategies have been designed for the effective control of microbial populations in different contexts, ranging from biomedicine to bioprocess engineering. Technically, such control strategies have been made possible by the advances made at the level of computational and synthetic biology combined with control theory. However, these control strategies have been applied mostly to synthetic gene circuits, impairing the applicability of the approach to natural circuits. In this review, we argue that it is possible to expand these control strategies to any cellular system and gene circuits based on a metric derived from this information theory, i.e., mutual information (MI). Indeed, based on this metric, it should be possible to characterize the natural frequency of any gene circuits and use it for controlling gene circuits within a population of cells.

Mathematical analysis of robustness of oscillations in models of the mammalian circadian clock

Circadian rhythms in a wide range of organisms are mediated by molecular mechanisms based on transcription-translation feedback. In this paper, we use bifurcation theory to explore mathematical models of genetic oscillators, based on Kim & Forger’s interpretation of the circadian clock in mammals. At the core of their models is a negative feedback loop whereby PER proteins (PER1 and PER2) bind to and inhibit their transcriptional activator, BMAL1. For oscillations to occur, the dissociation constant of the PER:BMAL1 complex, K^d, must be ≤ 0.04 nM, which is orders of magnitude smaller than a reasonable expectation of 1–10 nM for this protein complex. We relax this constraint by two modifications to Kim & Forger’s ‘single negative feedback’ (SNF) model: first, by introducing a multistep reaction chain for posttranscriptional modifications of Per mRNA and posttranslational phosphorylations of PER, and second, by replacing the first-order rate law for degradation of PER in the nucleus by a Michaelis-Menten rate law. These modifications increase the maximum allowable K^d to ~2 nM. In a third modification, we consider an alternative rate law for gene transcription to resolve an unrealistically large rate of Per2 transcription at very low concentrations of BMAL1. Additionally, we studied extensions of the SNF model to include a second negative feedback loop (involving REV-ERB) and a supplementary positive feedback loop (involving ROR). Contrary to Kim & Forger’s observations of these extended models, we find that, with our modifications, the supplementary positive feedback loop makes the oscillations more robust than observed in the models with one or two negative feedback loops. However, all three models are similarly robust when accounting for circadian rhythms (~24 h period) with K^d ≥ 1 nM. Our results provide testable predictions for future experimental studies.

The circadian rhythm aligns bodily functions to the day/night cycle and is important for our health. The rhythm originates from an intracellular molecular clock mechanism that mediates rhythmic gene expression. It is long understood that transcriptional negative feedback with sufficient time delay is key to generating circadian oscillations. However, some of the most widely cited mathematical models for the circadian clock suffer from problems of parameter ‘fragilities’. That is, sustained oscillations are possible only for physically unrealistic parameter values. A recent model by Kim & Forger nicely incorporates the inhibitory binding of PER proteins to their transcription activator BMAL1, but oscillations in the Kim-Forger model require a binding affinity between PER and BMAL1 that is orders of magnitude larger than observed binding affinities of protein complexes. To rectify this problem, we make several physiologically credible modifications to the Kim-Forger model, which allow oscillations to occur with more realistic binding affinities. The modified model is further extended to explore the potential roles of supplementary feedback loops in the mammalian clock mechanism. Ultimately, accurate models of the circadian clock will provide better predictive tools for chronotherapy and chrono-pharmacology studies.

Synthetic biology tools for engineering Goodwin oscillation in Trypanosoma brucei brucei

Kinetoplastid protozoa possess properties that are highly divergent from the mammalian, yeast and bacterial cells more commonly used in synthetic biology and represent a tantalisingly untapped source of bioengineering potential. Trypanosoma brucei brucei (T. b. brucei), an established model organism for studying the Kinetoplastida, is non-pathogenic to humans and provides an interesting test case for establishing synthetic biology in this phylogenetic class. To demonstrate further the tractability of Kinetoplastida to synthetic biology, we sought to construct and demonstrate a Goodwin oscillator, the simplest oscillatory gene network, in T. b. brucei for the first time. We report one completed iteration of the archetypal synthetic biology Design-Build-Test-Learn (DBTL) cycle; firstly, using Ab initio mathematical modelling of the behaviour a theoretical, oscillatory, trypanosomal synthetic gene network (SGN) to inform the design of a plasmid encoding that network. Once assembled, the plasmid was then used to generate a stable transfectant T. b. brucei cell line. To test the performance of the oscillatory SGN, a novel experimental setup was established to capture images of the fluorescent signal from motion-restricted live cells. Data captured were consistent with oscillatory behaviour of the SGN, with cellular fluorescence observed to oscillate with a period of 50 min, with varying amplitude and linear growth trend. This first DBTL cycle establishes a foundation for future cycles in which the SGN design and experimental monitoring setup can be further refined.

Mathematical Modeling in Circadian Rhythmicity

Circadian clocks are autonomous systems able to oscillate in a self-sustained manner in the absence of external cues, although such Zeitgebers are typically present. At the cellular level, the molecular clockwork consists of a complex network of interlocked feedback loops. This chapter discusses self-sustained circadian oscillators in the context of nonlinear dynamics theory. We suggest basic steps that can help in constructing a mathematical model and introduce how self-sustained generations can be modeled using ordinary differential equations. Moreover, we discuss how coupled oscillators synchronize among themselves or entrain to periodic signals. The development of mathematical models over the last years has helped to understand such complex network systems and to highlight the basic building blocks in which oscillating systems are built upon. We argue that, through theoretical predictions, the use of simple models can guide experimental research and is thus suitable to model biological systems qualitatively.

The Goodwin Oscillator and its Legacy

In the 1960's Brian Goodwin published a couple of mathematical models showing how feedback inhibition can lead to oscillations and discussed possible implications of this behaviour for the physiology of the cell. He also presented key ideas about the rich dynamics that may result from the coupling between such biochemical oscillators. Goodwin's work motivated a series of theoretical investigations aiming at identifying minimal mechanisms to generate limit cycle oscillations and deciphering design principles of biological oscillators. The three-variable Goodwin model (adapted by Griffith) can be seen as a core model for a large class of biological systems, ranging from ultradian to circadian clocks. We summarize here main ideas and results brought by Goodwin and review a couple of modeling works directly or indirectly inspired by Goodwin's findings.

Principles underlying the complex dynamics of temperature entrainment by a circadian clock

Autonomously oscillating circadian clocks resonate with daily environmental (zeitgeber) rhythms to organize physiology around the solar day. Although entrainment properties and mechanisms have been studied widely and in great detail for light-dark cycles, entrainment to daily temperature rhythms remains poorly understood despite that they are potent zeitgebers. Here we investigate the entrainment of the chronobiological model organism Neurospora crassa, subject to thermocycles of different periods and fractions of warm versus cold phases, mimicking seasonal variations. Depending on the properties of these thermocycles, regularly entrained rhythms, period-doubling (frequency demultiplication) but also irregular aperiodic behavior occurs. We demonstrate that the complex nonlinear phenomena of experimentally observed entrainment dynamics can be understood by molecular mathematical modeling.

A stochastic oscillator model simulates the entrainment of vertebrate cellular clocks by light

The circadian clock is a cellular mechanism that synchronizes various biological processes with respect to the time of the day. While much progress has been made characterizing the molecular mechanisms underlying this clock, it is less clear how external light cues influence the dynamics of the core clock mechanism and thereby entrain it with the light-dark cycle. Zebrafish-derived cell cultures possess clocks that are directly light-entrainable, thus providing an attractive laboratory model for circadian entrainment. Here, we have developed a stochastic oscillator model of the zebrafish circadian clock, which accounts for the core clock negative feedback loop, light input, and the proliferation of single-cell oscillator noise into population-level luminescence recordings. The model accurately predicts the entrainment dynamics observed in bioluminescent clock reporter assays upon exposure to a wide range of lighting conditions. Furthermore, we have applied the model to obtain refitted parameter sets for cell cultures exposed to a variety of pharmacological treatments and predict changes in single-cell oscillator parameters. Our work paves the way for model-based, large-scale screens for genetic or pharmacologically-induced modifications to the entrainment of circadian clock function.

Collective Oscillations in Coupled-Cell Systems

We investigate oscillations in coupled systems. The methodology is based on the Hopf bifurcation theorem and a condition extended from the Routh–Hurwitz criterion. Such a condition leads to locating the bifurcation values of the parameters. With such an approach, we analyze a single-cell system modeling the minimal genetic negative feedback loop and the coupled-cell system composed by these single-cell systems. We study the oscillatory properties for these systems and compare these properties between the model with Hill-type repression and the one with protein-sequestration-based repression. As the parameters move from the Hopf bifurcation value for single cells to the one for coupled cells, we compute the eigenvalues of the linearized systems to obtain the magnitude of the collective frequency when the periodic solution of the coupled-cell system is generated. Extending from this information on the parameter values, we further compute and compare the collective frequency for the coupled-cell system and the average frequency of the decoupled individual cells. To compare these scenarios with other biological oscillators, we perform parallel analysis and computations on a segmentation clock model.

An amplified derepression controller with multisite inhibition and positive feedback

How organisms are able to maintain robust homeostasis has in recent years received increased attention by the use of combined control engineering and kinetic concepts, which led to the discovery of robust controller motifs. While these motifs employ kinetic conditions showing integral feedback and homeostasis for step-wise perturbations, the motifs’ performance differ significantly when exposing them to time dependent perturbations. One type of controller motifs which are able to handle exponentially and even hyperbolically growing perturbations are based on derepression. In these controllers the compensatory reaction, which neutralizes the perturbation, is derepressed, i.e. its reaction rate is increased by the decrease of an inhibitor acting on the compensatory flux. While controllers in this category can deal well with different time-dependent perturbations they have the disadvantage that they break down once the concentration of the regulatory inhibitor becomes too low and the compensatory flux has gained its maximum value. We wondered whether it would be possible to bypass this restriction, while still keeping the advantages of derepression kinetics. In this paper we show how the inclusion of multisite inhibition and the presence of positive feedback loops lead to an amplified controller which is still based on derepression kinetics but without showing the breakdown due to low inhibitor concentrations. By searching for the amplified feedback motif in natural systems, we found it as a part of the plant circadian clock where it is highly interlocked with other feedback loops.

Robustness of a biomolecular oscillator to pulse perturbations

Biomolecular oscillators can function robustly in the presence of environmental perturbations, which can either be static or dynamic. While the effect of different circuit parameters and mechanisms on the robustness to steady perturbations has been investigated, the scenario for dynamic perturbations is relatively unclear. To address this, the authors use a benchmark three protein oscillator design - the repressilator - and investigate its robustness to pulse perturbations, computationally as well as use analytical tools of Floquet theory. They found that the metric provided by direct computations of the time it takes for the oscillator to settle after pulse perturbation is applied, correlates well with the metric provided by Floquet theory. They investigated the parametric dependence of the Floquet metric, finding that the parameters that increase the effective delay enhance robustness to pulse perturbation. They found that the structural changes such as increasing the number of proteins in a ring oscillator as well as adding positive feedback, both of which increase effective delay, facilitates such robustness. These results highlight such design principles, especially the role of delay, for designing an oscillator that is robust to pulse perturbation.

Hilbert transform-based time-series analysis of the circadian gene regulatory network

In this work, the authors propose the Hilbert transform (HT)‐based numerical method to analyse the time series of the circadian rhythms. They demonstrate the application of HT by taking both deterministic and stochastic time series that they get from the simulation of the fruit fly model Drosophila melanogaster and show how to extract the period, construct phase response curves, determine period sensitivity of the parameters to perturbations and build Arnold tongues to identify the regions of entrainment. They also derive a phase model that they numerically simulate to capture whether the circadian time series entrains to the forcing period completely (phase locking) or only partially (phase slips) or neither. They validate the phase model, and numerics with the experimental time series forced under different temperature cycles. Application of HT to the circadian time series appears to be a promising tool to extract the characteristic information about circadian rhythms.

Multi-scale modeling of the circadian modulation of learning and memory

We propose a multi-scale model to explain the time-of-day effects on learning and memory. We specifically model the circadian variation of hippocampus (HC) dependent long-term potentiation (LTP), depression (LTD), and the fear conditioning paradigm in amygdala. The model we built has both Goodwin type circadian gene regulatory network (GRN) and the conductance model of Morris-Lecar (ML) type to explain the spontaneous firing patterns (SFR) in suprachiasmatic nucleus (SCN). In the conductance model, we also include N-Methyl-D-aspartic acid receptor (NMDAR) to study the circadian dependent changes in LTP/LTD in hippocampus and include both NMDAR and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR) dynamics to explain the circadian modulation of fear conditioning paradigm in memory acquisition, recall, and extinction as seen in amygdala. Our multi-scale model captures the essential dynamics seen in the experiments and strongly supports the circadian time-of-the-day effects on learning and memory.

An Inactivation Switch Enables Rhythms in a Neurospora Clock Model

Autonomous endogenous time-keeping is ubiquitous across many living organisms, known as the circadian clock when it has a period of about 24 h. Interestingly, the fundamental design principle with a network of interconnected negative and positive feedback loops is conserved through evolution, although the molecular components differ. Filamentous fungus Neurospora crassa is a well-established chrono-genetics model organism to investigate the underlying mechanisms. The core negative feedback loop of the clock of Neurospora is composed of the transcription activator White Collar Complex (WCC) (heterodimer of WC1 and WC2) and the inhibitory element called FFC complex, which is made of FRQ (Frequency protein), FRH (Frequency interacting RNA Helicase) and CK1a (Casein kinase 1a). While exploring their temporal dynamics, we investigate how limit cycle oscillations arise and how molecular switches support self-sustained rhythms. We develop a mathematical model of 10 variables with 26 parameters to understand the interactions and feedback among WC1 and FFC elements in nuclear and cytoplasmic compartments. We performed control and bifurcation analysis to show that our novel model produces robust oscillations with a wild-type period of 22.5 h. Our model reveals a switch between WC1-induced transcription and FFC-assisted inactivation of WC1. Using the new model, we also study the possible mechanisms of glucose compensation. A fairly simple model with just three nonlinearities helps to elucidate clock dynamics, revealing a mechanism of rhythms' production. The model can further be utilized to study entrainment and temperature compensation.

Mathematical modeling of circadian rhythms

Circadian rhythms are endogenous ~24-hr oscillations usually entrained to daily environmental cycles of light/dark. Many biological processes and physiological functions including mammalian body temperature, the cell cycle, sleep/wake cycles, neurobehavioral performance, and a wide range of diseases including metabolic, cardiovascular, and psychiatric disorders are impacted by these rhythms. Circadian clocks are present within individual cells and at tissue and organismal levels as emergent properties from the interaction of cellular oscillators. Mathematical models of circadian rhythms have been proposed to provide a better understanding of and to predict aspects of this complex physiological system. These models can be used to: (a) manipulate the system in silico with specificity that cannot be easily achieved using in vivo and in vitro experimental methods and at lower cost, (b) resolve apparently contradictory empirical results, (c) generate hypotheses, (d) design new experiments, and (e) to design interventions for altering circadian rhythms. Mathematical models differ in structure, the underlying assumptions, the number of parameters and variables, and constraints on variables. Models representing circadian rhythms at different physiologic scales and in different species are reviewed to promote understanding of these models and facilitate their use. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Protein sequestration versus Hill-type repression in circadian clock models

Circadian (∼24 h) clocks are self-sustained endogenous oscillators with which organisms keep track of daily and seasonal time. Circadian clocks frequently rely on interlocked transcriptional-translational feedback loops to generate rhythms that are robust against intrinsic and extrinsic perturbations. To investigate the dynamics and mechanisms of the intracellular feedback loops in circadian clocks, a number of mathematical models have been developed. The majority of the models use Hill functions to describe transcriptional repression in a way that is similar to the Goodwin model. Recently, a new class of models with protein sequestration-based repression has been introduced. Here, the author discusses how this new class of models differs dramatically from those based on Hill-type repression in several fundamental aspects: conditions for rhythm generation, robust network designs and the periods of coupled oscillators. Consistently, these fundamental properties of circadian clocks also differ among Neurospora, Drosophila, and mammals depending on their key transcriptional repression mechanisms (Hill-type repression or protein sequestration). Based on both theoretical and experimental studies, this review highlights the importance of careful modelling of transcriptional repression mechanisms in molecular circadian clocks.

Modeling Reveals a Key Mechanism for Light-Dependent Phase Shifts of Neurospora Circadian Rhythms

Light shifts and synchronizes the phase of the circadian clock to daily environments, which is critical for maintaining the daily activities of an organism. It has been proposed that such light-dependent phase shifts are triggered by light-induced upregulation of a negative element of the core circadian clock (i.e., frq, Per1/2) in many organisms, including fungi. However, we find, using systematic mathematical modeling of the Neurospora crassa circadian clock, that the upregulation of the frq gene expression alone is unable to reproduce the observed light-dependent phase responses. Indeed, we find that the depression of the transcriptional activator white-collar-1, previously shown to be promoted by FRQ and VVD, is a key molecular mechanism for accurately simulating light-induced phase response curves for wild-type and mutant strains of Neurospora. Our findings elucidate specific molecular pathways that can be utilized to control phase resetting of circadian rhythms.

Asymmetry in Signal Oscillations Contributes to Efficiency of Periodic Systems

Oscillations are an important feature of cellular signaling that result from complex combinations of positive- and negative-feedback loops. The encoding and decoding mechanisms of oscillations based on amplitude and frequency have been extensively discussed in the literature in the context of intercellular and intracellular signaling. However, the fundamental questions of whether and how oscillatory signals offer any competitive advantages-and, if so, what-have not been fully answered. We investigated established oscillatory mechanisms and designed a study to analyze the oscillatory characteristics of signaling molecules and system output in an effort to answer these questions. Two classic oscillators, Goodwin and PER, were selected as the model systems, and corresponding no-feedback models were created for each oscillator to discover the advantage of oscillating signals. Through simulating the original oscillators and the matching no-feedback models, we show that oscillating systems have the capability to achieve better resource-to-output efficiency, and we identify oscillatory characteristics that lead to improved efficiency.

Hypothesis driven single cell dual oscillator mathematical model of circadian rhythms

Molecular mechanisms responsible for 24 h circadian oscillations, entrainment to external cues, encoding of day length and the time-of-day effects have been well studied experimentally. However, it is still debated from the molecular network point of view whether each cell in suprachiasmatic nuclei harbors two molecular oscillators, where one tracks dawn and the other tracks dusk activities. A single cell dual morning and evening oscillator was proposed by Daan et al., based on the molecular network that has two sets of similar non-redundant per1/cry1 and per2/cry2 circadian genes and each can independently maintain their endogenous oscillations. Understanding of dual oscillator dynamics in a single cell at molecular level may provide insight about the circadian mechanisms that encodes day length variations and its response to external zeitgebers. We present here a realistic dual oscillator model of circadian rhythms based on the series of hypotheses proposed by Daan et al., in which they conjectured that the circadian genes per1/cry1 track dawn while per2/cry2 tracks dusk and they together constitute the morning and evening oscillators (dual oscillator). Their hypothesis also provides explanations about the encoding of day length in terms of molecular mechanisms of per/cry expression. We frame a minimal mathematical model with the assumption that per1 acts a morning oscillator and per2 acts as an evening oscillator and to support and interpret this assumption we fit the model to the experimental data of per1/per2 circadian temporal dynamics, phase response curves (PRC's), and entrainment phenomena under various light-dark conditions. We also capture different patterns of splitting phenomena by coupling two single cell dual oscillators with neuropeptides vasoactive intestinal polypeptide (VIP) and arginine vasopressin (AVP) as the coupling agents and provide interpretation for the occurrence of splitting in terms of ME oscillators, though they are not required to explain the morning and evening oscillators. The proposed dual oscillator model based on Daan's hypothesis supports per1 and per2 playing the role of morning and evening oscillators respectively and this may be the first step towards the understanding of the core molecular mechanism responsible for encoding the day length.

Systems Biology-Derived Discoveries of Intrinsic Clocks

A systems approach to studying biology uses a variety of mathematical, computational, and engineering tools to holistically understand and model properties of cells, tissues, and organisms. Building from early biochemical, genetic, and physiological studies, systems biology became established through the development of genome-wide methods, high-throughput procedures, modern computational processing power, and bioinformatics. Here, we highlight a variety of systems approaches to the study of biological rhythms that occur with a 24-h period-circadian rhythms. We review how systems methods have helped to elucidate complex behaviors of the circadian clock including temperature compensation, rhythmicity, and robustness. Finally, we explain the contribution of systems biology to the transcription-translation feedback loop and posttranslational oscillator models of circadian rhythms and describe new technologies and "-omics" approaches to understand circadian timekeeping and neurophysiology.

Phase Sensitivity Analysis of Circadian Rhythm Entrainment

As a biological clock, circadian rhythms evolve to accomplish a stable (robust) entrainment to environmental cycles, of which light is the most obvious. The mechanism of photic entrainment is not known, but two models of entrainment have been proposed based on whether light has a continuous (parametric) or discrete (nonparametric) effect on the circadian pacemaker. A novel sensitivity analysis is developed to study the circadian entrainment in silico based on a limit cycle approach and applied to a model of Drosophila circadian rhythm. The comparative analyses of complete and skeleton photoperiods suggest a trade-off between the contribution of period modulation (parametric effect) and phase shift (nonparametric effect) in Drosophila circadian entrainment. The results also give suggestions for an experimental study to (in)validate the two models of entrainment.

Twilight, a Novel Circadian-Regulated Gene, Integrates Phototropism with Nutrient and Redox Homeostasis during Fungal Development

Phototropic regulation of circadian clock is important for environmental adaptation, organismal growth and differentiation. Light plays a critical role in fungal development and virulence. However, it is unclear what governs the intracellular metabolic response to such dark-light rhythms in fungi. Here, we describe a novel circadian-regulated Twilight (TWL) function essential for phototropic induction of asexual development and pathogenesis in the rice-blast fungus Magnaporthe oryzae. The TWL transcript oscillates during circadian cycles and peaks at subjective twilight. GFP-Twl remains acetylated and cytosolic in the dark, whereas light-induced phosphorylation (by the carbon sensor Snf1 kinase) drives it into the nucleus. The mRNA level of the transcription/repair factor TFB5, was significantly down regulated in the twl∆ mutant. Overexpression of TFB5 significantly suppressed the conidiation defects in the twl∆ mutant. Furthermore, Tfb5-GFP translocates to the nucleus during the phototropic response and under redox stress, while it failed to do so in the twl∆ mutant. Thus, we provide mechanistic insight into Twl-based regulation of nutrient and redox homeostasis in response to light during pathogen adaptation to the host milieu in the rice blast pathosystem.

The Goodwin model revisited: Hopf bifurcation, limit-cycle, and periodic entrainment

The three-variable Goodwin oscillator is a minimal model demonstrating the emergence of oscillations in simple biochemical feedback systems. As a prototypical oscillator, this model was extensively studied from a theoretical point of view and applied to various biological systems, including circadian clocks. Here, we reexamine this model, derive analytically the amplitude equation near the Hopf bifurcation and investigate the effect of a periodic modulation of the oscillator. In particular, we compare the entrainment performance when the free oscillator displays either self-sustained or damped oscillations. We discuss the results in the context of circadian oscillators.

Molecular mechanisms that regulate the coupled period of the mammalian circadian clock

In mammals, most cells in the brain and peripheral tissues generate circadian (∼24 h) rhythms autonomously. These self-sustained rhythms are coordinated and entrained by a master circadian clock in the suprachiasmatic nucleus (SCN). Within the SCN, the individual rhythms of each neuron are synchronized through intercellular signaling. One important feature of SCN is that the synchronized period is close to the population mean of cells' intrinsic periods. In this way, the synchronized period of the SCN stays close to the periods of cells in peripheral tissues. This is important because the SCN must entrain cells throughout the body. However, the mechanism that drives the period of the coupled SCN cells to the population mean is not known. We use mathematical modeling and analysis to show that the mechanism of transcription repression in the intracellular feedback loop plays a pivotal role in regulating the coupled period. Specifically, we use phase response curve analysis to show that the coupled period within the SCN stays near the population mean if transcriptional repression occurs via protein sequestration. In contrast, the coupled period is far from the mean if repression occurs through highly nonlinear Hill-type regulation (e.g., oligomer- or phosphorylation-based repression), as widely assumed in previous mathematical models. Furthermore, we find that the timescale of intercellular coupling needs to be fast compared to that of intracellular feedback to maintain the mean period. These findings reveal the important relationship between the intracellular transcriptional feedback loop and intercellular coupling. This relationship explains why transcriptional repression appears to occur via protein sequestration in multicellular organisms, mammals, and Drosophila, in contrast with the phosphorylation-based repression in unicellular organisms and syncytia. That is, transition to protein sequestration is essential for synchronizing multiple cells with a period close to the population mean (∼24 h).

The Goodwin model: behind the Hill function

The Goodwin model is a 3-variable model demonstrating the emergence of oscillations in a delayed negative feedback-based system at the molecular level. This prototypical model and its variants have been commonly used to model circadian and other genetic oscillators in biology. The only source of non-linearity in this model is a Hill function, characterizing the repression process. It was mathematically shown that to obtain limit-cycle oscillations, the Hill coefficient must be larger than 8, a value often considered unrealistic. It is indeed difficult to explain such a high coefficient with simple cooperative dynamics. We present here molecular models of the standard Goodwin model, based on single or multisite phosphorylation/dephosphorylation processes of a transcription factor, which have been previously shown to generate switch-like responses. We show that when the phosphorylation/dephosphorylation processes are fast enough, the limit-cycle obtained with a multisite phosphorylation-based mechanism is in very good quantitative agreement with the oscillations observed in the Goodwin model. Conditions in which the detailed mechanism is well approximated by the Goodwin model are given. A variant of the Goodwin model which displays sharp thresholds and relaxation oscillations is also explained by a double phosphorylation/dephosphorylation-based mechanism through a bistable behavior. These results not only provide rational support for the Goodwin model but also highlight the crucial role of the speed of post-translational processes, whose response curve are usually established at a steady state, in biochemical oscillators.

Mathematical modeling in chronobiology

Circadian clocks are autonomous oscillators entrained by external Zeitgebers such as light-dark and temperature cycles. On the cellular level, rhythms are generated by negative transcriptional feedback loops. In mammals, the suprachiasmatic nucleus (SCN) in the anterior part of the hypothalamus plays the role of the central circadian pacemaker. Coupling between individual neurons in the SCN leads to precise self-sustained oscillations even in the absence of external signals. These neuronal rhythms orchestrate the phasing of circadian oscillations in peripheral organs. Altogether, the mammalian circadian system can be regarded as a network of coupled oscillators. In order to understand the dynamic complexity of these rhythms, mathematical models successfully complement experimental investigations. Here we discuss basic ideas of modeling on three different levels (1) rhythm generation in single cells by delayed negative feedbacks, (2) synchronization of cells via external stimuli or cell-cell coupling, and (3) optimization of chronotherapy.

Entrainment of the suprachiasmatic nucleus network by a light-dark cycle

The synchronization of biological activity with the alternation of day and night (circadian rhythm) is performed in the brain by a group of neurons, constituting the suprachiasmatic nucleus (SCN). The SCN is divided into two subgroups of oscillating cells: the ventrolateral (VL) neurons, which are exposed to light (photic signal), and the dorsomedial (DM) neurons, which are coupled to the VL cells. When the coupling between these neurons is strong enough, the system synchronizes with the photic period. Upon increasing the cell coupling, the entrainment of the DM cells has been recently shown to occur via a very sharp (jumping) transition when the period of the photic input is larger than the intrinsic period of the cells. Here, we characterize this transition with a simple realistic model. We show that two bifurcations possibly lead to the disappearance of the endogenous mode. Using a mean-field model, we show that the jumping transition results from a supercritical Hopf-like bifurcation. This finding implies that both the period and strength of the stimulating photic signal, and the relative fraction of cells in the VL and DM compartments, are crucial in determining the synchronization of the system.

Fractional dynamics of globally slow transcription and its impact on deterministic genetic oscillation

In dynamical systems theory, a system which can be described by differential equations is called a continuous dynamical system. In studies on genetic oscillation, most deterministic models at early stage are usually built on ordinary differential equations (ODE). Therefore, gene transcription which is a vital part in genetic oscillation is presupposed to be a continuous dynamical system by default. However, recent studies argued that discontinuous transcription might be more common than continuous transcription. In this paper, by appending the inserted silent interval lying between two neighboring transcriptional events to the end of the preceding event, we established that the running time for an intact transcriptional event increases and gene transcription thus shows slow dynamics. By globally replacing the original time increment for each state increment by a larger one, we introduced fractional differential equations (FDE) to describe such globally slow transcription. The impact of fractionization on genetic oscillation was then studied in two early stage models – the Goodwin oscillator and the Rössler oscillator. By constructing a “dual memory” oscillator – the fractional delay Goodwin oscillator, we suggested that four general requirements for generating genetic oscillation should be revised to be negative feedback, sufficient nonlinearity, sufficient memory and proper balancing of timescale. The numerical study of the fractional Rössler oscillator implied that the globally slow transcription tends to lower the chance of a coupled or more complex nonlinear genetic oscillatory system behaving chaotically.

(Re)inventing the circadian feedback loop

For 20 years, researchers have thought that circadian clocks are defined by feedback loops of transcription and translation. The rediscovery of posttranslational circadian oscillators in diverse organisms forces us to rethink this paradigm. Meanwhile, the original "basic" feedback loops of canonical circadian clocks have swelled to include dozens of additional proteins acting in interlocked loops. We review several self-sustained clock mechanisms and propose that minimum requirements for diurnal timekeeping might be simpler than those of actual free-running circadian oscillators. Thus, complex mechanisms of circadian timekeeping might have evolved from random connections between unrelated feedback loops with independent but limited time-telling capability.

Photic desynchronization of two subgroups of circadian oscillators in a network model of the suprachiasmatic nucleus with dispersed coupling strengths

The suprachiasmatic nucleus (SCN) is the master circadian clock in mammals and is composed of thousands of neuronal oscillators expressing different intrinsic periods. These oscillators form a coupled network with a free-running period around 24 h in constant darkness and entrainable to the external light-dark cycle (T cycle). Coupling plays an important role in setting the period of the network and its range of entrainment. Experiments in rats have shown that two subgroups of oscillators within the SCN, a ventrolateral (VL) subgroup that receives photic input and a dorsomedial (DM) subgroup that is coupled to VL, can be desynchronized under a short (22-h) T cycle, with VL entrained to the cycle and DM free-running. We use a modified Goodwin model to understand how entrainment of the subgroups to short (22-h) and long (26-h) T cycles is influenced by light intensity, the proportion of neurons that receives photic input, and coupling heterogeneity. We find that the model's critical value for the proportion of photically-sensitive neurons is in accord with actual experimental estimates, while the model's inclusion of dispersed coupling can account for the experimental observation that VL and DM desynchronize more readily under the 22-h than under the 26-h T cycle. Heterogeneous intercellular coupling within the SCN is likely central to the generation of complex behavioral patterns.

A COMPARISON INVESTIGATION OF THE SIMPLEST MODELS OF CIRCADIAN RHYTHMS

We present three different simplest realistic models of the circadian pacemaker that can be interpreted in molecular terms: Goodwin model, piece-wise Goodwin model and delay model. The main purpose is to investigate how the common parameters in these three model affect the dynamical behavior. In particular, the existence and stability of periodic solutions of three models are studied, and the parameter space where the steady state becomes unstable is given which can be used to compare the effect of parameters on the oscillations of three models. Furthermore, the effects of parameters on the period of periodic solutions and phase response curves are investigated and discussed. Our results confirm that all three models can be used to simulate the realistic circadian rhythms. However, in practice, which model is chosen depends on the real data.

A simple negative interaction in the positive transcriptional feedback of a single gene is sufficient to produce reliable oscillations

Negative and positive transcriptional feedback loops are present in natural and synthetic genetic oscillators. A single gene with negative transcriptional feedback needs a time delay and sufficiently strong nonlinearity in the transmission of the feedback signal in order to produce biochemical rhythms. A single gene with only positive transcriptional feedback does not produce oscillations. Here, we demonstrate that this single-gene network in conjunction with a simple negative interaction can also easily produce rhythms. We examine a model comprised of two well-differentiated parts. The first is a positive feedback created by a protein that binds to the promoter of its own gene and activates the transcription. The second is a negative interaction in which a repressor molecule prevents this protein from binding to its promoter. A stochastic study shows that the system is robust to noise. A deterministic study identifies that the dynamics of the oscillator are mainly driven by two types of biomolecules: the protein, and the complex formed by the repressor and this protein. The main conclusion of this paper is that a simple and usual negative interaction, such as degradation, sequestration or inhibition, acting on the positive transcriptional feedback of a single gene is a sufficient condition to produce reliable oscillations. One gene is enough and the positive transcriptional feedback signal does not need to activate a second repressor gene. This means that at the genetic level an explicit negative feedback loop is not necessary. The model needs neither cooperative binding reactions nor the formation of protein multimers. Therefore, our findings could help to clarify the design principles of cellular clocks and constitute a new efficient tool for engineering synthetic genetic oscillators.

Modeling circadian clocks: From equations to oscillations

Circadian rhythms are generated at the cellular level by a small but tightly regulated genetic network. In higher eukaryotes, interlocked transcriptional-translational feedback loops form the core of this network, which ensures the activation of the right genes (proteins) at the right time of the day. Understanding how such a complex molecular network can generate robust, self-sustained oscillations and accurately responds to signals from the environment (such as light and temperature) is greatly helped by mathematical modeling. In the present paper we review some mathematical models for circadian clocks, ranging from abstract, phenomenological models to the most detailed molecular models. We explain how the equations are derived, highlighting the challenges for the modelers, and how the models are analyzed. We show how to compute bifurcation diagrams, entrainment, and phase response curves. In the subsequent paper, we discuss, through a selection of examples, how modeling efforts have contributed to a better understanding of the dynamics of the circadian regulatory network.

Modeling light adaptation in circadian clock: prediction of the response that stabilizes entrainment

Periods of biological clocks are close to but often different from the rotation period of the earth. Thus, the clocks of organisms must be adjusted to synchronize with day-night cycles. The primary signal that adjusts the clocks is light. In Neurospora, light transiently up-regulates the expression of specific clock genes. This molecular response to light is called light adaptation. Does light adaptation occur in other organisms? Using published experimental data, we first estimated the time course of the up-regulation rate of gene expression by light. Intriguingly, the estimated up-regulation rate was transient during light period in mice as well as Neurospora. Next, we constructed a computational model to consider how light adaptation had an effect on the entrainment of circadian oscillation to 24-h light-dark cycles. We found that cellular oscillations are more likely to be destabilized without light adaption especially when light intensity is very high. From the present results, we predict that the instability of circadian oscillations under 24-h light-dark cycles can be experimentally observed if light adaptation is altered. We conclude that the functional consequence of light adaptation is to increase the adjustability to 24-h light-dark cycles and then adapt to fluctuating environments in nature.

Mechanism of phase splitting in two coupled groups of suprachiasmatic-nucleus neurons

The phase-splitting behavior of coupled suprachiasmatic-nucleus neurons has been observed in many mammals, and its mechanism is still not completely understood. Based on our previous work [C. Gu, J. Wang, and Z. Liu, Phys. Rev. E 80, 030904(R) (2009)] on the free-running periods of neurons in the suprachiasmatic nucleus, we present here a modified Goodwin oscillator model to explain the mechanism of phase splitting. In contrast to the previous phase model, the modified Goodwin oscillator model contains the information on both the phase and amplitude and, thus, can show more features than the purely phase model, including all three behaviors of synchronization, phase splitting, and amplitude death and the distributed periodicity in the regions of synchronization and phase splitting, etc. An analytic phase model is extracted from the modified Goodwin oscillator model to explain the dependence of periodicity on the parameters. Moreover, both the modified Goodwin oscillator model and the analytic phase model show that the ensemble frequency can be enhanced or reduced by the time delay.

Circadian oscillators in eukaryotes

The biological clock, present in nearly all eukaryotes, has evolved such that organisms can adapt to our planet's rotation in order to anticipate the coming day or night as well as unfavorable seasons. As all modern high-precision chronometers, the biological clock uses oscillation as a timekeeping element. In this review, we describe briefly the discovery, historical development, and general properties of circadian oscillators. The issue of temperature compensation (TC) is discussed, and our present understanding of the underlying genetic and biochemical mechanisms in circadian oscillators are described with special emphasis on Neurospora crassa, mammals, and plants.

A design principle underlying the synchronization of oscillations in cellular systems

Biological oscillations are found ubiquitously in cells and are widely variable, with periods varying from milliseconds to months, and scales involving subcellular components to large groups of organisms. Interestingly, independent oscillators from different cells often show synchronization that is not the consequence of an external regulator. What is the underlying design principle of such synchronized oscillations, and can modeling show that the complex consequences arise from simple molecular or other interactions between oscillators? When biological oscillators are coupled with each other, we found that synchronization is induced when they are connected together through a positive feedback loop. Increasing the coupling strength of two independent oscillators shows a threshold beyond which synchronization occurs within a few cycles, and a second threshold where oscillation stops. The positive feedback loop can be composed of either double-positive (PP) or double-negative (NN) interactions between a node of each of the two oscillating networks. The different coupling structures have contrasting characteristics. In particular, PP coupling is advantageous with respect to stability of period and amplitude, when local oscillators are coupled with a short time delay, whereas NN coupling is advantageous for a long time delay. In addition, PP coupling results in more robust synchronized oscillations with respect to amplitude excursions but not period, with applied noise disturbances compared to NN coupling. However, PP coupling can induce a large fluctuation in the amplitude and period of the resulting synchronized oscillation depending on the coupling strength, whereas NN coupling ensures almost constant amplitude and period irrespective of the coupling strength. Intriguingly, we have also observed that artificial evolution of random digital oscillator circuits also follows this design principle. We conclude that a different coupling strategy might have been selected according to different evolutionary requirements.

How to achieve fast entrainment? The timescale to synchronization

Entrainment, where oscillators synchronize to an external signal, is ubiquitous in nature. The transient time leading to entrainment plays a major role in many biological processes. Our goal is to unveil the specific dynamics that leads to fast entrainment. By studying a generic model, we characterize the transient time to entrainment and show how it is governed by two basic properties of an oscillator: the radial relaxation time and the phase velocity distribution around the limit cycle. Those two basic properties are inherent in every oscillator. This concept can be applied to many biological systems to predict the average transient time to entrainment or to infer properties of the underlying oscillator from the observed transients. We found that both a sinusoidal oscillator with fast radial relaxation and a spike-like oscillator with slow radial relaxation give rise to fast entrainment. As an example, we discuss the jet-lag experiments in the mammalian circadian pacemaker.

Free-running period of neurons in the suprachiasmatic nucleus: Its dependence on the distribution of neuronal coupling strengths

The suprachiasmatic nucleus (SCN) pacemaker shows a free-running period ranging from 20 to 28 h for different species, which was usually explained from the angle of coupling strength. Based on the assumption of nonidentical coupling strengths in SCN, we find an alternative mechanism that the diversity of free-running period can be also caused by the distribution of coupling strengths. The free-running period is proportional to the average coupling strength and inverse proportional to the dispersion of couplings. Moreover, we present an analytic phase model to confirm our finding, which shows a solid foundation of our finding and opens a window to study the collective behaviors of SCN oscillators.

Entrainment of Circadian Programs

Of the three defining properties of circadian rhythmicity--persisting free-running rhythm, temperature compensation, and entrainment--the last is often poorly understood by many chronobiologists. This paper gives an overview of entrainment. Where have we been? Where are we now? Whence should we be going? Particular emphasis is given to a discussion of the Discrete vs. Continuous models for entrainment. We provide an integrated mechanism for entrainment from a limit-cycle perspective.

Simulating dark expressions and interactions of frq and wc-1 in the Neurospora circadian clock

Circadian rhythms are considered to play an essential part in the adaptation of organisms to their environments. The occurrence of circadian oscillations appears to be based on the presence of transcriptional-translational negative feedback loops. In Neurospora crassa, the protein FREQUENCY (FRQ) is part of such a negative feedback loop apparently by a direct interaction with its transcription factor WHITE COLLAR-1 (WC-1). Based on the observation that nuclear FRQ levels are significantly lower than nuclear WC-1 levels, it was suggested that FRQ would act more like a catalyst in inhibiting WC-1 rather than binding to WC-1 and making an inactive FRQ:WC-1 complex. Intrigued by this hypothesis, we constructed a model for the Neurospora circadian clock, which includes expression of the frq and the wc-1 genes and their possible interactions. The model suggests that even small amounts of nuclear FRQ-protein are capable of inhibiting frq transcription in a rhythmic manner by binding to WC-1 and promoting its degradation. Our model predicts the importance of a FRQ dependent degradation of WC-1 in closing the negative feedback loop. The model shows good agreement with experimental levels in nuclear and cytosolic FRQ and WC-1, their phase relationships, and several clock mutant phenotypes.

Robustness can evolve gradually in complex regulatory gene networks with varying topology

The topology of cellular circuits (the who-interacts-with-whom) is key to understand their robustness to both mutations and noise. The reason is that many biochemical parameters driving circuit behavior vary extensively and are thus not fine-tuned. Existing work in this area asks to what extent the function of any one given circuit is robust. But is high robustness truly remarkable, or would it be expected for many circuits of similar topology? And how can high robustness come about through gradual Darwinian evolution that changes circuit topology gradually, one interaction at a time? We here ask these questions for a model of transcriptional regulation networks, in which we explore millions of different network topologies. Robustness to mutations and noise are correlated in these networks. They show a skewed distribution, with a very small number of networks being vastly more robust than the rest. All networks that attain a given gene expression state can be organized into a graph whose nodes are networks that differ in their topology. Remarkably, this graph is connected and can be easily traversed by gradual changes of network topologies. Thus, robustness is an evolvable property. This connectedness and evolvability of robust networks may be a general organizational principle of biological networks. In addition, it exists also for RNA and protein structures, and may thus be a general organizational principle of all biological systems.

A molecular model for intercellular synchronization in the mammalian circadian clock

The mechanisms and consequences of synchrony among heterogeneous oscillators are poorly understood in biological systems. We present a multicellular, molecular model of the mammalian circadian clock that incorporates recent data implicating the neurotransmitter vasoactive intestinal polypeptide (VIP) as the key synchronizing agent. The model postulates that synchrony arises among circadian neurons because they release VIP rhythmically on a daily basis and in response to ambient light. Two basic cell types, intrinsically rhythmic pacemakers and damped oscillators, are assumed to arise from a distribution of Period gene transcription rates. Postsynaptic neurons show time-of-day dependent responses to VIP binding through a signaling cascade that activates Period mRNA transcription. The heterogeneous cell ensemble model self-synchronizes, entrains to ambient light-dark cycles, and desynchronizes in constant bright light or upon removal of VIP signaling. The degree of synchronicity observed depends on cell-specific features (e.g., mean and variability of parameters within the rhythm-generating loop), in addition to the more commonly studied effect of intercellular coupling strength. These simulations closely replicate experimental data and predict that heterogeneous oscillations (e.g., sustained, damped, and arrhythmic) arise from small differences in the molecular parameters between cells, that damped oscillators participate in entrainment and synchrony of the ensemble of cells, and that constant light desynchronizes oscillators by maximizing VIP release.