Blood flow distribution within the airway wall.

Altered perfusion of the bronchial mucosal plexus relative to the adventitial plexus may contribute to geometric changes in the airway wall and lumen. We studied bronchial perfusion distribution in sheep by using fluorescent microspheres at baseline and during intrabronchial artery challenge with methacholine chloride (MCh; n = 7). Additionally, we measured airway resistance (Raw) during MCh with control or increased perfusion (n = 9). Raw with MCh was significantly greater for high than control flow. Microspheres in histological sections lodged predominantly in the mucosa (60%), and this was not altered by MCh. However, more microspheres lodged in airways >1-mm in diameter during MCh and increased perfusion than MCh and control flow. In airways < or =1 mm in diameter, fewer microspheres lodged during control than increased flow. If the number of microspheres represents regional agonist access to airway smooth muscle, then the differences observed in Raw can be explained by the distribution of agonist. During challenge, there was greater MCh delivery to larger airways during increased flow and less delivery to smaller airways during control flow. The results demonstrate the effects of axial perfusion distribution on Raw.     

Radiographic visualization of airway wall movement during oscillatory flow in dogs

It has been suggested that radial movement of the central airway walls during oscillatory flow might contribute to the increased frequency dependence of compliance seen in chronic obstructive pulmonary disease (COPD) (J. Appl. Physiol. 26: 670-677, 1969). Radial airway wall motion has also been invoked to explain the frequency-dependent decreases in the efficiency of gas exchange during low-volume high-frequency ventilation (HFV) in histamine-bronchoconstricted dogs and in patients with respiratory insufficiency. To test the possibility that airway wall motion increases with bronchoconstriction, we measured central airway diameters using cinebronchoradiography in anesthetized tracheostomized dogs during oscillatory HFV [50 and 100 ml tidal volume (VT) at frequencies (f) of 2, 6, and 12 Hz], under control conditions, during electrical stimulation of the vagi, and after exposure to histamine aerosol. Cineradiobronchograms from two dogs were evaluated quantitatively for tracheal diameter and for lengths and diameters of a number of major airways. Under control conditions, the diameter of the airways fluctuated 7-9% of the mean with VT of 50 ml and 9-18% with VT of 100 ml in the range of frequencies studied. Bronchoconstriction produced by aerosolized histamine increased radial airway wall movement to 10-47% with VT of 50 ml, and during vagal stimulation diameters changed 7-20% at VT of 50 ml. After histamine, the central airways displayed large diameter changes during HFV, whereas more peripheral airways were markedly constricted and did not change in diameter.(ABSTRACT TRUNCATED AT 250 WORDS)     

Airway branching morphology of mature and immature rabbit lungs.

The scheme of Horsfield et al. for describing the pulmonary airway tree (J Appl Physiol 52: 21-26, 1982) catalogs each airway according to its order and the difference in order of its two daughters (denoted Delta). Although this scheme captures the natural asymmetry in the airway tree, it is still deterministic, because it assumes that all airways of a given order are the same; yet such variability is extremely important in determining the overall behavior of the lungs. We therefore analyzed complete lung lobes from three mature and two immature rabbits and determined the Horsfield order and Delta of every airway down to the terminal bronchioles. We also measured the diameter of each airway. This allowed us to determine the average structure of the rabbit airway tree, the variation about this average, and also how the structures of mature and immature airway trees compare. We found some variation in branching asymmetry and airway diameter at a given order between animals but no evidence of systematic differences in structure between mature and immature lungs. We found evidence of a difference in the branching structure of the peripheral vs. the central part of the airway tree (the break point being around order 20). We also determined the nature of the variation in Delta and diameter as a function of order, which should be valuable for the development of computer models seeking to encapsulate the naturally occurring regional variation in airway geometry in the normal rabbit lung.     

A comparison of the dose-response behavior of canine airways and parenchyma

We compared the histamine responsiveness of canine airways and parenchymal tissues in six anesthetized paralyzed open-chest mongrel dogs, partitioning total lung resistance (RL) into airway resistance (Raw) and tissue viscance (Vti). Pressure was measured during tidal breathing (frequency was 0.3 Hz) at the trachea and in three alveolar regions by use of alveolar capsules. Measurements were taken before and after the delivery of increasing concentrations of aerosolized histamine (0.1-30 mg/ml). We found that Vti accounted for 78 +/- 8% of RL under base-line conditions; this proportion remained relatively constant throughout the histamine concentration-response curve. There was a significant correlation between percent change in Vti and percent change in Raw at all levels of histamine-induced constriction (P less than 0.001). Moreover, the sensitivity of the tissues and airways (defined as the concentration of histamine required to double resistance) was remarkably similar. We conclude that, at this frequency of ventilation, Vti accounts for the major portion of RL both under base-line conditions and after histamine-induced constriction. Although increases in RL cannot be attributed solely to events occurring in the airways, the close correlation between changes in Raw and Vti and the similar sensitivities of the two support the use of indexes reflecting changes in airway caliber as an indicator of overall lung histamine responsiveness.     

Distribution of airway narrowing during hyperpnea-induced bronchoconstriction in guinea pigs

Increasing minute ventilation of dry gas shifts the principal burden of respiratory heat and water losses from more proximal airway to airways farther into the lung. If these local thermal transfers determine the local stimulus for bronchoconstriction, then increasing minute ventilation of dry gas might also extend the zone of airway narrowing farther into the lung during hyperpnea-induced bronchoconstriction (HIB). We tested this hypothesis by comparing tantalum bronchograms in tracheostomized guinea pigs before and during bronchoconstriction induced by dry gas hyperpnea, intravenous methacholine, and intravenous capsaicin. In eight animals subjected to 5 min of dry gas isocapnic hyperpnea [tidal volume (VT) = 2-5 ml, 150 breaths/min], there was little change in the diameter of the trachea or the main stem bronchi up to 0.75 cm past the main carina (zone 1). In contrast, bronchi from 0.75 to 1.50 cm past the main carina (zone 2) narrowed progressively at all minute ventilations greater than or equal to 300 ml/min (VT = 2 ml). More distal bronchi (1.50-3.10 cm past the main carina; zone 3) did not narrow significantly until minute ventilation was raised to 450 ml/min (VT = 3 ml). The estimated VT during hyperpnea needed to elicit a 50% reduction in airway diameter was significantly higher in zone 3 bronchi [4.3 +/- 0.8 (SD) ml] than in zone 2 bronchi (3.5 +/- 1.1 ml, P less than 0.012).(ABSTRACT TRUNCATED AT 250 WORDS)     

Airways impedance during single breaths of foreign gases.

The changes in airways resistance (Raw) and inertance (Iaw) during single inspirations of pure methane, helium, neon, and ethane at a flow of 0.1 l/s were measured in six healthy subjects by use of a forced-oscillation technique. Raw and Iaw were computed from respiratory transfer impedance obtained at a frequency of 20 Hz by applying pressure oscillations at the chest and measuring flow at the mouth with a bag-in-box system. Compared with the air data, the changes of Iaw after inhalation of 500 ml of gas averaged -41.1% with methane, -82.8% with helium, -25.8% with neon, and +4.8% with ethane. These changes were slightly less than the changes in gas density (-45%, -86%, -31%, and +5%, respectively). The inhaled volumes at which 50% of the changes had occurred (V50) did not differ significantly among gases and were approximately 100 ml. For Raw the data were more noisy than for Iaw; they were discarded in two subjects because of a strong and irreproducible volume dependence in air. Consistent differences were seen between the remaining subjects, one of whom exhibited a predominant viscosity dependence of Raw, one a predominant density dependence, and two an intermediate pattern. V50s were larger for Raw than for Iaw, indicating a more peripheral distribution of Raw. For Raw, V50s were lower with helium than with methane, in agreement with the notion that density-dependent resistance is located mainly in the large airways. The results suggest that some information on the serial distribution of Raw and Iaw may be derived from impedance measurements with foreign gases.     

Diameters, generations, and orders of branches in the bronchial tree

Studies of bronchial tree data by West et al. (J. Appl. Physiol. 60: 1089-1097, 1986) have shown that plots of mean diameter against generation, using log-log scales, can be represented by a power function with harmonic modulations. Other studies have shown that the mean diameter of the airways is exponentially related to order of branching. This paper demonstrates that both observations are compatible with a fractal model of branching, and because airway branching is fractal, this may explain why both are also true of the bronchial tree. Furthermore, the exponential relationship of mean diameter with generation in the larger airways, demonstrated by Weibel, is shown to result from the exponential relation of diameter with order in the fractal model.     

Numerical modelling and analysis of peripheral airway asymmetry and ventilation in the human adult lung

We present a new one-dimensional model of gas transport in the human adult lung. The model comprises asymmetrically branching airways, and heterogeneous interregional ventilation. Our model differs from previous models in that we consider the asymmetry in both the conducting and the acinar airways in detail. Another novelty of our model is that we use simple analytical relationships to produce physiologically realistic models of the conducting and acinar airway trees. With this new model, we investigate the effects of airway asymmetry and heterogeneous interregional ventilation on the phase III slope in multibreath washouts. The model predicts the experimental trend of the increase in the phase III slope with breath number in multibreath washout studies for nitrogen, SF(6) and helium. We confirm that asymmetrical branching in the acinus controls the magnitude of the first-breath phase III slope and find that heterogeneous interregional ventilation controls the way in which the slope changes with subsequent breaths. Asymmetry in the conducting airways appears to have little effect on the phase III slope. That the increase in slope appears to be largely controlled by interregional ventilation inhomogeneities should be of interest to those wishing to use multibreath washouts to detect the location of the structural abnormalities within the lung.     

A fractal analysis of the radial distribution of bronchial capillaries around large airways.

We analyzed published measurements of the bronchial circulation and airway wall (Anderson JC, Bernard SL, Luchtel DL, Babb AL, and Hlastala MP. Respir Physiol Neurobiol 132: 329-339, 2002) and determined that the radial distribution of bronchial capillary cross-sectional area was fractal. We limited our analysis to bronchial capillaries, diameter < or =10 mum, that resided between the epithelial basement membrane and adventitia-alveolar boundary, the airway wall tissue. Thirteen different radial distributions of capillary-to-tissue area were constructed simply by changing the number of annuli (i.e., the annular size) used to form each distribution. For the 13 distributions created, these annuli ranged in size from to of the size of the airway wall area. Radial distributions were excluded from the fractal analysis if the sectioning procedure resulted in an annulus with a radial thickness less than the diameter of a capillary. To determine the fractal dimension for a given airway, the coefficient of variation (CV) for each distribution was calculated, and ln(CV) was plotted against the logarithm of the relative piece area. For airways with diameter >2.4 mm, this relationship was linear, which indicated the radial distribution of bronchial capillary cross-sectional area was fractal with an average fractal dimension of 1.27. The radial distribution of bronchial capillary cross-sectional area was not fractal around airways with diameter <1.5 mm. We speculated on how the fractal nature of this circulation impacts the distribution of bronchial blood flow and the efficiency of mass transport during health and disease. A fractal analysis can be used as a tool to quantify and summarize investigations of the bronchial circulation.     

Mean airway pressure and mean alveolar pressure during high-frequency jet ventilation in rabbits

Mean airway pressure underestimates mean alveolar pressure during high-frequency oscillatory ventilation. We hypothesized that high inspiratory flows characteristic of high-frequency jet ventilation may generate greater inspiratory than expiratory pressure losses in the airways, thereby causing mean airway pressure to overestimate, rather than underestimate, mean alveolar pressure. To test this hypothesis, we ventilated anesthetized paralyzed rabbits with a jet ventilator at frequencies of 5, 10, and 15 Hz, constant inspiratory-to-expiratory time ratio of 0.5 and mean airway pressures of 5 and 10 cmH2O. We measured mean total airway pressure in the trachea with a modified Pitot probe, and we estimated mean alveolar pressure as the mean pressure corresponding in the static pressure-volume relationship to the mean volume of the respiratory system measured with a jacket plethysmograph. We found that mean airway pressure was similar to mean alveolar pressure at frequencies of 5 and 10 Hz but overestimated it by 1.1 and 1.4 cmH2O at mean airway pressures of 5 and 10 cmH2O, respectively, when frequency was increased to 15 Hz. We attribute this finding primarily to the combined effect of nonlinear pressure frictional losses in the airways and higher inspiratory than expiratory flows. Despite the nonlinearity of the pressure-flow relationship, inspiratory and expiratory net pressure losses decreased with respect to mean inspiratory and expiratory flows at the higher rates, suggesting rate dependence of flow distribution. Redistribution of tidal volume to a shunt airway compliance is thought to occur at high frequencies.(ABSTRACT TRUNCATED AT 250 WORDS)     

Identifying airways responsible for heterogeneous ventilation and mechanical dysfunction in asthma: an image functional modeling approach.

We present an image functional modeling approach, which synthesizes imaging and mechanical data with anatomically explicit computational models. This approach is utilized to identify the relative importance of small and large airways in the simultaneous deterioration of mechanical function and ventilation in asthma. Positron emission tomographic (PET) images provide the spatial distribution and relative extent of ventilation defects in asthmatic subjects postbronchoconstriction. We also measured lung resistance and elastance from 0.15 to 8 Hz. The first step in image functional modeling involves mapping ventilation three-dimensional images to the computational model and identifying the largest sized airways of the model that, if selectively constricted, could precisely match the size and anatomic location of ventilation defects imaged by PET. In data from six asthmatic subjects, these airways had diameters <2.39 mm and mostly <0.44 mm. After isolating and effectively closing airways in the model associated with these ventilation defects, we imposed constriction with various means and standard deviations to the remaining airways to match the measured lung resistance and elastance from the same subject. Our results show that matching both the degree of mechanical impairment and the size and location of the PET ventilation defects requires either constriction of airways <2.4 mm alone, or a simultaneous constriction of small and large airways, but not just large airways alone. Also, whereas larger airway constriction may contribute to mechanical dysfunction during asthma, degradation in ventilation function requires heterogeneous distribution of near closures confined to small airways.     

Zero-stress state of intra- and extraparenchymal airways from human, pig, rabbit, and sheep lung.

Alterations in airway wall anatomic properties and the consequential effects on airway narrowing have been assessed by use of computational models. In these models, it is generally assumed that at zero transmural pressure the airway wall exists in a zero-stress state. Many studies have shown that this is often not the case, as evidenced by a nonzero opening angle. In this study, we measured the opening angle of airway rings at zero transmural pressure to test this assumption. The airway tree was dissected from human, pig, sheep, and rabbit lungs. Airways were excised from the tree, and the opening angle was measured. There were obvious species and regional differences in opening angle. Rabbit airways from both extraparenchymal and intraparenchymal sites exhibited marked opening angles (7-82 degrees). Extraparenchymal airways from sheep had large opening angles (up to 50 degrees), but ovine intraparenchymal airways had small opening angles. Measurable opening angles were rarely observed in human and porcine airways of any size. The assumption of a stable zero-stress state at zero transmural pressure is therefore valid for human and porcine, but not rabbit and sheep, airways.     

Airway closure on imaging relates to airway hyperresponsiveness and peripheral airway disease in asthma

The regional pattern and extent of airway closure measured by three-dimensional ventilation imaging may relate to airway hyperresponsiveness (AHR) and peripheral airways disease in asthmatic subjects. We hypothesized that asthmatic airways are predisposed to closure during bronchoconstriction in the presence of ventilation heterogeneity and AHR. Fourteen asthmatic subjects (6 women) underwent combined ventilation single photon emission computed tomography/computed tomography scans before and after methacholine challenge. Regional airway closure was determined by complete loss of ventilation following methacholine challenge. Peripheral airway disease was measured by multiple-breath nitrogen washout from which S(cond) (index of peripheral conductive airway abnormality) was derived. Relationships between airway closure and lung function were examined by multiple-linear regression. Forced expiratory volume in 1 s was 87.5 ± 15.8% predicted, and seven subjects had AHR. Methacholine challenge decreased forced expiratory volume in 1 s by 23 ± 5% and increased nonventilated volume from 16 ± 4 to 29 ± 13% of computed tomography lung volume. The increase in airway closure measured by nonventilated volume correlated independently with both S(cond) (partial R(2) = 0.22) and with AHR (partial R(2) = 0.38). The extent of airway closure induced by methacholine inhalation in asthmatic subjects is greater with increasing peripheral airways disease, as measured by ventilation heterogeneity, and with worse AHR.     

Effect of PEEP on induced constriction is enhanced in decorin-deficient mice

Decorin (Dcn), a small leucine-rich proteoglycan, is present in the extracellular matrix of the airways and lung tissues, contributes to lung mechanical properties, and its deposition is altered in asthma. The effect of Dcn deficiency on airway parenchymal interdependence was examined during induced bronchoconstriction. Studies were performed in C57Bl/6 mice in which the Dcn gene was disrupted by targeted deletion (Dcn(-/-)) and in wild-type controls (Dcn(+/+)). Mice were mechanically ventilated, and respiratory system impedance was measured during in vivo ventilation at positive end-expiratory pressure (PEEP) = 2 and 10 cmH(2)0, before and after aerosol delivery of methacholine (MCh). Length vs. tension curves in isolated tracheal rings were measured in vitro. Dcn distribution in +/+ mice airways was characterized by immunofluorescence; differences in collagen structure in Dcn(+/+) and Dcn(-/-) mouse lungs was examined by electron microscopy. MCh caused similar increases in airway resistance (Raw) and tissue elastance (H) in Dcn(+/+) and Dcn(-/-) mice. During MCh-induced constriction, increasing PEEP caused a decrease in Raw that was greater in Dcn(-/-) mice and a decrease in H in Dcn(-/-) mice only. Tracheal ring compliance was greater in Dcn (-/-) mice. Imaging studies showed that Dcn was deposited primarily in the airway adventitial layer in Dcn(+/+) mice; in Dcn(-/-) mice, collagen had an irregular appearance, especially in the lung periphery. These results show that lack of Dcn alters the normal interaction between airways and lung parenchyma; in asthma, changes in Dcn could potentially contribute to abnormal airway physiology.     

Role of collateral paths in long-range diffusion in lungs.

The long-range apparent diffusion coefficient (LRADC) of (3)He gas in lungs, measured over times of several seconds and distances of 1-3 cm, probes the connections between the airways. Previous work has shown the LRADC to be small in health and substantially elevated in emphysema, reflecting tissue destruction, which is known to create collateral pathways. To better understand what controls LRADC, we report computer simulations and measurements of (3)He gas diffusion in healthy lungs. The lung is generated with a random algorithm using well-defined rules, yielding a three-dimensional set of nodes or junctions, each connected by airways to one parent node and two daughters; airway dimensions are taken from published values. Spin magnetization in the simulated lung is modulated sinusoidally, and the diffusion equation is solved to 1,000 s. The modulated magnetization decays with a time constant corresponding to an LRADC of approximately 0.001 cm(2)/s, which is smaller by a factor of approximately 20 than the values in healthy lungs measured here and previously in vivo and in explanted lungs. It appears that collateral gas pathways, not present in the simulations, are functional in healthy lungs; they provide additional and more direct routes for long-range motion than the canonical airway tree. This is surprising, inasmuch as collateral ventilation is believed to be physiologically insignificant in healthy lungs. We discuss the effect on LRADC of small collateral connections through airway walls and rule out other possible mechanisms. The role of collateral paths is supported by measurements of smaller LRADC in pigs, where collateral ventilation is known to be smaller.     

Airway anesthesia effects on hypercapnic breathing pattern in humans

Minute ventilation (VE) and breathing pattern during an abrupt increase in fractional CO2 were compared in 10 normal subjects before and after airway anesthesia. Subjects breathed 7% CO2-93% O2 for 5 min before and after inhaling aerosolized lidocaine. As a result of airway anesthesia, VE and tidal volume (VT) were greater during hypercapnia, but there was no effect on inspiratory time (TI). Therefore, airway anesthesia produced an increase in mean inspiratory flow (VT/TI) during hypercapnia. The increase in VT/TI was compatible with an increase in neuromuscular output. There was no effect of airway anesthesia on the inspiratory timing ratio or the shape and position of the curve relating VT and TI. We also compared airway resistance (Raw), thoracic gas volume, forced vital capacity, forced expired volume at 1s, and maximum midexpiratory flow rate before and after airway anesthesia. A small (0.18 cmH2O X l-1 X s) decrease in Raw occurred after airway anesthesia that did not correlate with the effect of airway anesthesia on VT/TI. We conclude that airway receptors accessible to airway anesthesia play a role in hypercapnic VE.     

Impedance, gas mixing, and bimodal ventilation in constricted lungs.

To evaluate the effect of increasing smooth muscle activation on the distribution of ventilation, lung impedance and expired gas concentrations were measured during a 16-breath He-washin maneuver in five nonasthmatic subjects at baseline and after each of three doses of aerosolized methacholine. Values of dynamic lung elastance (El,dyn), the curvature of washin plots, and the normalized slope of phase III (S(N)) were obtained. At the highest dose, El,dyn was 2.6 times the control value and S(N) for the 16th breath was 0.65 liter(-1). A previously described model of a constricted terminal airway was extended to include variable muscle activation, and the extended model was tested against these data. The model predicts that the constricted airway has two stable states. The impedances of the two stable states are independent of smooth muscle activation, but driving pressure and the number of airways in the high-resistance state increase with increasing muscle activation. Model predictions and experimental data agree well. We conclude that, as a result of the bistability of the terminal airways, the ventilation distribution in the constricted lung is bimodal.     

Contributions of diffusion jet flow and cardiac activity to regional ventilation in CFV.

The magnitude and regional distribution of local gas transport during constant-flow ventilation (CFV) were quantified by imaging the washout of nitrogen 13 (13NN) from anesthetized and paralyzed mongrel dogs with positron emission tomography. Equal jet flows, through two 2-mm-ID bronchial catheters 1 cm distal to the carina, were adjusted to provide eucapnic CFV (total flow = 57.6 ml.s-1.kg-1). Basal, midheart, and apical transverse sections were studied in supine and prone anesthetized dogs. The ventilation per unit volume (sV) of selected areas was computed from local 13NN concentration vs. time curves during washout. To separate the regional contributions of CFV and cardiogenic oscillation to enhanced molecular diffusion, additional supine dogs were also studied during unilateral CFV. In this protocol the CFV jet flow was delivered to a single lung while the contralateral lung was left apneic. For each lung, washout data were obtained under CFV and apnea both living and postmortem animals. The local contributions of diffusion, CFV jet effects, and cardiac activity to gas transport were evaluated and tested for additive and multiplicative synergistic interactions. The regional distribution of gas transport during CFV was found to be highly nonuniform and characterized by higher ventilation to regions located close to the main bronchi and those located in the direction in which the CFV jet pointed. No major differences were observed between supine and prone positions. This regional pattern of ventilation distribution was found to be the result of complementary and nearly multiplicative interaction between the regional effects of the CFV jet, concentrated in the central airways, and the preferential cardiogenic gas transport enhancement in ventral regions close to the heart. The data were also analyzed with a model that divides the regional diffusive gas transport resistance into a central component, affected by the CFV jet, and a peripheral component, affected only by cardiac activity. This analysis showed substantial regional heterogeneities in the effects of the different gas transport mechanisms, which are consistent with the geometry of the bronchial tree and the location of the heart in the dog. The results indicate that regional nonuniformities must be considered when modeling gas transport in CFV.     

A technique for quantitating airway size from bronchial casts

To develop a technique for quantitating the size of airways at various positions in the bronchial tree, we analyzed casts of formalin-fixed excised lungs of five mature male ferrets. The left lower lobe of each cast was dissected, the diameter and position of each terminal bronchiole were entered into a computer programmed to reconstruct the airway system, and the cross-sectional areas of 120 conducting airways were measured. The fraction of the lobe served by each measured airway was estimated by dividing the sum of the squared diameters of the terminal bronchioles subtended by that airway by the summed squared diameters of all terminal bronchioles in the lobe. In each cast the relationship between an airway's cross-sectional area (Y) and the fraction of the lobe it was estimated to subtend (X) was described (0.91 less than R2 less than 0.95) by the expression ln(Y) = A + B ln(X) + C [ln(X)]2. Linear regression of ln(Y) on ln(X) for 30-50 airways estimated to serve fractions of the lobe around each of three arbitrarily selected levels (airways serving 0.7, 2.2, and 9.5% of the lobe) was adequate to characterize the area of airways at each level in each of the five animals with 95% confidence intervals narrower than 8% of the estimated area. Variability of airway size at each level among the five casts was modest, suggesting that this technique identified analagous airways in the various animals. Interindividual variability did not increase when the data were reanalyzed with terminal units defined on the basis of airway diameters rather than on the morphological identification of terminal bronchioles.