Automated bidirectional coupling of multiscale models of aerosol dosimetry: validation with subject-specific deposition data.

Assessing the toxicity of airborne particulate matter or the efficacy of inhaled drug depends upon accurate estimates of deposited fraction of inhaled materials. In silico approaches can provide important insights into site- or airway-specific deposition of inhaled aerosols in the respiratory system. In this study, we improved on our recently developed 3D/1D model that simulate aerosol transport and deposition in the whole lung over multiple breath cycles (J. Aerosol Sci 151:105647). A subject-specific multiscale lung model of a healthy male subject using computational fluid-particle dynamics (CFPD) in a 3D model of the oral cavity through the large bronchial airways entering each lobe was bidirectionally coupled with a recently improved Multiple Path Particle Dosimetry (MPPD) model to predict aerosol deposition over the entire respiratory tract over multiple breaths for four conditions matching experimental aerosol exposures in the same subject from which the model was developed. These include two particle sizes (1 and 2.9 µm) and two subject-specific breathing rates of ~300 ml/s (slow breathing) and ~750 ml/s (fast breathing) at a target tidal volume of 1 L. In silico predictions of retained fraction were 0.31 and 0.29 for 1 µm and 0.66 and 0.62 for 2.9 µm during slow and fast breathing, respectively, and compared well with experimental data (1 µm: 0.31±0.01 (slow) and 0.27±0.01 (fast), 2.9 µm: 0.63±0.03 (slow) and 0.68±0.02 (fast)). These results provide a great deal of confidence in the validity and reliability of our approach.

Influence of alveolar mixing and multiple breaths of aerosol intake on particle deposition in the human lungs.

Predictive dosimetry models play an important role in assessing health effect of inhaled particulate matter and in optimizing delivery of inhaled pharmaceutical aerosols. In this study, the commonly used 1D Multiple-Path Particle Dosimetry model (MPPD) was improved by including a mechanistically based model component for alveolar mixing of particles and by extending the model capabilities to account for multiple breaths of aerosol intake. These modifications increased the retained fraction of particles and consequently particle deposition predictions in the deep lung during tidal breathing. Comparison with an existing dataset (J. Aerosol Sci., 99:27-39, 2016) obtained under two breathing conditions referred to as slow and fast breathing showed significant differences in 1 µm particle deposition between predictions based on subject-specific breathing patterns and lung volume (slow: 30 ± 1%, fast: 21 ± 1%, (average ± standard deviation), N = 7) and measurements (slow: 43 ± 9%, fast: 30 ± 5%) when the prior version of MPPD (single breath and no mixing, J. Aerosol Sci., 151:105647, 2021) was used. Adding a mixing model and multiple breaths moved the predictions (slow: 34 ± 2%, fast:25 ± 2%) closer to the range of deposition measurements. For 2.9 µm particles, predictions from both the original (slow: 70 ± 2%, fast: 57 ± 2%) and the revised MPPD model (slow: 71 ± 2%, fast: 59 ± 3%) compared well with experiments (slow: 67 ± 8%, fast: 58 ± 10%). This was expected as suspended fraction of 2.9 µm particles was small and thus the addition of alveolar mixing and multi breath capability only slightly increased the retained fraction for particles of this size and greater. The revised 1D model improves dose predictions in the deep lung and support human risk assessment from exposure to airborne particles.

Efficient bi-directional coupling of 3D Computational Fluid-Particle Dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry.

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. As a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFPD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. These methods can be used as the foundation for multi-scale modeling of the full respiratory system.

Computational fluid dynamics modeling of Bacillus anthracis spore deposition in rabbit and human respiratory airways.

Three-dimensional computational fluid dynamics and Lagrangian particle deposition models were developed to compare the deposition of aerosolized Bacillus anthracis spores in the respiratory airways of a human with that of the rabbit, a species commonly used in the study of anthrax disease. The respiratory airway geometries for each species were derived respectively from computed tomography (CT) and µCT images. Both models encompassed airways that extended from the external nose to the lung with a total of 272 outlets in the human model and 2878 outlets in the rabbit model. All simulations of spore deposition were conducted under transient, inhalation-exhalation breathing conditions using average species-specific minute volumes. Two different exposure scenarios were modeled in the rabbit based upon experimental inhalation studies. For comparison, human simulations were conducted at the highest exposure concentration used during the rabbit experimental exposures. Results demonstrated that regional spore deposition patterns were sensitive to airway geometry and ventilation profiles. Due to the complex airway geometries in the rabbit nose, higher spore deposition efficiency was predicted in the nasal sinus compared to the human at the same air concentration of anthrax spores. In contrast, higher spore deposition was predicted in the lower conducting airways of the human compared to the rabbit lung due to differences in airway branching pattern. This information can be used to refine published and ongoing biokinetic models of inhalation anthrax spore exposures, which currently estimate deposited spore concentrations based solely upon exposure concentrations and inhaled doses that do not factor in species-specific anatomy and physiology for deposition.

Extension of a PBPK model for ethylene glycol and glycolic acid to include the competitive formation and clearance of metabolites associated with kidney toxicity in rats and humans.

A previously developed PBPK model for ethylene glycol and glycolic acid was extended to include glyoxylic acid, oxalic acid, and the precipitation of calcium oxalate that is associated with kidney toxicity in rats and humans. The development and evaluation of the PBPK model was based upon previously published pharmacokinetic studies coupled with measured blood and tissue partition coefficients and rates of in vitro metabolism of glyoxylic acid to oxalic acid, glycine and other metabolites using primary hepatocytes isolated from male Wistar rats and humans. Precipitation of oxalic acid with calcium in the kidneys was assumed to occur only at concentrations exceeding the thermodynamic solubility product for calcium oxalate. This solubility product can be affected by local concentrations of calcium and other ions that are expressed in the model using an ion activity product estimated from toxicity studies such that calcium oxalate precipitation would be minimal at dietary exposures below the NOAEL for kidney toxicity in the sensitive male Wistar rat. The resulting integrated PBPK predicts that bolus oral or dietary exposures to ethylene glycol would result in typically 1.4-1.6-fold higher peak oxalate levels and 1.6-2-fold higher AUC's for calcium oxalate in kidneys of humans as compared with comparably exposed male Wistar rats over a dose range of 1-1000 mg/kg. The converse (male Wistar rats predicted to have greater oxalate levels in the kidneys than humans) was found for inhalation exposures although no accumulation of calcium oxalate is predicted to occur until exposures are well in excess of the theoretical saturated vapor concentration of 200mg/m(3). While the current model is capable of such cross-species, dose, and route-of-exposure comparisons, it also highlights several areas of potential research that will improve confidence in such predictions, especially at low doses relevant for most human exposures.

Magnetic resonance imaging and computational fluid dynamics (CFD) simulations of rabbit nasal airflows for the development of hybrid CFD/PBPK models.

The percentages of total airflows over the nasal respiratory and olfactory epithelium of female rabbits were calculated from computational fluid dynamics (CFD) simulations of steady-state inhalation. These airflow calculations, along with nasal airway geometry determinations, are critical parameters for hybrid CFD/physiologically based pharmacokinetic models that describe the nasal dosimetry of water-soluble or reactive gases and vapors in rabbits. CFD simulations were based upon three-dimensional computational meshes derived from magnetic resonance images of three adult female New Zealand White (NZW) rabbits. In the anterior portion of the nose, the maxillary turbinates of rabbits are considerably more complex than comparable regions in rats, mice, monkeys, or humans. This leads to a greater surface area to volume ratio in this region and thus the potential for increased extraction of water soluble or reactive gases and vapors in the anterior portion of the nose compared to many other species. Although there was considerable interanimal variability in the fine structures of the nasal turbinates and airflows in the anterior portions of the nose, there was remarkable consistency between rabbits in the percentage of total inspired airflows that reached the ethmoid turbinate region (approximately 50%) that is presumably lined with olfactory epithelium. These latter results (airflows reaching the ethmoid turbinate region) were higher than previous published estimates for the male F344 rat (19%) and human (7%). These differences in regional airflows can have significant implications in interspecies extrapolations of nasal dosimetry.

Incorporation of therapeutic interventions in physiologically based pharmacokinetic modeling of human clinical case reports of accidental or intentional overdosing with ethylene glycol.

Although occupational uses of the high production volume (HPV) chemical ethylene glycol (EG) have not been associated with adverse effects, there are case reports where humans have either intentionally or accidentally ingested large quantities of EG, primarily from antifreeze. The acute toxicity of EG can proceed through three stages, each associated with a different metabolite: central nervous system depression (ethylene glycol), cardiopulmonary effects associated with metabolic acidosis (glycolic acid), and ultimately renal toxicity (oxalic acid), depending on the total amounts consumed and the effectiveness of therapeutic interventions. A physiologically based pharmacokinetic (PBPK) model developed in a companion paper (Corley et al., 2005). Development of a physiologically based pharmacokinetic model for ethylene glycol and its metabolite, glycolic acid, in rats and humans. Toxicol. Sci., in press 2005) was refined in this study to include clinically relevant treatment regimens for EG poisoning such as hemodialysis or metabolic inhibition with either ethanol or fomepizole. Such modifications enabled the model to describe data from several human case reports, confirming the ability of the previous model to describe the pharmacokinetics of EG and its metabolite, glycolic acid, in humans across a broad range of doses and multiple exposure routes. By integrating the case report data sets with controlled studies in this PBPK model, it was demonstrated that fomepizole, if administered early enough in a clinical situation, can be more effective than ethanol or hemodialysis in preventing the metabolism of EG to more toxic metabolites. Hemodialysis remains an important option, however, if treatment is instituted after a significant amount of EG is metabolized or if renal toxicity has occurred.

Validation of human physiologically based pharmacokinetic model for vinyl acetate against human nasal dosimetry data.

Vinyl acetate has been shown to induce nasal lesions in rodents in inhalation bioassays. A physiologically based pharmacokinetic (PBPK) model for vinyl acetate has been used in human risk assessment, but previous in vivo validation was conducted only in rats. Controlled human exposures to vinyl acetate were conducted to provide validation data for the application of the model in humans. Five volunteers were exposed to 1, 5, and 10 ppm 13C1,13C2 vinyl acetate via inhalation. A probe inserted into the nasopharyngeal region sampled both 13C1,13C2 vinyl acetate and the major metabolite 13C1,13C2 acetaldehyde during rest and light exercise. Nasopharyngeal air concentrations were analyzed in real time by ion trap mass spectrometry (MS/MS). Experimental concentrations of both vinyl acetate and acetaldehyde were then compared to predicted concentrations calculated from the previously published human model. Model predictions of vinyl acetate nasal extraction compared favorably with measured values of vinyl acetate, as did predictions of nasopharyngeal acetaldehyde when compared to measured acetaldehyde. The results showed that the current PBPK model structure and parameterization are appropriate for vinyl acetate. These analyses were conducted from 1 to 10 ppm vinyl acetate, a range relevant to workplace exposure standards but which would not be expected to saturate vinyl acetate metabolism. Risk assessment based on this model further concluded that 24 h per day exposures up to 1 ppm do not present concern regarding cancer or non-cancer toxicity. Validation of the vinyl acetate human PBPK model provides support for these conclusions.

Derivation of a human equivalent concentration for n-butanol using a physiologically based pharmacokinetic model for n-butyl acetate and metabolites n-butanol and n-butyric acid.

The metabolic series approach for risk assessment uses a dosimetry-based analysis to develop toxicity information for a group of metabolically linked compounds using pharmacokinetic (PK) data for each compound and toxicity data for the parent compound. The metabolic series approach for n-butyl acetate and its subsequent metabolites, n-butanol and n-butyric acid (the butyl series), was first demonstrated using a provisional physiologically based pharmacokinetic (PBPK) model for the butyl series. The objective of this work was to complete development of the PBPK model for the butyl series. Rats were administered test compounds by iv bolus dose, iv infusion, or by inhalation in a recirculating closed chamber. Hepatic, vascular, and extravascular metabolic constants for metabolism were estimated by fitting the model to the blood time course data from these experiments. The respiratory bioavailability of n-butyl acetate (100% of alveolar ventilation) and n-butanol (50% of alveolar ventilation) was estimated from closed chamber inhalation studies and measured ventilation rates. The resulting butyl series PBPK model successfully reproduces the blood time course of these compounds following iv administration and inhalation exposure to n-butyl acetate and n-butanol in rats and arterial blood n-butanol kinetics following inhalation exposure to n-butanol in humans. These validated inhalation route models can be used to support species and dose-route extrapolations required for risk assessment of butyl series family of compounds. Human equivalent concentrations of 169 ppm and 1066 ppm n-butanol corresponding to the rat n-butyl acetate NOAELs of 500 and 3000 ppm were derived using the models.

Development of a physiologically based pharmacokinetic model for ethylene glycol and its metabolite, glycolic Acid, in rats and humans.

An extensive database on the toxicity and modes of action of ethylene glycol (EG) has been developed over the past several decades. Although renal toxicity has long been recognized as a potential outcome, in recent years developmental toxicity, an effect observed only in rats and mice, has become the subject of extensive research and regulatory reviews to establish guidelines for human exposures. The developmental toxicity of EG has been attributed to the intermediate metabolite, glycolic acid (GA), which can become a major metabolite when EG is administered to rats and mice at high doses and dose rates. Therefore, a physiologically based pharmacokinetic (PBPK) model was developed to integrate the extensive mode of action and pharmacokinetic data on EG and GA for use in developmental risk assessments. The resulting PBPK model includes inhalation, oral, dermal, intravenous, and subcutaneous routes of administration. Metabolism of EG and GA were described in the liver with elimination via the kidneys. Metabolic rate constants and partition coefficients for EG and GA were estimated from in vitro studies. Other biochemical constants were optimized from appropriate in vivo pharmacokinetic studies. Several controlled rat and human metabolism studies were used to validate the resulting PBPK model. When internal dose surrogates were compared in rats and humans over a broad range of exposures, it was concluded that humans are unlikely to achieve blood levels of GA that have been associated with developmental toxicity in rats following occupational or environmental exposures.

Development of a physiologically based pharmacokinetic model for propylene glycol monomethyl ether and its acetate in rats and humans.

Propylene glycol monomethyl ether (PM), along with its acetate, is the most widely used of the propylene glycol ether family of solvents. The most common toxic effects of PM observed in animal studies include sedation, very slight alpha(2u)-globulin mediated nephropathy (male rats only) and hepatomegally at high exposures (typically > 1000 ppm). Sedation in animal studies usually resolves within a few exposures to 3000 ppm (the highest concentration used in subchronic and chronic inhalation studies) due to the induction of metabolizing enzymes. Data from a variety of pharmacokinetic and mechanistic studies have been incorporated into a PBPK model for PM and its acetate in rats and mice. Published controlled exposure and workplace biomonitoring studies have also been included for comparisons of the internal dosimetry of PM and its acetate between laboratory animals and humans. PM acetate is rapidly hydrolyzed to PM, which is further metabolized to either glucuronide or sulfate conjugates (minor pathways) or propylene glycol (major pathway). In vitro half-lives for PM acetate range from 14 to 36 min depending upon the tissue and species. In vivo half-lives are considerably faster, reflecting the total contributions of esterases in the blood and tissues of the body, and are on the order of just a few minutes. Thus, very little PM acetate is found in vivo and, other than potential portal of entry irritation, the toxicity of PM acetate is related to PM. Regardless of the source for PM (either PM or its acetate), rats were predicted to have a higher Cmax and AUC for PM in blood than humans, especially at concentrations greater than the current ACGIH TLV of 100 ppm. This would indicate that the major systemic effects of PM would be expected to be less severe in humans than rats at comparable inhalation exposures.

Determination of age and gender differences in biochemical processes affecting the disposition of 2-butoxyethanol and its metabolites in mice and rats to improve PBPK modeling.

2-Butoxyethanol (BE) is the most widely used glycol ether solvent. BEs major metabolite, butoxyacetic acid (BAA), causes hemolysis with significant species differences in sensitivity. Several PBPK models have been developed over the past two decades to describe the disposition of BE and BAA in male rats and humans to refine health risk assessments. More recent efforts by Lee et al. [Lee, K.M., Dill, J.A., Chou, B.J., Roycroft, J.H., 1998. Physiologically based pharmacokinetic model for chronic inhalation of 2-butoxyethanol. Toxicol. Appl. Pharmacol. 153, 211-226] to describe the kinetics of BE and BAA in the National Toxicology Program (NTP) chronic inhalation studies required the use of several assumptions to extrapolate model parameters from earlier PBPK models developed for young male rats to include female F344 and both sexes of B6C3F1 mice and the effects of aging. To replace these assumptions, studies were conducted to determine the impact of age, gender and species on the metabolism of BE, and the tissue partitioning, renal acid transport and plasma protein binding of BAA. In the current study, the Lee et al. PBPK model was updated and expanded to include the further metabolism of BAA and the salivary excretion of BE and BAA which may contribute to the forestomach irritation observed in mice in the NTP study. The revised model predicted that peak blood concentrations of BAA achieved following 6 h inhalation exposures are greatest in young adult female rats at concentrations up to 300 ppm. This is not the case predicted for old (> or =18 months) animals, where peak blood concentrations of BAA in male and female mice were similar to or greater than female rats. The revised model serves as a quantitative tool for integrating an extensive pharmacokinetic and mechanistic database into a format that can readily be used to compare internal dosimetry across dose, route of exposure and species.

Physiologically based pharmacokinetic modeling of the temperature-dependent dermal absorption of chloroform by humans following bath water exposures.

The kinetics of chloroform in the exhaled breath of human volunteers exposed skin-only via bath water (concentrations < 100 ppb) were analyzed using a physiologically based pharmacokinetic (PBPK) model. Significant increases in exhaled chloroform (and thus bioavailability) were observed as exposure temperatures were increased from 30 to 40 degrees C. The blood flows to the skin and effective skin permeability coefficients (Kp) were both varied to reflect the temperature-dependent changes in physiology and exhalation kinetics. At 40 degrees C, no differences were observed between males and females. Therefore, Kps were determined (approximately 0.06 cm/hr) at a skin blood flow rate of 18% of the cardiac output. At 30 and 35 degrees C, males exhaled more chloroform than females, resulting in lower effective Kps calculated for females. At these lower temperatures, the blood flow to the skin was also reduced. Total amounts of chloroform absorbed averaged 41.9 and 43.6 microg for males and 11.5 and 39.9 microg for females exposed at 35 and 40 degrees C, respectively. At 30 degrees C, only 2/5 males and 1/5 females had detectable concentrations of chloroform in their exhaled breath. For perspective, the total intake of chloroform would have ranged from 79-194 microg if the volunteers had consumed 2 liters of water orally at the concentrations used in this study. Thus, the relative contribution of dermal uptake of chloroform to the total body burdens associated with bathing for 30 min and drinking 2 liters of water (ignoring contributions from inhalation exposures) was predicted to range from 1 to 28%, depending on the temperature of the bath.

Assessment of the percutaneous absorption of trichloroethylene in rats and humans using MS/MS real-time breath analysis and physiologically based pharmacokinetic modeling.

The development and validation of noninvasive techniques for estimating the dermal bioavailability of solvents in contaminated soil and water can facilitate the overall understanding of human health risk. To assess the dermal bioavailability of trichloroethylene (TCE), exhaled breath was monitored in real time using an ion trap mass spectrometer (MS/MS) to track the uptake and elimination of TCE from dermal exposures in rats and humans. A physiologically based pharmacokinetic (PBPK) model was used to estimate total bioavailability. Male F344 rats were exposed to TCE in water or soil under occluded or nonoccluded conditions by applying a patch to a clipper-shaved area of the back. Rats were placed in off-gassing chambers and chamber air TCE concentration was quantified for 3-5 h postdosing using the MS/MS. Human volunteers were exposed either by whole-hand immersion or by attaching patches containing TCE in soil or water on each forearm. Volunteers were provided breathing air via a face mask to eliminate inhalation exposure, and exhaled breath was analyzed using the MS/MS. The total TCE absorbed and the dermal permeability coefficient (K(P)) were estimated for each individual by optimization of the PBPK model to the exhaled breath data and the changing media and/or dermal patch concentrations. Rat skin was significantly more permeable than human skin. Estimates for K(P) in a water matrix were 0.31 +/- 0.01 cm/h and 0.015 +/- 0.003 cm/h in rats and humans, respectively. K(P) estimates were more than three times higher from water than soil matrices in both species. K(P) values calculated using the standard Fick's Law equation were strongly affected by exposure length and volatilization of TCE. In comparison, K(P) values estimated using noninvasive real-time breath analysis coupled with the PBPK model were consistent, regardless of volatilization, exposure concentration, or duration.

Utility of real time breath analysis and physiologically based pharmacokinetic modeling to determine the percutaneous absorption of methyl chloroform in rats and humans.

Due to the large surface area of the skin, percutaneous absorption has the potential to contribute significantly to the total bioavailability of some compounds. Breath elimination data, acquired in real-time using a novel MS/MS system, was assessed using a PBPK model with a dermal compartment to determine the percutaneous absorption of methyl chloroform (MC) in rats and humans from exposures to MC in non-occluded soil or occluded water matrices. Rats were exposed to MC using a dermal exposure cell attached to a clipper-shaved area on their back. The soil exposure cell was covered with a charcoal patch to capture volatilized MC and prevent contamination of exhaled breath. This technique allowed the determination of MC dermal absorption kinetics under realistic, non-occluded conditions. Human exposures were conducted by immersing one hand in 0.1% MC in water, or 0.75% MC in soil. The dermal PBPK model was used to estimate skin permeability (Kp) based on the fit of the exhaled breath data. Rat skin K(p)s were estimated to be 0.25 and 0.15 cm/h for MC in water and soil matrices, respectively. In comparison, human permeability coefficients for water matrix exposures were 40-fold lower at 0.006 cm/h. Due to evaporation and differences in apparent Kp, nearly twice as much MC was absorbed from the occluded water (61.3%) compared to the non-occluded soil (32.5%) system in the rat. The PBPK model was used to simulate dermal exposures to MC-contaminated water and soil in children and adults using worst-case EPA default assumptions. The simulations indicate that neither children nor adults will absorb significant amounts of MC from non-occluded exposures, independent of the length of exposure. The results from these simulations reiterate the importance of conducting dermal exposures under realistic conditions.

Proposed occupational exposure limits for select ethylene glycol ethers using PBPK models and Monte Carlo simulations.

Methoxyethanol (ethylene glycol monomethyl ether, EGME), ethoxyethanol (ethylene glycol monoethyl ether, EGEE), and ethoxyethyl acetate (ethylene glycol monoethyl ether acetate, EGEEA) are all developmental toxicants in laboratory animals. Due to the imprecise nature of the exposure data in epidemiology studies of these chemicals, we relied on human and animal pharmacokinetic data, as well as animal toxicity data, to derive 3 occupational exposure limits (OELs). Physiologically based pharmacokinetic (PBPK) models for EGME, EGEE, and EGEEA in pregnant rats and humans have been developed (M. L. Gargas et al., 2000, Toxicol. Appl. Pharmacol. 165, 53-62; M. L. Gargas et al., 2000, Toxicol. Appl. Pharmacol. 165, 63-73). These models were used to calculate estimated human-equivalent no adverse effect levels (NAELs), based upon internal concentrations in rats exposed to no observed effect levels (NOELs) for developmental toxicity. Estimated NAEL values of 25 ppm for EGEEA and EGEE and 12 ppm for EGME were derived using average values for physiological, thermodynamic, and metabolic parameters in the PBPK model. The uncertainties in the point estimates for the NOELs and NAELs were estimated from the distribution of internal dose estimates obtained by varying key parameter values over expected ranges and probability distributions. Key parameters were identified through sensitivity analysis. Distributions of the values of these parameters were sampled using Monte Carlo techniques and appropriate dose metrics calculated for 1600 parameter sets. The 95th percentile values were used to calculate interindividual pharmacokinetic uncertainty factors (UFs) to account for variability among humans (UF(h,pk)). These values of 1.8 for EGEEA/EGEE and 1.7 for EGME are less than the default value of 3 for this area of uncertainty. The estimated human equivalent NAELs were divided by UF(h,pk) and the default UFs for pharmacodynamic variability among animals and among humans to calculate the proposed OELs. This methodology indicates that OELs (8-h time-weighted average) that should protect workers from the most sensitive adverse effects of these chemicals are 2 ppm EGEEA and EGEE (11 mg/m(3) EGEEA, 7 mg/m(3) EGEE) and 0.9 ppm (3 mg/m(3)) EGME. These recommendations assume that dermal exposure will be minimal or nonexistent.

Pharmacokinetics of the trichothecene mycotoxin, T-2 toxin, in swine and cattle.

The pharmacokinetics of the trichothecene mycotoxin, T-2 toxin, were determined in growing gilts and heifers. Following intra-aortal administration in swine and intravenous administration in calves, the disappearance of the parent T-2 toxin followed a 2-compartment open model. Mean elimination phase half-lives were 13.8 and 17.4 min and mean apparent specific volumes of distribution were 0.366 and 0.376 l/kg in swine and calves, respectively. The fraction of T-2 toxin eliminated as parent compound in the urine was negligible. In spite of administration of a lethal oral dose in swine (2.4 mg/kg) and toxic oral doses (up to 3.6 mg/kg) in calves, no parent T-2 toxin was detected in plasma or urine. After intra-aortal administration in swine, tissue concentrations of T-2 toxin were consistently highest in lymphoid organs. Tissue residues of T-2 toxin were rapidly depleted such that, in spite of administration of a potentially lethal intra-aortal dose, no quantifiable T-2 toxin was present in any of the tissues collected at 4 hr after dosing. No T-2 toxin could be detected in liver, even at 1 hr after dosing.

A real-time in-vivo method for studying the percutaneous absorption of volatile chemicals.

Realistic estimates of percutaneous absorption following exposures to solvents in the workplace, or through contaminated soil and water, are critical to understanding human health risks. A method was developed to determine dermal uptake of solvents under non-steady-state conditions using real-time breath analysis in rats, monkeys, and humans. The exhaled breath was analyzed using an ion-trap mass spectrometer, which can quantitate chemicals in the exhaled breath stream in the 1-5 ppb range. The resulting data were evaluated using physiologically-based pharmacokinetic (PBPK) models to estimate dermal permeability constants (Kp) under various exposure conditions. The effects of exposure matrix (soil versus water), occlusion versus non-occlusion, and species differences on the absorption of methyl chloroform, trichloroethylene, and benzene were compared. Exposure concentrations were analyzed before and at 0.5-hour intervals throughout the exposures. The percentage of each chemical absorbed and the corresponding Kp were estimated by optimization of the PBPK model to the medium concentration and the exhaled-breath data. The method was found to be sufficiently sensitive for animal and human dermal studies at low exposure concentrations over small body surface areas, for short periods, using non-steady-state exposure conditions.

In vitro metabolism of benzo[a]pyrene and dibenzo[def,p]chrysene in rodent and human hepatic microsomes.

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous and often carcinogenic contaminants released into the environment during natural and anthropogenic combustion processes. Benzo[a]pyrene (B[a]P) is the prototypical carcinogenic PAH, and dibenzo[def,p]chrysene (DBC) is a less prevalent, but highly potent transplacental carcinogenic PAH. Both are metabolically activated by isoforms of the cytochrome P450 enzyme superfamily to form reactive carcinogenic and cytotoxic metabolites. Metabolism of B[a]P and DBC was studied in hepatic microsomes of male Sprague-Dawley rats, naïve and pregnant female B6129SF1/J mice, and female humans, corresponding to available pharmacokinetic data. Michaelis-Menten saturation kinetic parameters including maximum rates of metabolism (VMAX, nmol/min/mg microsomal protein), affinity constants (KM, µM), and rates of intrinsic clearance (CLINT, ml/min/kg body weight) were calculated from substrate depletion data. CLINT was also estimated from substrate depletion data using the alternative in vitro half-life method. VMAX and CLINT were higher for B[a]P than DBC, regardless of species. Clearance for both B[a]P and DBC was highest in naïve female mice and lowest in female humans. Clearance rates of B[a]P and DBC in male rat were more similar to female human than to female mice. Clearance of DBC in liver microsomes from pregnant mice was reduced compared to naïve mice, consistent with reduced active P450 protein levels and elevated tissue concentrations and residence times for DBC observed in previous in vivo pharmacokinetic studies. These findings suggest that rats are a more appropriate model organism for human PAH metabolism, and that pregnancy's effects on metabolism should be further explored.

A Bidirectional Coupling Procedure Applied to Multiscale Respiratory Modeling.

In this study, we present a novel multiscale computational framework for efficiently linking multiple lower-dimensional models describing the distal lung mechanics to imaging-based 3D computational fluid dynamics (CFD) models of the upper pulmonary airways in order to incorporate physiologically appropriate outlet boundary conditions. The framework is an extension of the Modified Newton's Method with nonlinear Krylov accelerator developed by Carlson and Miller [1, 2, 3]. Our extensions include the retention of subspace information over multiple timesteps, and a special correction at the end of a timestep that allows for corrections to be accepted with verified low residual with as little as a single residual evaluation per timestep on average. In the case of a single residual evaluation per timestep, the method has zero additional computational cost compared to uncoupled or unidirectionally coupled simulations. We expect these enhancements to be generally applicable to other multiscale coupling applications where timestepping occurs. In addition we have developed a "pressure-drop" residual which allows for stable coupling of flows between a 3D incompressible CFD application and another (lower-dimensional) fluid system. We expect this residual to also be useful for coupling non-respiratory incompressible fluid applications, such as multiscale simulations involving blood flow. The lower-dimensional models that are considered in this study are sets of simple ordinary differential equations (ODEs) representing the compliant mechanics of symmetric human pulmonary airway trees. To validate the method, we compare the predictions of hybrid CFD-ODE models against an ODE-only model of pulmonary airflow in an idealized geometry. Subsequently, we couple multiple sets of ODEs describing the distal lung to an imaging-based human lung geometry. Boundary conditions in these models consist of atmospheric pressure at the mouth and intrapleural pressure applied to the multiple sets of ODEs. In both the simplified geometry and in the imaging-based geometry, the performance of the method was comparable to that of monolithic schemes, in most cases requiring only a single CFD evaluation per time step. Thus, this new accelerator allows us to begin combining pulmonary CFD models with lower-dimensional models of pulmonary mechanics with little computational overhead. Moreover, because the CFD and lower-dimensional models are totally separate, this framework affords great flexibility in terms of the type and breadth of the adopted lower-dimensional model, allowing the biomedical researcher to appropriately focus on model design. Research funded by the National Heart and Blood Institute Award 1RO1HL073598.