Journal of Applied Physiology Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 102: 1348-1356, 2007. First published December 21, 2006; doi:10.1152/japplphysiol.00158.2006
8750-7587/07 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
102/4/1348    most recent
00158.2006v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by van de Vooren, H.
Right arrow Articles by Van der Wall, E. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by van de Vooren, H.
Right arrow Articles by Van der Wall, E. E.

Baroreflex sensitivity, blood pressure buffering, and resonance: what are the links? Computer simulation of healthy subjects and heart failure patients

Hedde van de Vooren,1 Maaike G. J. Gademan,1 Cees A. Swenne,1 Ben J. TenVoorde,2 Martin J. Schalij,1 and Ernst E. Van der Wall1

1Cardiology Department, Leiden University Medical Center, Leiden; and 2Department of Physics and Medical Technology, Vrije Universiteit Medical Center, Amsterdam, the Netherlands

Submitted 7 February 2006 ; accepted in final form 11 December 2006


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The arterial baroreflex buffers slow (<0.05 Hz) blood pressure (BP) fluctuations, mainly by controlling peripheral resistance. Baroreflex sensitivity (BRS), an important characteristic of baroreflex control, is often noninvasively assessed by relating heart rate (HR) fluctuations to BP fluctuations; more specifically, spectral BRS assessment techniques focus on the BP-to-HR transfer function around 0.1 Hz. Skepticism about the relevance of BRS to characterize baroreflex-mediated BP buffering is based on two considerations: 1) baroreflex-modulated peripheral vasomotor function is not necessarily related to baroreflex-HR transfer; and 2) although BP fluctuations around 0.1 Hz (Mayer waves) might be related to baroreflex BP buffering, they are merely a not-intended side effect of a closed-loop control system. To further investigate the relationship between BRS and baroreflex-mediated BP buffering, we set up a computer model of baroreflex BP control to simulate normal subjects and heart failure patients. Output variables for various randomly chosen combinations of feedback gains in the baroreflex arms were BP resonance, BP-buffering capacity, and BRS. Our results show that BP buffering and BP resonance are related expressions of baroreflex BP control and depend strongly on the sympathetic gain to the peripheral resistance. BRS is almost uniquely determined by the vagal baroreflex gain to the sinus node. In conclusion, BP buffering and BRS are unrelated unless coupled gains in all baroreflex limbs are assumed. Hence, the clinical benefit of a high BRS is most likely to be attributed to vagal effects on the heart instead of to effective BP buffering.

autonomic nervous system; cardiovascular variability; Mayer waves; spectral analysis; transfer function


IN DAILY LIFE, multiple processes perturb blood pressure. The duration of these challenges varies widely. For example, respiration makes blood pressure fluctuate with every breath (13), while physical or mental stress elevates blood pressure for minutes or even longer. The arterial baroreflex is a negative-feedback mechanism that effectively buffers such incidental blood pressure fluctuations (11, 20, 21, 23). In negative-feedback systems, feedback delay often causes resonance in a given frequency band; this is the price to be paid for effective buffering at other frequencies. Resonance in blood pressure (5, 8, 12, 31, 49) manifests in the form of the well-known Mayer (22, 33) waves (beat-to-beat blood pressure oscillations with a frequency of ~0.1 Hz/period of ~10 s). Effective baroreflex blood pressure buffering occurs below the Mayer frequency (10, 16).

Besides a sympathetic limb that modulates peripheral resistance, the baroreflex has also sympathetic and parasympathetic (vagal) limbs that influence cardiac contractility, venous return, and cardiac rhythm. Usually, baroreflex functioning is characterized by baroreflex sensitivity (BRS). This index of baroreflex vigor is defined as the reflex-induced change in interbeat interval in milliseconds per millimeter of Hg blood pressure change (14, 34, 36, 44). The prognostic value of BRS and the favorable consequences of successful interventions with BRS have amply been demonstrated (27, 28).

Little is known, however, about the representativeness of this index for the efficacy of blood pressure buffering. There are two reasons to be skeptical in this respect: 1) by definition—interbeat interval change per unit blood pressure change—BRS is bound to characterize baroreflex-mediated effects on the heart, while the baroreflex buffers blood pressure mainly by controlling peripheral resistance (2, 30); and 2) oftentimes being assessed in the Mayer frequency range of spontaneous heart rate and blood pressure fluctuations (15, 39), BRS might represent resonance rather than buffering baroreflex characteristics.

We addressed these skepticisms by simulations with a hybrid mathematical model of baroreflex blood pressure and heart rate control composed of hemodynamic elements that are evaluated on a beat-to-beat basis, linked to a time-continuous modeled neural control part. By changing some parameter settings, the model mimics physiological as well as pathological hemodynamic and autonomic conditions.

By simulating with various gain combination values, we quantified the role of the sympathetic and parasympathetic gains in the three baroreflex limbs for blood pressure variability (BPV) and heart rate variability (HRV) under physiological and pathological conditions. From the obtained systolic blood pressure and interbeat interval values, relations between BRS and blood pressure buffering and between blood pressure buffering and resonance were examined.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The simulation model we used for this study represents short-term human blood pressure control without breathing modulation. It is tuned for supine posture. This model, programmed in Matlab Simulink (MathWorks, Natick, MA), is, apart from some modifications, similar to the model as earlier designed and validated by TenVoorde and Kingma (46).

Model Description

A gross overview of the autonomically controlled model is given in Fig. 1 (see Tables 1 and 2 for abbreviations and model parameters). The model represents the systemic circulation and consists of three sections: a hemodynamic section, a baroreceptor section, and an autonomic control section. The model generates output in the form of time-dependent systolic blood pressure (SBP) values (mmHg) and interbeat interval (IBI) values (ms) by using a sinusoidal pressure probe (frequency adjustable, amplitude fixed at 1 mmHg) as an input signal. This apparently small perturbation at the input of the baroreflex produces reduced SBP fluctuations (amplitude <1 mmHg, buffering) but also Mayer waves (amplitude >1 mmHg, resonance), depending on the frequency of the pressure probe.


Figure 1
View larger version (22K):
[in this window]
[in a new window]

 
Fig. 1. Model of baroreflex arterial blood pressure control. IPFM, integral pulse frequency modulator. Model variables and model parameters are described in Tables 1 and 2. Model sections are discussed in METHODS. Adjustable parameters are placed in white boxes or triangles; gray boxes represent fixed model divisions.

 

View this table:
[in this window]
[in a new window]

 
Table 1. Model parameters and values under physiological and pathological conditions

 

View this table:
[in this window]
[in a new window]

 
Table 2. Model variables: sample and hold (n) or time dependent (t)

 
Hemodynamic section.   In the hemodynamic section, all signals are sample-and-hold signals: the beat-to-beat varying cardiovascular signals are modeled in elementary difference equations. All values are adapted when a new heartbeat emerges. Stroke volume Qn is produced by the one-chamber Starling heart. It depends on interbeat interval In, venous return volume Vn, and a contractility volume term Cn

Formula
where {delta}n is a left ventricle filling factor:

Formula
As this model will only be used to simulate different autonomic control states rather than different hemodynamic states (like standing posture), changes in cardiac contractility and venous return appear to generate only very small fluctuations in stroke volume (<5%). Therefore, we simplified the above relation into:

Formula
where Starling heart filling parameter Vref indicates the stroke volume when {delta}n = 1.

Stroke volume Qn, assuming a constant arterial compliance (CA), determines pulse pressure Pn by:

Formula
A systemic Windkessel simulates diastolic blood pressure Dn:

Formula
The Windkessel time constants Tn is controlled by the baroreflex (see autonomic control section) and is directly associated with total peripheral resistance. Although it is usual to compute diastolic pressure as the exponential decay of systolic pressure, we used this slightly modified formula to obtain more accurate systolic blood pressure values. Finally, systolic pressure Sn is computed by adding Pn and Dn.

Baroreceptor section.   The baroreceptors are modeled linearly within the range between a threshold of 90 mmHg and a saturation level of 150 mmHg. At the baroreceptors, the systolic blood pressure Sn is compared with a low-pass filtered systolic blood pressure reference value. This value functions as a dynamic blood pressure set point, mimicking the physiological process of baroreceptor resetting (47, 50, 51). The pressure variability source is added at the input of the baroreflex, rendering a sample-and-hold systolic blood pressure variability signal SBP, the first model output signal.

Autonomic control section.   In the time-continuous autonomic control section, SBP is converted into an afferent neural signal NE by factorizing this signal by the baroreceptor sensitivity coefficient {alpha}N. This signal serves as input for three effectors: vagal heart rate control (output: vagal signal n); sympathetic heart rate control (output: sympathetic signal m); and sympathetic peripheral resistance control (output: Windkessel time constant T). The vagal signal n represents vagal heart rate deceleration (0 < n < 1), while the sympathetic signal m represents sympathetic heart rate acceleration (m > 1).

The three effectors are modeled in frequency-dependent functional blocks, with specific sensitivity coefficients, time constants, time delays and by autonomic tones (N, MH, and MPR; see Table 1 for actual values). In addition to these model parameters, extra baroreflex gain multipliers (SH, V, and SPR) were added to strengthen or weaken the role of each baroreflex effector.

The neural time-continuous part and the hemodynamic beat-to-beat part are interconnected by an integral pulse frequency modulator (IPFM), which simulates cardiac pacemaker function (18). Rosenblueth and Simeone (40) have demonstrated that combined sympathetic and vagal influences on the sinus node contribute to the actual heart rate R according to the following relationship: R = R0·m·n, where R0 is the intrinsic heart rate. Integration of incoming neural activity results in the generation of the heart interval length In (18). This interbeat interval IBI is the second model output signal.

Adjustable Model Parameters

Thus the model is controlled by seven parameters: one (Vref) for stroke volume, three (MH, N, and MPR) for autonomic outflow, and three (SH, V, and SPR) multipliers for the gains in the three baroreflex limbs.

The first four parameters for stroke volume and autonomic outflow were set as two fixed combinations (Table 1) to represent either normal physiological or abnormal pathological resting conditions. With an increased sympathetic tone to the heart and to the peripheral resistance and decreased parasympathetic tone and reference stroke volume, the pathological parameter settings represent a serious pathological condition resembling congestive heart failure. Compared with the physiological conditions, the resting heart rate is higher (90 beats/min instead of 60 beats/min), and the average SBP is slightly lower (114 mmHg instead of 120 mmHg).

The last three parameters serve as potentiometers (multipliers) on the vagal and sympathetic baroreflex gains to the heart and to the peripheral resistance; V = SH = SPR = 1 is the reference value that is to represent a normally working baroreflex. When one of these parameter values equals 0, the corresponding limb of the baroreflex does not react to changes of SBP with respect to the reference value, and the corresponding effector output becomes the (fixed) tone. A value of 0.5 corresponds to weak involvement. The maximum value of these parameters is 3; this value corresponds to a strong involvement of a given baroreflex limb, e.g., as found in highly trained subjects.

Simulations and Frequency Characteristics

For a given combination of the seven model parameter values, 100 simulation runs were done. A single simulation run served to determine one SBP variability (BPV) frequency component, one IBI variability (HRV) frequency component, and the modulus of the SBP-to-IBI transfer function (TF, necessary to compute BRS) frequency component, at a given frequency of the sinusoid pressure probe. A single simulation run was executed as follows. First, the model was run until steady-state conditions were met. Cubic splines were then fitted through the resulting output signals to obtain the amplitudes of the SBP and IBI fluctuations caused by the pressure probe. Finally, the corresponding TF frequency component was computed by dividing HRV (the amplitude of the IBI fluctuations) by BPV (the amplitude of the SBP fluctuations). The 100 simulation runs were done to construct the complete frequency characteristics of BPV, HRV, and TF by computing all frequency components between 0.003 and 0.300 Hz (step 0.003 Hz).

A total of 162 frequency characteristics of HRV, BPV, and TF were made for both the physiological as well as for the pathophysiological conditions. These 162 frequency characteristics were made to represent 162 different combinations of baroreflex gain multiplier settings. One-hundred fifty gain multiplier combinations were randomly chosen to simulate uncoupled baroreflex gains (values between 0 and 3 from uniform distributions for V, SH, and SPR). In addition, 12 other V/SH/SPR combinations were made to simulate coupled baroreflex gains (0.5/0.5/0.5, 1/1/1, 1.5/1.5/1.5, 2/2/2, 2.5/2.5/2.5, 3/3/3). Besides these multiplier combinations, an extra set of simulation results (obtained with V/SH/SPR combinations 0/1/1, 3/1/1, 1/0/1, 1/3/1, 1/1/0, 1/1/3) was made for the generation of Fig. 2.


Figure 2
View larger version (36K):
[in this window]
[in a new window]

 
Fig. 2. Examples of simulated heart rate variability (HRV; A1, B1, C1), blood pressure variability (BPV; A2, B2, C2), and transfer function (TF; A3, B3, C3) frequency characteristics. Note that frequency axes are log scaled. For these simulations, the model parameters for stroke volume and autonomic outflow (Vref, MH, MPR, and N), were kept fixed at physiological values (see Table 1). The baroreflex gain multipliers SH, V, and SPR were varied around SH/V/SPR = 1/1/1 (solid line in all panels) in the following way. A1, A2, and A3: dotted line, SH/V/SPR = 0/1/1 (inactivated SH); dashed line, SH/V/SPR = 3/1/1 (strong SH). B1, B2, and B3: dotted line, SH/V/SPR = 1/0/1 (inactivated V); dashed line, SH/V/SPR = 1/3/1 (strong V). C1, C2, and C3: dotted line, SH/V/SPR = 1/1/0 (inactivated SPR); dashed line, SH/V/SPR = 1/1/3 (strong SPR). HRV and BPV amplitudes have to be related to the driving force of the sinusoidal pressure probe (1 mmHg). HRV = amplitude of interbeat interval (IBI) fluctuations; BPV = amplitude of systolic blood pressure (SBP) fluctuations; TF = modulus of the SBP-to-IBI transfer function.

 
Main Derived Simulation Variables: BRS, SBP-Buffering Capacity, SBP Resonance

After having computed a full BPV, HRV, and TF characteristic, we determined the following variables. BRS was computed as the averaged magnitude of TF in the low-frequency (LF, 0.05–0.15 Hz) band (15, 37, 39). This band incorporates the Mayer frequencies. SBP-buffering capacity was expressed as the amplitude of the original perturbation (the 1-mmHg sinusoidal pressure probe) divided by the BPV amplitude at the lowest simulated frequency component [0.003 Hz, which is still well above the baroreceptor resetting frequency (47)]; e.g., when the BPV at the lowest frequency component had an amplitude of 0.25 mmHg, the buffering capacity was 4. Maximal SBP resonance (in the LF band) was expressed as the maximal BPV divided by the amplitude of the original perturbation. To determine the relative importance of V, SH, and SPR for baroreflex sensitivity and blood pressure buffering and resonance, multiple linear regressions were done. For these regressions, only the simulations made with random generated baroreflex gain multipliers were used.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The simulation results obtained under physiological and pathological conditions (see Table 1) differ quantitatively (more outspoken characteristics under physiological conditions) rather than qualitatively: all frequency characteristics are smooth, and buffering occurs at the lowest frequencies while resonance occurs at the Mayer frequency around 0.1 Hz. Figure 2 displays examples of some HRV, BPV, and TF frequency characteristics obtained under physiological conditions. Figure 2 consists of three sets of HRV, BPV, and TF frequency characteristics, in each of which one of the three effectors was weakened or strengthened; i.e., baroreflex gain multiplier V, SH, or SPR was increased to 3 (strong) or reduced to 0 (inactivated) with respect to the default value of 1 (normal), while the other two baroreflex gain multipliers were kept at their default values of 1 (normal).

Figure 2A3 shows an unexpected influence of sympathetic heart rate control on IBI: the transfer function in the LF band (i.e., BRS) even decreases when control is strengthened (strong SH). Obviously, blood pressure buffering and resonance are completely insensitive for changes in the sympathetic gain to the heart (Fig. 2A2).


Figure 3
View larger version (28K):
[in this window]
[in a new window]

 
Fig. 3. Main graphical representations of the simulation results. A: SBP-buffering capacity increases with increasing sympathetic baroreflex gain to the peripheral resistance. B: baroreflex sensitivity (BRS) increases with increasing vagal baroreflex gain to the heart. C: SBP resonance increases with increasing SBP-buffering capacity. D: SBP-buffering capacity is only weakly related to baroreflex sensitivity. Filled (black) circles and filled (gray) squares, physiological conditions; open circles and squares, pathological conditions; circles (small), simulations with uncoupled (random generated) baroreflex gain multiplier combinations; squares (larger), simulations with coupled baroreflex gain multiplier combinations: SH/V/SPR = 0.5/0.5/0.5, 1/1/1, 1.5/1.5/1.5, 2/2/2, 2.5/2.5/2.5, 3/3/3. Dashed lines, linear regressions in the random data (hence, scouting simulations with coupled baroreflex gain multiplier combinations excluded). See text for further explanation.

 
Figure 2, B1–B3, shows how the HRV, BPV, and TF frequency characteristics react when the vagal heart rate control is weakened or strengthened (multiplier V assumes the value 0 or 3, respectively; multipliers SH and SPR are kept at a value of 1). Here, the impression arises that multiplier V strongly influences HRV and the BRS, while it does not affect the resonance and buffering behavior (relatively little differences in resonance and buffering are seen in Fig. 2B2).

Figure 2, C1–C2, shows the striking effect of a strengthened sympathetic peripheral resistance control (multiplier SPR assumes the value of 3; multipliers V and SH are kept at a value of 1) on the HRV and BPV frequency characteristics. Figure 2C2 shows that the original sinusoidal disturbance of SBP by the 1-mmHg pressure probe (see Fig. 1) is strongly weakened (buffered) for the lowest frequencies, is amplified (resonance) over nearly the whole LF band, and returns to ~1 mmHg for higher frequencies. The larger part of this effect, especially the resonance phenomenon, disappears under normal control (multiplier SPR assumes the value of 1). The frequency characteristic is almost flat when control is absent (multiplier SPR assumes the value of 0). The shapes of the HRV frequency characteristics in Fig. 2C1 grossly follow the BPV characteristics. As expected, the TF frequency characteristics (Fig. 2C3) are very much similar for all three SPR values, 0, 1, and 3. In summary, from the frequency characteristics shown in Fig. 2, C1–C3, the impression arises that sympathetic peripheral resistance control strongly influences resonance and buffering while it does not affect the TF or BRS.

Figure 3 depicts the strongest relations between vagal and sympathetic baroreflex gains, SBP-buffering capacity, SBP resonance, and BRS on the basis of the results of multiple linear regression analysis. It is pointed out that in a physiological setting, 83% of the variance in SBP buffering was attributable to sympathetic peripheral resistance control (multiplier SPR); under pathological conditions this was 78%. Also, 99% of the variance in BRS was attributable to vagal heart rate control (multiplier V); under pathological conditions this remained 99%.

The scatterplot of the SBP-buffering capacity as a function of SPR, together with linear fits for the physiological and the pathological data (Fig. 3A), shows close to perfect linear relationships. Also, there is little difference between the linear fits for the physiological and the pathological simulation results. Obviously, heart rate control, but also the settings of Vref and MPR, were of minor importance for blood pressure buffering.

The scatterplot of BRS as a function of V, together with linear fits for the physiological and the pathological data (Fig. 3B), shows nearly perfect linear relationships. Here, the physiological fit (slope 6.9 ms/mmHg) and the pathological fit (slope 4.0 ms/mmHg) differ considerably: with equal vagal gain multipliers, BRS is much larger in physiological conditions.

Figure 3C shows that SBP-buffering capacity and SBP resonance have a convex relationship and that the resonance phenomenon is much more prominent in physiological circumstances compared with pathological conditions. The strong link between buffering and resonance follows directly from regression analysis: also here, multiplier SPR attributes the most to variance in SBP resonance (95% under normal conditions, 91% under pathological conditions).

Finally, Fig. 3D shows that BRS was almost unrelated to SBP-buffering capacity, unless coupled baroreflex gains (simulation results represented by the open and solid squares) are assumed. The squared correlation coefficients of the linear regressions of SBP-buffering capacity on BRS were as low as 0.037 (physiological conditions) and 0.083 (pathological conditions).


    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We used a mathematical model to investigate the relation between baroreflex sensitivity (BRS, an index of baroreflex vigor) and baroreflex-mediated blood pressure-buffering capacity. This relation is not straightforward since the involved efferent baroreflex limbs (vagal and sympathetic pathways to the heart, and sympathetic pathways to the peripheral vasculature, respectively) differ. Moreover, baroreflex buffering occurs at lower frequencies than the Mayer frequency band in which BRS is noninvasively assessed and in which blood pressure resonates. Whether resonance disturbs the transfer function, thus precluding reliable BRS assessment in the Mayer frequency band, is not known. Also it is not clear what the relation is between, on the one hand, the "desired" phenomenon of blood pressure buffering at frequencies lower than the Mayer frequency (10, 16) and, on the other hand, the phenomenon of blood pressure resonance in the Mayer frequency band [nothing more than a by-product of baroreflex-mediated blood pressure control (10)].

Simulations with various combinations of baroreflex gains, under physiological as well as under pathological conditions (increased sympathetic tone, decreased vagal tone, reduced cardiac stroke volume), yielded frequency characteristics of the transfer function, of HRV, and of BPV, and values of BRS, blood pressure-buffering capacity and baroreflex resonance in a wide range of conditions that may be met in real life. All frequency characteristics had a smooth character, and even with striking resonance the transfer function did not show discontinuous or deviant behavior compared with its value below and above the frequency band of resonance (see Fig. 2). In the following, the simulation results will be discussed in the order they have been presented in Fig. 3.

Baroreflex Gains and Blood Pressure-Buffering Capacity

Our results suggest a predominant role for the sympathetic limb to the peripheral vasculature for the blood pressure-buffering capacity of the arterial baroreflex (Fig. 3A). There is almost no difference in buffering capacity between the physiological and the pathological conditions. This result clearly illustrates the fact that efficacy of baroreflex-mediated blood pressure control rests on dynamic control of the peripheral resistance. Modulation of heart rate by baroreflex-mediated modulation of the vagal and sympathetic tone to the heart is not very relevant for blood pressure control in the frequency range of interest for this study (0.05–0.3 Hz).

Obviously, the simulation results may not be interpreted in such a way that baroreflex-mediated blood pressure buffering in patients is not different from that in healthy subjects. The sympathetic feedback gain to the peripheral vasculature is the decisive factor here. We speculate that this gain will be lower in patients. Hence, it may have been somewhat unrealistic to extend the simulations in pathological conditions to a similar value of SPR than the simulations in physiological conditions. The consequence of our speculation would be that the blood pressure-buffering capacity in patients is smaller than that in healthy subjects.

Baroreflex Gains and BRS

Baroreflex sensitivity is linear with and depends almost exclusively on the vagal feedback gain to the heart (Fig. 3B). The slopes of the linear regressions (6.9 and 4.0 ms/mmHg with physiological and pathological conditions, respectively) are merely to be explained on the basis of heart rate differences between these two situations and the way the IPFM (18) reacts to fluctuations in vagal tone. The fact that BRS depends on heart rate has been recognized earlier (1), and proposals have been done to normalize BRS on heart rate, or, alternatively, to express BRS in beats per minute per millimeter Hg instead of in milliseconds per millimeter Hg. Such arithmetic operations would change the linear relationships in Fig. 3B to curved ones but leave the conclusions unaffected that BRS increases with increasing vagal feedback gain and that the vagal feedback gain almost uniquely determines BRS.

The predominant role of the vagal feedback gain on the baroreflex sensitivity (38) can also be formulated in a slightly different way: due to the differences in the time constants of the vagal and the sympathetic branches (in our model 0.1 and 4.0 s, respectively), the greater part of HRV is simply vagal transmission of blood pressure variability to the sinus node. This is easily perceived in Fig. 2, B1–B3 and C1–C3, and in accordance with the findings of Cevese et al. (9). When vagal feedback gain is zero (dotted lines in Fig. 2, B1–B3), there is almost no HRV (Fig. 2B1) while there still is appreciable BPV (Fig. 2B2). When there is appreciable vagal baroreflex feedback gain (solid and dashed lines in Fig. 2, B1–B3, and all lines in Fig. 2, C1–C3), the HRV frequency characteristics in Fig. 2, B1 and C1, have the same shape as the BPV frequency characteristics in Fig. 2, B2 and C2. In the case of overt (sympathetically mediated) blood pressure resonance, where the BPV frequency characteristic has a narrow peak (Fig. 2C2), a similar "monochrome" HRV frequency characteristic is seen in Fig. 2C1. Alternatively, when there is no outspoken resonance (Fig. 2B2), there is "broad band" HRV (Fig. 2B1).

Blood Pressure-Buffering Capacity and Mayer Waves

Figure 3C illustrates the principle that a negative-feedback control system with feedback delay buffers the controlled variable at certain frequencies at the cost of resonance at other frequencies. The baroreflex blood pressure control system as simulated in this study behaves in a way that is very similar to what was experimentally observed (10). Blood pressure buffering, a major function of the baroreflex, occurs at frequencies below the Mayer waves (resonance in the LF band, Fig. 2C2). Typically, the arterial baroreflex could dampen blood pressure and heart rate responses to stressors that last several minutes. On the one hand, neural control of blood pressure by sympathetically induced vasoconstriction is relatively fast (seconds). On the other hand, baroreceptor resetting (47, 50, 51) limits the maximal duration during which baroreflex mediated buffering of a stressor may continue. In our simulations, the BPV frequency characteristics in Fig. 2, A2, B2, and C2, shows that dampening (reduction of the sinusoidal pressure probe disturbance) is strongest for the lowest frequencies.

Although there still exists some controversy about the origin of the observed spontaneous blood pressure and heart rate variations around 0.1 Hz (32), we assume that this phenomenon is due to the dynamics of the closed-loop vasomotor control (arterial peripheral resistance), in which the time delay of a few seconds plays the major role. Building a baroreflex model with negative-feedback control, and with parameters estimated from physiologically known data, results in a model that simply shows such resonance behavior, without the need to postulate centrally driven oscillators or (strong) nonlinearities.

Resonance, the price to be paid for buffering, is likely to be useless in terms of homeostasis. At the same time, it may be an innocent phenomenon without any negative impact for the organism (22). The fact that Mayer waves, useless or not, exist, facilitates spectral BRS assessment in the LF band, by creating an input signal (blood pressure variability) for the baroreflex of which the output signal (heart rate variability) can easily be measured. There is no inherent signal analysis problem in measuring BRS by the transfer function around the resonance frequency. However, the 180° phase shift caused by the time lag in the sympathetic efferent baroreflex limb to the heart with respect to the phase shift in the efferent vagal limb, that has a much shorter time lag, may cause the sympathetic and vagal limbs to the heart to counteract in the LF band. This effect will become prominent with increased sympathetic gain to the heart (see, e.g., Fig. 2A1, dashed line). In this respect, lower TF frequencies would constitute a more realistic BRS estimate, because here vagal and sympathetic feedback to the heart is concordant (Fig. 2A3, dashed line). In general, TF values in the LF band are not too different at higher frequencies; TF values increase for lower frequencies (Fig. 2, A3, B3, and C3).

Baroreflex Sensitivity and Blood Pressure-Buffering Capacity

One of the major reasons to perform this study was the question of whether there is a relation between the primary function of the baroreflex, i.e., blood pressure buffering, and the generally accepted clinical index for baroreflex vigor, BRS. Figure 3D shows that this relation does almost not exist. The correlation coefficients of the regression lines of SBP-buffering capacity on BRS are very low, and the data are diffusely distributed.

Indeed, vagal control of heart rate (major cause of the BPV-to-HRV transfer and, hence, major determinant of BRS) and sympathetic modulation of the peripheral vasculature (major cause of peripheral resistance adaptations and, hence, a major determinant of blood pressure buffering) become effective via separate efferent pathways of the baroreflex. There should not necessarily be a fixed relationship between the feedback gains in both reflex limbs (43).

To our knowledge, there are no data regarding the relative strengths of the gains in the three baroreflex effector limbs. It might well be that subjects with a low gain in the vagal limb have also low gains in the sympathetic limbs, amongst others, because part of the origin of these gains is to be found in the common afferent pathway of the reflex starting at the baroreceptors in the arterial wall up to and including the nucleus tractus solitarii in the brain stem. Inspection of the simulation results obtained under coupled gains (closed and open squares in Fig. 3D) reveals that in such cases there is a seemingly linear relationship between BRS and blood pressure-buffering capacity in healthy subjects as well as in patients.

Limitations of the Model

Basic to our study is the representativeness of the mathematical model. The original model has extensively been validated (46), among others by comparing the results of modeled vagal blockade and of standing with real-world observations. The modified model, however, has a simplified hemodynamic structure. Since the simulations addressed blood pressure and heart rate control in the supine posture only, the dynamic control capabilities of cardiac contractility and venous return on cardiac output and hence, blood pressure, have completely been removed (obviously, such a simplification cannot be made in cases where the average IBI changes because of an altering circulatory load). Elimination of these feedforward mechanisms enabled us to concentrate on the role of the various baroreflex gains, especially in the LF band, rather than steady-state phenomena in the lower frequencies. As our simulation results are still comparable with the various spectra produced by the original model, we do believe that our model still produces relevant spectra.

The modified model that was used for our present study generates and explains some situations that are known from the clinic. It is obvious that the resonance phenomenon in the LF band, generally known as Mayer waves (33), is strongly under influence of the baroreflex. The only situations in which Mayer waves hardly appear is when the sympathetic baroreflex gain to the peripheral resistance is small (see Fig. 2C2, dotted line). This simulation result parallels studies in rats (19) and in humans (42, 45). The relevance of the model is underscored by the observation that SBP variability in the LF band decreases for a fixed heart rate (results not shown here). This phenomenon was described by Taylor and Eckberg (45) in a study in humans. The authors demonstrated that baroreflex-mediated heart rate control was not effective in reducing blood pressure variability that had a larger amplitude in sinus rhythm compared with fixed-interval atrial pacing.

Within the operating space constituted by the ranges of the parameters as given in Table 1, in combination with baroreflex gain multiplier values between 0 and 3, our model can be used without any difficulty. For example, as the baroreflex gain to the peripheral resistance (SPR) should not have any influence on the transfer function, Fig. 2C3 shows indeed that only varying SPR produces almost the same transfer functions. The minimal differences between those functions can be explained by nonlinearities in the model. Further expansion of the operating space may therefore be not allowed. Furthermore, higher baroreflex gains would no longer be realistic and lead to, e.g., unacceptably high blood pressure variability values.

BRS can be enhanced by training (7), and the beneficial effects of a thus increased BRS have convincingly been demonstrated (27). How this effect is accomplished remains uncertain. Inhibition of stressor-induced heart rate increases may be one reason; the vagal feedback gain in the cardiac efferent limb may predominantly cause this effect. Inhibition of stressor-induced blood pressure increases may be another reason; the sympathetic feedback gain in the baroreflex efferent limb to the peripheral vasculature may predominantly cause this effect. Both effects could help to inhibit a stressor-induced raise of myocardial oxygen consumption, which is proportional to the product of heart rate and systolic blood pressure (3, 29).

A final remark regards the phenomenon as seen in Fig. 2A3. It appears that BRS (the TF between 0.05 and 0.15 Hz) may lower with high sympathetic gain to the heart. This is caused by the differences in the latencies/time constants in the sympathetic (17, 41) and vagal (6, 48) limbs, bringing the vagal and the sympathetic effects in counterphase in the BRS frequency band. Hence, there are situations thinkable in which cancellation of vagal effects on heart rate by concurring sympathetic effects on heart rate in counterphase incorrectly suggest a deficient baroreflex. For higher frequencies, the influence of the sympathetic feedback gain weakens and disappears because of a low pass filtering effect caused by slow neurotransmitter diffusion at the synaptic clefts (17).

In conclusion, our simulation study suggests that, within the limitations of the model, BRS and baroreflex-mediated blood pressure buffering are unrelated baroreflex features unless there is a good physiological reason to assume a fixed relation between the baroreflex feedback gains in the efferent baroreflex limbs to the heart and peripheral vasculature.

Also, we conclude that baroreflex-mediated blood pressure-buffering capacity is almost uniquely determined by the sympathetic baroreflex feedback gain to the peripheral vasculature, while BRS is almost uniquely determined by the vagal feedback gain to the sinus node.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This study was in part supported by Netherlands Heart Foundation Grant 2003B094.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank Prof. Karel H. Wesseling for critically reading this manuscript and Sum-Che Man for help in preparing the figures.


    FOOTNOTES
 

Address for reprint requests and other correspondence: C. A. Swenne, Dept. of Cardiology, Leiden Univ. Medical Center, PO Box 9600, 2300 RC Leiden, Netherlands (e-mail: c.a.swenne{at}lumc.nl)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 

  1. Abrahamsson C, Ahlund C, Nordlander M, Lind L. A method for heart rate-corrected estimation of baroreflex sensitivity. J Hypertens 21: 2133–2140, 2003.[CrossRef][ISI][Medline]
  2. Aljuri N, Marini R, Cohen RJ. Test of dynamic closed-loop baroreflex and autoregulatory control of total peripheral resistance in intact and conscious sheep. Am J Physiol Heart Circ Physiol 287: H2274–H2286, 2004.[Abstract/Free Full Text]
  3. Baller D, Bretschneider HJ, Hellige G. A critical look at currently used indirect indices of myocardial oxygen consumption. Basic Res Cardiol 76: 163–181, 1981.[CrossRef][ISI][Medline]
  4. Berger RD, Saul JP, Cohen RJ. Transfer function analysis of autonomic regulation. I. Canine atrial rate response. Am J Physiol Heart Circ Physiol 256: H142–H152, 1989.[Abstract/Free Full Text]
  5. Bertram D, Barres C, Cuisinaud G, Julien C. The arterial baroreceptor reflex of the rat exhibits positive feedback properties at the frequency of Mayer waves. J Physiol 513: 251–261, 1998.[Abstract/Free Full Text]
  6. Borst C, Karemaker JM. Time delays in the human baroreceptor reflex. J Auton Nerv Syst 9: 399–409, 1983.[CrossRef][ISI][Medline]
  7. Buch AN, Coote JH, Townend JN. Mortality, cardiac vagal control and physical training—what's the link? Exp Physiol 87: 423–435, 2002.[Abstract]
  8. Burgess DE, Hundley JC, Li SG, Randall DC, Brown DR. First-order differential-delay equation for the baroreflex predicts the 0.4-Hz blood pressure rhythm in rats. Am J Physiol Regul Integr Comp Physiol 273: R1878–R1884, 1997.[Abstract/Free Full Text]
  9. Cevese A, Gulli G, Polati E, Gottin L, Grasso R. Baroreflex and oscillation of heart period at 0.1 Hz studied by alpha-blockade and cross-spectral analysis in healthy humans. J Physiol 531: 235–244, 2001.[Abstract/Free Full Text]
  10. Chapuis B, Vidal-Petiot E, Orea V, Barres C, Julien C. Linear modelling analysis of baroreflex control of arterial pressure variability in rats. J Physiol 559: 639–649, 2004.[Abstract/Free Full Text]
  11. Christou DD, Jones PP, Seals DR. Baroreflex buffering in sedentary and endurance exercise-trained healthy men. Hypertension 41: 1219–1222, 2003.[Abstract/Free Full Text]
  12. De Boer RW, Karemaker JM, Strackee J. Hemodynamic fluctuations and baroreflex sensitivity in humans: a beat-to-beat model. Am J Physiol Heart Circ Physiol 253: H680–H689, 1987.[Abstract/Free Full Text]
  13. Eckberg DL. The human respiratory gate. J Physiol 548.2: 339–352, 2003.
  14. Eckberg DL, Sleight P. Human Baroreflexes in Health and Disease. Oxford, UK: Clarendon Press, 1992.
  15. Frederiks J, Swenne CA, TenVoorde BJ, Honzíková N, Levert JV, Maan AC, Schalij MJ, Bruschke AVG. The importance of high-frequency paced breathing in spectral baroreflex sensitivity assessment. J Hypertens 11: 1635–1644, 2000.
  16. Hammer PE, Saul JP. Resonance in a mathematical model of baroreflex control: arterial blood pressure waves accompanying postural stress. Am J Physiol Regul Integr Comp Physiol 288: R1637–R1648, 2005.[Abstract/Free Full Text]
  17. Hill-Smith I, Purves RD. Synaptic delay in the heart: an iontophoretic study. J Physiol 279: 31–54, 1978.[Abstract/Free Full Text]
  18. Hyndman BW, Mohn RK. A model of the cardiac pacemaker and its use in decoding the information content of cardiac intervals. Automedica 1: 239–252, 1973.
  19. Japundzic N, Grichois ML, Zitoun P, Laude D, Elghozi JL. Spectral analysis of blood pressure and heart rate in conscious rats: effects of autonomic blockers. J Auton Nerv Syst 30: 91–100, 1990.[CrossRef][ISI][Medline]
  20. Jones PP, Christou DD, Jordan J, Seals DR. Baroreflex buffering is reduced with age in healthy men. Circulation 107: 1770–1774, 2003.[Abstract/Free Full Text]
  21. Jordan J, Tank J, Shannon JR, Diedrich A, Lipp A, Schroder C, Arnold G, Sharma AM, Biaggioni I, Robertson D, Luft FC. Baroreflex buffering and susceptibility to vasoactive drugs. Circulation 105: 1459–1464, 2002.[Abstract/Free Full Text]
  22. Julien C. The enigma of Mayer waves: Facts and models. Cardiovasc Res 70: 12–21, 2006.[Abstract/Free Full Text]
  23. Just A, Wittmann U, Nafz B, Wagner CD, Ehmke H, Kirchheim HR, Persson PB. The blood pressure buffering capacity of nitric oxide by comparison to the baroreceptor reflex. Am J Physiol Heart Circ Physiol 267: H521–H527, 1994.[Abstract/Free Full Text]
  24. Karemaker JM. Neurophysiology of the baroreceptor reflex. In: The Beat-to-Beat Investigation of Cardiovascular Function, edited by Kitney R and Rompelman O. Oxford: Clarendon Press, 1987, p. 27–49.
  25. Katona PG, McLean M, Dighton DH, Guz A. Sympathetic and parasympathetic cardiac control in athletes and nonathletes at rest. J Appl Physiol Respir Environ Exercise Physiol 52: 1652–1657, 1982.[Abstract/Free Full Text]
  26. Katona PG, Poitras JW, Barnett GO, Terry BS. Cardiac vagal efferent activity and heart period in the carotid sinus reflex. Am J Physiol 218: 1030–1037, 1970.[Free Full Text]
  27. La Rovere MT, Bersano C, Gnemmi M, Specchia G, Schwartz PJ. Exercise-induced increase in baroreflex sensitivity predicts improved prognosis after myocardial infarction. Circulation 106: 945–949, 2002.[Abstract/Free Full Text]
  28. La Rovere MT, Bigger JT Jr, Marcus FI, Mortara A, Schwartz PJ. Baroreflex sensitivity and heart-rate variability in prediction of total cardiac mortality after myocardial infarction. Lancet 351: 478–484, 1998.[CrossRef][ISI][Medline]
  29. Laurent D, Bolene-Williams C, Williams FL, Katz LN. Effects of heart rate on coronary flow and cardiac oxygen consumption. Am J Physiol 185: 355–364, 1956.[Abstract/Free Full Text]
  30. Liu HK, Guild SJ, Ringwood JV, Barrett CJ, Leonard BL, Nguang SK, Navakatikyan MA, Malpas SC. Dynamic baroreflex control of blood pressure: influence of the heart vs. peripheral resistance. Am J Physiol Regul Integr Comp Physiol 283: R533–R542, 2002.[Abstract/Free Full Text]
  31. Madwed JB, Albrecht P, Mark RG, Cohen RJ. Low-frequency oscillations in arterial pressure and heart rate: a simple computer model. Am J Physiol Heart Circ Physiol 256: H1573–H1579, 1989.[Abstract/Free Full Text]
  32. Malpas SC. Neural influences on cardiovascular variability: possibilities and pitfalls. Am J Physiol Heart Circ Physiol 282: H6–H20, 2002.[Abstract/Free Full Text]
  33. Mayer S. Studien zur Physiologie des Herzens und der Blutgefässe. V. über spontane Blutdruckschwankungen. Akad Wiss Math-Nat Kl (Wien) 74: 281–307, 1876.
  34. Parati G, Di Rienzo M, Mancia G. How to measure baroreflex sensitivity: from the cardiovascular laboratory to daily life. J Hypertens 18: 7–19, 2000.[ISI][Medline]
  35. Penáz J, Buriánek P. Zeitverlauf und Dynamik der durch Atmung ausgelösten Kreislaufänderungen beim Menschen. Pflügers Arch 267: 618–635, 1962.
  36. Pickering TG, Gribbin B, Sleight P. Comparison of the reflex heart rate response to rising and falling arterial pressure in man. Circ Res 6: 277–283, 1972.
  37. Pinna GD, Maestri R. Reliability of transfer function estimates in cardiovascular variability analysis. Med Biol Eng Comput 39: 338–347, 2001.[CrossRef][ISI][Medline]
  38. Pomeranz B, Macaulay RJ, Caudill MA, Kutz I, Adam D, Gordon D, Kilborn KM, Barger AC, Shannon DC, Cohen RJ. Assessment of autonomic function in humans by heart rate spectral analysis. Am J Physiol Heart Circ Physiol 248: H151–H153, 1985.[Abstract/Free Full Text]
  39. Robbe HWJ, Mulder LJM, Rüddel H, Langewitz WA, Veldman JBP, Mulder G. Assessment of baroreceptor reflex sensitivity by means of spectral analysis. Hypertension 10: 538–543, 1987.[Abstract/Free Full Text]
  40. Rosenblueth A, Simeone FA. The interrelations of vagal and accelerator effects on the cardiac rate. Am J Physiol 110: 42–45, 1934.[Free Full Text]
  41. Saul JP, Berger RD, Albrecht P, Stein SP, Chen MH, Cohen RJ. Transfer function analysis of the circulation: unique insights into cardiovascular regulation. Am J Physiol Heart Circ Physiol 261: H1231–H1245, 1991.[Abstract/Free Full Text]
  42. Scheffer GJ, TenVoorde BJ, Karemaker JM, Ros HH. Effects of epidural analgesia and atropine on heart rate and blood pressure variability: implications for the interpretation of beat-to-beat fluctuations. Eur J Anaesthesiol 11: 75–80, 1994.[ISI][Medline]
  43. Simms AE, Paton JF, Pickering AE. Disinhibition of the cardiac limb of the arterial baroreflex in rat: a role for metabotropic glutamate receptors in the nucleus tractus solitarii. J Physiol 575: 727–738, 2006.
  44. Smyth HS, Sleight P, Pickering GW. Reflex regulation of arterial pressure during sleep in man: a quantitative method of assessing baroreflex sensitivity. Circ Res 24: 109–121, 1969.[Abstract/Free Full Text]
  45. Taylor JA, Eckberg DL. Fundamental relations between short-term RR interval and arterial pressure oscillations in humans. Circulation 93: 1527–1532, 1996.[Abstract/Free Full Text]
  46. TenVoorde BJ, Kingma R. A baroreflex model of short term blood pressure and heart rate variability. Stud Health Technol Inform 71: 179–200, 2000.[Medline]
  47. Thrasher TN. Baroreceptors and the long-term control of blood pressure. Exp Physiol 89: 331–335, 2004.[Abstract/Free Full Text]
  48. Warner HR, Cox A. A mathematical model of heart rate control by sympathetic and vagus efferent information. J Appl Physiol 17: 349–355, 1962.[Free Full Text]
  49. Wesseling KH, Settels JJ. Baromodulation explains short-term blood pressure variability. In: Psychophysiology of Cardiovascular Control, edited by Orlebeke JF, Mulder G, and VanDoornen LJP. New York: Plenum, 1985, p. 69–97.
  50. Xavier-Neto J, Moreira E, Krieger E. Viscoelastic mechanisms of aortic baroreceptor resetting to hypotension and hypertension. Am J Physiol Heart Circ Physiol 271: H1407–H1415, 1996.[Abstract/Free Full Text]
  51. Xie P, McDowell T, Chapleau M, Hajduczok G, Abboud F. Rapid baroreceptor resetting in chronic hypertension. Implications for normalization of arterial pressure. Hypertension 17: 72–79, 1991.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
R. Kanbar, V. Orea, B. Chapuis, C. Barres, and C. Julien
A transfer function method for the continuous assessment of baroreflex control of renal sympathetic nerve activity in rats
Am J Physiol Regulatory Integrative Comp Physiol, November 1, 2007; 293(5): R1938 - R1946.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
G. D. Pinna
Assessing baroreflex sensitivity by the transfer function method: what are we really measuring?
J Appl Physiol, April 1, 2007; 102(4): 1310 - 1311.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
102/4/1348    most recent
00158.2006v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by van de Vooren, H.
Right arrow Articles by Van der Wall, E. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by van de Vooren, H.
Right arrow Articles by Van der Wall, E. E.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2007 by the American Physiological Society.