Modeling of autonomic control in sleep-disordered breathing

Unanticipated evolution of cardio-respiratory interactions with cognitive load during a Go-NoGo shooting task in virtual reality

The cardiovascular system interacts continuously with the respiratory system to maintain the vital balance of oxygen and carbon dioxide in our body. The interplay between the sympathetic and parasympathetic branches of the autonomic nervous system regulates the aforesaid involuntary functions. This study analyzes the dynamics of the cardio-respiratory (CR) interactions using RR Intervals (RRI), Systolic Blood Pressure (SBP), and Respiration signals after first-order differencing to make them stationary. It investigates their variation with cognitive load induced by a virtual reality (VR) based Go-NoGo shooting task with low and high levels of task difficulty. We use Pearson's correlation-based linear and mutual information-based nonlinear measures of association to indicate the reduction in RRI-SBP and RRI-Respiration interactions with cognitive load. However, no linear correlation difference was observed in SBP-Respiration interactions with cognitive load, but their mutual information increased. A couple of open-loop autoregressive models with exogenous input (ARX) are estimated using RRI and SBP, and one closed-loop ARX model is estimated using RRI, SBP, and Respiration. The impulse responses (IRs) are derived for each input-output pair, and a reduction in the positive and negative peak amplitude of all the IRs is observed with cognitive load. Some novel parameters are derived by representing the IR as a double exponential curve with cosine modulation and show significant differences with cognitive load compared to other measures, especially for the IR between SBP and Respiration.

Baroreflex Sensitivity and Blood Pressure Variability can Help in Understanding the Different Response to Therapy During Acute Phase of Septic Shock

BACKGROUND

Mean values of hemodynamic variables are poorly effective in evaluating an actual recovery of the short-term autonomic mechanisms for blood pressure (BP) and heart rate (HR) regulation. The aim of this work is to analyze the response to therapy in the early phase of septic shock to verify possible associations between BP recovery and BP autonomic control.

METHODS

This is an ancillary study from the multicenter prospective observational trial Shockomics (NCT02141607). A total of 21 septic shock patients were studied at two time points during the acute phase of shock and were classified according to changes in SOFA score. Time series of BP components and HR were analyzed in time and frequency domain. Baroreflex sensitivity (BRS) was assessed, and a mathematical model for the decomposition of diastolic arterial pressure (DAP) oscillations was used to understand the different contributions of BRS and HR on peripheral vascular resistance control.

RESULTS

Only those patients, who significantly improved organ function (responders, R), showed an increase of mean value and low frequency (LF) power in BP time series. Fluid accumulation was higher in the non-responders (NR). BRS increased in NR and the model of DAP variability showed that the contribution of HR was highly reduced in NR.

CONCLUSIONS

Although patients reached the mean BP target of 65 mmHg, our analyses highlighted important differences in terms of autonomic nervous system control. BP variability, HR variability and baroreflex trends can add information to individual vital sign measure such as mean BP, and can help in understanding the responsiveness to the combination of symphatomimetic drugs and fluid therapy.

Biophysical markers of the peripheral vasoconstriction response to pain in sickle cell disease

Painful vaso-occlusive crisis (VOC), a complication of sickle cell disease (SCD), occurs when sickled red blood cells obstruct flow in the microvasculature. We postulated that exaggerated sympathetically mediated vasoconstriction, endothelial dysfunction and the synergistic interaction between these two factors act together to reduce microvascular flow, promoting regional vaso-occlusions, setting the stage for VOC. We previously found that SCD subjects had stronger vasoconstriction response to pulses of heat-induced pain compared to controls but the relative degrees to which autonomic dysregulation, peripheral vascular dysfunction and their interaction are present in SCD remain unknown. In the present study, we employed a mathematical model to decompose the total vasoconstriction response to pain into: 1) the neurogenic component, 2) the vascular response to blood pressure, 3) respiratory coupling and 4) neurogenic-vascular interaction. The model allowed us to quantify the contribution of each component to the total vasoconstriction response. The most salient features of the components were extracted to represent biophysical markers of autonomic and vascular impairment in SCD and controls. These markers provide a means of phenotyping severity of disease in sickle-cell anemia that is based more on underlying physiology than on genotype. The marker of the vascular component (BM_v) showed stronger contribution to vasoconstriction in SCD than controls (p = 0.0409), suggesting a dominant myogenic response in the SCD subjects as a consequence of endothelial dysfunction. The marker of neurogenic-vascular interaction (BM_n-v) revealed that the interaction reinforced vasoconstriction in SCD but produced vasodilatory response in controls (p = 0.0167). This marked difference in BM_n-v suggests that it is the most sensitive marker for quantifying combined alterations in autonomic and vascular function in SCD in response to heat-induced pain.

The possible role of the vagal nervous system in the recovery of the blood pressure control after cardiac arrest: a porcine model study

Previous studies have proved that the baroreceptor reflex (baroreflex) control of heart rate can be used for stratification of post-infarction population and, in general, cardiovascular disease populations. Many methods have been proposed to estimate the so-called baroreflex sensitivity (BRS) expressed as ms mmhg^-1. Most of the studies that exploit BRS focus mainly on acute myocardial infarction (AMI) and there are no important works that investigate the role of BRS immediately after cardiac arrest (CA). The present work is a continuation of the published work of Ristagno et al (2014 Shock 41 72-8). In particular, the main objectives are: (1) to study the evolution of BRS after CA and following cardiopulmonary resuscitation (CPR); (2) to verify if the recovery of cardiovascular stability and arterial blood pressure is accompanied by a recovery of BR in a porcine model; (3) to investigate the possible causes of the BRS variations in response to CA and following cardiopulmonary resuscitation. All the BRS estimators adopted in this study show a significant decrease after CA. However, partial recovery is obtained in the last hours of post resuscitation. Analysis of impulse response showed a decrease in peak delay after CA and was significantly shorter 4 hours after CPR. This finding hints at a compensation mechanism: a faster response when baroreflex gain is not fully restored. The increase in the speed of baroreflex response is in line with the hypothesis of a key role of the parasympathetic nervous system, which is known to act at a higher firing rate.

A Model-Based Machine Learning Approach to Probing Autonomic Regulation From Nonstationary Vital-Sign Time Series

Physiological variables, such as heart rate (HR), blood pressure (BP) and respiration (RESP), are tightly regulated and coupled under healthy conditions, and a break-down in the coupling has been associated with aging and disease. We present an approach that incorporates physiological modeling within a switching linear dynamical systems (SLDS) framework to assess the various functional components of the autonomic regulation through transfer function analysis of nonstationary multivariate time series of vital signs. We validate our proposed SLDS-based transfer function analysis technique in automatically capturing 1) changes in baroreflex gain due to postural changes in a tilt-table study including ten subjects, and 2) the effect of aging on the autonomic control using HR/RESP recordings from 40 healthy adults. Next, using HR/BP time series of more than 450 adult ICU patients, we show that our technique can be used to reveal coupling changes associated with severe sepsis (AUC = 0.74, sensitivity = 0.74, specificity = 0.60). Our findings indicate that reduced HR/BP coupling is significantly associated with severe sepsis even after adjusting for clinical interventions (P  0.001). These results demonstrate the utility of our approach in phenotyping complex vital-sign dynamics, and in providing mechanistic hypotheses in terms of break-down of autoregulatory systems under healthy and disease conditions.

Neonatal hemodynamics: monitoring, data acquisition and analysis

Monitoring of cardiovascular function is critical to both clinical care and research as the use of sophisticated monitoring systems enable us to obtain accurate, reliable and real-time information on developmental hemodynamics in health and disease. Novel approaches to comprehensive hemodynamic monitoring and data acquisition will undoubtedly aid in developing a better understanding of developmental cardiovascular physiology in neonates. In addition, development and use of state-of-the-art, comprehensive hemodynamic monitoring systems enable the recognition of signs of cardiovascular compromise in its early stages, and provide information on the hemodynamic response to treatment in critically ill patients.

Autonomic-cardiorespiratory regulation: a physiology-based mathematical model

This paper presents a novel physiology-based mathematical model of autonomic-cardiorespiratory regulation described by a set of three nonlinear, coupled differential equations. We improved our previously proposed autonomic-cardiac regulation model by considering neuromechanical and mechanical coupling of cardiovascular and respiration systems including lung stretch-receptor reflex and venous return variation. We also introduced a differential equation describing respiration rate regulation which mainly originates in the medullary respiratory center. The results of simulation experiments suggest that the venous return variation generates a higher perturbation on heart rate and blood pressure than lung stretch-receptor reflex. The proposed model is also powerful in determining and removing direct respiratory impacts on parasympathetic activation tone to accurately extract parasympathetic activity caused by emotional states and environmental conditions.

A practical approach to parameter estimation applied to model predicting heart rate regulation

Mathematical models have long been used for prediction of dynamics in biological systems. Recently, several efforts have been made to render these models patient specific. One way to do so is to employ techniques to estimate parameters that enable model based prediction of observed quantities. Knowledge of variation in parameters within and between groups of subjects have potential to provide insight into biological function. Often it is not possible to estimate all parameters in a given model, in particular if the model is complex and the data is sparse. However, it may be possible to estimate a subset of model parameters reducing the complexity of the problem. In this study, we compare three methods that allow identification of parameter subsets that can be estimated given a model and a set of data. These methods will be used to estimate patient specific parameters in a model predicting baroreceptor feedback regulation of heart rate during head-up tilt. The three methods include: structured analysis of the correlation matrix, analysis via singular value decomposition followed by QR factorization, and identification of the subspace closest to the one spanned by eigenvectors of the model Hessian. Results showed that all three methods facilitate identification of a parameter subset. The "best" subset was obtained using the structured correlation method, though this method was also the most computationally intensive. Subsets obtained using the other two methods were easier to compute, but analysis revealed that the final subsets contained correlated parameters. In conclusion, to avoid lengthy computations, these three methods may be combined for efficient identification of parameter subsets.

Arterial blood pressure regulation following aorta clamping and declamping during surgery

In this paper, we propose the use of black box models for the system identification of the cardiopulmonary baroreflex control of arterial resistance and of ventricular contractility and of arterial baroreflex control of heart rate (HR) from invasive, continuous measurements of arterial blood pressure (ABP) and central venous pressure (CVP), and non invasive, continuous recordings of ECG and respiration. Two crucial phases of the abdominal aortic aneurism (AAA) repair were investigated: the clamping and declamping of aorta. The objective of the present work is to evaluate and to test the ability to monitor baroreflex responses to clamping and declamping maneuvers preceding and following aneurism removal.

Methods derived from nonlinear dynamics for analysing heart rate variability

Methods from nonlinear dynamics (NLD) have shown new insights into heart rate (HR) variability changes under various physiological and pathological conditions, providing additional prognostic information and complementing traditional time- and frequency-domain analyses. In this review, some of the most prominent indices of nonlinear and fractal dynamics are summarized and their algorithmic implementations and applications in clinical trials are discussed. Several of those indices have been proven to be of diagnostic relevance or have contributed to risk stratification. In particular, techniques based on mono- and multifractal analyses and symbolic dynamics have been successfully applied to clinical studies. Further advances in HR variability analysis are expected through multidimensional and multivariate assessments. Today, the question is no longer about whether or not methods from NLD should be applied; however, it is relevant to ask which of the methods should be selected and under which basic and standardized conditions should they be applied.