Effect of airways constriction on exhaled nitric oxide.

While airway constriction has been shown to affect exhaled nitric oxide (NO), the mechanisms and location of constricted airways most likely to affect exhaled NO remain obscure. We studied the effects of histamine-induced airway constriction and ventilation heterogeneity on exhaled NO at 50 ml/s (Fe(NO,50)) and combined this with model simulations of Fe(NO,50) changes due to constriction of airways at various depths of the lung model. In 20 normal subjects, histamine induced a 26 +/- 15(SD)% Fe(NO,50) decrease, a 9 +/- 6% forced expiratory volume in 1 s (FEV(1)) decrease, a 19 +/- 9% mean forced midexpiratory flow between 25% and 75% forced vital capacity (FEF(25-75)) decrease, and a 94 +/- 119% increase in conductive ventilation heterogeneity. There was a significant correlation of Fe(NO,50) decrease with FEF(25-75) decrease (P = 0.006) but not with FEV(1) decrease or with increased ventilation heterogeneity. Simulations confirmed the negligible effect of ventilation heterogeneity on Fe(NO,50) and showed that the histamine-induced Fe(NO,50) decrease was due to constriction, with associated reduction in NO flux, of airways located proximal to generation 15. The model also indicated that the most marked effect of airways constriction on Fe(NO,50) is situated in generations 10-15 and that airway constriction beyond generation 15 markedly increases Fe(NO,50) due to interference with the NO backdiffusion effect. These mechanical factors should be considered when interpreting exhaled NO in lung disease.     

Lung inflation does not increase maximal expiratory flow during induced obstruction in the dog

A deep inflation (DI) reverses induced bronchoconstriction in normal human subjects whether assessed by airway resistance before and after a DI or by isovolumic maximal expiratory flows (Vmax) from partial expiratory flow-volume (PEFV) vs. maximum expiratory flow-volume (MEFV) maneuvers. These observations suggest that with induced constriction the hysteresis of airways exceeds that of the parenchyma. In contrast with humans, a previous study of ours on dogs indicated that induced increases in airway resistance were unaffected by DI, suggesting that hysteresis of airways and parenchyma were equal. We hypothesized therefore that in constricted dog lungs, any differences that might arise in isovolumic Vmax between PEFV and MEFV maneuvers would not be due to changes in airway caliber but rather would be wholly determined by isovolumic differences in deflational recoil pressures. Recoil pressures were dynamically measured using six separate alveolar capsules in each of six dogs. At base line there were no significant differences between isovolumic recoil pressures or maximal flows with volume history, suggesting equal degrees of airway and parenchymal hysteresis. After histamine-induced constriction there were also no isovolumic differences in flows, but due to striking nonhomogeneities in dynamic recoil pressure among alveolar capsules, it was not possible to express a single meaningful recoil pressure pertinent to the lungs as a whole. These findings are consistent with the idea that isovolumic comparisons of Vmax serve as a reasonable indicator of changes in the relative degree of airway and parenchymal hysteresis.     

Effect of resting smooth muscle length on contractile response in resistance airways

We studied the effect of resting smooth muscle length on the contractile response of the major resistance airways (generations 0-5) in 18 mongrel dogs in vivo using tantalum bronchography. Dose-response curves to 10(-10) to 10(-7) mol/kg methacholine (MCh) were generated [at functional residual capacity (FRC)] by repeated intravenous bolus administration using tantalum bronchography after each dose. Airway constriction varied substantially with dose-equivalent stimulation and varied sequentially from trachea (8.8 +/- 2.2% change in airway diam) to fifth-generation bronchus (49.8 +/- 3.0%; P less than 0.001). Length-tension curves were generated for each airway to determine the airway diameter (i.e., resting in situ smooth muscle length) at which maximal constriction was elicited using bolus intravenous injection of 10(-8) mol/kg MCh. A Frank-Starling relationship was obtained for each airway; the transpulmonary pressure at which maximal constriction was elicited increased progressively from 2.50 +/- 1.12 cmH2O for trachea (approximately FRC) to 18.3 +/- 1.05 cmH2O for fifth-generation airways (approximately 50% TLC) (P less than 0.001). A similar relationship was obtained when change in airway diameter was plotted as a function of airway radius. We demonstrate substantial heterogeneity in the lung volumes at which maximal constriction is elicited and in distribution of parasympathomimetic constriction within the first few generations of resistance bronchi. Our data also suggest that lung hyperinflation may lead to augmented airway contractile responses by shifting resting smooth muscle length toward optimum resting smooth muscle length.     

Airway remodeling in asthma amplifies heterogeneities in smooth muscle shortening causing hyperresponsiveness.

Although airway remodeling and inflammation in asthma can amplify the constriction response of a single airway, their influence on the structural changes in the whole airway network is unknown. We present a morphometric model of the human lung that incorporates cross-sectional wall areas corresponding to the adventitia, airway smooth muscle (ASM), and mucosa for healthy and mildly and severely asthmatic airways and the influence of parenchymal tethering. A heterogeneous ASM percent shortening stimulus is imposed, causing distinct constriction patterns for healthy and asthmatic airways. We calculate lung resistance and elastance from 0.1 to 5 Hz. We show that, for a given ASM stimulus, the distribution of wall area in asthmatic subjects will amplify not only the mean but the heterogeneity of constriction in the lung periphery. Moreover, heterogeneous ASM shortening that would produce only mild changes in the healthy lung can cause hyperresponsive changes in lung resistance and elastance at typical breathing rates in the asthmatic lung, even with relatively small increases in airway resistance. This condition arises when airway closures occur randomly in the lung periphery. We suggest that heterogeneity is a crucial determinant of hyperresponsiveness in asthma and that acute asthma is more a consequence of extensive airway wall inflammation and remodeling, predisposing the lung to produce an acute pattern of heterogeneous constriction.     

Airway dynamics in transition between peak and maximal expiratory flow

Expiratory flow-volume curves with periodic interruption of flow showed flow transients exceeding maximal flow (Vmax) measured on the maximum expiratory flow-volume (MEFV) curve in a mechanical lung model and in five tracheotomized, vagotomized, open-chest, anesthetized dogs. Direct measurement of flow from the collapsing model airway showed that the volume of the flow transients in excess of the MEFV envelope was greater than that from the collapsing airway. Determination of wave-speed flows from local airway transmural pressure-area curves (J. Appl. Physiol. 52: 357-369, 1982) and photography of the airway led to the following conclusions. Flow transients exceeding Vmax are wave-speed flows determined by an initial and unstable configuration of the flow-limiting segment (FLS) with maximum compression in the midportion. The drop in flow from the peak to the following plateau is due to development of a more stable airway configuration with maximum compression at the mouthward end with a smaller area and a smaller maximal flow. When FLS jumps to a more peripheral position, the more distal airways may pass through similar configurational changes that are responsible for the sudden decrease of flow (the "knee") seen on most MEFV curves from dogs.     

Complex airway behavior and paradoxical responses to bronchoprovocation.

Heterogeneity of airway constriction and regional ventilation in asthma are commonly studied under the paradigm that each airway's response is independent from other airways. However, some paradoxical effects and contradictions in recent experimental and theoretical findings suggest that considering interactions among serial and parallel airways may be necessary. To examine airway behavior in a bronchial tree with 12 generations, we used an integrative model of bronchoconstriction, including for each airway the effects of pressure, tethering forces, and smooth muscle forces modulated by tidal stretching during breathing. We introduced a relative smooth muscle activation factor (T(r)) to simulate increasing and decreasing levels of activation. At low levels of T(r), the model exhibited uniform ventilation and homogeneous airway narrowing. But as T(r) reached a critical level, the airway behavior suddenly changed to a dual response with a combination of constriction and dilation. Ventilation decreased dramatically in a group of terminal units but increased in the rest. A local increase of T(r) in a single central airway resulted in full closure, while no central airway closed under global elevation of T(r). Lung volume affected the response to both local and global stimulation. Compared with imaging data for local and global stimuli, as well as for the time course of airway lumen caliber during bronchoconstriction recovery, the model predictions were similar. The results illustrate the relevance of dynamic interactions among serial and parallel pathways in airway interdependence, which may be critical for the understanding of pathological conditions in asthma.     

Airway stability and heterogeneity in the constricted lung.

The effect of bronchoconstriction on airway resistance is known to be spatially heterogeneous and dependent on tidal volume. We present a model of a single terminal airway that explains these features. The model describes a feedback between flow and airway resistance mediated by parenchymal interdependence and the mechanics of activated smooth muscle. The pressure-tidal volume relationship for a constricted terminal airway is computed and shown to be sigmoidal. Constricted terminal airways are predicted to have two stable states: one effectively open and one nearly closed. We argue that the heterogeneity of whole lung constriction is a consequence of this behavior. Airways are partitioned between the two states to accommodate total flow, and changes in tidal volume and end-expiratory pressure affect the number of airways in each state. Quantitative predictions for whole lung resistance and elastance agree with data from previously published studies on lung impedance.     

Interdependent regional lung emptying during forced expiration: a transistor model

We recognized similarities between isovolume pressure-flow curves of the lung and emitter-collector voltage-current characteristics of bipolar transistors, and used this analogy to model expiratory flow limitation in a two-generation branching network with parallel nonhomogeneity. In this model, each of two bronchi empty parenchymal compliances through a common trachea, and each branch includes resistances upstream and downstream of a flow-limiting site. Properties of each airway are specified independently, allowing simulation of differences between the tracheal and bronchial generations and between the parallel bronchial paths. Simulations of four types of parallel asymmetry were performed: unilateral peripheral bronchoconstriction; unilateral central bronchoconstriction; asymmetric redistribution of parenchymal compliance; and unilateral alteration of the bronchial area-transmural pressure characteristic. Our results indicate that multiple axial choke points can exist simultaneously in a symmetric lung when large airway opening-pleural pressure gradients exist; despite severe nonhomogeneity of regional lung emptying, flow interdependence among parallel branches tends to maintain a near normal configuration of the overall maximal expiratory flow-volume (MEFV) curve throughout a large fraction of the vital capacity; and sudden changes of slope of the MEFV curve ("knees" or "bumps") may reflect choking in one branch in a nonuniform lung, but need not be obvious even when severe heterogeneity of lung emptying exists.     

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.     

Maximum expiratory flow-volume curves during short periods of microgravity

To elucidate the effect of normal gravitation on the shape of the maximum expiratory flow-volume (MEFV) curve, we studied nine normal subjects in a National Aeronautics and Space Administration microgravity research aircraft. They performed multiple MEFV maneuvers at 0, 1, and approximately 2 G. The MEFV curves for each subject were filtered, aligned at residual volume, and ensemble averaged to produce an average MEFV curve for each state, allowing differences to be studied. Most subjects showed a decrease in the forced vital capacity at 0 G, which we attribute to an increased intrathoracic blood volume. In most of these subjects, the mean lung volume associated with a given flow was lower at 0 G over about the upper half of the vital capacity. This is similar to the change previously reported during headout immersion and is consistent with the known effect of engorgement of the lung with blood on elastic recoil. There were also consistent but highly individual changes in the position and magnitude of detailed features of the curve, the individual patterns being similar to those previously reported on transition from the erect to the supine position. This supports the idea that the location and motion of choke points that determine the detailed individual configuration of MEFV curves can be significantly influenced by gravitational forces, presumably via the effects of change in longitudinal tension on local airway pressure-diameter behavior and thus wave speed.     

Distribution of airway narrowing responses across generations and at branching points,assessed in vitro by anatomical optical coherence tomography

Background

Previous histological and imaging studies have shown the presence of variability in the degree of bronchoconstriction of airways sampled at different locations in the lung (i.e., heterogeneity). Heterogeneity can occur at different airway generations and at branching points in the bronchial tree. Whilst heterogeneity has been detected by previous experimental approaches, its spatial relationship either within or between airways is unknown.

Methods

In this study, distribution of airway narrowing responses across a portion of the porcine bronchial tree was determined in vitro. The portion comprised contiguous airways spanning bronchial generations (#3-11), including the associated side branches. We used a recent optical imaging technique, anatomical optical coherence tomography, to image the bronchial tree in three dimensions. Bronchoconstriction was produced by carbachol administered to either the adventitial or luminal surface of the airway. Luminal cross sectional area was measured before and at different time points after constriction to carbachol and airway narrowing calculated from the percent decrease in luminal cross sectional area.

Results

When administered to the adventitial surface, the degree of airway narrowing was progressively increased from proximal to distal generations (r = 0.80 to 0.98, P < 0.05 to 0.001). This ''serial heterogeneity'' was also apparent when carbachol was administered via the lumen, though it was less pronounced. In contrast, airway narrowing was not different at side branches, and was uniform both in the parent and daughter airways.

Conclusions

Our findings demonstrate that the bronchial tree expresses intrinsic serial heterogeneity, such that narrowing increases from proximal to distal airways, a relationship that is influenced by the route of drug administration but not by structural variations accompanying branching sites.     

Analysis of bronchial mechanics and density dependence of maximal expiratory flow

The computational model for expiratory flow in humans of Lambert and associates (J. Appl. Physiol. Respirat. Environ. Exercise Physiol. 52: 44-56, 1982) was used to investigate the effect of bronchial constrictions in three airway zones on the density dependence of maximal expiratory flow. It was found that constriction of the peripheral airways (less than 2 mm diam) reduced density dependence and increased the volume of isoflow. Constriction of the larger intraparenchymal airways resulted in increased density dependence at low lung volumes and essentially normal values at other volumes. The volume of isoflow was reduced. Extraparenchymal (but intrathoracic) airway constriction caused no change in the volume of isoflow but caused increased density dependence at the higher lung volumes. It was shown that in these model simulations the addition of extraparenchymal constriction to intraparenchymal constriction causes essentially no changes in density dependence. An increased volume of isoflow and significantly decreased density dependence at 50 and 25% vital capacity were produced by simulated constrictions only in the peripheral airways.     

Genetic influence on normal variability of maximum expiratory flow-volume curves.

To determine the importance of genetic influence on the variability of maximum expiratory flow-volume (MEFV) curves in normal individuals, MEFV curves breathing air and a mixture of 80% helium and 20% oxygen (He-O2), lung volumes, specific airway conductance, and closing capacity (CC) were obtained in 10 pairs of identical and 6 pairs of nonidentical twins, all nonsmokers and asymptomatic. For a given pair of identical twins, MEFV curves on air were more similar than those of a pair of nonidentical twins (P less than 0.02). The intrapair differences of identical twins were smaller than nonidentical twins of maximum expiratory flow (Vmax) at 60% of total lung capacity (TLC) on air (P less than 0.001) and on He-O2 (P less than 0.01). However, intrapair differences of Vmax at 40% TLC and CC were not significantly different in the two groups. Since Vmax at 60% TLC on air and He-O2 are dependent on the geometry of large airways these findings are suggestive that the geometry of large airways may be related to genetic factors. The relationship of the geometry of the peripheral airways and genetic factors has not been defined.     

Configuration of maximum expiratory flow-volume curve: model experiments with physiological implications

A two-compartment mechanical model of the lungs was constructed with two parallel peripheral and collapsible bronchi in series with one central and collapsible trachea. Maximal expiratory flow-volume (MEFV) curves similar to those obtained in most dogs and in some humans could be produced: a peak followed by a gently sloping plateau ending in a knee, where flow suddenly fell to a much smaller value approaching zero rather slowly over the last 25 to 50% of the expired vital capacity. It was shown that flow before the knee was limited in the trachea, and after the knee it was limited in the bronchi. Two patterns of changes in the configuration of the MEFV curve could be observed. Pattern of changes affecting the central airway, at a given volume, maximal flow during the first part of the expiration (i.e., before the knee) is decreased; the knee occurs at a lower lung volume; the flow at the beginning of the knee is decreased. This pattern was observed with the following interventions: decreased cross-sectional area of the trachea (partial obstruction); decreased axial tension of the trachea; and, increased frictional loss between the trachea and the bronchi. Pattern of changes affecting the airways in the periphery: the knee occurs at a higher lung volume; at a given volume, flow after the knee becomes smaller; the absolute flow at the start of the knee is almost unchanged. This pattern was observed with the following interventions: decreased cross-sectional area of the peripheral airways (partial obstruction); increased frictional loss upstream to the peripheral airways; and, decreased elastic recoil pressure.(ABSTRACT TRUNCATED AT 250 WORDS)     

How does airway inflammation modulate asthmatic airway constriction? An antigen challenge study.

During the late-phase (LP) response to inhaled allergen, mediators from neutrophils and eosinophils are released within the airways, resembling what occurs during an asthma attack. We compared the distribution of obstruction and degree of reversibility that follows a deep inspiration (DI) during early-phase (EP) and LP responses in nine asthmatic subjects challenged with allergen. Heterogeneity of constriction was assayed by determining frequency dependence of dynamic lung resistance and elastance, airway caliber by tracking airway resistance during a DI, and airway inflammation by measuring inflammatory cells in induced sputum postchallenge. Despite a paucity of eosinophils in the sputum at baseline (<1% of nonsquamous cells), asthmatic subjects showed a substantial EP response with highly heterogeneous constriction and reduced capacity to maximally dilate airways. The LP was associated with substantial airway inflammation in all subjects. However, five subjects showed only mild LP constriction, whereas four showed more marked LP constriction characterized by heterogeneous constriction similar to EP. Bronchoconstriction during LP was fully alleviated by administration of a bronchodilator. These findings, together with the impaired bronchodilatory response during a DI, indicate a physiological abnormality in asthma at the smooth muscle level and indicate that airway inflammation in asthma is associated with a highly nonuniform pattern of constriction. These data support the hypothesis that variability in responsiveness among asthmatic subjects derives from intrinsic differences in smooth muscle response to inflammation.     

Relationship between heterogeneous changes in airway morphometry and lung resistance and elastance

Lutchen, Kenneth R., and Heather Gillis. Relationshipbetween heterogeneous changes in airway morphometry and lung resistanceand elastance. J. Appl. Physiol.83(4): 1192-1201, 1997.We present a dog lung model to predictthe relation between inhomogeneous changes in airway morphometry andlung resistance (RL) andelastance (EL) for frequenciessurrounding typical breathing rates. TheRL andEL were sensitive in distinctways to two forms of peripheral constriction. First, when there is alarge and homogeneous constriction, theRL increases uniformly over thefrequency range. The EL israther unaffected below 1 Hz but then increases with frequencies up to5 Hz. This increase is caused by central airway wallshunting. Second, the RL andEL are extremely sensitive to mild inhomogeneous constriction in which a few highly constricted ornearly closed airways occur randomly throughout theperiphery. This results in extreme increases in the levelsand frequency dependence of RLand EL but predominantly attypical breathing rates (<1 Hz). Conversely, theRL andEL are insensitive to highly inhomogeneous airway constriction that does not produce any nearly closed airways. Similarly, alterations in theRL andEL due to central airway wallshunting are not likely until the preponderance of the peripheryconstricts substantially. The RLand EL spectra are far moresensitive to these two forms of peripheral constriction than toconstriction conditions known to occur in the central airways. On thebasis of these simulations, we derived a set of qualitative criteria toinfer airway constriction conditions from RL andEL spectra.

     

An analysis of pollutant gas transport and absorption in pulmonary airways

A mathematical model of ozone absorption, or for any soluble gas that has similar transport properties, is developed for a branching network of liquid-lined cylinders. In particular, we investigate specific flow regimes for finite length tubes where boundary layer phenomena and entrance effects exist in high Reynolds and Peclet (Pe) number airways. The smaller airways which have lower Reynolds and Peclet number flows are modelled by incorporating the detailed analysis found in [10] and modifying it for airways which have alveolated surfaces. We also consider a reacting gas and treat specific regimes where the reaction front is located at the air-liquid interface, within the liquid or at the liquid-tissue interface. Asymptotic methods are used in regions of the tracheobronchial tree where Pe much less than 1 and Pe much greater than 1. In addition, the fact that the radial transport parameter gamma much less than 1 for this toxin, and others such as nitrous oxides, is employed to simplify the analysis. The ozone concentrations, airway absorption and tissue dose are examined as a function of airway generation for several values of the governing parameters. The general result is a maximal dosing in airway generations 17 to 18 that is much larger (up to an order of magnitude) than the predictions of previous theories.     

Inverse modeling of dog airway and respiratory system impedances

Mechanical parameters of the respiratory system are often estimated from respiratory impedances using lumped-element inverse models. One such six-element model is composed of an airway branch [with a resistance (Raw) and inertance (Iaw)] separated from a tissue branch [with a resistance (Rt), inertance (It), and compliance (Ct)] by a shunt compliance representing alveolar gas compression (Cg). Even though the airways are known to have frequency-dependent resistance and inertance, these inverse models have been composed of linear frequency-independent elements. In this study we investigated the use of inverse models where the airway branch was represented by a frequency-independent Raw and Iaw, a Raw that is linearly related to frequency and an Iaw that is independent of frequency, and a system of identical parallel tubes the impedance of which was computed from the tube radius and length. These inverse models were used to analyze airway and respiratory impedances between 2 and 1,024 Hz that were predicted from an anatomically detailed forward model. The forward model represented the airways by an asymmetrically branched network with a terminal impedance representative of known Cg, Rt, It, and Ct. For respiratory impedances between 2 and 128 Hz, all models fit the data reasonably well, and reasonably accurate estimates of Cg, Rt, It, and Ct were extracted from these data. For data above 200 Hz, however, only the multiple-tube model accurately fitted respiratory impedances (Zrs). This model fitted the Zrs data best when composed of 27 tubes, each having a radius of 0.148 cm and a length of 16.5 cm.     

A Distribution-Moment Approximation for Coupled Dynamics of the Airway Wall and Airway Smooth Muscle

Asthma is fundamentally a disease of airway constriction. Due to a variety of experimental challenges, the dynamics of airways are poorly understood. Of specific interest is the narrowing of the airway due to forces produced by the airway smooth muscle wrapped around each airway. The interaction between the muscle and the airway wall is crucial for the airway constriction that occurs during an asthma attack. Although cross-bridge theory is a well-studied representation of complex smooth muscle dynamics, and these dynamics can be coupled to the airway wall, this comes at significant computational cost—even for isolated airways. Because many phenomena of interest in pulmonary physiology cannot be adequately understood by studying isolated airways, this presents a significant limitation. We present a distribution-moment approximation of this coupled system and study the validity of the approximation throughout the physiological range. We show that the distribution-moment approximation is valid in most conditions, and we explore the region of breakdown. These results show that in many situations, the distribution-moment approximation is a viable option that provides an orders-of-magnitude reduction in computational complexity; not only is this valuable for isolated airway studies, but it moreover offers the prospect that rich ASM dynamics might be incorporated into interacting airway models where previously this was precluded by computational cost.