Synchronization and entrainment of coupled circadian oscillators

Light-induced synchronization modulation: Enhanced in weak coupling and attenuated in strong coupling among suprachiasmatic nucleus neurons

Existing experiments demonstrated that constant light has either enhancing or diminishing effects on the behavioral rhythms of mammals, sparking our intense interest in the underlying mechanisms of this paradoxical phenomenon. The influence of constant light on behavioral rhythms involves the regulation of collective neuronal behavior. The robustness of behavioral rhythms stems from the synchronization of neurons. In mammals, the synchronization among neurons is regulated by the suprachiasmatic nucleus (SCN) located in the hypothalamus. Neurons within the SCN exhibit significant heterogeneity. The intrinsic frequency and coupling strength are two fundamental characteristics determining the internal dynamics of the SCN. In this study, the Poincaré model was employed to investigate the impact of constant light on SCN neuronal dynamics. We found that constant light can modulate neuronal synchronization, a phenomenon tightly linked to the critical threshold value of coupling strength among the neurons. Specifically, under weak coupling, constant light enhances neuronal synchronization. Under strong coupling, constant light weakens synchronization among oscillators. Furthermore, higher light intensity results in lengthened periods and reduced amplitudes. Our findings elucidate important underlying mechanisms by which constant light either enhances or diminishes mammalian behavioral rhythms, and provide a new perspective for understanding the complex regulation network of circadian rhythms.

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.

Nonmonotonic enhancement of diversity-induced resonance in systems of mobile oscillators

Diversity is omnipresent in natural and synthetic extended systems, the phenomenon of diversity-induced resonance (DIR), wherein a moderate degree of the diversity can provoke an optimal collective response, provides researchers a brand-new strategy to amplify and utilize the weak signal. As yet the relevant advances focus mostly on the ideal situations where the interactions among elements are uncorrelated with the physical proximity of agents. Such a consideration overlooks interactions mediated by the motion of agents in space. Here, we investigate the signal response of an ensemble of spatial mobile heterogeneous bistable oscillators with two canonical interacting modes: dynamic and preset. The oscillators are considered as mass points and perform random walks in a two-dimensional square plane. Under the dynamic scheme, the oscillators can only interact with other oscillators within a fixed vision radius. For the preset circumstance, the interaction among oscillators occurs only when all of them are in a predefined region at the same moment. We find that the DIR can be obtained in both situations. Additionally, the strength of resonance nonmonotonically rises with respect to the increase of moving speed, and the optimal resonance is acquired by an intermediate magnitude of speed. Finally, we propose reduced equations to guarantee the occurrence of such mobility-optimized DIR on the basis of the fast switching approximation theory and also examine the robustness of such phenomenon through the excitable FitzHugh-Nagumo model and a different spatial motion mechanism. Our results reveal for the first time that the DIR can be optimized by the spatial mobility and thus has promising potential application in the communication of mobile agents.

A minimal model of peripheral clocks reveals differential circadian re-entrainment in aging

The mammalian circadian system comprises a network of endogenous oscillators, spanning from the central clock in the brain to peripheral clocks in other organs. These clocks are tightly coordinated to orchestrate rhythmic physiological and behavioral functions. Dysregulation of these rhythms is a hallmark of aging, yet it remains unclear how age-related changes lead to more easily disrupted circadian rhythms. Using a two-population model of coupled oscillators that integrates the central clock and the peripheral clocks, we derive simple mean-field equations that can capture many aspects of the rich behavior found in the mammalian circadian system. We focus on three age-associated effects that have been posited to contribute to circadian misalignment: attenuated input from the sympathetic pathway, reduced responsiveness to light, and a decline in the expression of neurotransmitters. We find that the first two factors can significantly impede re-entrainment of the clocks following perturbation, while a weaker coupling within the central clock does not affect the recovery rate. Moreover, using our minimal model, we demonstrate the potential of using the feed-fast cycle as an effective intervention to accelerate circadian re-entrainment. These results highlight the importance of peripheral clocks in regulating the circadian rhythm and provide fresh insights into the complex interplay between aging and the resilience of the circadian system.

Stability of ecosystems under oscillatory driving with frequency modulation

Consumer-resource cycles are widespread in ecosystems, and seasonal forcing is known to influence them profoundly. Typically, seasonal forcing perturbs an ecosystem with time-varying frequency; however, previous studies have explored the dynamics of such systems under oscillatory forcing with constant frequency. Studies of the effect of time-varying frequency on ecosystem stability are lacking. Here we investigate isolated and network models of a cyclic consumer-resource ecosystem with oscillatory driving subjected to frequency modulation. We show that frequency modulation can induce stability in the system in the form of stable synchronized solutions, depending on intrinsic model parameters and extrinsic modulation strength. The stability of synchronous solutions is determined by calculating the maximal Lyapunov exponent, which determines that the fraction of stable synchronous solution increases with an increase in the modulation strength. We also uncover intermittent synchronization when synchronous dynamics are intermingled with episodes of asynchronous dynamics. Using the phase-reduction method for the network model, we reduce the system into a phase equation that clearly distinguishes synchronous, intermittently synchronous, and asynchronous solutions. While investigating the role of network topology, we find that variation in rewiring probability has a negligible effect on the stability of synchronous solutions. This study deepens our understanding of ecosystems under seasonal perturbations.

A synthetic population-level oscillator in non-microfluidic environments

Synthetic oscillators have become a research hotspot because of their complexity and importance. The construction and stable operation of oscillators in large-scale environments are important and challenging. Here, we introduce a synthetic population-level oscillator in Escherichia coli that operates stably during continuous culture in non-microfluidic environments without the addition of inducers or frequent dilution. Specifically, quorum-sensing components and protease regulating elements are employed, which form delayed negative feedback to trigger oscillation and accomplish the reset of signals through transcriptional and post-translational regulation. We test the circuit in devices with 1 mL, 50 mL, 400 mL of medium, and demonstrate that the circuit could maintain stable population-level oscillations. Finally, we explore potential applications of the circuit in regulating cellular morphology and metabolism. Our work contributes to the design and testing of synthetic biological clocks that function in large populations.

Limits of entrainment of circadian neuronal networks

Circadian rhythmicity lies at the center of various important physiological and behavioral processes in mammals, such as sleep, metabolism, homeostasis, mood changes, and more. Misalignment of intrinsic neuronal oscillations with the external day-night cycle can disrupt such processes and lead to numerous disorders. In this work, we computationally determine the limits of circadian synchronization to external light signals of different frequency, duty cycle, and simulated amplitude. Instead of modeling circadian dynamics with generic oscillator models (e.g., Kuramoto-type), we use a detailed computational neuroscience model, which integrates biomolecular dynamics, neuronal electrophysiology, and network effects. This allows us to investigate the effect of small drug molecules, such as Longdaysin, and connect our results with experimental findings. To combat the high dimensionality of such a detailed model, we employ a matrix-free approach, while our entire algorithmic pipeline enables numerical continuation and construction of bifurcation diagrams using only direct simulation. We, thus, computationally explore the effect of heterogeneity in the circadian neuronal network, as well as the effect of the corrective therapeutic intervention of Longdaysin. Last, we employ unsupervised learning to construct a data-driven embedding space for representing neuronal heterogeneity.

Effects of time-varying habitat connectivity on metacommunity persistence

Network structure or connectivity patterns are critical in determining collective dynamics among interacting species in ecosystems. Conventional research on species persistence in spatial populations has focused on static network structure, though most real network structures change in time, forming time-varying networks. This raises the question, in metacommunities, how does the pattern of synchrony vary with temporal evolution in the network structure. The synchronous dynamics among species are known to reduce metacommunity persistence. Here we consider a time-varying metacommunity small-world network consisting of a chaotic three-species food chain oscillator in each patch or node. The rate of change in the network connectivity is determined by the natural frequency or its subharmonics of the constituent oscillator to allow sufficient time for the evolution of species in between successive rewirings. We find that over a range of coupling strengths and rewiring periods, even higher rewiring probabilities drive a network from asynchrony towards synchrony. Moreover, in networks with a small rewiring period, an increase in average degree (more connected networks) pushes the asynchronous dynamics to synchrony. On the other hand, in networks with a low average degree, a higher rewiring period drives the synchronous dynamics to asynchrony resulting in increased species persistence. Our results also follow the calculation of synchronization time and are robust across other ecosystem models. Overall, our study opens the possibility of developing temporal connectivity strategies to increase species persistence in ecological networks.

The emergence of polyglot entrainment responses to periodic inputs in vicinities of Hopf bifurcations in slow-fast systems

Several distinct entrainment patterns can occur in the FitzHugh-Nagumo (FHN) model under external periodic forcing. Investigating the FHN model under different types of periodic forcing reveals the existence of multiple disconnected 1:1 entrainment segments for constant, low enough values of the input amplitude when the unforced system is in the vicinity of a Hopf bifurcation. This entrainment structure is termed polyglot to distinguish it from the single 1:1 entrainment region (monoglot) structure typically observed in Arnold tongue diagrams. The emergence of polyglot entrainment is then explained using phase-plane analysis and other dynamical system tools. Entrainment results are investigated for other slow-fast systems of neuronal, circadian, and glycolytic oscillations. Exploring these models, we found that polyglot entrainment structure (multiple 1:1 regions) is observed when the unforced system is in the vicinity of a Hopf bifurcation and the Hopf point is located near a knee of a cubic-like nullcline.

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.

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.

Modeling synchronization in globally coupled oscillatory systems using model order reduction

We construct reduced order models for two classes of globally coupled multi-component oscillatory systems, selected as prototype models that exhibit synchronization. These are the Kuramoto model, considered both in its original formulation and with a suitable change of coordinates, and a model for the circadian clock. The systems of interest possess strong reduction properties, as their dynamics can be efficiently described with a low-dimensional set of coordinates. Specifically, the solution and selected quantities of interest are well approximated at the reduced level, and the reduced models recover the expected transition to synchronized states as the coupling strengths vary. Assuming that the interactions depend only on the averages of the system variables, the surrogate models ensure a significant computational speedup for large systems.

Circadian rhythms and the HPA axis: A systems view

The circadian timing system comprises a network of time-keeping clocks distributed across a living host whose responsibility is to allocate resources and distribute functions temporally to optimize fitness. The molecular structures generating these rhythms have evolved to accommodate the rotation of the earth in an attempt to primarily match the light/dark periods during the 24-hr day. To maintain synchrony of timing across and within tissues, information from the central clock, located in the suprachiasmatic nucleus, is conveyed using systemic signals. Leading among those signals are endocrine hormones, and while the hypothalamic-pituitary-adrenal axis through the release of glucocorticoids is a major pacesetter. Interestingly, the fundamental units at the molecular and physiological scales that generate local and systemic signals share critical structural properties. These properties enable time-keeping systems to generate rhythmic signals and allow them to adopt specific properties as they interact with each other and the external environment. The purpose of this review is to provide a broad overview of these structures, discuss their functional characteristics, and describe some of their fundamental properties as these related to health and disease. This article is categorized under: Immune System Diseases > Computational Models Immune System Diseases > Biomedical Engineering.

Neither per, nor tim1, nor cry2 alone are essential components of the molecular circadian clockwork in the Madeira cockroach

Circadian clocks control rhythms in physiology and behavior entrained to 24 h light-dark cycles. Despite of conserved general schemes, molecular circadian clockworks differ between insect species. With RNA interference (RNAi) we examined an ancient circadian clockwork in a basic insect, the hemimetabolous Madeira cockroach Rhyparobia maderae. With injections of double-stranded RNA (dsRNA) of cockroach period (Rm´per), timeless 1 (Rm´tim1), or cryptochrome 2 (Rm´cry2) we searched for essential components of the clock´s core negative feedback loop. Single injections of dsRNA of each clock gene into adult cockroaches successfully and permanently knocked down respective mRNA levels within ~two weeks deleting daytime-dependent mRNA rhythms for Rm´per and Rm´cry2. Rm´perRNAi or Rm´cry2RNAi affected total mRNA levels of both genes, while Rm´tim1 transcription was independent of both, also keeping rhythmic expression. Unexpectedly, circadian locomotor activity of most cockroaches remained rhythmic for each clock gene knockdown employed. It expressed weakened rhythms and unchanged periods for Rm´perRNAi and shorter periods for Rm´tim1RNAi and Rm´cry2RNAi.As a hypothesis of the cockroach´s molecular clockwork, a basic network of switched differential equations was developed to model the oscillatory behavior of clock cells expressing respective clock genes. Data were consistent with two synchronized main groups of coupled oscillator cells, a leading (morning) oscillator, or a lagging (evening) oscillator that couple via mutual inhibition. The morning oscillators express shorter, the evening oscillators longer endogenous periods based on core feedback loops with either PER, TIM1, or CRY2/PER complexes as dominant negative feedback of the clockwork. We hypothesize that dominant morning oscillator cells with shorter periods express PER, but not CRY2, or TIM1 as suppressor of clock gene expression, while two groups of evening oscillator cells with longer periods either comprise TIM1 or CRY2/PER suppressing complexes. Modelling suggests that there is an additional negative feedback next to Rm´PER in cockroach morning oscillator cells.

Coupling protocol of interlocked feedback oscillators in circadian clocks

Circadian rhythms (approx. 24 h) show the robustness of key oscillatory features such as phase, period and amplitude against external and internal variations. The robustness of Drosophila circadian clocks can be generated by interlocked transcriptional-translational feedback loops, where two negative feedback loops are coupled through mutual activations. The mechanisms by which such coupling protocols have survived out of many possible protocols remain to be revealed. To address this question, we investigated two distinct coupling protocols: activator-coupled oscillators (ACO) and repressor-coupled oscillators (RCO). We focused on the two coupling parameters: coupling dissociation constant and coupling time-delay. Interestingly, the ACO was able to produce anti-phase or morning-evening cycles, whereas the RCO produced in-phase ones. Deterministic and stochastic analyses demonstrated that the anti-phase ACO provided greater fluctuations in amplitude not only with respect to changes in coupling parameters but also to random parameter perturbations than the in-phase RCO. Moreover, the ACO deteriorated the entrainability to the day-night master clock, whereas the RCO produced high entrainability. Considering that the real, interlocked feedback loops have evolved as the ACO, instead of the RCO, we first proposed a hypothesis that the morning-evening or anti-phase cycle is more essential for Drosophila than achieving robustness and entrainability.

Astrocytic Modulation of Neuronal Activity in the Suprachiasmatic Nucleus: Insights from Mathematical Modeling

The suprachiasmatic nucleus (SCN) of the hypothalamus consists of a highly heterogeneous neuronal population networked together to allow precise and robust circadian timekeeping in mammals. While the critical importance of SCN neurons in regulating circadian rhythms has been extensively studied, the roles of SCN astrocytes in circadian system function are not well understood. Recent experiments have demonstrated that SCN astrocytes are circadian oscillators with the same functional clock genes as SCN neurons. Astrocytes generate rhythmic outputs that are thought to modulate neuronal activity through pre- and postsynaptic interactions. In this study, we developed an in silico multicellular model of the SCN clock to investigate the impact of astrocytes in modulating neuronal activity and affecting key clock properties such as circadian rhythmicity, period, and synchronization. The model predicted that astrocytes could alter the rhythmic activity of neurons via bidirectional interactions at tripartite synapses. Specifically, astrocyte-regulated extracellular glutamate was predicted to increase neuropeptide signaling from neurons. Consistent with experimental results, we found that astrocytes could increase the circadian period and enhance neural synchronization according to their endogenous circadian period. The impact of astrocytic modulation of circadian rhythm amplitude, period, and synchronization was predicted to be strongest when astrocytes had periods between 0 and 2 h longer than neurons. Increasing the number of neurons coupled to the astrocyte also increased its impact on period modulation and synchrony. These computational results suggest that signals that modulate astrocytic rhythms or signaling (e.g., as a function of season, age, or treatment) could cause disruptions in circadian rhythm or serve as putative therapeutic targets.

Circadian gating of dark-induced increases in chloroplast- and cytosolic-free calcium in Arabidopsis

Changes in the spatiotemporal concentration of free Ca²⁺ ([Ca²⁺ ]) in different organelles of the cell contribute to responses of plants to physiological and environmental stimuli. One example are [Ca²⁺ ] increases in the stroma of chloroplasts during light-to-dark transitions; however, the function and mechanisms responsible are unknown, in part because there is a disagreement in the literature concerning whether corresponding dark-induced changes in cytosolic [Ca²⁺ ] ([Ca²⁺ ]_cyt ) can be detected. We have measured changes in [Ca²⁺ ]_cyt upon darkness in addition to the already known dark-induced increases in [Ca²⁺ ]_stroma in the aerial part of the Arabidopsis thaliana plant. These [Ca²⁺ ]_cyt transients depend on the photoperiod and time of day, peaking at anticipated dusk, and are superimposed on daily 24 h oscillations in [Ca²⁺ ]_cyt . We also find that the magnitude of the dark-induced increases in Ca²⁺ in both the cytosol and chloroplasts are gated by the nuclear circadian oscillator. The modulation of the magnitude of dark-induced increases in [Ca²⁺ ]_stroma and [Ca²⁺ ]_cyt by transcriptional regulators in the nucleus that are part of the circadian oscillator demonstrates a new role for the circadian system in subcellular Ca²⁺ signalling, in addition to its role in driving circadian oscillations of [Ca²⁺ ] in the cytosol and chloroplasts.

Long-range dispersal promotes species persistence in climate extremes

Anthropogenic global warming in this century can act as a leading factor for large scale species extinctions in the near future. Species, in order to survive, need to develop dispersal strategies depending upon their environmental niche. Based on empirical evidence only a few previous studies have addressed how dispersal can evolve with changing temperature. However, for the analytical tractability, there is a need to develop an explicit model to ask how the temperature-dependent dispersal alters ecological dynamics. We investigate the persistence of species in a spatial ecological model, where dispersal is considered as a function of temperature. Spatial persistence is of major concern and dispersal is reasonably an important factor for extinction risk in the context of promoting synchrony. Our study yields how the temperature influences species decision of dispersal, resulting in either short-range or long-range dispersal. We examine synchronous or asynchronous behavior of species under their thermal dependence of dispersal. Moreover, we also analyze the transients to study the collective behavior of species away from their final or asymptotic dynamics. One of the key findings is at the most unfavorable environmental conditions long-range dispersal works out as the driving force for the persistence of species.

Synchronization and entrainment of metapopulations: A trade-off among time-induced heterogeneity, dispersal, and seasonal force

Demographic and environmental heterogeneities are prevalent across many natural systems. Earlier studies on metapopulation models have mostly considered heterogeneities either in the demographic parameters or in the interaction strength and topology between the spatially separated patches. In contrast, here we study the dynamics of a metapopulation model where each of the uncoupled patches has different periods of oscillations (period mismatch). We show different synchronization dynamics governed by both period mismatch and dispersal in neighboring patches. Indeed, we find both appearance and disappearance of phase synchronization, quasiperiodic oscillations, and period doubling of limit cycle. We also quantify the effect of seasonal variation (entrainment) and dispersal on species synchrony using phase-response curve and a synchrony measure, which thereof identify the influence of stochasticity on species persistence through trade-off mechanisms. Our results show that trade-offs among period mismatch, dispersal, and external force can drive entrained oscillations as well as asynchronous population dynamics that structure ecological communities.

Coordination of robust single cell rhythms in the Arabidopsis circadian clock via spatial waves of gene expression

The Arabidopsis circadian clock orchestrates gene regulation across the day/night cycle. Although a multiple feedback loop circuit has been shown to generate the 24-hr rhythm, it remains unclear how robust the clock is in individual cells, or how clock timing is coordinated across the plant. Here we examine clock activity at the single cell level across Arabidopsis seedlings over several days under constant environmental conditions. Our data reveal robust single cell oscillations, albeit desynchronised. In particular, we observe two waves of clock activity; one going down, and one up the root. We also find evidence of cell-to-cell coupling of the clock, especially in the root tip. A simple model shows that cell-to-cell coupling and our measured period differences between cells can generate the observed waves. Our results reveal the spatial structure of the plant clock and suggest that unlike the centralised mammalian clock, the Arabidopsis clock has multiple coordination points.

Increased persistence via asynchrony in oscillating ecological populations with long-range interaction

Understanding the influence of the structure of a dispersal network on the species persistence and modeling a realistic species dispersal in nature are two central issues in spatial ecology. A realistic dispersal structure which favors the persistence of interacting ecological systems was studied [M. D. Holland and A. Hastings, Nature (London) 456, 792 (2008)NATUAS0028-083610.1038/nature07395], where it was shown that a randomization of the structure of a dispersal network in a metapopulation model of prey and predator increases the species persistence via clustering, prolonged transient dynamics, and amplitudes of population fluctuations. In this paper, by contrast, we show that a deterministic network topology in a metapopulation can also favor asynchrony and prolonged transient dynamics if species dispersal obeys a long-range interaction governed by a distance-dependent power law. To explore the effects of power-law coupling, we take a realistic ecological model, namely, the Rosenzweig-MacArthur model in each patch (node) of the network of oscillators, and show that the coupled system is driven from synchrony to asynchrony with an increase in the power-law exponent. Moreover, to understand the relationship between species persistence and variations in power-law exponent, we compute a correlation coefficient to characterize cluster formation, a synchrony order parameter, and median predator amplitude. We further show that smaller metapopulations with fewer patches are more vulnerable to extinction as compared to larger metapopulations with a higher number of patches. We believe that the present work improves our understanding of the interconnection between the random network and the deterministic network in theoretical ecology.

The asymmetry of the entrainment range induced by the difference in intrinsic frequencies between two subgroups within the suprachiasmatic nucleus

The rhythms of physiological and behavioral activities in mammals, which are regulated by the main clock suprachiasmatic nucleus (SCN) in the brain, can not be only synchronized to the natural 24 h light-dark cycle, but also to cycles with artificial periods. The range of the artificial periods that the animal can be synchronized to is called entrainment range. In the absence of the light-dark cycle, the animal can also maintain the circadian rhythm with an endogenous period close to 24 h. Experiments found that the entrainment range is not symmetrical with respect to the endogenous period. In the present study, an explanation is given for the asymmetry based on a Kuramoto model which describes the neuronal network of the SCN. Our numerical simulations and theoretical analysis show that the asymmetry results from the difference in the intrinsic frequencies between two subgroups of the SCN, as well as the entrainment range is affected by the difference.

Dissociation between two subgroups of the suprachiasmatic nucleus affected by the number of damped oscillated neurons

In mammals, the main clock located in the suprachiasmatic nucleus (SCN) of the brain synchronizes the body rhythms to the environmental light-dark cycle. The SCN is composed of about 2×10^{4} neurons which can be classified into three oscillatory phenotypes: self-sustained oscillators, damped oscillators, and arrhythmic neurons. Exposed to an artificial external light-dark cycle with a period of 22h instead of 24h, two subgroups of the SCN can become desynchronized (dissociated). The ventrolateral (VL) subgroup receives photic input and is entrained to the external cycle and a dorsomedial (DM) subgroup oscillates with its endogenous (i.e., free running) period and is synchronized to the external light-dark cycle through coupling from the VL. In the present study, we examined the effects of damped oscillatory neurons on the dissociation between VL and DM under an external 22h cycle. We found that, with increasing numbers of damped oscillatory neurons located in the VL, the dissociation between the VL and DM emerges, but if these neurons are increasingly present in the DM the dissociation disappears. Hence, the damped oscillatory neurons in different subregions of the SCN play distinct roles in the dissociation between the two subregions of the SCN. This shows that synchrony between SCN subregions is affected by the number of damped oscillatory neurons and the location of these cells. We suggest that more knowledge on the number and the location of these cells may explain why some species do show a dissociation between the subregions and others do not, as the distribution of oscillatory types of neurons offers a plausible and novel candidate mechanism to explain heterogeneity.

The circadian rhythm induced by the heterogeneous network structure of the suprachiasmatic nucleus

In mammals, the master clock is located in the suprachiasmatic nucleus (SCN), which is composed of about 20 000 nonidentical neuronal oscillators expressing different intrinsic periods. These neurons are coupled through neurotransmitters to form a network consisting of two subgroups, i.e., a ventrolateral (VL) subgroup and a dorsomedial (DM) subgroup. The VL contains about 25% SCN neurons that receive photic input from the retina, and the DM comprises the remaining 75% SCN neurons which are coupled to the VL. The synapses from the VL to the DM are evidently denser than that from the DM to the VL, in which the VL dominates the DM. Therefore, the SCN is a heterogeneous network where the neurons of the VL are linked with a large number of SCN neurons. In the present study, we mimicked the SCN network based on Goodwin model considering four types of networks including an all-to-all network, a Newman-Watts (NW) small world network, an Erdös-Rényi (ER) random network, and a Barabási-Albert (BA) scale free network. We found that the circadian rhythm was induced in the BA, ER, and NW networks, while the circadian rhythm was absent in the all-to-all network with weak cellular coupling, where the amplitude of the circadian rhythm is largest in the BA network which is most heterogeneous in the network structure. Our finding provides an alternative explanation for the induction or enhancement of circadian rhythm by the heterogeneity of the network structure.

Impact of dispersed coupling strength on the free running periods of circadian rhythms

The dominant endogenous clock, named the suprachiasmatic nucleus (SCN), regulates circadian rhythms of behavioral and physiological activity in mammals. One of the main characteristics of the SCN is that the animal maintains a circadian rhythm with a period close to 24 h in the absence of a daily light-dark cycle (called the free running period). The free running period varies among species due to heterogeneity of the SCN network. Previous studies have shown that the heterogeneity in cellular coupling as well as in intrinsic neuronal periods shortens the free running period. Furthermore, as derived from experiments, one neuron's coupling strength is negatively associated with its period. It is unknown what the effects of this association between coupling strength and period are on the free running period and how the heterogeneity in coupling strength influences this free running period. In the present study we found that in the presence of a negative relationship between one neuron's coupling strength and its period, surprisingly, the dispersion of coupling strengths increases the free running period. Our present finding may shed new light on the understanding of the heterogeneous SCN network and provides an alternative explanation for the diversity of free running periods between species.

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.

Coupling-induced synchronization in multicellular circadian oscillators of mammals

In mammals, circadian rhythms are controlled by the neurons located in the suprachiasmatic nucleus (SCN) of the hypothalamus. Each neuron in the SCN contains an autonomous molecular clock. The fundamental question is how the individual cellular oscillators, expressing a wide range of periods, interact and assemble to achieve phase synchronization. Most of the studies carried out so far emphasize the crucial role of the periodicity imposed by the light-dark cycle in neuronal synchronization. However, in natural conditions, the interaction between the SCN neurons is non-negligible and coupling between cells in the SCN is achieved partly by neurotransmitters. In this paper, we use a model of nonidentical, globally coupled cellular clocks considered as Goodwin oscillators. We mainly study the synchronization induced by coupling from an analytical way. Our results show that the role of the coupling is to enhance the synchronization to the external forcing. The conclusion of this paper can help us better understand the mechanism of circadian rhythm.

Enhanced entrainability of genetic oscillators by period mismatch

Biological oscillators coordinate individual cellular components so that they function coherently and collectively. They are typically composed of multiple feedback loops, and period mismatch is unavoidable in biological implementations. We investigated the advantageous effect of this period mismatch in terms of a synchronization response to external stimuli. Specifically, we considered two fundamental models of genetic circuits: smooth and relaxation oscillators. Using phase reduction and Floquet multipliers, we numerically analysed their entrainability under different coupling strengths and period ratios. We found that a period mismatch induces better entrainment in both types of oscillator; the enhancement occurs in the vicinity of the bifurcation on their limit cycles. In the smooth oscillator, the optimal period ratio for the enhancement coincides with the experimentally observed ratio, which suggests biological exploitation of the period mismatch. Although the origin of multiple feedback loops is often explained as a passive mechanism to ensure robustness against perturbation, we study the active benefits of the period mismatch, which include increasing the efficiency of the genetic oscillators. Our findings show a qualitatively different perspective for both the inherent advantages of multiple loops and their essentiality.

Diversity and noise effects in a model of homeostatic regulation of the sleep-wake cycle

Recent advances in sleep neurobiology have allowed development of physiologically based mathematical models of sleep regulation that account for the neuronal dynamics responsible for the regulation of sleep-wake cycles and allow detailed examination of the underlying mechanisms. Neuronal systems in general, and those involved in sleep regulation in particular, are noisy and heterogeneous by their nature. It has been shown in various systems that certain levels of noise and diversity can significantly improve signal encoding. However, these phenomena, especially the effects of diversity, are rarely considered in the models of sleep regulation. The present paper is focused on a neuron-based physiologically motivated model of sleep-wake cycles that proposes a novel mechanism of the homeostatic regulation of sleep based on the dynamics of a wake-promoting neuropeptide orexin. Here this model is generalized by the introduction of intrinsic diversity and noise in the orexin-producing neurons, in order to study the effect of their presence on the sleep-wake cycle. A simple quantitative measure of the quality of a sleep-wake cycle is introduced and used to systematically study the generalized model for different levels of noise and diversity. The model is shown to exhibit a clear diversity-induced resonance: that is, the best wake-sleep cycle turns out to correspond to an intermediate level of diversity at the synapses of the orexin-producing neurons. On the other hand, only a mild evidence of stochastic resonance is found, when the level of noise is varied. These results show that disorder, especially in the form of quenched diversity, can be a key-element for an efficient or optimal functioning of the homeostatic regulation of the sleep-wake cycle. Furthermore, this study provides an example of a constructive role of diversity in a neuronal system that can be extended beyond the system studied here.

Transcription-based oscillator model for light-induced splitting as antiphase circadian gene expression in the suprachiasmatic nuclei

Daily locomotor patterns of a variety of organisms have been interpreted as driven by dual circadian oscillators. Yet, in mammals, cellular data have revealed many circadian oscillators in the bilateral suprachiasmatic nucleus (SCN). To test how large numbers of oscillators could respond to environmental cues as a pair of oscillators, the authors developed a computational model composed of 2 groups of oscillators with strong local interactions and with weaker coupling between the 2 groups. Unlike previous models that assumed that light affects the timing or polarity of coupling between a pair of oscillators, this simulation assumed that light increased the transcription rate of a clock gene and consequently altered circadian properties of individual cells. In constant dark, weak local (within each of the 2 groups) and distant (between group) coupling established in-phase oscillations and a typical single bout of daily activity. In constant light, local synchrony developed only if coupling was strong and resulted in antiphase synchrony between the 2 groups and bimodal daily activity reminiscent of split behavior. These numerical simulations thus showed that splitting behavior can develop with increased light intensity without structural changes in the coupling topology or sign. Instead, the authors propose that light changes intrinsic oscillator properties through the increase of maximal transcription rate of a clock gene, so that as light intensity increases, the output of the coupled network transitions from a single bout of activity through irregular beating to 2 bouts and, in bright constant light, arrhythmicity.

Stochastic properties of the plant circadian clock

Circadian clocks are gene regulatory networks whose role is to help the organisms to cope with variations in environmental conditions such as the day/night cycle. In this work, we explored the effects of molecular noise in single cells on the behaviour of the circadian clock in the plant model species Arabidopsis thaliana. The computational modelling language Bio-PEPA enabled us to give a stochastic interpretation of an existing deterministic model of the clock, and to easily compare the results obtained via stochastic simulation and via numerical solution of the deterministic model. First, the introduction of stochasticity in the model allowed us to estimate the unknown size of the system. Moreover, stochasticity improved the description of the available experimental data in several light conditions: noise-induced fluctuations yield a faster entrainment of the plant clock under certain photoperiods and are able to explain the experimentally observed dampening of the oscillations in plants under constant light conditions. The model predicts that the desynchronization between noisy oscillations in single cells contributes to the observed damped oscillations at the level of the cell population. Analysis of the phase, period and amplitude distributions under various light conditions demonstrated robust entrainment of the plant clock to light/dark cycles which closely matched the available experimental data.

Effect of network architecture on synchronization and entrainment properties of the circadian oscillations in the suprachiasmatic nucleus

In mammals, the suprachiasmatic nucleus (SCN) of the hypothalamus constitutes the central circadian pacemaker. The SCN receives light signals from the retina and controls peripheral circadian clocks (located in the cortex, the pineal gland, the liver, the kidney, the heart, etc.). This hierarchical organization of the circadian system ensures the proper timing of physiological processes. In each SCN neuron, interconnected transcriptional and translational feedback loops enable the circadian expression of the clock genes. Although all the neurons have the same genotype, the oscillations of individual cells are highly heterogeneous in dispersed cell culture: many cells present damped oscillations and the period of the oscillations varies from cell to cell. In addition, the neurotransmitters that ensure the intercellular coupling, and thereby the synchronization of the cellular rhythms, differ between the two main regions of the SCN. In this work, a mathematical model that accounts for this heterogeneous organization of the SCN is presented and used to study the implication of the SCN network topology on synchronization and entrainment properties. The results show that oscillations with larger amplitude can be obtained with scale-free networks, in contrast to random and local connections. Networks with the small-world property such as the scale-free networks used in this work can adapt faster to a delay or advance in the light/dark cycle (jet lag). Interestingly a certain level of cellular heterogeneity is not detrimental to synchronization performances, but on the contrary helps resynchronization after jet lag. When coupling two networks with different topologies that mimic the two regions of the SCN, efficient filtering of pulse-like perturbations in the entrainment pattern is observed. These results suggest that the complex and heterogeneous architecture of the SCN decreases the sensitivity of the network to short entrainment perturbations while, at the same time, improving its adaptation abilities to long term changes.

In order to adapt to their cycling environment, virtually all living organisms have developed an internal timer, the circadian clock. In mammals, the circadian pacemaker is composed of about 20,000 neurons, called the suprachiasmatic nucleus (SCN) located in the hypothalamus. The SCN receives light signals from the retina and controls peripheral circadian clocks to ensure the proper timing of physiological processes. In each SCN neuron, a genetic regulatory network enables the circadian expression of the clock genes, but individual dynamics are highly heterogeneous in dispersed cell culture: many cells present damped oscillations and the period of the oscillations varies from cell to cell. In addition, the neurotransmitters that ensure the intercellular coupling, and thereby the synchronization of the cellular rhythms, differ between the two main regions of the SCN. We present here a mathematical model that accounts for this heterogeneous organization of the SCN and study the implication of the network topology on synchronization and entrainment properties. Our results show that cellular heterogeneity may help the resynchronization after jet lag and suggest that the complex architecture of the SCN decreases the sensitivity of the network to short entrainment perturbations while, at the same time, improving its adaptation abilities to long term changes.