Monitoring surgical nociception using multisensor physiological models.

Monitoring nociception, the flow of information associated with harmful stimuli through the nervous system even during unconsciousness, is critical for proper anesthesia care during surgery. Currently, this is done by tracking heart rate and blood pressure by eye. Monitoring objectively a patient's nociceptive state remains a challenge, causing drugs to often be over- or underdosed intraoperatively. Inefficient management of surgical nociception may lead to more complex postoperative pain management and side effects such as postoperative cognitive dysfunction, particularly in elderly patients. We collected a comprehensive and multisensor prospective observational dataset focused on surgical nociception (101 surgeries, 18,582 min, and 49,878 nociceptive stimuli), including annotations of all nociceptive stimuli occurring during surgery and medications administered. Using this dataset, we developed indices of autonomic nervous system activity based on physiologically and statistically rigorous point process representations of cardiac action potentials and sweat gland activity. Next, we constructed highly interpretable supervised and unsupervised models with appropriate inductive biases that quantify surgical nociception throughout surgery. Our models track nociceptive stimuli more accurately than existing nociception monitors. We also demonstrate that the characterizing signature of nociception learned by our models resembles the known physiology of the response to pain. Our work represents an important step toward objective multisensor physiology-based markers of surgical nociception. These markers are derived from an in-depth characterization of nociception as measured during surgery itself rather than using other experimental models as surrogates for surgical nociception.

Electrographic seizures during low-current thalamic deep brain stimulation in mice.

BACKGROUND: Deep brain stimulation of the central thalamus (CT-DBS) has potential for modulating states of consciousness, but it can also trigger electrographic seizures, including poly-spike-wave trains (PSWT). OBJECTIVES: To report the probability of inducing PSWTs during CT-DBS in awake, freely-moving mice. METHODS: Mice were implanted with electrodes to deliver unilateral and bilateral CT-DBS at different frequencies while recording electroencephalogram (EEG). We titrated stimulation current by gradually increasing it at each frequency until a PSWT appeared. Subsequent stimulations to test arousal modulation were performed at the current one step below the current that caused a PSWT during titration. RESULTS: In 2.21% of the test stimulations (10 out of 12 mice), CT-DBS caused PSWTs at currents lower than the titrated current, including currents as low as 20 µA. CONCLUSION: Our study found a small but significant probability of inducing PSWTs even after titration and at relatively low currents. EEG should be closely monitored for electrographic seizures when performing CT-DBS in both research and clinical settings.

Propofol anesthesia destabilizes neural dynamics across cortex.

Every day, hundreds of thousands of people undergo general anesthesia. One hypothesis is that anesthesia disrupts dynamic stability-the ability of the brain to balance excitability with the need to be stable and controllable. To test this hypothesis, we developed a method for quantifying changes in population-level dynamic stability in complex systems: delayed linear analysis for stability estimation (DeLASE). Propofol was used to transition animals between the awake state and anesthetized unconsciousness. DeLASE was applied to macaque cortex local field potentials (LFPs). We found that neural dynamics were more unstable in unconsciousness compared with the awake state. Cortical trajectories mirrored predictions from destabilized linear systems. We mimicked the effect of propofol in simulated neural networks by increasing inhibitory tone. This in turn destabilized the networks, as observed in the neural data. Our results suggest that anesthesia disrupts dynamical stability that is required for consciousness.

Dynamic Estimation of Cardiovascular State From Arterial Blood Pressure Recordings.

Real-time estimation of patient cardiovascular states, including cardiac output and systemic vascular resistance, is necessary for personalized hemodynamic monitoring and management. Highly invasive measurements enable reliable estimation of these states but increase patient risk. Prior methods using minimally invasive measurements reduce patient risk but have produced unreliable estimates limited due to trade-offs in accuracy and time resolution. Our objective was to develop an approach to estimate cardiac output and systemic vascular resistance with both a high time resolution and high accuracy from minimally invasive measurements. Using the two-element Windkessel model, we formulated a state-space method to estimate a dynamic time constant - the product of systemic vascular resistance and compliance - from arterial blood pressure measurements. From this time constant, we derived proportional estimates of systemic vascular resistance and cardiac output. We then validated our method with a swine cardiovascular dataset. Our estimates produced using arterial blood pressure measurements not only closely align with those using highly invasive measurements, but also closely align when derived from three separate locations on the arterial tree. Moreover, our estimates predictably change in response to standard cardiovascular drugs. Overall, our approach produces reliable, real-time estimates of cardiovascular states crucial for monitoring and control of the cardiovascular system.

Spike-wave discharges during low-current thalamic deep brain stimulation in mice.

BACKGROUND: Deep brain stimulation of central thalamus (CT-DBS) has potential for modulating states of consciousness, but it can also trigger spike-wave discharges (SWDs). OBJECTIVES: To report the probability of inducing SWDs during CT-DBS in awake mice. METHODS: Mice were implanted with electrodes to deliver unilateral and bilateral CT-DBS at different frequencies while recording EEG. We titrated stimulation current by gradually increasing it at each frequency until an SWD appeared. Subsequent stimulations to test arousal modulation were performed at the current one step below the current that caused an SWD during titration. RESULTS: In 2.21% of the test stimulations (10 out of 12 mice), CT-DBS caused SWDs at currents lower than the titrated current, at currents as low as 20 uA. CONCLUSION: Our study found a small but significant probability of inducing SWDs even after titration and at relatively low currents. EEG should be closely monitored for SWDs when performing CT-DBS in both research and clinical settings.

Ketamine can produce oscillatory dynamics by engaging mechanisms dependent on the kinetics of NMDA receptors.

Ketamine is an N-methyl-D-aspartate (NMDA)-receptor antagonist that produces sedation, analgesia, and dissociation at low doses and profound unconsciousness with antinociception at high doses. At high and low doses, ketamine can generate gamma oscillations (>25 Hz) in the electroencephalogram (EEG). The gamma oscillations are interrupted by slow-delta oscillations (0.1 to 4 Hz) at high doses. Ketamine's primary molecular targets and its oscillatory dynamics have been characterized. However, how the actions of ketamine at the subcellular level give rise to the oscillatory dynamics observed at the network level remains unknown. By developing a biophysical model of cortical circuits, we demonstrate how NMDA-receptor antagonism by ketamine can produce the oscillatory dynamics observed in human EEG recordings and nonhuman primate local field potential recordings. We have identified how impaired NMDA-receptor kinetics can cause disinhibition in neuronal circuits and how a disinhibited interaction between NMDA-receptor-mediated excitation and GABA-receptor-mediated inhibition can produce gamma oscillations at high and low doses, and slow-delta oscillations at high doses. Our work uncovers general mechanisms for generating oscillatory brain dynamics that differs from ones previously reported and provides important insights into ketamine's mechanisms of action as an anesthetic and as a therapy for treatment-resistant depression.

Multimodal MRI reveals brainstem connections that sustain wakefulness in human consciousness.

Consciousness is composed of arousal (i.e., wakefulness) and awareness. Substantial progress has been made in mapping the cortical networks that underlie awareness in the human brain, but knowledge about the subcortical networks that sustain arousal in humans is incomplete. Here, we aimed to map the connectivity of a proposed subcortical arousal network that sustains wakefulness in the human brain, analogous to the cortical default mode network (DMN) that has been shown to contribute to awareness. We integrated data from ex vivo diffusion magnetic resonance imaging (MRI) of three human brains, obtained at autopsy from neurologically normal individuals, with immunohistochemical staining of subcortical brain sections. We identified nodes of the proposed default ascending arousal network (dAAN) in the brainstem, hypothalamus, thalamus, and basal forebrain. Deterministic and probabilistic tractography analyses of the ex vivo diffusion MRI data revealed projection, association, and commissural pathways linking dAAN nodes with one another and with DMN nodes. Complementary analyses of in vivo 7-tesla resting-state functional MRI data from the Human Connectome Project identified the dopaminergic ventral tegmental area in the midbrain as a widely connected hub node at the nexus of the subcortical arousal and cortical awareness networks. Our network-based autopsy methods and connectivity data provide a putative neuroanatomic architecture for the integration of arousal and awareness in human consciousness.

Daytime dexmedetomidine sedation with closed-loop acoustic stimulation alters slow wave sleep homeostasis in healthy adults.

Background: The alpha-2 adrenergic agonist dexmedetomidine induces EEG patterns resembling those of non-rapid eye movement (NREM) sleep. Fulfilment of slow wave sleep (SWS) homeostatic needs would address the assumption that dexmedetomidine induces functional biomimetic sleep states. Methods: In-home sleep EEG recordings were obtained from 13 healthy participants before and after dexmedetomidine sedation. Dexmedetomidine target-controlled infusions and closed-loop acoustic stimulation were implemented to induce and enhance EEG slow waves, respectively. EEG recordings during sedation and sleep were staged using modified American Academy of Sleep Medicine criteria. Slow wave activity (EEG power from 0.5 to 4 Hz) was computed for NREM stage 2 (N2) and NREM stage 3 (N3/SWS) epochs, with the aggregate partitioned into quintiles by time. The first slow wave activity quintile served as a surrogate for slow wave pressure, and the difference between the first and fifth quintiles as a measure of slow wave pressure dissipation. Results: Compared with pre-sedation sleep, post-sedation sleep showed reduced N3 duration (mean difference of -17.1 min, 95% confidence interval -30.0 to -8.2, P=0.015). Dissipation of slow wave pressure was reduced (P=0.02). Changes in combined durations of N2 and N3 between pre- and post-sedation sleep correlated with total dexmedetomidine dose, (r=-0.61, P=0.03). Conclusions: Daytime dexmedetomidine sedation and closed-loop acoustic stimulation targeting EEG slow waves reduced N3/SWS duration and measures of slow wave pressure dissipation on the post-sedation night in healthy young adults. Thus, the paired intervention induces sleep-like states that fulfil certain homeostatic NREM sleep needs in healthy young adults. Clinical trial registration: ClinicalTrials.gov NCT04206059.

Convergent effects of different anesthetics are due to changes in phase alignment of cortical oscillations.

Many different anesthetics cause loss of responsiveness despite having diverse underlying molecular and circuit actions. To explore the convergent effects of these drugs, we examined how ketamine, an N-methyl-D-aspartate (NMDA) receptor antagonist, and dexmedetomidine, an α2 adrenergic receptor agonist, affected neural oscillations in the prefrontal cortex of nonhuman primates. Previous work has shown that anesthesia increases phase locking of low-frequency local field potential activity across cortex. We observed similar increases with anesthetic doses of ketamine and dexmedetomidine in the ventrolateral and dorsolateral prefrontal cortex, within and across hemispheres. However, the nature of the phase locking varied between regions. We found that oscillatory activity in different prefrontal subregions within each hemisphere became more anti-phase with both drugs. Local analyses within a region suggested that this finding could be explained by broad cortical distance-based effects, such as a large traveling wave. By contrast, homologous areas across hemispheres increased their phase alignment. Our results suggest that the drugs induce strong patterns of cortical phase alignment that are markedly different from those in the awake state, and that these patterns may be a common feature driving loss of responsiveness from different anesthetic drugs.

Ketamine can produce oscillatory dynamics by engaging mechanisms dependent on the kinetics of NMDA receptors.

Ketamine is an NMDA-receptor antagonist that produces sedation, analgesia and dissociation at low doses and profound unconsciousness with antinociception at high doses. At high and low doses, ketamine can generate gamma oscillations (>25 Hz) in the electroencephalogram (EEG). The gamma oscillations are interrupted by slow-delta oscillations (0.1-4 Hz) at high doses. Ketamine's primary molecular targets and its oscillatory dynamics have been characterized. However, how the actions of ketamine at the subcellular level give rise to the oscillatory dynamics observed at the network level remains unknown. By developing a biophysical model of cortical circuits, we demonstrate how NMDA-receptor antagonism by ketamine can produce the oscillatory dynamics observed in human EEG recordings and non-human primate local field potential recordings. We have discovered how impaired NMDA-receptor kinetics can cause disinhibition in neuronal circuits and how a disinhibited interaction between NMDA-receptor-mediated excitation and GABA-receptor-mediated inhibition can produce gamma oscillations at high and low doses, and slow-delta oscillations at high doses. Our work uncovers general mechanisms for generating oscillatory brain dynamics that differs from ones previously reported, and provides important insights into ketamine's mechanisms of action as an anesthetic and as a therapy for treatment-resistant depression.

Propofol-mediated Unconsciousness Disrupts Progression of Sensory Signals through the Cortical Hierarchy.

A critical component of anesthesia is the loss of sensory perception. Propofol is the most widely used drug for general anesthesia, but the neural mechanisms of how and when it disrupts sensory processing are not fully understood. We analyzed local field potential and spiking recorded from Utah arrays in auditory cortex, associative cortex, and cognitive cortex of nonhuman primates before and during propofol-mediated unconsciousness. Sensory stimuli elicited robust and decodable stimulus responses and triggered periods of stimulus-related synchronization between brain areas in the local field potential of Awake animals. By contrast, propofol-mediated unconsciousness eliminated stimulus-related synchrony and drastically weakened stimulus responses and information in all brain areas except for auditory cortex, where responses and information persisted. However, we found stimuli occurring during spiking Up states triggered weaker spiking responses than in Awake animals in auditory cortex, and little or no spiking responses in higher order areas. These results suggest that propofol's effect on sensory processing is not just because of asynchronous Down states. Rather, both Down states and Up states reflect disrupted dynamics.

Closed-loop control of anesthetic state in nonhuman primates.

Research in human volunteers and surgical patients has shown that unconsciousness under general anesthesia can be reliably tracked using real-time electroencephalogram processing. Hence, a closed-loop anesthesia delivery (CLAD) system that maintains precisely specified levels of unconsciousness is feasible and would greatly aid intraoperative patient management. The US Federal Drug Administration has approved no CLAD system for human use due partly to a lack of testing in appropriate animal models. To address this key roadblock, we implement a nonhuman primate (NHP) CLAD system that controls the level of unconsciousness using the anesthetic propofol. The key system components are a local field potential (LFP) recording system; propofol pharmacokinetics and pharmacodynamic models; the control variable (LFP power between 20 and 30 Hz), a programmable infusion system and a linear quadratic integral controller. Our CLAD system accurately controlled the level of unconsciousness along two different 125-min dynamic target trajectories for 18 h and 45 min in nine experiments in two NHPs. System performance measures were comparable or superior to those in previous CLAD reports. We demonstrate that an NHP CLAD system can reliably and accurately control in real-time unconsciousness maintained by anesthesia. Our findings establish critical steps for CLAD systems' design and testing prior to human testing.

Multifunctional fibers enable modulation of cortical and deep brain activity during cognitive behavior in macaques.

Recording and modulating neural activity in vivo enables investigations of the neurophysiology underlying behavior and disease. However, there is a dearth of translational tools for simultaneous recording and localized receptor-specific modulation. We address this limitation by translating multifunctional fiber neurotechnology previously only available for rodent studies to enable cortical and subcortical neural recording and modulation in macaques. We record single-neuron and broader oscillatory activity during intracranial GABA infusions in the premotor cortex and putamen. By applying state-space models to characterize changes in electrophysiology, we uncover that neural activity evoked by a working memory task is reshaped by even a modest local inhibition. The recordings provide detailed insight into the electrophysiological effect of neurotransmitter receptor modulation in both cortical and subcortical structures in an awake macaque. Our results demonstrate a first-time application of multifunctional fibers for causal studies of neuronal activity in behaving nonhuman primates and pave the way for clinical translation of fiber-based neurotechnology.

Propofol-mediated loss of consciousness disrupts predictive routing and local field phase modulation of neural activity.

Predictive coding is a fundamental function of the cortex. The predictive routing model proposes a neurophysiological implementation for predictive coding. Predictions are fed back from deep-layer cortex via alpha/beta (8-30Hz) oscillations. They inhibit the gamma (40-100Hz) and spiking that feed sensory inputs forward. Unpredicted inputs arrive in circuits unprepared by alpha/beta, resulting in enhanced gamma and spiking. To test the predictive routing model and its role in consciousness, we collected data from intracranial recordings of macaque monkeys during passive presentation of auditory oddballs (e.g., AAAAB) before and after propofol-mediated loss of consciousness (LOC). In line with the predictive routing model, alpha/beta oscillations in the awake state served to inhibit the processing of predictable stimuli. Propofol-mediated LOC eliminated alpha/beta modulation by a predictable stimulus in sensory cortex and alpha/beta coherence between sensory and frontal areas. As a result, oddball stimuli evoked enhanced gamma power, late (> 200 ms from stimulus onset) period spiking, and superficial layer sinks in sensory cortex. Therefore, auditory cortex was in a disinhibited state during propofol-mediated LOC. However, despite these enhanced feedforward responses in auditory cortex, there was a loss of differential spiking to oddballs in higher order cortex. This may be a consequence of a loss of within-area and inter-area spike-field coupling in the alpha/beta and gamma frequency bands. These results provide strong constraints for current theories of consciousness. Significance statement: Neurophysiology studies have found alpha/beta oscillations (8-30Hz), gamma oscillations (40-100Hz), and spiking activity during cognition. Alpha/beta power has an inverse relationship with gamma power/spiking. This inverse relationship suggests that gamma/spiking are under the inhibitory control of alpha/beta. The predictive routing model hypothesizes that alpha/beta oscillations selectively inhibit (and thereby control) cortical activity that is predictable. We tested whether this inhibitory control is a signature of consciousness. We used multi-area neurophysiology recordings in monkeys presented with tone sequences that varied in predictability. We recorded brain activity as the anesthetic propofol was administered to manipulate consciousness. Compared to conscious processing, propofol-mediated unconsciousness disrupted alpha/beta inhibitory control during predictive processing. This led to a disinhibition of gamma/spiking, consistent with the predictive routing model.

Modulatory dynamics mark the transition between anesthetic states of unconsciousness.

Unconsciousness maintained by GABAergic anesthetics, such as propofol and sevoflurane, is characterized by slow-delta oscillations (0.3 to 4 Hz) and alpha oscillations (8 to 14 Hz) that are readily visible in the electroencephalogram. At higher doses, these slow-delta-alpha (SDA) oscillations transition into burst suppression. This is a marker of a state of profound brain inactivation during which isoelectric (flatline) periods alternate with periods of the SDA patterns present at lower doses. While the SDA and burst suppression patterns have been analyzed separately, the transition from one to the other has not. Using state-space methods, we characterize the dynamic evolution of brain activity from SDA to burst suppression and back during unconsciousness maintained with propofol or sevoflurane in volunteer subjects and surgical patients. We uncover two dynamical processes that continuously modulate the SDA oscillations: alpha-wave amplitude and slow-wave frequency modulation. We present an alpha modulation index and a slow modulation index which characterize how these processes track the transition from SDA oscillations to burst suppression and back to SDA oscillations as a function of increasing and decreasing anesthetic doses, respectively. Our biophysical model reveals that these dynamics track the combined evolution of the neurophysiological and metabolic effects of a GABAergic anesthetic on brain circuits. Our characterization of the modulatory dynamics mediated by GABAergic anesthetics offers insights into the mechanisms of these agents and strategies for monitoring and precisely controlling the level of unconsciousness in patients under general anesthesia.

Sustaining wakefulness: Brainstem connectivity in human consciousness.

Consciousness is comprised of arousal (i.e., wakefulness) and awareness. Substantial progress has been made in mapping the cortical networks that modulate awareness in the human brain, but knowledge about the subcortical networks that sustain arousal is lacking. We integrated data from ex vivo diffusion MRI, immunohistochemistry, and in vivo 7 Tesla functional MRI to map the connectivity of a subcortical arousal network that we postulate sustains wakefulness in the resting, conscious human brain, analogous to the cortical default mode network (DMN) that is believed to sustain self-awareness. We identified nodes of the proposed default ascending arousal network (dAAN) in the brainstem, hypothalamus, thalamus, and basal forebrain by correlating ex vivo diffusion MRI with immunohistochemistry in three human brain specimens from neurologically normal individuals scanned at 600-750 µm resolution. We performed deterministic and probabilistic tractography analyses of the diffusion MRI data to map dAAN intra-network connections and dAAN-DMN internetwork connections. Using a newly developed network-based autopsy of the human brain that integrates ex vivo MRI and histopathology, we identified projection, association, and commissural pathways linking dAAN nodes with one another and with cortical DMN nodes, providing a structural architecture for the integration of arousal and awareness in human consciousness. We release the ex vivo diffusion MRI data, corresponding immunohistochemistry data, network-based autopsy methods, and a new brainstem dAAN atlas to support efforts to map the connectivity of human consciousness.

Recovery From Acute Respiratory Distress Syndrome Is Associated With Increasing Alpha Power in the Frontal Electroencephalogram During Propofol Sedation: A Case Report.

The effects of critical illness on electroencephalographic (EEG) signatures of sedatives have not been described, limiting the use of EEG-guided sedation in the intensive care unit (ICU). We report the case of a 36-year-old man recovering from acute respiratory distress syndrome (ARDS). Severe ARDS was characterized by slow-delta (0.1-4 Hz) and theta (4-8 Hz) oscillations but lacked the alpha (8-14 Hz) power expected during propofol sedation in a patient of this age. The alpha power emerged as ARDS resolved. This case raises the question of whether inflammatory states can alter EEG signatures during sedation.

Propofol Mediated Unconsciousness Disrupts Progression of Sensory Signals through the Cortical Hierarchy.

A critical component of anesthesia is the loss sensory perception. Propofol is the most widely used drug for general anesthesia, but the neural mechanisms of how and when it disrupts sensory processing are not fully understood. We analyzed local field potential (LFP) and spiking recorded from Utah arrays in auditory cortex, associative cortex, and cognitive cortex of non-human primates before and during propofol mediated unconsciousness. Sensory stimuli elicited robust and decodable stimulus responses and triggered periods of stimulus-induced coherence between brain areas in the LFP of awake animals. By contrast, propofol mediated unconsciousness eliminated stimulus-induced coherence and drastically weakened stimulus responses and information in all brain areas except for auditory cortex, where responses and information persisted. However, we found stimuli occurring during spiking Up states triggered weaker spiking responses than in awake animals in auditory cortex, and little or no spiking responses in higher order areas. These results suggest that propofol's effect on sensory processing is not just due to asynchronous down states. Rather, both Down states and Up states reflect disrupted dynamics.

Rapid thalamocortical network switching mediated by cortical synchronization underlies propofol-induced EEG signatures: a biophysical model.

Propofol-mediated unconsciousness elicits strong alpha/low-beta and slow oscillations in the electroencephalogram (EEG) of patients. As anesthetic dose increases, the EEG signal changes in ways that give clues to the level of unconsciousness; the network mechanisms of these changes are only partially understood. Here, we construct a biophysical thalamocortical network involving brain stem influences that reproduces transitions in dynamics seen in the EEG involving the evolution of the power and frequency of alpha/low-beta and slow rhythm, as well as their interactions. Our model suggests that propofol engages thalamic spindle and cortical sleep mechanisms to elicit persistent alpha/low-beta and slow rhythms, respectively. The thalamocortical network fluctuates between two mutually exclusive states on the timescale of seconds. One state is characterized by continuous alpha/low-beta-frequency spiking in thalamus (C-state), whereas in the other, thalamic alpha spiking is interrupted by periods of co-occurring thalamic and cortical silence (I-state). In the I-state, alpha colocalizes to the peak of the slow oscillation; in the C-state, there is a variable relationship between an alpha/beta rhythm and the slow oscillation. The C-state predominates near loss of consciousness; with increasing dose, the proportion of time spent in the I-state increases, recapitulating EEG phenomenology. Cortical synchrony drives the switch to the I-state by changing the nature of the thalamocortical feedback. Brain stem influence on the strength of thalamocortical feedback mediates the amount of cortical synchrony. Our model implicates loss of low-beta, cortical synchrony, and coordinated thalamocortical silent periods as contributing to the unconscious state.NEW & NOTEWORTHY GABAergic anesthetics induce alpha/low-beta and slow oscillations in the EEG, which interact in dose-dependent ways. We constructed a thalamocortical model to investigate how these interdependent oscillations change with propofol dose. We find two dynamic states of thalamocortical coordination, which change on the timescale of seconds and dose-dependently mirror known changes in EEG. Thalamocortical feedback determines the oscillatory coupling and power seen in each state, and this is primarily driven by cortical synchrony and brain stem neuromodulation.