|
|
||||||||
INNOVATIVE METHODOLOGY
1Cardiac MR Research Group, Division of Imaging Sciences, Kings College London; 2Department of Congenital Heart Disease, Guys Hospital; and 3Cardiothoracic Unit, Institute of Child Health and Great Ormond Street Hospital, London, United Kingdom
Submitted 13 September 2004 ; accepted in final form 18 April 2005
| ABSTRACT |
|---|
|
|
|---|
pulmonary hypertension; nitric oxide; magnetic resonance imaging
In pulmonary hypertensive disease, vessel wall distension and remodeling result in decreased arterial compliance (13, 26), causing an increase in pulse pressure, peak systolic pressure, and ventricular load. Morphological changes in the vessel wall, decreased segmental arterial compliance, and changes in compliance in response to vasodilators have been demonstrated in patients with pulmonary hypertensive disease (6, 22). Thus measuring total arterial compliance or elastance in the pulmonary arterial system in the clinical setting may provide additional diagnostic and prognostic information.
A simple approximation of compliance is the ratio of stroke volume to pulse pressure (Csv). However, it is thought to be an overestimate of true arterial compliance (23). The pulse pressure method (Cppm), based on the parameter optimization of a two-element windkessel model, is considered to be a more accurate estimate of true compliance (23, 27, 29). Both methods have been successfully used to calculate pulmonary arterial compliance in animals, and a close correlation between Csv and Cppm has been demonstrated (23). Thus Csv can be used to validate Cppm in vivo. Both methods require simultaneous measurement of arterial flow (stroke volume for Csv and flow curves for Cppm) and pulse pressure. In addition, the pulse pressure method requires calculation of vascular resistance. Such physiological information is difficult to acquire in humans in the clinical setting.
Velocity-encoded phase-contrast magnetic resonance (MR) enables quantification of volume flow in the pulmonary artery. Cardiac output and the pulmonary-to-systemic flow ratio (
p/
s) measured using this technique have been shown to be accurate (3, 4, 21). In addition, phase-contrast MR has been validated in numerous phantom experiments (3, 21). A new technique, MR-guided cardiac catheterization, allows simultaneous acquisition of phase-contrast MR information and invasive pressure measurements (18, 20). We previously demonstrated the feasibility of using this technique to calculate pulmonary vascular resistance in humans (18). We build on this previous work to evaluate the feasibility of calculating total arterial compliance in humans undergoing MR-guided cardiac catheterization. In this study, compliance was calculated using the pulse pressure method and the ratio of stroke volume to pulse pressure. In addition, measurement of other pulmonary hemodynamic parameters allowed a preliminary study of the relation between compliance and resistance and between compliance and pressure, particularly in response to NO.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Cardiac catheterization. Subjects underwent cardiac catheterization in an MR interventional suite (maximum gradient performance = 30 mT/m amplitude, slew rate = 150 T·m1·s1; 1.5-T Intera I/T MRI scanner, Philips Medical Systems) with X-ray backup (BV Pulsera cardiac X-ray unit, Philips Medical Systems) in the same room (XMR suite). General anesthesia was induced and maintained with propofol. The catheter was advanced under MR guidance or with a combination of X-ray and MR guidance, depending on ease of catheter manipulation. Angiographic catheters (4- to 6-Fr) with CO2-filled balloons were used for all procedures, allowing visualization in both imaging modalities. Hemodynamic variables were collected at 30% O2 (condition 1) and at 20 ppm NO + 30% O2 (condition 2).
MR catheter visualization. All imaging was performed using an adult or pediatric two-element receive-only phased-array radio-frequency coil (Flex L or Flex M, Philips Medical Systems); transmission was performed using the body coil. An interactive steady-state free precession sequence (TR = 2.9 ms, TE = 1.45 ms, matrix = 128 x 128, field of view = 250350 mm, 1014 frames/s, flip angle = 45°, slice thickness = 68 mm, free breathing acquisition) was used to visualize the catheter balloon during manipulation within the heart and great vessels. The imaging plane could be manipulated in real time into any plane to follow the catheter (17).
Flow and pressure quantification.
Pulmonary arterial flow data were acquired using a flow-sensitive, fast-field echo sequence (TR
9 ms, TE
5 ms, matrix = 128192 x 256, field of view = 250350 mm, flip angle = 15°, slice thickness = 57 mm, free breathing acquisition, 3 signal averages, retrospective gating, 40 phases, in-plane and through-plane resolution optimized to patient size). A dedicated nonlinear phase-correction filter, based on Chebichev polynomials (Philips Medical Systems), was used to correct for phase errors introduced by eddy currents and Maxwell terms. In addition, the image plane was centered in the bore of the magnet to further reduce the phase errors. Image planes were located at the midpoint of the pulmonary artery. Systolic pressure (Psystole), diastolic pressure (Pdiastole), and mean pressure (Pmean) were collected during MR flow data acquisition by means of standard fluid-filled catheters coupled to a computerized data-recording system, with a quoted measurement error of ±2 mmHg (S5 collect system, Datex Ohmeda). Pulse pressure (Ppulse) was calculated as Psystole Pdiastole.
Data analysis. Pulmonary blood flow curves and stroke volumes were generated from phase-contrast MR data using a semiautomatic vessel edge-detection algorithm (Flow, Medis) with operator correction. Pulmonary artery cardiac output is the product of pulmonary blood stroke volume and heart rate. Resistance (R) was calculated by dividing mean pulmonary arterial pressure by the pulmonary artery cardiac output and indexing to body surface area (23).
Calculation of total arterial compliance and elastance.
The equation defining the two-element windkessel model is as follows (16)
![]() |
(t)] with 50 cardiac cycles was constructed by repeating the single cardiac cycle phase-contrast MR flow curve.
(t) was used as an input to the windkessel model. R was a known variable for each patient at each condition and was entered into the model. For a given C, the equation was integrated using P(0) = 0 as an initial condition, generating a 50-cardiac cycle pressure curve. The modeled Pdiastole and Psystole are the minimum and maximum pressures between the 40th and 50th cycles, at which point the pressure curve has stabilized and is independent of initial conditions. The modeled Ppulse = modeled Psystole modeled Pdiastole. A series of pressure curves were generated using values of C between 0.001 and 7 ml/mmHg. The modulus of the difference between the modeled Ppulse and the invasively measured Ppulse was calculated, and an error plot was produced. The compliance that coincides with the global minimum of the error plot is taken as the best estimate of C, which was indexed to body surface area (which we call Cppm). All windkessel modeling was performed using Matlab (Mathworks). Elastance (E) is the inverse of C. Statistical analysis. Values are means (SD) unless otherwise stated. A two-tailed t-test was used to compare the hemodynamic responses to NO. A two-tailed t-test, correlation coefficients, linear regression analysis, and Bland-Altman analysis (7) were used to compare Cppm with Csv. Bias was the mean of the difference between the two methods, and agreement was the mean (2SD). P < 0.05 was taken as statistically significant. Correlation coefficients and linear regression analysis were also used to study the relation between steady-state hemodynamic variables and compliance. Statistical analysis was performed using Matlab.
| RESULTS |
|---|
|
|
|---|
|
, and R) in conditions 1 and 2 are shown in Table 1. At baseline, pulmonary arterial pressures and resistance were raised compared with the normal quoted range. There were statistically significant declines in R, Psystole, and Ppulse in response to 20 ppm NO.
|
|
|
10% was considered a response to NO (30). A response in Cppm was demonstrated in 7 of 17 patients and in Csv in 9 of 17 patients. The mean increase in Cppm was 16.9% (SD29.6), and the mean increase in Csv was 17.0% (SD30.1). There was no significant difference between percent increase in Cppm and Csv (P = 0.48). There was excellent correlation between the percent increase in Cppm and Csv (r = 0.89, P < 0.001, m = 0.87, c = 1.95). As a population, although mean Cppm and Csv increased in response to 20 ppm NO, statistical significance was not reached (P = 0.12 for both). Relation between steady-state hemodynamic variables and compliance. As pulmonary hypertensive disease worsens, compliance is thought to decrease. Thus R and Pmean were correlated with 1/Cppm = E, rather than Cppm (Figs. 4 and 5). There was good correlation between R and E in condition 1 (r = 0.86, P < 0.001, m = 0.54,c = 0.22), condition 2 (r = 0.92, P < 0.001, m = 0.56, c = 0.19), and overall (r = 0.89, P < 0.001, m = 0.55, c = 0.20). There was moderate correlation between Pmean and E in condition 1 (r = 0.72, P < 0.001, m = 12.6, c = 15.5), condition 2 (r = 0.73, P < 0.001, m = 14.0, c = 15.5), and overall (r = 0.72, P < 0.001, m = 13.1, c = 15.7). There was slight correlation between the percent change in E and R (r = 0.47, P = 0.06) and poor correlation between the percent change in E and Pmean (r = 0.02, P = 0.92).
|
|
| DISCUSSION |
|---|
|
|
|---|
Estimates of total arterial compliance. Several methods have been proposed to calculate total arterial compliance. The ratio of stroke volume to pulse pressure (Csv) is the simplest method and has been used to estimate pulmonary arterial compliance in humans. However, Csv is an overestimate of arterial compliance, inasmuch as it fails to take into consideration outflow from the arterial system during the ejection phase. Nevertheless, Csv has been shown to correlate strongly with other measures of compliance, and it also represents the upper limit of arterial compliance. It is therefore a useful method of validating other estimates of arterial compliance. Most other estimates of arterial compliance are based on the arterial windkessel model (16, 29). The diastolic decay method (19) and the area method (16) are widely used to estimate systemic arterial compliance based on the two-element windkessel model. Both require detailed pressure curves and are only applicable when there is zero flow in diastole (27, 29). These methods may not be suited to patients with pulmonary hypertension and pulmonary regurgitation. A method based on parameter optimization of a three-element windkessel model has also been used to calculate total pulmonary arterial compliance (23, 24), but this method overestimates pulmonary arterial compliance (23). The pulse pressure method is based on parameter estimation of a two-element windkessel model and has been found to be an accurate method of calculating the total arterial compliance in the pulmonary (23) and systemic vasculatures (2729). It does not require measurement of a detailed pressure curve or zero flow conditions in diastole and is well suited to patients with pulmonary hypertension (29). This method requires measurement of volume flow as a function of time (flow curves), pulse pressure, and vascular resistance (R). In the pulmonary artery, pressure can be measured at catheterization with invasive manometry. In animal studies, perivascular flow probes can be used to accurately measure flow curves (23); unfortunately, use of such methods is difficult in humans. In humans, blood flow velocity curves can be measured invasively (15) or noninvasively (14, 30). Combining this data with a measurement of cardiac output (using an indicator-dilution method) allows generation of flow curves (14, 15, 30). However, inaccuracies in measurement of cardiac output by indicator-dilution methods (10, 11) and measurement of blood flow velocities can lead to errors in measurement of volume flow.
MR-guided cardiac catheterization. Velocity-encoded phase-contrast MR is a validated technique that enables measurement of flow curves in major vessels (3, 4, 18). It is therefore well suited to generate flow curves for estimation of pulmonary arterial compliance by the pulse pressure method. MR-guided cardiac catheterization is a new technique that allows greater access to physiological and anatomic information (20). We previously demonstrated the feasibility of simultaneously acquiring phase-contrast MR and invasive pressure data in patients undergoing MR-guided cardiac catheterization to calculate pulmonary vascular resistance (18). In this study, we have used these data to calculate total arterial compliance by the pulse pressure method and as the ratio of stroke volume to pulse pressure.
Total arterial compliance. We have demonstrated the feasibility of calculating total pulmonary arterial compliance in patients undergoing MR-guided cardiac catheterization. As reported in other studies, Csv was consistently higher than Cppm (23). It is believed that Csv overestimates total arterial compliance, inasmuch as it assumes that the total stroke volume is buffered in the arteries (23). In fact, the volume buffered in the arteries is the stroke volume minus the arterial outflow volume during ejection. Our results demonstrated excellent correlation between Csv and Cppm in conditions 1 and 2 and overall, as well as good correlation between the percent change in Csv and Cppm, in agreement with previous studies (23). These results suggest that Cppm is a precise estimate of true arterial compliance, although we are unable to verify its accuracy, inasmuch as there is no in vivo "gold standard" method of measuring arterial compliance. However, previous work using models of the arterial system has shown Cppm to be an accurate estimate of total arterial compliance (27).
Relation to hemodynamic variables and response to NO. Our results show that, within our study population, R and Pmean are inversely related to Cppm, with good to moderate correlation. There are two possible explanations for these findings: 1) Vascular remodeling, which leads to an increase in R and, thus, in Pmean, causes vessel distension and reduced compliance. Studies have shown that reduced compliance in systemic arteries is secondary to distension. 2) Vascular remodeling occurs in the large conduit arteries as well as the resistance arteries as part of the pulmonary hypertensive disease process. Vessel wall thickening in large pulmonary arteries has been demonstrated using intravascular ultrasound and is associated with decreased vessel compliance (6, 22). It is probable that both factors contribute to the relation between R, Pmean, and Cppm. However, the strong correlation between Pmean and Cppm do suggest that distension of the conduit arteries is an important factor. To fully understand the relation between R, Pmean, and Cppm, this technique will need to be performed on a larger group of patients with pulmonary hypertensive disease.
Inhaled NO acts on pulmonary vascular smooth muscle, causing smooth muscle relaxation (12). NO causes a fall in resistance due to vasodilation of resistance arteries. In this study, we have demonstrated a hemodynamically significant increase in Cppm in response to NO in seven patients. However, the population change was small and not statistically significant. There was mild correlation between the percent change in R and 1/Cppm. A possible explanation for these findings is that a decrease in R leads to a decrease in Pmean, causing less distension and an increase in Cppm. In this study group, there was little correlation between the percent change in Pmean and 1/Cppm. Therefore, increase in compliance in response to NO cannot be explained purely in terms of reduced distension. Another explanation is that NO decreases resistance and increases compliance through its actions on smooth muscle distributed throughout the pulmonary vasculature. However, it should be noted that inhaled NO may not have significant bioavailability in the large conduit arteries.
Limitations. Phase-contrast MR flow is less accurate in patients with arrhythmias during acquisition and/or turbulent blood flow; the presence of these conditions is a general limitation of this technique, however, we did not encounter these events in our study. We did not repeat previous work internally validating phase-contrast MR in our patient population, inasmuch as we did not wish to prolong the procedure times. The study population was small and inhomogeneous. In addition, the change in hemodynamic parameters in response to NO was limited, possibly because of general anesthesia. This population was sufficient to prove the feasibility of this technique. However, study of the relation between the various hemodynamic parameters was hampered.
Conclusion. Using simultaneously acquired invasive pressure measurements and MR flow data, we have demonstrated the feasibility of quantifying total arterial compliance. This additional information can be combined with information regarding resistance to tailor patient management, particularly with the advent of new medical therapies. With proof of the feasibility of this technique, it will be possible to investigate the relation between compliance and other hemodynamic parameters in a larger more homogeneous population. In our facility, cardiac catheterization was accomplished by MR guidance with X-ray backup. Pulmonary artery catheterization can be accomplished without imaging before MR scanning, thus potentially allowing this technique to be used in a conventional cardiac MR scanner. Further investigation using this technique is required to fully elucidate the relation between compliance, pulmonary hypertensive disease, and medical intervention. In conclusion, we believe that measurement of pulmonary arterial compliance will give new insights into the pathophysiology of pulmonary hypertension and may be important in assessing suitability for long-term medical therapy.
| GRANTS |
|---|
|
|
|---|
| DISCLOSURES |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
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 |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J.-W. Lankhaar, N. Westerhof, T. J.C. Faes, C. Tji-Joong Gan, K. M. Marques, A. Boonstra, F. G. van den Berg, P. E. Postmus, and A. Vonk-Noordegraaf Pulmonary vascular resistance and compliance stay inversely related during treatment of pulmonary hypertension Eur. Heart J., July 1, 2008; 29(13): 1688 - 1695. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Morimont, B. Lambermont, A. Ghuysen, P. Gerard, P. Kolh, P. Lancellotti, V. Tchana-Sato, T. Desaive, and V. D'Orio Effective arterial elastance as an index of pulmonary vascular load Am J Physiol Heart Circ Physiol, June 1, 2008; 294(6): H2736 - H2742. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Vonk-Noordegraaf, J.-W. Lankhaar, M. J.W. Gotte, J. T. Marcus, P. E. Postmus, and N. Westerhof Magnetic resonance and nuclear imaging of the right ventricle in pulmonary arterial hypertension Eur. Heart J. Suppl., December 1, 2007; 9(suppl_H): H29 - H34. [Abstract] [Full Text] [PDF] |
||||
![]() |
G.-P. Diller, K. Dimopoulos, H. Kafka, S. Y. Ho, and M. A. Gatzoulis Model of chronic adaptation: right ventricular function in Eisenmenger syndrome Eur. Heart J. Suppl., December 1, 2007; 9(suppl_H): H54 - H60. [Abstract] [Full Text] [PDF] |
||||
![]() |
J.-W. Lankhaar, N. Westerhof, T. J. C. Faes, K. M. J. Marques, J. T. Marcus, P. E. Postmus, and A. Vonk-Noordegraaf Quantification of right ventricular afterload in patients with and without pulmonary hypertension Am J Physiol Heart Circ Physiol, October 1, 2006; 291(4): H1731 - H1737. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |