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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Potential technology for studying dosimetry and response to airborne chemical and biological pollutants

Advances in computational, and imaging techniques have enabled the rapid development of three-dimensional (3-D) models of biological systems in unprecedented detail. Using these advances, 3-D models of the lungs and nasal passages of the rat and human are being developed to ultimately improve predictions of airborne pollutant dosimetry. Techniques for imaging the respiratory tract by magnetic resonance imaging (MRI) were developed to improve the speed and accuracy of geometric data collection for mesh reconstruction. The MRI resolution is comparable to that obtained by manual measurements but at much greater speed and accuracy. Newly developed software (NWGrid) was utilized to translate imaging data from MR into 3-D mesh structures. Together, these approaches significantly reduced the time to develop a 3-D model. This more robust airway structure will ultimately facilitate modeling gas or vapor exchange between the respiratory tract and vasculature as well as enable linkages of dosimetry with cell response models. The 3-D, finite volume, viscoelastic mesh structures form the geometric basis for computational fluid dynamics modeling of inhalation, exhalation and the delivery of individual particles (or concentrations of gas or vapors) to discrete regions of the respiratory tract. The ability of these 3-D models to resolve dosimetry at such a high level of detail will require new techniques to measure regional airflows and particulate deposition for model validation.

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.

Development of a physiologically based pharmacokinetic model for hydroquinone

Hydroquinone (HQ) produces nephrotoxicity and renal tubular adenomas in male F344 rats following 2 years of oral dosing. Female F344 and SD rats are comparatively resistant to these effects. Nephrotoxicity and tumorigenicity have been associated with a minor glutathione conjugation pathway following the oxidation of HQ to benzoquinone (BQ). The majority of administered doses (90-99%) consists of glucuronide and sulfate conjugates of HQ. An initial physiologically based pharmacokinetic model was developed to characterize the role of kinetics in the strain differences observed in HQ-induced renal toxicity and tumorigenicity. Partition coefficients, protein-binding, and metabolic rate constants were determined directly or estimated from a series of in vivo and in vitro studies. Metabolism was confined to the liver and GI tract. The total flux through the glutathione pathway represented the "internal dose" of HQ for nephrotoxicity. Simulations were compared to a variety of data from male and female F344 rats, male SD rats, and a single male human volunteer. Simulations of intraperitoneal administration resulted in higher amounts of glutathione conjugates than comparable oral doses. This was consistent with protein-binding and toxicity studies and emphasized the importance of first-pass GI tract metabolism. In addition, male F344 rats were predicted to form more total glutathione conjugates than SD rats at equivalent dose levels, which was also consistent with the observed strain differences in renal toxicity. This model represents the first stage in the development of a biologically based dose-response model for improving the scientific basis for human health risk assessments of HQ.

A toxicokinetic study of inhaled ethylene glycol monomethyl ether (2-ME) and validation of a physiologically based pharmacokinetic model for the pregnant rat and human

Exposures to sufficiently high doses of ethylene glycol monomethyl ether (2-methoxyethanol, 2-ME) have been found to produce developmental effects in rodents and nonhuman primates. The acetic acid metabolite of 2-ME, 2-methoxyacetic acid (2-MAA), is the likely toxicant, and, as such, an understanding of the kinetics of 2-MAA is important when assessing the potential risks to humans associated with 2-ME. A previously described physiologically based pharmacokinetic (PBPK) model of 2-ME/2-MAA kinetics for rats exposed via oral or iv administration was extended and validated to inhalation exposures. Pregnant Sprague-Dawley rats were exposed for 5 days (gestation days 11-15), 6 h/day, to 2-ME vapor at 10 and 50 ppm. Validation consisted of comparing model output to maternal blood and fetal 2-ME and 2-MAA concentrations during and following 5 days of exposure (gestation days 11-15). These concentrations correspond to a known no observed effect level (NOEL) and a lowest observed effect level (LOEL) for developmental effects in rats. The rat PBPK model for 2-ME/2-MAA was scaled to humans and the model (without the pregnancy component) was used to predict data collected by other investigators on the kinetics of 2-MAA excretion in urine following exposures to 2-ME in human volunteers. The partially validated human model (with the pregnancy component) was used to predict equivalent human exposure concentrations based on 2-MAA dose measures (maximum blood concentration, C(max), and average daily area under the 2-MAA blood concentration curve, AUC, during pregnancy) that correspond to the concentrations measured at the rat NOEL and LOEL exposure concentrations. Using traditional PBPK scale-up techniques, it was calculated that pregnant women exposed for 8 h/day, 5 days/week, for the duration of pregnancy would need to be exposed to 12 or 60 ppm 2-ME to produce maternal 2-MAA blood concentrations (C(max) or average daily AUC) equivalent to those in rats exposed to the NOEL (10 ppm) or LOEL (50 ppm), respectively.

A toxicokinetic study of inhaled ethylene glycol ethyl ether acetate and validation of a physiologically based pharmacokinetic model for rat and human

The solvents ethylene glycol monoethyl ether acetate (EGEEA) and ethylene glycol monoethyl ether (EGEE), at sufficiently high doses, are known to be rodent developmental toxicants, exerting their toxic effects through the action of their metabolite 2-ethoxyacetic acid (2-EAA). Thus risks associated with exposure to these compounds are best evaluated based on a measure of the internal dose of 2-EAA. The goals of the work reported here were to develop physiologically based pharmacokinetic (PBPK) models of EGEEA and EGEE for pregnant rats and humans. These models were used to identify human exposure levels (ppm in air) equivalent to the rat no observed effect level (NOEL) and lowest observed effect level (LOEL) for developmental effects (Hanley et al., 1984). We exposed pregnant Sprague-Dawley rats to concentrations of EGEEA corresponding to the NOEL and LOEL. Maternal blood, urine, and fetal tissue concentrations of EGEE and 2-EAA measured in these experiments were used to validate the rat EGEEA and EGEE models. Data collected by other researchers were used to validate the capabilities of the rodent EGEEA and EGEE models to predict the kinetics in humans. The models for estimating circulating blood concentrations of 2-EAA were considered valid based on the ability of the model to accurately predict 2-EAA concentrations in rat blood, urine, and fetal tissue. The human inhaled concentration equivalent to the rat NOEL for EGEEA (50 ppm) was predicted to be 25 ppm using the maternal blood average daily area under the curve (AUC) and 40 ppm using the maximum concentration achieved in maternal blood (C(max)). The human inhaled concentration equivalent to the rat LOEL for EGEEA (100 ppm) was determined to be 55 ppm using the maternal blood average daily AUC and 80 ppm using the maternal blood C(max).

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.

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.

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.

Physiologically based pharmacokinetics and the dermal absorption of 2-butoxyethanol vapor by humans

It has generally been assumed that the skin contributes only minor amounts to the total uptake of solvent vapors, relative to the respiratory tract. Contrary to this assumption, the widely used glycol ether solvent, 2-butoxyethanol (BE), has been reported to be more effectively absorbed through the skin (75% of the total uptake) than through the lungs of humans (Johanson and Boman, 1991, Br. J. Ind. Med. 48, 788). The possibility that the finger prick blood sampling technique used in the Johanson and Boman study was confounded by locally high concentrations of BE at the site of absorption was suggested using a previously developed PBPK model (Corley et al., 1994, Toxicol. Appl. Pharmacol. 129, 61). The current study was conducted to verify the PBPK analysis and to determine whether or not the skin was the major site for absorption of BE vapor by exposing one arm from each of six human volunteers to 50 ppm 13C2-BE vapor for 2 hr. To evaluate the potential consequences of blood sampling techniques, samples were taken from both the unexposed arm (catheter; during and after exposure) and the exposed arm (finger prick; end of the exposure only) for analysis of both BE and its major metabolite, butoxyacetic acid (BAA). Butoxyacetic acid is responsible for the hemolysis observed in toxicity studies with laboratory animals. Humans, however, are significantly less sensitive to this effect. The concentration of BE in the finger prick blood samples averaged 1500 times higher than the corresponding concentration in venous blood sampled from a catheter installed in the unexposed arm at the end of the exposure. Blood BAA levels were generally within a factor of 4 of each other for the two techniques and, therefore, was considered a better indicator of systemic absorption. Urine was collected for 24 hr and analyzed for the following metabolites found in rat metabolism studies: free and conjugated BE, BAA, ethylene glycol (EG), and glycolic acid (GA), with only BAA detected in the human urine. More importantly, urinary BAA was found to be extensively conjugated ( approximately 67%) with glutamine, confirming recent reports. These results, coupled with PBPK modeling of worst-case exposure scenarios (no clothing, 100% of the body was exposed), demonstrated that no more than 15-27% (low-to-high relative temperatures and humidities), not 75%, of the total uptake of BE could be attributed to the skin of humans during simulated 8-hr exposures to the ACGIH TLV concentration of 25 ppm. Even less of the total uptake was attributed to the skin during simulations of exercise with whole-body exposures (5-9%) or by more realistic exposures of only the arms and head (1-8%). As a result, humans are unlikely to reach hemolytic concentrations of the metabolite BAA in blood following vapor exposures to BE.

Physiologically based pharmacokinetics of 2-butoxyethanol and its major metabolite, 2-butoxyacetic acid, in rats and humans

A physiologically based pharmacokinetic model was developed to describe the disposition of 2-butoxyethanol (CAS 111-76-2) and its major metabolite, 2-butoxyacetic acid, in rats and humans. A previous human inhalation model by Johanson (Toxicol. Lett. 34, 23 (1986)) was expanded to include additional routes of exposure, physiological descriptions for rats, competing pathways for metabolism of 2-butoxyethanol, and measured partition coefficients for 2-butoxyethanol and 2-butoxyacetic acid. Simulations were compared to data gathered from rats following either intravenous infusion or oral or inhalation exposure and from humans following either inhalation or dermal exposure to 2-butoxyethanol. It was necessary to add equations for both protein binding of 2-butoxyacetic acid in blood and saturable elimination of 2-butoxyacetic acid by the kidneys to consistently describe the data. While the model predicted that rats metabolize 2-butoxyethanol and eliminate the acid metabolite faster per kilogram body weight than humans, the balance of these two processes in addition to physiological differences between species resulted in higher predicted peak blood concentrations as well as total areas under the blood concentration time curves for 2-butoxyacetic acid for rats versus humans. These species differences in kinetics coupled with the fact that human blood is significantly less susceptible than rat blood to the hemolytic effects of 2-butoxyacetic acid indicate that there is considerably less risk for hemolysis in humans as a result of exposure to 2-butoxyethanol than would have been predicted solely from standard toxicity studies with rats.

Estimating the risk of liver cancer associated with human exposures to chloroform using physiologically based pharmacokinetic modeling

A physiologically based pharmacokinetic (PB-PK) model for CHCl3 has been used to prepare estimates of the probability that human populations exposed to low levels of CHCl3 will develop liver tumors similar to those seen in rodent bioassays. The PB-PK model for CHCl3 was based on a model reported earlier by Corley et al. (1990), but this model differed from that of Corley et al. in that it was also capable of describing a pharmacodynamic endpoint: induction of cytotoxicity in the liver of CHCl3-exposed animals produced by reactive metabolites of CHCl3. Pharmacodynamic descriptions in this model were derived from experimental measurements of cell replication ([3H]thymidine incorporation) as well as from quantitative histopathology in the liver of rats and mice. Two different approaches were used for hazard evaluation: (1) a "Safety Factor" approach based on no observed effect levels for liver tumors. and (2) calculation of lower confidence limits on risk-specific doses with the GLOBAL83 computer program. In each case, cytotoxicity produced by reactive CHCl3 metabolites was used as the measure of "dose" to the liver. The Safety Factor approach suggested that continuous exposure of human populations to concentrations of CHCl3 less than 2840 ppb in air or 13,900 ppb in water would not be likely to significantly increase the risk of developing liver tumors. The second approach suggested a "plausible upper 95% confidence limit" of 1 x 10(-5) for lifetime excess cancer risk for human populations continuously exposed to 2200 or 13,100 ppb CHCl3 in air or water, respectively.

Development of a physiologically based pharmacokinetic model for chloroform

A physiologically based pharmacokinetic model describing the disposition of chloroform in mice, rats, and humans was developed. This model was designed to facilitate extrapolations from high doses, such as those used in chronic rodent studies, to low doses that humans may be exposed to in the workplace or the environment. Kinetic constants for mice and rats were derived from in vivo experiments. Enzymatic studies conducted with samples of rodent and human tissues provided a rational basis for estimating human in vivo metabolic rate constants. Incorporation of physiological descriptions of the processes of absorption, distribution, metabolism, and excretion allowed extrapolation between different routes of exposure as well. The model was validated by comparing model predictions with experimental data gathered in mice, rats, and humans after inhalation, oral, or intraperitoneal administration of chloroform. Consistent with previous reports, the metabolic activation of chloroform to toxic intermediates was shown to occur most rapidly in the mouse, less rapidly in the rat, and most slowly in humans. Estimates of the "delivered dose" of chloroform metabolites to internal organs sensitive to chloroform toxicity were calculated. This model may be used to develop refined dose estimates for human populations exposed to low levels of chloroform in the environment.

Chlorpyrifos: a 13-week nose-only vapor inhalation study in Fischer 344 rats

Fischer 344 rats were exposed by the nose-only inhalation route to chlorpyrifos vapors at concentrations of 0, 5.2, 10.3, or 20.6 ppb, 6 hr/day, 5 days/week for 13 weeks. The exposure concentrations were limited by the low vapor pressure of chlorpyrifos (theoretical maximum vapor concentration of 25 ppb at 25 degrees C). No treatment-related signs of toxicity or changes in body weights were detected during the course of the study. Urinalysis, hematology, clinical chemistry, organ weights, gross pathologic, and histopathologic evaluations were performed at the end of the study with no treatment-related effects observed. In addition, no differences from controls were noted in plasma, red blood cell, or brain cholinesterase activities. The results of this study indicate that the no-observed-effect level for chlorpyrifos vapor was the highest attainable concentration, 20.6 ppb, in male and female Fischer 344 rats.

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.

Rapid thin layer chromatographic method for determination of zearalenone and zearalenol in grains and animal feeds

A rapid and inexpensive method has been developed for the analysis of zearalenone and zearalenol in grains and animal feeds. The method involves extraction with 75% methanol, precipitation of pigments with lead acetate, and defatting with petroleum ether. The mycotoxins are subsequently partitioned into toluene-ethyl acetate, chromatographed on high performance thin layer chromatographic plates, and detected after treatment with Fast Violet B salt solution. Sensitivity of the method is better than 80 ng/g for zearalenone and 200 ng/g for zearalenol. Ten samples can be completed in less than 2 h. The method is applicable for zearalenone in corn, wheat, barley, millet, and swine feeds.