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

Detailed comparative analysis of environmental microparticle deposition characteristics between human and monkey nasal cavities using a surface mapping technique

Inhaled particulate matter is associated with nasal diseases such as allergic rhinitis, rhinosinusitis and neural disorders. Its health risks on humans are usually evaluated by measurements on monkeys as they share close phylogenetic relationship. However, the reliability of cross-species toxicological extrapolation is in doubt due to physiological and anatomical variations, which greatly undermine the reliability of these expensive human surrogate models. This study numerically investigated in-depth microparticle transport and deposition characteristics on human and monkey (Macaca fuscata) nasal cavities that were reconstructed from CT-images. Deposition characteristics of 1-30μm particles were investigated under resting and active breathing conditions. Similar trends were observed for total deposition efficiencies and a single correlation using Stokes Number was fitted for both species and both breathing conditions, which is convenient for monkey-human extrapolation. Regional deposition patterns were carefully compared using the surface mapping technique. Deposition patterns of low, medium and high inertial particles, classified based on their total deposition efficiencies, were further analyzed in the 3D view and the mapped 2D view, which allows locating particle depositions on specific nasal regions. According to the particle intensity contours and regional deposition profiles, the major differences were observed at the vestibule and the floor of the nasal cavity, where higher deposition intensities of medium and high inertial particles were shown in the monkey case than the human case. Comparisons of airflow streamlines indicated that the cross-species variations of microparticle deposition patterns are mainly contributed by two factors. First, the more oblique directions of monkey nostrils result in a sharper airflow turn in the vestibule region. Second, the monkey's relatively narrower nasal valves lead to higher impaction of medium and high inertial particles on the nasal cavity floor. The methods and findings in this study would contribute to an improved cross-species toxicological extrapolation between human and monkey nasal cavities.

Incorporation of rapid association/dissociation processes in tissues into the monkey and human physiologically based pharmacokinetic models for manganese

In earlier physiologically based pharmacokinetic (PBPK) models for manganese (Mn), the kinetics of transport of Mn into and out of tissues were primarily driven by slow rates of association and dissociation of Mn with tissue binding sites. However, Mn is known to show rapidly reversible binding in tissues. An updated Mn model for primates, following similar work with rats, was developed that included rapid association/dissociation processes with tissue Mn-binding sites, accumulation of free Mn in tissues after saturation of these Mn-binding sites and rapid rates of entry into tissues. This alternative structure successfully described Mn kinetics in tissues in monkeys exposed to Mn via various routes including oral, inhalation, and intraperitoneal, subcutaneous, or intravenous injection and whole-body kinetics and tissue levels in humans. An important contribution of this effort is showing that the extension of the rate constants for binding and cellular uptake established in the monkey were also able to describe kinetic data from humans. With a consistent model structure for monkeys and humans, there is less need to rely on cadaver data and whole-body tracer studies alone to calibrate a human model. The increased biological relevance of the Mn model structure and parameters provides greater confidence in applying the Mn PBPK models to risk assessment. This model is also well-suited to explicitly incorporate emerging information on the role of transporters in tissue disposition, intestinal uptake, and hepatobiliary excretion of Mn.

Analyses on the influence of normal nasal morphological variations on odorant transport to the olfactory cleft

OBJECTIVE

Olfaction requires a combination of sensorineural components and conductive components, but conductive mechanisms have not typically received much attention. This study investigates the role of normal nasal vestibule morphological variations in ten healthy subjects on odorant flux in the olfactory cleft.

MATERIALS AND METHODS

Computed tomography images were used to create subject-specific nasal models. Each subject's unilateral nasal cavity was classified according to its nasal vestibule shape as Standard or Notched. Inspiratory airflow simulations were performed at 15 L/min, simulating resting inspiration using computational fluid dynamics modeling. Odorant transport simulations for three odorants (limonene, 2,4-dinitrotoluene, and acetaldehyde) were then performed at concentrations of 200 ppm for limonene and acetaldehyde, and 0.2 ppm for dinitrotoluene. Olfactory cleft odorant flux was computed for each simulation.

RESULTS AND DISCUSSION AND CONCLUSION

Simulated results showed airflow in the olfactory cleft was greater in the Standard phenotype compared to the Notched phenotype. For Standard, median airflow was greatest in the anterior region (0.5006 L/min) and lowest in the posterior region (0.1009 L/min). Median airflow in Notched was greatest in the medial region (0.3267 L/min) and lowest in the posterior region (0.0756 L/min). Median olfactory odorant flux for acetaldehyde and limonene was greater in Standard (Acetaldehyde: Standard = 140.45 pg/cm2-s; Notched = 122.20 pg/cm2-s. Limonene: Standard = 0.67 pg/cm2-s; Notched = 0.65 pg/cm2-s). Median dinitrotoluene flux was greater in Notched (Standard = 2.86 × 10-4pg/cm2-s; Notched = 4.29 × 10-4 pg/cm2-s). The impact of nasal vestibule morphological variations on odorant flux at the olfactory cleft may have implications on individual differences in olfaction, which should be investigated further.

Integration of Engineered Delivery with the Pharmacokinetics of Medical Candidates via Physiology-Based Pharmacokinetics

Physiologically based pharmacokinetic (PBPK) modeling is a mechanistic computational model that can be used to predict a drug product's ADME (absorption, distribution, metabolism, and excretion) and pharmacokinetics (PK). In recent years, PBPK modeling and simulation has been used increasingly to address many biopharmaceutics and clinical pharmacology questions, such as the effect of formulations, intrinsic factors (age, organ dysfunction, etc.), and extrinsic factors (comedications, food) on the PK of an investigational drug product. In this chapter, we will briefly introduce various PBPK models for ADME prediction and general procedures for PBPK modeling and simulations. The readers are encouraged to read updated literature on new applications of PBPK modeling and simulation which is still an emerging area in pharmaceutical development.

Direct nose to the brain nanomedicine delivery presents a formidable challenge

This advanced review describes the anatomical and physiological barriers and mechanisms impacting nanomedicine translocation from the nasal cavity directly to the brain. There are significant physiological and anatomical differences in the nasal cavity, olfactory area, and airflow reaching the olfactory epithelium between humans and experimentally studied species that should be considered when extrapolating experimental results to humans. Mucus, transporters, and tight junction proteins present barriers to material translocation across the olfactory epithelium. Uptake of nanoparticles through the olfactory mucosa and translocation to the brain can be intracellular via cranial nerves (intraneuronal) or other cells of the olfactory epithelium, or extracellular along cranial nerve pathways (perineural) and surrounding blood vessels (perivascular, the glymphatic system). Transport rates vary greatly among the nose to brain pathways. Nanomedicine physicochemical properties (size, surface charge, surface coating, and particle stability) can affect uptake efficiency, which is usually less than 5%. Incorporation of therapeutic agents in nanoparticles has been shown to produce pharmacokinetic and pharmacodynamic benefits. Assessment of adverse effects has included olfactory mucosa toxicity, ciliotoxicity, and olfactory bulb and brain neurotoxicity. The results have generally suggested the investigated nanomedicines do not present significant toxicity. Research needs to advance the understanding of nanomedicine translocation and its drug cargo after intranasal administration is presented. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Neurological Disease Therapeutic Approaches and Drug Discovery > Emerging Technologies Toxicology and Regulatory Issues in Nanomedicine > Toxicology of Nanomaterials.

Uniqueness of inspiratory airflow patterns in a realistic rat nasal cavity

In this study, we present a detailed flow analysis using an anatomically accurate rat nasal cavity model, in which the anatomy and physiology of the nasal airway was thoroughly examined. Special efforts were given to the swirling flow structures in the nasal vestibule (anterior section of the nose, lined by squamous epithelium), fractional flow patterns in the olfactory (posterior superior section of the rat nose, lined by olfactory epithelium), and a designated method to precisely quantify flow apportionment in the olfactory region was developed. Results revealed distinct inspiratory flow patterns in the anterior vestibule region, where the accelerated airflow undergoes two sharp turns as traveling through the tortuous airway, making a route in a shape of 8. Besides this, exceptionally large flow apportionment was observed at the interface of the olfactory recess, which can be as much as 15 times greater than that in the human nose. The thorough understanding of the airflow dynamics in the rat nasal cavity is necessary to avoid potential misinterpretation of rat-derived inhalation toxicity results. Research findings are expected to play a fundamental role in developing unbiased rat to human interspecies data extrapolation schemes.

New Approach Methodology for Assessing Inhalation Risks of a Contact Respiratory Cytotoxicant: Computational Fluid Dynamics-Based Aerosol Dosimetry Modeling for Cross-Species and In Vitro Comparisons

Regulatory agencies are considering alternative approaches to assessing inhalation toxicity that utilizes in vitro studies with human cells and in silico modeling in lieu of additional animal studies. In support of this goal, computational fluid-particle dynamics models were developed to estimate site-specific deposition of inhaled aerosols containing the fungicide, chlorothalonil, in the rat and human for comparisons to prior rat inhalation studies and new human in vitro studies. Under bioassay conditions, the deposition was predicted to be greatest at the front of the rat nose followed by the anterior transitional epithelium and larynx corresponding to regions most sensitive to local contact irritation and cytotoxicity. For humans, simulations of aerosol deposition covering potential occupational or residential exposures (1-50 µm diameter) were conducted using nasal and oral breathing. Aerosols in the 1-5 µm range readily penetrated the deep region of the human lung following both oral and nasal breathing. Under actual use conditions (aerosol formulations >10 µm), the majority of deposited doses were in the upper conducting airways. Beyond the nose or mouth, the greatest deposition in the pharynx, larynx, trachea, and bronchi was predicted for aerosols in the 10-20 µm size range. Only small amounts of aerosols >20 µm penetrated past the pharyngeal region. Using the ICRP clearance model, local retained tissue dose metrics including maximal concentrations and areas under the curve were calculated for each airway region following repeated occupational exposures. These results are directly comparable with benchmark doses from in vitro toxicity studies in human cells leading to estimated human equivalent concentrations that reduce the reliance on animals for risk assessments.

Recent Advances in Probabilistic Dose-Response Assessment to Inform Risk-Based Decision Making

Paradoxically, risk assessments for the majority of chemicals lack any quantitative characterization as to the likelihood, incidence, or severity of the risks involved. The relatively few cases where "risk" is truly quantified are based on either epidemiologic data or extrapolation of experimental animal cancer bioassay data. The paucity of chemicals and health endpoints for which such data are available severely limits the ability of decisionmakers to account for the impacts of chemical exposures on human health. The development by the World Health Organization International Programme on Chemical Safety (WHO/IPCS) in 2014 of a comprehensive framework for probabilistic dose-response assessment has opened the door to a myriad of potential advances to better support decision making. Building on the pioneering work of Evans, Hattis, and Slob from the 1990s, the WHO/IPCS framework provides both a firm conceptual foundation as well as practical implementation tools to simultaneously assess uncertainty, variability, and severity of effect as a function of exposure. Moreover, such approaches do not depend on the availability of epidemiologic data, nor are they limited to cancer endpoints. Recent work has demonstrated the broad feasibility of such approaches in order to estimate the functional relationship between exposure level and the incidence or severity of health effects. While challenges remain, such as better characterization of the relationship between endpoints observed in experimental animal or in vitro studies and human health effects, the WHO/IPCS framework provides a strong basis for expanding the breadth of risk management decision contexts supported by chemical risk assessment.

Application of a unified probabilistic framework to the dose-response assessment of acrolein

BACKGROUND

In quantitative chemical risk assessment, a reference value is an estimate of an exposure to a chemical that is "likely to be without appreciable risk." Because current "deterministic" approaches do not quantitatively characterize the likelihood or severity of harm, the National Academies has recommended using reference values derived from a risk-specific dose that are treated as random variables, with probability distributions characterizing uncertainty and variability.

OBJECTIVES

In order to build familiarity and address issues needed for routine and standardized derivation of probabilistic risk-specific dose distributions, a case example applying the unified probabilistic framework presented in Chiu and Slob (2015) is developed for acrolein. This case study is based on an updated systematic evidence map of literature (Keshava et al., 2020) identifying nasal lesions reported in Dorman et al. (2008) as the most appropriate endpoint and study for reference value derivation.

METHODS

The probability distribution was calculated for the risk-specific dose, which in this implementation of the approach was calculated for the dose at which 1% of the human population is estimated to experience minimal lesions, and a probabilistic reference value was computed as the 5th percentile of this distribution. A deterministic reference value was also derived for comparison, and a sensitivity analysis of the probabilistic reference value was conducted investigating alternative assumptions for the point of departure type and exposure duration.

RESULTS

The probabilistic reference value of 6 × 10-4 mg/m3 was slightly lower than the deterministic reference value of 8 × 10-4 mg/m3, and the risk-specific dose distribution had an uncertainty spanning a factor of 137 (95th-5th percentile ratio). Sensitivity analysis yielded slightly higher probabilistic reference values ranging between 9 × 10-4 mg/m3 and 2 × 10-3 mg/m3.

CONCLUSIONS

Using a probabilistic approach for deriving a reference value allows quantitative characterization of the severity, incidence, and uncertainty of effects at a given dose. The results can be used to inform risk management decisions and improve risk communication.

Application of systematic evidence mapping to assess the impact of new research when updating health reference values: A case example using acrolein

BACKGROUND

The environmental health community needs transparent, methodologically rigorous, and rapid approaches for updating human health risk assessments. These assessments often contain reference values for cancer and/or noncancer effects. Increasingly, the use of systematic review methods are preferred when developing these assessments. Systematic evidence maps are a type of analysis that has the potential to be very helpful in the update process, especially when combined with machine-learning software advances designed to expedite the process of conducting a review.

OBJECTIVES

To evaluate the applicability of evidence mapping to determine whether new evidence is likely to result in a change to an existing health reference value, using inhalation exposure to the air pollutant acrolein as a case example.

METHODS

New literature published since the 2008 California Environmental Protection Agency's Office of Environmental Health Hazard Assessment (OEHHA) Reference Exposure Level (REL) for acrolein was assessed. Systematic review methods were used to search the literature and screening included the use of machine-learning software. The Populations, Exposures, Comparators and Outcomes (PECO) criteria were kept broad to identify studies that characterized acute and chronic exposure and could be informative for hazard characterization. Studies that met the PECO criteria after full-text review were briefly summarized before their suitability for chronic point of departure (POD) derivation and calculation of a reference value was considered. Studies considered potentially suitable underwent a targeted evaluation to determine their suitability for use in dose-response analysis.

RESULTS

Over 15,000 studies were identified from scientific databases. Both machine-learning and manual screening processes were used to identify 60 studies considered PECO-relevant after full-text review. Most of these PECO-relevant studies were short-term exposure animal studies (acute or less than 1 month of exposure) and considered less suitable for deriving a chronic reference value when compared to the subchronic study in rats used in the 2008 OEHHA assessment. Thirteen epidemiological studies were identified but had limitations in the exposure assessment that made them less suitable for dose-response compared to the subchronic rat study. Among the 13 studies, there were four controlled trial studies that have the potential to be informative for future acute reference value derivation. Thus, the 2008 subchronic rat study used by OEHHA appears to still be the most appropriate study for chronic reference value derivation. In addition, advances in dosimetric modeling for gases, including new evidence pertinent to acrolein, could be considered when updating existing acrolein toxicity values.

CONCLUSIONS

Evidence mapping is a very useful tool to assess the need for updating an assessment based on understanding the potential impact of new studies on revising an existing health reference value. In this case example, the focus was to identify studies suitable for chronic exposure dose-response analysis, while also identifying studies that may be important to consider for acute exposure scenarios, hazard identification, or for future research. This allows the evidence map to be a useful resource for a range of decision-making contexts. Specialized systematic review software increased the efficiency of the process in terms of human resources and time to conduct the analysis.

A numerical investigation of the potential effects of e-cigarette smoking on local tissue dosimetry and the deterioration of indoor air quality

Electronic (e)-cigarette smoking is considered to be less harmful than traditional tobacco smoking because of the lack of a combustion process. However, e-cigarettes have the potential to release harmful chemicals depending on the constituents of the vapor. To date, there has been significant evidence on the adverse health effects of e-cigarette usage. However, what is less known are the impacts of the chemicals contained in exhaled air from an e-cigarette smoker on indoor air quality, the second-hand passive smoking of residents, and the toxicity of the exhaled air. In this study, we develop a comprehensive numerical model and computer-simulated person to investigate the potential effects of e-cigarette smoking on local tissue dosimetry and the deterioration of indoor air quality. We also conducted demonstrative numerical analyses for first-hand and second-hand e-cigarette smoking in an indoor environment. To investigate local tissue dosimetry, we used newly developed physiologically based pharmacokinetic/toxicokinetic models that reproduce inhalation exposure by way of the respiratory tract and dermal exposure through the human skin surface. These models were integrated into the computer-simulated person. Our numerical simulation results quantitatively demonstrated the potential impacts of e-cigarette smoking in enclosed spaces on indoor air quality.

Use of Nasal Pathology in the Derivation of Inhalation Toxicity Values for Hydrogen Sulfide

Nasal pathology can play an important role in the risk assessment process. For example, olfactory neuron loss (ONL) is one of the most sensitive end points seen in subchronic rodent hydrogen sulfide (H₂S) studies and has been used by several agencies to derive health-protective toxicity values. Alternative methods that rely on computational fluid dynamics (CFD) models to account for the influence of airflow on H₂S-induced ONL have been proposed. The use of CFD models result in toxicity values that are less conservative than those obtained using more traditional methods. These alternative approaches rely on anatomy-based CFD models. Model predictions of H₂S delivery (flux) to the olfactory mucosal wall are highly correlated with ONL in rodents. Three major areas of focus for this review include a brief description of nasal anatomy, H₂S-induced ONL in rodents, derivation of a chronic inhalation reference concentration for H₂S, and the use of CFD models to derive alternative toxicity values for this gas.

In silico approaches to respiratory nasal flows: A review

The engineering discipline of in silico fluid dynamics delivers quantitative information on airflow behaviour in the nasal regions with unprecedented detail, often beyond the reach of traditional experiments. The ability to provide visualisation and analysis of flow properties such as velocity and pressure fields, as well as wall shear stress, dynamically during the respiratory cycle may give significant insight to clinicians. Yet, there remains ongoing challenges to advance the state-of-the-art further, including for example the lack of comprehensive CFD modelling on varied cohorts of patients. The present article embodies a review of previous and current in silico approaches to simulating nasal airflows. The review discusses specific modelling techniques required to accommodate physiologically- and clinically-relevant findings. It also provides a critical summary of the reported results in the literature followed by an outlook on the challenges and topics anticipated to drive research into the future.

Correlation of regional deposition dosage for inhaled nanoparticles in human and rat olfactory

Background

Nose-to-brain transport of airborne ultrafine particles (UFPs) via the olfactory pathway has been verified as a possible route for particle translocation into the brain. The exact relationship between increased airborne toxicant exposure and neurological deterioration in the human central nervous system, is still unclear. However, the nasal olfactory is undoubtedly a critical junction where the time course and toxicant dose dependency might be inferred.

Method

Computational fluid-particle dynamics modeling of inhaled nanoparticles (1 to 100 nm) under low to moderate breathing conditions (5 to 14 L/min – human; and 0.14 to 0.40 L/min – rat) were performed in physiologically realistic human and rat nasal airways. The simulation emphasized olfactory deposition, and variations in airflow and particle flux caused by the inter-species airway geometry differences. Empirical equations were developed to predict regional deposition rates of inhaled nanoparticles on human and rat olfactory mucosa in sedentary breathing. Considering, breathing and geometric differences, quantified correlations between human and the rat olfactory deposition dose against a variety of metrics were proposed.

Results

Regional deposition of nanoparticles in human and the rat olfactory was extremely low, with the highest deposition (< 3.5 and 8.1%) occurring for high diffusivity particles of 1.5 nm and 5 nm, respectively. Due to significant filtering of extremely small particles (< 2 nm) by abrupt sharp turns at front of the rat nose, only small fractions of the inhaled nanoparticles (in this range) reached rat olfactory than that in human (1.25 to 45%); however, for larger sizes (> 3 nm), significantly higher percentage of the inhaled nanoparticles reached rat nasal olfactory than that in human (2 to 32 folds). Taking into account the physical and geometric features between human and rat, the total deposition rate (#/min) and deposition rate per unit surface area (#/min/mm²) were comparable for particles> 3 nm. However, when body mass was considered, the normalized deposition rate (#/min/kg) in the rat olfactory region exceeded that in the human. Nanoparticles < 1.5 nm were filtered out by rat anterior nasal cavity, and therefore deposition in human olfactory region exceeded that in the rat model.

Conclusion

Regional deposition dose of inhaled nanoparticles in a human and rat olfactory region was governed by particle size and the breathing rate. Interspecies correlation was determined by combining the effect of deposition dosage, physical\geometric features, and genetic differences. Developed empirical equations provided a tool to quantify inhaled nanoparticle dose in human and rat nasal olfactory regions, which lay the ground work for comprehensive interspecies correlation between the two species. Furthermore, this study contributes to the fields in toxicology, i.e., neurotoxicity evaluation and risk assessment of UFPs, in long-term and low-dose inhalation exposure scenarios.

Molecular Binding Contributes to Concentration Dependent Acrolein Deposition in Rat Upper Airways: CFD and Molecular Dynamics Analyses

Existing in vivo experiments show significantly decreased acrolein uptake in rats with increasing inhaled acrolein concentrations. Considering that high-polarity chemicals are prone to bond with each other, it is hypothesized that molecular binding between acrolein and water will contribute to the experimentally observed deposition decrease by decreasing the effective diffusivity. The objective of this study is to quantify the probability of molecular binding for acrolein, as well as its effects on acrolein deposition, using multiscale simulations. An image-based rat airway geometry was used to predict the transport and deposition of acrolein using the chemical species model. The low Reynolds number turbulence model was used to simulate the airflows. Molecular dynamic (MD) simulations were used to study the molecular binding of acrolein in different media and at different acrolein concentrations. MD results show that significant molecular binding can happen between acrolein and water molecules in human and rat airways. With 72 acrolein embedded in 800 water molecules, about 48% of acrolein compounds contain one hydrogen bond and 10% contain two hydrogen bonds, which agreed favorably with previous MD results. The percentage of hydrogen-bonded acrolein compounds is higher at higher acrolein concentrations or in a medium with higher polarity. Computational dosimetry results show that the size increase caused by the molecular binding reduces the effective diffusivity of acrolein and lowers the chemical deposition onto the airway surfaces. This result is consistent with the experimentally observed deposition decrease at higher concentrations. However, this size increase can only explain part of the concentration-dependent variation of the acrolein uptake and acts as a concurrent mechanism with the uptake-limiting tissue ration rate. Intermolecular interactions and associated variation in diffusivity should be considered in future dosimetry modeling of high-polarity chemicals such as acrolein.

Partitioning of dispersed nanoparticles in a realistic nasal passage for targeted drug delivery

The complex nasal structure poses obstacles for efficient nasal drug administration beyond the nasal valve, especially when targeting the olfactory region. This study numerically detailed the naturally inhaled nanoparticle transport process from the initial releasing locations to the final deposited sites using a realistic human nasal passage. Dispersed nanoparticles at different coronal cross-sections were partitioned into multiple groups according to their final deposited locations. Results showed inhaled nanoparticles are more likely to move along the septum. Olfactory deposited particles entered the nose through the inner superior corner of the nostril; the middle meatus deposited particles entered the nose through the top third of the nostril; the inferior deposited particles entered via the bottom floor regions of the nostril. Therefore, targeted nasal inhalation therapies that intentionally release therapeutic particles from these recognized regions at the nostril plane can considerably improve the resultant topical disposition doses. However, it remains challenging to completely prevent undesired particle depositions as particles coming from the same location may produce multiple-sites depositions due to partition overlapping. Nevertheless, the fraction of undesired particle deposition is anticipated to be reduced at a great extent compared to unplanned releasing approaches.

Detailed deposition analysis of inertial and diffusive particles in a rat nasal passage

Rats have been widely used as surrogates for evaluating the health effects of inhaled airborne particulate matter. To provide a thorough understanding of particle transport and deposition mechanisms in the rat nasal airway, this article presents a computational fluid dynamics (CFD) study of particle exposure in a realistic rat nasal passage under a resting flow condition. Particles covering a diameter range from 1 nm to 4 µm were passively released in front of the rat's breathing zone, and the Lagrangian particle tracking approach was used to calculate individual particle trajectories. Detailed particle deposition analysis shows the deposition of inertial particles >2 µm is high in the rat nasal vestibule and more than 70% of all inhaled inertial particles were trapped in this region. While for diffusive nanoparticles, the vestibule filtration effect is reduced, only less than 60% of inhaled nanoparticles were blocked by the anterior nasal structures. The particle exposure in the olfactory region only shows notable deposition for diffusive nanoparticles, which peaks at 9.4% for 5 nm particles. Despite the olfactory deposition remains at a low level, the ratio between the olfactory and the main passage is kept around 30-40% for 10-800 nm particles, which indicates a particle-size-independent distribution pattern in the main nasal passage and olfactory. This study provides a deep understanding of particles deposition features in a rat nasal passage, and the research findings can aid toxicologist in inter-species exposure-response extrapolation study.

Small Airway Absorption and Microdosimetry of Inhaled Corticosteroid Particles after Deposition

PURPOSE

To predict the cellular-level epithelial absorbed dose from deposited inhaled corticosteroid (ICS) particles in a model of an expanding and contracting small airway segment for different particle forms.

METHODS

A computational fluid dynamics (CFD)-based model of drug dissolution, absorption and clearance occurring in the surface liquid of a representative small airway generation (G13) was developed and used to evaluate epithelial dose for the same deposited drug mass of conventional microparticles, nanoaggregates and a true nanoaerosol. The ICS medications considered were budesonide (BD) and fluticasone propionate (FP). Within G13, total epithelial absorption efficiency (AE) and dose uniformity (microdosimetry) were evaluated.

RESULTS

Conventional microparticles resulted in very poor AE of FP (0.37%) and highly nonuniform epithelial absorption, such that <5% of cells received drug. Nanoaggregates improved AE of FP by a factor of 57-fold and improved dose delivery to reach approximately 40% of epithelial cells. True nanoaerosol resulted in near 100% AE for both drugs and more uniform drug delivery to all cells.

CONCLUSIONS

Current ICS therapies are absorbed by respiratory epithelial cells in a highly nonuniform manner that may partially explain poor clinical performance in the small airways. Both nanoaggregates and nanoaerosols can significantly improve ICS absorption efficiency and uniformity.

A tissue dose-based comparative exposure assessment of manganese using physiologically based pharmacokinetic modeling-The importance of homeostatic control for an essential metal

A physiologically-based pharmacokinetic (PBPK) model (Schroeter et al., 2011) was applied to simulate target tissue manganese (Mn) concentrations following occupational and environmental exposures. These estimates of target tissue Mn concentrations were compared to determine margins of safety (MOS) and to evaluate the biological relevance of applying safety factors to derive acceptable Mn air concentrations. Mn blood concentrations measured in occupational studies permitted verification of the human PBPK models, increasing confidence in the resulting estimates. Mn exposure was determined based on measured ambient air Mn concentrations and dietary data in Canada and the United States (US). Incorporating dietary and inhalation exposures into the models indicated that increases in target tissue concentrations above endogenous levels only begin to occur when humans are exposed to levels of Mn in ambient air (i.e. >10μg/m³) that are far higher than those currently measured in Canada or the US. A MOS greater than three orders of magnitude was observed, indicating that current Mn air concentrations are far below concentrations that would be required to produce the target tissue Mn concentrations associated with subclinical neurological effects. This application of PBPK modeling for an essential element clearly demonstrates that the conventional application of default factors to "convert" an occupational exposure to an equivalent continuous environmental exposure, followed by the application of safety factors, is not appropriate in the case of Mn. PBPK modeling demonstrates that the relationship between ambient Mn exposures and dose-to-target tissue is not linear due to normal tissue background levels and homeostatic controls.

Comparison of realistic and idealized breathing patterns in computational models of airflow and vapor dosimetry in the rodent upper respiratory tract

CONTEXT

Computational fluid dynamics (CFD) simulations of airflows coupled with physiologically based pharmacokinetic (PBPK) modeling of respiratory tissue doses of airborne materials have traditionally used either steady-state inhalation or a sinusoidal approximation of the breathing cycle for airflow simulations despite their differences from normal breathing patterns.

OBJECTIVE

Evaluate the impact of realistic breathing patterns, including sniffing, on predicted nasal tissue concentrations of a reactive vapor that targets the nose in rats as a case study.

MATERIALS AND METHODS

Whole-body plethysmography measurements from a free-breathing rat were used to produce profiles of normal breathing, sniffing and combinations of both as flow inputs to CFD/PBPK simulations of acetaldehyde exposure.

RESULTS

For the normal measured ventilation profile, modest reductions in time- and tissue depth-dependent areas under the curve (AUC) acetaldehyde concentrations were predicted in the wet squamous, respiratory and transitional epithelium along the main airflow path, while corresponding increases were predicted in the olfactory epithelium, especially the most distal regions of the ethmoid turbinates, versus the idealized profile. The higher amplitude/frequency sniffing profile produced greater AUC increases over the idealized profile in the olfactory epithelium, especially in the posterior region.

CONCLUSIONS

The differences in tissue AUCs at known lesion-forming regions for acetaldehyde between normal and idealized profiles were minimal, suggesting that sinusoidal profiles may be used for this chemical and exposure concentration. However, depending upon the chemical, exposure system and concentration and the time spent sniffing, the use of realistic breathing profiles, including sniffing, could become an important modulator for local tissue dose predictions.

Effects of Ambient Air Pollution Exposure on Olfaction: A Review

BACKGROUND

Olfactory dysfunction affects millions of people worldwide. This sensory impairment is associated with neurodegenerative disease and significantly decreased quality of life. Exposure to airborne pollutants has been implicated in olfactory decline, likely due to the anatomic susceptibility of the olfactory nerve to the environment. Historically, studies have focused on occupational exposures, but more recent studies have considered effects from exposure to ambient air pollutants.

OBJECTIVES

To examine all relevant human data evaluating a link between ambient pollution exposure and olfaction and to review supporting animal data in order to examine potential mechanisms for pollution-associated olfactory loss.

METHODS

We identified and reviewed relevant articles from 1950 to 2015 using PubMed and Web of Science and focusing on human epidemiologic and pathophysiologic studies. Animal studies were included only to support pertinent data on humans. We reviewed findings from these studies evaluating a relationship between environmental pollutant exposure and olfactory function.

RESULTS

We identified and reviewed 17 articles, with 1 additional article added from a bibliography search, for a total of 18 human studies. There is evidence in human epidemiologic and pathologic studies that increased exposure to ambient air pollutants is associated with olfactory dysfunction. However, most studies have used proxies for pollution exposure in small samples of convenience. Human pathologic studies, with supporting animal work, have also shown that air pollution can contact the olfactory epithelium, translocate to the olfactory bulb, and migrate to the olfactory cortex. Pollutants can deposit at each location, causing direct damage and disruption of tissue morphology or inducing local inflammation and cellular stress responses.

CONCLUSIONS

Ambient air pollution may impact human olfactory function. Additional studies are needed to examine air pollution-related olfactory impacts on the general population using measured pollution exposures and to link pollution exposure with olfactory dysfunction and related pathology. Citation: Ajmani GS, Suh HH, Pinto JM. 2016. Effects of ambient air pollution exposure on olfaction: a review. Environ Health Perspect 124:1683-1693; http://dx.doi.org/10.1289/EHP136.

From the Cover: Comparative Numerical Modeling of Inhaled Nanoparticle Deposition in Human and Rat Nasal Cavities

To gain a better understanding of nanoparticle exposure in human nasal cavities, laboratory animals (e.g. rat) are used for in vivo studies. However, due to anatomical differences between human and rodent nasal cavities, direct particle deposition comparisons between species are difficult. This paper presents a comparative nanoparticle (1 nm, 10 nm, and 100 nm) deposition study using anatomically realistic models of a human and rat nasal cavity. The particle deposition fraction was highest consistently in the main nasal passage, for all nanoparticles tested, in the human model; whereas this was only the case for 10 nm, and 100 nm particles for the rodent model, where greater deposition was found in the anterior nose for 1 nm particles. A deposition intensity (DI) term was introduced to represent the accumulated deposition fraction on cross-sectional slices. A common and preferential deposition site in the human model was found for all nanoparticles occurring at a distance of 3.5 cm inside the nasal passage. For the rodent model maximum DI occurred in the vestibule region at a distance of 0.3 cm, indicating that the rodent vestibule produces exceptionally high particle filtration capability. We also introduced a deposition flux which was a ratio of the regional deposition fraction relative to the region's surface area fraction. This value allowed direct comparison of deposition flux between species, and a regional extrapolation scaling factor was found (e.g. 1/10 scale for vestibule region for rat to human comparison). This study bridges the in vitro exposure experiments and in vivo nanomaterials toxicity studies, and can contribute towards improving inter-species exposure extrapolation studies in the future.

EPHECT III: Health risk assessment of exposure to household consumer products

In the framework of the EU EPHECT project (Emissions, Exposure Patterns and Health Effects of Consumer Products in the EU), irritative and respiratory effects were assessed in relation to acute (30-min) and long-term (24-h) inhalation exposure to key and emerging indoor air pollutants emitted during household use of selected consumer products. A detailed Health Risk Assessment (HRA) was performed for five selected pollutants of respiratory health relevance, namely acrolein, formaldehyde, naphthalene, d-limonene and α-pinene. For each pollutant, the Critical Exposure Limit (CEL) was compared to indoor air concentrations and exposure estimates for the use of 15 selected consumer products by two population groups (housekeepers and retired people) in the four geographical regions of Europe (North, West, South, East), which were derived previously based on microenvironmental modelling. For the present HRA, health-based CELs were derived for certain compounds in case indoor air quality guidelines were not available by the World Health Organization for end-points relevant to the current study. For each pollutant, the highest indoor air concentrations in each microenvironment and exposure estimates across home microenvironments during the day were lower than the corresponding acute and long-term CELs. However, considerable contributions, especially to acute exposures, were obtained in some cases, such as formaldehyde emissions resulting from single product use of a floor cleaning agent (82% CEL), a candle (10% CEL) and an electric air freshener (17% CEL). Regarding multiple product use, the case of 30-min formaldehyde exposure reaching 34% CEL when eight product classes were used across home microenvironments, i.e. all-purpose/kitchen/floor cleaning agents, furniture/floor polish, combustible/electric air fresheners, and perfume, needs to be highlighted. Such estimated values should be evaluated with caution, as these may be attributed to the exposure scenarios specifically constructed for the present study, following a 'most-representative worst-case scenario' approach for exposure and health risk assessment.

Comparative numerical modeling of inhaled micron-sized particle deposition in human and rat nasal cavities

Micron-sized particle deposition in anatomically realistic models of a rat and human nasal cavity was numerically investigated. A steady laminar inhalation flow rate was applied and particles were released from the outside air. Particles showing equivalent total particle deposition fractions were classified into low, medium and high inertial particle. Typical particle sizes are 2.5, 9 and 20 μm for the human model and 1, 2 and 3 μm for the rat model, respectively. Using a surface-mapping technique the 3D nasal cavity surface was "unwrapped" into a 2D domain and the particle deposition locations were plotted for complete visual coverage of the domain surface. The total surface area comparison showed that the surface area of the human nasal model was about ten times the size of the rat model. In contrast, the regional surface area percentage analysis revealed the olfactory region of the rat model was significantly larger than all other regions making up ∼55.6% of the total surface area, while that of the human nasal model only occupying 10.5%. Flow pattern comparisons showed rapid airflow acceleration was found at the nasopharynx region and the nostril region for the human and rat model, respectively. For the human model, the main passage is the major deposition region for micro-particles. While for the rat model, it is the vestibule. Through comparing the regional deposition flux between human and rat models, this study can contribute towards better extrapolation approach of inhalation exposure data between inter-subject species.

Comparative Risks of Aldehyde Constituents in Cigarette Smoke Using Transient Computational Fluid Dynamics/Physiologically Based Pharmacokinetic Models of the Rat and Human Respiratory Tracts

Computational fluid dynamics (CFD) modeling is well suited for addressing species-specific anatomy and physiology in calculating respiratory tissue exposures to inhaled materials. In this study, we overcame prior CFD model limitations to demonstrate the importance of realistic, transient breathing patterns for predicting site-specific tissue dose. Specifically, extended airway CFD models of the rat and human were coupled with airway region-specific physiologically based pharmacokinetic (PBPK) tissue models to describe the kinetics of 3 reactive constituents of cigarette smoke: acrolein, acetaldehyde and formaldehyde. Simulations of aldehyde no-observed-adverse-effect levels for nasal toxicity in the rat were conducted until breath-by-breath tissue concentration profiles reached steady state. Human oral breathing simulations were conducted using representative aldehyde yields from cigarette smoke, measured puff ventilation profiles and numbers of cigarettes smoked per day. As with prior steady-state CFD/PBPK simulations, the anterior respiratory nasal epithelial tissues received the greatest initial uptake rates for each aldehyde in the rat. However, integrated time- and tissue depth-dependent area under the curve (AUC) concentrations were typically greater in the anterior dorsal olfactory epithelium using the more realistic transient breathing profiles. For human simulations, oral and laryngeal tissues received the highest local tissue dose with greater penetration to pulmonary tissues than predicted in the rat. Based upon lifetime average daily dose comparisons of tissue hot-spot AUCs (top 2.5% of surface area-normalized AUCs in each region) and numbers of cigarettes smoked/day, the order of concern for human exposures was acrolein > formaldehyde > acetaldehyde even though acetaldehyde yields were 10-fold greater than formaldehyde and acrolein.

Advances in Inhalation Dosimetry Models and Methods for Occupational Risk Assessment and Exposure Limit Derivation

The purpose of this article is to provide an overview and practical guide to occupational health professionals concerning the derivation and use of dose estimates in risk assessment for development of occupational exposure limits (OELs) for inhaled substances. Dosimetry is the study and practice of measuring or estimating the internal dose of a substance in individuals or a population. Dosimetry thus provides an essential link to understanding the relationship between an external exposure and a biological response. Use of dosimetry principles and tools can improve the accuracy of risk assessment, and reduce the uncertainty, by providing reliable estimates of the internal dose at the target tissue. This is accomplished through specific measurement data or predictive models, when available, or the use of basic dosimetry principles for broad classes of materials. Accurate dose estimation is essential not only for dose-response assessment, but also for interspecies extrapolation and for risk characterization at given exposures. Inhalation dosimetry is the focus of this paper since it is a major route of exposure in the workplace. Practical examples of dose estimation and OEL derivation are provided for inhaled gases and particulates.

Tissue sensitivity of the rat upper and lower extrapulmonary airways to the inhaled electrophilic air pollutants diacetyl and acrolein

The target site for inhaled vapor-induced injury often differs in mouth-breathing humans compared with nose-breathing rats, thus complicating the use of rat inhalation toxicity data for assessment of human risk. We sought to examine sensitivity of respiratory/transitional nasal (RTM) and tracheobronchial (TBM) mucosa to two electrophilic irritant vapors: diacetyl and acrolein. Computational fluid dynamic physiologically based pharmacokinetic modeling was coupled with biomarker assessment to establish delivered dose-response relationships in RTM and TBM in male F344 rats following 6 h exposure to diacetyl or acrolein. Biomarkers included glutathione status, proinflammatory and antioxidant gene mRNA levels, and nuclear translocation of nuclear factor (erythroid-derived 2)-like 2 (Nrf2). Modeling revealed that 0.0094-0.1653 μg acrolein/min-cm(2) and 3.9-21.6 μg diacetyl/min-cm(2) were deposited into RTM/TBM. Results indicate RTM and TBM were generally of similar sensitivity to diacetyl and acrolein. For instance, both tissues displayed induction of antioxidant and proinflammatory genes, and nuclear accumulation of Nrf2 after electrophile exposure. Hierarchical cellular response patterns were similar in RTM and TBM but differed between vapors. Specifically, diacetyl exposure induced proinflammatory and antioxidant genes concomitantly at low exposure levels, whereas acrolein induced antioxidant genes at much lower exposure levels than that required to induce proinflammatory genes. Generally, diacetyl was less potent than acrolein, as measured by maximal induction of transcripts. In conclusion, the upper and lower extrapulmonary airways are of similar sensitivity to inhaled electrophilic vapors. Dosimetrically based extrapolation of nasal responses in nose-breathing rodents may provide an approach to predict risk to the lower airways of humans during mouth-breathing.

Acrolein and asthma attack prevalence in a representative sample of the United States adult population 2000-2009

Background

Acrolein is an air toxic and highly potent respiratory irritant. There is little epidemiology available, but US EPA estimates that outdoor acrolein is responsible for about 75 percent of non-cancer respiratory health effects attributable to air toxics in the United States, based on the Agency's 2005 NATA (National-Scale Air Toxics Assessment) and acrolein's comparatively potent inhalation reference concentration of 0.02 µg/m³.

Objectives

Assess the association between estimated outdoor acrolein exposure and asthma attack reported by a representative cross-sectional sample of the adult United States population.

Methods

NATA 2005 chronic outdoor acrolein exposure estimates at the census tract were linked with residences oif adults (≥18 years old) in the NHIS (National Health Interview Survey) 2000 – 2009 (n = 271,348 subjects). A sample-weighted logistic regression model characterized the association between the prevalence of reporting at least one asthma attack in the 12 months prior to survey interview and quintiles of exposure to outdoor acrolein, controlling for potential confounders.

Results

In the highest quintile of outdoor acrolein exposure (0.05 – 0.46 µg/m³), there was a marginally significant increase in the asthma attack pOR (prevalence-odds ratio [95% CI]  = 1.08 [0.98∶1.19]) relative to the lowest quintile. The highest quintile was also associated with a marginally significant increase in prevalence-odds (1.13 [0.98∶1.29]) in a model limited to never smokers (n = 153,820).

Conclusions

Chronic exposure to outdoor acrolein of 0.05 – 0.46 µg/m³ appears to increase the prevalence-odds of having at least one asthma attack in the previous year by 8 percent in a representative cross-sectional sample of the adult United States population.

Nanoparticle mass transfer from lung airways to systemic regions--Part II: Multi-compartmental modeling

This is the second article of a two-part paper, combining high-resolution computer simulation results of inhaled nanoparticle deposition in a human airway model (Kolanjiyil and Kleinstreuer, 2013, "Nanoparticle Mass Transfer From Lung Airways to Systemic Regions--Part I: Whole-Lung Aerosol Dynamics," ASME J. Biomech. Eng., 135(12), p. 121003) with a new multicompartmental model for insoluble nanoparticle barrier mass transfer into systemic regions. Specifically, it allows for the prediction of temporal nanoparticle accumulation in the blood and lymphatic systems and in organs. The multicompartmental model parameters were determined from experimental retention and clearance data in rat lungs and then the validated model was applied to humans based on pharmacokinetic cross-species extrapolation. This hybrid simulator is a computationally efficient tool to predict the nanoparticle kinetics in the human body. The study provides critical insight into nanomaterial deposition and distribution from the lungs to systemic regions. The quantitative results are useful in diverse fields such as toxicology for exposure-risk analysis of ubiquitous nanomaterial and pharmacology for nanodrug development and targeting.

Effects of endogenous formaldehyde in nasal tissues on inhaled formaldehyde dosimetry predictions in the rat, monkey, and human nasal passages

Formaldehyde is a nasal carcinogen in rodents at high doses and is an endogenous compound that is present in all living cells. Due to its high solubility and reactivity, quantitative risk estimates for inhaled formaldehyde have relied on internal dose estimates in the upper respiratory tract. Dosimetry calculations are complicated by the presence of endogenous formaldehyde concentrations in the respiratory mucosa. Anatomically accurate computational fluid dynamics (CFD) models of the rat, monkey, and human nasal passages were used to simulate uptake of inhaled formaldehyde. An epithelial structure was implemented in the nasal CFD models to estimate formaldehyde absorption from air:tissue partitioning, species-specific metabolism, first-order clearance, DNA binding, and endogenous formaldehyde production. At an exposure concentration of 1 ppm, predicted formaldehyde nasal uptake was 99.4, 86.5, and 85.3% in the rat, monkey, and human, respectively. Endogenous formaldehyde in nasal tissues did not significantly affect wall mass flux or nasal uptake predictions at exposure concentrations > 500 ppb; however, reduced nasal uptake was predicted at lower exposure concentrations. At an exposure concentration of 1 ppb, predicted nasal uptake was 17.5 and 42.8% in the rat and monkey; net desorption of formaldehyde was predicted in the human model. The nonlinear behavior of formaldehyde nasal absorption will affect the dose-response analysis and subsequent risk estimates at low exposure concentrations. Updated surface area partitioning of nonsquamous epithelium and average flux values in regions where DNA-protein cross-links and cell proliferation rates were measured in rats and monkeys are reported for use in formaldehyde risk models of carcinogenesis.

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

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.

Dynamic multiscale boundary conditions for 4D CT of healthy and emphysematous rats

Changes in the shape of the lung during breathing determine the movement of airways and alveoli, and thus impact airflow dynamics. Modeling airflow dynamics in health and disease is a key goal for predictive multiscale models of respiration. Past efforts to model changes in lung shape during breathing have measured shape at multiple breath-holds. However, breath-holds do not capture hysteretic differences between inspiration and expiration resulting from the additional energy required for inspiration. Alternatively, imaging dynamically--without breath-holds--allows measurement of hysteretic differences. In this study, we acquire multiple micro-CT images per breath (4DCT) in live rats, and from these images we develop, for the first time, dynamic volume maps. These maps show changes in local volume across the entire lung throughout the breathing cycle and accurately predict the global pressure-volume (PV) hysteresis. Male Sprague-Dawley rats were given either a full- or partial-lung dose of elastase or saline as a control. After three weeks, 4DCT images of the mechanically ventilated rats under anesthesia were acquired dynamically over the breathing cycle (11 time points, ≤100 ms temporal resolution, 8 cmH2O peak pressure). Non-rigid image registration was applied to determine the deformation gradient--a numerical description of changes to lung shape--at each time point. The registration accuracy was evaluated by landmark identification. Of 67 landmarks, one was determined misregistered by all three observers, and 11 were determined misregistered by two observers. Volume change maps were calculated on a voxel-by-voxel basis at all time points using both the Jacobian of the deformation gradient and the inhaled air fraction. The calculated lung PV hysteresis agrees with pressure-volume curves measured by the ventilator. Volume maps in diseased rats show increased compliance and ventilation heterogeneity. Future predictive multiscale models of rodent respiration may leverage such volume maps as boundary conditions.

Biologically-based modeling insights in inhaled vapor absorption and dosimetry

The lung is a route of entry and also a target site for inhaled vapors, therefore, knowledge of the total absorbed dose and/or the dose absorbed in each airway during inhalation exposure is essential. Vapor absorption characteristics result primarily from the fact that vapors demonstrate equilibrium/saturation behavior in fluids. Thus, during inhalation exposures blood and airway tissue vapor concentrations increase to a steady state value and increase no further no matter how long the exposure. High tissue concentrations can be obtained with highly soluble vapors, thus solubility, as measured by blood:air partition coefficient, is a fundamentally important physical/chemical characteristic of vapors. While it is classically thought that vapor absorption occurs only in the alveoli it is now understood that this is not the case. Soluble vapors can be efficiently absorbed in the airways themselves and do not necessarily penetrate to the alveolar level. Such vapors are more likely to injure the proximal than distal airways because that is the site of the greatest delivered dose. There are substantial species differences in airway vapor absorption between laboratory animals and humans making interpretation of laboratory animal inhalation toxicity data difficult. Airway absorption is dependent on vapor solubility and is enhanced by local metabolism and/or direct reaction within airway tissues. Modern simulation models that incorporate terms for solubility, metabolism, and reaction rate accurately predict vapor absorption patterns in both animals and humans and have become essential tools for understanding the pharmacology and toxicology of airborne vapors.

Phase-contrast MRI and CFD modeling of apparent ³He gas flow in rat pulmonary airways

Phase-contrast (PC) magnetic resonance imaging (MRI) with hyperpolarized ³He is potentially useful for developing and testing patient-specific models of pulmonary airflow. One challenge, however, is that PC-MRI provides apparent values of local ³He velocity that not only depend on actual airflow but also on gas diffusion. This not only blurs laminar flow patterns in narrow airways but also introduces anomalous airflow structure that reflects gas-wall interactions. Here, both effects are predicted in a live rat using computational fluid dynamics (CFD), and for the first time, simulated patterns of apparent ³He gas velocity are compared with in vivo PC-MRI. Results show (1) that correlations (R²) between measured and simulated airflow patterns increase from 0.23 to 0.79 simply by accounting for apparent ³He transport, and (2) that remaining differences are mainly due to uncertain airway segmentation and partial volume effects stemming from relatively coarse MRI resolution. Higher-fidelity testing of pulmonary airflow predictions should therefore be possible with future imaging improvements.

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.

Update on a Pharmacokinetic-Centric Alternative Tier II Program for MMT-Part II: Physiologically Based Pharmacokinetic Modeling and Manganese Risk Assessment

Recently, a variety of physiologically based pharmacokinetic (PBPK) models have been developed for the essential element manganese. This paper reviews the development of PBPK models (e.g., adult, pregnant, lactating, and neonatal rats, nonhuman primates, and adult, pregnant, lactating, and neonatal humans) and relevant risk assessment applications. Each PBPK model incorporates critical features including dose-dependent saturable tissue capacities and asymmetrical diffusional flux of manganese into brain and other tissues. Varied influx and efflux diffusion rate and binding constants for different brain regions account for the differential increases in regional brain manganese concentrations observed experimentally. We also present novel PBPK simulations to predict manganese tissue concentrations in fetal, neonatal, pregnant, or aged individuals, as well as individuals with liver disease or chronic manganese inhalation. The results of these simulations could help guide risk assessors in the application of uncertainty factors as they establish exposure guidelines for the general public or workers.

A validated hybrid computational fluid dynamics-physiologically based pharmacokinetic model for respiratory tract vapor absorption in the human and rat and its application to inhalation dosimetry of diacetyl

Diacetyl vapor is associated with bronchiolar injury in man but primarily large airway injury in the rat. The goal of this study was to develop a physiologically based pharmacokinetic model for inspired vapor dosimetry and to apply the model to diacetyl. The respiratory tract was modeled as a series of airways: nose, trachea, main bronchi, large bronchi, small bronchi, bronchioles, and alveoli with tissue dimensions obtained from the literature. Airborne vapor was allowed to absorb (or desorb) from tissues based on mass transfer coefficients. Transfer of vapor within tissues was based on molecular diffusivity with direct reaction with tissue substrates and/or metabolism being allowed in each tissue compartment. In vitro studies were performed to provide measures of diacetyl metabolism kinetics and direct reaction rates allowing for the development of a model with no unassigned variables. Respiratory tract uptake of halothane, acetone, ethanol and diacetyl was measured in male F344 rat to obtain data for model validation. The human model was validated against published values for inspired vapor uptake. For both the human and rat models, a close concordance of model estimates with experimental measurements was observed, validating the model. The model estimates that limited amounts of inspired diacetyl penetrate to the bronchioles of the rat (<2%), whereas in the lightly exercising human, 24% penetration to the bronchioles is estimated. Bronchiolar tissue concentrations of diacetyl in the human are estimated to exceed those in the rat by 40-fold. These inhalation dosimetric differences may contribute to the human-rat differences in diacetyl-induced airway injury.

Effects of Surface Smoothness on Inertial Particle Deposition in Human Nasal Models

Computational fluid dynamics (CFD) predictions of inertial particle deposition have not compared well with data from nasal replicas due to effects of surface texture and the resolution of tomographic images. To study effects of geometric differences between CFD models and nasal replicas, nasal CFD models with different levels of surface smoothness were reconstructed from the same MRI data used to construct the nasal replica used by Kelly et al. (2004) [Aerosol Sci. Technol. 38:1063-1071]. One CFD model in particular was reconstructed without any surface smoothing to preserve the detailed topology present in the nasal replica. Steady-state inspiratory airflow and Lagrangian particle tracking were simulated using Fluent software. Particle deposition estimates from the smoother models under-predicted nasal deposition from replica casts, which was consistent with previous findings. These discrepancies were overcome by including surface artifacts that were not present in the reduced models and by plotting deposition efficiency versus the Stokes number, where the characteristic diameter was defined in terms of the pressure-flow relationship to account for changes in airflow resistance due to wall roughness. These results indicate that even slight geometric differences have significant effects on nasal deposition and that this information should be taken into account when comparing particle deposition data from CFD models with experimental data from nasal replica casts.

Analysis of manganese tracer kinetics and target tissue dosimetry in monkeys and humans with multi-route physiologically based pharmacokinetic models

Manganese (Mn) is an essential nutrient with the capacity for toxicity from excessive exposure. Accumulation of Mn in the striatum, globus pallidus, and other midbrain regions is associated with neurotoxicity following high-dose Mn inhalation. Physiologically based pharmacokinetic (PBPK) models for ingested and inhaled Mn in rats and nonhuman primates were previously developed. The models contained saturable Mn tissue-binding capacities, preferential fluxes of Mn in specific tissues, and homeostatic control processes such as inducible biliary excretion of Mn. In this study, a nonhuman primate model was scaled to humans and was further extended to include iv, ip, and sc exposure routes so that past studies regarding radiolabeled carrier-free (54)MnCl(2) tracer kinetics could be evaluated. Simulation results accurately recapitulated the biphasic elimination behavior for all exposure routes. The PBPK models also provided consistent cross-species descriptions of Mn tracer kinetics across multiple exposure routes. These results indicate that PBPK models can accurately simulate the overall kinetic behavior of Mn and predict conditions where exposures will increase free Mn in various tissues throughout the body. Simulations with the human model indicate that globus pallidus Mn concentrations are unaffected by air concentrations < 10 μg/m(3) Mn. The use of this human Mn PBPK model can become a key component of future human health risk assessment of Mn, allowing the consideration of various exposure routes, natural tissue background levels, and homeostatic controls to explore exposure conditions that lead to increased target tissue levels resulting from Mn overexposure.

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

Human exposure to hydrogen sulfide (H(2)S) gas occurs from natural and industrial sources and can result in dose-related neurological, respiratory, and cardiovascular effects. Olfactory neuronal loss in H(2)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(2)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(2)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(2)S dosimetry in the olfactory region was predicted to be similar among individuals. Differences in the 99 th 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(2)S dosimetry in the olfactory region.

Multi-dose-route, multi-species pharmacokinetic models for manganese and their use in risk assessment

Manganese (Mn) is an essential element that may be toxic in conditions of overexposure. Nearly 10 years ago, some of the authors of this article published a proposed methodology to perform a tissue-dose-based risk assessment and a detailed list of data needs necessary to perform the assessment. Since that time, a substantial body of Mn pharmacokinetic (PK) data has been generated in rats and nonhuman primates, allowing for the construction of physiologically based pharmacokinetic (PBPK) models for Mn. This study reviews the development of the Mn PBPK models, reassesses the previously identified data needs, and details potential uses of these models in risk assessment of Mn. Based upon numerous animal experiments, pharmacokinetic (PK) models have effectively simulated tissue kinetics of Mn from both inhaled and oral Mn intake. PK models achieve this by incorporating homeostatic control processes, saturable tissue binding capacities, and preferential fluxes in various tissue regions. While minor data gaps still exist, the models captured the main dose-dependent characteristics of Mn disposition in rodents and monkeys and provide a structure to parameterize an equivalent PK description in humans. These models are organized to contribute to a tissue-dose based risk assessment of Mn that simultaneously considers ingestion and inhalation kinetics of Mn along with homeostatic control of Mn.

Evaluating placental transfer and tissue concentrations of manganese in the pregnant rat and fetuses after inhalation exposures with a PBPK model

A Physiologically Based Pharmaco Kinetic (PBPK) model, based on a published description of manganese (Mn) kinetics in adult rats, has been developed to describe Mn uptake and tissue distribution in the pregnant dam and fetus during dietary and inhalation exposures. This extension incorporated key physiological processes controlling Mn pharmacokinetics during pregnancy and fetal development. After calibration against tissue Mn concentrations observed during late gestation, the model accurately simulated Mn tissue distribution in the dam and fetus following both diet and inhalation exposures to the pregnant rat. Maternal to fetal transfer of Mn through placenta was described using two pathways: a saturable active transport with high affinity and a simple diffusion. The active transport dominates at basal and lower Mn exposure, whereas at higher Mn exposure, the relative contribution of the diffusion pathway increases. To simulate fetal tissue Mn, tissue-binding parameters and preferential influx/efflux rates in fetal brain were adjusted from the adult model based on differential developmental processes and varying tissue demands for Mn in early life. Model simulations were consistent with observed tissue Mn concentrations in fetal tissues, including brain for diet alone and for combined diet and inhalation. Simulations of Mn in placenta and other maternal tissues in late gestation correlated well with measured tissue concentrations. This model, together with our published models for Mn kinetics during lactation and postnatal development, will help to address concerns about Mn neurotoxicity in potentially sensitive human subpopulation, such as infants and children by providing an estimate of Mn exposure in the population of interest.

Upper respiratory tract uptake of naphthalene

Naphthalene is a nasal toxicant and carcinogen in the rat. Upper respiratory tract (URT) uptake of naphthalene was measured in the male and female F344 rat at exposure concentrations of 1, 4, 10, or 30 ppm at inspiratory flow rates of 150 or 300 ml/min. To assess the potential importance of nasal cytochrome P450 (CYP) metabolism, groups of rats were pretreated with the CYP inhibitor 5-phenyl-1-pentyne (PP) (100 mg/kg, ip). In vitro metabolism of naphthalene was similar in nasal tissues from both genders and was reduced by 80% by the inhibitor. URT uptake in female rats was concentration dependent with uptake efficiencies (flow 150 ml/min) of 56, 40, 34, and 28% being observed at inspired concentrations of 1, 4, 10, and 30 ppm, respectively. A similar effect was observed in male rats (flow 150 ml/min) with uptake efficiencies of 57, 49, 37, and 36% being observed. Uptake was more efficient in the male than female rat, likely due to their larger size (226 vs. 144 g). Uptake of naphthalene was significantly reduced by inhibitor pretreatment with the effect being greater at the lower inspired concentrations. Specifically, in pretreated female rats (150 ml/min), URT uptake averaged 25, 29, and 26% at inspired concentrations of 4, 10, and 30 ppm, respectively. Thus, the concentration dependence of uptake was virtually abolished by PP pretreatment. These results provide evidence that nasal CYP metabolism of naphthalene contributes to URT scrubbing of this vapor and is also involved in the concentration dependence of uptake that is observed.