Simulation of the uptake of a reactive gas in a rat respiratory tract model with an asymmetric tracheobronchial region patterned on complete conducting airway cast data

Comparative Simulation of Gas Transport in Airway Models of Rat, Dog, and Human

Although a number of animal studies have been conducted to investigate the toxic effects of gaseous pollutants on human airways, the anatomical and physiological differences between animals and humans represent a challenge in extrapolating the animal data to humans. The aim of this study was to examine how interspecies anatomical and physiological differences influence the transport of the inhaled gases throughout the airways and alveoli. We designed mathematical airway models of three mammalian species, rats, dogs, and humans, in which interspecies differences in airway dimensions and respiratory patterns were taken into account. We then simulated the bulk flow of three gases (ozone [O(3)], nitrogen dioxide [NO(2)], and sulfur dioxide [SO(2)]) and obtained the intra-airway concentrations of the gases and the amount absorbed using these models. For all three gases, both real-time and mean concentrations in the upper and lower airways were higher in humans when compared with rats and dogs. For example, the mean concentration of O(3) in the 5th bronchi of humans was 3 and 12 times higher than in rats and dogs, respectively. Similarly, the amount of absorbed gases corrected for airway surface area was again higher in the upper and lower airways of humans than the other two species. Sensitivity analysis indicated that tidal volume, respiratory rate, and surface area of the upper and lower airways had significant impact on the results. In conclusion, kinetics of inhaled gaseous substances vary substantially among animals and humans, and such variations are, at least partially, the result of anatomical and physiological differences in their airways.

Three-Dimensional Simulations of Reactive Gas Uptake in Single Airway Bifurcations

The pattern of lung injury induced by the inhalation of ozone (O(3)) depends on the dose delivered to different tissues in the airways. This study examined the distribution of O(3) uptake in a single, symmetrically branched airway bifurcation. Reaction in the epithelial lining fluid was assumed to be so rapid that O(3) concentration was negligible along the entire surface of the bifurcation wall. Three-dimensional numerical solutions of the continuity, Navier-Stokes and convection-diffusion equations were obtained for steady inspiratory and expiratory flows at Reynolds numbers ranging from 100 to 500. The total rate of O(3) uptake was found to increase with increasing flow rate during both inspiration and expiration. Hot spots of O(3) flux appeared at the carina of the bifurcation for virtually all inspiratory and expiratory Reynolds numbers considered in the simulations. At the lowest expiratory Reynolds number, however, the location of the maximum flux was shifted to the outer wall of the daughter branch. For expiratory flow, additional hot spots of flux were found on the parent branch wall just downstream of the branching region. In all cases, O(3) uptake in the single bifurcation was larger than that in a straight tube of equal inlet radius and wall surface area. This study provides insight into the effect of flow conditions on O(3) uptake and dose distribution in individual bifurcations.

Dose response for formaldehyde-induced cytotoxicity in the human respiratory tract

Human studies of the sensory irritant effects of formaldehyde are complicated by the subjective nature of some clinical endpoints. This limits the usefulness of such studies for quantitative noncancer risk assessment of airborne formaldehyde. Objective measures of the noncancer effects of formaldehyde, such as the rate of regenerative cellular proliferation (RCP) secondary to cytolethality, can be obtained from laboratory animals but present the challenge of interspecies extrapolation of the data. To the extent that uncertainties associated with this extrapolation can be reduced, however, dose-response data obtained in laboratory animals are a viable alternative to clinical studies. Here, we describe the extrapolation of dose-response data for RCP from F344 rats to humans. Rats inhaled formaldehyde (0, 0.7, 2.0, 6.0, 10, and 15 ppm) 6 h/day, 5 days/week for up to 2 years. The dose response for RCP was J-shaped, with the rates of RCP at 0.7 and 2.0 ppm below but not statistically different from control, while the rates at the higher concentrations were significantly greater than control. Both the raw J-shaped data and a hockey-stick-shaped curve fitted to the raw data were used for predicting the human dose response for RCP. Cells lining the nasal airways of F344 rats and rhesus monkeys are comparably sensitive to the cytolethal effects of inhaled formaldehyde, suggesting that the equivalent human cells are also likely to be comparably sensitive. Using this assumption, the challenge of rat-to-human extrapolation was reduced to accurate prediction of site-specific flux of formaldehyde from inhaled air into the tissue lining the human respiratory tract. A computational fluid dynamics model of air flow and gas transport in the human nasal airways was linked to a typical path model of the human lung to provide site-specific flux predictions throughout the respiratory tract. Since breathing rate affects formaldehyde dosimetry, cytotoxicity dose-response curves were predicted for three standard working levels. With the most vigorous working level, the lowest concentrations of formaldehyde predicted to exert any cytotoxic effects in humans were 1.0 and 0.6 ppm, for the J-shaped and hockey-stick-shaped RCP curves, respectively. The predicted levels of response at the lowest effect concentrations are smaller than can be measured clinically. Published literature showing that the cytotoxicity of inhaled formaldehyde is related to exposure level rather than to duration of exposure suggests that the present analysis is a reasonable basis for derivation of standards for continuous human exposure.

NO2 interfacial transfer is reduced by phospholipid monolayers

Nitrogen dioxide (NO2) is a ubiquitous, pollutant gas that produces a broad range of pathological and physiological effects on the lung. Absorption of inhaled NO2 is coupled to near-interfacial reactions between the solute gas and constituents of the airway and alveolar epithelial lining fluid. Although alveolar surfactant imparts limited resistance to respiratory gas exchange compared with that contributed by either the pulmonary membrane or uptake in red blood cells, resistance to NO2 flux could have a significant effect on NO2 absorption kinetics. To investigate the effect of interfacial surfactant on NO2 absorption, we designed an apparatus permitting exposure of variably compressed monolayers. Our results suggest that compressed monolayers enriched in 1,2-dipalmitoyl-sn-3-glycero-phosphocholine present significant resistance to NO2 absorption even at surface tensions greater than those achieved in vivo. However, monolayers composed of pure unsaturated phospholipids failed to alter NO2 absorption significantly when compressed, in spite of similar reductions in surface tension. The results demonstrate that phospholipid monolayers appreciably limit NO2 absorption and further that monolayer-induced resistance to NO2 flux is related to physicochemical properties of the film itself rather than alterations within the aqueous and gas phases. On the basis of these findings, we propose that pulmonary surfactant may influence the intrapulmonary gas phase distribution of inhaled NO2.

Dosimetry modeling of inhaled formaldehyde: the human respiratory tract

Formaldehyde (HCHO), which has been shown to be a nasal carcinogen in rats and mice, is used widely and extensively in various manufacturing processes. Studies in rhesus monkeys suggest that the lower respiratory tract may be at risk and some epidemiologic studies have reported an increase in lung cancer associated with HCHO; other studies have not. Thus, an assessment of possible human risk to HCHO exposure based on dosimetry information throughout the respiratory tract (RT) is desirable. To obtain dosimetry estimates for a risk assessment, two types of models were used. The first model (which is the subject of another investigation) used computational fluid dynamics (CFD) to estimate local fluxes in a 3-dimensional model of the nasal region. The subject of the present investigation (the second model) applied a 1-dimensional equation of mass transport to each generation of an adult human symmetric, bifurcating Weibel-type RT anatomical model, augmented by an upper respiratory tract. The two types of modeling approaches were made consistent by requiring that the 1-dimensional version of the nasal passages have the same inspiratory air-flow rate and uptake during inspiration as the CFD simulations for 4 daily human activity levels. Results obtained include the following: (1) More than 95% of the inhaled HCHO is predicted to be retained by the RT. (2) The CFD predictions for inspiration, modified to account for the difference in inspiration and complete breath times, are a good approximation to uptake in the nasal airways during a single breath. (3) In the lower respiratory tract, flux is predicted to increase for several generations and then decrease rapidly. (4) Compared to first pulmonary region generation fluxes, the first few tracheobronchial generations fluxes are over 1000 times larger. Further, there is essentially no flux in the alveolar sacs. (5) Predicted fluxes based on the 1-dimensional model are presented that can be used in a biologically based dose-response model for human carcinogenesis. Use of these fluxes will reduce uncertainty in a risk assessment for formaldehyde carcinogenicity.

Photochemical oxidants: state of the science

Atmospheric photochemical processes resulting in the production of tropospheric ozone (O(3)) and other oxidants are described. The spatial and temporal variabilities in the occurrence of surface level oxidants and their relationships to air pollution meteorology are discussed. Models of photooxidant formation are reviewed in the context of control strategies and comparisons are provided of the air concentrations of O(3) at select geographic locations around the world. This overall oxidant (O(3)) climatology is coupled to human health and ecological effects. The discussion of the effects includes both acute and chronic responses, mechanisms of action, human epidemiological and plant population studies and briefly, efforts to establish cause-effect relationships through numerical modeling. A short synopsis is provided of the interactive effects of O(3) with other abiotic and biotic factors. The overall emphasis of the paper is on identifying the current uncertainties and gaps in our understanding of the state of the science and some suggestions as to how they may be addressed.

O3-induced formation of bioactive lipids: estimated surface concentrations and lining layer effects

Recent evidence suggests that inhaled ozone (O3) does not induce toxicity via direct epithelial interactions. Reactions with epithelial lining fluid (ELF) constituents limit cellular contact and generate products, including lipid ozonation products, postulated to initiate pathophysiological cascades. To delineate specific aspects of lipid ozonation product formation and to estimate in situ surface concentrations, we studied the O3 absorption characteristics of ELF constituent mixtures and measured hexanal, heptanal, and nonanal yields as a function of ascorbic acid (AH2) concentration. Exposures of isolated rat lungs, bronchoalveolar lavage fluid (BALF) and egg phosphatidylcholine (PC) liposomes were conducted. 1) O3 absorption by AH2, uric acid, and albumin exceeded that by egg PC and glutathione. O3 reaction with egg PC occurred when AH2 concentrations were reduced. 2) Aldehydes were produced in low yield during lung and BALF exposures in a time- and O3 concentration-dependent manner. 3) Diminishing BALF AH2 content lowered O3 uptake but increased aldehyde yields. Conversely, AH2 addition to egg PC increased O3 uptake but reduced aldehyde yields. Estimations of bioactive ozonation and autoxidation product accumulation within the ELF suggested possible nanomolar to low micromolar concentrations. The use of reaction products as metrics of O3 exposure may have intrinsic sensitivity and specificity limitations. Moreover, due to the heterogenous nature of O3 reactions within the ELF, dose-response relationships may not be linear with respect to O3 absorption.