- Original Paper
- Published:
Open-loop dynamic and static characteristics of the carotid sinus baroreflex in rats with chronic heart failure after myocardial infarction
The Journal of Physiological Sciences volume 60, pages 283–298 (2010)
Abstract
We estimated open-loop dynamic characteristics of the carotid sinus baroreflex in normal control rats and chronic heart failure (CHF) rats after myocardial infarction. First, the neural arc transfer function from carotid sinus pressure to splanchnic sympathetic nerve activity (SNA) and its corresponding step response were examined. Although the steady-state response was attenuated in CHF, the negative peak response and the time to peak did not change significantly, suggesting preserved neural arc dynamic characteristics. Next, the peripheral arc transfer function from SNA to arterial pressure (AP) and its corresponding step response were examined. The steady-state response and the initial slope were reduced in CHF, suggesting impaired end-organ responses. In a simulation study based on the dynamic and static characteristics, the percent recovery of AP was reduced progressively as the size of disturbance increased in CHF, suggesting that a reserve for AP buffering is lost in CHF despite relatively maintained baseline AP.
Introduction
The arterial baroreflex is an important negative feedback system that stabilizes systemic arterial pressure (AP) against exogenous disturbances. The rapidness of AP regulation may be best described by the system dynamic characteristics. With respect to the sympathetic arterial baroreflex system, previous studies in rabbits [1] and rats [2] have indicated that the transfer function of the baroreflex neural arc from baroreceptor pressure input to efferent sympathetic nerve activity (SNA) exhibits “derivative” characteristics, which means that the dynamic gain of the SNA response to pressure perturbation becomes greater as the modulation frequency increases. On the other hand, the transfer function of the baroreflex peripheral arc from SNA to AP shows “low-pass” characteristics, which means that the dynamic gain of the AP response to SNA variation becomes smaller as the modulation frequency increases. In short, the neural arc provides an accelerating mechanism of the dynamic AP response in the arterial baroreflex system [1].
Although a number of studies have indicated that the baroreflex function is depressed in heart failure [3–6], dynamic characteristics of the arterial baroreflex in heart failure have not been fully described. The aim of the present study was to identify the open-loop dynamic characteristics of the carotid sinus baroreflex in a rat model of chronic heart failure (CHF) after myocardial infarction. To obtain a total picture of the AP regulation in CHF, we also estimated the open-loop static characteristics of the carotid sinus baroreflex.
Materials and methods
Animals were cared for in strict accordance with the Guiding Principles for the Care and Use of Animals in the Field of Physiological Sciences, which has been approved by the Physiological Society of Japan. All experimental protocols were reviewed and approved by the Animal Subjects Committee at the National Cerebral and Cardiovascular Center.
Myocardial infarction
Coronary artery ligation was performed under sterile conditions in 8-week-old male Sprague-Dawley rats according to a previously established procedure [7]. After inducing anesthesia by halothane inhalation, the rat was intubated and mechanically ventilated. The left chest was opened at the fourth intercostal space and the left coronary artery was ligated with a 5-0 polypropylene suture (PROLENE, Ethicon, GA, USA). An electrocardiogram was monitored for 1 h after the coronary ligation, and the heart was defibrillated as necessary by manual prodding. Thereafter air was evacuated from the thoracic cavity and the incision was closed. The rats were allowed to recover from anesthesia, and were fed ad libitum with a standard laboratory chow and given free access to water. Although we initially planned to wait until the rats showed an objective sign of advanced heart failure such as facial edema or labored breathing, the rats with such severe heart failure were too weak to survive the acute baroreflex study described in the following section. Accordingly, the rats that survived for 100–200 days (156 ± 18 days) after myocardial infarction were used without solid criteria for starting the acute baroreflex study. Instead, the rats were regarded as experiencing CHF when the central venous pressure was higher than 2.5 mmHg or the biventricular weight was greater than 2.5 g/kg body weight at the time of the acute baroreflex study.
Acute baroreflex study
Experiments were performed in normal control (n = 12) and CHF (n = 7) rats. Sham operation was not performed in the control rats. Among the control rats, eight rats were matched with the CHF rats based on body weight, but they were younger in age due to the retardation of growth in the CHF rats. The remaining four rats were age-matched with the CHF rats, but they were heavier in body weight. Because none of the parameters of the baroreflex dynamic and static characteristics differed statistically between the two subgroups of the control rats, we report pooled data obtained from the 12 control rats.
Each rat was anesthetized with an intraperitoneal injection (2 ml/kg) of a mixture of urethane (250 mg/ml) and α-chloralose (40 mg/ml), and mechanically ventilated with oxygen-enriched room air. A venous catheter was inserted into the right femoral vein, and 20-fold diluted anesthetic mixture was administered continuously (2–3 ml kg−1 h−1). An arterial catheter was inserted into the right femoral artery to measure AP. Heart rate (HR) was obtained from AP through a cardiotachometer. Another venous catheter was inserted into the left femoral vein and advanced into the inferior vena cava to measure central venous pressure and to supply Ringer solution (6 ml kg−1 h−1).
A postganglionic branch from the splanchnic sympathetic nerve was exposed through a left flank incision, and a pair of stainless steel wire electrodes (Bioflex wire AS633, Cooner Wire, CA, USA) was attached to record SNA. The nerve and electrodes were covered with silicone glue (Kwik-Sil, World Precision Instruments, FL, USA) for insulation and fixation. To quantify the nerve activity, the preamplified nerve signal was band-pass filtered at 150–1,000 Hz, and then full-wave rectified and low-pass filtered with a cut-off frequency of 30 Hz. Pancuronium bromide (0.4 mg kg−1 h−1) was administered to prevent muscular activity from contaminating the SNA recording. At the end of the experiment, we confirmed the disappearance of SNA in response to an intravenous bolus injection of a ganglionic blocker, hexamethonium bromide (60 mg kg−1), and recorded the noise level.
Bilateral vagal and aortic depressor nerves were sectioned at the neck to avoid reflexes from the cardiopulmonary region and aortic arch. The carotid sinus regions were isolated from the systemic circulation bilaterally according to previously reported procedures [8, 9]. Briefly, a 7-0 polypropylene suture with a fine needle (PROLENE, Ethicon, GA, USA) was passed through the tissue between the external and internal carotid arteries, and the external carotid artery was ligated close to the carotid bifurcation. The internal carotid artery was embolized using two to three steel balls (0.8 mm in diameter, Tsubaki Nakashima, Nara, Japan) injected from the common carotid artery. Under these conditions, the brain stem area was perfused by patent bilateral vertebral arteries. The isolated carotid sinuses were filled with warmed Ringer solution through catheters inserted into the common carotid arteries. Carotid sinus pressure (CSP) was controlled using a servo-controlled piston pump. Heparin sodium (100 U kg−1) was given intravenously to prevent blood coagulation. Body temperature was maintained at approximately 38°C with a heating pad.
Protocols
After the above surgical procedures were completed, reflex responses in SNA, AP, and HR to CSP input were monitored for more than 30 min. The rat was excluded from further study and analysis in the event that the reflex responses became smaller within this period. Possible causes for the deterioration in the reflex responses include surgical damage to the carotid sinus nerves and brain ischemia due to the bilateral carotid occlusion.
To estimate the dynamic input-output relationship of the carotid sinus baroreflex, CSP was perturbed for 20 min using a Gaussian white noise (GWN) signal with the mean of 120 mmHg and standard deviation of 20 mmHg. The selection of the mean pressure and amplitude of GWN does not significantly affect the estimation of the system dynamic characteristics except for a factor of proportionality (see Appendix 1 for details). The switching interval of GWN was set at 500 ms. The input power spectral density was relatively constant up to 1 Hz, which covered the upper frequency range of interest with respect to the sympathetic arterial baroreflex in rats [2].
To estimate the static input-output relationship of the carotid sinus baroreflex, CSP was decreased to 60 mmHg for 4–6 min, and then increased stepwise from 60 to 180 mmHg in increments of 20 mmHg every minute [10].
Data analysis
Data were sampled at 200 Hz using a 16-bit analog-to-digital converter and stored on a dedicated laboratory computer system. In each rat, the noise level of SNA recorded after the administration of hexamethonium bromide was treated as zero. Because the absolute voltage of SNA varied among animals depending on the recording conditions, SNA averaged during the last 10 s at CSP of 60 mmHg in the stepwise input protocol was defined as 100%. The same normalization factor was used for the analysis of the baroreflex dynamic characteristics.
Dynamic characteristics of the baroreflex neural arc, peripheral arc, total baroreflex, and HR control were estimated by an open-loop transfer function analysis as follows [11]. Data were analyzed from 90 s after the initiation of the GWN input. The input-output data pairs were resampled at 10 Hz and segmented into 50% overlapping bins of 1,024 points each. For each segment, a linear trend was removed, and a Hanning window was applied. Fast Fourier transform was performed to obtain the frequency spectra of the input and output signals. The ensemble averages of the input power spectral density [S XX (f)], output power spectral density [S YY (f)], and cross spectral density between the input and output signals [S YX (f)] were calculated over 12 segments, where f denotes frequency. Finally, the transfer function [H(f)] from input to output was estimated as:
The transfer function is a complex-valued function that can be expressed by the modulus and phase at each frequency. In the present study, we refer to the modulus of the transfer function as the dynamic gain. To quantify the linear dependence between the input and output signals, a magnitude squared coherence function [Coh(f)] was calculated as:
The coherence function is a real-valued function ranging from zero to unity. When the output signal is perfectly explained by the linear dynamics with the input signal, the coherence value becomes unity. When the output signal is totally independent of the input signal, the coherence value becomes zero.
To facilitate understanding of the transfer function, the step response was also calculated as follows. The system impulse response was derived from the inverse Fourier transform of H(f). The step response was then obtained from the time integral of the impulse response.
To quantify the open-loop static characteristics of the carotid sinus baroreflex, mean SNA, AP, and HR were obtained during the last 10 s at each CSP level of the stepwise input protocol. In each rat, data from two consecutive step cycles were averaged. The static characteristics of the baroreflex neural arc (the CSP-SNA relationship), the total baroreflex (the CSP-AP relationship), and the HR control (the CSP-HR relationship) were described by fitting four-parameter logistic functions to the input-output data as follows [10, 12]:
where x and y denote the input (CSP) and output (SNA, AP, or HR), respectively; P 1 is the response range of output; P 2 is the slope coefficient; P 3 is the midpoint pressure of input; and P 4 is the minimum value of output. For convenience, the maximum slope or the maximum gain of the logistic function is reported by a positive value as P 1 P 2/4.
The static characteristics of the baroreflex peripheral arc (the SNA-AP relationship) were quantified by a linear regression analysis as follows [10, 13]:
where a and b represent the slope and intercept, respectively, of the regression line.
Statistical analysis
All data are presented in mean and SE values. Unpaired t tests were used to compare the parameters of the baroreflex dynamic and static characteristics between the control and CHF groups [13]. To compare the transfer functions between the two groups, we arbitrarily selected the dynamic gain values at 0.01, 0.1, and 1 Hz (G 0.01, G 0.1, and G 1). For the step response relating to the neural arc transfer function, the negative peak response (S peak), the time to the negative peak (T peak), the value of the step response at 10 s (S 10), the steady-state response at 50 s (S 50), and the ratio of the peak response to the steady-state response (S peak/S 50) were calculated. For the step response relating to the peripheral arc, total baroreflex, or HR control, an initial slope of the response was calculated (see Appendix 2) in addition to S50. Differences were considered significant when P < 0.05.
Simulation study
Once the open-loop dynamic and static characteristics of a system are both identified, closed-loop system responses can be simulated [14–16]. To compare closed-loop behavior of the carotid sinus baroreflex between control and CHF conditions, step inputs ranging from −10 to −60 mmHg were applied as exogenous disturbances, and resulting AP responses were simulated. Percent recovery was calculated as the magnitude of steady-state AP recovery relative to the size of step disturbance. The initial slope for the recovery response was also calculated (see Appendix 2).
Results
Postmortem examination confirmed that the left ventricular free wall was reduced to a membrane-like scar in the CHF group. Biventricular weights, both absolute and relative to body weight, were significantly greater in the CHF than in the control group (Table 1). The central venous pressure was significantly higher and the baseline AP and HR were significantly lower in the CHF group. Although the duration after myocardial infarction ranged from 100 to 200 days in the CHF group, there was no significant correlation between the duration and biventricular weight (y = 0.0012x + 2.53, r 2 = 0.035, P = 0.72, x duration in days, y biventricular weight in grams per kilogram) or between the duration and central venous pressure (y = −0.013x + 7.5, r 2 = 0.07, P = 0.58, x duration in days, y central venous pressure in millmeters of mercury).
Typical experimental recordings obtained from a control rat are shown in Fig. 1a. m-SNA indicates a 2-s moving averaged signal of SNA. In the dynamic input protocol, CSP was changed dynamically according to a GWN signal. m-SNA varied dynamically in response to the CSP perturbation. Although AP changed dynamically, the AP variation seemed more sluggish than the SNA variation. Changes in HR were less obvious from the time series data. In the static input protocol, a stepwise increase in CSP decreased m-SNA, AP, and HR. The noise level of SNA obtained after the intravenous administration of hexamethonium bromide was set at zero. Because of the normalization procedure, m-SNA during the last 10 s at CSP of 60 mmHg approximated 100%. Figure 1b represents CSP and SNA signals sampled at 200 Hz during the dynamic input protocol. High CSP inputs suppressed SNA to a noise level.
Typical experimental recordings obtained from a CHF rat are shown in Fig. 2a. In the dynamic input protocol, CSP was changed dynamically according to a GWN signal. Although m-SNA varied dynamically in response to the CSP perturbation, changes in AP and HR were not obvious from the time series data. In the static input protocol, a stepwise increase in CSP decreased m-SNA, AP, and HR. The magnitudes of the responses in m-SNA, AP, and HR seem smaller than those in the control rat. Figure 2b represents CSP and SNA signals sampled at 200 Hz during the dynamic input protocol. High CSP inputs suppressed SNA to a noise level.
Dynamic characteristics of the carotid sinus baroreflex
Figure 3 summarizes the open-loop dynamic characteristics of the carotid sinus baroreflex averaged for the control and CHF rats. The neural arc transfer function from CSP to SNA showed derivative characteristics (Fig. 3a). G 0.01 tended to be lower in the CHF rats, whereas G 0.1 and G 1 did not differ significantly between the two groups (Table 2). The phase plot showed an out-of-phase relationship in the frequency range from 0.01 to 1 Hz in both groups. The coherence plot with values less than unity suggests that SNA contained an unknown noise signal unrelated to the baroreflex and/or a nonlinear system response to the CSP input. The bottom panel of Fig. 3a represents the step responses of SNA for a unit increase in CSP. Although S 10 and S 50 were significantly attenuated in the CHF rats, S peak, T peak, and S peak/S 50 did not differ significantly between the two groups (Table 2).
The peripheral arc transfer function from SNA to AP showed low-pass characteristics (Fig. 3b). G 0.01 and G 0.1 were significantly smaller in the CHF rats, whereas G 1 did not differ significantly between the two groups (Table 2). The phase was close to zero radians at the lowest frequency and delayed with an increase in frequency in both groups. The coherence plot with values less than unity indicates that a part of AP variation was not explained by the linear dynamics with the SNA variation. In the step response of AP for a unit increase in SNA (Fig. 3b, bottom), S 50 was significantly smaller and the initial slope (the dotted line) was significantly gentler in the CHF than in the control rats (Table 2).
The total baroreflex transfer function from CSP to AP showed low-pass characteristics (Fig. 3c). G 0.01 and G 0.1 were significantly smaller in the CHF than in the control rats (Table 2), but G 1 did not differ significantly. The phase was close to −π radians at the lowest frequency in both groups, reflecting the negative feedback operation attained by the total baroreflex. The phase delayed with an increase in frequency. The coherence values seem lower than those in the peripheral arc transfer function. In the step response of AP for a unit increase in CSP (Fig. 3c, bottom), both S50 and the initial slope (the dotted line) were significantly attenuated in the CHF in comparison with the control rats (Table 2).
The transfer function from CSP to HR also displayed low-pass characteristics (Fig. 3d). G 0.01 and G 0.1 were significantly smaller in the CHF than in the control rats (Table 2). We did not compare G 1 because the coherence values close to zero and the irregular change in the phase shift above 0.8 Hz suggested poor reliability of the estimated transfer function. In the step response of HR for a unit increase in CSP (Fig. 3d, bottom), both S 50 and the initial slope (the dotted line) were significantly attenuated in the CHF in comparison with the control rats (Table 2).
Static characteristics of the carotid sinus baroreflex
Figure 4 summarizes the open-loop static characteristics of the carotid sinus baroreflex obtained from the control and CHF rats. The baroreflex neural arc showed a decreasing SNA response with an increase in CSP (Fig. 4a). The response range of SNA was significantly narrower in the CHF rats (Table 3). Consequently, the minimum SNA was significantly higher in the CHF rats. The midpoint pressure on the CSP axis was significantly lower in the CHF rats. Despite the significant attenuation in the response range of SNA, the maximum slope of the neural arc was not reduced in the CHF rats compared with the control rats.
The peripheral arc from SNA to AP approximated a straight line in both the control and CHF rats (Fig. 4b). The slope of the regression line was significantly less steep in the CHF rats, whereas the AP intercept did not differ significantly between the two groups (Table 3).
The total baroreflex function from CSP to AP approximated an inverse sigmoidal curve (Fig. 4c). The response range of AP was significantly narrower in the CHF rats (Table 3). The slope coefficient, the midpoint pressure on the CSP axis, and the minimum AP did not differ significantly between the two groups. The maximum gain was significantly smaller in the CHF rats compared with the control rats.
The static CSP-HR relationship also approximated an inverse sigmoidal curve (Fig. 4d). The response range of HR was significantly narrower, and the minimum HR was significantly lower in the CHF rats (Table 3). The midpoint pressure on the CSP axis was significantly lower in the CHF rats. The slope coefficient and the maximum slope did not differ significantly between the two groups.
The baroreflex equilibrium diagram is obtained by plotting the neural and peripheral arcs on a pressure-SNA plane (Fig. 4e). The ordinate is either CSP (for the neural arc) or AP (for the peripheral arc). The intersection between the neural and peripheral arcs gives the closed-loop operating point [17, 18]. The operating-point AP (the horizontal arrow) was significantly lower in the CHF rats, whereas the operating-point SNA (the vertical arrow) did not differ significantly between the two groups (Table 3). The total baroreflex gain at the operating point, calculated from the product of the tangential slope of the neural arc (ΔSNA/ΔCSP) and the slope of the peripheral arc (ΔAP/ΔSNA), did not differ significantly between the two groups (Table 3): ΔCSP, ΔSNA, and ΔAP indicate small changes in CSP, SNA, and AP, respectively, at the intersection point. Figure 4f depicts a putative baroreflex equilibrium diagram where the SNA axis is scaled so that the maximum value of SNA in the CHF group becomes two times higher than that in the control group (see Discussion).
Simulation study
A block diagram of the simulation study is shown in Fig. 5a [14–16]. We used impulse responses derived from the group-averaged neural and peripheral arc transfer functions (H N and H P) to calculate dynamic responses of the carotid sinus baroreflex. The steady-state gains of H N and H P were normalized to unity in order that absolute values of the steady-state gains could be determined by their corresponding static characteristics. The static characteristics of the neural and peripheral arcs were modeled as a logistic function and a regression line, respectively, using group-averaged parameter values (Table 3). A sinusoidal wave with an amplitude of 15 mmHg (peak-to-peak pressure of 30 mmHg) and a frequency of 5 Hz (corresponding to 300 beats/min) was added to the output from the peripheral arc to mimic pulsatile pressure. The AP signal was fed back into the neural arc to achieve a closed-loop simulation. After the AP signal reached steady state, step disturbances ranging from −10 to −60 mmHg were imposed.
Typical examples of the transient AP response to a step disturbance of −40 mmHg under control and CHF conditions are shown in Fig. 5b. In each panel, the pulsatile pressure is shown in gray, and the mean AP signal is shown as a solid bold line. The horizontal dashed lines represent the mean AP values immediately before and after the onset of the step disturbance. The AP signal decreased abruptly by 40 mmHg at time zero and recovered gradually thereafter. The upward arrow indicates the magnitude of steady-state AP recovery. The recovery was greater in the control than in the CHF simulation.
Figure 5c depicts the percent recovery of AP relative to the size of the step disturbance. The steady-state AP recovery was approximately 50% in the control and approximately 32% in the CHF simulation for a step disturbance of −10 mmHg, indicating that the feedback AP regulation in the CHF simulation was approximately 64% as effective as that in the control simulation. On the other hand, the AP recovery was approximately 35% in the control and approximately 14% in the CHF simulation for a step disturbance of −60 mmHg, indicating that the efficiency of feedback AP regulation in the CHF simulation reduced to only 40% of that in the control simulation.
Figure 5d shows the initial slope of the AP recovery. The initial slope increased as the size of step disturbance increased in the control simulation. Although the initial slope of the AP recovery in the CHF simulation was comparable to that in the control simulation for a step disturbance of −10 mmHg, it did not increase significantly with an increase in the size of step disturbance.
Discussion
The major findings of the present study are (1) the dynamic characteristics of the baroreflex neural arc were preserved in CHF, whereas those of the baroreflex peripheral arc were significantly depressed (Fig. 3), and (2) the total baroreflex gain at the closed-loop operating point seemed preserved in CHF, whereas the range of baroreflex operation was significantly narrowed (Fig. 4). Because of these modulations in the baroreflex characteristics, the AP regulation was less robust against exogenous disturbances in CHF (Fig. 5), which may partly explain the incidence of acute decompensation in stable CHF patients caused by noncompliance with salt and water restriction [19].
Dynamic characteristics of the carotid sinus baroreflex in CHF rats
The derivative characteristics of the neural arc are preserved in the CHF rats (Fig. 3a), being consistent with a previous study in heart failure rabbits [20]. Moreover, the present results indicate that the dynamic AP response to SNA is significantly depressed in CHF (Fig. 3b), suggesting impaired end-organ responses to SNA. The total baroreflex function in terms of the AP regulation was more sluggish in the CHF rats (Fig. 3c, bottom) despite the preserved neural arc derivative characteristics.
In the transfer function from CSP to HR, dynamic gain values were significantly attenuated in the CHF rats (Fig. 3d). Because the dynamic gain of the neural arc transfer function did not reduce significantly (G 0.1), the attenuation of the dynamic gain of the HR control may be attributable to the reduced HR response to SNA such as that related to the downregulation of β-adrenergic receptors [21].
Static characteristics of the carotid sinus baroreflex in CHF rats
In both the control and CHF rats, the input-output relationship of the peripheral arc approximated a straight line (Fig. 4b), and the sigmoidal nonlinearity is primarily attributed to the neural arc (Fig. 4a). The impairment of the total baroreflex in the CHF rats, characterized by a significant attenuation of the response range of AP and a reduction of the maximum gain (Fig. 4c, Table 3), is in essence similar to that obtained previously in canine models of heart failure [3, 5, 6].
Wang et al. [4, 5] reported poor end-organ responses and normal central control of renal SNA in the heart failure dogs. The attenuated slope of the peripheral arc in the CHF rats (Fig. 4b) may be consistent with the poor end-organ responses. Although the response range of SNA in the neural arc was narrowed in the CHF rats (Fig. 4a), the maximum slope of the neural arc did not change significantly, which should contribute to the maintained total baroreflex gain within a small range around the operating point (Fig. 4e, Table 3).
The response range of HR was significantly reduced in the CHF rats (Fig. 4d) in agreement with the depressed dynamic HR response to CSP (Fig. 3d). Although 24-h averaged HR was higher in the CHF rats in a previous study [7], the minimum HR was significantly lower in the CHF than in the control rats. Anesthesia and vagotomy might have affected the results, and further studies are needed to reconcile some of the discrepancies between the present and previous results.
Equilibrium diagram and simulation study
A baroreflex equilibrium diagram provides information on the closed-loop operating point (Fig. 4e) [10, 17, 18, 22–26]. Because absolute SNA is considered to be higher in CHF [7], if we can use the absolute value for SNA, the SNA axis of the equilibrium diagram should be scaled in a manner that renders the peripheral arc much shallower in the CHF than in the control rats (Fig. 4f). The operating-point AP was decreased only by 10 mmHg in the CHF rats and the total baroreflex gain at the operating point did not differ between the control and CHF rats. The baroreflex function may therefore seem preserved in the CHF rats when it is assessed only within a small range around the closed-loop operating point.
We used splanchnic SNA to construct a baroreflex equilibrium diagram. The analysis based on SNA of only one limb of the sympathetic nervous system can be an oversimplification in light of the well-known regional differences in SNAs [27–29]. One rationale for using splanchnic SNA is that the splanchnic region has been regarded as a major site for the blood flow redistribution [17, 30]. Furthermore, we have continued to observe insignificant differences in the steady-state responses to CSP between left and right cardiac SNAs [31], between cardiac and renal SNAs [32], and among cardiac, renal, and muscle SNAs [33] in anesthetized rabbits. It seems that the carotid sinus baroreflex brings about common activity in addition to regional activity among sympathetic nerves innervating several districts (see Appendix 3). Whether the similarity among SNAs to different districts will hold in other experimental conditions awaits future investigations.
The simulation study is useful and necessary for the integrated understanding of the estimated dynamic and static characteristics of the carotid sinus baroreflex (Fig. 5a). The feedback AP regulation becomes much weaker in CHF as the size of step disturbance increases (Fig. 5c, d), suggesting that a reserve for AP buffering is lost in CHF. Because not only the baroreceptor input but also loading conditions to the heart are changed during orthostatic tilt, the results of the present simulation can not be directly extrapolated to the actual tilting conditions. However, it may be of notice that patients with severe congestive heart failure show a significant drop in systolic AP upon orthostatic tilt regardless of relatively maintained baseline AP [34]. The exogenous disturbance in the opposite direction may also be hazardous. Even though baseline AP appears normal in stable CHF patients, noncompliance with salt and water restriction may easily collapse the AP regulation and induce acute decompensation [19].
Study perspective
The progression of CHF is closely associated with the autonomic imbalance between sympathetic and vagal nerve activities. Although pharmacological treatments such as β-adrenergic blockers [35], angiotensin converting enzyme inhibitors [36], and angiotensin receptor blockers [37] have been developed, the therapeutic outcome is not always satisfactory and the mortality rate for CHF remains high. Novel therapeutic strategies for CHF beyond conventional treatments such as electrical vagal nerve stimulation are expected to be developed [7, 38].
Although the mechanisms for the autonomic imbalance in CHF are not completely understood, the loss of baroreflex control of the sympathetic and vagal systems has been considered an important factor [39, 40]. Zucker et al. [41] demonstrated that chronic baroreceptor activation improved the survival of dogs with pacing-induced heart failure. The study seems rudimentary, however, in that the stimulation intensity of the carotid sinus baroreceptors was adjusted intermittently to obtain a prescribed sympathoinhibitory response. The stimulation could be too strong or too weak during the intervals of adjustments. If the stimulation intensity is feedback-controlled continuously, the therapeutic effect and safety could be enhanced.
Quantification of the dynamic characteristics of a target system is the first step for developing a robust controller [42–44]. Next, the controller should be designed in such a way that the variation of the target system dynamics does not affect the controlling result much [45]. The simulation of the closed-loop AP response in CHF, such as that shown in Fig. 5, will be utilized in designing a robust controller system of the carotid sinus activation.
Limitations
First, we used normal rats without sham operation as the control group. Sham-operated rats may serve as a more appropriate control. Although the surgical operation of thoracotomy itself could have affected the baroreflex function in the CHF group, because over 100 days had elapsed after the surgical operation, the effect of surgical operation might have been limited. Actually, two rats that underwent coronary ligation but did not meet the criteria for CHF retained the maximum gain of the total baroreflex close to the control group (rat A: central venous pressure = 2.28 mmHg, biventricular weight = 1.94 g kg−1, total baroreflex gain = 1.61; rat B: central venous pressure = 1.69 mmHg, biventricular weight = 2.33 g kg−1, total baroreflex gain = 1.67).
Second, the data were obtained under anesthetic conditions. Because the anesthesia affects autonomic nervous activities, the results should be carefully interpreted. Vagotomy and lower perfusion of the brain due to carotid artery occlusion might have also affected the baroreflex function. Nevertheless, the present results would provide a unique clue to integrated understanding of the impaired baroreflex function in CHF.
Third, we did not measure cardiac output or peripheral vascular resistance. Subdividing the peripheral arc into the cardiac and vascular components will be necessary to identify the mechanisms for an approximately linear relationship between SNA and AP in both the control and CHF rats.
Conclusions
The dynamic and static characteristics of the neural arc, peripheral arc, and total baroreflex were analyzed comprehensively in rats with CHF after myocardial infarction. Although the derivative characteristics of the baroreflex neural arc are preserved in the CHF rats, the dynamic AP regulation is depressed in both the magnitude and response speed. The equilibrium diagram indicates that the baroreflex gain may seem preserved in CHF in a small range around the closed-loop operating point. However, the percent recovery of AP and the speed of recovery are reduced progressively as the size of exogenous disturbance increases in the CHF rats. The reserve for AP buffering may be lost in CHF despite the relatively maintained baseline AP.
References
Ikeda Y, Kawada T, Sugimachi M, Kawaguchi O, Shishido T, Sato T, Miyano H, Matsuura W, Alexander J Jr, Sunagawa K (1996) Neural arc of baroreflex optimizes dynamic pressure regulation in achieving both stability and quickness. Am J Physiol 271:H882–H890
Sato T, Kawada T, Inagaki M, Shishido T, Sugimachi M, Sunagawa K (2003) Dynamics of sympathetic baroreflex control of arterial pressure in rats. Am J Physiol Regul Integr Comp Physiol 285:R262–R270
White CW (1981) Abnormalities in baroreflex control of heart rate in canine heart failure. Am J Physiol 240:H793–H799
Wang W, Chen JS, Zucker IH (1990) Carotid sinus baroreceptor sensitivity in experimental heart failure. Circulation 81:1959–1966
Wang W, Chen JS, Zucker IH (1991) Carotid sinus baroreceptor reflex in dogs with experimental heart failure. Circ Res 68:1294–1301
Wang W, Brändle M, Zucker IH (1993) Influence of vagotomy on the baroreflex sensitivity in anesthetized dogs with experimental heart failure. Am J Physiol 265:H1310–H1317
Li M, Zheng C, Sato T, Kawada T, Sugimachi M, Sunagawa K (2004) Vagal nerve stimulation markedly improves long-term survival after chronic heart failure in rats. Circulation 109:120–124
Shoukas AA, Callahan CA, Lash JM, Haase EB (1991) New technique to completely isolate carotid sinus baroreceptor regions in rats. Am J Physiol 260:H300–H303
Sato T, Kawada T, Miyano H, Shishido T, Inagaki M, Yoshimura R, Tatewaki T, Sugimachi M, Alexander J Jr, Sunagawa K (1999) New simple methods for isolating baroreceptor regions of carotid sinus and aortic depressor nerves in rats. Am J Physiol 276:H326–H332
Kawada T, Kamiya A, Li M, Shimizu S, Uemura K, Yamamoto H, Sugimachi M (2009) High levels of circulating angiotensin II shift the open-loop baroreflex control of splanchnic sympathetic nerve activity, heart rate and arterial pressure in anesthetized rats. J Physiol Sci 59:447–455
Marmarelis PZ, Marmarelis VZ (1978) The white noise method in system identification. In: Analysis of physiological systems. Plenum, New York, pp 131–221
Kent BB, Drane JW, Blumenstein B, Manning JW (1972) A mathematical model to assess changes in the baroreceptor reflex. Cardiology 57:295–310
Glantz SA (2002) Primer of biostatistics, 5th edn. McGraw-Hill, New York
Kawada T, Yanagiya Y, Uemura K, Miyamoto T, Zheng C, Li M, Sugimachi M, Sunagawa K (2003) Input-size dependence of the baroreflex neural arc transfer characteristics. Am J Physiol Heart Circ Physiol 284:H404–H415
Yamamoto K, Kawada T, Kamiya A, Takaki H, Shishido T, Sunagawa K, Sugimachi M (2008) Muscle mechanoreflex augments arterial baroreflex-mediated dynamic sympathetic response to carotid sinus pressure. Am J Physiol Heart Circ Physiol 295:H1081–H1089
Kamiya A, Kawada T, Yamamoto K, Mizuno M, Shimizu S, Sugimachi M (2008) Upright tilt resets dynamic transfer function of baroreflex neural arc to minify the pressure disturbance in total baroreflex control. J Physiol Sci 58:189–198
Mohrman DE, Heller LJ (2006) Cardiovascular physiology, 6th edn. McGraw Hill, New York
Sato T, Kawada T, Inagaki M, Shishido T, Takaki H, Sugimachi M, Sunagawa K (1999) New analytic framework for understanding sympathetic baroreflex control of arterial pressure. Am J Physiol 276:H2251–H2261
Lepage S (2008) Acute decompensated heart failure. Can J Cardiol 24(Suppl B):6B–8B
Masaki H, Imaizumi T, Harasawa Y, Takeshita A (1994) Dynamic arterial baroreflex in rabbits with heart failure induced by rapid pacing. Am J Physiol 267:H92–H99
Kawai H, Mohan A, Hagen J, Dong E, Armstrong J, Stevens SY, Liang CS (2000) Alterations in cardiac adrenergic terminal function and β-adrenoceptor density in pacing-induced heart failure. Am J Physiol Heart Circ Physiol 278:H1708–H1716
Kawada T, Shishido T, Inagaki M, Zheng C, Yanagiya Y, Uemura K, Sugimachi M, Sunagawa K (2002) Estimation of baroreflex gain using a baroreflex equilibrium diagram. Jpn J Physiol 52:21–29
Kashihara K, Kawada T, Li M, Sugimachi M, Sunagawa K (2004) Bezold-Jarisch reflex blunts arterial baroreflex via the shift of neural arc toward lower sympathetic nerve activity. Jpn J Physiol 54:395–404
Yamamoto K, Kawada T, Kamiya A, Takaki H, Miyamoto T, Sugimachi M, Sunagawa K (2004) Muscle mechanoreflex induces the pressor response by resetting the arterial baroreflex neural arc. Am J Physiol Heart Circ Physiol 286:H1382–H1388
Kamiya A, Kawada T, Yamamoto K, Michikami D, Ariumi H, Uemura K, Zheng C, Shimizu S, Aiba T, Miyamoto T, Sugimachi M, Sunagawa K (2005) Resetting of the arterial baroreflex increases orthostatic sympathetic activation and prevents postural hypotension in rabbits. J Physiol 566:237–246
Michikami D, Kamiya A, Kawada T, Inagaki M, Shishido T, Yamamoto K, Ariumi H, Iwase S, Sugenoya J, Sunagawa K, Sugimachi M (2006) Short-term electroacupuncture at Zusanli resets the arterial baroreflex neural arc toward lower sympathetic nerve activity. Am J Physiol Heart Circ Physiol 291:H318–H326
Ninomiya I, Nisimaru N, Irisawa H (1971) Sympathetic nerve activity to the spleen, kidney, and heart in response to baroceptor input. Am J Physiol 221:1346–1351
Matsukawa K, Ninomiya I, Nishiura N (1993) Effects of anesthesia on cardiac and renal sympathetic nerve activities and plasma catecholamines. Am J Physiol 265:R792–R797
Yamamoto H, Kawada T, Kamiya A, Kita T, Sugimachi M (2008) Electroacupuncture changes the relationship between cardiac and renal sympathetic nerve activities in anesthetized cats. Auton Neurosci 144:43–49
Rowell LB (1974) Human cardiovascular adjustments to exercise and thermal stress. Physiol Rev 54:75–159
Kawada T, Uemura K, Kashihara K, Jin Y, Li M, Zheng C, Sugimachi M, Sunagawa K (2003) Uniformity in dynamic baroreflex regulation of left and right cardiac sympathetic nerve activities. Am J Physiol Regul Integr Comp Physiol 284:R1506–R1512
Kawada T, Shishido T, Inagaki M, Tatewaki T, Zheng C, Yanagiya Y, Sugimachi M, Sunagawa K (2001) Differential dynamic baroreflex regulation of cardiac and renal sympathetic nerve activities. Am J Physiol Heart Circ Physiol 280:H1581–H1590
Kamiya A, Kawada T, Yamamoto K, Michikami D, Ariumi H, Miyamoto T, Shimizu S, Uemura K, Aiba T, Sunagawa K, Sugimachi M (2005) Dynamic and static baroreflex control of muscle sympathetic nerve activity (SNA) parallels that of renal and cardiac SNA during physiological change in pressure. Am J Physiol Heart Circ Physiol 289:H2641–H2648
Kassis E (1987) Cardiovascular response to orthostatic tilt in patients with severe congestive heart failure. Cardiovasc Res 21:362–368
Packer M, Fowler MB, Roecker EB, Coats AJ, Katus HA, Krum H, Mohacsi P, Rouleau JL, Tendera M, Staiger C, Holcslaw TL, Amann-Zalan I, DeMets DL (2002) Effect of carvedilol on the morbidity of patients with severe chronic heart failure: results of the Carvedilol Prospective Randomized Cumulative Survival (COPERNICUS) study. Circulation 106:2194–2199
Yusuf S, Sleight P, Pogue J, Bosch J, Davies R, Dagenais G (2000) Effects of an angiotensin-converting-enzyme inhibitor, ramipril, on cardiovascular events in high-risk patients. The heart outcomes prevention evaluation study investigators. N Engl J Med 342:145–153
Telmisartan Randomised AssessmeNt Study in ACE iNtolerant subjects with cardiovascular Disease (TRANSCEND) Investigators, Yusuf S, Teo K, Anderson C, Pogue J, Dyal L, Copland I, Schumacher H, Dagenais G, Sleight P (2008) Effects of the angiotensin-receptor blocker telmisartan on cardiovascular events in high-risk patients intolerant to angiotensin-converting enzyme inhibitors: a randomised controlled trial. Lancet 372:1174–1183
Schwartz PJ, De Ferrari GM, Sanzo A, Landolina M, Rordorf R, Raineri C, Campana C, Revera M, Ajmone-Marsan N, Tavazzi L, Odero A (2008) Long term vagal stimulation in patients with advanced heart failure: first experience in man. Eur J Heart Fail 10:884–891
Mancia G, Seravalle G, Giannattasio C, Bossi M, Preti L, Cattaneo BM, Grassi G (1992) Reflex cardiovascular control in congestive heart failure. Am J Cardiol 69:17G–23G
Eckberg DL, Drabinsky M, Braunwald E (1971) Defective cardiac parasympathetic control in patients with heart disease. N Engl J Med 285:877–883
Zucker IH, Hackley JF, Cornish KG, Hiser BA, Anderson NR, Kieval R, Irwin ED, Serdar DJ, Peuler JD, Rossing MA (2007) Chronic baroreceptor activation enhances survival in dogs with pacing-induced heart failure. Hypertension 50:904–910
Sato T, Kawada T, Sugimachi M, Sunagawa K (2002) Bionic technology revitalizes native baroreflex function in rats with baroreflex failure. Circulation 106:730–734
Gotoh TM, Tanaka K, Morita H (2005) Controlling arterial blood pressure using a computer-brain interface. Neuroreport 16:343–347
Kawada T, Shimizu S, Yamamoto H, Shishido T, Kamiya A, Miyamoto T, Sunagawa K, Sugimachi M (2009) Servo-controlled hind-limb electrical stimulation for short-term arterial pressure control. Circ J 73:851–859
Åström K, Hägglund T (1995) PID controllers: theory, design, and tuning, 2nd edn. International Society of Automation, Research Triangle Park, NC
Hunter IW, Korenberg MJ (1986) The identification of nonlinear biological systems: Wiener and Hammerstein cascade models. Biol Cybern 55:135–144
Acknowledgments
This study was supported by Health and Labour Sciences Research Grants (H18-nano-Ippan-003, H19-nano-Ippan-009, H20-katsudo-Shitei-007, and H21-nano-Ippan-005) from the Ministry of Health, Labour and Welfare of Japan; by a Grant-in-Aid for Scientific Research (No. 20390462) from the Ministry of Education, Culture, Sports, Science and Technology of Japan; and by the Industrial Technology Research Grant Program from the New Energy and Industrial Technology Development Organization (NEDO) of Japan.
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendix 1: feature of a Gaussian white noise input
The neural arc of the arterial baroreflex may be roughly modeled by a cascade system, which consists of a dynamic linear subsystem followed by a static nonlinear subsystem [14]. In this type of system, the selection of the mean input pressure and the amplitude of input does not significantly affect the estimation of the system linear dynamics except for a factor of proportionality when the system is tested with a Gaussian white noise (GWN) signal [46]. In order to demonstrate the above notion, we performed a simulation study. In reference to Fig. 6a, the dynamic linear subsystem for the baroreflex neural arc was modeled by the following equation:
where f C1 and f C2 are the corner frequencies determining the frequency-dependent changes in the dynamic gain, L denotes the pure dead time, f and j indicate frequency and an imaginary unit, respectively. We set f C1 = 0.1 Hz, f C2 = 0.9 Hz, and L = 0.1 s to mimic the neural arc transfer function. The static nonlinear subsystem of the neural arc was modeled by the logistic function (Eq. 3) with parameters of P 1 = 60%, P 2 = 0.1 mmHg−1, P 3 = 120 mmHg, and P 4 = 40%. Using this model, the linear input-output relationship was estimated by a GWN input or a binary white noise (BIN) input. Mean input pressure was changed among 120, 90, and 150 mmHg (P 3, P 3 – 30, and P 3 + 30 mmHg).
Figure 6b through 6d shows the estimation results. In each panel, thin smooth curves in the gain and phase plots indicate the transfer function of the dynamic linear subsystem described by Eq. 5. The solid bold curves, dotted thin curves, and dashed thin curves represent the estimation results obtained by mean input pressures of 120, 90, and 150 mmHg, respectively. Figure 6b and c represents the estimation results obtained by GWNs with standard deviations of 10 and 20 mmHg, respectively. As can be seen, the estimated gain plots are nearly parallel to the gain plot of the given dynamic linear subsystem. The estimated phase plots are superimposable on the phase plot of the given dynamic linear subsystem.
Figure 6d shows the estimation results obtained by BINs with an amplitude of 10 mmHg (i.e., the peak-to-peak pressure of 20 mmHg). None of the estimated gain plots are parallel to the gain plot of the given dynamic linear subsystem. The estimated phase plots also deviated from the phase plot of the given dynamic linear subsystem. When the system was tested by the BIN input with the mean input pressure of 120 mmHg, the coherence value became close to unity, suggesting that the BIN input had caused on-off behavior in the system output.
The above simulation results confirm that the selection of the mean input pressure and the amplitude of GWN does not significantly affect the estimation of the dynamic characteristics of the system except for a factor of proportionality when the mean input pressure and the standard deviation of GWN are not far beyond the system operating range.
Appendix 2: estimation of the initial slope of the step response or recovery response
To characterize the speed of the response, the initial slope of the response was calculated as follows. First, a response threshold is determined as the 5% value of the steady-state step response. In the case of the recovery response, 5% of the maximum recovery is used as the threshold. Next, an initial time point at which the response exceeds the threshold value is determined. Starting from this initial time point, a linear regression analysis giving a slope and an intercept is performed repeatedly while increasing the number of analyzed data points. The steepest slope thus obtained is used as the initial slope of the response. By setting the 5% threshold, the linear regression can be performed excluding the initial data points constituting the dead time.
Appendix 3: comparison of splanchnic and renal SNAs under the present experimental conditions
We compared splanchnic and renal SNAs under the present experimental conditions in an additional three rats. A branch of the right renal sympathetic nerve was exposed through a right flank incision, and renal SNA was recorded simultaneously with splanchnic SNA. Figure 7a presents time series of 10-Hz decimated CSP, splanchnic SNA, and renal SNA obtained from one rat. White lines in splanchnic and renal SNAs indicate 2-s moving averaged signals. In each of the splanchnic and renal SNAs, a 10-s averaged value for the last 10 s at CSP of 60 mmHg was defined as 100%, and a 10-s averaged value after the administration of hexamethonium bromide was defined as 0%, in the same manner as the main study. Figure 7b shows the relationships of SNAs versus CSP. Although the renal SNA values were slightly greater than the splanchnic SNA values at CSPs of 160 and 180 mmHg in this rat, the general sigmoidal relationships were similar between the two SNAs. The regression line for renal SNA (y) versus splanchnic SNA (x) was y = 0.94x + 5.8 (r 2 = 0.998) in this rat (Fig. 7c). The other two rats showed regression lines of y = 1.01x + 0.3 (r 2 = 0.999) and y = 0.99x + 3.7 (r 2 = 0.995), suggesting that both splanchnic and renal SNAs can represent common activity for the AP regulation during the stepwise input protocol under the present experimental conditions.
Rights and permissions
About this article
Cite this article
Kawada, T., Li, M., Kamiya, A. et al. Open-loop dynamic and static characteristics of the carotid sinus baroreflex in rats with chronic heart failure after myocardial infarction. J Physiol Sci 60, 283–298 (2010). https://doi.org/10.1007/s12576-010-0096-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s12576-010-0096-9