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

Amplified signal response by cluster synchronization competition in rings with short-distance couplings

Topological resonance has been revealed in degree-heterogeneous scale-free networks for weak signal amplification, but not in degree-homogeneous all-to-all networks [Acebrón et al., Phys. Rev. Lett. 99, 128701 (2007)0031-900710.1103/PhysRevLett.99.128701]. Here, we show that when the coupling distance of the all-to-all networks is reduced from global to local, i.e., converting all-to-all networks into rings, we can observe a resonant response to a weak signal similar to scale-free networks. We find that such a resonance effect induced by ring topology is robust across a wide range of ring sizes and signal frequencies. We further show that at intermediate coupling strength, oscillators in the rings can form separate synchronous clusters that compete with each other, resulting in large amplitude oscillations of boundary nodes between clusters and thus giving rise to the resonant signal amplification. Finally, we propose a structure of a three-node feed-forward motif simplified from the observed cluster synchronization competition to analyze the mechanism underlying the resonance behavior, which corresponds well with the numerical results.

Diversity-induced decoherence

We analyze the effect of small-amplitude noise and heterogeneity in a network of coupled excitable oscillators with strong timescale separation. Using mean-field analysis, we uncover the mechanism of a nontrivial effect-diversity-induced decoherence (DIDC)-in which heterogeneity modulates the mechanism of self-induced stochastic resonance to inhibit the coherence of oscillations. We argue that DIDC may offer one possible mechanism via which, in excitable neural systems, generic heterogeneity and background noise can synergistically prevent unwanted resonances that may be related to hyperkinetic movement disorders.

Stochastic Spatial Heterogeneity in Activities of H+-ATP-Ases in Electrically Connected Plant Cells Decreases Threshold for Cooling-Induced Electrical Responses

H⁺-ATP-ases, which support proton efflux through the plasma membrane, are key molecular transporters for electrogenesis in cells of higher plants. Initial activities of the transporters can influence the thresholds of generation of electrical responses induced by stressors and modify other parameters of these responses. Previously, it was theoretically shown that the stochastic heterogeneity of individual cell thresholds for electrical responses in a system of electrically connected neuronal cells can decrease the total threshold of the system (“diversity-induced resonance”, DIR). In the current work, we tested a hypothesis about decreasing the thresholds of generation of cooling-induced electrical responses in a system of electrically connected plant cells with increasing stochastic spatial heterogeny in the initial activities of H⁺-ATP-ases in these cells. A two-dimensional model of the system of electrically connected excitable cells (simple imitation of plant leaf), which was based on a model previously developed in our works, was used for the present investigation. Simulation showed that increasing dispersion in the distribution of initial activities of H⁺-ATP-ases between cells decreased the thresholds of generation of cooling-induced electrical responses. In addition, the increasing weakly influenced the amplitudes of electrical responses. Additional analysis showed two different mechanisms of the revealed effect. The increasing spatial heterogeneity in activities of H⁺-ATP-ases induced a weak positive shift of the membrane potential at rest. The shift decreased the threshold of electrical response generation. However, the decreased threshold induced by increasing the H⁺-ATP-ase activity heterogeneity was also observed after the elimination of the positive shift. The result showed that the “DIR-like” mechanism also participated in the revealed effect. Finally, we showed that the standard deviation of the membrane potentials before the induction of action potentials could be used for the estimation of thresholds of cooling-induced plant electrical responses. Thus, spatial heterogeneity in the initial activities of H⁺-ATP-ases can be a new regulatory mechanism influencing the generation of electrical responses in plants under actions of stressors.

Hubs, diversity, and synchronization in FitzHugh-Nagumo oscillator networks: Resonance effects and biophysical implications

Using the FitzHugh-Nagumo equations to represent the oscillatory electrical behavior of β-cells, we develop a coupled oscillator network model with cubic lattice topology, showing that the emergence of pacemakers or hubs in the system can be viewed as a natural consequence of oscillator population diversity. The optimal hub to nonhub ratio is determined by the position of the diversity-induced resonance maximum for a given set of FitzHugh-Nagumo equation parameters and is predicted by the model to be in a range that is fully consistent with experimental observations. The model also suggests that hubs in a β-cell network should have the ability to "switch on" and "off" their pacemaker function. As a consequence, their relative amount in the population can vary in order to ensure an optimal oscillatory performance of the network in response to environmental changes, such as variations of an external stimulus.

Optimal Self-Induced Stochastic Resonance in Multiplex Neural Networks: Electrical vs. Chemical Synapses

Electrical and chemical synapses shape the dynamics of neural networks, and their functional roles in information processing have been a longstanding question in neurobiology. In this paper, we investigate the role of synapses on the optimization of the phenomenon of self-induced stochastic resonance in a delayed multiplex neural network by using analytical and numerical methods. We consider a two-layer multiplex network in which, at the intra-layer level, neurons are coupled either by electrical synapses or by inhibitory chemical synapses. For each isolated layer, computations indicate that weaker electrical and chemical synaptic couplings are better optimizers of self-induced stochastic resonance. In addition, regardless of the synaptic strengths, shorter electrical synaptic delays are found to be better optimizers of the phenomenon than shorter chemical synaptic delays, while longer chemical synaptic delays are better optimizers than longer electrical synaptic delays; in both cases, the poorer optimizers are, in fact, worst. It is found that electrical, inhibitory, or excitatory chemical multiplexing of the two layers having only electrical synapses at the intra-layer levels can each optimize the phenomenon. Additionally, only excitatory chemical multiplexing of the two layers having only inhibitory chemical synapses at the intra-layer levels can optimize the phenomenon. These results may guide experiments aimed at establishing or confirming to the mechanism of self-induced stochastic resonance in networks of artificial neural circuits as well as in real biological neural networks.

Sleep Modelling across Physiological Levels

Sleep and circadian rhythms are regulated across multiple functional, spatial and temporal levels: from genes to networks of coupled neurons and glial cells, to large scale brain dynamics and behaviour. The dynamics at each of these levels are complex and the interaction between the levels is even more so, so research have mostly focused on interactions within the levels to understand the underlying mechanisms—the so-called reductionist approach. Mathematical models were developed to test theories of sleep regulation and guide new experiments at each of these levels and have become an integral part of the field. The advantage of modelling, however, is that it allows us to simulate and test the dynamics of complex biological systems and thus provides a tool to investigate the connections between the different levels and study the system as a whole. In this paper I review key models of sleep developed at different physiological levels and discuss the potential for an integrated systems biology approach for sleep regulation across these levels. I also highlight the necessity of building mechanistic connections between models of sleep and circadian rhythms across these levels.

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.

System-size expansion of the moments of a master equation

We study an expansion method of the general multidimensional master equation, based on a system-size expansion of the time evolution equations of the moments. The method turns out to be more accurate than the traditional van Kampen expansion for the first and second moments, with an error that scales with system-size similar to an alternative expansion, also applied to the equations of the moments, called Gaussian approximation, with the advantage that it has less systematic errors. Besides, we analyze a procedure to find the solution of the expansion method and we show different cases where it greatly simplifies. This includes the analytical solution of the average value and fluctuations in the number of infected nodes of the SIS epidemic model in complex networks, under the degree-based approximation.

Spiking patterns and synchronization of thalamic neurons along the sleep-wake cycle

Spiking patterns and synchronization dynamics of thalamic neurons along the sleep-wake cycle are studied in a minimal model of four coupled conductance-based neurons. The model simulates two thalamic neurons coupled via a gap junction and driven by a synaptic input from a two-neuron model of sleep regulation by the hypothalamus. In accord with experimental data, the model shows that during sleep, when hypothalamic wake-active neurons are silent, the thalamic neurons discharge bursts of spikes. During wake, the excitatory synaptic input from the hypothalamus drives the coupled thalamic neurons to a state of tonic firing (single spikes). In the deterministic case, the thalamic neurons synchronize in-phase in the bursting regime but demonstrate multi-stability of out-of-phase, in-phase, and asynchronous states in the tonic firing. However, along the sleep-wake cycle, once the neurons synchronize in-phase during sleep (bursting), they stay synchronized in wake (tonic firing). It is thus found that noise is needed to reproduce the experimentally observed transitions between synchronized bursting during sleep and asynchronous tonic firing during wake. Overall, synchronization of bursting is found to be more robust to noise than synchronization of tonic firing, where a small disturbance is sufficient to desynchronize the thalamic neurons. The model predicts that the transitions between sleep and wake happen via chaos because a single thalamic neuron exhibits chaos between regular bursting and tonic activity. The results of this study suggest that the sleep- and wake-related dynamics in the thalamus may be generated at a level of gap junction-coupled clusters of thalamic neurons driven from the hypothalamus which would then propagate throughout the thalamus and cortex via axonal long-range connections.

Activity patterns with silent states in a heterogeneous network of gap-junction coupled Huber-Braun model neurons

A mathematical model of a network of nearest neighbor gap-junction coupled neurons has been used to examine the impact of neuronal heterogeneity on the networks' activity during increasing coupling strength. Heterogeneity has been introduced by Huber-Braun model neurons with randomization of the temperature as a scaling factor. This leads to neurons of an enormous diversity of impulse pattern, including burst discharges, chaotic activity, and two different types of tonic firing-all of them experimentally observed in the peripheral as well as central nervous system. When the network is composed of all these types of neurons, randomly selected, a particular phenomenon can be observed. At a certain coupling strength, the network goes into a completely silent state. Examination of voltage traces and inter-spike intervals of individual neurons suggests that all neurons, irrespective of their original pattern, go through a well-known bifurcation scenario, resembling those of single neurons especially on external current injection. All the originally spontaneously firing neurons can achieve constant membrane potentials at which all intrinsic and gap-junction currents are balanced. With limited diversity, i.e., taking out neurons of specific patterns from the lower and upper temperature range, spontaneous firing can be reinstalled with further increasing coupling strength, especially when the tonic firing regimes are missing. Reinstalled firing develops from slowly increasing subthreshold oscillations leading to tonic firing activity with already fairly well synchronized action potentials, while the subthreshold potentials can still be significantly different. Full in phase synchronization is not achieved. Additional studies are needed elucidating the underlying mechanisms and the conditions under which such particular transitions can appear.

Neural circuit interactions between the dorsal raphe nucleus and the lateral hypothalamus: an experimental and computational study

Orexinergic/hypocretinergic (Ox) neurotransmission plays an important role in regulating sleep, as well as in anxiety and depression, for which the serotonergic (5-HT) system is also involved in. However, little is known regarding the direct and indirect interactions between 5-HT in the dorsal raphe nucleus (DRN) and Ox neurons in the lateral hypothalamus (LHA). In this study, we report the additional presence of 5-HT1BR, 5-HT2AR, 5-HT2CR and fast ligand-gated 5-HT3AR subtypes on the Ox neurons of transgenic Ox-enhanced green fluorescent protein (Ox-EGFP) and wild type C57Bl/6 mice using single and double immunofluorescence (IF) staining, respectively, and quantify the colocalization for each 5-HT receptor subtype. We further reveal the presence of 5-HT3AR and 5-HT1AR on GABAergic neurons in LHA. We also identify NMDAR1, OX1R and OX2R on Ox neurons, but none on adjacent GABAergic neurons. This suggests a one-way relationship between LHA's GABAergic and Ox neurons, wherein GABAergic neurons exerts an inhibitory effect on Ox neurons under partial DRN's 5-HT control. We also show that Ox axonal projections receive glutamatergic (PSD-95 immunopositive) and GABAergic (Gephyrin immunopositive) inputs in the DRN. We consider these and other available findings into our computational model to explore possible effects of neural circuit connection types and timescales on the DRN-LHA system's dynamics. We find that if the connections from 5-HT to LHA's GABAergic neurons are weakly excitatory or inhibitory, the network exhibits slow oscillations; not observed when the connection is strongly excitatory. Furthermore, if Ox directly excites 5-HT neurons at a fast timescale, phasic Ox activation can lead to an increase in 5-HT activity; no significant effect with slower timescale. Overall, our experimental and computational approaches provide insights towards a more complete understanding of the complex relationship between 5-HT in the DRN and Ox in the LHA.