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.

Artifact characterization and mitigation techniques during concurrent sensing and stimulation using bidirectional deep brain stimulation platforms.

Bidirectional deep brain stimulation (DBS) platforms have enabled a surge in hours of recordings in naturalistic environments, allowing further insight into neurological and psychiatric disease states. However, high amplitude, high frequency stimulation generates artifacts that contaminate neural signals and hinder our ability to interpret the data. This is especially true in psychiatric disorders, for which high amplitude stimulation is commonly applied to deep brain structures where the native neural activity is miniscule in comparison. Here, we characterized artifact sources in recordings from a bidirectional DBS platform, the Medtronic Summit RC + S, with the goal of optimizing recording configurations to improve signal to noise ratio (SNR). Data were collected from three subjects in a clinical trial of DBS for obsessive-compulsive disorder. Stimulation was provided bilaterally to the ventral capsule/ventral striatum (VC/VS) using two independent implantable neurostimulators. We first manipulated DBS amplitude within safe limits (2-5.3 mA) to characterize the impact of stimulation artifacts on neural recordings. We found that high amplitude stimulation produces slew overflow, defined as exceeding the rate of change that the analog to digital converter can accurately measure. Overflow led to expanded spectral distortion of the stimulation artifact, with a six fold increase in the bandwidth of the 150.6 Hz stimulation artifact from 147-153 to 140-180 Hz. By increasing sense blank values during high amplitude stimulation, we reduced overflow by as much as 30% and improved artifact distortion, reducing the bandwidth from 140-180 Hz artifact to 147-153 Hz. We also identified artifacts that shifted in frequency through modulation of telemetry parameters. We found that telemetry ratio changes led to predictable shifts in the center-frequencies of the associated artifacts, allowing us to proactively shift the artifacts outside of our frequency range of interest. Overall, the artifact characterization methods and results described here enable increased data interpretability and unconstrained biomarker exploration using data collected from bidirectional DBS devices.

Neural effects of propofol-induced unconsciousness and its reversal using thalamic stimulation.

The specific circuit mechanisms through which anesthetics induce unconsciousness have not been completely characterized. We recorded neural activity from the frontal, parietal, and temporal cortices and thalamus while maintaining unconsciousness in non-human primates (NHPs) with the anesthetic propofol. Unconsciousness was marked by slow frequency (~1 Hz) oscillations in local field potentials, entrainment of local spiking to Up states alternating with Down states of little or no spiking activity, and decreased coherence in frequencies above 4 Hz. Thalamic stimulation 'awakened' anesthetized NHPs and reversed the electrophysiologic features of unconsciousness. Unconsciousness is linked to cortical and thalamic slow frequency synchrony coupled with decreased spiking, and loss of higher-frequency dynamics. This may disrupt cortical communication/integration.

Layer and rhythm specificity for predictive routing.

In predictive coding, experience generates predictions that attenuate the feeding forward of predicted stimuli while passing forward unpredicted "errors." Different models have suggested distinct cortical layers, and rhythms implement predictive coding. We recorded spikes and local field potentials from laminar electrodes in five cortical areas (visual area 4 [V4], lateral intraparietal [LIP], posterior parietal area 7A, frontal eye field [FEF], and prefrontal cortex [PFC]) while monkeys performed a task that modulated visual stimulus predictability. During predictable blocks, there was enhanced alpha (8 to 14 Hz) or beta (15 to 30 Hz) power in all areas during stimulus processing and prestimulus beta (15 to 30 Hz) functional connectivity in deep layers of PFC to the other areas. Unpredictable stimuli were associated with increases in spiking and in gamma-band (40 to 90 Hz) power/connectivity that fed forward up the cortical hierarchy via superficial-layer cortex. Power and spiking modulation by predictability was stimulus specific. Alpha/beta power in LIP, FEF, and PFC inhibited spiking in deep layers of V4. Area 7A uniquely showed increases in high-beta (∼22 to 28 Hz) power/connectivity to unpredictable stimuli. These results motivate a conceptual model, predictive routing. It suggests that predictive coding may be implemented via lower-frequency alpha/beta rhythms that "prepare" pathways processing-predicted inputs by inhibiting feedforward gamma rhythms and associated spiking.

A simulation-based comparative analysis of PID and LQG control for closed-loop anesthesia delivery.

Closed loop anesthesia delivery (CLAD) systems can help anesthesiologists efficiently achieve and maintain desired anesthetic depth over an extended period of time. A typical CLAD system would use an anesthetic marker, calculated from physiological signals, as real-time feedback to adjust anesthetic dosage towards achieving a desired set-point of the marker. Since control strategies for CLAD vary across the systems reported in recent literature, a comparative analysis of common control strategies can be useful. For a nonlinear plant model based on well-established models of compartmental pharmacokinetics and sigmoid-Emax pharmacodynamics, we numerically analyze the set-point tracking performance of three output-feedback linear control strategies: proportional-integral-derivative (PID) control, linear quadratic Gaussian (LQG) control, and an LQG with integral action (ILQG). Specifically, we numerically simulate multiple CLAD sessions for the scenario where the plant model parameters are unavailable for a patient and the controller is designed based on a nominal model and controller gains are held constant throughout a session. Based on the numerical analyses performed here, conditioned on our choice of model and controllers, we infer that in terms of accuracy and bias PID control performs better than ILQG which in turn performs better than LQG. In the case of noisy observations, ILQG can be tuned to provide a smoother infusion rate while achieving comparable steady state response with respect to PID. The numerical analysis framework and findings reported here can help CLAD developers in their choice of control strategies. This paper may also serve as a tutorial paper for teaching control theory for CLAD.