Am J Physiol Heart Circ Physiol 288: H1874-H1884, 2005.
First published December 16, 2004; doi:10.1152/ajpheart.00754.2004
0363-6135/05 $8.00
A ventricular-vascular coupling model in presence of aortic stenosis
Damien Garcia,1
Paul J. C. Barenbrug,2
Philippe Pibarot,3
André L. A. J. Dekker,2
Frederik H. van der Veen,2
Jos G. Maessen,2
Jean G. Dumesnil,3 and
Louis-Gilles Durand1
1Clinical Research Institute of Montreal, Montreal, Quebec, Canada; 2Department of Cardiothoracic Surgery, Academic Hospital, Maastricht, The Netherlands; and 3Quebec Heart Institute, Laval Hospital, Laval University, Quebec, Canada
Submitted 27 July 2004
; accepted in final form 8 December 2004
 |
ABSTRACT
|
|---|
In patients with aortic stenosis, the left ventricular afterload is determined by the degree of valvular obstruction and the systemic arterial system. We developed an explicit mathematical model formulated with a limited number of independent parameters that describes the interaction among the left ventricle, an aortic stenosis, and the arterial system. This ventricular-valvular-vascular (V3) model consists of the combination of the time-varying elastance model for the left ventricle, the instantaneous transvalvular pressure-flow relationship for the aortic valve, and the three-element windkessel representation of the vascular system. The objective of this study was to validate the V3 model by using pressure-volume loop data obtained in six patients with severe aortic stenosis before and after aortic valve replacement. There was very good agreement between the estimated and the measured left ventricular and aortic pressure waveforms. The total relative error between estimated and measured pressures was on average (standard deviation) 7.5% (SD 2.3) and the equation of the corresponding regression line was y = 0.99x 2.36 with a coefficient of determination r2 = 0.98. There was also very good agreement between estimated and measured stroke volumes (y = 1.03x + 2.2, r2 = 0.96, SEE = 2.8 ml). Hence, this mathematical V3 model can be used to describe the hemodynamic interaction among the left ventricle, the aortic valve, and the systemic arterial system.
mathematical modeling; cardiovascular system; cardiac catheterization; left ventricle
LEFT VENTRICULAR (LV) pressure, aortic pressure, and cardiac output result from a matching among the oxygen demand of the body, the LV performance, and the LV afterload. The presence of an aortic valve stenosis causes an obstruction to LV outflow, thus resulting in an increase in LV afterload (20). Patients with aortic stenosis also often have concomitant diseases, including hypertension, hyperlipidemia, diabetes, and atherosclerosis (1, 6, 22, 25). These diseases have been shown to alter the structural and functional properties of the systemic arterial system (8, 40). More specifically, they may reduce arterial elasticity and/or increase arteriolar resistance. Hence in patients with aortic stenosis, the left ventricle is often facing a double load: a valvular load imposed by the aortic stenosis and an arterial load caused by a decrease in systemic arterial compliance (C) and/or an increase in systemic vascular resistance (R). Briand et al. (3) have recently shown that reduced C is a frequent occurrence in patients with aortic stenosis where it independently contributes to increase afterload and decrease LV function. In addition, Antonini-Canterin et al. (1) have reported that symptoms of aortic stenosis develop at a lesser degree of valvular obstruction in hypertensive compared with normotensive patients. It is thus important to assess the respective contributions of the aortic valve and the systemic arterial system to the LV workload in such patients. This information could be useful to preferentially target the treatment on the valve (i.e., valve replacement), on the systemic arterial system (i.e., medical therapy), or on both, given the type and magnitude of afterload augmentation. We thus developed a mathematical ventricular-valvular-vascular model (V3 model) to study the coupled system formed by the left ventricle, a normal or stenotic aortic valve, and the systemic arterial system. A simple mathematical model describing the interaction between the left ventricle and the arterial system has been recently proposed by Segers et al. (2628) and has been proven to correctly predict some clinical observations. We (13) completed the model of Segers with a newly validated equation describing the instantaneous transvalvular pressure-flow relationship. The objective of the present study was to assess the descriptive validity of the resulting V3 model using catheterization data obtained in six patients immediately before and after aortic valve replacement (9).
 |
MATERIALS AND METHODS
|
|---|
Analytic model.
Transvalvular flow rate was simulated using the V3 model (Fig. 1). The latter consists on the interaction among 1) the time-varying elastance model for the left ventricle (35), 2) the instantaneous net pressure gradient across the aortic valve (13), and 3) the three-element windkessel model for the systemic arterial load (38). The LV time-varying elastance (E) model relates the ventricular pressure (PV) to the ventricular volume (V) and the unloaded volume, defined as the intercept of the end-systolic pressure-volume relationship with the volume axis (V0) (35) as follows:
 | (1) |

View larger version (21K):
[in this window]
[in a new window]
|
Fig. 1. Schema of the ventricular-valvular-vascular interaction mathematical model. PV, left ventricular (LV) pressure; V, LV volume; EOA, effective orifice area of the aortic valve; AA, cross-sectional area of the aorta at the sinotubular junction; PA, pressure in the ascending aorta; Z0, aortic characteristic impedance; C, total arterial compliance; R, systemic vascular resistance; PVE, central venous pressure.
|
|
With regard to the aortic valve, we (13) have recently demonstrated that the instantaneous net pressure gradient (TPGnet) across the valve during LV ejection, that is, the difference between the LV pressure and the pressure in the ascending aorta (PA), can be written as:
 | (2) |
where
is the fluid density, Q is the transvalvular flow rate, and ELCo is the energy loss coefficient of the aortic valve (12, 13). ELCo is a function of the valve effective orifice area (EOA) and the aortic cross-sectional area (AA) and is defined as ELCo = EOA AA/(AA EOA). The so-called aortic valve EOA is the minimal cross-sectional area of the flow jet, i.e., the cross-sectional area of the vena contracta, downstream of a native or bioprosthetic aortic heart valve (Fig. 1). Finally, the three-element windkessel model relates the PA to the transvalvular flow rate (Q) as follows:
 | (3) |
where Z0 is the aortic characteristic impedance, R is the total vascular resistance, C is the total arterial compliance, and PVE is the central venous pressure. Figure 1 depicts a schema of the complete coupling model. The sum
(Eq. 2)/
t + (Eq. 2)/(RC) yields:
 | (4) |
The second term of Eq. 3 is substituted for the expression
PA/
t + PA/RC in the first term of Eq. 4, and PV is replaced using Eq. 1. Also, during ejection, transvalvular flow rate can be written as: Q(t) =
V(t)/
t. Therefore, Eq. 4 becomes:
 | (5) |
where
At the onset of the ejection (t = t0), V is equal to the LV end-diastolic volume (LVEDV) and Q = 0. Ejection begins when LV pressure reaches aortic pressure, i.e., when TPGnet is zero. According to Eq. 2,
Q/
t is thus also equal to zero at the ejection onset. Because Q(t) =
V(t)/
t, the initial conditions are therefore:
 | (6) |
The third-order nonlinear differential Eq. 5 and the corresponding initial conditions (Eq. 6) completely describe the LV volume during ejection in the presence of aortic stenosis under the conditions that the ventricular contractility, the arterial properties, and the aortic stenosis severity are known. The transvalvular flow rate is simply calculated using the negative derivative of V(t). It has been shown that the normalized elastance (EN) is quasi-independent of the preload, the afterload, and the cardiac inputs (30, 35). Elastance E(t) can thus be defined from its peak value (Emax) and time to Emax (TEmax) as follows: E(t) = Emax EN(t/TEmax) (30). Table 1 summarizes the independent parameters necessary for simulating LV volume during ejection using the V3 model.
Numerical computation.
Time reference (t = 0) is fixed at the onset of the isovolumic contraction. The computational algorithm is somewhat equivalent to the one used by Segers et al. (26). An arbitrary diastolic pressure (DPA) is chosen. Time corresponding to the ejection onset (t0) is determined from Eq. 1 such that t0 satisfies the following equation: E(t0) = DPA/(LVEDV V0). LV volume during ejection is then calculated from Eqs. 5 and 6 using an explicit Runge-Kutta method. The end of LV ejection is reached when Q(t) becomes zero and it is further assumed that Q = 0 during diastole. PA is then deduced using Eq. 3 with the calculated transvalvular flow rate (see APPENDIX). A second DPA is therefore obtained. If the relative difference between the two DPA is >1%, a new iteration is performed using the second DPA value and so on until the desired relative precision is reached. PV during isovolumic contraction, ejection, and isovolumic relaxation is then calculated from Eq. 1 and is extrapolated during the ventricular filling period using an exponential function (see APPENDIX).
Clinical protocol.
Clinical materials and methods are described in detail by Dekker et al. (9). Briefly, LV and aortic pressures, LV volume, and ECG were measured in six symptomatic patients with severe aortic valve stenosis just before and just after the aortic valve replacement. All patients provided their consent, and the study was approved by the institutional committee of the Academic Hospital of Maastricht. LV volume and pressure were obtained by inserting a conductance catheter with a pressure sensor (Sentron ANP-223N, Roden, The Netherlands) into the left ventricle. The conductance catheter was connected to a Leycom Sigma 5DF system (CD Leycom, Zoetermeer, The Netherlands) used in the dual-frequency mode. PA was measured using a pressure tip catheter (Sentron ANP-530N, Roden, The Netherlands) positioned in the aortic root.
Hemodynamic data (LV and aortic pressures and LV volume) were acquired over a few cardiac cycles before cardiopulmonary bypass when the hemodynamics conditions were stable (dataset 1), and then a series of pressure-volume loops were acquired during inferior caval vein occlusion (dataset 2). The same protocol was repeated 15 min after successful weaning from the cardiopulmonary bypass. Pressure-volume loops during caval vein occlusion could not be obtained preoperatively in two patients due to the presence of ectopic beats. The onset of the isovolumic contraction was detected from the peak of the QRS wave of the ECG signal.
Estimation of cardiovascular parameters in patients.
For each patient, the normalized elastance was determined before and after valve replacement from the series of pressure-volume loops contained in dataset 2. Emax and V0 were first calculated according to the method described by Kono et al. (16). With V0 known, elastance was then determined from Eq. 1 for each loop, and the average value was computed and used thereafter. EN was obtained by dividing elastance by Emax. TEmax was calculated as the elapsed time between the onset of the isovolumic contraction and time when elastance is maximal; i.e., when E(t) = Emax. LVEDV was directly obtained from the LV volume waveforms measured during the steady-state condition (dataset 1). The other parameters Emax, V0, R, C, Z0, PVE, and ELCo were estimated using the Nelder-Mead minimization method (17). The function to be minimized was the sum of the relative errors between measured and simulated LV and aortic pressures, respectively. Typical values encountered in the literature were used as initial estimates. We chose to estimate Emax and V0 with the use of the minimization method rather than to use those calculated from the LV pressure-volume loops of dataset 2, because LV contractility may vary significantly during operation (4, 9). Indeed, according to our preliminary calculations, we observed that Emax and V0 measured during caval vein occlusion were different from those calculated during the steady-state condition. Once the minimization was completed, PV and PA and stroke volume (SV) estimated by the model were compared with the averaged experimental data to validate the V3 model.
Comparison of estimated and measured pressures.
Estimated PV and PA were compared with measured pressures by calculating the total relative error defined as:
where |...| represents the vectorial 1-norm. The goodness of fit was evaluated by performing a linear regression using all the data points from the pressure waveforms, and the coefficient of determination (r2), as well as the root mean square error, were determined. We also calculated the mean of the residuals, defined as the differences between the measured pressures and the fitted pressures, and we studied their correlation with the measured pressure data to test whether the residuals displayed a systematic pattern. A similar analysis was performed with LV volumes measured and simulated during LV ejection.
Statistical analysis.
Standard cardiovascular parameters (heart rate; stroke volume; cardiac output; LV end-diastolic and end-systolic volumes; LV ejection fraction; LV end-diastolic, peak, and end-systolic pressures; aortic systolic, diastolic and mean pressures; mean transvalvular pressure gradient; and Gorlin EOA) were measured for each patient. The data measured before valve replacement were compared with those measured after valve replacement with the use of a two-tailed paired t-test. A similar statistical analysis was also performed with the seven cardiovascular parameters estimated by the minimization method (Emax, V0, R, C, Z0, PVE, and ELCo).
 |
RESULTS
|
|---|
Normalized elastance before and after aortic valve replacement.
Aortic valve replacement did not significantly alter EN for a given patient, as shown in Fig. 2. Because it was impossible to measure EN in patients 3 and 5 before valve replacement and given that the EN remained unchanged during surgery in all other patients, we chose to use the EN determined postoperatively to simulate preoperative LV and aortic pressures in these two patients.

View larger version (18K):
[in this window]
[in a new window]
|
Fig. 2. Normalized elastance before (solid lines) and after (dotted lines) aortic valve replacement (AVR). Comparison is not possible in patients 3 and 5 because P-V loop data could not be obtained before surgery.
|
|
Comparison of cardiovascular parameters before and after aortic valve replacement.
The mean values and standard deviation of cardiovascular parameters measured before and after valve replacement are shown in Table 2. Among these parameters, only LV peak systolic pressure, Gorlin EOA, and mean transvalvular pressure gradient were significantly different before versus after valve replacement. Consistently, among the cardiovascular parameters determined from the minimization method, only ELCo changed significantly after replacement (Table 3). Hence, only valvular characteristics were significantly and markedly improved after valve replacement. There was no significant improvement in the vascular and ventricular parameters (except LV peak systolic pressure) immediately after valve replacement. Nonetheless, the parameters of LV systolic function (LV end-systolic volume, LV ejection fraction, and SV) tended to improve after valve replacement, but the differences did not reach statistical significance, likely because of the small number of patients included in this study. Emax decreased in all patients after valve replacement (Table 3). The difference between pre- and postreplacement Emax values was not statistically significant (P = 0.07) when analyzed with a two-sided paired t-test but reached statistical significance (P = 0.03) when analyzed with a one-sided paired t-test. The decrease in Emax after valve replacement may be related to the acute depression of myocardial contractility that is commonly observed after cardiopulmonary bypass. In addition, it is also possible that the reduction in LV systolic pressure resulting from valve replacement may have contributed to the decrease in Emax, because it has been suggested that Emax is not completely independent of LV afterload (34).
View this table:
[in this window]
[in a new window]
|
Table 3. Cardiovascular parameters estimated by the minimization method before and after valve replacement in each patient
|
|
It is also interesting to note that patient 6 had a relatively high residual transvalvular pressure gradient (28 mmHg) after valve replacement (see Fig. 7). The ELCo was only 0.75 cm2 (see Fig. 10), and the energy loss index (ELCo indexed for body surface area: 1.98 m2 in this patient) was 0.38 cm2/m2 after replacement, which is consistent with severe patient-prosthesis mismatch (24) and definitely explains the high gradient in this patient.

View larger version (28K):
[in this window]
[in a new window]
|
Fig. 7. Comparison between measured (solid lines) and estimated (dotted lines) LV and aortic pressure waveforms in the six patients after AVR.
|
|

View larger version (16K):
[in this window]
[in a new window]
|
Fig. 10. Comparison between EOA (in cm2) measured using the Gorlin formula and energy loss coefficient (ELCo) estimated from the minimization method. Dotted line, identity line. Solid and shaded circles correspond to data before and after surgery, respectively. Note the severe patient-prosthesis mismatch in patient 6 (isolated gray circle).
|
|
Estimated versus measured LV volumes and SV.
Overall, there was a very good agreement between estimated and measured LV volume waveforms as shown in Figs. 3 and 4. Total relative error ranged between 1.2% (patient 4 before surgery, Fig. 3) and 8.9% (patient 6 after surgery, Fig. 4) and was on average (standard deviation) 5.3% (SD 2.7). Slope and intercept (with 95% confidence interval) of the regression line between estimated and measured volume data during ejection were 0.95 (0.947, 0.955) and 3.84 (3.49, 4.18) ml, respectively, and r2 was 0.978. The root-mean-square error was 4.4 ml. The mean (SD) of the residuals was 0.00 (SD 4.4) ml, and no correlation was found (r < 103) between the residuals and the volume data, suggesting that the residuals behaved randomly. These data confirm that there was a good fit between the LV volume waveforms simulated by the V3 model and those directly measured in patients. Measured LV SV ranged between 43 and 88 ml. There was an excellent agreement between estimated and measured SV (y = 1.03x + 2.2, r2 = 0.96, SEE = 2.8 ml, Fig. 5). The mean (SD) absolute error was 3.8 (SD 2.7) ml.

View larger version (21K):
[in this window]
[in a new window]
|
Fig. 3. Comparison between measured (solid lines) and estimated (dotted lines) LV volume waveforms during ejection in the six patients before AVR.
|
|

View larger version (20K):
[in this window]
[in a new window]
|
Fig. 4. Comparison between measured (solid lines) and estimated (dotted lines) LV volume waveforms during ejection in the six patients after AVR.
|
|

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 5. Relationship between measured stroke volume (in ml) and that estimated from the minimization method. Dotted line, identity line. Black and grey dots correspond to data before and after surgery, respectively.
|
|
Estimated versus measured pressure waveforms.
Overall, there was a very good agreement between estimated and measured pressure waveforms as shown in Figs. 6 and 7. Total relative error ranged between 3.7% (patient 3 after surgery, Fig. 7) and 10.5% (patient 4 after surgery, Fig. 7) and was on average 7.5% (SD 2.3). Slope and intercept of the regression line between the simulated and the measured pressure data were 0.99 (95% confidence interval: 0.9981, 1.002) and 2.36 (2.50, 2.22) mmHg, respectively, and r2 was 0.98. The root-mean-square error was 6.0 mmHg. The mean (SD) of the residuals was 0.0 (SD 6.0) mmHg, and no correlation was found (r <0.001; P = not significant) between the residuals and the pressure data. These results show that there was an excellent fit between the pressure waveforms generated by the V3 model and the pressure waveforms directly measured in patients.

View larger version (27K):
[in this window]
[in a new window]
|
Fig. 6. Comparison between measured (solid lines) and estimated (dotted lines) LV and aortic pressure waveforms in the six patients before AVR.
|
|
Comparison between estimated and measured cardiovascular parameters.
In this study, seven among the nine parameters included in the model (Table 1) were estimated by the Nelder-Mead minimization method. The individual and average data of these parameters estimated by minimization are shown in Table 3. The minimization method usually provides a local minimum and not a global minimum of the function to be minimized around the chosen starting point. For a complex system, such as the one described by the V3 model, it is uncertain whether there is only one unique set of values for these parameters, which leads to a good concordance between the pressure waveforms. To confirm that the detected local minimum corresponds to the actual value of the cardiovascular parameters, we therefore compared three of these parameters, i.e., R, C, and ELCo, obtained from the minimization method with those directly obtained from the experimental measurements. R and C were calculated with the equations used in clinical practice and were thus represented by the abbreviations Rclin and Cclin, respectively:
 | (7) |
and
 | (8) |
where MAP, PP, and T represent the mean aortic pressure, the pulse pressure, and the cardiac period, respectively. We previously demonstrated that ELCo is equivalent to the valve EOA calculated using the standard Gorlin formula (in cm2): EOAGorlin = Qmean/(44.3
), where Qmean (in ml/s) is the mean transvalvular flow rate and TPGnet (in mmHg) is the mean net transvalvular pressure gradient, which includes the pressure recovery occurring downstream of the vena contracta (39). For the purpose of this analysis, ELCo was thus compared with EOAGorlin.
As shown in Figs. 810, there was a strong correlation between the parameters estimated by the optimization method and the parameters obtained from experimental measurements (estimated vs. measured: resistance, y = 0.91x 0.04, r2 = 0.98; compliance, y = 1.12x 0.15, r2 = 0.94; energy loss coefficient, y = 1.00x + 0.04, r2 = 0.95). This suggests that the minimization algorithm likely converges toward a consistent minimum.

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 8. Relationship between resistance (in mmHg·s·ml1) measured using the standard clinical formula (Eq. 7) and that estimated from the minimization method. Dotted line, identity line. Solid and shaded circles correspond to data before and after surgery, respectively.
|
|
 |
DISCUSSION
|
|---|
V3 model.
A few LV-arterial coupling models have been proposed to analyze the effect of vascular properties on LV function (10, 27, 32, 36, 41). The model proposed by Segers et al. (26, 32) is relatively simple and easily applicable in the clinical setting. This model accurately predicted the SV and the systolic and diastolic aortic pressures in animals and in humans. The model of Segers et al. (26, 32) is based on the combination of the time-varying elastance representation of the left ventricle with a four-element windkessel model. However, this model assumes that the aortic valve is normal and it is thus not applicable in patients with an aortic stenosis. Computer models also including the effect of aortic stenosis have been proposed by Li et al. (18) and Smith et al. (31). In these models, the aortic stenosis was represented by the valve resistance defined by the Ohm's law as the ratio of TPG to Q. However, other investigators (2, 5, 37) have demonstrated that the resistance is not appropriate to describe the LV burden imposed by the stenotic valve. In the present study, we utilized the model proposed by Segers et al. in which we incorporated an analytic model of the transvalvular pressure gradient. This mathematical model of aortic stenosis is based on the energy loss concept as described by Eq. 2, and it has been recently validated in vitro (13). The resulting V3 model is thus represented by a third-order nonlinear differential equation whose resolution can be performed using appropriate initial conditions (see Eqs. 5 and 6). The main advantage of the V3 model is that it may explicitly and accurately describe the behavior of LV volume during ejection if a few (nine) cardiovascular parameters are known.
LV normalized elastance.
The LV elastance normalized with respect to its amplitude and time to peak value is an essential element of the V3 model. Although normalized elastance was similar before and after aortic valve replacement in a given patient, it was different from one patient to the other (Fig. 2). When every patient is considered, the mean absolute difference between the individual normalized elastances and the averaged normalized elastance was 0.053 (maximum = 0.372). The interindividual variability observed in the present study was higher than that reported in the study of Senzaki et al. (30). In their study that included 52 patients, the variability of normalized elastance was minimal despite the presence of interindividual differences with regard to etiology of heart disease, LV myocardial contractility, and loading conditions. However, their study did not include patients with aortic stenosis. Moreover, their results are also in disagreement with those recently published by Kjørstad et al. (15), where a considerable variability of the normalized elastance was observed among the patients. Further studies in a larger number of patients with aortic stenosis are thus needed to assess the interindividual variability of normalized elastance in these patients.
Cardiovascular parameters.
As expected, the systemic vascular resistance estimated by the minimization method was systematically lower than that calculated from experimental data (Fig. 8). Indeed, it can be shown from Eqs. 3 and 7 that Rclin is related to R as follows: Rclin = R + Z0. The aortic characteristic impedance was not measured in the present study. Nonetheless, the value estimated by the minimization method (0.071 ± 0.014 mmHg·s·ml1) was very consistent with that observed in patients by other investigators (19).
Although Segers et al. (29) reported in an animal study that the stroke volume-to-aortic pulse pressure ratio (Eq. 8) systematically underestimates the compliance determined by the area method, Chemla et al. (7) showed that it provides an accurate estimation of the total arterial compliance in patients with a normal aortic valve. Our results (Fig. 9) are consistent with the results of Chemla et al., and they suggest that the stroke volume-to-aortic pulse pressure ratio may also be valid to estimate compliance in patients with aortic stenosis as well as in patients with aortic valve prosthesis. This observation, however, remains to be confirmed in a larger number of patients with aortic stenosis.

View larger version (16K):
[in this window]
[in a new window]
|
Fig. 9. Relationship between compliance (in ml/mmHg) measured using the standard clinical formula (Eq. 8) and that estimated from the minimization method. Dotted line, identity line. Solid and shaded circles correspond to data before and after surgery, respectively.
|
|
Garcia et al. (11) demonstrated that there is a strong agreement between the ELCo and the EOA determined by the Gorlin formula in animals with experimentally induced supravalvular aortic stenosis as well as in patients with aortic valve stenosis. They also showed that ELCo and EOAGorlin are equivalent parameters that accurately reflect the energy loss caused by the valvular stenosis. The good concordance between the ELCo estimated by the minimization method in the present study and the EOA determined by the Gorlin formula (Fig. 10) is consistent with these previous results.
Limitations of the mathematical model.
In the V3 model, we chose a three-element windkessel model even though it has been recently shown by Stergiopulos et al. (33) that a four-element model, including the so-called arterial inertance, is superior in describing the vascular properties. Inserting this fourth element in our coupling model would have resulted in a fourth-order differential equation for the LV volume, which would have therefore necessitated an additional initial condition (see Eq. 6), which is a priori unavailable. This would have therefore resulted in an implicit coupling model that is difficult to apply in the clinical situation. On the other hand, the utilization of the three-element windkessel model provided an explicit model of the LV volume during ejection, which is easily applicable in the context of clinical studies. According to our results, this simplification of our model was not a major limitation, because the estimated compliance and resistance were very consistent with the ones directly obtained from experimental data.
The estimated unloaded volume (V0) was
0 in seven among the twelve datasets (Table 3). This parameter was estimated by the linear approximation of the pressure-volume relation given by Eq. 1. However, the pressure-volume relationship is not linear when LV volume becomes small (14, 21). Hence, V0 has no physical or physiological meaning but is rather a virtual parameter computed by extrapolation. Despite this limitation, the time-varying elastance model remains one of the best simple models of the left ventricle dynamics, and its linear characteristic is appropriate over a large physiological range (14).
Potential clinical implications.
It is often difficult to perform a comprehensive analysis of the interaction between different physiological systems in the context of a clinical study. Indeed, this approach usually requires the realization of a large number of physiological measurements in a large cohort of patients, and these measurements are often difficult or even impossible to perform in patients. The utilization of mathematical models may overcome this problem. Hence, the V3 model validated in the present study may be useful to investigate the interaction between the systemic arterial system, the aortic valve, and the left ventricle in different types of cardiovascular diseases. Moreover, as opposed to other descriptive models previously reported in the literature, the V3 model should be potentially applicable in the clinical setting.
Several recent studies suggested that the so-called degenerative aortic stenosis is probably due to an atherosclerotic process (22). Hence, in the vast majority of the patients, aortic stenosis is probably not a disease solely limited to the valve but rather one manifestation of an atherosclerotic process involving various components of the vascular system including the aorta. Accordingly, Briand et al. (3) have recently reported that 41% of the patients with aortic stenosis concomitantly have an abnormally low systemic arterial compliance. In addition, previous studies (1, 23) have reported that systemic hypertension is a frequent occurrence (3040%) in patients with aortic stenosis. These findings underline that patients with aortic stenosis often have abnormal vascular properties. In such patients, the V3 model may be useful to better assess the respective contributions of the aortic valve and the systemic arterial system to the LV workload in patients having concomitantly aortic stenosis and vascular diseases (i.e., hypertension, atherosclerosis). This could contribute to determine whether the treatment should be targeted on the reduction of the valvular load by aortic valve replacement, on the reduction of the arterial load by medical treatment, or on both aspects. And in this context, the V3 model may be superior to existing simple methods for characterizing the two components of the LV load. For example, the transvalvular pressure gradient and systolic arterial pressure are often used to assess the magnitude of the valvular and vascular loads, respectively. However, these parameters are highly flow dependent and are therefore subject to pseudonormalization in patients with reduced cardiac output (3). This implies that, in such patients, these parameters may underestimate the impact of the aortic valve and/or of the systemic arterial system on the LV workload. The utilization of parameters that are less flow dependent such as the Gorlin EOA or the ELCo for the characterization of valvular properties and the arterial resistance and compliance for the characterization of vascular properties may improve the ability to assess the two components of the LV load. However, when these parameters are analyzed separately, it is difficult to estimate what is the independent contribution of each parameter to the overall LV workload. The most frequently used approach to overcome this limitation in clinical studies is to perform multivariate regression analysis. However, this approach assumes a linear relationship between the different parameters, which is obviously not the case in the ventricular-valvular-vascular system. The utilization of a mathematical model such as the V3 model that relates the aforementioned vascular and valvular parameters with an appropriate formulation may be superior to conventional approaches for the assessment of the respective contributions of the aortic valve and systemic arterial system to the LV workload. This potential clinical usefulness of the V3 model, however, remains to be confirmed in a larger cohort of patients.
 |
APPENDIX
|
|---|
Determination of aortic pressure.
We note t0 as the time at which ejection begins and T as the cardiac period. The solution of Eq. 3 gives the aortic pressure:
where Q* is related to Q as follows:
and the diastolic aortic blood pressure is found using the periodic property of aortic pressure:
Ventricular pressure during the ventricular filling period.
Ventricular pressure during the ventricular filling period is extrapolated using the following expression:
where time t* corresponds to the end of the isovolumic relaxation. In our study, according to the results of Senzaki et al. (30), we postulated that t* = 1.8TEmax. Parameters a and n were calculated such that PV is continuously differentiable. Reminding that time reference (t = 0) is fixed at the onset of the isovolumic contraction, it can be demonstrated that parameters a and n are thus given by:
 |
FOOTNOTES
|
|---|
Address for reprint requests and other correspondence: D. Garcia, Biomedical Engineering Laboratory, IRCM, 110 Pine West Ave., Montreal, Quebec, Canada H2W 1R7 (E-mail: Damien.Garcia{at}ircm.qc.ca)
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
|
|---|
- Antonini-Canterin F, Huang G, Cervesato E, Faggiano P, Pavan D, Piazza R, and Nicolosi GL. Symptomatic aortic stenosis: does systemic hypertension play an additional role? Hypertension 41: 12681272, 2003.[Abstract/Free Full Text]
- Blais C, Pibarot P, Dumesnil JG, Garcia D, Chen D, and Durand LG. Comparison of valve resistance with effective orifice area regarding flow dependence. Am J Cardiol 88: 4552, 2001.[CrossRef][Web of Science][Medline]
- Briand M, Dumesnil JG, Kadem L, Tongue AG, Rieu R, Garcia D, and Pibarot P. Reduced systemic arterial compliance impacts significantly on LV afterload and function in aortic stenosis: implications for diagnosis and treatment. J Am Coll Cardiol. In press.
- Brookes CI, White PA, Bishop AJ, Oldershaw PJ, Redington AN, and Moat NE. Validation of a new intraoperative technique to evaluate load-independent indices of right ventricular performance in patients undergoing cardiac operations. J Thorac Cardiovasc Surg 116: 468476, 1998.[Abstract/Free Full Text]
- Burwash IG, Hay KM, and Chan KL. Hemodynamic stability of valve area, valve resistance, and stroke work loss in aortic stenosis: a comparative analysis. J Am Soc Echocardiogr 15: 814822, 2002.[CrossRef][Web of Science][Medline]
- Chan KL. Is aortic stenosis a preventable disease? J Am Coll Cardiol 42: 593599, 2003.[Abstract/Free Full Text]
- Chemla D, Hebert JL, Coirault C, Zamani K, Suard I, Colin P, and Lecarpentier Y. Total arterial compliance estimated by stroke volume-to-aortic pulse pressure ratio in humans. Am J Physiol Heart Circ Physiol 274: H500H505, 1998.[Abstract/Free Full Text]
- Cohn JN. Arterial compliance to stratify cardiovascular risk: more precision in therapeutic decision making. Am J Hypertens 14: 258S263S, 2001.[CrossRef][Web of Science][Medline]
- Dekker AL, Barenbrug PJ, van der Veen FH, Roekaerts P, Mochtar B, and Maessen JG. Pressure-volume loops in patients with aortic stenosis. J Heart Valve Dis 12: 325332, 2003.[Web of Science][Medline]
- Fitchett DH. LV-arterial coupling: interactive model to predict effect of wave reflections on LV energetics. Am J Physiol Heart Circ Physiol 261: H1026H1033, 1991.[Abstract/Free Full Text]
- Garcia D, Dumesnil JG, Durand LG, Kadem L, and Pibarot P. Discrepancies between catheter and Doppler estimates of valve effective orifice area can be predicted from the pressure recovery phenomenon: practical implications with regard to quantification of aortic stenosis severity. J Am Coll Cardiol 41: 435442, 2003.[Abstract/Free Full Text]
- Garcia D, Pibarot P, Dumesnil JG, Sakr F, and Durand LG. Assessment of aortic valve stenosis severity: a new index based on the energy loss concept. Circulation 101: 765771, 2000.[Abstract/Free Full Text]
- Garcia D, Pibarot P, and Durand LG. Analytical modeling of the instantaneous pressure gradient across the aortic valve. J Biomech doi:101016/jjbiomech.2004.06.018: 2005.
- Kass DA and Maughan WL. From "Emax" to pressure-volume relations: a broader view. Circulation 77: 12031212, 1988.[Free Full Text]
- Kjorstad KE, Korvald C, and Myrmel T. Pressure-volume-based single-beat estimations cannot predict left ventricular contractility in vivo. Am J Physiol Heart Circ Physiol 282: H1739H1750, 2002.[Abstract/Free Full Text]
- Kono A, Maughan WL, Sunagawa K, Hamilton K, Sagawa K, and Weisfeldt ML. The use of left ventricular end-ejection pressure and peak pressure in the estimation of the end-systolic pressure-volume relationship. Circulation 70: 10571065, 1984.[Abstract/Free Full Text]
- Lagarias J, Reeds J, Wright M, and Wright P. Convergence properties of the Nelder-Mead simplex method in low dimensions. SIAM J Optimiz 9: 112147, 1998.[CrossRef]
- Li JK, Zhu JY, and Nanna M. Computer modeling of the effects of aortic valve stenosis and arterial system afterload on left ventricular hypertrophy. Comput Biol Med 27: 477485, 1997.[CrossRef][Web of Science][Medline]
- Merillon JP, Fontenier GJ, Lerallut JF, Jaffrin MY, Motte GA, Genain CP, and Gourgon RR. Aortic input impedance in normal man and arterial hypertension: its modification during changes in aortic pressure. Cardiovasc Res 16: 646656, 1982.[Web of Science][Medline]
- Nishimura RA. Cardiology patient pages. Aortic valve disease. Circulation 106: 770772, 2002.[Free Full Text]
- Noda T, Cheng CP, De Tombe PP, and Little WC. Curvilinearity of LV end-systolic pressure-volume and dP/dtmax-end-diastolic volume relations. Am J Physiol Heart Circ Physiol 265: H910H917, 1993.[Abstract/Free Full Text]
- Novaro GM and Griffin BP. Calcific aortic stenosis: another face of atherosclerosis? Cleve Clin J Med 70: 471477, 2003.[Abstract/Free Full Text]
- Pate GE. Association between aortic stenosis and hypertension. J Heart Valve Dis 11: 612614, 2002.[Web of Science][Medline]
- Pibarot P and Dumesnil JG. Hemodynamic and clinical impact of prosthesis-patient mismatch in the aortic valve position and its prevention. J Am Coll Cardiol 36: 11311141, 2000.[Abstract/Free Full Text]
- Rossebo AB and Pedersen TR. Hyperlipidaemia and aortic valve disease. Curr Opin Lipidol 15: 447451, 2004.[CrossRef][Web of Science][Medline]
- Segers P, Steendijk P, Stergiopulos N, and Westerhof N. Predicting systolic and diastolic aortic blood pressure and stroke volume in the intact sheep. J Biomech 34: 4150, 2001.[CrossRef][Web of Science][Medline]
- Segers P, Stergiopulos N, and Westerhof N. Quantification of the contribution of cardiac and arterial remodeling to hypertension. Hypertension 36: 760765, 2000.[Abstract/Free Full Text]
- Segers P, Stergiopulos N, and Westerhof N. Relation of effective arterial elastance to arterial system properties. Am J Physiol Heart Circ Physiol 282: H1041H1046, 2002.[Abstract/Free Full Text]
- Segers P, Verdonck P, Deryck Y, Brimioulle S, Naeije R, Carlier S, and Stergiopulos N. Pulse pressure method and the area method for the estimation of total arterial compliance in dogs: sensitivity to wave reflection intensity. Ann Biomed Eng 27: 480485, 1999.[CrossRef][Web of Science][Medline]
- Senzaki H, Chen CH, and Kass DA. Single-beat estimation of end-systolic pressure-volume relation in humans. A new method with the potential for noninvasive application. Circulation 94: 24972506, 1996.[Abstract/Free Full Text]
- Smith BW, Chase JG, Nokes RI, Shaw GM, and Wake G. Minimal haemodynamic system model including ventricular interaction and valve dynamics. Med Eng Phys 26: 131139, 2004.[CrossRef][Web of Science][Medline]
- Stergiopulos N, Meister JJ, and Westerhof N. Determinants of stroke volume and systolic and diastolic aortic pressure. Am J Physiol Heart Circ Physiol 270: H2050H2059, 1996.[Abstract/Free Full Text]
- Stergiopulos N, Westerhof BE, and Westerhof N. Total arterial inertance as the fourth element of the windkessel model. Am J Physiol Heart Circ Physiol 276: H81H88, 1999.[Abstract/Free Full Text]
- Suga H. Cardiac energetics: from E(max) to pressure-volume area. Clin Exp Pharmacol Physiol 30: 580585, 2003.[CrossRef][Web of Science][Medline]
- 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 on the ratio. Circ Res 32: 314322, 1973.[Abstract/Free Full Text]
- Ursino M. Interaction between carotid baroregulation and the pulsating heart: a mathematical model. Am J Physiol Heart Circ Physiol 275: H1733H1747, 1998.[Abstract/Free Full Text]
- Voelker W, Reul H, Nienhaus G, Stelzer T, Schmitz B, Steegers A, and Karsch KR. Comparison of valvular resistance, stroke work loss, and Gorlin valve area for quantification of aortic stenosis. An in vitro study in a pulsatile aortic flow model. Circulation 91: 11961204, 1995.[Abstract/Free Full Text]
- Westerhof N, Elzinga G, and Sipkema P. An artificial arterial system for pumping hearts. J Appl Physiol 31: 776781, 1971.[Free Full Text]
- Wilkinson JL. Haemodynamic calculations in the catheter laboratory. Heart 85: 113120, 2001.[Free Full Text]
- Willens HJ, Davis W, Herrington DM, Wade K, Kesler K, Mallon S, Brown WV, Reiber JH, and Raines JK. Relationship of peripheral arterial compliance and standard cardiovascular risk factors. Vasc Endovascular Surg 37: 197206, 2003.[Abstract/Free Full Text]
- Zacek M and Krause E. Numerical simulation of the blood flow in the human cardiovascular system. J Biomech 29: 1320, 1996.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:

|
 |

|
 |
 
W. K. Laskey, W. G. Kussmaul III, and A. Noordergraaf
Systemic Arterial Response to Exercise in Patients With Aortic Valve Stenosis
Circulation,
February 24, 2009;
119(7):
996 - 1004.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
D. Garcia, P. G. Camici, L.-G. Durand, K. Rajappan, E. Gaillard, O. E. Rimoldi, and P. Pibarot
Impairment of coronary flow reserve in aortic stenosis
J Appl Physiol,
January 1, 2009;
106(1):
113 - 121.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
D. Tanne, L. Kadem, R. Rieu, and P. Pibarot
Hemodynamic impact of mitral prosthesis-patient mismatch on pulmonary hypertension: an in silico study
J Appl Physiol,
December 1, 2008;
105(6):
1916 - 1926.
[Abstract]
[Full Text]
[PDF]
|
 |
|
Copyright © 2005 by the American Physiological Society.