# Probability of ventricular fibrillation: allometric model based on the ST deviation

- Maria P Bonomini
^{1}, - Pedro D Arini
^{1, 2}and - Max E Valentinuzzi
^{1}Email author

**10**:2

https://doi.org/10.1186/1475-925X-10-2

© Bonomini et al; licensee BioMed Central Ltd. 2011

**Received: **23 November 2010

**Accepted: **13 January 2011

**Published: **13 January 2011

## Abstract

### Background

Allometry, in general biology, measures the relative growth of a part in relation to the whole living organism. Using reported clinical data, we apply this concept for evaluating the probability of ventricular fibrillation based on the electrocardiographic ST-segment deviation values.

### Methods

Data collected by previous reports were used to fit an allometric model in order to estimate ventricular fibrillation probability. Patients presenting either with death, myocardial infarction or unstable angina were included to calculate such probability as, *VF*
_{
p
}
*= δ + β (ST)*, for three different ST deviations. The coefficients *δ* and *β* were obtained as the best fit to the clinical data extended over observational periods of 1, 6, 12 and 48 months from occurrence of the first reported chest pain accompanied by ST deviation.

### Results

By application of the above equation in log-log representation, the fitting procedure produced the following overall coefficients: Average *β* = 0.46, with a maximum = 0.62 and a minimum = 0.42; Average *δ* = 1.28, with a maximum = 1.79 and a minimum = 0.92. For a 2 mm ST-deviation, the full range of predicted ventricular fibrillation probability extended from about 13% at 1 month up to 86% at 4 years after the original cardiac event.

### Conclusions

These results, at least preliminarily, appear acceptable and still call for full clinical test. The model seems promising, especially if other parameters were taken into account, such as blood cardiac enzyme concentrations, ischemic or infarcted epicardial areas or ejection fraction. It is concluded, considering these results and a few references found in the literature, that the allometric model shows good predictive practical value to aid medical decisions.

## Background

Ventricular fibrillation can be viewed as a probabilistic event that appears biased under certain pathophysiological and daily life situations. Physicians, in their practice, try to predict as close as possible how high such probability is. Since most of cardiac deaths are due to ventricular fibrillation, it can be said that such death, in the end, would be an event that counts as a ventricular fibrillation. Valentinuzzi, in 2010, has reviewed at large such arrhythmia and its countermeasures [1]. Empirical tests, as possible quantitative criteria to screen out patients of high risk (that is, searching for a better answer to the question *shall we confine the patient to the coronary unit*?) have been attempted with moderate success, but always the degree of uncertainty is rather large. In such endeavor, we might try an appealing and old universal scaling, the allometric law, although in principle apparently not related to the cardiac risk overall concept, it might find a place in it and at least deserves to be reminded bringing about first a nice and well carried out paper by Noujaim *et al*, in 2004 [2]. In it, it is recalled that from mouse to whale the electrocardiographic PR interval increases 10^{1} times whereas body mass (*B*
_{
M
}) augments 10^{6}. This is the first use we found of the allometric law in cardiology encouraging us to proceed further up.

Scaling of many biological processes can be described by the allometric equation, *Y* = *a*(*B*
_{
M
})^{
b
}, where *Y* is the biological process and *a* and *b* are scaling constants. In general, the weights of most individual organs scale as a constant fraction of body mass (i.e., the body mass exponent, *b* equals 1.0). Biological rates (heart rate, respiratory rate) scale as *b* close to 0.25. Finally, volume rates (the product of volume and rate) such as cardiac output, ventilation and oxygen uptake vary as *b* around 0.75. These emergent patterns provide insights into body-size dependent 'principles of design' that seem to dictate several blueprint aspects and function across species among all mammals [3].

Noujaim *et al* [2] assumed that the heart behaves as a set of "fractal-like" networks tending to minimize propagation time across the conducting system while ensuring a hemodynamically optimal atrioventricular activation sequence. With the potential relationship given above and, subsequently, based on previously published values of PR interval, heart rate, and body masses of 541 mammals, they reported as best fit the equation *PR* = 53(*B*
_{
M
})^{0.24}.

Inspired in the latter report, the following question seems pertinent: Would a relationship similar to the allometric equation be conceivable, say, between the probability of cardiac risk (or cardiac event or episode, all equivalent terms) and heart weight, or perhaps other parameter somehow related to the latter, as for example, the number of cardiac diseased fibers or the ST shift seen in the ECG? The objective of this communication tries to find an answer to such question. The mentioned ECG deviation appears as a good candidate because well-known is the fact that the larger the ST step, in either direction (upward or downward), the larger the compromised myocardial mass. In some cases, such change includes inverted T-wave, which tends to complicate the wave-pattern. Many reports confirm this concept, such as Klootwijk, in 1998 [4], Kléber, in 2000 [5], or Balian *et al*, in 2006 [6], among others, where often the ST shift is defined as a change of ST amplitude in one or more leads of at least ± 100 μV from the baseline ST level, developing within a 10 minute period and persisting for at least 1 minute. However, differences among authors regarding these criteria are frequent.

## Methods

### Theoretical background

*y*is related to the growth velocity of another component (or the whole organism)

*x*in a constant way. This was clearly described by von Bertalanffy in 1957 [8]. Thus, the relative rate of change of a given event

*y*is proportional to the relative rate of change of body mass or body weight

*x*, i.e.,

*y*was the weight of an organ (heart, stomach, other) and

*x*was body weight or mass. The parameters

*A*and

*B*require numerical estimation by an appropriate procedure usually using empirical information. By the same token, let us say that the probability of fibrillation

*P*

_{ F }(and we use

*P*

_{ F }because, as stated above, most of the cardiac episodes end up in ventricular fibrillation) follows a relationship with the number of ventricular diseased fibers (

*N*

_{ DF }) formally equal to (2), i.e.,

*y*in equation (3) is replaced by

*P*

_{ F }in (4), and

*N*

_{ DF }in the latter takes the place of

*x*in the former. After all, the number of diseased cardiac fibers (ischemic or infarcted or both) are part of the cardiac mass. Besides, since the electrocardiographic ST-segment deviation (Δ

_{ ST }) is a traditional estimator of cardiac injury, it sounds sensible to state that,

*VF*

_{ P }as ventricular fibrillation probability, where

Hence, equation (8) in log-log plot would represent the probability of fibrillation as function of the ECG ST-depression or elevation.

### Numerical procedure

To calculate out the two constants *δ* and *β* of equation (8) and later on apply the mathematical expression for predictive purposes, the probability of the data having occurred can be estimated by, (a) simply assuming an arbitrary and theoretical set of coupled pairs of numbers, as for example, a quadratic law of the type w = K z^{2}, that is, the *VF*
_{
P
}would be accepted as being proportional to the square of the ST, or (b) using a particular hypothesis, say, based on clinical data.

*et al*[9] and another from Kaul

*et al*[10], as follows:

- (i)
In the first one, 642 patients had been admitted to coronary care unit with prolonged chest pain. Due to the exclusion criteria applied by these authors, 469 were removed leaving a net number of 173 for their study. Besides, they reported survival rates at 1 and 4 years after the first admission.

- (ii)
In the second paper (PARAGON-A trial), out of 2,282 patients with chest discomfort within the previous 12 hours, there was a screen out of 694 due to either missing or not clear enough records leaving a net of 1,588 cases. They were evaluated at 1 month and 6 months.

- (iii)
Besides, the latter authors had 8,001 patients (GUSTO-IIb trial) comparing hirudin and heparin therapy when unstable angina or acute myocardial infarction was present without ST-segment. Out of this total, only 6,301 were evaluated at 1 month and 6 months.

In Hyde *et al* [9], patients with ≥ 0.5 mm ST-segment depression were classified as "true depression". This deviation was subclassified as 0.5 mm, 1 mm or ≥ 2 mm. In their own words, "ST segment depression was measured using calipers 80 ms after the J point in intervals of 0.5 mm. ECGs were analyzed blinded to the clinical outcome".

The ST-segment criteria in Kaul *et al* [10], instead, rounded out the depression of 0.5 mm to 1 mm, of 1.5 mm to 2 mm, including in the latter larger deviations, distinguishing three groups: No ST-segment depression, 1 mm ST-segment depression in two contiguous leads, and ST-segment depression of 2 mm in two contiguous leads. The 12-lead ECGs were recorded at a paper speed of 25 mm/s. ST segment depression was judged to be present if the J point was depressed by ≥1 mm and was followed by a horizontal or downsloping ST segment for at least 0.08s in one or more of the 12 leads, except for the aVR lead.

*β*and

*δ*and goodness of fit

*r*

^{ 2 }were computed by linear regression and all quadratic fits used values within the 0.5-2 mm range (see Table 1). It should be recalled that the standard ECG calibration of 10 mm = 1 mV is used in all the paper.

Coefficients *β* and *δ*, in equation (6), and adjusted r-square for all 9 curves

Months | β | δ | Adjusted r-square | |
---|---|---|---|---|

| 1 | 0.42 | 1.09 | 0.9884 |

| 6 | 0.45 | 1.23 | 0.9917 |

| 1 | 0.62 | 0.92 | 0.9245 |

| 6 | 0.54 | 1.09 | 0.9427 |

| 1 | 0.50 | 1.01 | 0.9615 |

| 6 | 0.49 | 1.16 | 0.9734 |

12 | 0.54 | 1.56 | 0.9844 | |

| 48 | 0.48 | 1.79 | 0.9477 |

| 1 to 48 | 0.46 | 1.28 | 0.9959 |

## Results

*et al*[10]. As example, for 1.5 mm shift, the predicted probability at one month ranges from slightly below 11% (GUSTO-IIb) to 15% (PARAGON-A). For the same ST shift selected above, the probability values at 6 months span from 15% (GUSTO-IIb) to a 20% (PARAGON-A). Finally, the upper two curves describe the behavior at 1 (open squares) and 4 years (dark squares) after the event, according to Hyde

*et al*[9]. For the same previous ST deviation, the foreseen probability range goes from 45% to about 75%. The fitted adjustments pass essentially through the depicted points.

*et al*[10], respectively, at 1 and 6 months, from bottom to top, showing also one half of the standard deviation for each data point. Notice the spread decrease comparing the bottom with the upper curve at any ST value. Besides, an inverse relationship between ST deviation and standard deviation is manifest, which speaks of the gradual nature of ST changes and, therefore, points out to the importance of such amplitude.

## Discussion

This report has developed an allometric equation simply based on the electrocardiographic ST-segment deviation. The calculated coefficients permit predictions at different times after the first cardiac episode or, with a much wider spread, as an overall quantitative evaluation applying the relationship given in Figure 3. Obviously, the model must be tested in the clinical environment to better assess its accuracy and predictive power. In Figure 2, dispersion increases at 1 month after the cardiac episode, which might be interpreted as the patient still traversing a period of dangerous instability. Conversely, the upper curve, after a longer time, shows a marked spread decrease. We read this fact as a stable condition because of compensation.

With the aging process, along perhaps with an increase in the ischemic areas or deterioration of the myocardial scar tissue, it seems quite acceptable a consequent increase, too, in the probability of an arrhythmic event, as well depicted in the three figures, throughout longer observational periods.

It is convenient to underline that the best fit quadratic equation (see Methods) supplies the numerical information needed to estimate the *β* and *δ* constants of the allometric law.

This model uses only the ST-segment as criterion, which obviously leaves out other possible parameters, such as myocardial enzymes (CPK, for example), quantitatively obtainable by blood sample analysis, or ischemic or infarcted epicardial surface, from appropriate imaging procedures, or ejection fraction as evaluated by echocardiography. Other anthropological data, such as patient's sex and age, could also be included. Any of these criteria would lead to allometric equations as the one herein reported. One tempting and difficult approach would try to combine all the mentioned parameters in a single mathematical model.

The results herein presented foresee a direct application in the clinical environment to better predict the evaluation of a cardiac patient. However, this kind of validation remains to be carried out.

## Conclusions

The allometric statement seems to maintain interest, especially in general mammalian biology [11, 12] and the results reported here would indicate an attractive line of research with their consequent clinical tests. Once more, it should be underlined the proportion basis of the allometric statement, since it numerically links here a specific number of compromised fibers (ischemic or even dead) with the concept of cardiac risk.

## Declarations

### Acknowledgements

Supported by grants PIP #538 from *Consejo Nacional de Investigaciones Científicas y Técnicas* (CONICET) and *Fundación Florencio Fiori* ni, both research institutions from Argentina. Details of the numerical calculations can be requested via email from the authors by any interested reader. This paper was partially communicated to the IEEE/32nd International Conference on Engineering in Medicine and Biology, September 1-4, Buenos Aires, Argentina.

## Authors’ Affiliations

## References

- Valentinuzzi ME:
**Cardiac Fibrillation-Defibrillation: Clinical and Engineering Aspects.**In*Series on Bioengineering and Biomedical Engineering*. World Scientific Publishers, Singapore; 2010.Google Scholar - Noujaim SF, Lucca E, Muñoz V, Persaud D, Berenfeld O, Meijler FL, Jalife J:
**From mouse to whale: A universal scaling relation for the PR interval of the electrocardiogram of mammals.***Circulation*2004,**110:**2801–2808. 10.1161/01.CIR.0000146785.15995.67View ArticleGoogle Scholar - Lindstedt SL, Schaeffer PJ:
**Use of allometry in predicting anatomical and physiological parameters of mammals.***Laboratory Animals*2002,**36:**1–19. 10.1258/0023677021911731View ArticleGoogle Scholar - Klootwijk APJ:
**Dynamic computer-assisted ST segment monitoring in patients with acute coronary syndromes.**In*PhD Thesis*. Erasmus Universiteit, Rótterdam; 1998.Google Scholar - Kléber AG:
**ST-segment elevation in the electrocardiogram: a sign of myocardial ischemia.***Cardiovasc Res*2000,**45**(1):111–118.View ArticleGoogle Scholar - Balian V, Galli M, Marcassa C, Cecchin G, Child M, Barlocco F, Petrucci E, Filippini G, Michi R, Onofri M:
**Intracoronary ST-segment shift soon after elective percutaneous coronary intervention accurately predicts periprocedural myocardial injury.***Circulation*2006,**114:**1948–1954. 10.1161/CIRCULATIONAHA.106.620476View ArticleGoogle Scholar - Snell O:
**Die Abhängigkeit des Hirngewichtes von dem Körpergewicht und den geistigen Fähigkeiten (In German, Dependence of brain weight on body weight and the intellectual capacity).***Archiv für Psychiatrie und Nervenkrankheiten*1891,**110:**2801–2808.Google Scholar - Von Bertalanffy L:
**Quantitative laws in metabolism and growth.***Quaterly Review Biology*1957,**32**(3):217–231. 10.1086/401873View ArticleGoogle Scholar - Hyde TA, French JK, Wong CK, Strasnicky IT, Whitlock RML, White HD:
**Four-year survival of patients with acute coronary syndromes without ST-segment elevation and prognostic significance of 0.5 mm ST-segment depression.***American J Cardiol*1999,**84:**379–385. 10.1016/S0002-9149(99)00319-7View ArticleGoogle Scholar - Kaul P, Fu Y, Chang WC, Harrington RA, Wagner GS, Goodman SG, Granger CB, Moliterno DJ, Van der Werf F, Califf RM, Topol EJ, Armstrong PW:
**Prognostic value of ST-segment depression in acute coronary syndromes: Insights from PARAGON-A applied to GUSTO-IIb.***J American College Cardiol*2001,**38**(1):64–71. 10.1016/S0735-1097(01)01307-9View ArticleGoogle Scholar - Adler CP, Costabel U:
**Cell number in human heart in atrophy, hypertrophy, and under the influence of cytostatics.***Recent Adv and Studies in Cardiac Structure and Metabolism*1975,**6:**343–355.Google Scholar - Martin RD, Genoud M, Hemelrijk CK:
**Problems of allometric scaling analysis: Examples from mammalian reproductive biology.***J Experimental Biology*2005,**208:**1731–1747. 10.1242/jeb.01566View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.