|
|
||||||||
Departments of 1 Physiology and Biophysics and 2 Medicine, University of Washington, Seattle, Washington 98195
| |
ABSTRACT |
|---|
|
|
|---|
Although several factors are known to influence nonuniformity of ventilation, including lung mechanical properties (regional structure and compliance), external factors (chest wall, pleural pressure, heart), and ventilatory parameters (tidal and preinspiratory volume, flow rate), their relative contributions are poorly understood. We studied five excised, unperfused, canine right-middle lobes under varied levels of tidal volume (VT), thus eliminating many factors affecting heterogeneity. Multiple-breath washouts of N2 were analyzed for anatomic dead space volume (VDanat), nonuniformity of N2 washout, and nonuniformity between joined acinar regions vs. that occurring between larger joined regions. Approximately 80% of ventilation heterogeneity was found among joined acinar regions at resting levels of VT, but increasing VT reduced intra-acinar heterogeneity to about 25% of that found at resting levels. Increasing VT had essentially no effect on VDanat and heterogeneity among larger joined regions. The results indicate that the magnitude of VT is a major influence on the dominant intra-acinar component of ventilation heterogeneity and that VT effects on VDanat are likely due to perfusion and/or influences normally external to the lobar structure.
intraregional ventilation distribution; phase III slope; multiple breath inert-gas washout; normalized phase III slope
| |
INTRODUCTION |
|---|
|
|
|---|
INCREASING TIDAL VOLUME (VT) has long been known to increase exhaled anatomic dead space volume (VDanat) and alter the uniformity of ventilation distribution (1, 6, 15, 24, 27). Human studies suggest that increasing VT tends to increase the uniformity of gas distribution among very small joined airway regions on the scale of acini (6, 24), whereas ventilation between larger regions becomes more uneven (6). Quantitative and mechanistic understanding of these VT-dependent changes remains incomplete, in part because most evidence is from studies of intact animals and humans in which many influences on distribution of ventilation act simultaneously. In addition, methodological constraints have limited measurement of intra-acinar ventilation heterogeneity, and this component may dominate the nonuniformity of gas distribution in the whole lung (8). Also, although it is generally accepted that VT-dependent changes in VDanat are due to enlargement of proximal airways (1), computer models designed to test mechanisms of gas mixing within peripheral or proximal airways also predict VT effects on VDanat (9, 17, 25), and the relative contributions of these mechanisms have not been clearly demonstrated.
We undertook a reductionist approach to this problem by measuring ventilation heterogeneity in excised, suspended canine right-middle lobes under tightly controlled conditions of negative-pressure breathing and varied VT. Other ventilatory parameters, including preinspiratory volume, inspiratory time (TI), expiratory time (TE), and volume history, were held constant. This preparation removed most extralobar influences on ventilation distribution, including chest wall and diaphragm, the beating heart, and the normal gravity-determined pleural pressure gradient. Presently, the most sensitive and advanced method of assessing the uniformity of ventilation distribution is normalized phase III slope (SnIII) analysis of exhaled inert gas from multiple breath washout (MBW) (12, 13). A biexponential curve-fit approach to SnIII analysis was developed for these studies (see APPENDIX A) and used to separate diffusive-convective-dependent inhomogeneity (dcdi) of gas mixing that occurs among joined acinar airways (intraregional) from that which occurs between airways separating larger regions (interregional) due to only convective-dependent inhomogeneity (cdi). VDanat was assessed from the same data. The results clearly demonstrate that, in this preparation, increasing VT over the normal physiological range largely eliminates the dominant intraregional component of ventilation heterogeneity found at resting levels of VT but has little or no influence on interregional heterogeneity and VDanat.
| |
METHODS |
|---|
|
|
|---|
Apparatus and General Procedures
All animal procedures were reviewed and approved by the University of Washington's Animal Care Committee.Anesthesia was induced in five mongrel dogs (20-30 kg) of either sex by bolus injection of thiopental sodium (20 mg/kg iv) followed, as needed, by bolus injections of pentobarbital sodium (30-120 mg iv) to maintain a deep surgical plane. After induction, the animals were intubated and ventilated by piston pump ventilator (model 608, Harvard Apparatus, South Natick, MA) at a fixed VT of 15 ml/kg with zero end-expiratory pressure (EEP). Ventilatory rate was varied to maintain normal end-expired CO2, arterial O2, and arterial pH status, and an occasional "sigh" was given by double inflation. Concentrations of inspired and expired gases were continuously monitored with a respiratory mass spectrometer (model 1100, Perkin-Elmer, Pomona, CA), sampling gas midstream in a connector between the endotracheal tube and ventilator at a rate of 1 ml/s. The animals were anticoagulated with ~6,000 units of sodium heparin and exsanguinated, and the thorax was then accessed via sternotomy. EEP was increased to +5 cmH2O to maintain the lung near normal end-expiratory volume, except when partial deflation was needed during the excision procedure to avoid pleural damage. The right-middle lobar bronchus was exposed by dissecting away surrounding blood vessels and connective tissue, carefully sparing the lobar pleura. The lobar bronchus was cross-clamped near the hilum at the end of a tidal inspiration, transected central to the cross-clamp, and mounted on a semirigid plastic tube (Tygon). The lung surface was occasionally misted with normal saline throughout the excision procedure.
The lobe was then suspended in a sealed plethysmograph (constructed for
these experiments) that allowed air from outside of the box to flow
into and out of the lobe through the mounting tube (Fig.
1). Rigid, closed-cell foam was used to
fill much of the chamber volume not occupied by the lobe. The lobe was
inflated by a negative box pressure (Pbox) of
5
cmH2O [i.e., to a transpulmonary pressure (Ptp) of +5
cmH2O] by a small vacuum pump adjusted to maintain a
stable EEP. A computer-controlled linear-motor ventilator (custom
designed and locally built) was used to remove and replace a volume of
the box gas to produce the desired swing of Ptp and ventilate the lobe.
This pressure-controlled system allowed VT, TI,
TE, and breath-hold time (after inspiration and/or
expiration) to be varied independently. To achieve constant
instantaneous gas flow rates during inspiration and expiration given
the nonlinear compliance characteristics of the system (box, gas,
lobe), a custom software program was developed and used to modify
instantaneous computer voltage output controlling the ventilator.
Negative Pbox was measured by a differential pressure transducer
(±2.5 psi, Statham Instruments, Hato Rey, PR) and delivered to a
chart recorder (DMS 1000, Western Graphtech, Irvine, CA). The signal
was also fed to a separate computer for digital recording in a data
file and for visual monitoring of pressure levels on a screen display during the experiments. The pressure transducer voltage signal was
calibrated over the range of
55 cmH2O to +5
cmH2O.
|
Initially, air was removed from the box with a large syringe to produce
several Ptp swings of 30 cmH2O (Pbox reduced from
5
to
35 cmH2O) to apply a standard volume history. A
resting-level tidal breath was defined as a Ptp swing of 12 cmH2O (Pbox from
5 to
17 cmH2O).
TI and TE were set at 2 s and 4 s,
respectively. The mass spectrometer sample inlet was positioned to
access the center of the bronchial tube where it exited the box. A
Krogh pediatric spirometer (locally built) was attached to this
mounting tube when measurement of VT was needed. The
spirometer was then removed and replaced on the mounting tube by a
pediatric T-piece two-way valve (2230 series, Hans Rudolph, Kansas
City, MO) arranged to deliver the marker and washout gases. Lobes were
ventilated at 10 breaths/min with a gas mixture of 1% He-1%
SF6 (balance air). Signals for concentrations of He,
N2, SF6, and Pbox were stored in a computer
file as digital values sampled at 36 Hz for each data-acquisition channel.
Initial gas concentrations were considered to be uniform throughout the lobe when concentrations of gases at the sampling site were invariant over the entire ventilatory cycle. After several minutes of ventilation with the chosen pattern, the inspired gas mixture was changed to 100% O2, and MBW of indicator gases was continued and recorded for at least 30 breaths. The lobar surface was maintained in moist condition by injecting ~50 ml of room-temperature normal saline through a misting system surrounding the lobe between each MBW.
Six end-inspiratory pressures (EIP) were chosen in a random manner to
produce tidal Ptp swings from 5 to 35 cmH2O in
5-cmH2O increments. Apparatus dead space volume
(VDapp) was ~6 ml, TI and
TE remained at 2 and 4 s, respectively, and EEP remained at
5 cmH2O. After each washout, VT was
returned to the resting level, a double inflation "sigh" was
given, the next EIP was selected, VT was measured, and
another MBW was performed. The final (seventh) MBW was carried out
under the same conditions as the first to assess reproducibility and
stability of measurements over the course of the experimental procedures.
Volume at end-expiration [lobe volume (VL)] for each lobe
was estimated only from N2 MBW. With each breath of
O2, the alveolar N2 is diluted by the fraction
VL/(VL + VA) = 1/(1 + x),
where alveolar volume (VA) equals VT
total dead space volume (VD), and x equals
VA/VL. Thus after n
breaths the ratio of the final phase III N2 concentration
(FnN2) to the initial value (F0N2) is
|
(1) |
n log(1 + x). The slope of a graph of
log(Fn/F0) vs. n has slope =
log(1 + x). Hence, x and then VL can
be derived from washout data. In the presence of substantial
heterogeneity of ventilation, this method is weighted toward the
compartments with early washouts (i.e., less than the total
VL).
Measurement of Exhaled Gas Dead Space Volume and Heterogeneity
VDanat.
VD was quantified from the N2 concentration
profile of the second washout breath (Fig.
2A) by a method similar to that
described by Young (36) but replacing CO2 by
N2, in which volume is assessed to the point at which
exhaled indicator gas concentration reaches 50% of the alveolar level
(defined at 0.75 VT). VDanat was
then obtained by subtracting VDapp from
VD (VDanat = VD
VDapp).
|
Early and late phase III slopes. MBW results in a series of diminishing exhaled marker gas concentration profiles that each display the three phases classically labeled as I, II, and III (Fig. 2B). Slopes of two portions of the phase III plateau were determined from each washout breath as indicators of alveolar gas concentration uniformity (Fig. 2A). Early emptying units were represented by marker gas measured from 0.5 to 0.75 of the exhaled VT (phase IIIE), and late emptying units were represented by 0.75 to 0.95 of VT (phase IIIL). Slopes were determined as a function of volume by best fitting a linear equation to data points obtained from phase IIIE and phase IIIL.
Normalization and further analysis of phase III slopes. To discern contributions of only cdi and dcdi to nonuniformity of exhaled gas, a normalization procedure was employed in a manner similar to that described by Crawford et. al. (8), with several modifications (see Ref. 11 and APPENDIX A). To remove the effect of progressive gas dilution on slopes, the normalized slopes of phase IIIE (SnIIIE) and phase IIIL (SnIIIL) were found by dividing best fit slopes by the average concentration of marker gas in the respective intervals (i.e., E or L).
Measured values of SnIII (SnIIIE observed and SnIIIL observed) from each washout breath were plotted as a function of breath number (Fig. 2C). Contributions of convective-dependent inhomogeneity (SnIII cdi) and diffusive-convective-dependent inhomogeneity (SnIII dcdi) to the total calculated inhomogeneity of each washout breath (SnIII total) were estimated by a curve-fitting procedure (Fig. 2C). For each breath of a MBW, cdi is assumed to add a slightly diminishing, but nearly continuous, additional interregional inhomogeneity of marker gas dilution to that present in each previous breath (23). Therefore, SnIII cdi as a function of n is approximated by an equation of the form
|
(2) |
|
(3) |
|
(4) |
Statistical Analysis
Paired t-tests and the Wilcoxon signed-rank test were used for comparison of the following: 1) VT, VL, and SnIIIL total measurements recorded under identical conditions of EIP at the beginning and end of each experiment, 2) early (E) and late (L) values of SnIII 1, Snf, and cdimax, obtained under the smallest vs. largest conditions of tidal volume (EIP =
10 cmH2O vs.
35
cmH2O), 3) simultaneously obtained E vs. L values
of SnIII 1, Snf, and
cdimax obtained with EIP =
10 cmH2O and
35 cmH2O, and 4) the difference between
simultaneously obtained E and L values of
SnIII 1, normalized by the value of that
difference found for EIP =
10 cmH2O, for EIP =
10 cmH2O vs.
35
cmH2O.
The effect of altering EIP on VDanat, SnIII total, SnIII dcdi, and SnIII cdi was described by the slope of the best fit linear and single exponential regression equations, but the results are reported for only the best fit solutions (linear or single-exponential; see RESULTS). For this analysis, the six individual values of SnIII total, SnIII dcdi, and SnIII cdi per lobe were normalized by the average of those measurements for that lobe to control for average differences between lobes. The statistical significance of any EIP influence on VDanat and SnIII results (normalized) was also tested by paired t-test comparison of best fit slope values to a zero slope predicted by the null hypothesis.
Except where indicated, all results are expressed as means ± SE in the text and Figs. 4-6.
| |
RESULTS |
|---|
|
|
|---|
General Lobar Measurements
None of the lobes exhibited visible increases in gas trapping at the end of the experimental procedures. The average amount of time from the first to the last MBW was 1.7 ± 0.5 h (mean ± SD).Table 1 lists the values of VT,
VL, and N2
SnIIIL total obtained for each lobe from the
first and last MBW runs performed under equal conditions of
EIP. Compared with the initial washout, the same EIP used
for the final washout produced an ~4% increase in VT,
which was not found to be statistically significant. An ~18% greater
VL measured in the final vs. initial washout was found to
be significant by the Wilcoxon signed-rank test (P < 0.05)
but not the paired t-test. N2
SnIIIL total measured from the final washout
averaged 0.76 ± 0.21 (mean ± SD) of that obtained from the initial
washout, but this potential trend of decreasing N2
SnIIIL total over the course of the procedures
was not found to be statistically significant.
|
Exhaled Dead Space Volume
VDanat was not markedly altered by increasing EIP as VT increased by approximately twofold (Fig. 3). Regressions of VDanat were calculated against both EIP and ln(EIP) (best correlations found against EIP, R2 range; 0.38 to 0.66). Slopes of the best fit regression equations ranged from
0.07 to +0.07 (ml/cmH2O), not
significantly different from zero for the group (P = 0.66).
|
Heterogeneity of Exhaled N2
Measured values of SnIII.
SnIIIE 1 was greater than
SnIIIL 1 when measured with EIP =
35
cmH2O and
10 cmH2O, in every instance
(average values, Fig. 4).
SnIIIE 1 was found to be significantly greater
than SnIIIL 1 for EIP =
35
cmH2O by the two-tailed paired t-test
(P < 0.04) and Wilcoxon signed-rank test (P < 0.05), but a significant difference was only found with EIP =
10 cmH2O by the Wilcoxon signed-rank test
(P < 0.05). When normalized by the difference
between SnIIIE 1 and
SnIIIL 1 found with EIP =
10
cmH2O, the difference value was found to be significantly
reduced with EIP increased to
35 cmH2O, when
analyzed by the paired t-test (P < 0.04) but not the Wilcoxon signed-rank test.
|
35
cmH2O vs.
10 cmH2O, as tested by the
paired t-test, was precluded by the wide range of
SnIII 1 values found for the group with EIP =
10 cmH2O. However, the consistent reductions in
SnIIIE 1 and SnIIIL 1
with increased VT were found to be significant by the
Wilcoxon signed-rank test (P < 0.05).
When compared with L Snf, simultaneously measured
E Snf values were found to be significantly
larger for EIP =
10 cmH2O by the Wilcoxon
signed-rank test (P < 0.05) but not for EIP =
35 cmH2O (Fig. 5A).
Snf was not found to be significantly larger for
the E vs. L portions of the phase III slopes with EIP =
10 cmH2O or
35 cmH2O, when tested by the
paired t-test. Snf for the E and L
portions of the phase III slope were not found to be significantly
different when measured with EIP =
10 cmH2O vs.
35 cmH2O, as tested by the paired t-test and
Wilcoxon signed-rank test.
|
Curve-fit values of SnIII and comparison with measured
values.
In regard to goodness of fit, curve-fit results for N2
SnIIIE observed vs. breath number were very
similar to results found for N2
SnIIIL observed vs. breath number. In the
interest of brevity, parameters that produced the best fit solutions of
Eq. 4 to measured values of N2
SnIII observed are only presented for the more
commonly studied L portion of the phase III slope (Table
2). Observed SnIIIE 1 and SnIIIL 1
values were highly correlated with the corresponding best fit
SnIIIE total and
SnIIIL total values, respectively
(R2 = 1.0), and weighted correlation of the best
fit solutions to each series of weighted
SnIIIL observed measurements from each MBW was
excellent in every instance (R2
0.98). This high
degree of correlation reflects the emphasis on measured values from
earlier vs. later washout breaths imposed by the weighting procedure
that strongly influenced the best fit solutions. Correlations between
curve-fit values and the corresponding unweighted values of
SnIIIE observed and
SnIIIL observed were also very good in most
instances (R2 = 0.78 ± 0.04 and 0.84 ± 0.03, respectively; individual correlations given in Table 2 for
SnIIIL observed).
|
10
cmH2O, but not with EIP =
35 cmH2O (Fig.
5B). The E portion of cdimax was significantly
larger when measured with EIP =
10 cmH2O vs.
35 cmH2O (Wilcoxon signed-rank test, P < 0.05). L values of cdimax were not found to be different
when measured with EIP =
10 cmH2O vs.
35
cmH2O, as tested by the Wilcoxon signed-rank test.
Significant differences were not found when these comparisons of
cdimax were made by the paired t-test.
As suggested by the grouped mean best fit values from the first washout
breath (Figs. 6),
SnIIIE total was greater than simultaneously measured SnIIIL total in every
instance (individual data not shown, VDapp/VT = 0). Increasing
EIP resulted in a large decline in N2 SnIIIE total and
SnIIIL total primarily because of the dcdi
component. The dominant SnIII dcdi component
of SnIII total varied widely among the five
lobes at the lowest levels of EIP.
|
0.064 ± 0.023 l/[ln(cmH2O)] (R2 = 0.72 ± 0.08). For normalized N2
SnIIIL total vs. ln(EIP), slopes averaged
0.047 ± 0.013 l/[ln(cmH2O)]
(R2 = 0.78 ± 0.04). For the group of five lobes,
slopes were significantly greater than zero as tested by the paired
t-test (P < 0.048 and P < 0.024 for E and L values, respectively).
N2 SnIIIE dcdi and
SnIIIL dcdi results mirrored those of
SnIIIE total and
SnIIIL total. Fig.
7 contains the within-lobe normalized values of SnIIIE dcdi and
SnIIIL dcdi vs. EIP for the five lobes tested
(VDapp/VT = 0), with lines
representing the best fit solutions. For normalized N2
SnIIIE dcdi vs. ln(EIP), slopes averaged
0.068 ± 0.022 l/ [ln(cmH2O)]
(R2 = 0.67 ± 0.13, however, for one lobe
R2 = 0.16 and for the other four lobes
R2 = 0.80 ± 0.06). For normalized N2
SnIIIL dcdi vs. EIP, slopes averaged
0.054 ± 0.013 l/[ln(cmH2O)]
(R2 = 0.79 ± 0.04). For the group of five lobes,
slopes were significantly greater than zero as tested by the paired
t-test (P < 0.039 and P < 0.014 for E and L values, respectively).
|
0.038 ± 0.017 l/[ln(cmH2O)] (R2 = 0.44 ± 0.18). Slopes of normalized SnIIIL cdi vs.
ln(EIP) averaged
0.021 ± 0.011 l/[ln(cmH2O)] (R2 = 0.27 ± 0.14). Slopes of SnIIIE cdi and
SnIIIL cdi vs. ln(EIP) were not found to
be significantly different from zero by the paired t-test
(P = 0.084 and P = 0.442, respectively).
| |
DISCUSSION |
|---|
|
|
|---|
Overview
The present study was designed to improve understanding of how intralobar VDanat and heterogeneity of gas distribution are influenced by altering VT. Increasing VT by varying EIP was found to have little influence on VDanat (Fig. 3), but significantly decreased E and L exhaled SnIII (Figs. 4, 6, and 7), primarily because of decreasing dcdi (Figs. 6 and 7). Maintaining EEP fixed and controlling EIP, as opposed to predefining levels of VT, facilitates interpretation of results from lobes of varied size. Assuming that compliance characteristics of lobes are similar, equal conditions of Ptp will favor 1) the same degree of stress applied to all lobes regardless of size, 2) equal specific ventilation, and 3) diffusion-convection fronts in similarly sized airways. Ventilatory patterns were delivered to the lobes with a very high degree of repeatability and accuracy. The effect of time on the preparation was to increase parenchymal compliance, as evidenced by generally increased VT and VL when measured under equal conditions of EIP for the initial and final washouts. Although there was an apparent temporal trend in four out of five lobes toward decreased SnIIIL total when the conditions of the first MBW were repeated for the last, this was not found to be statistically significant. The randomized nature in which EIP test conditions were applied also largely removes concerns that trends in the results were influenced by time. However, it remains a possibility that variability in the measurements was increased by aging of the preparation and thus potentially reduced or obscured minor trends in the data.Tidal Volume Influences on Exhaled Dead Space Volume
The absence of a measured change in exhaled VDanat with increasing VT is an important finding of this study and is at odds with the results and conclusions of most other investigations of VDanat vs. VT (9, 17-19, 31). Increasing VT is usually found to affect VDanat regardless of the method of measurement and has generally been explained by expansion of proximal airways (1) and various mechanisms that influence gas mixing within peripheral or proximal airways (9, 18, 25). VDanat derived from clearance curves increased with increasing VT in two human studies (9, 25), in agreement with predictions of mathematical models representing intra-airway and cardiogenic gas mixing mechanisms that were sensitive to peripheral airway geometry (25). VDanat was also found to increase with VT when estimated from human exhaled gas profiles by use of methods similar to this study (18), in agreement with model predictions that included gas dispersion in relatively proximal airway regions (19). However, the same modeling study did not find the VT influence on VDanat to depend on distal airway geometry (19). At least one study found no change in Fowler VDanat when VT was varied in humans, but no explanation was offered for the discrepancy with other studies (26). The present results from excised canine lobes are sufficient to allow hypotheses of intra-airway gas-mixing mechanisms to be discounted when addressing the influence of variable VT on VDanat. Canine and human lungs appear to be similar in regard to heterogeneity of distal airway branching geometry (E. R. Weibel, personal communication), and properties of more proximal airways that are predicted to influence gas dispersion are probably also similar.It is unlikely that this result is due to errors in the procedures of the experiment, in which the average SD of VDanat was ±0.6 ml (allowing for a resolution of about 1 ml). The second washout breath was used for measurement of VDanat to avoid variability in first-breath values due to manual valving. VDanat is known to vary with progressive washout breaths to a small degree (8), but this effect is not expected to alter our overall conclusions. The method used to estimate VDanat in these studies (36) is similar to the method proposed by Fowler that has often been used to determine this value (14), and differences are not expected to affect the overall conclusions of this study.
Along with every rise in VT, inspiratory flow rate also increased proportionally in this study because of fixed TI. Increasing inspiratory flow rate is predicted to move the diffusion-convection front position toward the lung periphery (12, 23). VDanat represents the summed average appearance on expiration of gas from the vicinity of the diffusion-convection front of all parallel airway regions (13), and therefore increasing inspiratory flow rate is predicted to increase VDanat (12, 22, 23). However, increasing airway cross-sectional area during inspiration tends to decrease convective velocity and to move the front position proximally (22). In this study, the average volume from the gas-sampling site in the lobar bronchus to the diffusion-convection front position appears to have remained unaltered. This finding implies that the average front position must have moved mouthward because increasing VT increases the volume of the airways. Therefore, under the conditions of these experiments, the increase in airway cross-sectional area was a greater influence on the average position of the diffusion-convection front than was the increase in convective flow rate. In partial support of this result, VDanat was found to remain unchanged when inspiratory flow rate was varied over a wide range in intact dogs mechanically ventilated with a constant VT (2). In that study, molecular diffusivity was found to be a greater determinant of the diffusion-convection front position than was convection.
The lack of VT influence on VDanat in this study of excised, unperfused lobes leads to the conclusion that, when such an influence is found in more intact preparations and humans, it is due to extralobar mechanisms and/or perfusion. Cardiogenic mixing and the nonuniform intrapleural pressure gradient are two extralobar mechanisms that are known to influence ventilation distribution (13), but it is not clear how those influences would produce the VT effect on VDanat. A need for the presence of these two extralobar mechanisms to produce the reported correlation between VT and VDanat is contradicted by the results of one experiment in which an excised, perfused lung preparation was used to demonstrate the VT influence on dead space ventilation (31).
One explanation for the discrepancy between our present results and those of previous studies where VDanat correlated with VT may be found in the known effects of gas exchange on alveolar gas concentrations. As pointed out by Cormier and Bélanger (4), the exhaled N2 profile during a washout with 100% O2 in perfused lungs is determined by several factors involving gas exchange. Alveolar N2 is initially diluted slightly in the beginning of O2 inspiration because of increased alveolar CO2 transfer, a fairly constant O2 consumption, and the low rate of N2 transfer between the liquid and gas phases. For a given CO2 production and a given TI and TE, the time required for the respiratory quotient to become <1.0 during a breath is longer for a large- compared with a small-volume inspiration, presumably resulting in a larger dilution of early-exhaled N2. Further investigation is required to address whether this mechanism contributes to the increase in VDanat found when VT is increased in perfused lungs.
Uniformity of Exhaled Alveolar Gas and Influences of Tidal Volume
Contributions of dcdi and cdi to overall heterogeneity. SnIII dcdi was a far greater influence on the uniformity of exhaled alveolar gas than was SnIII cdi in this preparation. The dcdi component was responsible for ~85% of N2 SnIIIL total when values for the lowest two conditions of EIP and VT were averaged and for ~70% of the total at the two highest levels of EIP and VT (corrected to VDapp/VT = 0, Fig. 6B). Human MBW studies have found the dcdi component of SnIII to be the major contributor to the phase III slope (6, 8). In a previous study of upright and awake humans breathing 1-liter volumes at a moderate rate, the dcdi contribution to the total N2 SnIII (measured for late-exhaled gas) was found to be ~62% after gas exchange was taken into account (8). The relatively greater contribution of dcdi in this excised lobar preparation was expected because of the absence of most external influences on ventilation distribution.
Evidence for diffusive-convective-dependent gas mixing during exhalation. A decrease in expired gas heterogeneity as exhalation proceeds (curvilinearity of the phase III slope) has been predicted by computer models of intrapulmonary ventilation distribution (24, 25) and is considered to be evidence for continued homogenization of pulmonary gas due to diffusive-convective interaction (6, 25). In agreement with those predictions and a previous demonstration in humans (6), we found that values representing early-exhaled gas (SnIIIE 1, SnIIIE total, and SnIIIE dcdi) were greater than corresponding values representing late-exhaled gas (SnIIIL 1, SnIIIL total, and SnIIIL dcdi) in every instance. Results from multiple 1.5-liter breath washouts in humans analyzed for the slope of phase III from 0.75 to 1.0 liters and from 1.0 to 1.5 liters suggested this result, when slopes from the later expirate were found to be consistently less than from the earlier expirate (6). A single-breath, inspiratory-capacity washout study in humans failed to detect a difference in heterogeneity of early- vs. late-exhaled alveolar gas (24); however, the disparity may be due to added interregional heterogeneities imposed by vital capacity maneuvers and/or the decrease in intraregional heterogeneity expected from large volume breathing. The ability to resolve this decreasing heterogeneity of exhaled gas for both small and large tidal breaths in these experiments was probably enhanced by the lack of change in VDanat (an effect that would influence early-exhaled alveolar gas), and because heterogeneity of exhaled gas was not influenced by gas exchange (an effect that increases heterogeneity of later exhaled alveolar gas; Ref. 4).
VT-dependent alteration of total heterogeneity and the
dcdi and cdi components.
The measured values of SnIII 1 were found to
be decreased when tested with the largest vs. the smallest level of
VT in every instance, for both early- and late-exhaled
phase III gas. However, the wide variability in this index among the
lobes with EIP =
10 cmH2O limited the statistical
significance of this result (Fig. 4). Variation in
SnIII 1 among the various lobes was greatly reduced with EIP =
35 cmH2O. Best fit values of
SnIII total were essentially identical to the
measured SnIII 1 values and also displayed
wide variability among the lobes (Fig. 6). Most human studies have
found heterogeneity of exhaled alveolar gas concentrations to decrease
with increased VT, beginning with early investigations of
exhaled gas profiles (15, 27). More recently, single-breath washouts in
humans revealed phase III slopes that were larger for 1-liter
inspirations than for an inspiratory capacity maneuver (24), in general
agreement with human N2 clearance studies (3, 28).
Hyperventilation (increased VT and frequency) decreased N2 clearance in one human study (38), but comparison with
the present results is limited because voluntary hyperventilation is
known to alter chest wall deformation patterns and interregional sequencing in addition to gas flow rate (13). On the whole, increasing
VT has been found to result in more uniform mixing of tidal
and resident gas as measured by N2 clearance or exhaled gas uniformity.
Potential mechanisms for VT-dependent effects on dcdi and cdi. Decreased intraregional heterogeneity with increased VT may be the result of mouthward movement of the diffusion-convection front, decreased resistance to diffusive mixing, and/or increased uniformity of acinar inflation. Heterogeneity of airway branching geometry and nonuniform airway and parenchymal mechanical characteristics are inherent properties of the normal human lung structure that increase among progressively smaller regions and are generally thought to influence heterogeneity of intraregional ventilation (20, 23, 29, 33, 35). As previously discussed, the VDanat results argue for a mouthward movement of the average diffusion-convection front position as VT was increased. This shift in the front position is predicted to allow increased diffusive mixing of parallel airway regions and increase the uniformity of exhaled alveolar gas concentrations (22, 23). Resistance to diffusive gas transfer is expected to be nonuniform among parallel airways with nonuniform geometry due to differences in path length for molecular travel (24), and this nonuniformity is expected to decrease with increased tidal volume due to increased airway dimensions (6).
Increasing EIP in this preparation increases stretching force uniformly throughout the lung, but heterogeneous airway and parenchymal compliance characteristics must produce nonuniform expansion, perhaps resulting in more uniform ventilation with greater VT. Modeling studies that utilize only documented nonuniformity of peripheral airway geometry generally underestimate the magnitude of exhaled gas nonuniformity and its VT dependence in humans (24, 33). Sufficient modification of model geometry can bring results of simulations closer to those actually observed experimentally (33), but the lung structure clearly possesses heterogeneous mechanical properties that should also be taken into account when interpreting the results of these experiments (29, 35). Heterogeneity of intraregional compliance may be able to account for discrepancies between model predictions and experimental results, but the complexity of modeling mechanical heterogeneities in the fine airways still precludes their use in detailed computer simulations of intrapulmonary gas mixing. The configuration of the excised lobe preparation may have influenced the cdi results to a small degree, in particular the gravity-determined nonuniform tissue stretch that is expected due to hanging lobes by their bronchi. The weight of the blood-free lobe is supported by tissue adjacent to the bronchus, and this local stress progressively diminishes peripherally. Distortion-dependent changes in ventilation distribution are expected to be manifest primarily as altered cdi, reflecting regional differences in lung expansion between top and bottom. The small decrease in cdi found for early-exhaled alveolar gas when VT was increased is consistent with the presence of this potential distortion effect. Increasing ventilatory flow rate has been shown to have mixed influences on the uniformity of exhaled indicator gas (15, 16, 25, 37, 38). Alteration of interregional sequencing is probably responsible for much of the flow rate effect on ventilation heterogeneity in spontaneously breathing humans (37), but peripheral airway effects of flow rate have also been supported by human investigations (25). In contrast to increasing heterogeneity found when inspiratory and expiratory flows were voluntarily increased during vital capacity maneuvers (25), these experiments demonstrate the opposite effect in response to increased VT and flow rate. Furthermore, in related experiments using this preparation, we found that altering inspiratory time from
to 8 s (fixed VT) failed
to significantly alter SnIII dcdi or SnIII total (11). Results from the present
study argue against a role for the unperfused lobar structure in
producing increased ventilation heterogeneity often associated with
increasing convective flow rate (over the range tested).
Unlike healthy young humans, dogs possess extensive collateral
ventilation that has been implicated in homogenizing gas distribution between pulmonary regions (10) and may account for some of the improved
uniformity of exhaled gas when VT is increased. However, several factors decrease the possible importance of this mechanism in
producing the VT effect. Communications in the walls of
acinar and alveolar regions, as represented by pores of Kohn, are
expected to enlarge with increased lung stretch, but the normal liquid lining probably reduces the free communication of gases. Also, the
N2 diffusion-convection front position is predicted to be near the entrance of acinar regions, resulting in uniform mixing within
more distal regions, regardless of the path (22). Interregional communications exist at higher levels of the canine bronchial tree
(10), and these pathways may increase diffusive- and/or convective-dependent mixing of gases with increased VT.
Although interregional communications may have contributed to these
results to some degree, the generally similar improvement of gas mixing in human and canine lungs with increased VT argues against
the necessity of such communications for the VT effect.
Conclusions
These studies of the effects of VT on ventilation heterogeneity and VDanat in excised canine lobes have demonstrated that larger VT significantly decreases intraregional heterogeneity of ventilation within and among acinar regions. In contrast, VDanat was found to be unaltered by changes in VT, and also interregional heterogeneity of ventilation that occurs among more proximally branching regions was essentially unchanged or slightly decreased. These results indicate that the magnitude of VT is an important consideration when interpreting results of other investigations in which uniformity of exhaled gas is a measured parameter. Furthermore, the ability to detect changes in SnIII is reduced with large VT maneuvers because of substantial VT effects on the dominant dcdi component of ventilation heterogeneity. Thus large-volume maneuvers cannot accurately represent the heterogeneity of ventilation present under conditions of resting VT. These results also suggest that VT-dependent influences on VDanat and interregional heterogeneity require the presence of normal extralobar influences on the lung parenchyma and/or perfusion.| |
APPENDIX A |
|---|
|
|
|---|
Modified SnIII Analysis
Procedures used in this study for determination of SnIII and its cdi and dcdi components contain several modifications of the method previously introduced by Crawford et al. (8), which had evolved from work by Paiva (21) and Paiva and Engel (22). Normalizing the phase III slope by the mean gas concentration of the specific volume of interest (e.g., E or L in this study, Fig. 2A) is a variation of the original method by which the mean gas concentration of the entire exhaled volume was utilized (8, 11) and essentially removes the assumption that early- and late-exhaled marker gases are washed out at similar rates. The present procedure for determining cdi and dcdi components of SnIII total is the first to utilize a biexponential equation best fit to measured slope values from MBW. In the original method, only the cdi process is represented by a best fit function, and here we add a second single-exponential equation representing the dcdi process in the complete equation. The form of this second equation is appropriate when considering previous descriptions and simulations of the dcdi contribution to SnIII during MBW (7, 8, 22).The addition of a second equation in this version of SnIII analysis, as well as the weighting of data in the curve-fit procedure, imposes certain restrictions on the possible solutions. In addition to requiring a solution with more components, the 1/n weighting procedure gives earlier washout breaths more influence over the best fit solution of Eq. 4, thus reducing the influence of increased error in values measured from later washout breaths due to falling marker gas concentrations (Fig. 2C). This SnIII analysis procedure probably gives slightly smaller values of SnIII cdi and greater values of SnIII dcdi than would have been found previously. Model analysis predicts that the dcdi process does not actually reach dynamic equilibration before the fifth breath (7, 8), and the exact number of breaths required to reach dynamic equilibrium of dcdi in lungs of humans and dogs is not known. Left unaccounted for, a residual influence of dcdi on the rate of increase in SnIII vs. n beyond breath 4 likely leads to an overestimated rate of increasing SnIII cdi.
By Eq. 4, the best fit solution is required to pass through the origin. This requirement essentially formalizes previously stated assumptions that all influences on the normalized slope values are generated by dcdi and/or cdi (excluding gas exchange; Ref. 4) and that marker gas concentrations are uniformly distributed in the lung before washout (8, 21, 34). Although some objections may be raised concerning the effect of this restriction on the possible best fit solutions, this method conveys some theoretical and practical advantages over earlier methods (8, 34). By the original method, a single-exponential fit to slope values obtained after breath 4 represents the cdi process, and the solution was back-extrapolated to the y-axis (in which the model function is zero) (8). SnIII cdi for breath 1 was the difference between the back-extrapolated values of SnIII from breath 1 and the y-intercept, and SnIII dcdi was the difference between the first-breath measured SnIII and the back-extrapolated value of SnIII cdi (8). By the original method and recent variations, indexes representing ventilation heterogeneity in peripheral airways are subject to all error present in single measurement of alveolar slope from the first washout breath (8, 34). We have found that, in many experimental situations, slope values from the first washout breath are often prone to errors for a variety of reasons. An inherent advantage of this fitting procedure lies in the use of all data points from the washout to determine both peripheral and proximal airway contributions to ventilation heterogeneity, reducing the influence of single-measurement errors. In addition, required passage of solutions through the origin prevents generation of solutions with negative y-intercepts. Such solutions lead to either subjective judgment that they are not appropriate or difficult theoretical explanation for potentially negative values of SnIII cdi and/or SnIII dcdi. Earlier detailed comparison of this, earlier, and other SnIII methodologies found that solutions for most sample MBW data sets were very similar (11). However, methods that employed the 1/n weighting procedure and/or required passage through the origin were superior for allowing computer-generated solutions to be found for difficult data sets (i.e., data sets containing outlying values or wide variations in measured breath-to-breath values of SnIII). In our experience, the present methodology allows nonsubjective SnIII analysis to be applied more successfully in situations in which such variations normally occur (e.g., spontaneously breathing humans, manual valving, low marker gas concentrations, etc.).
Quantitative comparison of the dcdi vs. cdi contribution to total ventilation heterogeneity is maximal for values determined for the first washout breath, because starting conditions for both processes are equal only before MBW. The dcdi dominates first-breath SnIII values (8, 12, 32), and later washout breaths contain increased information regarding cdi (23). The components of SnIII have been quantified in other ways that allow for more detailed examination of the cdi process (5, 7). The change in SnIII beyond five or six normal tidal breaths (or approximately >1.5 lung volume turnovers) is predicted to be almost entirely determined by the cdi process (8, 34), and SnIII of some or all of those breaths can be averaged as other indications of cdi influence on SnIII (5-7). SnIII cdi is usually a minor contributor to SnIII total for the first washout breath, but small changes in cdi measured for washout breath 1 can produce large changes in SnIII of later washout breaths because of compounded influence of cdi. We have found that SnIII cdi and Snf are only adequate for detecting relatively large changes in cdi. The values of SnIII cdi were usually a small fraction of the already small values of SnIII total, resulting in a small signal-to-noise ratio. Similarly, SnIII observed from later washout breaths normally exhibits wide variability due to increasing measurement error of decreasing marker gas concentrations and phase III slope values, thus limiting resolution of Snf. Therefore, we quantified cdimax as the value of SnIII cdi determined by curve fit for breath 25 of the MBW. In theory, cdimax contains the same cdi information found in Snf but excludes the influence of dcdi and is less prone to single measurement errors.
| |
APPENDIX B |
|---|
|
|
|---|
Influences of Apparatus Dead Space on SnIII
Even small changes in the ratio of VDapp to VT can have large influences on SnIII. Aside from any true alteration in the distribution of inspired gas, an alteration in VDapp and/or VT can introduce variability in VDapp/VT and SnIII. To control for the influence of VDapp, a subset of He, N2, and SF6 MBWs was performed under varied conditions of VDapp. SnIIIL results were then evaluated to create correction factors that were applied to other SnIII results from these studies.Methods
Right middle lobes were removed from four mongrel dogs and prepared for negative-pressure ventilation as described in the text. The data acquisition system sampled each channel at 24 Hz. The minimum VDapp (~10 ml) consisted of the mounting tube and pediatric one-way valve. VDapp was then varied in a random manner from 10 to 45 ml in 5-ml increments by altering the length of tubing between the gas sampling site and the gas source/exit. All other parameters of ventilation were kept constant.