The synchronization of neuronal oscillators determined by the directed network structure of the suprachiasmatic nucleus under different photoperiods

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.

The Physiological and Pharmacological Significance of the Circadian Timing of the HPA Axis: A Mathematical Modeling Approach

As a potent endogenous regulator of homeostasis, the circadian time-keeping system synchronizes internal physiology to periodic changes in the external environment to enhance survival. Adapting endogenous rhythms to the external time is accomplished hierarchically with the central pacemaker located in the suprachiasmatic nucleus (SCN) signaling the hypothalamus-pituitary-adrenal (HPA) axis to release hormones, notably cortisol, which help maintain the body's circadian rhythm. Given the essential role of HPA-releasing hormones in regulating physiological functions, including immune response, cell cycle, and energy metabolism, their daily variation is critical for the proper function of the circadian timing system. In this review, we focus on cortisol and key fundamental properties of the HPA axis and highlight their importance in controlling circadian dynamics. We demonstrate how systems-driven, mathematical modeling of the HPA axis complements experimental findings, enhances our understanding of complex physiological systems, helps predict potential mechanisms of action, and elucidates the consequences of circadian disruption. Finally, we outline the implications of circadian regulation in the context of personalized chronotherapy. Focusing on the chrono-pharmacology of synthetic glucocorticoids, we review the challenges and opportunities associated with moving toward personalized therapies that capitalize on circadian rhythms.

Optimal phase-selective entrainment of heterogeneous oscillator ensembles

We develop a framework to design optimal entrainment signals that entrain an ensemble of heterogeneous nonlinear oscillators, described by phase models, at desired phases. We explicitly take into account heterogeneity in both oscillation frequency and the type of oscillators characterized by different Phase Response Curves. The central idea is to leverage the Fourier series representation of periodic functions to decode a phase-selective entrainment task into a quadratic program. We demonstrate our approach using a variety of phase models, where we entrain the oscillators into distinct phase patterns. Also, we show how the generalizability gained from our formulation enables us to meet a wide range of design objectives and constraints, such as minimum-power, fast entrainment, and charge-balanced controls.

Generation and Disruption of Circadian Rhythms in the Suprachiasmatic Nucleus: A Core-Shell Model

We focus our research on how the core-shell organization controls behavior of the suprachiasmatic nucleus (SCN), how the core and shell are synchronized to the environment, what impact they have on the behavior of the SCN under different lighting conditions, and what mechanisms disrupt synchronization. To this end, we use a reduced Kuramoto model, with parameters inferred from experimental observations and calibrated for mice, and perform a detailed comparison between the model and experimental data under light-dark (LD), dark-dark (DD), and light-light (LL) conditions. The operating limits of free-running and entrained SCN activity under symmetric LD cycles are analyzed, with particular focus on the phenomena of anticipation and dissociation. Results reveal that the core-shell organization of the SCN enables anticipation of future events over circadian cycles. The model predicts the emergence of a second (dissociated) rhythm for large and small LD periods. Our results are in good qualitative and quantitative agreement with experimental observations of circadian dissociation. We further describe SCN activity under LL conditions and show that our model satisfies Aschoff’s first rule, according to which the endogenous free-running circadian period observed under complete darkness will shorten in diurnal animals and lengthen in nocturnal animals under constant light. Our results strongly suggest that the Kuramoto model captures essential features of synchronization and entrainment in the SCN. Moreover, our approach is easily extendible to an arbitrary number of groups, with dynamics described by explicit equations for the group phase and synchronization index. Viewed together, the reduced Kuramoto model presents itself as a useful tool for exploring open problems in the study of circadian rhythms, one that can account for evolving views of the circadian system’s organization, including peripheral clocks and inter-hemispheric interaction, and can be translated to other nocturnal and diurnal animals, including humans.

Light-induced synchronization of the SCN coupled oscillators and implications for entraining the HPA axis

The suprachiasmatic nucleus (SCN) synchronizes the physiological rhythms to the external light-dark cycle and tunes the dynamics of circadian rhythms to photoperiod fluctuations. Changes in the neuronal network topologies are suggested to cause adaptation of the SCN in different photoperiods, resulting in the broader phase distribution of neuron activities in long photoperiods (LP) compared to short photoperiods (SP). Regulated by the SCN output, the level of glucocorticoids is elevated in short photoperiod, which is associated with peak disease incidence. The underlying coupling mechanisms of the SCN and the interplay between the SCN and the HPA axis have yet to be fully elucidated. In this work, we propose a mathematical model including a multiple-cellular SCN compartment and the HPA axis to investigate the properties of the circadian timing system under photoperiod changes. Our model predicts that the probability-dependent network is more energy-efficient than the distance-dependent network. Coupling the SCN network by intra-subpopulation and inter-subpopulation forces, we identified the negative correlation between robustness and plasticity of the oscillatory network. The HPA rhythms were predicted to be strongly entrained to the SCN rhythms with a pro-inflammatory high-amplitude glucocorticoid profile under SP. The fast temporal topology switch of the SCN network was predicted to enhance synchronization when the synchronization is not complete. These synchronization and circadian dynamics alterations might govern the seasonal variation of disease incidence and its symptom severity.

Network Structure of the Master Clock Is Important for Its Primary Function

A master clock located in the suprachiasmatic nucleus (SCN) regulates the circadian rhythm of physiological and behavioral activities in mammals. The SCN has two main functions in the regulation: an endogenous clock produces the endogenous rhythmic signal in body rhythms, and a calibrator synchronizes the body rhythms to the external light-dark cycle. These two functions have been determined to depend on either the dynamic behaviors of individual neurons or the whole SCN neuronal network. In this review, we first introduce possible network structures for the SCN, as revealed by time series analysis from real experimental data. It was found that the SCN network is heterogeneous and sparse, that is, the average shortest path length is very short, some nodes are hubs with large node degrees but most nodes have small node degrees, and the average node degree of the network is small. Secondly, the effects of the SCN network structure on the SCN function are reviewed based on mathematical models of the SCN network. It was found that robust rhythms with large amplitudes, a high synchronization between SCN neurons and a large entrainment ability exists mainly in small-world and scale-free type networks, but not other types. We conclude that the SCN most probably is an efficient small-world type or scale-free type network, which drives SCN function.

Modelling the functional roles of synaptic and extra-synaptic γ-aminobutyric acid receptor dynamics in circadian timekeeping

In the suprachiasmatic nucleus (SCN), γ-aminobutyric acid (GABA) is a primary neurotransmitter. GABA can signal through two types of GABAA receptor subunits, often referred to as synaptic GABAA (gamma subunit) and extra-synaptic GABAA (delta subunit). To test the functional roles of these distinct GABAA in regulating circadian rhythms, we developed a multicellular SCN model where we could separately compare the effects of manipulating GABA neurotransmitter or receptor dynamics. Our model predicted that blocking GABA signalling modestly increased synchrony among circadian cells, consistent with published SCN pharmacology. Conversely, the model predicted that lowering GABAA receptor density reduced firing rate, circadian cell fraction, amplitude and synchrony among individual neurons. When we tested these predictions, we found that the knockdown of delta GABAA reduced the amplitude and synchrony of clock gene expression among cells in SCN explants. The model further predicted that increasing gamma GABAA densities could enhance synchrony, as opposed to increasing delta GABAA densities. Overall, our model reveals how blocking GABAA receptors can modestly increase synchrony, while increasing the relative density of gamma over delta subunits can dramatically increase synchrony. We hypothesize that increased gamma GABAA density in the winter could underlie the tighter phase relationships among SCN cells.

Universal principles governing multiple random searchers on complex networks: The logarithmic growth pattern and the harmonic law

We propose a unified framework to evaluate and quantify the search time of multiple random searchers traversing independently and concurrently on complex networks. We find that the intriguing behaviors of multiple random searchers are governed by two basic principles-the logarithmic growth pattern and the harmonic law. Specifically, the logarithmic growth pattern characterizes how the search time increases with the number of targets, while the harmonic law explores how the search time of multiple random searchers varies relative to that needed by individual searchers. Numerical and theoretical results demonstrate these two universal principles established across a broad range of random search processes, including generic random walks, maximal entropy random walks, intermittent strategies, and persistent random walks. Our results reveal two fundamental principles governing the search time of multiple random searchers, which are expected to facilitate investigation of diverse dynamical processes like synchronization and spreading.

Globally fixed-time synchronization of coupled neutral-type neural network with mixed time-varying delays

This paper mainly studies the globally fixed-time synchronization of a class of coupled neutral-type neural networks with mixed time-varying delays via discontinuous feedback controllers. Compared with the traditional neutral-type neural network model, the model in this paper is more general. A class of general discontinuous feedback controllers are designed. With the help of the definition of fixed-time synchronization, the upper right-hand derivative and a defined simple Lyapunov function, some easily verifiable and extensible synchronization criteria are derived to guarantee the fixed-time synchronization between the drive and response systems. Finally, two numerical simulations are given to verify the correctness of the results.

Differences in intrinsic amplitudes of neuronal oscillators improve synchronization in the suprachiasmatic nucleus

In mammals, a main clock located in the suprachiasmatic nucleus (SCN) regulates the ∼24 h rhythms of behavioral and physiological activities exposed to a natural 24 light-dark cycle or even under constant darkness. The rhythms originate from self-sustained oscillations of the SCN neurons, which differ in both intrinsic periods and intrinsic amplitudes. The intrinsic periods and the intrinsic amplitudes were found to be bound to specific regions in the previous experiments. In particular, neurons of smaller amplitudes and larger periods are located in a ventrolateral part, and neurons of larger amplitudes and smaller periods are in a dorsomedial part. In the present study, we examined the effects of the differences in the intrinsic frequencies and the differences in the intrinsic amplitudes of neuronal oscillators on the synchronization, respectively. We found that the differences in the intrinsic frequencies weaken the synchronization, whereas the differences in the intrinsic amplitudes strengthen the synchronization. Our finding may shed light on the effects of the heterogenous properties of individual neurons on the collective behaviors of the SCN network and provide a way to enhance the synchronization.

Entrainment range of the suprachiasmatic nucleus affected by the difference in the neuronal amplitudes between the light-sensitive and light-insensitive regions

Mammals not only can be synchronized to the natural 24-h light-dark cycle, but also to a cycle with a non-24-h period. The range of the period of the external cycle, for which the animals can be entrained to, is called the entrainment range, which differs among species. The entrainment range as a characteristic of the animal is determined by the main circadian clock, i.e., the suprachiasmatic nucleus (SCN) in the brain. The SCN is composed of ∼10000 heterogeneous neurons, which can be divided into two subgroups, i.e., the ventrolateral subgroup (VL) directly receiving the light information from the retina and relaying the information to the dorsomedial subgroup (DM). Among the SCN neurons, the amplitudes are different; however, it is unclear that the amplitude is related to the location of the neurons in experiments. In the present study, we examined the effect of the difference in the neuronal amplitude between the VL and the DM on the entrainment range of the SCN, based on a mathematical model, i.e., the Poincaré model, which is used to describe the circadian clock. We find that the maximal entrainment range is obtained when the difference is equal to a critical point. If the difference of the amplitudes of the VL neurons to the amplitudes of the DM neurons is smaller than a critical point, with the increase of the difference, the entrainment range of the SCN increases, while if the difference is larger than the critical point, the entrainment range decreases with the increase of the difference. Our finding may give a potential explanation for the diversity of the entrainment range among species.