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.

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.

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.

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.

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.

Cortical-subcortical interactions in goal-directed behavior.

Flexibly selecting appropriate actions in response to complex, ever-changing environments requires both cortical and subcortical regions, which are typically described as participating in a strict hierarchy. In this traditional view, highly specialized subcortical circuits allow for efficient responses to salient stimuli, at the cost of adaptability and context specificity, which are attributed to the neocortex. Their interactions are often described as the cortex providing top-down command signals for subcortical structures to implement; however, as available technologies develop, studies increasingly demonstrate that behavior is represented by brainwide activity and that even subcortical structures contain early signals of choice, suggesting that behavioral functions emerge as a result of different regions interacting as truly collaborative networks. In this review, we discuss the field's evolving understanding of how cortical and subcortical regions in placental mammals interact cooperatively, not only via top-down cortical-subcortical inputs but through bottom-up interactions, especially via the thalamus. We describe our current understanding of the circuitry of both the cortex and two exemplar subcortical structures, the superior colliculus and striatum, to identify which information is prioritized by which regions. We then describe the functional circuits these regions form with one another, and the thalamus, to create parallel loops and complex networks for brainwide information flow. Finally, we challenge the classic view that functional modules are contained within specific brain regions; instead, we propose that certain regions prioritize specific types of information over others, but the subnetworks they form, defined by their anatomical connections and functional dynamics, are the basis of true specialization.

Algebraic approach for subspace decomposition and clustering of neural activity.

We developed an approach to decompose neuronal signals into disjoint components, corresponding to task- or event-based epochs. This protocol describes how to project behavioral templates onto a low-dimensional subspace of neuronal responses to derive neuronal templates, then how to decompose and cluster neuronal responses using these derived templates. We outline these steps on complementary datasets of calcium imaging and spiking activity. Our approach relies on fundamental, linear algebraic principles and is adaptive to the temporal structure of the neural data. For complete details on the use and execution of this protocol, please refer to Adam et al. (2022).1.

Dynamic control of visually guided locomotion through corticosubthalamic projections.

Goal-directed locomotion requires control signals that propagate from higher order areas to regulate spinal mechanisms. The corticosubthalamic hyperdirect pathway offers a short route for cortical information to reach locomotor centers in the brainstem. We developed a task in which head-fixed mice run to a visual landmark and then stop and wait to collect the reward and examined the role of secondary motor cortex (M2) projections to the subthalamic nucleus (STN) in controlling locomotion. Our behavioral modeling, calcium imaging, and optogenetics manipulation results suggest that the M2-STN pathway can be recruited during visually guided locomotion to rapidly and precisely control the pedunculopontine nucleus (PPN) of the mesencephalic locomotor region through the basal ganglia. By capturing the physiological dynamics through a feedback control model and analyzing neuronal signals in M2, PPN, and STN, we find that the corticosubthalamic projections potentially control PPN activity by differentiating an M2 error signal to ensure fast input-output dynamics.

Deep brain stimulation in the subthalamic nucleus for Parkinson's disease can restore dynamics of striatal networks.

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is highly effective in alleviating movement disability in patients with Parkinson's disease (PD). However, its therapeutic mechanism of action is unknown. The healthy striatum exhibits rich dynamics resulting from an interaction of beta, gamma, and theta oscillations. These rhythms are essential to selection and execution of motor programs, and their loss or exaggeration due to dopamine (DA) depletion in PD is a major source of behavioral deficits. Restoring the natural rhythms may then be instrumental in the therapeutic action of DBS. We develop a biophysical networked model of a BG pathway to study how abnormal beta oscillations can emerge throughout the BG in PD and how DBS can restore normal beta, gamma, and theta striatal rhythms. Our model incorporates STN projections to the striatum, long known but understudied, found to preferentially target fast-spiking interneurons (FSI). We find that DBS in STN can normalize striatal medium spiny neuron activity by recruiting FSI dynamics and restoring the inhibitory potency of FSIs observed in normal conditions. We also find that DBS allows the reexpression of gamma and theta rhythms, thought to be dependent on high DA levels and thus lost in PD, through cortical noise control. Our study highlights that DBS effects can go beyond regularizing BG output dynamics to restoring normal internal BG dynamics and the ability to regulate them. It also suggests how gamma and theta oscillations can be leveraged to supplement DBS treatment and enhance its effectiveness.

Distinct prefrontal top-down circuits differentially modulate sensorimotor behavior.

Sensorimotor behaviors require processing of behaviorally relevant sensory cues and the ability to select appropriate responses from a vast behavioral repertoire. Modulation by the prefrontal cortex (PFC) is thought to be key for both processes, but the precise role of specific circuits remains unclear. We examined the sensorimotor function of anatomically distinct outputs from a subdivision of the mouse PFC, the anterior cingulate cortex (ACC). Using a visually guided two-choice behavioral paradigm with multiple cue-response mappings, we dissociated the sensory and motor response components of sensorimotor control. Projection-specific two-photon calcium imaging and optogenetic manipulations show that ACC outputs to the superior colliculus, a key midbrain structure for response selection, principally coordinate specific motor responses. Importantly, ACC outputs exert control by reducing the innate response bias of the superior colliculus. In contrast, ACC outputs to the visual cortex facilitate sensory processing of visual cues. Our results ascribe motor and sensory roles to ACC projections to the superior colliculus and the visual cortex and demonstrate for the first time a circuit motif for PFC function wherein anatomically non-overlapping output pathways coordinate complementary but distinct aspects of visual sensorimotor behavior.