AJP - Heart  AJP: Regulatory, Integrative and Comparative Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 288: H1740-H1746, 2005. First published November 18, 2004; doi:10.1152/ajpheart.00584.2004
0363-6135/05 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
288/4/H1740    most recent
00584.2004v1
Right arrow Alert me when this article is cited
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 Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Mohri, S.
Right arrow Articles by Suga, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Mohri, S.
Right arrow Articles by Suga, H.

Normal distribution of ventricular pressure-volume area of arrhythmic beats under atrial fibrillation in canine heart

Satoshi Mohri,1 Juichiro Shimizu,1 Gentaro Iribe,1 Haruo Ito,1 Terumasa Morita,2 Hiroki Yamaguchi,2 Shunji Sano,2 Fumihiko Kajiya,1 and Hiroyuki Suga1,3

Departments of 1Cardiovascular Physiology and 2Cardiovascular Surgery, Okayama University Graduate School of Medicine and Dentistry, Okayama; and 3National Cardiovascular Center Research Institute, Osaka, Japan

Submitted 14 June 2004 ; accepted in final form 16 November 2004


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
We previously found the frequency distribution of the left ventricular (LV) effective afterload elastance (Ea) of arrhythmic beats to be nonnormal or non-Gaussian in contrast to the normal distribution of the LV end-systolic elastance (Emax) in canine in situ LVs during electrically induced atrial fibrillation (AF). These two mechanical variables determine the total mechanical energy [systolic pressure-volume area (PVA)] generated by LV contraction when the LV end-diastolic volume is given on a per-beat basis. PVA and Emax are the two key determinants of the LV O2 consumption per beat. In the present study, we analyzed the frequency distribution of PVA during AF by its {chi}2, significance level, skewness, and kurtosis and compared them with those of other major cardiodynamic variables including Ea and Emax. We assumed the volume intercept (V0) of the end-systolic pressure-volume relation needed for Emax determination to be stable during arrhythmia. We found that PVA distributed much more normally than Ea and slightly more so than Emax during AF. We compared the {chi}2, significance level, skewness, and kurtosis of all the complex terms of the PVA formula. We found that the complexity of the PVA formula attenuated the effect of the considerably nonnormal distribution of Ea on the distribution of PVA along the central limit theorem. We conclude that mean (SD) of PVA can reliably characterize the distribution of PVA of arrhythmic beats during AF, at least in canine hearts.

end-systolic elastance; effective afterload elastance; frequency distribution; normality


MANY GROUPS, including ours, investigated the statistical characteristics of cardiodynamic variables during ventricular arrhythmia under atrial fibrillation (AF) in basic and clinical studies (15, 713, 15, 23, 24). We recently found that the frequency distribution of the left ventricular (LV) effective afterload elastance (Ea; Fig. 1) (16, 21) is nonnormal or non-Gaussian, with considerable skewness and kurtosis in canine hearts under AF, in contrast to the LV contractility index (Emax; Fig. 1), with relatively normal or Gaussian frequency distribution (11, 12, 24). Therefore, the mean (SD) of Ea alone could not reliably represent its nonnormal distribution, with considerable skewness and kurtosis, and may lead to an incorrect evaluation of the LV mechanical performance during AF (11, 12).



View larger version (63K):
[in this window]
[in a new window]
 
Fig. 1. Left ventricular (LV) pressure-volume (P-V) diagram showing relationships of P-V area (PVA) to other major cardiodynamic variables. ESP, end-systolic pressure; Ea, effective afterload elastance; Emax, end-systolic ventricular elastance; PE, potential energy (hatched area); EW, external work (vertically shaded area); PVA, systolic PVA as a measure of total mechanical energy generated by LV contraction; SV, stroke volume; V0, dead volume at which systolic pressure is zero; ESV, end-systolic volume; EDV, end-diastolic volume; {Delta}PVA (white triangular area), 3-ms increment ({Delta}t) in PVA.

 
We previously found that Emax and Ea are the two major determinants of the total mechanical energy of LV contraction that is measurable as the systolic pressure-volume area (PVA; Fig. 1) (18, 19). The PVA is the sum of the LV external mechanical work (EW) and the elastic potential energy (PE; Fig. 1). Then we wondered whether normally or nonnormally the PVA determined by the combination of the nonnormally distributing Ea and the normally distributing Emax would distribute. If PVA distributes normally, its mean (SD) could fully characterize its distribution. However, if PVA distributes nonnormally, its mean (SD) alone could not fully characterize its distribution.

In the present study, we investigated the frequency distribution of PVA and compared it with that of other major cardiodynamic variables including Emax (16, 18, 19) and Ea (16, 21) during AF in the canine in situ ejecting heart. The present experimental data showed that PVA distributed much more normally, with less skewness and kurtosis, than Ea and even slightly more normally than Emax. The PVA is determined predominantly by Emax and Ea when the end-diastolic volume (EDV) of each arrhythmic beat is given (16, 18). Here, Emax is primarily determined by the preceding beat intervals (24), Ea is determined by the ventricular-aortic coupling in each arrhythmic beat (16), and EDV is determined by the beat intervals, Emax, and Ea of the preceding beats. We found the normal distribution of PVA to be attributable to the complexity of the PVA formula as a function of Emax and Ea. The present result indicates that the mean (SD) of PVA of arrhythmic beats of the ventricle could reasonably characterize the distribution of PVA in the same manner as Emax.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Surgical preparation. We performed the canine experiments in conformity with the guiding principles for the care and use of animals in the field of physiological sciences of the Japan Physiological Society and the American Physiological Society. Six mongrel dogs (6–12 kg) were anesthetized with pentobarbital sodium (25 mg/kg iv) after premedication with ketamine hydrochloride (50 mg/kg im) and intubated in each experiment. The anesthesia was maintained by fentanyl (100 µg/h per dog iv) as usual in our laboratory (11, 12, 24).

The chest of the dog was opened midsternally. A 3-F catheter-tipped micromanometer was inserted into the LV from the apex to measure the LV pressure (LVP). A 7-F eight-electrode conductance catheter (Webster Laboratories, Baldwin Park, CA) was introduced into the LV through an apical stab and placed along the ventricular long axis to measure LV volume (LVV) continuously. The method for measuring LVV with this catheter was described in detail previously (11, 12, 24).

Briefly, the catheter measured continuously the time-varying electrical conductance [Gi(t)] of the five segments (i = 1–5) of blood in the LV cavity. The LV total blood volume was then continuously calculated from the five segmental Gi(t) (i = 1–5) after calibration of blood conductivity in the sampling cuvette. Our custom-made signal conditioner-processor (SI Medicotech) was used to convert the segmental Gi(t) to LV conductance volume. The parallel conductance (Gp) due to the conductance of the LV wall and the surrounding tissues and fluid was obtained by the standard hypertonic saline-dilution method. A constant offset volume (Vc) was calculated from Gp. The absolute LVV was obtained by subtracting Vc from the LV conductance converted from the segmental Gi(t).

We determined LV volume-axis intercept (V0) of the end-systolic pressure-volume (P-V) relation (ESPVR) by clamping the inferior vena cava during 10–20 regular beats, as described elsewhere (11, 12, 14, 16, 18, 21, 24). The clamp gradually decreased LVP and LVV and shifted the shrinking LV P-V loop toward the origin of the P-V diagram. We drew a straight line tangential to the left-upper or end-systolic corners of these P-V loops and extrapolated the line to the volume axis to obtain its V0 (11, 12, 24).

We assumed that this V0 was reasonably stable, even in individual arrhythmic beats, because there was no direct method to obtain V0 in every arrhythmic beat with a variable Emax. We also knew that V0 decreased and increased by a few milliliters with considerable increases and decreases, respectively, in Emax in regular beats (20). We had also confirmed that these changes in V0 could not significantly affect Emax and PVA values to the extent that they mislead the cardiac mechanoenergetics (16).

We attached a pair of stimulation electrodes to the left atrial appendage. Suprathreshold electrical stimulation at 20 Hz via these electrodes induced and maintained AF (11, 12, 24). AF started on and stopped off the stimulation. We produced AF for 2 min and recorded LVP and LVV during the latter 1 min at sampling intervals of 3 ms in a computer.

PVA. The systolic PVA quantifies the total mechanical energy generated by a ventricular contraction (18, 19). Figure 1 illustrates PVA in the ventricular P-V diagram. PVA consists of two areas: the EW (rectangular area) within the P-V loop and the elastic PE (triangular area) between the end-systolic P-V or Emax line and the end-diastolic P-V curve on the origin side of the P-V loop.

We calculated PVA of each contraction by summing all the PVA increments per 3 s from end diastole to end systole. Each 3-ms ({Delta}t) PVA increment ({Delta}PVA) corresponds to the narrow triangular area scanned by a line connecting the instantaneous P-V data point drawing the P-V loop and V0 on the volume axis (16) (Fig. 1). We obtained EW by summing {Delta}PVA from the onset of the present contraction to the onset of the next contraction (16). We obtained PE by subtracting EW from PVA (16).

We further obtained the contractility index, Emax, of each contraction as the slope of the ESPVR (Fig. 1) (16, 18, 19). We also obtained the effective Ea of each contraction as the slope of the line connecting the end-systolic P-V point on the left upper corner of the P-V loop and the EDV on the LVV axis (Fig. 1) (16, 1821).

We know that PVA is formulated approximately by Eq. A10 as explained in the APPENDIX. Equation A10 indicates that PVA is a complex function of Emax, Ea, and EDV.

Statistics. We performed the basic statistical analyses of PVA as well as Emax, Ea, beat interval (RR), EW, PE, EDV, and stroke volume (SV) of all the arrhythmic beats sampled in an arbitrary 1 min in each heart. We studied their frequency distributions and performed {chi}2 test of their normality, evaluating it by its P value (Table 1; see Figs. 4 and 5) (6, 17).


View this table:
[in this window]
[in a new window]
 
Table 1. Basic statistics of PVA during atrial fibrillation in canine hearts

 


View larger version (26K):
[in this window]
[in a new window]
 
Fig. 4. A representative set of histograms of frequency distributions of PVA (A) and Ea (B) in heart 1 (Table 1). Normal distribution curves are the best-fit curves with the least squares. Note closeness of the histogram to the normal curve for PVA and their separation for Ea.

 


View larger version (31K):
[in this window]
[in a new window]
 
Fig. 5. Means ± SDs of {chi}2 (A), its P value (B), skewness (C), and kurtosis (D) of major cardiodynamic variables in 6 hearts. RR, beat interval. Solid arrows indicate significance by Bonferroni's t-test and least significant difference (LSD) method. Dotted arrows indicate significance by LSD method only. Values in inset represent P < 0.05 for significance by the respective multiple comparisons.

 
We also obtained two dimensionless measures of departure from the normality or Gaussianity of the frequency distributions, i.e., skewness and kurtosis, as in our previous study (Table 1; see Figs. 4 and 5) (11, 12). The normal distribution has zero skewness and zero kurtosis. A positive skewness indicates a leftward shift of the peak frequency, and a negative skewness indicates its rightward shift. A positive kurtosis indicates a sharper peak, and a negative kurtosis indicates a dull peak or even two separate peaks.

We analyzed correlations among Emax, Ea, EDV, and PVA and obtained their correlation coefficients (r) and coefficients of determination (r2; Table 2). We also obtained autocorrelation coefficients of EDV, Emax, Ea, and PVA with a lag of only one beat interval, e.g., correlation coefficient between EDVi and EDVi + 1 where i was increased from 1 to the one beat before the last beat in the sampled 1 min (Table 3).


View this table:
[in this window]
[in a new window]
 
Table 2. Correlation coefficients and coefficients of determination among Emax, Ea, EDV, and PVA in canine hearts

 

View this table:
[in this window]
[in a new window]
 
Table 3. Autocorrelation coefficients of Emax, Ea, EDV, and PVA between two adjacent arrhythmic beats in 1 min in six canine hearts

 
We compared the basic statistics of EmaxEa2, Emax2Ea, Emax2Ea + EmaxEa2/2, (EDV – V0)2, (Emax + Ea)2, Emax + Ea, and (Emax2Ea + EmaxEa2/2)/(Emax + Ea)2, which are the variously complex terms of Eq. A10, to calculate PVA from EDV, Emax, and Ea (see APPENDIX). We also compared the basic statistics between PVA calculated by Eq. A10 and the measured PVA (see Fig. 6).



View larger version (30K):
[in this window]
[in a new window]
 
Fig. 6. Means ± SDs of {chi}2 (A), its P value (B), skewness (C), and kurtosis (D) of terms in Eq. A10 to calculate PVA from Emax, Ea, and EDV. Abscissas show terms 1–8 arranged in decreasing order of {chi}2, as shown in inset. Solid arrows indicate significance by Bonferroni's t-test and LSD method. Dotted arrows indicate significance by LSD method only. Values in inset represent P < 0.05 for significance by the respective multiple comparisons.

 
We obtained moving-averaged PVA values by increasing the number of moving-averaged beats from 1 to 50 (see Fig. 7).



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 7. Decreasing percent coefficient of variation (%CV = 100 x SD/mean) of moving-averaged PVA with an increasing number of moving-averaged beats in 6 hearts. Ordinate shows %CV of moving-averaged values. Middle (solid) curve, mean; top and bottom (dashed) curves, mean ± SD. Horizontal dotted lines indicate 5% and 10% CV measured PVA. Values in inset show P < 0.05 for significance. Solid arrows indicate significance by Bonferroni's t-test and LSD method. Dotted arrows indicate significance by LSD method only.

 
In these statistical analyses, we used Microsoft Office 2001 Excel and StatView version 4.5. Values are means (SD).


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Figure 2 shows a representative 3-s segment of the continuous traces of LVP and LVV during AF. It contains nine arrhythmic beats separated by various beat intervals and featured by all varying peak LVPs and end-diastolic and end-systolic LVVs. All six hearts showed similar arrhythmic beats during AF.



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 2. An arbitrarily chosen 3-s segment of the continuous traces of LV pressure (LVP, thick traces) and volume (LVV, thin traces) of several arrhythmic beats during electrically induced atrial fibrillation (AF) in a canine in situ heart.

 
Figure 3 shows the same set of LV P-V loops of the four consecutive arrhythmic beats during AF and their PVA as well as Emax and Ea (16, 1821). These correspond to the first four beats in Fig. 2. Figure 3A illustrates PVA of the arrhythmic beat with the highest contractility (largest Emax) and largest SV of the four beats. PVA consists of the sum of the rectangular (shaded) area within the P-V loop and the triangular (hatched) area between the ESPVR line and the EDPVR curve on the origin side of the P-V loop of this contraction. The former area represents EW performed by the LV during this contraction. The latter area represents PE that was generated by this contraction in the LV wall but was not converted to EW (18, 19). Figure 3B illustrates PVA of the arrhythmic beat with the lowest contractility (smallest Emax) and the smallest SV of the four beats. In this beat, EW was nearly zero.



View larger version (30K):
[in this window]
[in a new window]
 
Fig. 3. LV P-V loops of 4 consecutive arrhythmic beats during AF. Beats correspond to the first 4 beats in Fig. 2. A: PVA, Emax, and Ea of the strongest contraction of the 4 beats. B: PVA, Emax, and Ea of the weakest contraction of the same 4 beats. ESPVR, end-systolic P-V relation, which has slope Emax; EW, external work or stroke work; PE, potential energy (PVA = PE + EW); SV, stroke volume; V0, volume axis intercept of ESPVR, where end-systolic volume is zero. Note considerable changes in all cardiodynamic variables between beats.

 
Figure 4 shows a representative set of the frequency distributions of PVA and Ea in one heart (heart 1 in Table 1). PVA distributed much more normally than Ea. We chose Ea to compare with PVA, because we had found Ea to have a considerable nonnormal distribution (11). The {chi}2, its P value, skewness, and kurtosis quantify the difference of the frequency distribution between PVA and Ea. The mean, median, and mode scattered less for PVA than for Ea. The other hearts showed similar results (Table 1). For PVA, no P value of {chi}2 test was <0.05; the smallest P value was 0.162. Neither skewness (P > 0.661) nor kurtosis (P > 0.116) was significantly different from 0.

Figure 5 compares means (SD) of {chi}2, its P value, skewness, and kurtosis of not only Ea and PVA but also other cardiodynamic variables (RR, Emax, EW, PE, SV, and EDV) observed in all six hearts. The order of the cardiodynamic variables on the abscissas was the increasing order of {chi}2 values from left to right. Ea had the largest {chi}2, skewness, and kurtosis and the smallest P value for {chi}2 test. On the other hand, PVA as well as EDV and SV had much smaller {chi}2, skewness, and kurtosis and larger P value for {chi}2 test. For Emax, {chi}2 and P values were intermediate but skewness and kurtosis were small.

The one-way repeated-measures analysis of variance (ANOVA) showed statistical significance of F-test for {chi}2, its P value, skewness, and kurtosis among the cardiodynamic variables on the abscissas. The inset P values for {chi}2 test on the solid horizontal arrows in Fig. 5 indicate that Bonferroni's t-test and least significant difference (LSD) method showed statistical significance between the variables bridged by the arrows. The P values on the dotted arrows indicate statistical significance only by LSD.

The results in Fig. 5, A and B, indicate that PVA, as well as Emax, PE, EW, SV, and EDV, had significantly smaller {chi}2 values with significantly higher P values for {chi}2 test than Ea. RR had intermediate {chi}2 and P values between Ea and either Emax or PVA. PVA as well as EDV, SV, and PE had significantly smaller {chi}2 values with significantly higher P values for {chi}2 test than Ea and RR. The P values for {chi}2 test of PVA was >0.4 in five hearts and 0.1 in one heart (Table 1), negating the null hypothesis of the nonnormal distribution of PVA and, hence, supporting its normal distribution. However, the P values for {chi}2 test of Ea were <0.05 in three hearts, supporting the null hypothesis of the nonnormal distribution of Ea and, hence, negating its normal distribution. The other variables had P values for {chi}2 test >0.1. This negated the null hypothesis of their nonnormal distributions.

Moreover, the results in Fig. 5, C and D, indicate that PVA, as well as EDV, SV, EW, Emax, and RR, had much smaller skewness and kurtosis than Ea and PE. PVA, as well as EW and Emax, had virtually zero skewness on average. PVA, as well as SV, EW, and Emax, had small negative kurtosis on average. Although we did not include end-systolic pressure (ESP) and end-systolic volume (ESV) in Fig. 5, their small {chi}2, skewness, and kurtosis were comparable to those of PVA and EDV and much smaller than those of Ea.

Figure 6 compares means (SD) of {chi}2, its P value, skewness, and kurtosis of the various complex terms of Eq. A10 to determine PVA from EDV, Emax, and Ea (see APPENDIX) in all six hearts. Terms 1–7 (Fig. 6) are arranged in decreasing order of {chi}2 values from left to right. Term 8 is the PVA calculated by Eq. A10, and term 9 is the measured PVA.

The calculated PVA (term 8) had relatively small {chi}2 and large P (>0.155, all insignificant). However, its skewness and kurtosis were closest to zero. These statistical values were comparable to those of the measured PVA (term 9). The measured PVA data used here are the same as those used in Fig. 5.

The results in Fig. 6 indicate that EmaxEa2 (term 1) had the most nonnormal distribution with the largest {chi}2, the smallest P value (<0.05 only in 2 hearts), and the largest skewness and kurtosis. However, term 7 (Emax2Ea + EmaxEa2/2)/(Emax + Ea)2 had the smallest {chi}2, the largest P value (>0.772), and small, but not the smallest, skewness and kurtosis, although it contains multiple Ea values.

Table 2 shows r and r2 values among Emax, Ea, and EDV in Eq. A10 and the calculated PVA. They were positively or negatively correlated significantly in all six hearts. The r2 values indicate that the variation of any of Emax, Ea, EDV, and PVA was attributable to those of the other three variables by 30–80%. This result indicates that Emax, Ea, EDV, and PVA were variably, although not fully, independent of each other.

Table 3 shows the autocorrelation coefficients (r) of EDV, Emax, Ea, and PVA with a lag of only one beat. Most r values were insignificant. Therefore, these variables changed on a per-arrhythmic-beat basis, rather randomly, independent of their respective values of the last beat in most of the six hearts.

Figure 7 shows that the percent coefficient of variation (CV) of moving-averaged PVA decreased with an increasing number of moving-averaged beats in the six hearts. The percent CV of ~35% on average of raw PVA values during AF decreased rapidly to ~10% and further to ~5% as a result of an increase in the number of moving-averaged beats to 20 and 40 beats, respectively.


    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
The present results have evidently shown for the first time the normality or Gaussianity of the frequency distribution of PVA of arrhythmic beats in canine LVs under electrically induced AF (Figs. 4 and 5). This normality existed despite the significantly nonnormal or non-Gaussian frequency distribution of Ea of the same arrhythmic beats. Here, Ea is one of the major determinants of PVA. Because PVA is a measure of the total mechanical energy of each ventricular contraction (18, 19), its normal distribution indicates that the ventricular total mechanical energy of each arrhythmic beat distributes normally during AF, despite the considerably nonnormally distributing Ea (11, 12, 24).

The present results have also shown that although Ea (Fig. 4) and EmaxEa2 (Fig. 6) in the formula (Eq. A10) to calculate PVA distributed nonnormally, the other terms of the formula distributed normally to the comparable extent of the calculated PVA as well as the measured PVA. Of these, Emax and EDV have been determined by the onset of each contraction as the results of the cardiodynamic variables of the preceding arrhythmic beats (22, 24). Ea is defined as ESP/SV in each contraction, where ESP and SV are determined as the result of the ventricular-aortic interaction during each contraction (14, 16, 21).

Therefore, PVA is determined as the consequence of a complex interaction of various cardiodynamic variables with normal or nonnormal frequency distributions (3, 4, 7, 913, 24). The normal distribution of PVA, therefore, seems principally attributable to the central limit theorem or the law of large numbers (6). This theorem assumes all the component variables to be mutually independent. However, the component variables, i.e., EDV, Emax, and Ea, of PVA correlated mutually by 30–80% in terms of r2 (Table 2). This may partly account for some, although not significant, residual {chi}2 of PVA, despite its virtually zero skewness and kurtosis (Table 1, Fig. 6). The present study for the first time revealed the normal or Gaussian frequency distribution of PVA during AF and its underlying mechanism, although in the LV of normal canine hearts.

The present study has also shown that the mean of PVA per beat during AF could be estimated reliably with a reasonably small (10–5%) CV of the mean by averaging PVA values of 20–40 beats, respectively (Fig. 7). This indicates that the mean (SD) of PVA of only 20–40 beats could give us a reasonably reliable estimation of the frequency distribution of PVA of arrhythmic beats during AF, at least in normal canine hearts. This suggests that only 7–15 s are necessary to estimate reasonable mean (SD) of PVA under AF. This rather rapid convergence of the SD of moving-averaged PVA seems partly attributable to the largely insignificant autocorrelation of PVA between the present and preceding beats as the result of the similarly insignificant autocorrelation of the components (EDV, Emax, and Ea) of PVA (Table 3).

Although we analyzed neither LV myocardial Ca2+ handling nor cross-bridge cycling in the present study, we know that an increased EDV increases Ca2+ handling and cross-bridge cycling according to the Starling law of the heart (16, 18). An increased Emax manifests an LV contractility enhanced with increased Ca2+ handling and recruitability of cross-bridge cycling (16, 18). Here, myocardial Ca2+ handling is determined by its restitution and potentiation over the last two preceding beats (24), whereas cross-bridge cycling is primarily determined by EDV, Emax, and Ea (16, 18). Ea can directly and simply be related to neither myocardial Ca2+ handling nor cross-bridge cycling, because it is defined as ESP/SV (11, 21), where ESP and SV are determined by EDV, Emax, and total peripheral resistance (16, 21). These cardiac system-element interrelations suggest that PVA of an arrhythmic beat is determined in a complex manner by not only element, but also system, factors. Although we do not know the statistical characteristics of all the system-element factors involved in PVA determination, we could speculate that the central limit theorem may serve for the normal distribution of PVA.

An obvious limitation of the present study is the use of normal canine hearts with electrically induced AF. It is possible that statistical characteristics of PVA during AF in human diseased hearts are different from those we have found in the present study (14, 10). This warrants future studies to investigate the statistical characteristics of PVA in various pathological hearts under artificial or spontaneous AF. Their statistical characteristics may reflect the pathophysiological conditions of the heart.

One more limitation of the present study is the assumption of a constant V0 under arrhythmia with a widely varying Emax. Fortunately, our previous studies showed that changes in V0, if any, with a varying Emax would be relatively small, e.g., a few milliliters within a physiologically working range of Emax (16, 20). Any change in Emax affects PVA via PE (Figs. 1 and 3, Eqs. A1 and A3 in APPENDIX). For example, even a 2-ml decrease in V0 with an increase in Emax from 2.5 to 10 mmHg/ml could increase PE only by an additional triangular area bound by the end-systolic P-V point and the new 2-ml base on the volume axis in the P-V diagram, though not shown graphically. When the end-systolic pressure increases representatively from 50 to 100 mmHg with the increased Emax (Figs. 2 and 3), the decreased V0 with the increased Emax would increase PVA by only 50–100 mmHg/ml. This increase in PVA seems too small to deviate the frequency distribution of PVA from normality (e.g., Fig. 4A).

We therefore conclude that the total mechanical energy generated by arrhythmic beats during AF can reliably be characterized by the mean (SD) of PVA in the normal canine LV. The normal distribution of PVA during AF seems theoretically attributable to the complexity of the PVA formula as a function of Ea, Emax, and EDV of individual arrhythmic beats.


    APPENDIX
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Figure 1 helps in our understanding of the following relations of PVA to the other cardiodynamic variables, including Ea and Emax, in the LV P-V diagram (16).

First, by definition,

(A1)
These EW and PE are formulated as follows

(A2)

(A3)
Here, we assumed that LVP during ejection was constant at ESP and the end-diastolic P-V relation was virtually zero below ESV. This assumption does not disadvantageously affect the magnitude of PVA, except in the presence of aortic stenosis or regurgitation.

ESV and ESP are formulated by the following simultaneous equations

(A4)

(A5)
Solving simultaneous Eqs. A4 and A5 yields

(A6)

(A7)
Substituting Eqs. A6 and A7 into Eqs. A2 and A3 yields

(A8)

(A9)
Substituting Eqs. A8 and A9 into Eq. A1 yields

(A10)
Equation A10 indicates that PVA is a complex function of Emax, Ea, and EDV, all of which change beat-by-beat with variable magnitudes of SD, skewness, and kurtosis (Fig. 5). The complexity of the PVA equation (Eq. A10) includes additions, multiplications, squares, and divisions of Emax, Ea, and EDV.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
This work was partly supported by Ministry of Education, Culture, Sports, Science, and Technology (Japan) Scientific Research Grants 13558113, 13770350, 13854030, 13878185, 13878192, 14380405, 15650095, 15659186, 15700330, and 16659057 and Ministry of Health, Labour, and Welfare (Japan) Cardiovascular Diseases Research Grant 14A-1.


    FOOTNOTES
 

Address for reprint requests and other correspondence: S. Mohri, Dept. of Cardiovascular Physiology, Okayama Univ. Graduate School of Medicine and Dentistry, 2-5-1 Shikatacho, Okayama 700-8558, Japan (E-mail: smohri{at}md.okayama-u.ac.jp)

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
 APPENDIX
 GRANTS
 REFERENCES
 

  1. Atwood JE, Myers J, Sandhu S, Lachterman B, Friis R, Oshita A, Forbes S, Walsh D, and Froelicher V. Optimal sampling interval to estimate heart rate at rest and during exercise in atrial fibrillation. Am J Cardiol 63: 45–48, 1989.[CrossRef][Web of Science][Medline]
  2. Bootsma BK, Hoelen AJ, Strackerr J, and Meijler FL. Analysis of RR intervals in patients with atrial fibrillation at rest and during exercise. Circulation 41: 783–793, 1970.[Abstract/Free Full Text]
  3. Brookes CI, White PA, Staples M, Oldershaw PJ, Redington AN, Collins PD, and Noble MI. Myocardial contractility is not constant during spontaneous atrial fibrillation in patients. Circulation 98: 1762–1768, 1998.[Abstract/Free Full Text]
  4. Daoud EG, Weiss R, Bahu M, Knight BP, Bogun F, Goval R, Harvey M, Strickberger SA, Man KC, and Morady F. Effect of an irregular ventricular rhythm on cardiac output. Am J Cardiol 78: 1433–1436, 1996.[CrossRef][Web of Science][Medline]
  5. Dubrey SW and Falk RH. Optimal number of beats for the Doppler measurement of cardiac output in atrial fibrillation. J Am Soc Echocardiogr 10: 67–71, 1997.[CrossRef][Web of Science][Medline]
  6. Glantz SA. Primer of Biostatistics. New York: McGraw-Hill, 1981, p. 217–221, 284–292.
  7. Hardman SMC, Noble MIM, and Seed WA. Postextrasystolic potentiation and its contribution to the beat-to-beat variation of the pulse during atrial fibrillation. Circulation 86: 1223–1232, 1992.[Abstract/Free Full Text]
  8. Hoopen MT and Bongaarts JPM. Probabilistic characterization of R-R intervals. Cardiovasc Res 3: 218–226, 1969.[Web of Science][Medline]
  9. Horan LG and Kistler JC. Study of ventricular response in atrial fibrillation. Circ Res 9: 305–312, 1961.[Abstract/Free Full Text]
  10. Kerr AJ, Simmonds MB, and Steward RA. Influence of heart rate on stroke volume variability in atrial fibrillation in patients with normal and impaired left ventricular function. Am J Cardiol 82: 1496–1500, 1998.[CrossRef][Web of Science][Medline]
  11. Mitani H, Yamaguchi H, Morita T, Oshima Y, Shimizu J, Ito H, Araki J, Takaki M, Sano S, and Suga H. Effective arterial elastance of irregular beats during atrial fibrillation in canine left ventricle. Jpn J Physiol 50: 77–89, 2000.[CrossRef][Web of Science][Medline]
  12. Morita T, Araki J, Oshima Y, Mitani H, Iribe G, Mohri S, Shimizu J, Sano S, Kajiya F, and Suga H. Frequency distribution, variance, and moving average of left ventricular rhythm and contractility during atrial fibrillation in dog. Jpn J Physiol 52: 41–49, 2002.[CrossRef][Web of Science][Medline]
  13. Muntinga HJ, Gosselink AT, Bianksma PK, De Kam PJ, Van Der Wall EE, and Crijns HJ. Left ventricular beat-to-beat performance in atrial fibrillation: dependence on contractility, preload, and afterload. Heart 82: 575–580, 1999.[Abstract/Free Full Text]
  14. Nozawa T, Yasumura Y, Futaki S, Tanaka N, Uenishi M, and Suga H. Efficiency of energy transfer from pressure-volume area to external mechanical work increases with contractile state and decreases with afterload in the left ventricle of the anesthetized closed-chest dog. Circulation 77: 1116–1124, 1988.[Abstract/Free Full Text]
  15. Rawles JM and Rowland E. Is the pulse in atrial fibrillation irregularly irregular? Br Heart J 56: 4–11, 1986.[Abstract/Free Full Text]
  16. Sagawa K, Maughan L, Suga H, and Sunagawa K. Cardiac Contraction and the Pressure-Volume Relationship. New York: Oxford University Press, 1988, p. 232–298.
  17. Snedecor GW and Cochran WG. Statistical Methods (6th ed.). Ames: Iowa State University Press, 1971.
  18. Suga H. Ventricular energetics. Physiol Rev 70: 247–277, 1990.[Free Full Text]
  19. Suga H. Cardiac mechanics, and energetics—from Emax to PVA. Front Med Biol Eng 2: 3–22, 1990.[Medline]
  20. Suga H, Yasumura Y, Nozawa T, Futaki S, and Tanaka N. Pressure-volume relation around zero transmural pressure in excised cross-circulated dog left ventricle. Circ Res 63: 361–372, 1988.[Abstract/Free Full Text]
  21. Sunagawa K, Maughan WL, Burkhoff D, and Sagawa K. Left ventricular interaction with arterial load studied in isolated canine heart. Am J Physiol Heart Circ Physiol 245: H773–H780, 1983.[Abstract/Free Full Text]
  22. Suzuki S, Araki J, Morita T, Mohri S, Mikane T, Yamaguchi H, Sano S, Ohe T, Hirakawa M, and Suga H. Ventricular contractility in atrial fibrillation is predictable by mechanical restitution and potentiation. Am J Physiol Heart Circ Physiol 275: H1513–H1519, 1998.[Abstract/Free Full Text]
  23. Watt JH, Donner AP, McKinney CM, and Klein GJ. Atrial fibrillation: minimal sampling interval to estimate average rate. J Electrocardiol 17: 153–156, 1984.[Web of Science][Medline]
  24. Yamaguchi H, Takaki M, Ito H, Tachibana H, Lee S, and Suga H. Pressure-interval relationship characterizes left ventricular irregular beat contractilities and their mean level during atrial fibrillation. Jpn J Physiol 47: 101–110, 1997.[CrossRef][Web of Science][Medline]




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
288/4/H1740    most recent
00584.2004v1
Right arrow Alert me when this article is cited
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 Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Mohri, S.
Right arrow Articles by Suga, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Mohri, S.
Right arrow Articles by Suga, H.


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