Power Law versus Exponential State Transition Dynamics: Application to Sleep-Wake Architecture

BACKGROUND: Despite the common experience that interrupted sleep has a negative impact on waking function, the features of human sleep-wake architecture that best distinguish sleep continuity versus fragmentation remain elusive. In this regard, there is growing interest in characterizing sleep architecture using models of the temporal dynamics of sleep-wake stage transitions. In humans and other mammals, the state transitions defining sleep and wake bout durations have been described with exponential and power law models, respectively. However, sleep-wake stage distributions are often complex, and distinguishing between exponential and power law processes is not always straightforward. Although mono-exponential distributions are distinct from power law distributions, multi-exponential distributions may in fact resemble power laws by appearing linear on a log-log plot. METHODOLOGY/PRINCIPAL FINDINGS: To characterize the parameters that may allow these distributions to mimic one another, we systematically fitted multi-exponential-generated distributions with a power law model, and power law-generated distributions with multi-exponential models. We used the Kolmogorov-Smirnov method to investigate goodness of fit for the "incorrect" model over a range of parameters. The "zone of mimicry" of parameters that increased the risk of mistakenly accepting power law fitting resembled empiric time constants obtained in human sleep and wake bout distributions. CONCLUSIONS/SIGNIFICANCE: Recognizing this uncertainty in model distinction impacts interpretation of transition dynamics (self-organizing versus probabilistic), and the generation of predictive models for clinical classification of normal and pathological sleep architecture

An alerting structure: Human orexin receptor 1

[No abstract available

Role of Homer Proteins in the Maintenance of Sleep-Wake States

Sleep is an evolutionarily conserved process that is linked to diurnal cycles and normal daytime wakefulness. Healthy sleep and wakefulness are integral to a healthy lifestyle; this occurs when an organism is able to maintain long bouts of both sleep and wake. Homer proteins, which function as adaptors for group 1 metabotropic glutamate receptors, have been implicated in genetic studies of sleep in both Drosophila and mouse. Drosophila express a single Homer gene product that is upregulated during sleep. By contrast, vertebrates express Homer as both constitutive and immediate early gene (H1a) forms, and H1a is up-regulated during wakefulness. Genetic deletion of Homer in Drosophila results in fragmented sleep and in failure to sustain long bouts of sleep, even under increased sleep drive. However, deletion of Homer1a in mouse results in failure to sustain long bouts of wakefulness. Further evidence for the role of Homer1a in the maintenance of wake comes from the CREB alpha delta mutant mouse, which displays a reduced wake phenotype similar to the Homer1a knockout and fails to up-regulate Homer1a upon sleep loss. Homer1a is a gene whose expression is induced by CREB. Sustained behaviors of the sleep/wake cycle are created by molecular pathways that are distinct from those for arousal or short bouts, and implicate an evolutionarily-conserved role for Homer in sustaining these behaviors

Obstructive Sleep Apnea Alters Sleep Stage Transition Dynamics

Enhanced characterization of sleep architecture, compared with routine polysomnographic metrics such as stage percentages and sleep efficiency, may improve the predictive phenotyping of fragmented sleep. One approach involves using stage transition analysis to characterize sleep continuity.We analyzed hypnograms from Sleep Heart Health Study (SHHS) participants using the following stage designations: wake after sleep onset (WASO), non-rapid eye movement (NREM) sleep, and REM sleep. We show that individual patient hypnograms contain insufficient number of bouts to adequately describe the transition kinetics, necessitating pooling of data. We compared a control group of individuals free of medications, obstructive sleep apnea (OSA), medical co-morbidities, or sleepiness (n = 374) with mild (n = 496) or severe OSA (n = 338). WASO, REM sleep, and NREM sleep bout durations exhibited multi-exponential temporal dynamics. The presence of OSA accelerated the "decay" rate of NREM and REM sleep bouts, resulting in instability manifesting as shorter bouts and increased number of stage transitions. For WASO bouts, previously attributed to a power law process, a multi-exponential decay described the data well. Simulations demonstrated that a multi-exponential process can mimic a power law distribution.OSA alters sleep architecture dynamics by decreasing the temporal stability of NREM and REM sleep bouts. Multi-exponential fitting is superior to routine mono-exponential fitting, and may thus provide improved predictive metrics of sleep continuity. However, because a single night of sleep contains insufficient transitions to characterize these dynamics, extended monitoring of sleep, probably at home, would be necessary for individualized clinical application

Cholinergic Modulation of Narcoleptic Attacks in Double Orexin Receptor Knockout Mice

To investigate how cholinergic systems regulate aspects of the sleep disorder narcolepsy, we video-monitored mice lacking both orexin (hypocretin) receptors (double knockout; DKO mice) while pharmacologically altering cholinergic transmission. Spontaneous behavioral arrests in DKO mice were highly similar to those reported in orexin-deficient mice and were never observed in wild-type (WT) mice. A survival analysis revealed that arrest lifetimes were exponentially distributed indicating that random, Markovian processes determine arrest lifetime. Low doses (0.01, 0.03 mg/kg, IP), but not a high dose (0.08 mg/kg, IP) of the cholinesterase inhibitor physostigmine increased the number of arrests but did not alter arrest lifetimes. The muscarinic antagonist atropine (0.5 mg/kg, IP) decreased the number of arrests, also without altering arrest lifetimes. To determine if muscarinic transmission in pontine areas linked to REM sleep control also influences behavioral arrests, we microinjected neostigmine (50 nl, 62.5 µM) or neostigmine + atropine (62.5 µM and 111 µM respectively) into the nucleus pontis oralis and caudalis. Neostigmine increased the number of arrests in DKO mice without altering arrest lifetimes but did not provoke arrests in WT mice. Co-injection of atropine abolished this effect. Collectively, our findings establish that behavioral arrests in DKO mice are similar to those in orexin deficient mice and that arrests have exponentially distributed lifetimes. We also show, for the first time in a rodent narcolepsy model, that cholinergic systems can regulate arrest dynamics. Since perturbations of muscarinic transmission altered arrest frequency but not lifetime, our findings suggest cholinergic systems influence arrest initiation without influencing circuits that determine arrest duration

The Neuronal Transition Probability (NTP) Model for the Dynamic Progression of Non-REM Sleep EEG: The Role of the Suprachiasmatic Nucleus

Little attention has gone into linking to its neuronal substrates the dynamic structure of non-rapid-eye-movement (NREM) sleep, defined as the pattern of time-course power in all frequency bands across an entire episode. Using the spectral power time-courses in the sleep electroencephalogram (EEG), we showed in the typical first episode, several moves towards-and-away from deep sleep, each having an identical pattern linking the major frequency bands beta, sigma and delta. The neuronal transition probability model (NTP) – in fitting the data well – successfully explained the pattern as resulting from stochastic transitions of the firing-rates of the thalamically-projecting brainstem-activating neurons, alternating between two steady dynamic-states (towards-and-away from deep sleep) each initiated by a so-far unidentified flip-flop. The aims here are to identify this flip-flop and to demonstrate that the model fits well all NREM episodes, not just the first. Using published data on suprachiasmatic nucleus (SCN) activity we show that the SCN has the information required to provide a threshold-triggered flip-flop for timing the towards-and-away alternations, information provided by sleep-relevant feedback to the SCN. NTP then determines the pattern of spectral power within each dynamic-state. NTP was fitted to individual NREM episodes 1–4, using data from 30 healthy subjects aged 20–30 years, and the quality of fit for each NREM measured. We show that the model fits well all NREM episodes and the best-fit probability-set is found to be effectively the same in fitting all subject data. The significant model-data agreement, the constant probability parameter and the proposed role of the SCN add considerable strength to the model. With it we link for the first time findings at cellular level and detailed time-course data at EEG level, to give a coherent picture of NREM dynamics over the entire night and over hierarchic brain levels all the way from the SCN to the EEG

A Thalamocortical Neural Mass Model of the EEG during NREM Sleep and Its Response to Auditory Stimulation

Few models exist that accurately reproduce the complex rhythms of the thalamocortical system that are apparent in measured scalp EEG and at the same time, are suitable for large-scale simulations of brain activity. Here, we present a neural mass model of the thalamocortical system during natural non-REM sleep, which is able to generate fast sleep spindles (12–15 Hz), slow oscillations (<1 Hz) and K-complexes, as well as their distinct temporal relations, and response to auditory stimuli. We show that with the inclusion of detailed calcium currents, the thalamic neural mass model is able to generate different firing modes, and validate the model with EEG-data from a recent sleep study in humans, where closed-loop auditory stimulation was applied. The model output relates directly to the EEG, which makes it a useful basis to develop new stimulation protocols