|
|
||||||||
1 Program in Biomedical Engineering, 2 Department of Physiology, University of Arizona, Tucson, Arizona 85724
| |
ABSTRACT |
|---|
|
|
|---|
Oxygen transport from capillaries to exercising skeletal muscle is studied by use of a Krogh-type cylinder model. The goal is to predict oxygen consumption under conditions of high demand, on the basis of a consideration of transport processes occurring at the microvascular level. Effects of the decline in oxygen content of blood flowing along capillaries, intravascular resistance to oxygen diffusion, and myoglobin-facilitated diffusion are included. Parameter values are based on human skeletal muscle. The dependence of oxygen consumption on oxygen demand, perfusion, and capillary density are examined. When demand is moderate, the tissue is well oxygenated and consumption is slightly less than demand. When demand is high, capillary oxygen content declines rapidly with axial distance and radial oxygen transport is limited by diffusion resistance within the capillary and the tissue. Under these conditions, much of the tissue is hypoxic, consumption is substantially less than demand, and consumption is strongly dependent on capillary density. Predicted consumption rates are comparable with experimentally observed maximal rates of oxygen consumption.
oxygen diffusion; Krogh cylinder model; maximal oxygen consumption
| |
INTRODUCTION |
|---|
|
|
|---|
OXYGEN TRANSPORT FROM
BLOOD into tissue occurs by passive diffusion. The maximum
distance that oxygen can diffuse from a blood microvessel into
surrounding tissue decreases with increasing oxygen consumption rate
and is a few tens of micrometers in tissues with high oxygen demand,
such as heavily working skeletal muscle. For a tissue with a
given capillary density, the diffusion of oxygen to the mitochondria
where it is consumed is one of the factors limiting the maximal rate of
oxygen consumption (
O2 max). Observations of maximally working single muscles in humans have shown a
substantial oxygen saturation in the venous blood, in the range
15-30% (1, 28, 24, 25), showing that oxygen extraction is incomplete, reflecting diffusive limitation of oxygen transport (33). More complete extraction of the oxygen
available in the blood would result in very low values for the partial
pressure of oxygen (PO2) at the venous end of
capillaries, which would not provide a sufficient
PO2 gradient for oxygen diffusion. Clearly, both convective and diffusive limitations of oxygen transport are
important in determining
O2 max.
A further factor that may limit the oxygen consumption of strongly
stimulated muscle is the maximal rate of turnover of mitochondrial enzymes in the Krebs cycle, which limits oxygen demand. Blomstrand et
al. (3) showed a correlation between the oxygen
consumption rate of a maximally working human leg muscle and the
maximal activity of the enzyme oxoglutarate dehydrogenase. Thus limited
oxygen demand, perfusion limitation, and diffusion limitation may all play a part in determining
O2 max in
skeletal muscle.
At the microvascular level, the distribution of
PO2 in tissue surrounding microvessels in
heavily working muscle has been debated. In some analyses of oxygen
diffusion from blood to tissue (see below), it has been assumed that
the decline in PO2 occurs mainly in the tissue.
However, studies of dog gracilis muscle at
O2 max (11) showed low
values and small gradients of PO2 throughout
most of the tissue. The results were attributed to effects of the
particulate nature of blood on oxygen transport within capillaries and
to the facilitation of oxygen transport by myoglobin. Richardson et al.
(26) used magnetic resonance spectroscopy to measure
average myoglobin saturation in exercising human skeletal muscle and
concluded that average tissue PO2 is low during
maximal exercise. These studies suggested that
PO2 declines steeply in the radial direction
inside capillaries and within the first few micrometers outside
capillaries, with small gradients in the bulk of the tissue.
Since the classical work of Krogh (19), many investigations of oxygen transport to skeletal muscle have used theoretical models (22). Krogh's model is based on the following assumptions: each capillary is the sole oxygen supply for a surrounding cylindrical region of tissue, the PO2 at the vessel wall is assumed to equal that of the blood, the decline of PO2 along a capillary is neglected, oxygen diffuses radially from the capillary, and consumption is uniform in the tissue. In subsequent models (21), several of these assumptions were relaxed. For example, Blum (4) included axial decline in blood PO2 and considered both constant and linear oxygen uptake kinetics. The importance of intravascular resistance to oxygen transport from blood to tissue was shown by Hellums (14) and Hellums et al. (15). Theoretical estimates of the minimum PO2 required to achieve complete tissue oxygenation at maximal consumption in a variety of muscles and species were obtained by Roy and Popel (29).
Federspiel (9) developed a two-dimensional model of oxygen delivery to a muscle fiber with a prescribed distribution of PO2 at its surface. Groebe and Thews (12) considered a single cylindrical muscle fiber with several adjacent capillaries in an effort to evaluate the effects of model geometry on predicted PO2 profiles in heavily working muscle. Both studies predicted relatively shallow gradients of PO2 within muscle fibers, as observed by Gayeski and Honig (11). Secomb and Hsu (30) computed the PO2 at each point in a finite tissue domain, taking into account the effects of all microvessels within the domain. Their model showed that PO2 levels in capillaries in resting skeletal muscle can be strongly influenced by oxygen diffusion from arterioles at distances of 100 µm or more. As oxygen consumption rate increases, such relatively long-range effects become less important, and when consumption is very high, the assumption of the Krogh model, that each point in the tissue receives oxygen only from the nearest capillary, becomes increasingly justified. Under such conditions, the use of a model based on the Krogh cylinder concept is appropriate.
The goal of the present study is to use a theoretical model for oxygen transport from capillaries to exercising skeletal muscle to predict the dependence of oxygen consumption rate on oxygen demand, on the basis of a consideration of transport processes occurring at the microvascular level. A wide range of oxygen demand is considered, including conditions of high oxygen demand in which a significant fraction of the tissue becomes hypoxic (PO2 < 1 Torr). None of the previous models cited above has analyzed the dependence of oxygen consumption on demand under such conditions. Oxygen consumption is assumed to depend on PO2 according to Michaelis-Menten kinetics, and results are compared with those obtained when zero-order consumption kinetics is assumed. Nonuniform oxygen consumption due to mitochondrial clustering around the capillary is considered. Effects of the decline of oxygen content in blood flowing through capillaries, intravascular resistance to oxygen diffusion, and myoglobin-facilitated diffusion are included. As in the Krogh cylinder model, the capillaries are assumed to be parallel and evenly spaced, so that each capillary supplies oxygen to a surrounding uniform cylindrical region of tissue. A range of capillary densities is considered.
This model is used to investigate the roles played by oxygen demand,
perfusion limitation, and diffusion limitation in determining
O2 max of skeletal muscle. Profiles of
PO2 in the tissue surrounding the capillary are
determined and are compared with measurements of tissue
PO2 in maximally stimulated skeletal muscle. The overall average oxygen consumption rate as a function of demand is
calculated and compared with experimental
O2 max values in human skeletal muscle.
Previous studies (29) have explored the relationship
between oxygen transport at the capillary level and
O2 max. However, no theoretical results
including all the factors considered here have previously been compared with experimentally determined oxygen consumption rates. The results are used to test the hypothesis that oxygen consumption rates in
maximally stimulated skeletal muscle can be predicted on the basis of
observed values of capillary density, blood perfusion, and oxygen
transport parameters in blood and in tissue.
| |
METHODS |
|---|
|
|
|---|
Governing equations.
The tissue is represented as an array of uniformly spaced cylinders
with capillaries along the axes. Each tissue cylinder is assumed to be
supplied with oxygen exclusively by the capillary within it (Fig.
1). Oxygen diffusion in the axial
direction is neglected, as evidenced by the fact that the
PO2 gradients are much steeper in the radial
direction than in the axial direction. Then Fick's law of diffusion
and conservation of mass lead to the following equation for oxygen
diffusion in muscle tissue
|
(1) |
|
= Mt
[Pb
P(Rc)] where
is the rate of oxygen diffusion from the capillary per unit length,
Mt is a mass transfer coefficient,
Pb is the average partial pressure of oxygen within the
blood, and P(Rc) is the partial pressure of
oxygen at the capillary wall, r = Rc. Equating this to the diffusive flux into the
tissue gives the following boundary condition at the capillary wall
|
(2) |
|
(3) |
|
(4) |
, i.e.
|
(5) |
is blood flow rate and z is distance
along the capillary. The oxygen content of blood is given by C = CB SHb(Pb), where CB is
the carrying capacity of blood at 100% saturation and the
oxyhemoglobin saturation is described using the Hill equation (22)
|
(6) |
Parameter values.
Parameter values are chosen to represent blood flow in human skeletal
muscle under physiological conditions and are summarized in Table
1. At high flow rates, such as assumed
here, the decline in PO2 in arterioles is
small, and the PO2 entering the capillary is
assumed to equal that of the arterial blood, i.e., Pb = 100 Torr when z = 0. Experimentally determined values
for K in skeletal muscle range from 5 × 10
10 to 10 × 10
10
(cm2/s)(cm3
O2 · cm
3 · Torr
1).
Here, K = 9.4 × 10
10
(cm2/s)(cm3
O2 · cm
3 · Torr
1)
is assumed (2). To model the effects of myoglobin, the
following parameter values are assumed: DMb = 1.73 × 10
7 cm2/s (18);
CMb = 3.83 × 10
7
mol/cm3 (20); Vm = 2.24 × 104 cm3. The oxygen-myoglobin saturation
SMb(P) is represented by a Michaelis-Menten equation, with
half-saturation at a PO2 of
PMb50 = 3.2 Torr (26).
|
5 cm3
O2 · cm
3 · Torr
1,
the amount of free dissolved oxygen is <2% of the total amount of
oxygen in the blood when PO2 is between 15 and
100 Torr. Therefore, only hemoglobin-bound oxygen is considered.
M0 can vary over a wide range in skeletal
muscle. Experiments based on mitochondrial enzyme turnover imply that
M0 = 40 cm3
O2 · 100 cm
3 · min
1 in a maximally working
muscle (3), but observed consumption rates
(25) imply higher values. To show the effects of demand on
consumption, a range of M0 from 0 to 80 cm3 O2 · 100 cm
3 · min
1 is considered, and a
reference value of 40 cm3 O2 · 100 cm
3 · min
1 is used. The dependence
of oxygen consumption on PO2 is not precisely known. Here, Michaelis-Menten kinetics are assumed, with half-maximal consumption at a PO2 of P0 = 1 Torr. This value is in the range indicated by experimental studies
(6, 27). According to this assumption, oxygen consumption
is less than 50% of demand in hypoxic regions, where
PO2 <1 Torr. To assess the sensitivity of
the results to the Michaelis-Menten form of oxygen uptake kinetics and
the value chosen for P0, calculations were carried out
using zero-order kinetics, which corresponds to the limiting case of
Michaelis-Menten kinetics when P0
0. To model possible
effects of mitochondrial clustering around capillaries, calculations
were carried with an oxygen demand of 80 cm3
O2 · 100 cm
3 · min
1 in a cylindrical region
around the capillary representing half of the total tissue volume and
zero outside this region, so that the average demand is equal to the
reference level, i.e., 40 cm3 O2 · 100 cm
3 · min
1.
In mammalian skeletal muscle, typical capillary diameters are in the
range of 4 to 8 µm. Here, a diameter of 5 µm is assumed, so
Rc = 2.5 µm. Capillary lengths,
L, and average flow velocities,
, have not
been determined in exercising human skeletal muscle. However, the model
results depend only on their ratio, as shown by Eq. 5, which
is unchanged if z and
are altered by the same factor. The ratio
/L may be estimated from
the perfusion rate
|
(7) |
3 · s
1. For the
assumed capillary density and radius, this gives
/L = 4.5 s
1. For purposes of
illustration, we choose L = 0.5 mm and
= 2.25 mm/s, which are in the range of observed values in mammalian skeletal muscle. In the standard case, the tissue cylinder radius, Rt, is based on a capillary density in human
skeletal muscle tissue of 468 capillaries per square millimeter
(7). The Krogh tissue cylinder radius that gives the same
capillary density is Rt = 26 µm. However,
the training level of the subjects has been shown to have a significant
effect on capillary density. To compare model results with experimental
O2 max values, we consider a range of
capillary densities from 468 to 1,000 capillaries/mm2.
Intravascular resistance to oxygen diffusion depends on the mass
transfer coefficient, which is given by Mt =
Kpl Sh, where Kpl is the Krogh diffusion coefficient in plasma
and Sh is the Sherwood number, a nondimensional constant
that depends on the oxygen transport process occurring within the
vessel and reflects the particulate nature of blood. Hellums et al.
(15) compiled theoretical and experimental data showing
the dependence of Sh on vessel diameter. For a diameter of 5 µm, Sh ranging from 1.5 to 3 are consistent with available
data. Here, Sh = 2.5 is assumed, except where noted,
and Kpl = 8.3 × 10
10
(cm2/s) (cm3
O2 · cm
3 · Torr
1)
(15), giving Mt = 6.52 × 10
9 (cm2/s) (cm3
O2 · cm
3 · Torr
1).
Numerical procedures.
The capillary is discretized into 100 points along its length. At each
point, Eqs. 1, 3, and 4 are solved
numerically to determine the radial profile of
PO2. The total consumption per unit length is
found by numerically integrating the consumption per unit volume over
the cross section of the tissue cylinder, i.e.
|
(8) |
| |
RESULTS |
|---|
|
|
|---|
Figure 2 shows the predicted
variation of PO2 with position in the tissue
cylinder, including effects of intravascular resistance and
myoglobin-facilitated diffusion. The decline in the average blood
PO2 as it flows down the length of the
capillary is shown in Fig. 2A. Radial
PO2 profiles at three distances along the
capillary are shown in Fig. 2B. At each location along the
cylinder, PO2 declines with increasing distance
from the capillary. The PO2 in the tissue
adjacent to the capillary is substantially lower than the mean
PO2 within the capillary as a result of
intravascular resistance to radial oxygen transport. In reality,
intracapillary PO2 is continuous with the
tissue PO2 but varies with both position and
time as red blood cells traverse the capillary. For simplicity, only
the average intracapillary level is shown in Fig. 2B. At the
upstream end of the capillary (z = 0), the
PO2 remains relatively high throughout the
tissue, dropping from 78.9 Torr at the capillary wall to 35.5 Torr at
the outer boundary of the cylinder. However, PO2 declines rapidly with distance along the
capillary. At the midpoint (z = L/2), the
mean intravascular PO2 is 34.7 Torr, and the
tissue PO2 drops to 0.5 Torr at the outer edge
of the cylinder. At this point, much of the tissue is hypoxic, leading
to a reduced rate of oxygen consumption, according to Michaelis-Menten
kinetics. Beyond this point, intravascular PO2
declines more slowly, because less oxygen is being consumed per unit
capillary length.
|
Figure 3 shows cumulative frequency
distributions of tissue PO2 and myoglobin
saturation, when demand is 40 cm3
O2 · 100 cm
3 · min
1. Figure 3A
indicates that a large fraction of tissue has low tissue
PO2. With 468 capillaries/mm2 and
P50 = 26 Torr, 37% of the tissue is hypoxic
(PO2 < 1 Torr) and more than 90% of the
tissue has a PO2 less than 12.5 Torr. Increasing capillary density to 600 capillaries/mm2 reduces
the fraction of tissue that is hypoxic to 20%. If a right shift in the
oxyhemoglobin dissociation curve is considered (by linearly increasing
P50 from 26 Torr at the arterial end to 39 Torr at the
venous end of the capillary) in addition to an increased capillary
density of 600 capillaries/mm2, 1.4% of the tissue is
hypoxic. The corresponding distributions of myoglobin
saturation are shown in Fig. 3B, indicating that saturation
is fairly evenly distributed over a wide range.
|
The relationship between oxygen supply, demand, and consumption is
shown in Fig. 4. The rate of convective
oxygen supply to the tissue cylinder is 49.2 cm3
O2 · 100 cm
3 · min
1, and this represents an
upper limit to consumption. The demand is the rate at which oxygen
would be consumed if all the tissue was well oxygenated. At low levels
of demand (M0 < 10 cm3
O2 · 100 cm
3 · min
1), the consumption is
equal to the demand, but at higher levels, consumption falls short of
demand. The limited ability of oxygen to diffuse into the tissue
restricts the consumption that is achieved. Predicted consumption
decreases with increasing intravascular resistance, i.e., with
decreasing Sh. In the absence of intravascular resistance to
oxygen diffusion, a demand of 40 cm3
O2 · 100 cm
3 · min
1 results in a predicted
consumption rate of 30.3 cm3 O2 · 100 cm
3 · min
1. Inclusion of the
effects of intravascular resistance to oxygen transport further reduces
the predicted rate of consumption. For example, when the demand is 40 cm3 O2 · 100 cm
3 · min
1, the predicted
consumption is 24.7 cm3 O2 · 100 cm
3 · min
1 for the standard value
of intravascular resistance (Sh = 2.5). The effect of
varying the intravascular resistance is also shown in Fig. 4. The limit
Sh
corresponds to the case of no intravascular resistance.
|
To assess the sensitivity of the results to the assumed
Michaelis-Menten oxygen uptake kinetics, further calculations were carried out with the assumption that local consumption equals demand
until the PO2 falls to zero (zero-order
kinetics). The region of tissue that receives oxygen is smaller with
zero-order kinetics than with Michaelis-Menten, but the consumption
rate is higher. Overall, zero-order kinetics increase oxygen
consumption by 3% (from 24.7 to 25.5 cm3
O2 · 100 cm
3 · min
1) relative to
Michaelis-Menten kinetics when M0 = 40 cm3 O2 · 100 cm
3 · min
1.
The diffusive resistance to oxygen transport depends on the capillary
density, N. For a given level of oxygen demand, increasing N decreases the amount of oxygen that each capillary must
deliver, and so a smaller PO2 gradient is
needed to produce the necessary flux. Furthermore, higher N
corresponds to a smaller Rt and therefore a
shorter distance over which oxygen must diffuse. Consequently, predicted oxygen consumption increases substantially with increasing N, as shown in Fig. 5.
Increasing N from 468 to 800 capillaries/mm2
increases predicted consumption by 33%, to 32.8 cm3
O2 · 100 cm
3 · min
1 when
M0 = 40 cm3
O2 · 100 cm
3 · min
1.
|
A conceptual model has previously been proposed (33) in
which
O2 max is determined by the
condition that the PO2 at the venous end of
each capillary equals the minimum PO2 required for complete oxygenation of the adjacent tissue. This concept is
explored in Fig. 6, which shows the
variation of predicted mean intravascular PO2
at the venous end of capillaries with oxygen consumption rate (solid
line). Because the perfusion rate is assumed to be fixed, venous
PO2 drops with increasing consumption rate. The
shape of this curve is dictated by conservation of mass and does not
depend on the distribution of tissue PO2. The
predicted consumption that can be achieved when demand is 40 cm3 O2 · 100 cm
3 · min
1 corresponds to a venous
PO2 of 25.5 Torr (or 24.2 Torr if mitochondrial clustering is assumed). The dashed curve in Fig. 6 shows the minimum blood PO2 that is required to supply oxygen to
the outer boundary of the tissue cylinder (r = Rt) if a uniform consumption rate is assumed.
The intersection of the two curves represents the point at which the
blood PO2 at the downstream end of the
capillary is equal to the minimum PO2 required
to meet MO. This implies that the maximum consumption rate
that could be achieved assuming a uniform consumption rate with no
hypoxia throughout the tissue cylinder is 19.0 cm3
O2 · 100 cm
3 · min
1. According to the
present model, significantly higher overall consumption rates can be
achieved if demand is increased and a portion of the tissue is
permitted to become hypoxic.
|
Predicted rates of oxygen consumption, under several different sets of
assumptions, are compared in Fig. 7 with
experimentally measured
O2 max values.
Experimental results are indicated by solid bars, and the other bars
represent theoretical predictions. For each, the shaded region shows
consumption when M0 = 40 cm3
O2 · 100 cm
3 · min
1, the hatched region
shows the additional consumption when mitochondrial clustering near the
capillary is assumed, and the white region shows the additional
consumption when M0 is increased to 80 cm3 O2 · 100 cm
3 · min
1 throughout the tissue.
The leftmost bars correspond to the case in which effects of
intravascular resistance and myoglobin-facilitated diffusion are
neglected. Subsequent bars show the cumulative results of successively
including intravascular resistance, myoglobin-facilitated diffusion, a
right shift in the oxyhemoglobin dissociation curve (P50
increasing linearly from 26 Torr at the entrance to 39 Torr at the
venous end of the capillary), and increased capillary density (as
discussed below).
|
The predictions shown in Fig. 7A correspond to the results
of Andersen and Saltin (1) for the human quadriceps
muscle, with a perfusion rate of 0.041 cm3
blood · cm
3 · s
1,
corresponding to an oxygen supply of 49.2 cm3
O2 · 100 cm
3 · min
1. The third bar (+Myo) in
Fig. 7A, including effects of intravascular resistance and
myoglobin-facilitated diffusion, corresponds to the standard case shown
in Figs. 2, 4, and 5. Inclusion of intravascular resistance reduces
predicted consumption significantly, from 29.9 to 24.3 cm3
O2 · 100 cm
3 · min
1 when
M0 = 40 cm3
O2 · 100 cm
3 · min
1 (model and +IVR bars).
Myoglobin-facilitated diffusion leads to only a small (<2%) increase
in consumption to 24.7 cm3 O2 · 100 cm
3 · min
1 for
M0 = 40 cm3
O2 · 100 cm
3 · min
1 (+Myo bar). This level
of consumption is substantially lower than the reported
O2 max value of 35.0 cm3
O2 · 100 cm
3 · min
1 (1). A
right shift in the oxyhemoglobin dissociation curve (+RS bar) increases
consumption by almost 14%. Predicted consumption with
M0 = 80 cm3
O2 · 100 cm
3 · min
1 is then 13% less than
the observed
O2 max. To obtain consumption levels close to the observed
O2 max, a capillary density higher than
468 capillaries/mm2 must be assumed. For example,
increasing capillary density to 600 capillaries/mm2 (+HD
bar) leads to predicted consumption rates of 32.1 cm3
O2 · 100 cm
3 · min
1 for
M0 = 40 cm3
O2 · 100 cm
3 · min
1, 33.8 cm3
O2 · 100 cm
3 · min
1 with mitochondrial
clustering, and 36.1 cm3 O2 · 100 cm
3 · min
1 when
M0 = 80 cm3
O2 · 100 cm
3 · min
1, which are close to the
observed level.
In the study by Richardson et al. (25), the
perfusion in the quadriceps muscle was 0.0642 cm3
blood · cm
3 · s
1,
corresponding to an oxygen supply of 75.7 cm3
O2 · 100 cm
3 · min
1. Predictions for this
higher perfusion level are shown in Fig. 7B. Relative to
Fig. 7A, oxygen supply is increased by ~54%, but the
increase in consumption is less, ~12% for the case including effects
of intravascular resistance and myoglobin (gray +Myo bar). Inclusion of
intravascular resistance, myoglobin-facilitated diffusion, and right
shifting the oxyhemoglobin dissociation curve lead to similar effects
to those seen in Fig. 7A. In this case, even when a right
shift is included, the predicted consumption rates (31.4 cm3 O2 · 100 cm
3 · min
1 when
M0 = 40 cm3
O2 · 100 cm
3 · min
1, 33.4 cm3
O2 · 100 cm
3 · min
1 with clustering, and
34.9 cm3 O2 · 100 cm
3 · min
1 when
M0 = 80 cm3
O2 · 100 cm
3 · min
1) fall far short of the
reported
O2 max, which is 60.2 cm3 O2 · 100 cm
3 · min
1 (25).
According to the model, such high oxygen consumption can only be
achieved by a combination of high demand and high N.
Clearly, this reported consumption rate cannot be achieved with a
demand of 40 cm3 O2 · 100 cm
3 · min
1. If N is
increased from 468 to 1,000 capillaries/mm2 and a demand of
40 cm3 O2 · 100 cm
3 · min
1 is considered, nearly
all of the tissue is well oxygenated, and the predicted consumption,
38.8 cm3 O2 · 100 cm
3 · min
1, nearly equals demand.
As seen in the +HD bar, under these conditions, mitochondrial
clustering has practically no effect on oxygen consumption, and
predicted consumption still falls well short of the reported
O2 max value. However, if an increased
N is coupled with a demand of 80 cm3
O2 · 100 cm
3 · min
1, the resulting
consumption is 57.4 cm3 O2 · 100 cm
3 · min
1, close to the observed
O2 max.
| |
DISCUSSION |
|---|
|
|
|---|
Factors determining oxygen consumption rate.
When M0 is high in skeletal muscle, the rate of
oxygen consumption depends not only on the demand, but also on the rate
of convective delivery and on the diffusive processes occurring within capillaries and in the surrounding tissue. The present analysis shows
that all these factors have significant effects on the consumption rate
that is achieved. Clearly, consumption cannot exceed either the demand
of the muscle or the convective supply (Fig. 4). When demand is low,
all tissue has sufficient oxygen, and predicted consumption equals
demand. As demand increases, the limited rate of oxygen diffusion leads
to low oxygen levels in the outer regions of the tissue cylinder, as
shown schematically in Fig. 1. This oxygen deficiency causes the
average consumption within the cylinder to fall below the demand.
Further increases in demand, beyond ~40 cm3
O2 · 100 cm
3 · min
1, lead to only slight
increases in consumption (Fig. 4).
3 · min
1 before leveling off
(Fig. 5).
The role of myoglobin in oxygen delivery has been controversial.
Gayeski and Honig (11) considered that myoglobin was an important factor in oxygen transport within muscle fibers. Conley et
al. (5) used measurements of myoglobin saturation in vivo made with magnetic resonance spectroscopy to support the opinion that
myoglobin may play a role in facilitating oxygen transport in skeletal
muscle but did not quantify the effects of myoglobin on oxygen
transport. In a recent review, Jürgens et al. (17) concluded that myoglobin plays at most a minor role in oxygen transport
under physiological conditions. A similar conclusion was reached by Roy
and Popel (29). According to the present analysis,
inclusion of myoglobin-facilitated diffusion has relatively small
effects on overall oxygen transport. It leads to shallower PO2 gradients, lower intravascular
PO2 values, higher PO2
levels in the outer part of the cylinder where tissue is hypoxic, and a
slight (<2%) net increase in consumption (Fig. 7, +Myo). To understand this result, it is helpful to consider the facilitation pressure (Pf) = DmbCmbVm/K,
which represents the contribution of myoglobin to oxygen transport in
terms of an effective additional partial pressure of oxygen available
to drive diffusive transport (10). For the parameter
values used here (DMb = 1.73 × 10
7 cm2/s), Pf = 1.6 Torr,
which is much less than the capillary PO2 even
at the venous end. This result reflects the fact that the diffusivity
of myoglobin in muscle is two orders of magnitude less than the
diffusivity of oxygen. Therefore, myoglobin-facilitated oxygen
transport would not be expected to have much effect. Some earlier
models were based on the diffusivity of myoglobin measured in protein
solutions (18), DMb = 8 × 10
7 cm2/s. Assuming this value gives
Pf = 7.3 Torr and leads to a 5% increase in the oxygen
consumption rate when M0 = 40 cm3 O2 · 100 cm
3 · min
1.
The effects of a right shift of the oxyhemoglobin dissociation curve
and mitochondrial clustering can also increase consumption, as
indicated in Fig. 7. A right shift of the oxyhemoglobin dissociation curve increases blood PO2 at a given saturation
and therefore increases the PO2 gradient for
oxygen diffusion into tissue. Simulating a right shift by linearly
increasing P50 from 26 to 39 Torr along the length of the
capillary leads to a significantly increased consumption (Fig. 7, +RS).
Additional simulations in which P50 varied smoothly from 26 to 39 Torr, but with nonlinear dependence on distance traveled, showed
very similar results. Effects of mitochondrial clustering are
represented here by assuming that consumption occurs only in a region
surrounding each capillary representing half of the tissue volume. For
an average consumption rate 40 cm3
O2 · 100 cm
3 · min
1, clustering leads to a
significant increase in consumption (Fig. 7, hatched areas) whenever
hypoxic regions are present. When N is increased to 1,000 capillaries/mm2 (Fig. 7B, +HD), virtually all
the tissue is well oxygenated at a uniform demand of 40 cm3
O2 · 100 cm
3 · min
1 and clustering has
almost no effect.
Model predictions are sensitive to N, which is equal to the
capillary-to-fiber ratio divided by the average fiber cross-sectional area. Changes in N can be caused by changes in the number of
capillaries per fiber or in fiber size. Because the model assumes that
each capillary supplies oxygen exclusively to the surrounding cylinder of tissue, regardless of the number of fibers intersected by that cylinder, the results are not sensitive to the capillary-to-fiber ratio
or the fiber size per se, but only to N (i.e.,
capillaries/mm2).
Some effects that are not considered in this model may also influence
consumption rates. Only oxygen diffusion from capillaries is
considered. In resting muscle, a significant fraction of the M0 can be met by oxygen diffusion from
arterioles. However, this fraction decreases with increasing
M0 (30). In the model, the contribution of oxygen dissolved in the plasma to convective delivery is neglected. Additional calculations including this component were
carried out, showing that predicted consumption increased by <1.1%.
Capillaries are lined with a glycocalyx or endothelial surface layer
0.5-1 µm thick that can exclude flowing red blood cells
(23). This layer may increase intravascular
resistance to oxygen diffusion, but its effects are unknown. The
assumed Sherwood number, Sh = 2.5, was based on
theoretical estimates that neglect such effects (15).
Effects of nonuniform capillary spacing and flow.
The Krogh-type model assumes that all capillaries in the tissue are
identical, with the same flow in each and uniform spacing between them.
In reality, capillaries are not evenly spaced, and this heterogeneity
may lead to reduced oxygen delivery, because capillaries in more
densely perfused regions may flow through regions of higher
PO2 where less oxygen is extracted. Under
conditions of very high demand, however, even relatively densely
perfused areas can be assumed to contain a significant amount of
hypoxic tissue, particularly near the downstream end of each capillary, so that extraction from each capillary is not greatly affected by the
presence of near neighbors. According to this argument, the effect of
nonuniform spacing on average consumption is expected to be small when
demand is very high. Heterogeneity in capillary length or flow rate
could also affect the results. To estimate the effects of variations in
capillary lengths, two cylinders were considered, one 



Relationship to venous PO2.
According to the conceptual model of Wagner (33), for a
given perfusion rate, the PO2 at the venous end
of capillaries declines with increasing consumption, whereas the
PO2 required to achieve tissue oxygenation
increases. Wagner suggested that
O2 max is reached when venous PO2 falls to the minimum
needed for complete tissue oxygenation. This condition is indicated
graphically in Fig. 6 by the intersection of the two curves. Wagner
defines
O2 max as the point at which
venous blood PO2 is equal to the
PO2 required to fully meet the
M0 of the tissue, i.e., no hypoxic regions are present. However, as shown in Fig. 6, significantly higher overall oxygen consumption rates can be achieved if oxygen demand increases beyond the level at which hypoxia first appears. With further increases
in demand, consumption in the well-oxygenated regions increases,
whereas the hypoxic regions spread and consumption decreases there
according to Michaelis-Menten kinetics. The decrease in consumption in
hypoxic regions is more than compensated for by the increased
consumption in well-oxygenated regions, and the net effect is an
increase in consumption. For example, the curves in Fig. 6 intersect
when demand is 19 cm3 O2 · 100 cm
3 · min
1 and consumption is 17.9 cm3 O2 · 100 cm
3 · min
1. Increasing demand to 40 cm3 O2 · 100 cm
3 · min
1 results in a
significantly higher consumption rate, 24.7 cm3
O2 · 100 cm
3 · min
1. Therefore, Wagner's
conceptual model, if applied quantitatively, may underestimate the
maximal oxygen consumption rate. Similarly, the predictions of
end-capillary PO2 at
O2 max by Roy and Popel
(29), which were mainly in the range 20-30 Torr, were based on the assumption that no hypoxia is present and correspond to
the intersection of the two curves in Fig. 6. According to the present
model, their approach may overestimate end-capillary PO2 at
O2 max
by several millimeters of mercury.
Comparison with observations of tissue oxygen levels.
Gayeski and Honig (11) used cryospectrophotometry to
measure myoglobin saturation of dog gracilis muscle rapidly frozen at
O2 max. Measured myoglobin saturation
levels were between 8% and 40%, corresponding to
PO2 values of 0.5 Torr and 3.5 Torr, based on
PMb50 = 5.3 Torr for dog muscle. They reported shallow
PO2 gradients of <0.1 Torr/µm and low,
nonzero PO2 values throughout the tissue at
distances greater than 3 µm from the nearest capillary. They
concluded that most of the drop in PO2 between
red blood cells in a capillary and an adjacent muscle fiber occurs
outside the fiber or within a few micrometers inside the fiber. These
results were stated to contradict classical theories of oxygen
transport to tissue such as the Krogh model and were attributed to the
effects of intravascular resistance and myoglobin facilitation. In
subsequent studies, Voter and Gayeski (32) showed that the
sampling region of the earlier studies was larger than assumed by
Gayeski and Honig (11), implying that the oxygen gradients
may have been larger than originally estimated. However, these studies
supported the previous conclusion that oxygen gradients are shallow
with low PO2 throughout the tissue at maximal consumption.
Comparison with measured values of
O2 max.
Oxygen consumption depends on M0. On the basis
of the experimental work of Blomstrand et al. (3) and
Richardson et al. (25), a range of
M0 from 0-80 cm3
O2 · 100 cm
3 · min
1 was considered. In the
absence of a measured value for M0, no unique
prediction of
O2 max can be made.
However, predicted consumption increases slowly with
M0 in the range from 40 to 80 cm3
O2 · 100 cm
3 · min
1 (Fig. 4). Therefore,
consumption rates for M0 values in this range
can be considered as approximate predictions of
O2 max.
O2 max, independent of systemic
limitations on oxygen delivery. Results using this approach (1,
25) are compared in Fig. 7 with model predictions. With standard
parameter values (N = 468 capillaries/mm2,
P50 = 26 Torr, including intravascular resistance and
myoglobin), the model underestimates the measured
O2 max values. However, the values of
some key parameters were not measured in the experimental studies, and
it is of interest to examine what values of these parameters would lead
to consumption rates close to measured
O2 max values.
Andersen and Saltin (1) observed an average
O2 max of 35 cm3
O2 · 100 cm
3 · min
1 in subjects who ranged
in training from sedentary to endurance trained. In the corresponding
model results (Fig. 7A), when a right shift in the
oxyhemoglobin dissociation curve is included by linearly increasing
P50 from 26 Torr at the arterial end to 39 Torr at the
venous end of the capillary, consumption with
M0 = 40 cm3
O2 · 100 cm
3 · min
1 increases from 24.7 to
28.0 cm3 O2 · 100 cm
3 · min
1 but is still less than
the measured value. With mitochondrial clustering and a right shift,
the consumption is 29.9 cm3 O2 · 100 cm
3 · min
1, and when
M0 = 80 cm3
O2 · 100 cm
3 · min
1 the consumption with the
right shift is 30.6 cm3 O2 · 100 cm
3 · min
1. Assuming a higher
capillary density, N = 600 capillaries/mm2,
leads to a predicted consumption of 32.1 cm3
O2 · 100 cm
3 · min
1 when
M0 = 40 cm3
O2 · 100 cm
3 · min
1, close to the measured
value. Therefore, the observed
O2 max can be accounted for by modifying the standard model to include a right
shift of the oxyhemoglobin dissociation curve and a slightly higher
capillary density.
The subjects in the experiments of Richardson et al. (25)
were competitive endurance bicyclists, and
O2 max values with a mean of 60.2 cm3 O2 · 100 cm
3 · min
1 were obtained, much
higher than any model predictions with a capillary density of 468 capillaries/mm2. However, if capillary density is increased
to 1,000 capillaries/mm2 and a right shift is included,
with M0 = 80 cm3
O2 · 100 cm
3 · min
1, a consumption rate of
57.4 cm3 O2 · 100 cm
3 · min
1 is predicted, close to
the measured value. Endurance training has been shown to have many
effects on skeletal muscle, including increasing N and
mitochondrial number and causing an even more pronounced right shift in
the oxyhemoglobin dissociation curve during exercise when compared with
untrained subjects. The highest N in human skeletal muscle
found in the literature was 810 capillaries/mm2
(8). N estimated by staining methods may
substantially underestimate the true N (13).
According to the model results, a combination of these effects,
including a marked increase in N, must be present in the
muscles of the subjects studied by Richardson et al. (25) to account for the very high muscle-specific
O2 max levels observed.
In conclusion, a theoretical model for oxygen transport to tissue has
been developed on the basis of the tissue cylinder concept of Krogh
(19) but incorporating recent information on oxygen transport mechanisms and parameters. According to the model, oxygen consumption rates under conditions of high demand depend on both convective and diffusive limitations on oxygen delivery. Intravascular resistance to oxygen diffusion significantly restricts delivery, whereas myoglobin-facilitated diffusion in tissue plays a negligible role. A right shift of the oxyhemoglobin dissociation curve can significantly enhance delivery, and consumption rates are sensitively dependent on N. When demand is high, much of the tissue is
predicted to be hypoxic, with low PO2
gradients, as observed experimentally.
Predicted
O2 max are comparable in
magnitude to experimentally observed
O2 max values. Andersen and Saltin (1) considered subjects whose level of physical activity
ranged from sedentary to endurance trained. The average value
determined for these subjects can be predicted with this model using
capillary densities found in the literature. Richardson et al.
(25) reported a much higher average
O2 max value for endurance trained athletes. This value can only be explained with the model by assuming Ns near 1,000 capillaries/mm2.
| |
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: T. W. Secomb, Dept. of Physiology, Univ. of Arizona, Tucson, AZ 85724-5051 (E-mail: secomb{at}u.arizona.edu).
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.
Received 2 November 2000; accepted in final form 16 July 2001.
| |
REFERENCES |
|---|
|
|
|---|
1.
Andersen, P,
and
Saltin B.
Maximal perfusion of skeletal muscle in man.
J Physiol (Lond)
366:
233-249,
1985
2.
Bentley, TB,
Meng H,
and
Pittman RN.
Temperature dependence of oxygen diffusion and consumption in mammalian striated muscle.
Am J Physiol Heart Circ Physiol
264:
H1825-H1830,
1993
3.
Blomstrand, E,
Radegran G,
and
Saltin B.
Maximum rate of oxygen uptake by human skeletal muscle in relation to maximal activities of enzymes in the Krebs cycle.
J Physiol (Lond)
501:
455-460,
1997
4.
Blum, J.
Concentration profiles in and around capillaries.
Am J Physiol
198:
991-998,
1960.
5.
Conley, KE,
Ordway GA,
and
Richardson RS.
Deciphering the mysteries of myoglobin in striated muscle.
Acta Physiol Scand
168:
623-634,
2000[Web of Science][Medline].
6.
Costa, LE,
Mendez G,
and
Boveris A.
Oxygen dependence of mitochondrial function measured by high-resolution respirometry in long-term hypoxic rats.
Am J Physiol Cell Physiol
273:
C852-C858,
1997
7.
Coyle, EF,
Martin WH,
Sinacore DR,
Joyner MJ,
Hagberg JM,
and
Holloszy JO.
Time course of loss of adaptations after stopping prolonged intense endurance training.
J Appl Physiol
57:
1857-1864,
1984
8.
Crenshaw, AG,
Friden J,
Thornell LE,
and
Hargens AR.
Extreme endurance training: evidence of capillary and mitochondria compartmentalization in human skeletal muscle.
Eur J Appl Physiol
63:
173-178,
1991[Web of Science].
9.
Federspiel, WJ.
A model study of intracellular oxygen gradients in a myoglobin-containing skeletal muscle fiber.
Biophys J
49:
857-868,
1986[Web of Science][Medline].
10.
Fletcher, JE.
On facilitated diffusion in muscle tissues.
Biophys J
29:
437-458,
1980[Web of Science][Medline].
11.
Gayeski, TE,
and
Honig CR.
O2 gradients from sarcolemma to cell interior in red muscle at maximal
O2.
Am J Physiol Heart Circ Physiol
251:
H789-H799,
1986
12.
Groebe, K,
and
Thews G.
Role of geometry and anisotropic diffusion for modeling PO2 profiles in working red muscle.
Respir Physiol
79:
225-278,
1990.
13.
Hansen-Smith, FM,
Blackwell LH,
and
Joswiak GR.
Expression of muscle capillary alkaline phosphatase is affected by hypoxia.
J Appl Physiol
73:
776-780,
1992
14.
Hellums, JD.
The resistance to oxygen transport in the capillaries relative to that in the surrounding tissue.
Microvasc Res
13:
131-136,
1977[Web of Science][Medline].
15.
Hellums, JD,
Nair PK,
Huang NS,
and
Ohshima N.
Simulation of intraluminal gas transport processes in the microcirculation.
Ann Biomed Eng
24:
1-24,
1996[Web of Science][Medline].
16.
Hoppeler, H,
Mathieu O,
Krauer R,
Classen H,
Armstrong RB,
and
Weibel ER.
Design of the mammalian respiratory system. VI. Distribution of mitochondria and capillaries in various muscles.
Respir Physiol
44:
87-111,
1981[Web of Science][Medline].
17.
Jürgens, KD,
Papadopoulos S,
Peters T,
and
Gros G.
Myoglobin: just an oxygen store or also an oxygen transporter?
News Physiol Sci
15:
269-274,
2000
18.
Jürgens, KD,
Peters T,
and
Gros G.
Diffusivity of myoglobin in intact skeletal muscle cells.
Proc Natl Acad Sci USA
91:
3829-3833,
1994
19.
Krogh, A.
The number and the distribution of capillaries in muscle with the calculation of the oxygen pressure necessary for supplying the tissue.
J Physiol (Lond)
52:
409-515,
1919.
20.
Meng, H,
Bentley TB,
and
Pittman RN.
Myoglobin content of hamster skeletal muscle.
J Appl Physiol
74:
2194-2197,
1993
21.
Middleman, S.
Transport Phenomena in the Cardiovascular System. New York: Wiley-Interscience, 1972.
22.
Popel, AS.
Theory of oxygen transport to tissue.
Crit Rev Biomed Eng
17:
257-321,
1989[Web of Science][Medline].
23.
Pries, AR,
Secomb TW,
and
Gaehtgens P.
The endothelial surface layer.
Pflügers Arch
440:
653-666,
2000[Web of Science][Medline].
24.
Richardson, R, S,
Poole DC,
Knight DR,
Kurdak SS,
Hogan MC,
Grassi B,
Johnson EC,
Kendrick KF,
Erickson BK,
and
Wagner PD.
High muscle blood flow in man: is maximal O2 extraction compromised?
J Appl Physiol
75:
1911-1916,
1993
25.
Richardson, RS,
Knight DR,
Poole DC,
Kurdak SS,
Hogan MC,
Grassi B,
and
Wagner PD.
Determinants of maximal exercise VO2 during single leg knee-extensor exercise in humans.
Am J Physiol Heart Circ Physiol
268:
H1453-H1461,
1995
26.
Richardson, RS,
Noyszewski EA,
Kendrick KF,
Leigh JS,
and
Wagner PD.
Myoglobin O2 desaturation during exercise. Evidence of limited O2 transport.
J Clin Invest
96:
1916-1926,
1995.
27.
Richmond, KN,
Shonat RD,
Lynch RM,
and
Johnson PC.
Critical PO2 of skeletal muscle in vivo.
Am J Physiol Heart Circ Physiol
277:
H1831-H1840,
1999
28.
Rowell, LB,
Saltin B,
Kiens B,
and
Christensen NJ.
Is peak quadriceps blood flow in humans even higher during exercise with hypoxemia?
Am J Physiol Heart Circ Physiol
251:
H1038-H1044,
1986
29.
Roy, TK,
and
Popel AS.
Theoretical predictions of end-capillary PO2 in muscles of athletic and nonathletic animals at
O2 max.
Am J Physiol Heart Circ Physiol
271:
H721-H737,
1996
30.
Secomb, TW,
and
Hsu R.
Simulation of O2 transport in skeletal muscle: diffusive exchange between arterioles and capillaries.
Am J Physiol Heart Circ Physiol
267:
H1214-H1221,
1994
31.
Thomson, JM,
Dempsey JA,
Chosy LW,
Shahidi NT,
and
Reddan WG.
Oxygen transport and oxyhemoglobin dissociation during prolonged muscular work.
J Appl Physiol
37:
658-664,
1974
32.
Voter, WA,
and
Gayeski TE.
Determination of myoglobin saturation of frozen specimens using a reflecting cryospectrophotometer.
Am J Physiol Heart Circ Physiol
269:
H1328-H1341,
1995
33.
Wagner, PD.
An integrated view of the determinants of maximum oxygen uptake.
In: Oxygen Transfer from Atmosphere to Tissues, , edited by Gonzalez N,
and Fedde MR.. New York: Plenum, 1988, p. 245-256.
This article has been cited by other articles:
![]() |
J. C. Arciero, B. E. Carlson, and T. W. Secomb Theoretical model of metabolic blood flow regulation: roles of ATP release by red blood cells and conducted responses Am J Physiol Heart Circ Physiol, October 1, 2008; 295(4): H1562 - H1571. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. I. Levy, E. L. Schiffrin, J.-J. Mourad, D. Agostini, E. Vicaut, M. E. Safar, and H. A.J. Struijker-Boudier Impaired Tissue Perfusion: A Pathology Common to Hypertension, Obesity, and Diabetes Mellitus Circulation, August 26, 2008; 118(9): 968 - 976. [Full Text] [PDF] |
||||
![]() |
M. J. Lyon, L. M. Steer, and L. T. Malmgren Stereological estimates indicate that aging does not alter the capillary length density in the human posterior cricoarytenoid muscle J Appl Physiol, November 1, 2007; 103(5): 1815 - 1823. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Liguzinski and B. Korzeniewski Oxygen delivery by blood determines the maximal VO2 and work rate during whole body exercise in humans: in silico studies Am J Physiol Heart Circ Physiol, July 1, 2007; 293(1): H343 - H353. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Goldman, R. M. Bateman, and C. G. Ellis Effect of decreased O2 supply on skeletal muscle oxygenation and O2 consumption during sepsis: role of heterogeneous capillary spacing and blood flow Am J Physiol Heart Circ Physiol, June 1, 2006; 290(6): H2277 - H2285. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. D. Noakes, P. E. di Prampero, C. Capelli, T. Zaobornyj, L. B Valdez, A. Boveris, M. Ashenden, T. W. Secomb, S. Dufour, E. Ponsot, et al. Comments on Point:Counterpoint "Positive effects of intermittent hypoxia (live high:train low) on exercise performance are/are not mediated primarily by augmented red cell volume" J Appl Physiol, December 1, 2005; 99(6): 2453 - 2462. [Full Text] [PDF] |
||||
![]() |
W. J. van der Laarse, A. L. des Tombe, B. J. van Beek-Harmsen, M. B. E. Lee-de Groot, and R. T. Jaspers Krogh's diffusion coefficient for oxygen in isolated Xenopus skeletal muscle fibers and rat myocardial trabeculae at maximum rates of oxygen consumption J Appl Physiol, December 1, 2005; 99(6): 2173 - 2180. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Frisbee Reduced nitric oxide bioavailability contributes to skeletal muscle microvessel rarefaction in the metabolic syndrome Am J Physiol Regulatory Integrative Comp Physiol, August 1, 2005; 289(2): R307 - R316. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. J. McGuire and T. W. Secomb Estimation of capillary density in human skeletal muscle based on maximal oxygen consumption rates Am J Physiol Heart Circ Physiol, December 1, 2003; 285(6): H2382 - H2391. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. N. Bruce and M. C. Bruce A multicompartment model of carboxyhemoglobin and carboxymyoglobin responses to inhalation of carbon monoxide J Appl Physiol, September 1, 2003; 95(3): 1235 - 1247. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Lo, A. J. Fuglevand, and T. W. Secomb Oxygen delivery to skeletal muscle fibers: effects of microvascular unit structure and control mechanisms Am J Physiol Heart Circ Physiol, August 7, 2003; 285(3): H955 - H963. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Takahashi and K. Asano Mitochondrial respiratory control can compensate for intracellular O2 gradients in cardiomyocytes at low PO2 Am J Physiol Heart Circ Physiol, September 1, 2002; 283(3): H871 - H878. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |