Linear and Nonlinear Analysis of the Carotid Sinus Baroreflex Dynamic Characteristics

Effects of nitric oxide inhalation on pulmonary arterial impedance: differences between normal and pulmonary hypertension male rats

Nitric oxide (NO) inhalation improves pulmonary hemodynamics in participants with pulmonary arterial hypertension (PAH). Although it can reduce pulmonary vascular resistance (PVR) in PAH, its impact on the dynamic mechanics of pulmonary arteries and its potential difference between control and participants with PAH remain unclear. PA impedance provides a comprehensive description of PA mechanics. With an arterial model, PA impedance can be parameterized into peripheral pulmonary resistance (Rp), arterial compliance (Cp), characteristic impedance of the proximal arteries (Zc), and transmission time from the main PA to the reflection site. This study investigated the effects of inhaled NO on PA impedance and its associated parameters in control and monocrotaline-induced pulmonary arterial hypertension (MCT-PAH) male rats (6/group). Measurements were obtained at baseline and during NO inhalation at 40 and 80 ppm. In both groups, NO inhalation decreased PVR and increased the left atrial pressure. Notably, its impact on PA impedance was frequency dependent, as revealed by reduced PA impedance modulus in the low-frequency range below 10 Hz, with little effect on the high-frequency range. Furthermore, NO inhalation attenuated Rp, increased Cp, and prolonged transmission time without affecting Zc. It reduced Rp more pronouncedly in MCT-PAH rats, whereas it increased Cp and delayed transmission time more effectively in control rats. In conclusion, the therapeutic effects of inhaled NO on PA impedance were frequency dependent and may differ between the control and MCT-PAH groups, suggesting that the effect on the mechanics differs depending on the pathological state.NEW & NOTEWORTHY Nitric oxide inhalation decreased pulmonary arterial impedance in the low-frequency range (<10 Hz) with little impact on the high-frequency range. It reduced peripheral pulmonary resistance more pronouncedly in pulmonary hypertension rats, whereas it increased arterial compliance and transmission time in control rats. Its effect on the mechanics of the pulmonary arteries may differ depending on the pathological status.

Short-term dynamic characteristics of diuresis during exogenous pressure perturbations with and without arterial baroreflex control

Although body fluid volume control by the kidneys may be classified as a long-term arterial pressure (AP) control system, it does not necessarily follow that the urine flow (UF) response to changes in AP is slow. We quantified the dynamic characteristics of the UF response to short-term AP changes by changing mean AP between 60 mmHg and 100 mmHg every 10 s according to a binary white noise sequence in anesthetized rats (n = 8 animals). In a baro-on trial (the carotid sinus baroreflex was enabled), the UF response represented the combined synergistic effects of pressure diuresis (PD) and neurally mediated antidiuresis (NMA). In a baro-fix trial (the carotid sinus pressure was fixed at 100 mmHg), the UF response mainly reflected the effect of PD. The UF step response was quantified using the sum of two exponential decay functions. The fast and slow components had time constants of 6.5 ± 3.6 s and 102 ± 85 s (means ± SD), respectively, in the baro-on trial. Although the gain of the fast component did not differ between the two trials (0.49 ± 0.21 vs. 0.66 ± 0.22 µL·min-1·kg-1·mmHg-1), the gain of the slow component was greater in the baro-on than in the baro-fix trial (0.51 ± 0.14 vs. 0.09 ± 0.39 µL·min-1·kg-1·mmHg-1, P = 0.023). The magnitude of NMA relative to PD was calculated to be 32.2 ± 29.8%. In conclusion, NMA contributed to the slow component, and its magnitude was approximately one-third of that of the effect of PD.NEW & NOTEWORTHY We quantified short-term dynamic characteristics of the urine flow (UF) response to arterial pressure (AP) changes using white noise analysis. The UF step response approximated the sum of two exponential decay functions with time constants of ∼6.5 s and 102 s. The neurally mediated antidiuretic (NMA) effect contributed to the slow component of the UF step response, with the magnitude of approximately one-third of that of the pressure diuresis (PD) effect.

Input-size dependence of the baroreflex neural arc transfer characteristics during Gaussian white noise inputs

Although Gaussian white noise (GWN) inputs offer a theoretical framework for identifying higher-order nonlinearity, an actual application to the data of the neural arc of the carotid sinus baroreflex did not succeed in fully predicting the well-known sigmoidal nonlinearity. In the present study, we assumed that the neural arc can be approximated by a cascade of a linear dynamic (LD) component and a nonlinear static (NS) component. We analyzed the data obtained using GWN inputs with a mean of 120 mmHg and standard deviations (SDs) of 10, 20, and 30 mmHg for 15 min each in anesthetized rats (n = 7). We first estimated the linear transfer function from carotid sinus pressure to sympathetic nerve activity (SNA) and then plotted the measured SNA against the linearly predicted SNA. The predicted and measured data pairs exhibited an inverse sigmoidal distribution when grouped into 10 bins based on the size of the linearly predicted SNA. The sigmoidal nonlinearity estimated via the LD-NS model showed a midpoint pressure (104.1 ± 4.4 mmHg for SD of 30 mmHg) lower than that estimated by a conventional stepwise input (135.8 ± 3.9 mmHg, P < 0.001). This suggests that the NS component is more likely to reflect the nonlinearity observed during pulsatile inputs that are physiological to baroreceptors. Furthermore, the LD-NS model yielded higher R2 values compared with the linear model and the previously suggested second-order Uryson model in the testing dataset.NEW & NOTEWORTHY We examined the input-size dependence of the baroreflex neural arc transfer characteristics during Gaussian white noise inputs. A linear dynamic-static nonlinear model yielded higher R2 values compared with a linear model and captured the well-known sigmoidal nonlinearity of the neural arc, indicating that the nonlinear dynamics contributed to determining sympathetic nerve activity. Ignoring such nonlinear dynamics might reduce our ability to explain underlying physiology and significantly limit the interpretation of experimental data.

Dynamic accentuated antagonism of heart rate control during different levels of vagal nerve stimulation intensity in rats

Accentuated antagonism refers to a phenomenon in which the vagal effect on heart rate (HR) is augmented by background sympathetic tone. The dynamic aspect of accentuated antagonism remains to be elucidated during different levels of vagal nerve stimulation (VNS) intensity. We performed VNS on anesthetized rats (n = 8) according to a binary white noise signal with a switching interval of 500 ms at three different stimulation rates (low-intensity: 0-10 Hz, moderate-intensity: 0-20 Hz, and high-intensity: 0-40 Hz). The transfer function from VNS to HR was estimated with and without concomitant tonic sympathetic nerve stimulation (SNS) at 5 Hz. The asymptotic low-frequency (LF) gain (in beats/min/Hz) of the transfer function increased with SNS regardless of the VNS rate [low-intensity: 3.93 ± 0.70 vs. 5.82 ± 0.65 (P = 0.021), moderate-intensity: 3.87 ± 0.62 vs. 5.36 ± 0.53 (P = 0.018), high-intensity: 4.77 ± 0.85 vs. 7.39 ± 1.36 (P = 0.011)]. Moreover, SNS slightly increased the ratio of high-frequency (HF) gain to the LF gain. These effects of SNS were canceled by the pretreatment of ivabradine, an inhibitor of hyperpolarization-activated cyclic nucleotide-gated channels, in another group of rats (n = 6). Although background sympathetic tone antagonizes the vagal effect on mean HR, it enables finer HR control by increasing the dynamic gain of the vagal HR transfer function regardless of VNS intensity. When interpreting the HF component of HR variability, the augmenting effect from background sympathetic tone needs to be considered.

Linear and nonlinear identification of the carotid sinus baroreflex in the very low-frequency range

Since the arterial baroreflex system is classified as an immediate control system, the focus has been on analyzing its dynamic characteristics in the frequency range between 0.01 and 1 Hz. Although the dynamic characteristics in the frequency range below 0.01 Hz are not expected to be large, actual experimental data are scant. The aim was to identify the dynamic characteristics of the carotid sinus baroreflex in the frequency range down to 0.001 Hz. The carotid sinus baroreceptor regions were isolated from the systemic circulation, and carotid sinus pressure (CSP) was changed every 10 s according to Gaussian white noise with a mean of 120 mmHg and standard deviation of 20 mmHg for 90 min in anesthetized Wistar‐Kyoto rats (n = 8). The dynamic gain of the linear transfer function relating CSP to arterial pressure (AP) at 0.001 Hz tended to be greater than that at 0.01 Hz (1.060 ± 0.197 vs. 0.625 ± 0.067, p = 0.080), suggesting that baroreflex control was largely maintained at 0.001 Hz. Regarding nonlinear analysis, a second‐order Uryson model predicted AP with a higher R² value (0.645 ± 0.053) than a linear model (R² = 0.543 ± 0.057, p = 0.025) or a second‐order Volterra model (R² = 0.589 ± 0.055, p = 0.045) in testing data. These pieces of information may be used to create baroreflex models that can add a component of autonomic control to a cardiovascular digital twin for predicting acute hemodynamic responses to treatments and tailoring individual treatment strategies.

Determinants of orthostatic cerebral oxygenation assessed using near-infrared spectroscopy

BACKGROUND

To understand the relationship between blood pressure changes during standing up and clinical outcome, cerebral oxygenation needs to be measured, which may be performed using near-infrared spectroscopy (NIRS). However, the role of potential determinants of NIRS-derived orthostatic cerebral oxygenation, i.e., age, sex, type of postural change (i.e., standing up from sitting versus supine position), blood pressure (BP) and baroreflex sensitivity (BRS) is still unknown and needed to better interpret findings from studies using orthostatic NIRS measurements.

METHODS

34 younger (median age 25 years, inter quartile range (IQR) 22-45) and 31 older adults (median age 77 years, IQR 72-81) underwent BP, BRS and NIRS measurements during standing up from sitting and supine position. Linear regression models were used to assess the potential determinant role of age, sex, type of postural change, BP and BRS in orthostatic cerebral oxygenation drop and recovery. Orthostatic cerebral oxygenation test-retest reliability was assessed using intra class correlations.

RESULTS

Younger age, male sex and standing up from supine compared to sitting position were positively associated with cerebral oxygenation drop; older age and standing up from sitting compared to supine position were associated with higher cerebral oxygenation recovery. Test-retest reliability was highest (ICC > 0.83) during standing up from supine position.

CONCLUSION

Based on the findings of this study, age, sex and type of postural change are significant determinants of NIRS-derived orthostatic cerebral oxygenation and should be taken into account in the interpretation of NIRS measurements. In the design of new studies, standing up from supine position is preferable (higher reliability) over standing up from sitting position.

Ivabradine increases the high frequency gain ratio in the vagal heart rate transfer function via an interaction with muscarinic potassium channels

Muscarinic potassium channels (IK,ACh ) are thought to contribute to the high frequency (HF) dynamic heart rate (HR) response to vagal nerve stimulation (VNS) because they act faster than the pathway mediated by hyperpolarization-activated cyclic nucleotide-gated (HCN) channels. However, the interactions between the two pathways have not yet been fully elucidated. We previously demonstrated that HCN channel blockade by ivabradine (IVA) increased the HF gain ratio of the transfer function from VNS to HR. To test the hypothesis that IVA increases the HF gain ratio via an interaction with IK,ACh , we examined the dynamic HR response to VNS under conditions of control (CNT), IK,ACh blockade by tertiapin-Q (TQ, 50 nM/kg), and TQ plus IVA (2 mg/kg) (TQ + IVA) in anesthetized rats (n = 8). In each condition, the right vagal nerve was stimulated for 10 min with binary white noise signals between 0-10, 0-20, and 0-40 Hz. On multiple regression analysis, the HF gain ratio positively correlated with the VNS rate with a coefficient of 1.691 ± 0.151 (×0.01) (p < 0.001). TQ had a negative effect on the HF gain ratio with a coefficient of -1.170 ± 0.214 (×0.01) (p < 0.001). IVA did not significantly increase the HF gain ratio in the presence of TQ. The HF gain ratio remained low under the TQ + IVA condition compared to controls. These results affirm that the IVA-induced increase in the HF gain ratio is dependent on the untethering of the hyperpolarizing effect of IK,ACh .

Closed-Loop Identification of Baroreflex Properties in the Frequency Domain

The arterial baroreflex system plays a key role in maintaining the homeostasis of arterial pressure (AP). Changes in AP affect autonomic nervous activities through the baroreflex neural arc, whereas changes in the autonomic nervous activities, in turn, alter AP through the baroreflex peripheral arc. This closed-loop negative feedback operation makes it difficult to identify open-loop dynamic characteristics of the neural and peripheral arcs. Regarding sympathetic AP controls, we examined the applicability of a nonparametric frequency-domain closed-loop identification method to the carotid sinus baroreflex system in anesthetized rabbits. This article compares the results of an open-loop analysis applied to open-loop data, an open-loop analysis erroneously applied to closed-loop data, and a closed-loop analysis applied to closed-loop data. To facilitate the understanding of the analytical method, sample data files and sample analytical codes were provided. In the closed-loop identification, properties of the unknown central noise that modulated the sympathetic nerve activity and the unknown peripheral noise that fluctuated AP affected the accuracy of the estimation results. A priori knowledge about the open-loop dynamic characteristics of the arterial baroreflex system may be used to advance the assessment of baroreflex function under closed-loop conditions in the future.

Ivabradine augments high-frequency dynamic gain of the heart rate response to low- and moderate-intensity vagal nerve stimulation under β-blockade

Our previous study indicated that intravenously administered ivabradine (IVA) augmented the dynamic heart rate (HR) response to moderate-intensity vagal nerve stimulation (VNS). Considering an accentuated antagonism, the results were somewhat paradoxical; i.e., the accentuated antagonism indicates that an activation of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels via the accumulation of intracellular cyclic adenosine monophosphate (cAMP) augments the HR response to VNS, whereas the inhibition of HCN channels by IVA also augmented the HR response to VNS. To remove the possible influence from the accentuated antagonism, we examined the effects of IVA on the dynamic vagal control of HR under β-blockade. In anesthetized rats (n = 7), the right vagal nerve was stimulated for 10 min according to binary white noise signals between 0 and 10 Hz (V_0-10), between 0 and 20 Hz (V_0-20), and between 0 and 40 Hz (V_0-40). The transfer function from VNS to HR was estimated. Under β-blockade (propranolol, 2 mg/kg iv), IVA (2 mg/kg iv) did not augment the asymptotic low-frequency gain but increased the asymptotic high-frequency gain in V_0-10 (0.53 ± 0.10 vs. 1.74 ± 0.40 beats/min/Hz, P < 0.01) and V_0-20 (0.79 ± 0.14 vs. 2.06 ± 0.47 beats/min/Hz, P < 0.001). These changes, which were observed under a minimal influence from sympathetic background tone, may reflect an increased contribution of the acetylcholine-sensitive potassium channel (I_K,ACh) pathway after IVA, because the HR control via the I_K,ACh pathway is faster and acts in the frequency range higher than the cAMP-mediated pathway.NEW & NOTEWORTHY Since ivabradine (IVA) inhibits hyperpolarization-activated cyclic nucleotide-gated channels, interactions among the sympathetic effect, vagal effect, and IVA can occur in the control of heart rate (HR). To remove the sympathetic effect, we estimated the transfer function from vagal nerve stimulation to HR under β-blockade in anesthetized rats. IVA augmented the high-frequency dynamic gain during low- and moderate-intensity vagal nerve stimulation. Untethering the hyperpolarizing effect of acetylcholine-sensitive potassium channels after IVA may be a possible underlying mechanism.