Journal of Applied Physiology Watch the video to learn how APS reaches out to developing nations.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 101: 92-101, 2006. First published March 30, 2006; doi:10.1152/japplphysiol.01305.2005
8750-7587/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
101/1/92    most recent
01305.2005v1
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 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 Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Shmuylovich, L.
Right arrow Articles by Kovács, S. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Shmuylovich, L.
Right arrow Articles by Kovács, S. J.

Load-independent index of diastolic filling: model-based derivation with in vivo validation in control and diastolic dysfunction subjects

Leonid Shmuylovich and Sándor J. Kovács

Cardiovascular Biophysics Laboratory, Washington University School of Medicine, St. Louis, Missouri

Submitted 12 October 2005 ; accepted in final form 20 March 2006


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Maximum elastance is an experimentally validated, load-independent systolic function index stemming from the time-varying elastance paradigm that decoupled extrinsic load from (intrinsic) contractility. Although Doppler echocardiography is the preferred method of diastolic function (DF) assessment, all echo-derived indexes are load dependent, and no invasive or noninvasive load-independent index of filling (LIIF) exists. In this study, we derived and experimentally validated a LIIF. We used a kinematic filling paradigm (the parameterized diastolic filling formalism) to predict and derive the (dimensionless) dynamic diastolic efficiency M, defined by the slope of the peak driving force [maximum driving force (kxo) {propto} peak atrioventricular (AV) gradient] to maximum viscoelastic resistive force [peak resistive force (cEpeak)] relation. To validate load independence, we analyzed E-waves recorded while load was varied via tilt table (head up, horizontal, and head down) in 16 healthy volunteers. For the group, linear regression of E-wave derived kxo vs. cEpeak yielded kxo = M (cEpeak) + B, r2 = 0.98; where M = 1.27 ± 0.09 and B = 5.69 ± 1.70. Effects of diastolic dysfunction (DD) on M were assessed by analysis of preexisting simultaneous cath-echo data in six DD vs. five control subjects. Average M for the DD group (M = 0.98 ± 0.07) was significantly lower than controls (M = 1.17 ± 0.05, P < 0.001). We conclude that M is a LIIF because it uncouples intrinsic DF (i.e., the pressure-flow relation) from extrinsic load (left ventricular end-diastolic pressure). Larger M values imply better DF in that increasing AV pressure gradient results in relatively smaller increases in peak resistive losses (cEpeak). Conversely, lower M implies that increasing AV gradient leads to larger increases in resistive losses. Further prospective validation characterizing M in well-defined pathological states is warranted.

diastolic function; Doppler echocardiography; mathematical modeling


THE ANALYSIS OF PRESSURE-VOLUME data obtained from the in vivo heart for the purpose of determining intrinsic pump function is complex, because pressures and volumes are dependent on both intrinsic variables (such as contractility and muscle stiffness) and extrinsic variables (such as preload and afterload). By viewing cardiac function mechanically, as a system having time-varying elastance, Suga et al. (39) successfully uncoupled these variables. They showed experimentally that the instantaneous pressure-to-volume ratio defines a time-varying elastance that attains the same maximum value at a fixed contractile state, regardless of preload changes. This established maximum elastance (Emax) as an experimentally validated, load-independent index of systolic function (36). Additional conceptual validation of Emax as a load-independent index has been achieved using a kinematic, forced harmonic oscillator-based argument showing that the slope of the maximum force-displacement relationship (Emax analog) depends only on the intrinsic oscillator parameters rather than the initial (load) conditions (29). While Emax is a chamber property that is uncoupled from the effects of load on systolic function, no noninvasive, load-independent attribute for diastole has been theoretically proposed and validated, or empirically observed in the course of experiment.

Doppler echocardiography is the preferred method for noninvasive diastolic function (DF) assessment. Doppler-derived indexes have been used to characterize DF in numerous cardiac disorders, including heart failure, myocardial infarction, hypertrophic cardiomyopathy, and hypertension (1, 37). In current practice, most DF indexes are derived by visual inspection of transmitral E- and A-waves. These shape-derived indexes include peak velocity of the E-wave (Epeak), duration of the E-wave, acceleration (AT) and deceleration (DT) times of the E-wave, and area under the E-wave [velocity-time integral (VTI)]. Additional indexes include the peak velocity of the A wave (Apeak) and the ratio of the Epeak and Apeak velocities (E/A). Most clinically relevant Doppler-derived DF indexes have proved to be load dependent in animals and humans, both in health and in disease (3, 8, 11, 14, 17, 26, 30, 32, 33, 38, 4042, 45, 48). Several newer load-independent indexes have been empirically suggested, from velocity of propagation (Vp), to annular velocities derived from Doppler tissue imaging (E'/A', E', E/E'), but there is no consensus as to whether these indexes are truly load independent (9, 15, 24, 27, 28, 31, 47).

To quantify DF, we modeled filling in kinematic terms via the parameterized diastolic filling (PDF) formalism (1820). The formalism is a lumped parameter, predictive rather than accommodative (23), model that characterizes transmitral flow according to the motion of a damped simple harmonic oscillator (SHO). This model characterizes transmitral blood flow velocity in terms of elastic, inertial, and damping forces. During filling, the elastic driving force generates both inertial forces, causing acceleration, and resistive (damping) forces, opposing acceleration. The three (mathematically) independent model parameters, spring constant k, damping constant c, and initial spring displacement xo, fully characterize the velocity of the SHO (i.e., E-wave velocity contour). These parameters are determined by solving the "inverse" problem of filling, using the clinical E-wave contour as the input, and the model parameters as the best-fit determined output (12). The model-predicted velocity provides an excellent fit to the Doppler E-wave contour recorded clinically (12). In accordance with first principles-based predictions and causality, the model parameters have physiological analogs that have been experimentally validated in vivo (Table 1). Specifically, the chamber stiffness (dP/dV) is linearly related to k (21); the initial peak force (kxo) that drives the oscillator is the analog of the peak atrioventricular (AV) pressure gradient generating transmitral blood flow (2); 1/2kxFormula is the energy (ergs) available before valve opening (20); and xo is linearly related to the volumetric load, i.e., the VTI of the E-wave (20). In addition, the PDF formalism has been extensively tested and validated in subjects with a wide range of cardiac pathologies and loads, including hypertension (22), heart failure (34), and diabetes (35).


View this table:
[in this window]
[in a new window]
 
Table 1. Physiological analogs of PDF parameters and parameter-derived indexes

 
The load-independent index of diastolic filling problem has remained theoretically and experimentally unsolved (5). To tackle the problem theoretically, we used the kinematic paradigm of filling to derive a dimensionless ratio, with the previously invoked requirement (29) (showing that Emax is indeed load independent) that the algebraic expression for the load-independent index of filling (LIIF) must be independent of load in the terms of the model, i.e., independent of initial spring displacement xo. As noted in Table 1, mathematical expressions for the peak AV gradient and the resistive forces opposing blood flow at peak velocity can be derived from the E-wave contour using the PDF model. Mathematical analysis reveals that these expressions are linearly related, regardless of load. Therefore, a simple load-independent expression having clear physiological meaning relates maximum driving force (analog of peak AV gradient, kxo) to resistive (viscoelastic) force at peak flow (cEpeak). Although superficially similar in appearance, the relationship of kxo to cEpeak differs importantly from previous work by Hasegawa et al. (13), who compared Epeak to peak AV gradient (Epeak vs. kxo using our terminology) at different loading conditions. First, unlike the slope of their Epeak-to-peak AV gradient relationship, the slope of our peak AV gradient-to-peak resistive force relationship is dimensionless. Second and more importantly, our relation accounts for the resistive (viscoelastic) component of the filling process, whereas the Epeak-to-peak AV gradient relation has no explicit term accounting for resistive forces. The kxo (peak AV gradient kxo)-to-cEpeak relationship presented in this study can, therefore, be viewed as complementary to, but independent of, the prior work of Hasegawa et al.

To test the load independence of the predicted ratio, we altered load in healthy human volunteer subjects by changing tilt table position while recording transmitral flow. We investigated changes in traditional echo-derived indexes, changes in PDF parameters, and the load dependence of the predicted dimensionless ratio. Results of other studies have unambiguously shown that, relative to horizontal, preload increases with head-down tilt and decreases with head-up tilt (7, 11). In addition, many studies have altered preload by performing tilt table position changes (8, 11, 17, 30, 32, 33, 40, 48).

To assess the ability of slope M to differentiate between control and diastolic dysfunction (DD) states, we analyzed preexisting data from patients undergoing simultaneous catheterization-echocardiography using micromanometric (Millar) catheters. We calculated the kxo-to-cEpeak relationship from E-waves as a function of load, which was defined by beat-to-beat left ventricular (LV) end-diastolic pressure (LVEDP) respiratory variation. We determined the linearity of the kxo-to-cEpeak relationship as a function of LVEDP in this data set, as well as the ability of the slope (M) of the linear regression relationship to differentiate between DD and control states.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Theory

Modeling physiology: decoupling the kinematics of filling from load.   E-wave velocity contours are known to be load dependent, and, because the PDF formalism provides an excellent fit to all observed E-waves, the PDF parameters will vary as E-waves vary in response to changes in load.

To derive a load-independent index, we sought to decouple load (xo) from filling by employing our kinematic model of filling (the PDF model). Specifically, we note that, regardless of how E-wave contours respond to variations in load, the equation of motion (Newton's law), which includes inertial, damping, and recoil forces for the SHO, remains valid. The equation that approximates the E-wave is:

Formula 1(1)
where m, c, and k denote inertia, damping, and spring constants, respectively. This relationship applies at all times, including at time t = tpeak when the peak of the E-wave is inscribed. The acceleration term d2x/dt2 is zero at t = tpeak. We change notation so that dx/dt (velocity of the damped spring) is expressed as Epeak, and Eq. 1 becomes:

Formula 2(2)
Because the peak velocity is attained quickly, the instantaneous driving force exerted by the spring at peak velocity, kx (tpeak), can be linearly approximated by the initial, peak driving force kxo, the model's analog of the peak AV gradient. We can approximate this linear relationship as:

Formula 3(3)
and therefore approximate and rewrite Eq. 2 as

Formula 4(4)
where kxo is the mathematical analog for the maximum AV gradient and cEpeak represents the resistive (viscoelastic) force opposing filling. Hence we predict that the relation (slope) between maximum driving force and peak viscous force (kxo vs. cEpeak) should be linear and should be load independent. This predicted linearity between kxo and cEpeak should be load independent because it is derived from equations and approximations that are load independent.

Experiment

Data acquisition.   Pulsed Doppler echocardiography was used to acquire continuous transmitral blood flow data from 16 subjects (9 men, 7 women, ages 20–30 yr), while subjects were positioned on a tilt table (Trex Medical, Danbury, CT). Doppler data were acquired in the apical four-chamber views with the sample volume gated at 1.5–2.5 mm and directed between the tips of the mitral valve leaflets orthogonal to the mitral valve plane. The subjects were healthy medical and graduate students on no prescribed medications. None of the subjects had a history of heart disease, coronary artery disease, hypertension, or diabetes. Before participation in the study, all subjects gave informed consent for participation in accordance with a protocol reviewed and approved by the Washington University Medical Center Human Studies Committee (IRB).

Doppler data were obtained with a clinical echocardiographic imaging system (Acuson Sequia 256, Mountain View, CA) equipped with a 2-MHz transducer. Heart rate was recorded simultaneously via ECG limb lead II and displayed on the E- and A-wave images. Blood pressure was monitored via a digital blood pressure cuff (Medtronic LifePak 12, Minneapolis, MN).

The data acquisition protocol consisted of initial (baseline) E- and A-wave recording with the subject supine and the tilt table in the horizontal position for 5 min, such that the heart rate was in steady state. Baseline blood pressure was also obtained. After 5 min in the horizontal position, the table, including appropriate padding and straps to ensure safety, was gradually tilted to a 90° head-down position. Since the heart/diaphragm shifted during tilt, the transducer was suitably reoriented to obtain transmitral flow. After resolution of heart rate transients, typically after several minutes, transmitral flow was recorded. After several minutes, the tilt table was returned to the horizontal position. Once heart rate and blood pressure returned to baseline levels, transmitral flow was again recorded. After 5 min horizontally, the tilt table was moved to 90° head-up position. The transducer was adjusted to account for the heart shift during tilt, and, after heart rate and blood pressure transients resolved, transmitral flow was recorded. In 2 of the 16 subjects, the persistent increase in heart rate during head-up tilt resulted in significant E- and A-wave merging. With significant merging, the E-wave deceleration portion is lost, and thus it is not possible to reliably fit merged waves via the PDF formalism. Hence for these two subjects, head-up E-waves were not analyzed, but the head-down and horizontal data were included in the overall analysis. After 5 min upright, subjects were tilted to the horizontal position, and a concluding set of transmitral flow images was recorded. All data were recorded on VHS tape for off-line analysis using a custom-editing station.

Data analysis.   For each subject at each stage of tilt, five beats were captured and analyzed by model-based image processing (MBIP). The MBIP method by which PDF parameters are obtained from digitized recordings of transmitral flow has been previously described (1921). Briefly, transmitral waves for individual beats were analyzed from still images of transmitral velocity profiles, and conventional Doppler-derived indexes were calculated. In addition, E- and A-wave contours were fit via the PDF formalism, which permits the contour of the E-wave to be expressed in a closed-form mathematical expression. Minimizing the error between the model-predicted contour and the clinical E-wave contour via the Levenberg-Marquardt algorithm provides best-fit (unique) model parameters (c, k, xo) for each E-wave. Since directly measured and model calculated values of Epeak and AT show nearly perfect agreement, Epeak was calculated directly from the model-predicted contour for the E-wave.

Additional validation in pathological cases.   In addition to testing the LIIF noninvasively in healthy normal subjects via tilt table, we tested the load dependence of M in subjects with and without DD undergoing diagnostic catheterization. This additional analysis used existing data from previous studies (2, 25), utilizing the Cardiovascular Biophysics Laboratory database of simultaneous Doppler echocardiographic transmitral flow and micromanometric (Millar) catheter-derived intraventricular (LV) pressure. We note that the original intent of the Cardiovascular Biophysics Laboratory database was not explicit testing of load dependence. However, because some subjects manifested significant beat-to-beat respiratory variation of LVEDP (load) and simultaneously recorded E-wave contours, the data were suitable for determination of the kxo-to-cEpeak relationship as a function of load, i.e., varying LVEDP. We, therefore, selected subjects who had good quality E-waves as well as significant load variation (LVEDP variation > 10 mmHg) in response to the respiratory cycle. Selection criteria for inclusion in the DD group required normal sinus rhythm, no evidence of valvular disease, no active ischemia, normal ejection fraction (EF) > 60%, and elevated LVEDP (>19 mmHg). Subjects in the control group had normal sinus rhythm, no valvular disease, no active ischemia, normal EF > 60%, and normal LVEDP. Because subjects are referred for catheterization to establish the presence of coronary artery disease, variable degrees of coronary artery disease were present in both the control and DD groups. However, no subject in either group had ongoing or active ischemia. Demographics are presented in Table 2. For each subject, 25 consecutive E-waves and simultaneous LV pressures were analyzed. For each subject, good quality E-waves were selected and analyzed using the MBIP procedure described above. LVEDP was determined from the simultaneous LV pressure data.


View this table:
[in this window]
[in a new window]
 
Table 2. Demographics of simultaneous catheterization-echocardiography subjects analyzed retrospectively

 

    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Variation of E-wave Shape With Tilt Table Position

In agreement with other studies (8, 11, 18, 30, 32, 33, 40), both E- and A-wave shapes varied in response to changes in load generated by changes in tilt table position. Because the PDF parameters are determined from the contour of the waves, changes in load resulted in changes in the PDF parameters. Figure 1 shows representative Doppler waves from a subject in the head-up, horizontal, and head-down position, as well as the PDF model-predicted fit to the E-wave contours. Table 3 summarizes the average values for both traditional Doppler indexes and PDF parameters under different load states i.e., tilt table positions.


Figure 1
View larger version (94K):
[in this window]
[in a new window]
 
Fig. 1. Pulsed wave transmitral flow-velocity images from 1 subject at 3 different preload (tilt table positions) states. Single diastolic interval at each tilt table position is shown. Parameterized diastolic filling (PDF) model-predicted fit to each E-wave is shown in bottom panels.

 

View this table:
[in this window]
[in a new window]
 
Table 3. Average values of Doppler indexes and PDF parameters in one typical tilt table subject at three different preload states

 
Since subjects had a distribution of values for Doppler indexes and PDF parameter values in the horizontal tilt table position, determining an average over all subjects required reference to a common baseline. Horizontal was chosen as the baseline tilt table position, and all values were calculated as percent change from horizontal. Table 4 presents both PDF and traditional Doppler data as percent change from horizontal for both head-up and head-down tilt.


View this table:
[in this window]
[in a new window]
 
Table 4. Averaged, normalized percent increase or decrease relative to baseline (horizontal tilt table position) for all tilt table subjects

 
Determination of Predicted Load Independence

In accordance with the prediction (see METHODS) that the maximum driving force must be linearly related to peak viscous (resistive) force, a linear regression via least mean square error of cEpeak to kxo was performed. Figure 2A shows a representative kxo vs. cEpeak plot for one subject. Head-up data from two subjects were not included, but their head-down and horizontal data fit the kxo vs. cEpeak regression with high r2. Among all subjects, the average slope of the kxo vs. cEpeak plot was 1.27 (SD 0.09), and the average intercept was 5.69 (SD 1.70). The average r2 value for each subjects' linear regression was r2 = 0.95. The combined data from all subjects plotted together yielded a linear regression with r2 = 0.98 (Fig. 2B).


Figure 2
View larger version (18K):
[in this window]
[in a new window]
 
Fig. 2. A: maximum driving force [kxo, peak atrioventricular (AV) gradient] vs. peak resistive force (cEpeak) for 1 subject at 3 different preload states. Note slope of best linear fit is independent of tilt table position. B: kxo vs. cEpeak for all (n = 16) subjects at different preload states. Reported values represent 5-beat average for kxo and cEpeak for each subject at each preload state. See text for details.

 
Analysis of the kxo vs. cEpeak Relation: The Role of Randomness

Our experimental tilt table-derived results show the existence of a load-independent linear relationship between kxo and cEpeak. One issue is whether the damped-SHO model itself requires that the kxo vs. cEpeak relation must always be linear. To address this issue, we selected random values of k, xo, and c to create E-waves, with the requirement that their amplitudes and durations fall within the range of values that has been previously observed and used to fit Doppler data via the PDF formalism. If our observations were a tautology, these E-waves would show a linear relationship between k, xo, and c. The graph of kxo vs. cEpeak using such random values (within the physiological range) is shown in Fig. 3.


Figure 3
View larger version (23K):
[in this window]
[in a new window]
 
Fig. 3. A: 3 typical model-generated E-waves created by randomly picking values for initial spring displacement xo, damping constant c, and spring constant k, known to be in the physiological range. B: kxo vs. cEpeak for the 3 random E-waves shown in A. Note deterioration of r2. C: increase in randomly generated E-wave sample size to n = 10 indicates further, substantial deterioration (r2 = 0.01) of the observed, highly linear, kxo-to-cEpeak relationship. See text for details.

 
Figure 3 shows that not all values of k, xo, and c generate a linear relationship between kxo and cEpeak. Therefore, the observed linear relation between kxo and cEpeak is not a general property of a damped spring model assembled from arbitrary components. Furthermore, the range of values of k, xo, and c that satisfy the kxo vs. cEpeak relationship for a particular heart represents the physiological regime in which that heart's filling attributes reside.

It should be noted that k, c, and xo are, from a mathematical standpoint, independent parameters. There is therefore no a priori correlation between k, c, xo, or any combination of the parameters that can be predicted. These parameters become physiologically coupled, and their magnitudes are constrained, once they are fit to clinical data (the E-wave). Once determined by fitting to actual E-waves, several correlations can be seen, such as Epeak to kxo, k to cEpeak, and c to kxo. However, kxo to cEpeak is the strongest observed correlation and is the only one with a dimensionless slope amenable to a simple physiological interpretation.

Effect of DD on M: Results From Available Data

Using existing simultaneous cath-echo data, we analyzed the kxo vs. cEpeak relationship in subjects with significant respiratory (LVEDP) variation. The results, including statistics, are summarized in Fig. 4 and Table 5. The kxo vs. cEpeak relation was highly linear for both normal (average r2 = 0.96 ± 0.02) and DD groups (average r2 = 0.90 ± 0.05). The average slope for the normal group was M = 1.17 ± 0.05, and the average value for the DD group was M = 0.98 ± 0.07 (P < 0.001 by ANOVA). Additionally, the slope intercept of the normal cath-echo group was B = 6.69 ± 0.91, whereas the slope intercept of the DD cath-echo group was B = 10.67 ± 3.35 (P = 0.03 by ANOVA). The difference in M between the cath-echo DD subjects and the healthy tilt table subjects showed stronger statistical significance (P < 0.00001) compared with the difference in M between the cath-echo normal subjects and the healthy tilt table subjects (P = 0.02). The difference in intercept B is also statistically significant between DD subjects and tilt table subjects (P < 0.0001). There was no statistically significant difference, however, between the intercept of the normal cath-echo subjects and the tilt table subjects.


Figure 4
View larger version (9K):
[in this window]
[in a new window]
 
Fig. 4. Summary comparison of slope M, intercept B, and left ventricular end-diastolic pressure (LVEDP) for all (n = 27) subjects comprising the tilt table group (n = 16) and the catheterization-echocardiography group (n = 11). A: thick lines denote average values, thin lines above and below denote 1 SD relative to the mean LVEDP. B: regression slope (M) comparison between groups. Thick line is the average value from the tilt table study. Dotted lines are 1 SD relative to the mean value. C: intercept (B) comparison between groups. Thick line is average value from the tilt table study. Dotted lines are 1 SD relative to the mean intercept value. See text for details.

 

View this table:
[in this window]
[in a new window]
 
Table 5. Comparison of the slope (M), intercept (B), and LVEDP for all (n = 27) subjects comprised of the tilt table group (n = 16) and the catheterization-echocardiography group (n = 11)

 

    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Background on Doppler Echo

A major challenge facing interpretation of echo-derived DF indexes is that the indexes in current clinical use (Epeak, AT, DT, E/A, VTI) are load dependent (3, 8, 11, 14, 17, 26, 32, 33, 4042, 44, 45, 48). Despite this uncertainty, it is often assumed that changes observed in the indexes represent changes in pathological processes. Studies have clearly established, however, that traditional transmitral flow-derived diastolic indexes are load dependent, and, therefore, changes in DF indexes could be due to load variability rather than to pathology.

Using blood pressure cuffs to noninvasively reduce preload in 12 normal subjects, Triulzi et al. (44), for example, found Epeak, E/A, VTI, DT, and AT to change significantly with preload and concluded that preload reduction in normal subjects produces a "pattern that mimics changes in left ventricular diastolic dysfunction." By altering preload with postural changes, Downes et al. (8) found that, although E/A changed, it did not change enough for normal values to migrate into the abnormal range and concluded that "simple changes in venous return do not ‘normalize’ an abnormal pattern, nor do they ‘abnormalize’ a normal pattern." Kmetzo et al. (17), however, found that 80° head-up tilt produced abnormal relaxation patterns in 22 normal heart-healthy volunteers.

Newer load-independent Doppler-derived indexes of DF have been proposed (Vp, E', E/E'), but all have been shown to have load dependence. E' has been noted to be load independent in patients with chronic ischemic syndrome (47) but has shown variability with changing preload in normal subjects (30, 31), in animals (15), and in patients undergoing hemodialysis (28). Vp has shown load independence and strong correlation to the time constant of isovolumic relaxation in patients undergoing cardiac surgery (9) but proven to be load dependent in normal volunteers (30). Even E/E', which many studies find to be load independent (25, 28, 30), has proved in animal experiments to be preload dependent (15).

Despite studies that have shown DF indexes to be load dependent, no studies have derived diastolic indexes that compensate for load. Recent studies (13) have shown, however, that intrinsic properties of the ventricle may correlate better with load-induced changes in the Doppler indexes rather than the Doppler indexes themselves. For example, Tanabe et al. (43) found that the change in E/A following a decrease in preload correlated strongly with an invasively derived intrinsic property of the ventricle.

In contrast, Hasegawa et al. (13) observed Epeak (ml/s) volumetric flow to be linearly related to the peak AV pressure gradient (r2 = 0.94) in accordance with the Bernoulli relation relating pressure gradient to flow velocity. Although load was not explicitly altered, the observed change in peak AV gradient is consistent with a variable load. Interestingly, the slope of the observed relation is highly linear and is nearly indistinguishable between the control and heart failure groups. Because the slope of the Epeak vs. peak AV gradient relation discussed by Hasegawa et al. remains linear and does not change in the setting of heart failure, although load may have changed, in our view it reflects the applicability of the Bernoulli equation rather than fulfilling the criterion for a LIIF.

Applying the PDF Model

The PDF model helps elucidate the connection between the experimentally determined Epeak vs. peak AV gradient relation by Hasegawa et al. (13) and the present study. In agreement with Hasegawa et al., a plot (not shown) of Epeak vs. kxo (the peak AV gradient equivalent) using our data also generates a linear relation. The slope has units of seconds per kilogram; the linear regression relation is Epeak = 0.017(kxo) + 0.39, r2 = 0.78. The fact that an approximately linear relation is observed is predictable in part from the nonsteady Bernoulli relation, which relates pressure gradient to flow velocity. Moreover, the Epeak vs. kxo relation does not explicitly account for the effect of damping/viscosity as can be achieved by inclusion of the parameter c.

A stronger and more meaningful correlation is achieved by the kxo (peak AV gradient kxo) vs. cEpeak relation, which attains r2 = 0.98 compared with the r2 = 0.78 for the Epeak vs. peak AV gradient (kxo) relation. In addition to attaining a stronger correlation, the slope of the kxo (peak AV gradient kxo) vs. cEpeak relation is dimensionless. It complements the Hasegawa et al. (13) observation by explicitly including the effect of damping. A high value for cEpeak implies a large cEpeak must be generated to achieve the observed flow. This higher resistive force indicates that conversion of (maximum AV gradient) potential energy into kinetic energy (flow) is less efficient. The PDF (viscoelastic) parameter c is important because it has been shown to be significantly different from normal in pathophysiological states with delayed relaxation patterns such as diabetes (6, 35). A higher value of the parameter c in diabetic (human and rat) hearts implies that diabetic hearts, compared with normal hearts, are less efficient in being able to convert potential energy (pressure gradient) into kinetic (flow) energy. Therefore, we expect that, unlike the Epeak vs. peak AV gradient (kxo) relationship, the kxo (peak AV gradient kxo) vs. cEpeak relationship will have quantifiable differences between normal and abnormal LV function states.

Although not the primary intent of the present study, Fig. 4 (derived from cath-echo data acquired for another purpose) provides compelling evidence that both slope M and intercept B, of the kxo vs. cEpeak linear regression, change in LV DD states. The DD group had lower M values and higher B values compared with the normal cath-echo group and the healthy tilt table group. The normal cath-echo group had the same B values as the healthy tilt table group and slightly lower M values. This difference in M between the normal cath-echo subjects and the tilt table subjects may be due to the higher average age of the normal cath-echo subjects (44 yr) relative to the healthy tilt table group (25 yr).

These preliminary data, particularly because they were not acquired with the specific intent of determining load independence, strongly support the view that M and B can differentiate between normal and pathological states. Because M is the slope of a regression relation, it provides relative information regarding the efficiency with which the ventricle adapts to load. In contrast, the intercept correlates with LVEDP, an absolute rather than relative index. New (cath-echo) studies, specifically including maneuvers that vary load, are indicated to more fully elucidate the relationship between M and B, and phenotypically well-characterized DD states.

Interestingly, only certain regimes of the kxo vs. cEpeak graph can exist physiologically. Energy conservation makes it impossible for cEpeak, the peak resistive force opposing filling, to exceed kxo, the initial (and maximum) driving force initiating filling. If M is low, specifically <1, as preload increases, the kxo vs. cEpeak regression line approaches the kxo = cEpeak boundary, eventually crossing it and moving into the nonphysiological regime. A ventricle having M < 1 would, therefore, not be able to maintain filling function in the face of increased preload (Fig. 5).


Figure 5
View larger version (29K):
[in this window]
[in a new window]
 
Fig. 5. Equation governing filling (Eq. 1) constrains the physiology to the upper half of the plot because the peak resistive force can never exceed the maximum driving force. A chamber with M < 1 operates on a regression line that may eventually reach the prohibited regime for sufficiently elevated peak driving force (AV gradient) values. In contrast, a chamber having M ≥ 1 is not similarly constrained. See text for details.

 
Limitations

A minor limitation is the absence of end-diastolic volume (EDV) data as a correlate of preload alteration at different stages of tilt. EDV could not be reliably calculated, because only a four-chamber view was archived. Transthoracic echocardiographic studies have shown that EDV calculated from four-chamber views alone achieve only an r =0.61 correlation with true EDV (46). Although exact preload change with tilt alteration is not reported, many studies have altered preload using tilt table methodology (8, 11, 17, 30, 32, 33, 40, 47), and it is certain that preload is highest in head-down tilt and lowest in head-up tilt (7, 11).

A second limitation is that the PDF formalism is most applicable to E-waves having both ascending and descending portions. E- and A-waves become difficult to separate and discern when the A-wave merges with the E-wave and covers more than two-thirds of the E-wave deceleration (4), which typically occurs at heart rate >100 beats/min. Thus M values for subjects with high heart rates (>100 beats/min) were not computed. In the present study, only 2 of 16 subjects during head-up tilt had limited E-wave deceleration portions. Hence head-up data from these two subjects were not fit by the PDF formalism.

Finally, in our tilt table study, all data were acquired from heart-healthy volunteers, and so a detailed comparison of the kxo vs. cEpeak relation between normal and pathological states cannot be made. We predict, and results from analysis of previously acquired data obtained for other purposes in a modest sample size show, that the slope and intercept of the kxo vs. cEpeak relation changes (decreases) in the setting of dysfunction. The full potential of the index resides in performance of follow-up echocardiographic studies, where each person is his own control, showing alteration of E-wave morphology in response to pharmacological therapy. Whether observed E-wave changes are due to load, or due to intrinsic changes in the LV in response to therapy (LV remodeling), could be addressed via slope M and intercept B. Because no prediction regarding the existence of any LIIF or its noninvasive validation in normal subjects has been previously achieved, the proposed approach is a reasonable first step.

In conclusion, it is known that changing tilt table position alters load in a predictable manner and affects E-waves in a predictable way (3, 8, 11, 14, 17, 26, 32, 33, 38, 4042, 44, 45, 48). Changes in load thereby affect the E-wave determined PDF parameters (k, xo, and c). The observed E-wave variation reflects the physiological mechanisms that respond to changes in load. Because the slope M of the maximum driving force vs. peak resistive force relationship remains constant in response to alteration of tilt table-induced load (Fig. 2), we conclude that M represents a LIIF in normal subjects. To determine whether M differentiates between control and DD states, we analyzed previously acquired cath-echo data. Results show M to be significantly lower in DD states (6 subjects with increased LVEDP and normal EF), compared with control states (5 subjects with normal LVEDP and normal EF). We emphasize that additional prospective studies, specifically designed to test load independence in well-defined pathophysiological subsets, are required to fully characterize and establish the utility of our proposed noninvasive LIIF.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This study was supported in part by the Whitaker Foundation, the National Heart, Lung, and Blood Institute (HL-54179, HL-04023), the American Heart Association, the Alan A. and Edith L. Wolff Charitable Trust, and the Barnes-Jewish Hospital Foundation.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We acknowledge the assistance of anonymous reviewers in revising this manuscript.


    FOOTNOTES
 

Address for reprint requests and other correspondence: S. J. Kovács, Cardiovascular Biophysics Laboratory, Washington Univ. Medical Center, 660 South Euclid Ave., Box 8086, St. Louis, MO 63110 (e-mail: sjk{at}wuphys.wustl.edu)

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. Appleton CP, Firstenberg MS, Garcia MJ, and Thomas JD. The echo-Doppler evaluation of left ventricular diastolic function. Cardiol Clin 18: 513–546, 2000.[CrossRef][Medline]
  2. Bauman L, Chung CS, Karamanoglu M, and Kovács S. The peak atrioventricular pressure gradient to transmitral flow relation: kinematic model prediction with in-vivo validation. J Am Soc Echocardiogr 8: 839–844, 2004.
  3. Choong CY, Herrmann HC, Weyman AE, and Fifer MA. Preload dependence of Doppler-derived indexes of left ventricular diastolic function in humans. J Am Coll Cardiol 10: 800–808, 1987.[Abstract]
  4. Chung CS and Kovács SJ. Consequences of increasing heart rate on deceleration time, the velocity time integral, and E/A. Am J Cardiol 97: 130–136, 2006.[CrossRef][Web of Science][Medline]
  5. Courtois M, Ludbrook PA, and Kovacs SJ. Unsolved problems in diastole. Cardiol Clin 18: 653–667, 2000.[CrossRef][Medline]
  6. Dent CL, Bowman AW, Scott MJ, Allen JS, Lisauskas JB, Janif M, Wickline SA, and Kovács SJ. Echocardiographic characterization of fundamental mechanisms of abnormal diastolic filling in diabetic rats with a parametrized diastolic filling formalism. J Am Soc Echocardiogr 14: 1166–1172, 2001.[CrossRef][Web of Science][Medline]
  7. Dirschedl P, Gregull A, and Lollgen H. Volume loading of the heart by "leg up" position and head down tilting (–6 degrees) (HDT). Acta Astronaut 27: 41–43, 1992.[CrossRef][Web of Science][Medline]
  8. Downes TR, Nomeir AM, Stewart K, Mumma M, Kerensky R, and Little WC. Effect of alteration in loading conditions on both normal and abnormal patterns of left ventricular filling in healthy individuals. Am J Cardiol 65: 377–382, 1990.[CrossRef][Web of Science][Medline]
  9. Garcia MJ, Smedira NG, Greenberg NL, Main M, Firstenberg MS, Odabashian J, and Thomas JD. Color M-mode Doppler flow propagation velocity is a preload insensitive index of left ventricular relaxation: animal and human validation. J Am Coll Cardiol 35: 201–208, 2000.[Abstract/Free Full Text]
  10. Granzier HL and Labeit S. The giant protein titin: a major player in myocardial mechanics, signaling, and disease. Circ Res 94: 284–295, 2004.[Abstract/Free Full Text]
  11. Guazzi M, Maltagliati A, Tamborini G, Celeste F, Pepi M, Muratori M, Berti M, and Guazzi MD. How the left and right sides of the heart, as well as pulmonary venous drainage, adapt to an increasing degree of head-up tilting in hypertophic cardiomyopathy: differences from the normal heart. J Am Coll Cardiol 36: 195–193, 2000.
  12. Hall AF and Kovács SJ. Automated method for characterization of diastolic transmitral Doppler velocity contours: early rapid filling. Ultrasound Med Biol 20: 107–116, 1994.[CrossRef][Web of Science][Medline]
  13. Hasegawa H, Little WC, Ohno M, Brucks S, Morimoto A, Cheng HJ, and Cheng CP. Diastolic mitral annular velocity during the development of heart failure. J Am Coll Cardiol 41: 1590–1597, 2003.[Abstract/Free Full Text]
  14. Hurrell DG, Nishimura RA, Ilstrup DM, and Appleton CP. Utility of preload alteration in assessment of left ventricular filling pressure by Doppler echocardiography: a simultaneous catheterization and Doppler echocardiographic study. J Am Coll Cardiol 30: 459–467, 1997.[Abstract]
  15. Jacques DC, Pinsky MR, Severyn D, and Gorcsan J. Influence of alterations in loading on mitral annular velocity by tissue doppler echocardiography and its associated ability to predict filling pressures. Chest 126: 1910–1918, 2004.[Abstract/Free Full Text]
  16. Kass DA, Bronzwaer JGF, and Paulus WJ. What mechanisms underlie diastolic dysfunction in heart failure? Circ Res 94: 1533–1542, 2004.[Abstract/Free Full Text]
  17. Kmetzo JJ, Plotnick GD, and Gottdiener JS. Effect of postural changes and isometric exercise on Doppler-derived measurements of diastolic function in normal subjects. Chest 100: 357–363, 1991.[Abstract/Free Full Text]
  18. Kovács SJ, Meisner JS, and Yellin EL. Modeling of diastole. Cardiol Clin 18: 459–487, 2000.[CrossRef][Medline]
  19. Kovács SJ, McQueen MD, and Peskin CS. Modeling cardiac fluid dynamics and diastolic function. Phil Trans R Soc Lond A Math Phys Eng Sci 359: 1299–1314, 2001.[CrossRef]
  20. Kovács SJ, Barzilai B, and Perez JE. Evaluation of diastolic function with Doppler echocardiography: the PDF formalism. Am J Physiol Heart Circ Physiol 252: H178–H187, 1987.[Abstract/Free Full Text]
  21. Kovács SJ, Setser R, and Hall AF. Left ventricular chamber stiffness from model based image processing of transmitral Doppler E-waves. Coron Artery Dis 8: 179–187, 1997.[Web of Science][Medline]
  22. Kovács SJ, Rosado J, McGuire ALM, and Hall AF. Can transmitral Doppler E-waves differentiate hypertensive hearts from normal? Hypertension 30: 788–795, 1997.[Abstract/Free Full Text]
  23. Lipton P. Testing hypotheses: prediction and prejudice. Science 307: 219–221, 2005.[Abstract/Free Full Text]
  24. Lisauskas JB, Singh J, Courtois MR, and Kovács SJ Jr. The relation of the peak Doppler E-wave to peak mitral annulus velocity ratio to diastolic function. Ultrasound Med Biol 27: 499–507, 2001.[CrossRef][Web of Science][Medline]
  25. Lisauskas JB, Singh J, Bowman AW, and Kovács SJ. Chamber properties from transmitral flow: prediction of average and passive left ventricular diastolic stiffness. J Appl Physiol 91: 154–162, 2001.[Abstract/Free Full Text]
  26. Masuyama T, St Goar FG, Alderman EL, and Popp RL. Effects of nitroprusside on transmitral flow velocity patterns in extreme heart failure: a combined hemodynamic and Doppler echocardiographic study of varying loading conditions. J Am Coll Cardiol 16: 1175–1185, 1990.[Abstract]
  27. Nagueh SF, Lakkis NM, Middleton KJ, Spencer WH, Zoghbi WA, and Quinones MA. Doppler estimation of left ventricular filling pressures in patients with hypertophic cardiomyopathy. Circulation 99: 254–261, 2003.
  28. Oguzhan A, Arinc H, Abaci A, Topsakal R, Eryol NK, Ozdogru I, Basar E, and Ergin A. Preload dependence of Doppler tissue imaging derived indexes of left ventricular diastolic function. Echocardiography 22: 320–325, 2005.[CrossRef][Web of Science][Medline]
  29. Oomen B, Karamanoglu M, and Kovács SJ. Modeling time varying elastance: the meaning of "load-independence". Cardiovasc Eng 3: 123–130, 2003.[CrossRef]
  30. Paelinck BP, van Eck JWM, De Hert SG, and Gillebert TC. Effects of postural changes on cardiac function in healthy subjects. Eur J Echocardiogr 4: 196–201, 2003.[Abstract/Free Full Text]
  31. Pela G, Regolisti G, Coghi P, Cabassi A, Basile A, Cavatorta A, Manca C, and Borghetti A. Effects of the reduction of preload on left and right ventricular myocardial velocities analyzed by Doppler tissue echocardiography in healthy subjects. Eur J Echocardiogr 5: 262–271, 2004.[Medline]
  32. Pepi M, Guazzi M, Maltagliati A, Berna G, and Tamborini G. Diastolic ventricular interaction in normal and dilated heart during head-up tilting. Clin Cardiol 23: 665–672, 2000.[Web of Science][Medline]
  33. Plotnick GD, Kmetzo JJ, and Gottdiener JS. Effect of autonomic blockade, postural changes and isometric exercise on Doppler indexes of diastolic left ventricular function. Am J Cardiol 67: 1284–1290, 1991.[CrossRef][Web of Science][Medline]
  34. Rich MW, Stitziel NO, and Kovács SJ. Prognostic value of diastolic filling parameters derived using novel image processing technique in patients ≥ 70 years of age with congestive heart failure. Am J Cardiol 84: 82–86, 1999.[CrossRef][Web of Science][Medline]
  35. Riordan MM, Chung CS, and Kovács SJ. Diabetes and diastolic function: stiffness and relaxation from transmitral flow. Ultrasound Med Biol 31: 1589–1596, 2005.[CrossRef][Web of Science][Medline]
  36. Sagawa K. Editorial: the end-systolic pressure-volume relation of the ventricle: definition, modifications and clinical use. Circulation 63: 1223–1227, 1981.[Free Full Text]
  37. Soto PF and Kovacs SJ. The diagnosis and treatment of diastolic heart failure. Cardiovasc Rev Rep 23: 493–500, 2002.
  38. Stoddard MF, Pearson AC, Kern MJ, Ratcliff J, Mrosek DG, and Labovitz AJ. Influence of alteration in preload on the pattern of left ventricular diastolic filling as assessed by Doppler echocardiography in humans. Circulation 79: 1226–1236, 1989.[Abstract/Free Full Text]
  39. Suga H, Sagawa K, and Shoukas AA. Load independence of the instantaneous pressure-volume ratio of the canine left ventricle and effects of epinephrine and heart rate of the ratio. Circ Res 32: 314–322, 1973.[Abstract/Free Full Text]
  40. Suzuki T, Sato K, and Aoki K. Influence of postural change on transmitral flow velocity profile assessed by pulsed Doppler echocardiography in normal individuals and in patients with myocardial infarction. Am Heart J 120: 110–115, 1990.[CrossRef][Web of Science][Medline]
  41. Sztajzel J, Ruedin P, Monin C, Stoermann C, Leski M, Rutishauser W, and Lerch R. Effect of altered loading conditions during haemodialysis on left ventricular filling pattern. Eur Heart J 14: 655–661, 1993.[Abstract/Free Full Text]
  42. Takahashi T, Iizuka M, Sato H, Serizawa T, Momomura S, Mochizuki T, Kohmoto O, Aoyagi T, Matsui H, Ikenouchi H, et al. Doppler echocardiographic-determined changes in left ventricular diastolic filling flow velocity during the lower body positive and negative pressure method. Am J Cardiol 65: 237–241, 1990.[CrossRef][Web of Science][Medline]
  43. Tanabe M, Onishi K, Dohi K, Funabiki Koji T, Kitamura T, Ito M, Nobori T, and Nakano T. Change in filling pattern with preload reduction reflects left ventricular relaxation. Int J Cardiol 98: 67–72, 2005.[CrossRef][Web of Science][Medline]
  44. Triulzi MO, Castini D, Ornaghi M, and Vitolo E. Effects of preload reduction on mitral flow velocity pattern in normal subjects. Am J Cardiol 66: 995–1001, 1990.[CrossRef][Web of Science][Medline]
  45. Voller H, Uhrig A, Spielberg C, von Ameln H, Schroder K, Bruggemann T, and Schroder R. Acute alterations of pre- and afterload: are Doppler-derived diastolic filling patterns able to differentiate the loading condition? Int J Card Imaging 9: 231–240, 1993.[CrossRef][Medline]
  46. Weyman AE. Principles and Practice of Echocardiography (2nd ed.). Philadelphia, PA: Lea & Febiger, 1994, p. 575–624.
  47. Yalcin F, Kaftan A, Muderrisoglu H, Korkmaz ME, Flachskampf F, Garcia M, and Thomas JD. Is Doppler tissue velocity during early left ventricular filling preload independent? Heart 87: 336–339, 2002.[Abstract/Free Full Text]
  48. Yamamoto K, Masuyama T, Tanouchi J, Uematsu M, Doi Y, Mano T, Hori M, Tada M, and Kamada T. Peak early diastolic filling velocity may decrease with preload augmentation: effect of concomitant increase in the rate of left atrial pressure drop in early diastole. J Am Soc Echocardiogr 6: 245–254, 1993.[Medline]



This article has been cited by other articles:


Home page
Eur J Heart FailHome page
E.-Y. Choi, B. W. Choi, S.-A. Kim, S. J. Rhee, C. Y. Shim, Y. J. Kim, S.-M. Kang, J.-W. Ha, and N. Chung
Patterns of late gadolinium enhancement are associated with ventricular stiffness in patients with advanced non-ischaemic dilated cardiomyopathy
Eur J Heart Fail, June 1, 2009; 11(6): 573 - 580.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
L. Shmuylovich and S. J. Kovacs
Stiffness and relaxation components of the exponential and logistic time constants may be used to derive a load-independent index of isovolumic pressure decay
Am J Physiol Heart Circ Physiol, December 1, 2008; 295(6): H2551 - H2559.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
C. S. Chung and S. J. Kovacs
Physical determinants of left ventricular isovolumic pressure decline: model prediction with in vivo validation
Am J Physiol Heart Circ Physiol, April 1, 2008; 294(4): H1589 - H1596.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
M. M. Riordan, E. P. Weiss, T. E. Meyer, A. A. Ehsani, S. B. Racette, D. T. Villareal, L. Fontana, J. O. Holloszy, and S. J. Kovacs
The effects of caloric restriction- and exercise-induced weight loss on left ventricular diastolic function
Am J Physiol Heart Circ Physiol, March 1, 2008; 294(3): H1174 - H1182.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
L. Shmuylovich and S. J. Kovacs
E-wave deceleration time may not provide an accurate determination of LV chamber stiffness if LV relaxation/viscoelasticity is unknown
Am J Physiol Heart Circ Physiol, June 1, 2007; 292(6): H2712 - H2720.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
C. S. Chung, A. Strunc, R. Oliver, and S. J. Kovacs
Diastolic ventricular-vascular stiffness and relaxation relation: elucidation of coupling via pressure phase plane-derived indexes
Am J Physiol Heart Circ Physiol, November 1, 2006; 291(5): H2415 - H2423.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
101/1/92    most recent
01305.2005v1
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 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 Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Shmuylovich, L.
Right arrow Articles by Kovács, S. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Shmuylovich, L.
Right arrow Articles by Kovács, S. J.


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