A probabilistic framework for a physiological representation of dynamically evolving sleep state

Unsupervised Estimation of Mouse Sleep Scores and Dynamics Using a Graphical Model of Electrophysiological Measurements

The proportion, number of bouts, and mean bout duration of different vigilance states (Wake, NREM, REM) are useful indices of dynamics in experimental sleep research. These metrics are estimated by first scoring state, sometimes using an algorithm, based on electrophysiological measurements such as the electroencephalogram (EEG) and electromyogram (EMG), and computing their values from the score sequence. Isolated errors in the scores can lead to large discrepancies in the estimated sleep metrics. But most algorithms score sleep by classifying the state from EEG/EMG features independently in each time epoch without considering the dynamics across epochs, which could provide contextual information. The objective here is to improve estimation of sleep metrics by fitting a probabilistic dynamical model to mouse EEG/EMG data and then predicting the metrics from the model parameters. Hidden Markov models (HMMs) with multivariate Gaussian observations and Markov state transitions were fitted to unlabeled 24-h EEG/EMG feature time series from 20 mice to model transitions between the latent vigilance states; a similar model with unbiased transition probabilities served as a reference. Sleep metrics predicted from the HMM parameters did not deviate significantly from manual estimates except for rapid eye movement sleep (REM) ([Formula: see text]; Wilcoxon signed-rank test). Changes in value from Light to Dark conditions correlated well with manually estimated differences (Spearman’s rho 0.43–0.84) except for REM. HMMs also scored vigilance state with over 90% accuracy. HMMs of EEG/EMG features can therefore characterize sleep dynamics from EEG/EMG measurements, a prerequisite for characterizing the effects of perturbation in sleep monitoring and control applications.

Closed-loop feedback control and bifurcation analysis of epileptiform activity via optogenetic stimulation in a mathematical model of human cortex

Optogenetics provides a method of neuron stimulation that has high spatial, temporal, and cell-type specificity. Here we present a model of optogenetic feedback control that targets the inhibitory population, which expresses light-sensitive channelrhodopsin-2 channels, in a mean-field model of undifferentiated cortex that is driven to seizures. The inhibitory population is illuminated with an intensity that is a function of electrode measurements obtained via the cortical model. We test the efficacy of this control method on seizurelike activity observed in two parameter spaces of the cortical model that most closely correspond to seizures observed in patients. We also compare the effect of closed-loop and open-loop control on seizurelike activity using a less-complicated ordinary differential equation model of the undifferentiated cortex in parameter space. Seizurelike activity is successfully suppressed in both parameter planes using optimal illumination intensities less likely to have adverse effects on cortical tissue.

A Pharmacokinetics-Neural Mass Model (PK-NMM) for the Simulation of EEG Activity during Propofol Anesthesia

Modeling the effects of anesthetic drugs on brain activity is very helpful in understanding anesthesia mechanisms. The aim of this study was to set up a combined model to relate actual drug levels to EEG dynamics and behavioral states during propofol-induced anesthesia. We proposed a new combined theoretical model based on a pharmacokinetics (PK) model and a neural mass model (NMM), which we termed PK-NMM--with the aim of simulating electroencephalogram (EEG) activity during propofol-induced general anesthesia. The PK model was used to derive propofol effect-site drug concentrations (C(eff)) based on the actual drug infusion regimen. The NMM model took C(eff) as the control parameter to produce simulated EEG-like (sEEG) data. For comparison, we used real prefrontal EEG (rEEG) data of nine volunteers undergoing propofol anesthesia from a previous experiment. To see how well the sEEG could describe the dynamic changes of neural activity during anesthesia, the rEEG data and the sEEG data were compared with respect to: power-frequency plots; nonlinear exponent (permutation entropy (PE)); and bispectral SynchFastSlow (SFS) parameters. We found that the PK-NMM model was able to reproduce anesthesia EEG-like signals based on the estimated drug concentration and patients' condition. The frequency spectrum indicated that the frequency power peak of the sEEG moved towards the low frequency band as anesthesia deepened. Different anesthetic states could be differentiated by the PE index. The correlation coefficient of PE was 0.80 ± 0.13 (mean ± standard deviation) between rEEG and sEEG for all subjects. Additionally, SFS could track the depth of anesthesia and the SFS of rEEG and sEEG were highly correlated with a correlation coefficient of 0.77 ± 0.13. The PK-NMM model could simulate EEG activity and might be a useful tool for understanding the action of propofol on brain activity.

Epoch-based Entropy for Early Screening of Alzheimer’s Disease

In this paper, we introduce a novel entropy measure, termed epoch-based entropy. This measure quantifies disorder of EEG signals both at the time level and spatial level, using local density estimation by a Hidden Markov Model on inter-channel stationary epochs. The investigation is led on a multi-centric EEG database recorded from patients at an early stage of Alzheimer’s disease (AD) and age-matched healthy subjects. We investigate the classification performances of this method, its robustness to noise, and its sensitivity to sampling frequency and to variations of hyperparameters. The measure is compared to two alternative complexity measures, Shannon’s entropy and correlation dimension. The classification accuracies for the discrimination of AD patients from healthy subjects were estimated using a linear classifier designed on a development dataset, and subsequently tested on an independent test set. Epoch-based entropy reached a classification accuracy of 83% on the test dataset (specificity = 83.3%, sensitivity = 82.3%), outperforming the two other complexity measures. Furthermore, it was shown to be more stable to hyperparameter variations, and less sensitive to noise and sampling frequency disturbances than the other two complexity measures.

Real-time automated EEG tracking of brain states using neural field theory

A real-time fitting system is developed and used to fit the predictions of an established physiologically-based neural field model to electroencephalographic spectra, yielding a trajectory in a physiological parameter space that parametrizes intracortical, intrathalamic, and corticothalamic feedbacks as the arousal state evolves continuously over time. This avoids traditional sleep/wake staging (e.g., using Rechtschaffen-Kales stages), which is fundamentally limited because it forces classification of continuous dynamics into a few discrete categories that are neither physiologically informative nor individualized. The classification is also subject to substantial interobserver disagreement because traditional staging relies in part on subjective evaluations. The fitting routine objectively and robustly tracks arousal parameters over the course of a full night of sleep, and runs in real-time on a desktop computer. The system developed here supersedes discrete staging systems by representing arousal states in terms of physiology, and provides an objective measure of arousal state which solves the problem of interobserver disagreement. Discrete stages from traditional schemes can be expressed in terms of model parameters for backward compatibility with prior studies.

Emergence from general anesthesia and the sleep-manifold

The electroencephalogram (EEG) during the re-establishment of consciousness after general anesthesia and surgery varies starkly between patients. Can the EEG during this emergence period provide a means of estimating the underlying biological processes underpinning the return of consciousness? Can we use a model to infer these biological processes from the EEG patterns? A frontal EEG was recorded from 84 patients. Ten patients were chosen for state-space analysis. Five showed archetypal emergences; which consisted of a progressive decrease in alpha power and increase peak alpha frequency before return of responsiveness. The five non-archetypal emergences showed almost no spectral EEG changes (even as the volatile general anesthetic decreased) and then an abrupt return of responsiveness. We used Bayesian methods to estimate the likelihood of an EEG pattern corresponding to the position of the patient on a 2-dimensional manifold in a state space of excitatory connection strength vs. change in intrinsic resting neuronal membrane conductivity. We could thus visualize the trajectory of each patient in the state-space during their emergence period. The patients who followed an archetypal emergence displayed a very consistent pattern; consisting of progressive increase in conductivity, and a temporary period of increased connection strength before return of responsiveness. The non-archetypal emergence trajectories remained fixed in a region of phase space characterized by a relatively high conductivity and low connection strength throughout emergence. This unexpected progressive increase in conductivity during archetypal emergence may be due to an abating of the surgical stimulus during this period. Periods of high connection strength could represent forays into dissociated consciousness, but the model suggests all patients reposition near the fold in the state space to take advantage of bi-stable cortical dynamics before transitioning to consciousness.