USE OF COMPUTATIONAL FLUID DYNAMICS MODELS FOR DOSIMETRY OF INHALED GASES IN THE NASAL PASSAGES

Computational fluid dynamics (CFD) models of the nasal passages of a rat, monkey, and human are being used (1) to determine important factors affecting nasal uptake, (2) to make interspecies dosimetric comparisons, (3) to provide detailed anatomical information for the rat, monkey, and human nasal passages, and (4) to provide estimates of regional air-phase mass transport coefficients (a measure of the resistance to gas transport from inhaled air to airway walls) in the nasal passages of all three species. For many inhaled materials, lesion location in the nose follows patterns that are both site and species specific. For reactive, water-soluble (Category 1) gases, regional uptake can be a major factor in determining lesion location. Since direct measurement of airflow and uptake is experimentally difficult, CFD models are used here to predict uptake patterns quantitatively in three-dimensional reconstructions of the F344 rat, rhesus monkey, and human nasal passages. In formaldehyde uptake simulations, absorption processes were assumed to be as rapid as possible, and regional flux (transport rate) of inhaled formaldehyde to airway walls was calculated for rats, primates, and humans. For uptake of gases like vinyl acetate and acrylic acid vapors, physiologically based pharmacokinetic uptake models incorporating anatomical and physical information from the CFD models were developed to estimate nasal tissue dose in animals and humans. The use of biologically based models in risk assessment makes sources of uncertainty explicit and, in doing so, allows quantification of uncertainty through sensitivity analyses. Limited resources can then be focused on reduction of important sources of uncertainty to make risk estimates more accurate.     

Dosimetry modeling of inhaled formaldehyde: binning nasal flux predictions for quantitative risk assessment.

Interspecies extrapolations of tissue dose and tumor response have been a significant source of uncertainty in formaldehyde cancer risk assessment. The ability to account for species-specific variation of dose within the nasal passages would reduce this uncertainty. Three-dimensional, anatomically realistic, computational fluid dynamics (CFD) models of nasal airflow and formaldehyde gas transport in the F344 rat, rhesus monkey, and human were used to predict local patterns of wall mass flux (pmol/[mm(2)-h-ppm]). The nasal surface of each species was partitioned by flux into smaller regions (flux bins), each characterized by surface area and an average flux value. Rat and monkey flux bins were predicted for steady-state inspiratory airflow rates corresponding to the estimated minute volume for each species. Human flux bins were predicted for steady-state inspiratory airflow at 7.4, 15, 18, 25.8, 31.8, and 37 l/min and were extrapolated to 46 and 50 l/min. Flux values higher than half the maximum flux value (flux median) were predicted for nearly 20% of human nasal surfaces at 15 l/min, whereas only 5% of rat and less than 1% of monkey nasal surfaces were associated with fluxes higher than flux medians at 0.576 l/min and 4.8 l/min, respectively. Human nasal flux patterns shifted distally and uptake percentage decreased as inspiratory flow rate increased. Flux binning captures anatomical effects on flux and is thereby a basis for describing the effects of anatomy and airflow on local tissue disposition and distributions of tissue response. Formaldehyde risk models that incorporate flux binning derived from anatomically realistic CFD models will have significantly reduced uncertainty compared with risk estimates based on default methods.     

Dosimetry modeling of inhaled formaldehyde: comparisons of local flux predictions in the rat, monkey, and human nasal passages.

Formaldehyde-induced nasal squamous cell carcinomas in rats and squamous metaplasia in rats and rhesus monkeys occur in specific regions of the nose with species-specific distribution patterns. Experimental approaches addressing local differences in formaldehyde uptake patterns and dose are limited by the resolution of dissection techniques used to obtain tissue samples and the rapid metabolism of absorbed formaldehyde in the nasal mucosa. Anatomically accurate, 3-dimensional computational fluid dynamics models of F344 rat, rhesus monkey, and human nasal passages were used to estimate and compare regional inhaled formaldehyde uptake patterns predicted among these species. Maximum flux values, averaged over a breath, in nonsquamous epithelium were estimated to be 2620, 4492, and 2082 pmol/(mm(2)-h-ppm) in the rat, monkey, and human respectively. Flux values predicted in sites where cell proliferation rates were measured as similar in rats and monkeys were also similar, as were fluxes predicted in a region of high tumor incidence in the rat nose and the anterior portion of the human nose. Regional formaldehyde flux estimates are directly applicable to clonal growth modeling of formaldehyde carcinogenesis to help reduce uncertainty in human cancer risk estimates.     

Dosimetry of nasal uptake of water-soluble and reactive gases: A first study of interhuman variability

Certain inhaled chemicals, such as reactive, water-soluble gases, are readily absorbed by the nasal mucosa upon inhalation and may cause damage to the nasal epithelium. Comparisons of the spatial distribution of nasal lesions in laboratory animals exposed to formaldehyde with gas uptake rates predicted by computational models reveal that lesions usually occur in regions of the susceptible epithelium where gas absorption is highest. Since the uptake patterns are influenced by air currents in the nose, interindividual variability in nasal anatomy and ventilation rates due to age, body size, and gender will affect the patterns of gas absorption in humans, potentially putting some age groups at higher risk when exposed to toxic gases. In this study, interhuman variability in the nasal dosimetry of reactive, water-soluble gases was investigated by means of computational fluid dynamics (CFD) models in 5 adults and 2 children, aged 7 and 8 years old. Airflow patterns were investigated for allometrically scaled inhalation rates corresponding to resting breathing. The spatial distribution of uptake at the airway walls was predicted to be nonuniform, with most of the gas being absorbed in the anterior portion of the nasal passages. Under the conditions of these simulations, interhuman variability in dose to the whole nose (mass per time per nasal surface area) due to differences in anatomy and ventilation was predicted to be 1.6-fold among the 7 individuals studied. Children and adults displayed very similar patterns of nasal gas uptake; no significant differences were noted between the two age groups.     

Comparative computational modeling of airflows and vapor dosimetry in the respiratory tracts of rat, monkey, and human

Computational fluid dynamics (CFD) models are useful for predicting site-specific dosimetry of airborne materials in the respiratory tract and elucidating the importance of species differences in anatomy, physiology, and breathing patterns. We improved the imaging and model development methods to the point where CFD models for the rat, monkey, and human now encompass airways from the nose or mouth to the lung. A total of 1272, 2172, and 135 pulmonary airways representing 17±7, 19±9, or 9±2 airway generations were included in the rat, monkey and human models, respectively. A CFD/physiologically based pharmacokinetic model previously developed for acrolein was adapted for these anatomically correct extended airway models. Model parameters were obtained from the literature or measured directly. Airflow and acrolein uptake patterns were determined under steady-state inhalation conditions to provide direct comparisons with prior data and nasal-only simulations. Results confirmed that regional uptake was sensitive to airway geometry, airflow rates, acrolein concentrations, air:tissue partition coefficients, tissue thickness, and the maximum rate of metabolism. Nasal extraction efficiencies were predicted to be greatest in the rat, followed by the monkey, and then the human. For both nasal and oral breathing modes in humans, higher uptake rates were predicted for lower tracheobronchial tissues than either the rat or monkey. These extended airway models provide a unique foundation for comparing material transport and site-specific tissue uptake across a significantly greater range of conducting airways in the rat, monkey, and human than prior CFD models.     

Inhalation dosimetry of hexamethylene diisocyanate vapor in the rat and human respiratory tracts

Abstract

Hexamethylene diisocyanate (HDI) is a reactive chemical used in the commercial production of polyurethanes. Toxic effects in rodents exposed to HDI vapor primarily occur in the nasal passages, yet some individuals exposed occupationally to concentrations exceeding current regulatory limits may experience temporary reduction in lung function and asthma-like symptoms. Knowledge of interspecies differences in respiratory tract dosimetry of inhaled HDI would improve our understanding of human health risks to this compound. HDI uptake was measured in the upper respiratory tract of anesthetized Fischer-344 rats. Nasal uptake of HDI was >90% in rats at unidirectional flow rates of 150 and 300?ml/min and a target air concentration of 200?ppb. Uptake data was used to calibrate nasal and lung dosimetry models of HDI absorption in rats and humans. Computational fluid dynamics (CFD) models of the nasal passages were used to simulate inspiratory airflow and HDI absorption. Transport of HDI through lung airways was simulated using convection-diffusion based mass transport models. HDI nasal uptake of 90% and 78% was predicted using the rat and human nasal CFD models, respectively. Total respiratory tract uptake was estimated to be 99% in rats and 97% in humans under nasal breathing. Predicted human respiratory uptake decreased to 87% under oral breathing conditions. Absorption rates of inhaled HDI in human lung airways were estimated to be higher than the rat due to lower uptake in head airways. Model predictions demonstrated significant penetration of HDI to human bronchial airways, although absorption rates were sensitive to breathing style.     

Studies of inspiratory airflow patterns in the nasal passages of the F344 rat and rhesus monkey using nasal molds: relevance to formaldehyde toxicity

For highly water soluble and reactive gases, such as formaldehyde, the reported distribution of nasal lesions in rats and rhesus monkeys following inhalation exposure may be attributable, at least in part, to regional gas uptake patterns that are a consequence of nasal airflow characteristics. Inspiratory nasal airflow was studied at flow rates across the physiologic range using a unidirectional dynamically similar water-dye siphon system in clear acrylic molds of the nasal airways of F344 rats and rhesus monkeys. In both species there were complex and inspiratory flow streams, exhibiting regions of simple laminar, complex secondary (vortices, eddies, swirling), and turbulent flows, with only minor effects of the volumetric flow rates studied on these flow patterns. There was a precise association between points of dye intake at the nostril with complex but generally coherent streaklines throughout the nose, indicating the potential for sensitive dependence of nasal airflow on nostril geometry. On the basis of these studies, a classification for the major airways (meatuses) in the nasal passages of rats and rhesus monkeys was proposed. The spiral shape of the anterior nasal airway of the rat was considered to play an important role in local mixing of inspired airstreams. In the rhesus monkey, the complex geometry of the nasal vestibule contributed to the formation of secondary flows and turbulence in the anterior nose, which represents a potentially important difference between rheusus monkeys and humans. There was a good correlation between routes of flow, regional secondary flows, turbulence, and impaction of airstreams on the airway wall, with the reported distribution of formaldehyde-induced nasal lesions in rats and rhesus monkeys. These studies support the proposal that nasal airflow patterns play an important role in the distribution of lesions induced by formaldehyde.     

A distributed-parameter model for formaldehyde uptake and disposition in the rat nasal lining

DNA-protein cross-links (DPX) serve as a dosimeter for inhaled formaldehyde and are associated with tumor induction in rat nasal passages after chronic exposure to 6 ppm and above. To determine the role of epithelium-specific morphometry in formaldehyde-induced patterns of injury, we developed a mathematical model that links airflow-driven formaldehyde uptake with DPX formation in regions of the rat nose with high and low tumor incidence. A three-dimensional, anatomically accurate computational fluid dynamics model of rat nasal airflow and inhaled gas uptake was integrated with a physiologically based mathematical model incorporating tissue thickness, formaldehyde diffusion, its removal by enzymatic and nonenzymatic processes, and DNA distribution in the nasal mucosa to predict DPX formation. The model implicitly incorporates the reversible conversion of formaldehyde to methylene glycol. Where possible, parameter values were taken from the literature or estimated using published correlations. The Michaelis-Menten kinetic constants Vmax and Km, as well as a first-order constant for formaldehyde removal, were left as fitted parameters. The resultant model fit to the experimentally measured DPX in the high- and low-tumor-incidence regions of the rat nasal passages was very good. Sensitivity analysis indicates that among the fitted parameters, model fits are most sensitive to Vmax and that predictions were sensitive to changes in tissue thickness when all other parameters are held constant. The model structure facilitates extrapolation to primates and humans and application to other soluble, reactive gases.     

A Distributed-Parameter Model for Formaldehyde Uptake and Disposition in the Rat Nasal Lining

DNA–protein cross-links (DPX) serve as a dosimeter for inhaled formaldehyde and are associated with tumor induction in rat nasal passages after chronic exposure to 6 ppm and above. To determine the role of epithelium-specific morphometry in formaldehyde-induced patterns of injury, we developed a mathematical model that links airflow-driven formaldehyde uptake with DPX formation in regions of the rat nose with high and low tumor incidence. A three-dimensional, anatomically accurate computational fluid dynamics model of rat nasal airflow and inhaled gas uptake was integrated with a physiologically based mathematical model incorporating tissue thickness, formaldehyde diffusion, its removal by enzymatic and nonenzymatic processes, and DNA distribution in the nasal mucosa to predict DPX formation. The model implicitly incorporates the reversible conversion of formaldehyde to methylene glycol. Where possible, parameter values were taken from the literature or estimated using published correlations. The Michaelis–Menten kinetic constants V_max and K_m, as well as a first-order constant for formaldehyde removal, were left as fitted parameters. The resultant model fit to the experimentally measured DPX in the high- and low-tumor-incidence regions of the rat nasal passages was very good. Sensitivity analysis indicates that among the fitted parameters, model fits are most sensitive to V_max and that predictions were sensitive to changes in tissue thickness when all other parameters are held constant. The model structure facilitates extrapolation to primates and humans and application to other soluble, reactive gases.     

Computational fluid dynamics simulations of inhaled nano- and microparticle deposition in the rhesus monkey nasal passages

Abstract

Anatomically accurate computational fluid dynamics (CFD) models of the nasal passages of an infant (6 months old, 1.3?kg) and adult (7 years old, 11.9?kg) rhesus monkey were used to predict nasal deposition of inhaled nano- and microparticles. Steady-state, inspiratory airflow simulations were conducted at flow rates equal to 100, 200 and 300% of the estimated minute volume for resting breathing in each model. Particle transport and deposition simulations were conducted using the Lagrangian method to track the motion of inhaled particles. Nasal deposition fractions were higher in the infant model than the adult model at equivalent physiologic flow rates. Deposition curves collapsed when differences in nasal geometry were accounted for by plotting microparticle deposition versus the Stokes number and nanoparticle deposition as a function of the Schmidt number and diffusion parameter. Particle deposition was also quantified on major nasal epithelial types. Maximum olfactory deposition ranged from 5 to 14% for 1–2?nm particles in the adult and infant models, depending on flow rate. For these particle sizes, deposition on respiratory/transitional epithelia ranged from 40 to 50%. Increased deposition was also predicted for olfactory and respiratory/transitional epithelia for particle sizes >5?µm in the infant model and >8?µm in the adult model. Semi-empirical curves were developed based on the CFD simulation results to allow for simplified calculations of age-based deposition in the rhesus monkey nasal passages that can be implemented into lung dosimetry models.     

Application of physiological computational fluid dynamics models to predict interspecies nasal dosimetry of inhaled acrolein

Acrolein is a highly soluble and reactive aldehyde and is a potent upper-respiratory-tract irritant. Acrolein-induced nasal lesions in rodents include olfactory epithelial atrophy and inflammation, epithelial hyperplasia, and squamous metaplasia of the respiratory epithelium. Nasal uptake of inhaled acrolein in rats is moderate to high, and depends on inspiratory flow rate, exposure duration, and concentration. In this study, anatomically accurate three-dimensional computational fluid dynamics (CFD) models were used to simulate steady-state inspiratory airflow and to quantitatively predict acrolein tissue dose in rat and human nasal passages. A multilayered epithelial structure was included in the CFD models to incorporate clearance of inhaled acrolein by diffusion, blood flow, and first-order and saturable metabolic pathways. Kinetic parameters for these pathways were initially estimated by fitting a pharmacokinetic model with a similar epithelial structure to time-averaged acrolein nasal extraction data and were then further adjusted using the CFD model. Predicted air:tissue flux from the rat nasal CFD model compared well with the distribution of acrolein-induced nasal lesions from a subchronic acrolein inhalation study. These correlations were used to estimate a tissue dose-based no-observed-adverse-effect level (NOAEL) for inhaled acrolein. A human nasal CFD model was used to extrapolate effects in laboratory animals to human exposure conditions on the basis of localized tissue dose and tissue responses. Assuming that equivalent tissue dose will induce similar effects across species, a NOAEL human equivalent concentration for inhaled acrolein was estimated to be 8 ppb.     

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.     

Incorporation of tissue reaction kinetics in a computational fluid dynamics model for nasal extraction of inhaled hydrogen sulfide in rats.

Rodents exposed to hydrogen sulfide (H2S) develop olfactory neuronal loss. This lesion has been used by the risk assessment community to develop occupational and environmental exposure standards. A correlation between lesion locations and areas of high H2S flux to airway walls has been previously demonstrated, but a quantitative dose assessment is needed to extrapolate dose at lesion sites to humans. In this study, nasal extraction (NE) of 10, 80, and 200 ppm H2S was measured in the isolated upper respiratory tract of anesthetized rats under constant unidirectional inspiratory flow rates of 75, 150, and 300 ml/min. NE was dependent on inspired H2S concentration and air flow rate: increased NE was observed when H2S exposure concentrations or inspiratory air flow rates were low. An anatomically accurate, three-dimensional computational fluid dynamics (CFD) model of rat nasal passages was used to predict NE of inhaled H2S. To account for the observed dependence of NE on H2S exposure concentration, the boundary condition used at airway walls incorporated first-order and saturable kinetics in nasal tissue to govern mass flux at the air:tissue interface. Since the kinetic parameters cannot be obtained using the CFD model, they were estimated independently by fitting a well-mixed, two-compartment pharmacokinetic (PK) model to the NE data. Predicted extraction values using this PK-motivated CFD approach were in good agreement with the experimental measurements. The CFD model provides estimates of localized H2S flux to airway walls and can be used to calibrate lesion sites by dose.     

Application of magnetic resonance (MR) imaging for the development and validation of computational fluid dynamic (CFD) models of the rat respiratory system

Computational fluid dynamic (CFD) models of the respiratory system provide a quantitative basis for extrapolating the localized dose of inhaled materials and improving human health risk assessments based upon inhalation studies conducted in animals. Nevertheless, model development and validation have historically been tedious and time-consuming tasks. In recognition of this, we previously reported on the use of proton (1H) magnetic resonance (MR) imaging for visualizing nasal-sinus passages in the rat, and for speeding computational mesh generation. Here, the generation and refinement of meshes for rat nasal airways are described in more detail and simulated airflows are presented. To extend the CFD models to the complete respiratory tract, three-dimensional (3D) 1H MR imaging of rat pulmonary casts was also utilized to construct pulmonary airway meshes using procedures developed for the nasal airways. Furthermore, the feasibility of validating CFD predictions with MR was tested by imaging hyperpolarized 3He gas at physiological flow rates in a straight pipe with a diameter comparable to the rat trachea. Results from these diverse studies highlight the potential utility of MR imaging not only for speeding CFD development but also possibly for model validation.     

Application of Magnetic Resonance (MR) Imaging for the Development and Validation of Computational Fluid Dynamic (CFD) Models of the Rat Respiratory System

Computational fluid dynamic (CFD) models of the respiratory system provide a quantitative basis for extrapolating the localized dose of inhaled materials and improving human health risk assessments based upon inhalation studies conducted in animals. Nevertheless, model development and validation have historically been tedious and time-consuming tasks. In recognition of this, we previously reported on the use of proton (¹H) magnetic resonance (MR) imaging for visualizing nasal-sinus passages in the rat, and for speeding computational mesh generation. Here, the generation and refinement of meshes for rat nasal airways are described in more detail and simulated airflows are presented. To extend the CFD models to the complete respiratory tract, three-dimensional (3D) ¹H MR imaging of rat pulmonary casts was also utilized to construct pulmonary airway meshes using procedures developed for the nasal airways. Furthermore, the feasibility of validating CFD predictions with MR was tested by imaging hyperpolarized ³He gas at physiological flow rates in a straight pipe with a diameter comparable to the rat trachea. Results from these diverse studies highlight the potential utility of MR imaging not only for speeding CFD development but also possibly for model validation.     

Human respiratory tract cancer risks of inhaled formaldehyde: dose-response predictions derived from biologically-motivated computational modeling of a combined rodent and human dataset.

Formaldehyde inhalation at 6 ppm and above causes nasal squamous cell carcinoma (SCC) in F344 rats. The quantitative implications of the rat tumors for human cancer risk are of interest, since epidemiological studies have provided only equivocal evidence that formaldehyde is a human carcinogen. Conolly et al. (Toxicol. Sci. 75, 432-447, 2003) analyzed the rat tumor dose-response assuming that both DNA-reactive and cytotoxic effects of formaldehyde contribute to SCC development. The key elements of their approach were: (1) use of a three-dimensional computer reconstruction of the rat nasal passages and computational fluid dynamics (CFD) modeling to predict regional dosimetry of formaldehyde; (2) association of the flux of formaldehyde into the nasal mucosa, as predicted by the CFD model, with formation of DNA-protein cross-links (DPX) and with cytolethality/regenerative cellular proliferation (CRCP); and (3) use of a two-stage clonal growth model to link DPX and CRCP with tumor formation. With this structure, the prediction of the tumor dose response was extremely sensitive to cell kinetics. The raw dose-response data for CRCP are J-shaped, and use of these data led to a predicted J-shaped dose response for tumors, notwithstanding a concurrent low-dose-linear, directly mutagenic effect of formaldehyde mediated by DPX. In the present work the modeling approach used by Conolly et al. (ibid.) was extended to humans. Regional dosimetry predictions for the entire respiratory tract were obtained by merging a three-dimensional CFD model for the human nose with a one-dimensional typical path model for the lower respiratory tract. In other respects, the human model was structurally identical to the rat model. The predicted human dose response for DPX was obtained by scale-up of a computational model for DPX calibrated against rat and rhesus monkey data. The rat dose response for CRCP was used "as is" for the human model, since no preferable alternative was identified. Three sets of baseline parameter values for the human clonal growth model were obtained through separate calibrations against respiratory tract cancer incidence data for nonsmokers, smokers, and a mixed population of nonsmokers and smokers, respectively. Additional risks of respiratory tract cancer were predicted to be negative up to about one ppm for all three cases when the raw CRCP data from the rat were used. When a hockey-stick-shaped model was fit to the rat CRCP data and used in place of the raw data, positive maximum likelihood estimates (MLE) of additional risk were obtained. These MLE estimates were lower, for some comparisons by as much as 1,000-fold, than MLE estimates from previous cancer dose-response assessments for formaldehyde. Breathing rate variations associated with different physical activity levels did not make large changes in predicted additional risks. In summary, this analysis of the human implications of the rat SCC data indicates that (1) cancer risks associated with inhaled formaldehyde are de minimis (10(-6) or less) at relevant human exposure levels, and (2) protection from the noncancer effects of formaldehyde should be sufficient to protect from its potential carcinogenic effects.     

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.     

A computational fluid dynamics approach to assess interhuman variability in hydrogen sulfide nasal dosimetry

Human exposure to hydrogen sulfide (H₂S) gas occurs from natural and industrial sources and can result in dose-related neurological, respiratory, and cardiovascular effects. Olfactory neuronal loss in H₂S-exposed rats has been used to develop occupational and environmental exposure limits. Using nasal computational fluid dynamics (CFD) models, a correlation was found between wall mass flux and olfactory neuronal loss in rodents, suggesting an influence of airflow patterns on lesion locations that may affect interspecies extrapolation of inhaled dose. Human nasal anatomy varies considerably within a population, potentially affecting airflow patterns and dosimetry of inhaled gases. This study investigates interhuman variability of H₂S nasal dosimetry using anatomically accurate CFD models of the nasal passages of five adults and two children generated from magnetic resonance imaging (MRI) or computed tomography (CT) scan data. Using allometrically equivalent breathing rates, steady-state inspiratory airflow and H₂S uptake were simulated. Approximate locations of olfactory epithelium were mapped in each model to compare air:tissue flux in the olfactory region among individuals. The fraction of total airflow to the olfactory region ranged from 2% to 16%. Despite this wide range in olfactory airflow, H₂S dosimetry in the olfactory region was predicted to be similar among individuals. Differences in the 99th percentile and average flux values were <1.2-fold at inhaled concentrations of 1, 5, and 10 ppm. These preliminary results suggest that differences in nasal anatomy and ventilation among adults and children do not have a significant effect on H₂S dosimetry in the olfactory region.     

A model for the uptake of inhaled vapors in the nose of the dog during cyclic breathing

A model was developed to simulate the uptake of inhaled vapors in the nasal airway of the Beagle dog during cyclic breathing. Input data to the model were morphological and physiological data for the dog, and physiocochemical data for the vapors. The model simulates the nasal airway as a slit-like duct, where air passes between the two airway walls in an ideal plug flow. The thickness of the walls corresponds to the distance between the air interface and the average position where vapor molecules are removed into the capillary blood. All resistance to radial mass transfer is assumed to arise on the liquid side in the diffusion of vapors through the air/blood tissue barrier and in transport by the blood. The model results agreed reasonably well with experimental data. The nasal absorption of vapors on inhalation increased from 1% for a compound with a blood/air partition coefficient of 1 to around 95% uptake for a compound with a partition coefficient of 2000. Desorption from the nasal tissues on exhalation increased from approximately 1% to approximately 30% over the same range of partition coefficients. The nasal uptake over one complete breathing cycle, defined as absorption on inhalation minus desorption on exhalation, was almost zero for low partition coefficient compounds and plateaued at around 65% for high partition coefficient compounds. The model indicates that diffusional resistance and inertia of the nasal tissues result in temporary storage of absorbed vapors upon inhalation, followed by desorption of vapors back into the airstream upon exhalation. An important consequence of this phenomenon is a shift in exposure to inhaled vapors from the lungs to the nasal airway during cyclic flow compared with predictions from experiments and models based on unidirectional flow.     

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₃], nitrogen dioxide [NO₂], and sulfur dioxide [SO₂]) 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₃ 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.