Allometric scaling of the maximum metabolic rate of mammals: oxygen transport from the lungs to the heart is a limiting step
© Painter; licensee BioMed Central Ltd. 2005
Received: 22 March 2005
Accepted: 11 August 2005
Published: 11 August 2005
The maximum metabolic rate (MMR) of mammals is approximately proportional to M 0.9 , where M is the mammal's body weight. Therefore, MMR increases with body weight faster than does the basal metabolic rate (BMR), which is approximately proportional to M 0.7 . MMR is strongly associated with the capacity of the cardiovascular system to deliver blood to capillaries in the systemic circulation, but properties of this vascular system have not produced an explanation for the scaling of MMR.
Here we focus on the pulmonary circulation where resistance to blood flow (impedance) places a limit on the rate that blood can be pumped through the lungs before pulmonary edema occurs. The maximum pressure gradient that does not produce edema determines the maximum rate that blood can flow through the pulmonary veins without compromising the diffusing capacity of oxygen. We show that modeling the pulmonary venous tree as a fractal-like vascular network leads to a scaling equation for maximum cardiac output that predicts MMR as a function of M as well as the conventional power function aM b does and that least-squares regression estimates of the equation's slope-determining parameter correspond closely to the value of the parameter calculated directly from Murray's law.
The assumption that cardiac output at the MMR is limited by pulmonary capillary pressures that produce edema leads to a model that is in agreement with experimental measurements of MMR scaling, and the rate of blood flow in pulmonary veins may be rate-limiting for the pathway of oxygen.
The maximum metabolic rate (MMR) of mammals is measured as the rate of oxygen consumption during the maximum sustainable rate of exercise . Unlike the basal metabolic rate (BMR), which consumes oxygen at rates far below the delivery capacity of the cardiovascular system [1, 2], the MMR is largely determined by the maximal rate that the cardiovascular system can deliver oxygen to mitochondria in muscle tissue .
MMR has been measured in mammals ranging in size, M, from 0.007 kg (pygmy mice) to 575 kg (cattle). Regression of the logarithm of MMR (denoted Q) on the logarithm of M gives a maximum-likelihood estimate (MLE) of the exponent b in the allometric expression
Q = aM b (1)
of 0.872 with a 95% confidence interval (CI) of 0.812–0.931 for MMR data from 32 mammalian species . In contrast, regression analysis of BMR data from 619 mammalian species gives a MLE of the slope, 0.69, with 95% CI, 0.68–0.70 
To explain the scaling of the metabolic rate in mammals, West et al.  and Bengtson and Eden  model the arterial network as a structure that starts with a single tube (aorta) that repeatedly branches into two (or more) smaller tubes. Branching continues until a tube (small arteriole) that supplies capillaries is reached. They assume that all paths from the heart to capillaries pass through n tubes and that the arterial network is a truncated self-similar fractal (i.e., a fractal-like network). The smallest vessels of the circulatory system have dimensions that vary little with body size, whereas the dimensions of the aorta and other great vessels are highly dependent on size. For convenience, we define level 1 of the arterial tree (or venous tree) as the smallest arterioles (or venules). These have radius r 1 and length l 1 . Each level 2 vascular tube with radius r 2 and length l 2 is connected to η 1 level 1 structures. In general, each level i+1 tube of radius ri+1and length li+1is connected to η i level i tubes. It follows from the assumption of a self-similar fractal that the branching ratio is a constant (denoted η) and that the ratio of tube lengths, li+1/l i , is also a constant (denoted L) throughout the network.
The theory of West et al. minimizes the (pressure) × (volume) work of the heart that is required to pump one unit of blood against a difference in pressure equal to the pressure in the aorta minus the pressure in capillaries. This work per unit of blood flow is proportional to the impedance in the arterial network. Minimization of this energy cost for pulsatile flow in arteries is claimed to require area-preserving branching of the network (i.e., the ratio ri+1/r i , termed R, is equal to η1/2) and, as a consequence, to require that the density of capillaries in tissues is proportional to M-1/4(assuming that the diameter of the aorta scales proportionally to M3/8or that arterial blood volume scales proportionally to M). The theory's 3/4-power scaling prediction for metabolic rate follows from the assumption that metabolic rate is proportional to the total number of capillaries calculated as tissue capillary density multiplied by M, an assumption that is reasonable for MMR but not for BMR . The theory of Bengtson and Eden assumes that energy dissipation per endothelial surface area is constant, leading to the conclusions that R is equal to η2/5and that the total number of capillaries is proportional to M15/17if the volume of blood in arteries scales proportionally to M. If it is assumed that the diameter of the aorta scales proportionally to M3/8, the number of capillaries is proportional to M15/16.
The scaling of the total number of capillaries in skeletal muscle, where over 90% of energy metabolism occurs during MMR exercise, is nearly identical to the scaling of MMR , and, as noted above, this scaling is not proportional to M3/4. The 95% CI for the scaling exponent for total capillary volume, 0.909 – 1.0559, contains 15/16 but not 3/4. Moreover, if either of these theories is adequate for predicting capillary density, it should correctly predict the scaling exponent for capillaries in the lung, which is 1.00 with 95% CI of 0.912 – 1.087 . This CI contains 15/16 but not 3/4.
A model for the maximum metabolic rate
While minimization of impedance does not by itself lead to a correct prediction of capillary density in muscle and lung tissue, it is clearly an important principle for design of mammalian vascular systems [7, 8]. The potential importance of impedance is most apparent in the pulmonary venous circulation, where the entire output of the heart's right ventricle flows before blood enters the left atrium of the heart. The driving force for pulmonary venous return to the heart is the pressure at the venous end of pulmonary capillaries minus the diastolic pressure in the left atrium (denoted P LA ).
The output of oxygen by the left ventricle of the heart into the aorta is equal to the input of oxygen from the lungs to the heart. This is equal to the cardiac blood output rate multiplied by the maximum amount of oxygen per ml of blood multiplied by the percent saturation of blood with oxygen. Pressure in the model is strictly increasing with flow. However, as pressure rises above oncotic pressure, interstitial edema increases and then more and more fluid accumulates within alveoli. Therefore, oxygen saturation is strictly decreasing as a consequence of the increasing barrier to oxygen diffusion from pulmonary air into capillaries. As a result, there is a blood flow rate, denoted F max , that produces the maximum uptake of oxygen in the lungs, which is also the maximum output of oxygen to the body. The pressure near the venous end of alveolar capillaries at F max is denoted Π max . Consequently, the pressure gradient that drives the return of blood in pulmonary capillaries back to the heart is
ΔP max = F max I p (2)
where ΔP max = Π max - P LA and I p is the impedance of the pulmonary venous network. It is assumed that Π max is proportional to the oncotic pressure of blood, denoted Π o . The value of Π max is assumed to be approximately the same in mammals of different sizes because Π o appears to be nearly invariant in mammalian species, being approximately 20 mm Hg [9–11] and P LA is approximately 1 mm Hg. (All pressures in this article are measured relative to ambient pressure.) Therefore, the scaling of F max with body size depends largely on the scaling of I p .
The impedance of the pulmonary venous network is a consequence of its physical structure and the viscosity of blood (termed ν). The pulmonary arteries and veins form parallel fractal-like networks in each lung with arteries and veins of the same level having similar dimensions [12, 13]. Small venules have dimensions that are body-size-invariant (r 1 approximately 10-5 m and l 1 approximately 10-4 m). These vascular tubes receive blood from the capillaries in pulmonary acini, the structures that comprise approximately 10,000 alveoli and that appear to be body-size-invariant in mammals .
The impedance of a fractal-like network is the sum of impedances contributed by each level of the network. We assume that the impedance I i due to level i is the value calculated from the Poiseuille theory for non-turbulent fluid flow, , where N i is the number of level i vessels . Consequently, Ii+1is equal to . The observation that dimensions within acini are size-invariant leads to the conclusion that η, R and L are size-invariant in acini. We assume that these ratios remain constant throughout the network. Therefore, the factor η L/R 4 (denoted α) is assumed to be size-invariant, and the expression for I p is a geometric series (when α ≠ 1) that simplifies to
Substitution of this formula into Equation (2) gives
where C is a constant. Note that ζ depends on the base used to define the logarithm. The base 10 is used in the following regression analysis. When α = 1, we have
Testing model predictions
Maximum metabolic rates (VO2 max) of mammals from Weibel et al..
VO2 max (ml/min)
Maximum metabolic rates of mammals adjusted to a standard hematocrit of 0.42 from Weibel et al..
Body mass (kg)
VO2 max (ml/min)
LSLR using the data in Table 2 and the model of Equation (1) gives the value of 0.957 for b (Rc2 = .9697) and SSR = 0.5890) when the SSR is minimized. LSLR using Equation (5) finds that the SSR is minimized when ζ equals 0.801 (SSR = 0.5833). LSLR of predicted values of cardiac output from Equation (5) using values of M from Table 2 and the estimate for ζ of 0.801 gives b = 0.958 and Rc2 = 0.9991. Clearly the predictions from Equation (5) are again nearly indistinguishable from those of Equation (1), and Equation (5) fits these data as well as Equation (1) does.
Comparison with Murray's law
L = η 1/3 , (7)
and this equation for L leads to the formula
Equation (8) is remarkably similar to Murray's law for the scaling of radii of arterial or venous networks, which states that flow rate is proportional to the third power of vessel radius . For our network model, Murray's law implies R 3 = η, and this equation together with the condition L = η 1/3 implies α = 1. With this value of α, the slope of the logarithm of Equation (6) depends only on the estimate of M 1 . For M 1 = 0.00001 kg, Equation (6) is nearly identical to a power function with b = 0.916. Therefore, Murray's law and the fractal length scaling relationship lead to the constrained PVFCP model and predict that the slope parameter of the scaling function is in the range of observed values.
The PVFCP model predicts that the logarithm of maximum oxygen uptake in mammals is approximately proportional to the logarithm of body mass. If the radii of veins in the pulmonary venous tree obey Murray's law, then the constant of proportionality is in the range of experimentally observed values for MMR. The PVFCP model, like other published explanations for MMR scaling, focuses on the supply of oxygen to the tissues. However, the PVFCP model differs from other explanations for MMR scaling because it focuses on pulmonary blood flow.
The PVFCP model and the model of Bengtson and Eden  use the same mathematical description of pressure-flow relationships in a vascular tree. While the model of Bengtson and Eden  is consistent with current data on MMR, the model's assumption of energy dissipation that is proportional to vascular surface area is questionable as a principle of mammalian design. For example, a hypothetical mammalian species that replaces the R = η 2/5 requirement of their theory with the R = η 1/3 relationship of Murray's law would reduce total energy dissipation in arteries. This replacement would also give a higher predicted capillary density and consequently a higher MMR.
It is instructive to compare the number of independent parameters and assumptions in the PVFCP model with the number of parameters and assumptions in the two fractal-like models of the arterial network that predict metabolic scaling [4, 5]. All three models describe the vascular network as a self-similar fractal-like tubular structure with pressure gradients determined by Poiseuille's law. All assume that the size of terminal (smallest) network tubes is the same in mammals of different size and that blood viscosity does not vary with body size. All contain the branching ratio parameter η and the network length parameter n. In the PVFCP model, a relationship between η, n and body mass is derived from the assumption that the number of terminal segments is proportional to body mass, an assumption that is supported by observations. In the other two models, a relation between these parameters is derived from the assumption that arterial blood volume is proportional to body mass, an assumption without direct observational support. Network structure is related to metabolic rate in the PVFCP model by Equation (5), which specifies the maximum rate of blood flow that does not compromise pulmonary function. In the other models, such a relation is derived from the assumption that metabolic rate is proportional to the number of capillaries in the systemic circulation. In the PVFCP model, there is one more independent parameter, α, which is defined by fitting experimental data. The other models have two additional parameters, L and R. Both models specify L indirectly using the assumption of Equation (7). The parameter R is specified by an energy minimization principle in one model  and by an energy dissipation principle in the other . While the number of parameters and assumptions in the PVFCP model is relatively large, it is less than the number in the fractal-like network models previously published. Another recent mathematical description of metabolic scaling, the "Allometric Cascade" model , is not discussed here because it is not a mechanistic model. Indeed, the two models appear compatible because the PVFCP model could be integrated into the "Allometric Cascade" model to provide a mechanism-based scaling term for the maximum rate of blood flow.
Weibel et al.  argue that it is the volume of mitochondria in muscle tissue and the blood supply in capillaries in muscle tissue that determine the scaling of MMR. This view is supported by their demonstration that MMR is remarkably correlated with and is proportional to mitochondrial volume (b = 1.09, R c 2 = 0.9939) and to estimated capillary blood volume in muscle tissue (b = 0.975, R c 2 = 0.9846). However, total mitochondrial volume and blood volume in muscle capillaries can be increased by exercise conditioning, and the correlation between capillary surface area and MMR or between mitochondrial volume and MMR may arise from such conditioning.
In the formulation of the PVFCP model, the role of gravity in facilitating or impeding the return of pulmonary blood to the heart has been ignored. Blood that is one inch higher than the left atrium has potential energy to facilitate its return to the heart that is approximately equivalent to a 2 mm Hg pressure gradient. For small mammals (e.g., mice), gravitational effects would be small compared with the approximately 20 mm Hg pressure gradient that we assume drives blood return during MMR exercise. However, for large mammals (e.g., elephants and whales), the effects of gravity will significantly increase blood return from regions of lung above the heart, but decrease blood return from regions below the heart. Therefore, Equation (5) may not adequately describe MMR blood flow in large mammals.
A second reason for doubting the validity of Equation (5) for large mammals is that intervals of the heart cycle increase with body size. The minimum length of the heart cycle (at maximum heart rate) is largely composed of the time required for the ventricles to fill plus the time required for the ventricles to eject blood into the pulmonary artery and aorta. At maximal heart rate, ventricular filling time is nearly equal to the PR interval, which is approximately proportional to the 1/4-power of body mass . If the sum of the QRS interval and the ST segment, which is nearly equal to the time required to eject blood from the ventricles, has similar scaling, then the scaling exponent for maximum heart rate is less than the scaling exponent for the MMR divided by body mass, i.e., the specific maximum metabolic rate (SMMR). Thus, maximum heart rate, not the limitation posed by pulmonary venous impedance, may limit MMR for very large mammals.
The biological plausibility of the relation between MMR and I p proposed in the PVFCP model depends on whether pressures in lung capillaries approach the oncotic pressure of blood during periods of maximal exertion. In healthy humans at rest, the pressure difference between pulmonary capillaries and the left atrium ranges from approximately 5 to 11 mm Hg . Assuming that the value of 5 mm Hg occurs when pulmonary veins are dilated, this pressure difference is predicted to increase by a factor of approximately 4 during heavy exercise in a trained athlete when cardiac output increases by a factor of 4 (assuming that the pulmonary veins are in a comparable state of dilation). This would require the capillary pressure to rise to approximately 21 mm Hg. It is noteworthy that signs of pressure stress are sometimes observed in pulmonary tissue from trained endurance athletes .
Studies of human patients with narrowing of the mitral valve, the valve between the left atrium and left ventricle, are consistent with the hypothesis that I p limits maximum metabolic rate. This condition, termed mitral stenosis, causes an increase in P LA . Patients with a P LA below 20 mm Hg usually do not have pulmonary edema at rest but may develop it with exercise. Furthermore, women with a P LA between 18 and 20 mm Hg are at risk for developing pulmonary edema during pregnancy where the cardiac output at rest increases on average by approximately 50% [20–22].
Additional support for the proposed role of pulmonary impedance in determining MMR comes from studies of horses, which have an MMR well above the value predicted by the allometric equation fitted to the data in Table 1. Horses at rest have pulmonary capillary blood pressures that are above those in humans with mitral stenosis and pulmonary edema with exercise. Horses are apparently able to exercise without developing pulmonary edema because they are able to "concentrate" their blood during periods of exertion. The concentration of erythrocytes (measured as the hematocrit) is increased during exercise . This requires a preferential loss of water that likely occurs in capillaries of the systemic circulation. As a result, the concentration of albumin in blood is increased and the oncotic pressure of blood is increased. This adaptation enables a horse at a gallop to tolerate pulmonary capillary pressures as high as 38 mm Hg .
Horses possess a second adaptation that allows them to increase their SMMR. Their ratio of lung volume to body mass is approximately 20% greater than the average value for mammals . To pump blood through their large lungs at an unusually high rate per unit lung volume, horses possess a heart that is larger (as a fraction of body mass) than the average value for mammals . This enables them to achieve a SMMR that is more than twice that of a cow of similar size. However, even with its remarkable adaptations, no horse can sustain the SMMR that pygmy mice and other small mammals can achieve .
I thank Charles Salocks and Danielle Ketchum for their careful reviews and helpful comments.
- Weibel ER, Bacigalupe LD, Schmitt B, Hoppeler H: Allometric scaling of maximal metabolic rate in mammals: muscle aerobic capacity as determinant factor. Respir Physiol Neurobiol. 2004, 140: 115-132. 10.1016/j.resp.2004.01.006.View ArticlePubMedGoogle Scholar
- Darveau C-A, Suarez RK, Andrews RD, Hochachka PW: Allometric cascade as a unifying principle of body mass effects on metabolism. Nature. 2002, 417: 166-170. 10.1038/417166a.View ArticlePubMedGoogle Scholar
- White CR, Seymour RS: Mammalian basal metabolic rate is proportional to body mass2/3. Proc Natl Acad Sci USA. 2003, 100: 4046-4049. 10.1073/pnas.0436428100.PubMed CentralView ArticlePubMedGoogle Scholar
- West GB, Brown JH, Enquist BJ: A general model for the origin of allometric scaling laws in biology. Science. 1997, 276: 122-126. 10.1126/science.276.5309.122.View ArticlePubMedGoogle Scholar
- Bengtson H-U, Eden P: A simple model for the arterial system. J Theoret Biol. 2003, 221: 437-443. 10.1006/jtbi.2003.3198.View ArticleGoogle Scholar
- Gehr P, Mwangi DK, Ammann A, Maloiy GM, Taylor CR, Weibel ER: Design of the mammalian respiratoey system. V. Scaling morphometric pulmonary diffusing capacity to body mass: wild and domestic animals. Respir Physiol. 1981, 44: 61-86. 10.1016/0034-5687(81)90077-3.View ArticlePubMedGoogle Scholar
- Murray CD: The physiological principle of minimum work. I. The vascular system and the cost of blood volume. Proc Natl Acad Sci USA. 1926, 12: 207-214.PubMed CentralView ArticlePubMedGoogle Scholar
- LaBarbera M: Principles of design of fluid transport systems in zoology. Science. 1990, 249: 992-1000.View ArticlePubMedGoogle Scholar
- Stohrer M, Boucher Y, Stangassinger M, Jain RK: Oncotic pressure in solid tumors is elevated. Cancer Res. 2000, 60: 4251-4255.PubMedGoogle Scholar
- Gabel JC, Scott RL, Adair TH, Drake RE, Traber DL: Errors in calculated oncotic pressure of dog plasma. Am J Physiol. 1980, 239: H810-H812.PubMedGoogle Scholar
- Madigan JE, Rahal CJ: Measurement of Plasma Colloid Osmotic Pressure in Normal Thoroughbred Neonatal Foals. 2003, Davis California: University of CaliforniaGoogle Scholar
- Agur AM, Ming JL, Grant JC: Grant's Atlas of Anatomy. 2004, Baltimore: Lippincott, Williams and WilliamsGoogle Scholar
- Bergman RA, Afifi AK, Heidger PM: Atlas of Microscopic Anatomy. 1996, Philadelphia: WB SaundersGoogle Scholar
- Weibel ER: Fractal geometry: a design principle for living organisms. Am J Physiol. 1991, 261: L361-L369.PubMedGoogle Scholar
- Turcotte DL, Pelletier JD, Newman WI: Networks with side branching in biology. J Theoret Biol. 1998, 193: 577-592. 10.1006/jtbi.1998.0723.View ArticleGoogle Scholar
- Fung YC: Biomechanics. 1990, New York: Springer-VerlagView ArticleGoogle Scholar
- Noujaim SF, Lucca E, Munoz 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: 2802-2808. 10.1161/01.CIR.0000146785.15995.67.View ArticlePubMedGoogle Scholar
- Tabata T, Oki T, Fukuda N, Iuchi A, Manabe K, Kageji Y, Sasaki M, Yamada H, Ito S: Influence of left atrial pressure on left atrial appendage flow velocity patterns in patients in sinus rhythm. J Am Soc Echocardiogr. 1996, 9: 857-864. 10.1016/S0894-7317(96)90478-2.View ArticlePubMedGoogle Scholar
- West JB: Vulnerability of pulmonary capillaries during exercise. Exerc Sport Sci Rev. 2004, 32: 24-30. 10.1097/00003677-200401000-00006.View ArticlePubMedGoogle Scholar
- Reis G, Motta MS, Barbosa MM, Esteves WA, Souza SF, Bocchi EA: Dobutamine stress echocardiography for noninvasive assessment and risk stratification of patients with rheumatic mitral stenosis. J Am Col Cardiol. 2004, 43: 393-401. 10.1016/j.jacc.2003.09.037.View ArticleGoogle Scholar
- Desai DK, Adanlawo M, Naidoo DP, Moodley J, Kleinschmidt I: Mitral stenosis in pregnancy: a four-year experience at King Edward VIII Hospital, Durban, South Africa. BJOG. 2000, 107: 953-958.View ArticlePubMedGoogle Scholar
- van Oppen AC, Stigter RH, Bruinse HW: Cardiac output in normal pregnancy: a critical review. Obstet Gynecol. 1996, 87: 310-318. 10.1016/0029-7844(95)00348-7.View ArticlePubMedGoogle Scholar
- Weber J-M, Dobson GP, Parkhouse WS, Wheeldon D, Harman JC, Snow DH, Hochachka PW: Cardiac output and oxygen consumption in exercising Thoroughbred horses. Am J Physiol. 1987, 253: R890-R895.PubMedGoogle Scholar
- Hackett RP, Ducharme NG, Gleed RD, Mitchell L, Soderholm LV, Erickson BK, Erb HN: Do Thoroughbred and Standardbred horses have similar increases in pulmonary vascular pressures during exercise?. Can J Vet Res. 2003, 67: 291-296.PubMed CentralPubMedGoogle Scholar
- Hoppeler H, Lindstedt SL, Claassen H, Taylor CR, Mathieu O, Weibel ER: Scaling mitochondrial volume in heart to body mass. Respir Physiol. 1984, 55: 131-137. 10.1016/0034-5687(84)90018-5.View ArticlePubMedGoogle Scholar
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.