Data-based modeling of cerebral hemodynamics quantifies impairment of cerebral blood flow regulation in type-2 diabetes.

We studied the regulation dynamics of cerebral blood velocity (CBv) at middle cerebral arteries (MCA) in response to spontaneous changes of arterial blood pressure (ABP), termed dynamic cerebral autoregulation (dCA), and end-tidal CO2 as proxy for blood CO2 tension, termed dynamic vasomotor reactivity (DVR), by analyzing time-series data collected at supine rest from 36 patients with Type-2 Diabetes Mellitus (T2DM) and 22 age/sex-matched non-diabetic controls without arterial hypertension. Our analysis employed a robust dynamic modeling methodology that utilizes Principal Dynamic Modes (PDM) to estimate subject-specific dynamic transformations of spontaneous changes in ABP and end-tidal CO2 (viewed as two "inputs") into changes of CBv at MCA measured via Transcranial Doppler ultrasound (viewed as the "output"). The quantitative results of PDM analysis indicate significant alterations in T2DM of both DVR and dCA in terms of two specific PDM contributions that rise to significance (p < 0.05). Our results further suggest that the observed DVR and dCA alterations may be due to reduction of cholinergic activity (based on previously published results from cholinergic blockade data) that may disturb the sympatho-vagal balance in T2DM. Combination of these two model-based "physio-markers" differentiated T2DM patients from controls (p = 0.0007), indicating diabetes-related alteration of cerebrovascular regulation, with possible diagnostic implications.

Impaired dynamic cerebral autoregulation measured in the middle cerebral artery in patients with vertebrobasilar ischemia is associated with autonomic failure.

OBJECTIVES: To assess whether vertebrobasilar artery ischemia (VBI) affects cortical cerebral blood flow (CBF) regulation. MATERIAL AND METHODS: 107 consecutive patients (mean age 65 ± 15 years; women 21) with VBI underwent structured stroke care with assessment of dynamic cerebral autoregulation (dCA) in both middle cerebral arteries (MCAs) by transfer function analysis using spontaneous oscillations of blood pressure (BP) and CBF velocity that yields by extraction of phase and gain information in the very low (0.02-0.07 Hz), low (0.07-0.15 Hz) and high frequency (0.15-0.5 Hz) ranges. Additionally, power spectrum analysis of BP and heart rate variability (HRV) was performed. The control group consists of 29 age- and sex-matched healthy persons. RESULTS: Compared to controls, phase in the VBI patients was significantly reduced and gain increased in the very low frequencies (VLF), in the low (LF), phase was significantly reduced only ipsilaterally. In the high frequencies (HF), phase reduction was only marginally significant. BP power spectral density (PSD) was much higher in the patients than in the controls across all frequencies. In the PSD of heart rate variability the controls but not the patients exhibited a strong peak around 0.11Hz, while the patients, but not the controls, exhibit a strong peak around 0.36 Hz. In regression analysis, patient's phase and gain results were not related to age, sex, arterial hypertension, diabetes mellitus, renal dysfunction, heart failure as indicated by left ventricular ejection fraction, stroke subtype, presence or absence of cerebral small vessel disease. CONCLUSION: Patients with VBI exhibit bilateral cortical autoregulation impairment in association with an autonomic nervous system disbalance. GOV IDENTIFIER: NCT04611672.

Dynamic effects of cholinergic blockade upon cerebral blood flow autoregulation in healthy adults.

Background: Cerebral flow autoregulation (CFA) is a homeostatic mechanism critical for survival. The autonomic nervous system (ANS) plays a key role in maintaining proper CFA function. More quantitative studies of how the ANS influences CFA are desirable. Objective: To discover and quantify the dynamic effects of cholinergic blockade upon CFA in response to changes of arterial blood pressure and blood CO2 tension in healthy adults. Methods: We analyzed time-series data of spontaneous beat-to-beat mean arterial blood pressure (ABP) and cerebral blood flow velocity in the middle cerebral arteries (CFV), as well as breath-to-breath end-tidal CO2 (CO2), collected in 9 adults before and after cholinergic blockade, in order to obtain subject-specific predictive input-output models of the dynamic effects of changes in ABP and CO2 (inputs) upon CFV (output). These models are defined in convolutional form using "kernel" functions (or, equivalently, Transfer Functions in the frequency domain) that are estimated via the robust method of Laguerre expansions. Results: Cholinergic blockade caused statistically significant changes in the obtained kernel estimates (and the corresponding Transfer Functions) that define the linear dynamics of the ABP-to-CFV and CO2-to-CFV causal relations. The kernel changes due to cholinergic blockade reflect the effects of the cholinergic mechanism and exhibited, in the frequency domain, resonant peaks at 0.22 Hz and 0.06 Hz for the ABP-to-CFV and CO2-to-CFV dynamics, respectively. Conclusion: Quantitative estimates of the dynamics of the cholinergic component in CFA are found as average changes of the ABP-to-CFV and CO2-to-CFV kernels, and corresponding Transfer Functions, before and after cholinergic blockade.

The Dynamic Relationship Between Cortical Oxygenation and End-Tidal CO 2 Transient Changes Is Impaired in Mild Cognitive Impairment Patients.

Background: Recent studies have utilized data-based dynamic modeling to establish strong association between dysregulation of cerebral perfusion and Mild Cognitive Impairment (MCI), expressed in terms of impaired CO2 dynamic vasomotor reactivity in the cerebral vasculature. This raises the question of whether this is due to dysregulation of central mechanisms (baroreflex and chemoreflex) or mechanisms of cortical tissue oxygenation (CTO) in MCI patients. We seek to answer this question using data-based input-output predictive dynamic models. Objective: To use subject-specific data-based multivariate input-output dynamic models to quantify the effects of systemic hemodynamic and blood CO2 changes upon CTO and to examine possible differences in CTO regulation in MCI patients versus age-matched controls, after the dynamic effects of central regulatory mechanisms have been accounted for by using cerebral flow measurements as another input. Methods: The employed model-based approach utilized the general dynamic modeling methodology of Laguerre expansions of kernels to analyze spontaneous time-series data in order to quantify the dynamic effects upon CTO (an index of relative capillary hemoglobin saturation distribution measured via near-infrared spectroscopy) of contemporaneous changes in end-tidal CO2 (proxy for arterial CO2), arterial blood pressure and cerebral blood flow velocity in the middle cerebral arteries (measured via transcranial Doppler). Model-based indices (physio-markers) were computed for these distinct dynamic relationships. Results: The obtained model-based indices revealed significant statistical differences of CO2 dynamic vasomotor reactivity in cortical tissue, combined with "perfusivity" that quantifies the dynamic relationship between flow velocity in cerebral arteries and CTO in MCI patients versus age-matched controls (p = 0.006). Significant difference between MCI patients and age-matched controls was also found in the respective model-prediction accuracy (p = 0.0001). Combination of these model-based indices via the Fisher Discriminant achieved even smaller p-value (p = 5 × 10-5) when comparing MCI patients with controls. The differences in dynamics of CTO in MCI patients are in lower frequencies (<0.05 Hz), suggesting impairment in endocrine/metabolic (rather than neural) mechanisms. Conclusion: The presented model-based approach elucidates the multivariate dynamic connectivity in the regulation of cerebral perfusion and yields model-based indices that may serve as physio-markers of possible dysregulation of CTO during transient CO2 changes in MCI patients.

Dysregulation of CO2-Driven Heart-Rate Chemoreflex Is Related Closely to Impaired CO2 Dynamic Vasomotor Reactivity in Mild Cognitive Impairment Patients.

BACKGROUND: Significant reduction of dynamic vasomotor reactivity (DVR) was recently reported in patients with amnestic mild cognitive impairment (MCI) relative to age-matched controls. These results were obtained via a novel approach that utilizes data-based predictive dynamic models to quantify DVR. OBJECTIVE: Using the same methodological approach, we seek to quantify the dynamic effects of the CO2-driven chemoreflex and baroreflex upon heart-rate in order to examine their possible correlation with the observed DVR impairment in each MCI patient. METHODS: The employed approach utilizes time-series data to obtain subject-specific predictive input-output models of the dynamic effects of changes in arterial blood pressure and end-tidal CO2 (putative "inputs") upon cerebral blood flow velocity in large cerebral arteries, cortical tissue oxygenation, and heart-rate (putative "outputs"). RESULTS: There was significant dysregulation of CO2-driven heart-rate chemoreflex (p = 0.0031), but not of baroreflex (p = 0.5061), in MCI patients relative to age-matched controls. The model-based index of CO2-driven heart-rate chemoreflex gain (CRG) correlated significantly with the DVR index in large cerebral arteries (p = 0.0146), but not with the DVR index in small/micro-cortical vessels (p = 0.1066). This suggests that DVR impairment in small/micro-cortical vessels is not mainly due to CO2-driven heart-rate chemoreflex dysregulation, but to other factors (possibly dysfunction of neurovascular coupling). CONCLUSION: Improved delineation between MCI patients and controls is achieved by combining the DVR index for small/micro-cortical vessels with the CRG index (p = 2×10-5). There is significant correlation (p < 0.01) between neuropsychological test scores and model-based DVR indices. Combining neuropsychological scores with DVR indices reduces the composite diagnostic index p-value (p∼10-10).

Assessment of dynamic cerebral autoregulation in humans: Is reproducibility dependent on blood pressure variability?

We tested the influence of blood pressure variability on the reproducibility of dynamic cerebral autoregulation (DCA) estimates. Data were analyzed from the 2nd CARNet bootstrap initiative, where mean arterial blood pressure (MABP), cerebral blood flow velocity (CBFV) and end tidal CO2 were measured twice in 75 healthy subjects. DCA was analyzed by 14 different centers with a variety of different analysis methods. Intraclass Correlation (ICC) values increased significantly when subjects with low power spectral density MABP (PSD-MABP) values were removed from the analysis for all gain, phase and autoregulation index (ARI) parameters. Gain in the low frequency band (LF) had the highest ICC, followed by phase LF and gain in the very low frequency band. No significant differences were found between analysis methods for gain parameters, but for phase and ARI parameters, significant differences between the analysis methods were found. Alternatively, the Spearman-Brown prediction formula indicated that prolongation of the measurement duration up to 35 minutes may be needed to achieve good reproducibility for some DCA parameters. We conclude that poor DCA reproducibility (ICC<0.4) can improve to good (ICC > 0.6) values when cases with low PSD-MABP are removed, and probably also when measurement duration is increased.

Quantification of dynamic cerebral autoregulation and CO2 dynamic vasomotor reactivity impairment in essential hypertension.

The study of dynamic cerebral autoregulation (DCA) in essential hypertension has received considerable attention because of its clinical importance. Several studies have examined the dynamic relationship between spontaneous beat-to-beat arterial blood pressure data and contemporaneous cerebral blood flow velocity measurements (obtained via transcranial Doppler at the middle cerebral arteries) in the form of a linear input-output model using transfer function analysis. This analysis is more reliable when the contemporaneous effects of changes in blood CO2 tension are also taken into account, because of the significant effects of CO2 dynamic vasomotor reactivity (DVR) upon cerebral flow. In this article, we extract such input-output predictive models from spontaneous time series hemodynamic data of 24 patients with essential hypertension and 20 normotensive control subjects under resting conditions, using the novel methodology of principal dynamic modes (PDMs) that achieves improved estimation accuracy over previous methods for relatively short and noisy data. The obtained data-based models are subsequently used to compute indexes and markers that quantify DCA and DVR in each subject or patient and therefore can be used to assess the effects of essential hypertension. These model-based DCA and DVR indexes were properly defined to capture the observed effects of DCA and VR and found to be significantly different (P < 0.05) in the hypertensive patients. We also found significant differences between patients and control subjects in the relative contribution of three PDMs to the model output prediction, a finding that offers the prospect of identifying the physiological mechanisms affected by essential hypertension when the PDMs are interpreted in terms of specific physiological mechanisms.NEW & NOTEWORTHY This article presents novel model-based methodology for obtaining diagnostic indexes of dynamic cerebral autoregulation and dynamic vasomotor reactivity in hypertension.

Dynamic Cerebral Autoregulation Reproducibility Is Affected by Physiological Variability.

Parameters describing dynamic cerebral autoregulation (DCA) have limited reproducibility. In an international, multi-center study, we evaluated the influence of multiple analytical methods on the reproducibility of DCA. Fourteen participating centers analyzed repeated measurements from 75 healthy subjects, consisting of 5 min of spontaneous fluctuations in blood pressure and cerebral blood flow velocity signals, based on their usual methods of analysis. DCA methods were grouped into three broad categories, depending on output types: (1) transfer function analysis (TFA); (2) autoregulation index (ARI); and (3) correlation coefficient. Only TFA gain in the low frequency (LF) band showed good reproducibility in approximately half of the estimates of gain, defined as an intraclass correlation coefficient (ICC) of >0.6. None of the other DCA metrics had good reproducibility. For TFA-like and ARI-like methods, ICCs were lower than values obtained with surrogate data (p < 0.05). For TFA-like methods, ICCs were lower for the very LF band (gain 0.38 ± 0.057, phase 0.17 ± 0.13) than for LF band (gain 0.59 ± 0.078, phase 0.39 ± 0.11, p ≤ 0.001 for both gain and phase). For ARI-like methods, the mean ICC was 0.30 ± 0.12 and for the correlation methods 0.24 ± 0.23. Based on comparisons with ICC estimates obtained from surrogate data, we conclude that physiological variability or non-stationarity is likely to be the main reason for the poor reproducibility of DCA parameters.

Comparing model-based cerebrovascular physiomarkers with DTI biomarkers in MCI patients.

OBJECTIVE: To compare the novel model-based hemodynamic physiomarker of Dynamic Vasomotor Reactivity (DVR) with biomarkers based on Diffusion Tensor Imaging (DTI) and some widely used neurocognitive scores in terms of their ability to delineate patients with amnestic Mild Cognitive Impairment (MCI) from age-matched cognitively normal controls. MATERIALS & METHODS: The model-based DVR and MRI-based DTI markers were obtained from 36 patients with amnestic MCI and 16 age-matched controls without cognitive impairment, for whom widely used neurocognitive scores were available. These markers and scores were subsequently compared in terms of statistical delineation between patients and controls. RESULTS: It was found that statistically significant delineation between MCI patients and controls was comparable for DVR or DTI markers (p < 0.01). The performance of both types of markers was consistent with the scores of some (but not all) widely used neurocognitive tests. CONCLUSION: Since DTI offers a measure of cerebral white matter integrity, the results suggest that the model-based hemodynamic marker of DVR may correlate with cognitive impairment due to white matter lesions. This finding is consistent with the hypothesis that dysregulation of cerebral microcirculation may be an early cause of cognitive impairment, which has been recently corroborated by several studies.

Multi-Input, Multi-Output Neuronal Mode Network Approach to Modeling the Encoding Dynamics and Functional Connectivity of Neural Systems.

This letter proposes a novel method, multi-input, multi-output neuronal mode network (MIMO-NMN), for modeling encoding dynamics and functional connectivity in neural ensembles such as the hippocampus. Compared with conventional approaches such as the Volterra-Wiener model, linear-nonlinear-cascade (LNC) model, and generalized linear model (GLM), the NMN has several advantages in terms of estimation accuracy, model interpretation, and functional connectivity analysis. We point out the limitations of current neural spike modeling methods, especially the estimation biases caused by the imbalanced class problem when the number of zeros is significantly larger than ones in the spike data. We use synthetic data to test the performance of NMN with a comparison of the traditional methods, and the results indicate the NMN approach could reduce the imbalanced class problem and achieve better predictions. Subsequently, we apply the MIMO-NMN method to analyze data from the human hippocampus. The results indicate that the MIMO-NMN method is a promising approach to modeling neural dynamics and analyzing functional connectivity of multi-neuronal data.

Closed-loop modeling of the heart-rate reflex for improved diagnosis and monitoring of Mild Cognitive Impairment.

Analysis of beat-to-beat spontaneous cerebral hemodynamic data has yielded predictive dynamic models of cerebral hemodynamics and has shown previously that patients with Mild Cognitive Impairment (MCI) exhibit significantly reduced cerebral vasomotor reactivity to CO2 relative to cognitively normal control subjects [1]. The present work examines the heart-rate reflex (HRR) dynamics of 46 MCI patients compared to 20 control subjects, using closed-loop modeling of HRR under resting conditions of spontaneous variations of arterial blood pressure (baroreflex) and end-tidal CO2 (chemoreflex). These subject-specific predictive dynamic models are obtained via the methodology of Principal Dynamic Modes [2] and allow the computation of model-based markers of baroreflex and chemoreflex function. We found that the chemoreflex gain is significantly weakened in MCI patients relative to controls (p=0.0086), while the baroreflex is not significantly affected. These findings offer another tool for diagnosis and monitoring of MCI (via model-based markers), when used in conjunction with current methods.

Reproducibility of dynamic cerebral autoregulation parameters: a multi-centre, multi-method study.

OBJECTIVE: Different methods to calculate dynamic cerebral autoregulation (dCA) parameters are available. However, most of these methods demonstrate poor reproducibility that limit their reliability for clinical use. Inter-centre differences in study protocols, modelling approaches and default parameter settings have all led to a lack of standardisation and comparability between studies. We evaluated reproducibility of dCA parameters by assessing systematic errors in surrogate data resulting from different modelling techniques. APPROACH: Fourteen centres analysed 22 datasets consisting of two repeated physiological blood pressure measurements with surrogate cerebral blood flow velocity signals, generated using Tiecks curves (autoregulation index, ARI 0-9) and added noise. For reproducibility, dCA methods were grouped in three broad categories: 1. Transfer function analysis (TFA)-like output; 2. ARI-like output; 3. Correlation coefficient-like output. For all methods, reproducibility was determined by one-way intraclass correlation coefficient analysis (ICC). MAIN RESULTS: For TFA-like methods the mean (SD; [range]) ICC gain was 0.71 (0.10; [0.49-0.86]) and 0.80 (0.17; [0.36-0.94]) for VLF and LF (p = 0.003) respectively. For phase, ICC values were 0.53 (0.21; [0.09-0.80]) for VLF, and 0.92 (0.13; [0.44-1.00]) for LF (p < 0.001). Finally, ICC for ARI-like methods was equal to 0.84 (0.19; [0.41-0.94]), and for correlation-like methods, ICC was 0.21 (0.21; [0.056-0.35]). SIGNIFICANCE: When applied to realistic surrogate data, free from the additional exogenous influences of physiological variability on cerebral blood flow, most methods of dCA modelling showed ICC values considerably higher than what has been reported for physiological data. This finding suggests that the poor reproducibility reported by previous studies may be mainly due to the inherent physiological variability of cerebral blood flow regulatory mechanisms rather than related to (stationary) random noise and the signal analysis methods.

Mechanism-Based and Input-Output Modeling of the Key Neuronal Connections and Signal Transformations in the CA3-CA1 Regions of the Hippocampus.

This letter examines the results of input-output (nonparametric) modeling based on the analysis of data generated by a mechanism-based (parametric) model of CA3-CA1 neuronal connections in the hippocampus. The motivation is to obtain biological insight into the interpretation of such input-output (Volterra-equivalent) models estimated from synthetic data. The insights obtained may be subsequently used to interpretat input-output models extracted from actual experimental data. Specifically, we found that a simplified parametric model may serve as a useful tool to study the signal transformations in the hippocampal CA3-CA1 regions. Input-output modeling of model-based synthetic data show that GABAergic interneurons are responsible for regulating neuronal excitation, controlling the precision of spike timing, and maintaining network oscillations, in a manner consistent with previous studies. The input-output model obtained from real data exhibits intriguing similarities with its synthetic-data counterpart, demonstrating the importance of a dynamic resonance in the system/model response around 2 Hz to 3 Hz. Using the input-output model from real data as a guide, we may be able to amend the parametric model by incorporating more mechanisms in order to yield better-matching input-output model. The approach we present can also be applied to the study of other neural systems and pathways.

Compartmental and Data-Based Modeling of Cerebral Hemodynamics: Nonlinear Analysis.

OBJECTIVE: As an extension to our study comparing a putative compartmental and data-based model of linear dynamic cerebral autoregulation (CA) and CO2-vasomotor reactivity (VR), we study the CA-VR process in a nonlinear context. METHODS: We use the concept of principal dynamic modes (PDM) in order to obtain a compact and more easily interpretable input-output model. This in silico study permits the use of input data with a dynamic range large enough to simulate the classic homeostatic CA and VR curves using a putative structural model of the regulatory control of the cerebral circulation. The PDM model obtained using theoretical and experimental data are compared. RESULTS: It was found that the PDM model was able to reflect accurately both the simulated static CA and VR curves in the associated nonlinear functions (ANFs). Similar to experimental observations, the PDM model essentially separates the pressure-flow relationship into a linear component with fast dynamics and nonlinear components with slow dynamics. In addition, we found good qualitative agreement between the PDMs representing the dynamic theoretical and experimental CO2-flow relationship. CONCLUSION: Under the modeling assumption and in light of other experimental findings, we hypothesize that PDMs obtained from experimental data correspond with passive fluid dynamical and active regulatory mechanisms. SIGNIFICANCE: Both hypothesis-based and data-based modeling approaches can be combined to offer some insight into the physiological basis of PDM model obtained from human experimental data. The PDM modeling approach potentially offers a practical way to quantify the status of specific regulatory mechanisms in the CA-VR process.

Comparison of Model-Based Indices of Cerebral Autoregulation and Vasomotor Reactivity Using Transcranial Doppler versus Near-Infrared Spectroscopy in Patients with Amnestic Mild Cognitive Impairment.

We recently introduced model-based "physiomarkers" of dynamic cerebral autoregulation and CO2 vasomotor reactivity as an aid for diagnosis of early-stage Alzheimer's disease (AD) [1], where significant impairment of dynamic vasomotor reactivity (DVR) was observed in early-stage AD patients relative to age-matched controls. Milder impairment of DVR was shown in patients with amnestic mild cognitive impairment (MCI) using the same approach in a subsequent study [2]. The advocated approach utilizes subject-specific data-based models of cerebral hemodynamics to quantify the dynamic effects of resting-state changes in arterial blood pressure and end-tidal CO2 (the putative inputs) upon cerebral blood flow velocity (the putative output) measured at the middle cerebral artery via transcranial Doppler (TCD). The obtained input-output models are then used to compute model-based indices of DCA and DVR from model-predicted responses to an input pressure pulse or an input CO2 pulse, respectively. In this paper, we compare these model-based indices of DVR and DCA in 46 amnestic MCI patients, relative to 20 age-matched controls, using TCD measurements with their counterparts using Near-Infrared Spectroscopy (NIRS) measurements of blood oxygenation at the lateral prefrontal cortex in 43 patients and 22 age-matched controls. The goal of the study is to assess whether NIRS measurements can be used instead of TCD measurements to obtain model-based physiomarkers with comparable diagnostic utility. The results corroborate this view in terms of the ability of either output to yield model-based physiomarkers that can differentiate the group of aMCI patients from age-matched healthy controls.

Assessment of cerebrovascular reactivity during resting state breathing and its correlation with cognitive function in hypertension.

BACKGROUND: Hypertension is associated with cognitive deficits, particularly executive function, and decreased cerebral microvascular responsiveness to CO2 (CO2 vasoreactivity). The relation between CO2 vasoreactivity and executive function is not known. Protocols to assess CO2 vasoreactivity are cumbersome and require inhaling a CO2-enriched gas. We explored the ability to measure CO2 vasoreactivity using end-tidal CO2 fluctuations during normal breathing and the association of this measure with cognitive function in hypertension. METHODS: Executive function (Trail-Making Test parts A/B), memory, attention and blood flow velocity (BFV) in the middle cerebral artery using transcranial Doppler were measured in hypertensive subjects who were tapered off their treatment for 3 weeks. BFV was measured while sitting and normally breathing for 5 min, followed by breathing 5% CO2 gas and hyperventilation for 2 min each. We calculated CO2 vasoreactivity as the rate of BFV change from hypoventilation to hyperventilation, and as a model-derived measure using the normal breathing data. The latter was derived using nonlinear principal dynamic modes (PDM), which modelled the dynamic effect of fluctuations in end-tidal CO2 and blood pressure upon BFV during normal room-air respiration. Multiple regression analyses were used to correlate cerebral hemodynamics with cognitive measures. RESULTS: Data were collected from 41 individuals with hypertension (mean age 71 years, 24% African Americans, 61% women, off antihypertensive therapy). Lower CO2 vasoreactivity was associated with a worse executive function test score using both calculation methods: p value using the hyper/hypoventilation data was 0.04 and from the PDM analysis was 0.009. PDM calculations showed a stronger correlation with executive function (0.41 vs. 0.21 using the hyper/hypoventilation data). There were no associations with memory or attention measures. There was a weak but statistically significant correlation between the two calculation methods of CO2 vasoreactivity (R(2) = 14%, p = 0.02). CONCLUSION: This study suggests that the decrease in CO2 vasoreactivity in hypertension is associated with lower executive function. This may offer new insight into the vascular underpinning of cognitive decline in hypertension. We demonstrate that calculating CO2 vasoreactivity is possible during normal breathing. If replicated in future studies, this may offer a more convenient clinical way to assess CO2 vasoreactivity in hypertension and cognitive disorders.

Time-varying modeling of cerebral hemodynamics.

The scientific and clinical importance of cerebral hemodynamics has generated considerable interest in their quantitative understanding via computational modeling. In particular, two aspects of cerebral hemodynamics, cerebral flow autoregulation (CFA) and CO2 vasomotor reactivity (CVR), have attracted much attention because they are implicated in many important clinical conditions and pathologies (orthostatic intolerance, syncope, hypertension, stroke, vascular dementia, mild cognitive impairment, Alzheimer's disease, and other neurodegenerative diseases with cerebrovascular components). Both CFA and CVR are dynamic physiological processes by which cerebral blood flow is regulated in response to fluctuations in cerebral perfusion pressure and blood CO2 tension. Several modeling studies to date have analyzed beat-to-beat hemodynamic data in order to advance our quantitative understanding of CFA-CVR dynamics. A confounding factor in these studies is the fact that the dynamics of the CFA-CVR processes appear to vary with time (i.e., changes in cerebrovascular characteristics) due to neural, endocrine, and metabolic effects. This paper seeks to address this issue by tracking the changes in linear time-invariant models obtained from short successive segments of data from ten healthy human subjects. The results suggest that systemic variations exist but have stationary statistics and, therefore, the use of time-invariant modeling yields "time-averaged models" of physiological and clinical utility.

Facilitation of memory encoding in primate hippocampus by a neuroprosthesis that promotes task-specific neural firing.

OBJECTIVE: Memory accuracy is a major problem in human disease and is the primary factor that defines Alzheimer's, ageing and dementia resulting from impaired hippocampal function in the medial temporal lobe. Development of a hippocampal memory neuroprosthesis that facilitates normal memory encoding in nonhuman primates (NHPs) could provide the basis for improving memory in human disease states. APPROACH: NHPs trained to perform a short-term delayed match-to-sample (DMS) memory task were examined with multi-neuron recordings from synaptically connected hippocampal cell fields, CA1 and CA3. Recordings were analyzed utilizing a previously developed nonlinear multi-input multi-output (MIMO) neuroprosthetic model, capable of extracting CA3-to-CA1 spatiotemporal firing patterns during DMS performance. MAIN RESULTS: The MIMO model verified that specific CA3-to-CA1 firing patterns were critical for the successful encoding of sample phase information on more difficult DMS trials. This was validated by the delivery of successful MIMO-derived encoding patterns via electrical stimulation to the same CA1 recording locations during the sample phase which facilitated task performance in the subsequent, delayed match phase, on difficult trials that required more precise encoding of sample information. SIGNIFICANCE: These findings provide the first successful application of a neuroprosthesis designed to enhance and/or repair memory encoding in primate brain.

Analysis of intravenous glucose tolerance test data using parametric and nonparametric modeling: application to a population at risk for diabetes.

BACKGROUND: Modeling studies of the insulin-glucose relationship have mainly utilized parametric models, most notably the minimal model (MM) of glucose disappearance. This article presents results from the comparative analysis of the parametric MM and a nonparametric Laguerre based Volterra Model (LVM) applied to the analysis of insulin modified (IM) intravenous glucose tolerance test (IVGTT) data from a clinical study of gestational diabetes mellitus (GDM). METHODS: An IM IVGTT study was performed 8 to 10 weeks postpartum in 125 women who were diagnosed with GDM during their pregnancy [population at risk of developing diabetes (PRD)] and in 39 control women with normal pregnancies (control subjects). The measured plasma glucose and insulin from the IM IVGTT in each group were analyzed via a population analysis approach to estimate the insulin sensitivity parameter of the parametric MM. In the nonparametric LVM analysis, the glucose and insulin data were used to calculate the first-order kernel, from which a diagnostic scalar index representing the integrated effect of insulin on glucose was derived. RESULTS: Both the parametric MM and nonparametric LVM describe the glucose concentration data in each group with good fidelity, with an improved measured versus predicted r² value for the LVM of 0.99 versus 0.97 for the MM analysis in the PRD. However, application of the respective diagnostic indices of the two methods does result in a different classification of 20% of the individuals in the PRD. CONCLUSIONS: It was found that the data based nonparametric LVM revealed additional insights about the manner in which infused insulin affects blood glucose concentration.

Nonlinear modeling of dynamic interactions within neuronal ensembles using Principal Dynamic Modes.

A methodology for nonlinear modeling of multi-input multi-output (MIMO) neuronal systems is presented that utilizes the concept of Principal Dynamic Modes (PDM). The efficacy of this new methodology is demonstrated in the study of the dynamic interactions between neuronal ensembles in the Pre-Frontal Cortex (PFC) of a behaving non-human primate (NHP) performing a Delayed Match-to-Sample task. Recorded spike trains from Layer-2 and Layer-5 neurons were viewed as the "inputs" and "outputs", respectively, of a putative MIMO system/model that quantifies the dynamic transformation of multi-unit neuronal activity between Layer-2 and Layer-5 of the PFC. Model prediction performance was evaluated by means of computed Receiver Operating Characteristic (ROC) curves. The PDM-based approach seeks to reduce the complexity of MIMO models of neuronal ensembles in order to enable the practicable modeling of large-scale neural systems incorporating hundreds or thousands of neurons, which is emerging as a preeminent issue in the study of neural function. The "scaling-up" issue has attained critical importance as multi-electrode recordings are increasingly used to probe neural systems and advance our understanding of integrated neural function. The initial results indicate that the PDM-based modeling methodology may greatly reduce the complexity of the MIMO model without significant degradation of performance. Furthermore, the PDM-based approach offers the prospect of improved biological/physiological interpretation of the obtained MIMO models.