How the cortico-thalamic feedback affects the EEG power spectrum over frontal and occipital regions during propofol-induced sedation

Induced synchronization by endogenous noise modulation in finite-size random neural networks: A stochastic mean-field study

Event-related synchronization and desynchronization (ERS/ERD) are well-known features found experimentally in brain signals during cognitive tasks. Their understanding promises to have much better insights into neural information processes in cognition. Under the hypothesis that neural information affects the endogenous neural noise level in populations, we propose to employ a stochastic mean-field model to explain ERS/ERD in the γ-frequency range. The work extends previous mean-field studies by deriving novel effects from finite network size. Moreover, numerical simulations of ERS/ERD and their analytical explanation by the mean-field model suggest several endogenous noise modulation schemes, which may modulate the system's synchronization.

Propofol modulates neural dynamics of thalamo-cortical system associated with anesthetic levels in rats

The thalamocortical system plays an important role in consciousness. How anesthesia modulates the thalamocortical interactions is not completely known.  We simultaneously recorded local field potentials(LFPs) in thalamic reticular nucleus(TRN) and ventroposteromedial thalamic nucleus(VPM), and electrocorticographic(ECoG) activities in frontal and occipital cortices in freely moving rats (n = 11). We analyzed the changes in thalamic and cortical local spectral power and connectivities, which were measured with phase-amplitude coupling (PAC), coherence and multivariate Granger causality, at the states of baseline, intravenous infusion of propofol 20, 40, 80 mg/kg/h and after recovery of righting reflex. We found that propofol-induced burst-suppression results in a synchronous decrease of spectral power in thalamus and cortex (p < 0.001 for all frequency bands). The cross-frequency PAC increased by propofol, characterized by gradually stronger 'trough-max' pattern in TRN and stronger 'peak-max' pattern in cortex. The cross-region PAC increased in the phase of TRN modulating the amplitude of cortex. The functional connectivity (FC) between TRN and cortex for α/β bands also significantly increased (p < 0.040), with increased directional connectivity from TRN to cortex under propofol anesthesia. In contrast, the corticocortical FC significantly decreased (p < 0.047), with decreased directional connectivity from frontal cortex to occipital cortex. However, the thalamothalamic functional and directional connectivities remained largely unchanged by propofol anesthesia.  The spectral powers and connectivities are differentially modulated with the changes of propofol doses, suggesting the changes in neural dynamics in thalamocortical system could be used for distinguishing different vigilance levels caused by propofol.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-022-09912-0.

Delayed closed-loop neurostimulation for the treatment of pathological brain rhythms in mental disorders: a computational study

Mental disorders are among the top most demanding challenges in world-wide health. A large number of mental disorders exhibit pathological rhythms, which serve as the disorders characteristic biomarkers. These rhythms are the targets for neurostimulation techniques. Open-loop neurostimulation employs stimulation protocols, which are rather independent of the patients health and brain state in the moment of treatment. Most alternative closed-loop stimulation protocols consider real-time brain activity observations but appear as adaptive open-loop protocols, where e.g., pre-defined stimulation sets in if observations fulfil pre-defined criteria. The present theoretical work proposes a fully-adaptive closed-loop neurostimulation setup, that tunes the brain activities power spectral density (PSD) according to a user-defined PSD. The utilized brain model is non-parametric and estimated from the observations via magnitude fitting in a pre-stimulus setup phase. Moreover, the algorithm takes into account possible conduction delays in the feedback connection between observation and stimulation electrode. All involved features are illustrated on pathological α- and γ-rhythms known from psychosis. To this end, we simulate numerically a linear neural population brain model and a non-linear cortico-thalamic feedback loop model recently derived to explain brain activity in psychosis.

Amortized Bayesian inference on generative dynamical network models of epilepsy using deep neural density estimators

Whole-brain modeling of epilepsy combines personalized anatomical data with dynamical models of abnormal activities to generate spatio-temporal seizure patterns as observed in brain imaging data. Such a parametric simulator is equipped with a stochastic generative process, which itself provides the basis for inference and prediction of the local and global brain dynamics affected by disorders. However, the calculation of likelihood function at whole-brain scale is often intractable. Thus, likelihood-free algorithms are required to efficiently estimate the parameters pertaining to the hypothetical areas, ideally including the uncertainty. In this study, we introduce the simulation-based inference for the virtual epileptic patient model (SBI-VEP), enabling us to amortize the approximate posterior of the generative process from a low-dimensional representation of whole-brain epileptic patterns. The state-of-the-art deep learning algorithms for conditional density estimation are used to readily retrieve the statistical relationships between parameters and observations through a sequence of invertible transformations. We show that the SBI-VEP is able to efficiently estimate the posterior distribution of parameters linked to the extent of the epileptogenic and propagation zones from sparse intracranial electroencephalography recordings. The presented Bayesian methodology can deal with non-linear latent dynamics and parameter degeneracy, paving the way for fast and reliable inference on brain disorders from neuroimaging modalities.

Mathematical Model Insights into EEG Origin under Transcranial Direct Current Stimulation (tDCS) in the Context of Psychosis

Schizophrenia is a psychotic disease that develops progressively over years with a transition from prodromal to psychotic state associated with a disruption in brain activity. Transcranial Direct Current Stimulation (tDCS), known to alleviate pharmaco-resistant symptoms in patients suffering from schizophrenia, promises to prevent such a psychotic transition. To understand better how tDCS affects brain activity, we propose a neural cortico-thalamo-cortical (CTC) circuit model involving the Ascending Reticular Arousal System (ARAS) that permits to describe major impact features of tDCS, such as excitability for short-duration stimulation and electroencephalography (EEG) power modulation for long-duration stimulation. To this end, the mathematical model relates stimulus duration and Long-Term Plasticity (LTP) effect, in addition to describing the temporal LTP decay after stimulus offset. This new relation promises to optimize future stimulation protocols. Moreover, we reproduce successfully EEG-power modulation under tDCS in a ketamine-induced psychosis model and confirm the N-methyl-d-aspartate (NMDA) receptor hypofunction hypothesis in the etiopathophysiology of schizophrenia. The model description points to an important role of the ARAS and the δ-rhythm synchronicity in CTC circuit in early-stage psychosis.

Arousal Fluctuations Govern Oscillatory Transitions Between Dominant [Formula: see text] and [Formula: see text] Occipital Activity During Eyes Open/Closed Conditions

Arousal results in widespread activation of brain areas to increase their response in task and behavior relevant ways. Mediated by the Ascending Reticular Arousal System (ARAS), arousal-dependent inputs interact with neural circuitry to shape their dynamics. In the occipital cortex, such inputs may trigger shifts between dominant oscillations, where [Formula: see text] activity is replaced by [Formula: see text] activity, or vice versa. A salient example of this are spectral power alternations observed while eyes are opened and/or closed. These transitions closely follow fluctuations in arousal, suggesting a common origin. To better understand the mechanisms at play, we developed and analyzed a computational model composed of two modules: a thalamocortical feedback circuit coupled with a superficial cortical network. Upon activation by noise-like inputs originating from the ARAS, our model is able to demonstrate that noise-driven non-linear interactions mediate transitions in dominant peak frequency, resulting in the simultaneous suppression of [Formula: see text] limit cycle activity and the emergence of [Formula: see text] oscillations through coherence resonance. Reduction in input provoked the reverse effect - leading to anticorrelated transitions between [Formula: see text] and [Formula: see text] power. Taken together, these results shed a new light on how arousal shapes oscillatory brain activity.

Reconfiguration of αmplitude driven dominant coupling modes (DoCM) mediated by α-band in adolescents with schizophrenia spectrum disorders

Electroencephalography (EEG) based biomarkers have been shown to correlate with the presence of psychotic disorders. Increased delta and decreased alpha power in psychosis indicate an abnormal arousal state. We investigated brain activity across the basic EEG frequencies and also dynamic functional connectivity of both intra and cross-frequency coupling that could reveal a neurophysiological biomarker linked to an aberrant modulating role of alpha frequency in adolescents with schizophrenia spectrum disorders (SSDs). A dynamic functional connectivity graph (DFCG) has been estimated using the imaginary part of phase lag value (iPLV) and correlation of the envelope (corrEnv). We analyzed DFCG profiles of electroencephalographic resting state (eyes closed) recordings of healthy controls (HC) (n = 39) and SSDs subjects (n = 45) in basic frequency bands {δ,θ,α₁,α₂,β₁,β₂,γ}. In our analysis, we incorporated both intra and cross-frequency coupling modes. Adopting our recent Dominant Coupling Mode (DοCM) model leads to the construction of an integrated DFCG (iDFCG) that encapsulates the functional strength and the DοCM of every pair of brain areas. We revealed significantly higher ratios of delta/alpha₁,₂ power spectrum in SSDs subjects versus HC. The probability distribution (PD) of amplitude driven DoCM mediated by alpha frequency differentiated SSDs from HC with absolute accuracy (100%). The network Flexibility Index (FI) was significantly lower for subjects with SSDs compared to the HC group. Our analysis supports the central role of alpha frequency alterations in the neurophysiological mechanisms of SSDs. Currents findings open up new diagnostic pathways to clinical detection of SSDs and support the design of rational neurofeedback training.

Computational Models in Electroencephalography

Computational models lie at the intersection of basic neuroscience and healthcare applications because they allow researchers to test hypotheses in silico and predict the outcome of experiments and interactions that are very hard to test in reality. Yet, what is meant by “computational model” is understood in many different ways by researchers in different fields of neuroscience and psychology, hindering communication and collaboration. In this review, we point out the state of the art of computational modeling in Electroencephalography (EEG) and outline how these models can be used to integrate findings from electrophysiology, network-level models, and behavior. On the one hand, computational models serve to investigate the mechanisms that generate brain activity, for example measured with EEG, such as the transient emergence of oscillations at different frequency bands and/or with different spatial topographies. On the other hand, computational models serve to design experiments and test hypotheses in silico. The final purpose of computational models of EEG is to obtain a comprehensive understanding of the mechanisms that underlie the EEG signal. This is crucial for an accurate interpretation of EEG measurements that may ultimately serve in the development of novel clinical applications.

Serotonergic neurons in the dorsal raphe nucleus mediate the arousal-promoting effect of orexin during isoflurane anesthesia in male rats

Previous studies have demonstrated that the activation of orexinergic neurons facilitates the recovery of animals from general anesthesia. Moreover, serotonergic neurons that receive projections from orexin neurons have also been shown to participate in sleep-wakefulness regulation. In the present study, we aimed to explore whether orexinergic neurons facilitate emergence from isoflurane anesthesia in rats by activating serotonergic neurons. Orexin A (30 or 100 pmol), orexin B (30 or 100 pmol), and their respective antagonists SB-334867 and TCS-OX2-29 (5 or 20 μg) were microinjected into the dorsal raphe nucleus (DRN) of rats, and their effects on induction and emergence times were analyzed. Electroencephalogram (EEG) changes were also recorded and analyzed to illuminate the effect of orexin injection into the DRN on cortical excitability under isoflurane anesthesia. Activation of serotonergic neurons was detected via immunohistochemical analysis of c-Fos expression following orexin administration. Our results indicated that injection of neither orexins nor orexin antagonists into the rat DRN exerted an impact on induction time, whereas orexin-A injection (100 pmol) enhanced arousal when compared with the saline group. In contrast, administration of orexin receptor type 1 antagonist SB-334867 (20 μg) prolonged emergence time from isoflurane anesthesia. Microinjection of orexin-A induced an arousal pattern on EEG, and decreased the burst suppression ratio under isoflurane anesthesia. Isoflurane anesthesia inhibited the activity of serotonergic neurons, as shown by decrease in the number of c-Fos-immunoreactive serotonergic neurons when compared with the sham group. This inhibitory effect was partially reversed by administration of orexin-A. Taken together, our findings suggest that orexinergic signals facilitate emergence from isoflurane anesthesia, at least partially, by reversing the effects of isoflurane on serotonergic neurons of the DRN.

Neonatal exposure to propofol affects interneuron development in the piriform cortex and causes neurobehavioral deficits in adult mice

RATIONALE

Animal studies have shown that early postnatal propofol administration is involved in neurobehavioral alterations in adults. However, the underlying mechanism is not clear.

METHODS

We used c-Fos immunohistochemistry to identify activated neurons in brain regions of neonatal mice under propofol exposure and performed behavioral tests to observe the long-term consequences.

RESULTS

Exposure to propofol (30g or 60 mg/kg) on P7 produced significant c-Fos expression in the deep layers of the piriform cortex on P8. Double immunofluorescence of c-Fos with interneuron markers in the piriform cortex revealed that c-Fos was specifically induced in calbindin (CB)-positive interneurons. Repeated propofol exposure from P7 to P9 induced behavioral deficits in adult mice, such as olfactory function deficit in a buried food test, decreased sociability in a three-chambered choice task, and impaired recognitive ability of learning and memory in novel object recognition tests. However, locomotor activity in the open-field test was not generally affected. Propofol treatment also significantly decreased the number of CB-positive interneurons in the piriform cortex of mice on P21 and adulthood.

CONCLUSIONS

These results suggest that CB-positive interneurons in the piriform cortex are vulnerable to propofol exposure during the neonatal period, and these neurons are involved in the damage effects of propofol on behavior changes. These data provide a new target of propofol neurotoxicity and may elucidate the mechanism of neurobehavioral deficits in adulthood.

Optimal Model Parameter Estimation from EEG Power Spectrum Features Observed during General Anesthesia

Mathematical modeling is a powerful tool that enables researchers to describe the experimentally observed dynamics of complex systems. Starting with a robust model including model parameters, it is necessary to choose an appropriate set of model parameters to reproduce experimental data. However, estimating an optimal solution of the inverse problem, i.e., finding a set of model parameters that yields the best possible fit to the experimental data, is a very challenging problem. In the present work, we use different optimization algorithms based on a frequentist approach, as well as Monte Carlo Markov Chain methods based on Bayesian inference techniques to solve the considered inverse problems. We first probe two case studies with synthetic data and study models described by a stochastic non-delayed linear second-order differential equation and a stochastic linear delay differential equation. In a third case study, a thalamo-cortical neural mass model is fitted to the EEG spectral power measured during general anesthesia induced by anesthetics propofol and desflurane. We show that the proposed neural mass model fits very well to the observed EEG power spectra, particularly to the power spectral peaks within δ - (0 - 4 Hz) and α - (8 - 13 Hz) frequency ranges. Furthermore, for each case study, we perform a practical identifiability analysis by estimating the confidence regions of the parameter estimates and interpret the corresponding correlation and sensitivity matrices. Our results indicate that estimating the model parameters from analytically computed spectral power, we are able to accurately estimate the unknown parameters while avoiding the computational costs due to numerical integration of the model equations.

Suppression of underlying neuronal fluctuations mediates EEG slowing during general anaesthesia

The physiological mechanisms by which anaesthetic drugs modulate oscillatory brain activity remain poorly understood. Combining human data, mathematical and computational analysis of both spiking and mean-field models, we investigated the spectral dynamics of encephalographic (EEG) beta-alpha oscillations, observed in human patients undergoing general anaesthesia. The effect of anaesthetics can be modelled as a reduction of neural fluctuation intensity, and/or an increase in inhibitory synaptic gain in the thalamo-cortical circuit. Unlike previous work, which suggested the primary importance of gamma-amino-butryic-acid (GABA) augmentation in causing a shift to low EEG frequencies, our analysis demonstrates that a non-linear transition, triggered by a simple decrease in neural fluctuation intensity, is sufficient to explain the clinically-observed appearance - and subsequent slowing - of the beta-alpha narrowband EEG peak. In our model, increased synaptic inhibition alone, did not correlate with the clinically-observed encephalographic spectral changes, but did cause the anaesthetic-induced decrease in neuronal firing rate. Taken together, our results show that such a non-linear transition results in functional fragmentation of cortical and thalamic populations; highly correlated intra-population dynamics triggered by anaesthesia decouple and isolate neural populations. Our results are able to parsimoniously unify and replicate the observed anaesthetic effects on both the EEG spectra and inter-regional connectivity, and further highlight the importance of neural activity fluctuations in the genesis of altered brain states.

Electroencephalographic delta/alpha frequency activity differentiates psychotic disorders: a study of schizophrenia, bipolar disorder and methamphetamine-induced psychotic disorder

Electroencephalography (EEG) has been proposed as a neurophysiological biomarker to delineate psychotic disorders. It is known that increased delta and decreased alpha, which are apparent in psychosis, are indicative of inappropriate arousal state, which leads to reduced ability to attend to relevant information. On this premise, we investigated delta/alpha frequency activity, as this ratio of frequency activity may serve as an effective neurophysiological biomarker. The current study investigated differences in delta/alpha frequency activity, in schizophrenia (SCZ), bipolar I disorder with psychotic features and methamphetamine-induced psychosis. One hundred and nine participants, including individuals with SCZ (n = 28), bipolar I disorder with psychotic features (n = 28), methamphetamine-induced psychotic disorder (MPD) (n = 24) and healthy controls (CON, n = 29). Diagnosis was ascertained with the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders, 4th Edition disorders and current medication was recorded. EEG was undertaken in three testing conditions: resting eyes open, resting eyes closed and during completion of a simple cognitive task (visual continuous performance task). EEG delta/alpha frequency activity was investigated across these conditions. First, delta/alpha frequency activity during resting eyes closed was higher in SCZ and MPD globally, when compared to CON, then lower for bipolar disorder (BPD) than MPD for right hemisphere. Second, delta/alpha frequency activity during resting eyes open was higher in SCZ, BPD and MPD for all electrodes, except left frontal, when compared to CON. Third, delta/alpha frequency activity during the cognitive task was higher in BPD and MPD for all electrodes, except left frontal, when compared to CON. Assessment of EEG delta/alpha frequency activity supports the delineation of underlying neurophysiological mechanisms present in psychotic disorders, which are likely related to dysfunctional thalamo-cortical connectivity. Delta/alpha frequency activity may provide a useful neurophysiological biomarker to delineate psychotic disorders.

Noise-driven neuromorphic tuned amplifier

We study a simple stochastic model of neuronal excitatory and inhibitory interactions. The model is defined on a directed lattice and internodes couplings are modulated by a nonlinear function that mimics the process of synaptic activation. We prove that such a system behaves as a fully tunable amplifier: the endogenous component of noise, stemming from finite size effects, seeds a coherent (exponential) amplification across the chain generating giant oscillations with tunable frequencies, a process that the brain could exploit to enhance, and eventually encode, different signals. On a wider perspective, the characterized amplification process could provide a reliable pacemaking mechanism for biological systems. The device extracts energy from the finite size bath and operates as an out of equilibrium thermal machine, under stationary conditions.

Network Properties in Transitions of Consciousness during Propofol-induced Sedation

Reliable electroencephalography (EEG) signatures of transitions between consciousness and unconsciousness under anaesthesia have not yet been identified. Herein we examined network changes using graph theoretical analysis of high-density EEG during patient-titrated propofol-induced sedation. Responsiveness was used as a surrogate for consciousness. We divided the data into five states: baseline, transition into unresponsiveness, unresponsiveness, transition into responsiveness, and recovery. Power spectral analysis showed that delta power increased from responsiveness to unresponsiveness. In unresponsiveness, delta waves propagated from frontal to parietal regions as a traveling wave. Local increases in delta connectivity were evident in parietal but not frontal regions. Graph theory analysis showed that increased local efficiency could differentiate the levels of responsiveness. Interestingly, during transitions of responsive states, increased beta connectivity was noted relative to consciousness and unconsciousness, again with increased local efficiency. Abrupt network changes are evident in the transitions in responsiveness, with increased beta band power/connectivity marking transitions between responsive states, while the delta power/connectivity changes were consistent with the fading of consciousness using its surrogate responsiveness. These results provide novel insights into the neural correlates of these behavioural transitions and EEG signatures for monitoring the levels of consciousness under sedation.

Anesthetic action on the transmission delay between cortex and thalamus explains the beta-buzz observed under propofol anesthesia

In recent years, more and more surgeries under general anesthesia have been performed with the assistance of electroencephalogram (EEG) monitors. An increase in anesthetic concentration leads to characteristic changes in the power spectra of the EEG. Although tracking the anesthetic-induced changes in EEG rhythms can be employed to estimate the depth of anesthesia, their precise underlying mechanisms are still unknown. A prominent feature in the EEG of some patients is the emergence of a strong power peak in the β-frequency band, which moves to the α-frequency band while increasing the anesthetic concentration. This feature is called the beta-buzz. In the present study, we use a thalamo-cortical neural population feedback model to reproduce observed characteristic features in frontal EEG power obtained experimentally during propofol general anesthesia, such as this beta-buzz. First, we find that the spectral power peak in the α- and δ-frequency ranges depend on the decay rate constant of excitatory and inhibitory synapses, but the anesthetic action on synapses does not explain the beta-buzz. Moreover, considering the action of propofol on the transmission delay between cortex and thalamus, the model reveals that the beta-buzz may result from a prolongation of the transmission delay by increasing propofol concentration. A corresponding relationship between transmission delay and anesthetic blood concentration is derived. Finally, an analytical stability study demonstrates that increasing propofol concentration moves the systems resting state towards its stability threshold.

Lack of Responsiveness during the Onset and Offset of Sevoflurane Anesthesia Is Associated with Decreased Awake-Alpha Oscillation Power

Anesthetic drugs are typically administered to induce altered states of arousal that range from sedation to general anesthesia (GA). Systems neuroscience studies are currently being used to investigate the neural circuit mechanisms of anesthesia-induced altered arousal states. These studies suggest that by disrupting the oscillatory dynamics that are associated with arousal states, anesthesia-induced oscillations are a putative mechanism through which anesthetic drugs produce altered states of arousal. However, an empirical clinical observation is that even at relatively stable anesthetic doses, patients are sometimes intermittently responsive to verbal commands during states of light sedation. During these periods, prominent anesthesia-induced neural oscillations such as slow-delta (0.1–4 Hz) oscillations are notably absent. Neural correlates of intermittent responsiveness during light sedation have been insufficiently investigated. A principled understanding of the neural correlates of intermittent responsiveness may fundamentally advance our understanding of neural dynamics that are essential for maintaining arousal states, and how they are disrupted by anesthetics. Therefore, we performed a high-density (128 channels) electroencephalogram (EEG) study (n = 8) of sevoflurane-induced altered arousal in healthy volunteers. We administered temporally precise behavioral stimuli every 5 s to assess responsiveness. Here, we show that decreased eyes-closed, awake-alpha (8–12 Hz) oscillation power is associated with lack of responsiveness during sevoflurane effect-onset and -offset. We also show that anteriorization—the transition from occipitally dominant awake-alpha oscillations to frontally dominant anesthesia induced-alpha oscillations—is not a binary phenomenon. Rather, we suggest that periods, which were defined by lack of responsiveness, represent an intermediate brain state. We conclude that awake-alpha oscillation, previously thought to be an idling rhythm, is associated with responsiveness to behavioral stimuli.

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.