Screening-level risk assessment for styrene-acrylonitrile (SAN) trimer detected in soil and groundwater.

A screening-level risk assessment was conducted for styrene-acrylonitrile (SAN) Trimer detected at the Reich Farm Superfund site in Toms River, NJ. Consistent with a screening-level approach, on-site and off-site exposure scenarios were evaluated using assumptions that are expected to overestimate actual exposures and hazards at the site. Environmental sampling data collected for soil and groundwater were used to estimate exposure point concentrations. Several exposure scenarios were evaluated to assess potential on-site and off-site exposures, using parameter values for exposures to soil (oral, inhalation of particulates, and dermal contact) and groundwater (oral, dermal contact) to reflect central tendency exposure (CTE) and reasonable maximum exposure (RME) conditions. Three reference dose (RfD) values were derived for SAN Trimer for short-term, subchronic, and chronic exposures, based upon its effects on the liver in exposed rats. Benchmark (BMD) methods were used to assess the relationship between exposure and response, and to characterize appropriate points of departure (POD) for each RfD. An uncertainty factor of 300 was applied to each POD to yield RfD values of 0.1, 0.04, and 0.03 mg/kg-d for short-term, subchronic, and chronic exposures, respectively. Because a chronic cancer bioassay for SAN Trimer in rats (NTP 2011a) does not provide evidence of carcinogenicity, a cancer risk assessment is not appropriate for this chemical. Potential health hazards to human health were assessed using a hazard index (HI) approach, which considers the ratio of exposure dose (i.e., average daily dose, mg/kg-d) to toxicity dose (RfD, mg/kg-d) for each scenario. All CTE and RME HI values are well below 1 (where the average daily dose is equivalent to the RfD), indicating that there is no concern for potential noncancer effects in exposed populations even under the conservative assumptions of this screening-level assessment.

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.

Acrylamide: Consideration of species differences and nonlinear processes in estimating risk and safety for human ingestion.

Acrylamide in cooked foods results in wide-spread, low-level human exposure. Potential risks from dietary intake remain unclear due to apparent conflicting results from cancer bioassays conducted in rats that reported tumors and epidemiology studies that are suggestive but provide little or no evidence of increased cancer. Risk estimation often includes two common assumptions: (1) tumor response rates in test species can be extrapolated systematically to estimate human response rates and (2) tumor rates observed following high-dose exposures can be linearly extrapolated to predict response rates following low-dose exposures. The validity of these assumptions was evaluated for acrylamide based upon the examination of relevant toxicokinetic and toxicodynamic differences between humans and rats, including sources of nonlinearity that modify high to low dose extrapolation of cancer incidence. Important species differences and sources of nonlinearity are identified, and recommendations for addressing them within the quantitative framework of a PBTK/TD model are discussed. These differences are likely to estimate risk levels up to several orders of magnitude lower in humans than in rats. Quantitative inclusion of these TK/TD factors will more closely estimate actual human cancer risk derived from high-dose rodent studies, since detoxification processes for acrylamide and glycidamide appear adequately protective against toxicity from human dietary doses.

Derivation of noncancer reference values for acrylonitrile.

Dose-response assessments were conducted for the noncancer effects of acrylonitrile (AN) for the purposes of deriving subchronic and chronic oral reference dose (RfD) and inhalation reference concentration (RfC) values. Based upon an evaluation of available toxicity data, the irritation and neurological effects of AN were determined to be appropriate bases for deriving reference values. A PBPK model, which describes the toxicokinetics of AN and its metabolite 2-cyanoethylene oxide (CEO) in both rats and humans, was used to assess the dose-response data in terms of an internal dose measure for the oral RfD values, but could not be used in deriving the inhalation RfC values. Benchmark dose (BMD) methods were used to derive all reference values. Where sufficient information was available, data-derived uncertainty factors were applied to the points of departure determined by BMD methods. From this assessment, subchronic and chronic oral RfD values of 0.5 and 0.05 mg/kg/day, respectively, were derived. Similarly, subchronic and chronic inhalation RfC values of 0.1 and 0.06 mg/m(3), respectively, were derived. Confidence in the reference values derived for AN was considered to be medium to high, based upon a consideration of the confidence in the key studies, the toxicity database, dosimetry, and dose-response modeling.

Physiologically based modeling of the inhalation pharmacokinetics of ethylbenzene in B6C3F1 mice.

A physiologically based pharmacokinetic (PBPK) model was developed for inhaled ethylbenzene (EB) in B6C3F1 mice. The mouse physiological parameters were obtained from the literature, but the blood:air and tissue:air partition coefficients were determined by vial equilibration technique. The maximal velocity for hepatic metabolism (Vmax) obtained from a previously published rat study was increased by a factor of approximately 3 to account for enzyme induction during repeated exposures. The Michaelis affinity constant (Km) for hepatic metabolism of EB, obtained from a previously published rat PBPK modeling study, was kept unchanged during single and repeated exposure scenarios. Hepatic metabolism alone could not adequately describe the clearance of EB from mouse blood. Additional metabolism was assumed to be localized in the lung. The parameters for pulmonary metabolism were obtained by optimization of PBPK model fits to kinetic data collected following exposures to 75-1000 ppm. The PBPK model successfully predicted all available blood and tissue concentration data in mice exposed to 75 or 750 ppm EB. Overall, the results indicate that the rate of EB clearance is markedly higher in B6C3F1 mice than rats or humans and exceeds the hepatic metabolism capacity. Available biochemical evidence is consistent with a significant role for pulmonary metabolism; however, the extent to which the extrahepatic metabolism is localized in the lung is unclear. Overall, the PBPK model developed for the mouse adequately simulated the blood and tissue kinetics of EB by accounting for its high rate of clearance.

Cancer dose--response assessment for acrylonitrile based upon rodent brain tumor incidence: use of epidemiologic, mechanistic, and pharmacokinetic support for nonlinearity.

A cancer dose-response assessment was conducted for acrylonitrile (AN) using updated information on mechanism of action, epidemiology, toxicity, and pharmacokinetics. Although more than 10 chronic bioassays indicate that AN produces multiple tumors in rats and mice, a number of large, well-conducted epidemiology studies provide no evidence of a causal association between AN exposure and cancer mortality of any type. The epidemiological data include early industry exposures that are far higher than occur today and that approach or exceed levels found to be tumorigenic in animals. Despite the absence of positive findings in the epidemiology data, a dose-response assessment was conducted for AN based on brain tumors in rats. Mechanistic studies implicate the involvement of oxidative stress in rat brain due to a metabolite (2-cyanoethylene oxide or CEO, cyanide), but do not conclusively rule out a potential role for the direct genotoxicity of CEO. A PBPK model was used to predict internal doses (peak CEO in brain) for 12 data sets, which were pooled together to provide a consistent characterization of the dose-response relationship for brain tumor incidence in the rat. The internal dose corresponding to a 5% increase in extra risk (ED 05=0.017 mg/L brain) and its lower confidence limit (LED 05=0.014 mg/L brain) was used as the point of departure. The weight-of-evidence supports the use of a nonlinear extrapolation for the cancer dose-response assessment. A quantitative comparison of the epidemiology exposure-response data (lung and brain cancer mortality) to the rat brain tumor data in terms of internal dose adds to the confidence in the nonlinear extrapolation. Uncertainty factors of 200 and 220 (for the oral and inhalation routes, respectively) were applied to the LED 05 to account for interspecies variation, intraspecies variation, and the severity of the response measure. Accordingly, oral doses below 0.009 mg/kg-day and air concentrations below 0.1mg/m(3) are not expected to pose an appreciable risk to human populations exposed to AN.

Using physiologically-based pharmacokinetic modeling to address nonlinear kinetics and changes in rodent physiology and metabolism due to aging and adaptation in deriving reference values for propylene glycol methyl ether and propylene glycol methyl ether acetate.

Reference values, including an oral reference dose (RfD) and an inhalation reference concentration (RfC), were derived for propylene glycol methyl ether (PGME), and an oral RfD was derived for its acetate (PGMEA). These values were based on transient sedation observed in F344 rats and B6C3F1 mice during a two-year inhalation study. The dose-response relationship for sedation was characterized using internal dose measures as predicted by a physiologically-based pharmacokinetic (PBPK) model for PGME and its acetate. PBPK modeling was used to account for changes in rodent physiology and metabolism due to aging and adaptation, based on data collected during Weeks 1, 2, 26, 52, and 78 of a chronic inhalation study. The peak concentration of PGME in richly perfused tissues (i.e., brain) was selected as the most appropriate internal dose measure based on a consideration of the mode of action for sedation and similarities in tissue partitioning between brain and other richly perfused tissues. Internal doses (peak tissue concentrations of PGME) were designated as either no-observed-adverse-effect levels (NOAELs) or lowest-observed-adverse-effect levels (LOAELs) based on the presence or the absence of sedation at each time point, species, and sex in the two-year study. Distributions of the NOAEL and LOAEL values expressed in terms of internal dose were characterized using an arithmetic mean and standard deviation, with the mean internal NOAEL serving as the basis for the reference values, which was then divided by appropriate uncertainty factors. Where data were permitting, chemical-specific adjustment factors were derived to replace default uncertainty factor values of 10. Nonlinear kinetics, which was predicted by the model in all species at PGME concentrations exceeding 100 ppm, complicate interspecies, and low-dose extrapolations. To address this complication, reference values were derived using two approaches that differ with respect to the order in which these extrapolations were performed: (1) default approach of interspecies extrapolation to determine the human equivalent concentration (PBPK modeling) followed by uncertainty factor application, and (2) uncertainty factor application followed by interspecies extrapolation (PBPK modeling). The resulting reference values for these two approaches are substantially different, with values from the latter approach being seven-fold higher than those from the former approach. Such a striking difference between the two approaches reveals an underlying issue that has received little attention in the literature regarding the application of uncertainty factors and interspecies extrapolations to compounds where saturable kinetics occur in the range of the NOAEL. Until such discussions have taken place, reference values based on the former approach are recommended for risk assessments involving human exposures to PGME and PGMEA.

Addressing nonlinearity in the exposure-response relationship for a genotoxic carcinogen: cancer potency estimates for ethylene oxide.

Ethylene oxide (EO) has been identified as a carcinogen in laboratory animals. Although the precise mechanism of action is not known, tumors in animals exposed to EO are presumed to result from its genotoxicity. The overall weight of evidence for carcinogenicity from a large body of epidemiological data in the published literature remains limited. There is some evidence for an association between EO exposure and lympho/hematopoietic cancer mortality. Of these cancers, the evidence provided by two large cohorts with the longest follow-up is most consistent for leukemia. Together with what is known about human leukemia and EO at the molecular level, there is a body of evidence that supports a plausible mode of action for EO as a potential leukemogen. Based on a consideration of the mode of action, the events leading from EO exposure to the development of leukemia (and therefore risk) are expected to be proportional to the square of the dose. In support of this hypothesis, a quadratic dose-response model provided the best overall fit to the epidemiology data in the range of observation. Cancer dose-response assessments based on human and animal data are presented using three different assumptions for extrapolating to low doses: (1) risk is linearly proportionate to dose; (2) there is no appreciable risk at low doses (margin-of-exposure or reference dose approach); and (3) risk below the point of departure continues to be proportionate to the square of the dose. The weight of evidence for EO supports the use of a nonlinear assessment. Therefore, exposures to concentrations below 37 microg/m3 are not likely to pose an appreciable risk of leukemia in human populations. However, if quantitative estimates of risk at low doses are desired and the mode of action for EO is considered, these risks are best quantified using the quadratic estimates of cancer potency, which are approximately 3.2- to 32-fold lower, using alternative points of departure, than the linear estimates of cancer potency for EO. An approach is described for linking the selection of an appropriate point of departure to the confidence in the proposed mode of action. Despite high confidence in the proposed mode of action, a small linear component for the dose-response relationship at low concentrations cannot be ruled out conclusively. Accordingly, a unit risk value of 4.5 x 10(-8) (microg/m3)(-1) was derived for EO, with a range of unit risk values of 1.4 x 10(-8) to 1.4 x 10(-7) (microg/m3)(-1) reflecting the uncertainty associated with a theoretical linear term at low concentrations.

Assessing the dose-dependency of allometric scaling performance using physiologically based pharmacokinetic modeling.

The performance of allometric scaling of dose as a power of body weight under a variety of extrapolation conditions with respect to species, route, exposure intensity, and mechanism/mode of action, remains untested in many cases. In this paper, animal-human internal dose ratio comparisons have been developed for 12 chemicals (benzene, carbon tetrachloride, chloroform, diisopropylfluorophosphate, ethanol, ethylene oxide, methylene chloride, methylmercury, styrene, tetrachloroethene, trichloroethene, and vinyl chloride). This group of predominantly volatile and lipophilic chemicals was selected on the basis that their kinetics have been well-studied and can be predicted in mice, rats, and humans using physiologically based pharmacokinetic (PBPK) models. PBPK model predictions were compared to the allometric scaling predictions for interspecies extrapolation. Recommendations for the application of the allometric scaling are made with reference to internal dose measure (mode of action) and concentration level. The results of this assessment generally support the use of scaling factors recommended in the published literature, which includes scaling factors of 1.0 for risk assessments in which toxicity is attributed to the parent chemical or stable metabolite, and -0.75 for dose-response assessments in which toxicity is attributed to the formation of a reactive metabolite from an inhaled compound. A scaling factor of 0.75 is recommended for dose-response assessments of orally administered compounds in which toxicity is attributed to the parent chemical or stable metabolite and 1.0 for risk assessments in which toxicity is attributed to the formation of a reactive metabolite from a compound administered by the oral route. A dose-dependency in the results suggests that the scaling factors appropriate at high exposures may differ from those at low exposures, primarily due to the impact of saturable metabolism.

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.

Development of a preliminary physiologically based toxicokinetic (PBTK) model for 1,3-butadiene risk assessment.

Potential health effects of human exposure to 1,3-butadiene (BD) are of concern due to the use of BD in industry and its low-level presence throughout the environment. Physiologically based toxicokinetic (PBTK) models of BD in rodents have been developed by multiple research groups in an effort to explain species differences in toxicity (especially carcinogenic potency) through toxicokinetics. PBTK modeling of dose metrics related to a non-cancer endpoint, ovotoxicity in experimental animals, was conducted. The cumulative area under the blood concentration vs. time curve (AUC) for the metabolite diepoxybutane (butadiene diepoxide, DEB) was found to be consistent with ovotoxicity in mice and rats exposed to BD by inhalation or epoxybutene (butadiene monoepoxide, EB) or DEB by intraperitoneal injection. This suggests that cumulative DEB AUC may also be an appropriate metric for possible human risk. A preliminary human PBTK model was assembled for the eventual assessment of reproductive risk to humans and for prioritizing the determination of model parameters. The preliminary model accurately predicted published data on exhaled breath BD concentrations in a human volunteer exposed to BD by inhalation. The fit was relatively insensitive to the rate constant for BD epoxidation. Sensitivity analyses were conducted on this human PBTK model. Using a range of published rate constants, human blood DEB was found to be sensitive to rates of epoxidation of EB to DEB and hydrolysis of EB and DEB, but not BD epoxidation. Because of the large ranges of rates measured in vitro for these reactions, different combinations of in-vitro rates produce varying predictions of blood DEB concentration. Thus, validation of a human PBTK model with human biomonitoring data will be essential to produce a PBTK model that can be applied to risk assessment.

Improving cancer dose-response characterization by using physiologically based pharmacokinetic modeling: an analysis of pooled data for acrylonitrile-induced brain tumors to assess cancer potency in the rat.

Historically, U.S. regulators have derived cancer slope factors by using applied dose and tumor response data from a single key bioassay or by averaging the cancer slope factors of several key bioassays. Recent changes in U.S. Environmental Protection Agency (EPA) guidelines for cancer risk assessment have acknowledged the value of better use of mechanistic data and better dose-response characterization. However, agency guidelines may benefit from additional considerations presented in this paper. An exploratory study was conducted by using rat brain tumor data for acrylonitrile (AN) to investigate the use of physiologically based pharmacokinetic (PBPK) modeling along with pooling of dose-response data across routes of exposure as a means for improving carcinogen risk assessment methods. In this study, two contrasting assessments were conducted for AN-induced brain tumors in the rat on the basis of (1) the EPA's approach, the dose-response relationship was characterized by using administered dose/concentration for each of the key studies assessed individually; and (2) an analysis of the pooled data, the dose-response relationship was characterized by using PBPK-derived internal dose measures for a combined database of ten bioassays. The cancer potencies predicted for AN by the contrasting assessments are remarkably different (i.e., risk-specific doses differ by as much as two to four orders of magnitude), with the pooled data assessments yielding lower values. This result suggests that current carcinogen risk assessment practices overestimate AN cancer potency. This methodology should be equally applicable to other data-rich chemicals in identifying (1) a useful dose measure, (2) an appropriate dose-response model, (3) an acceptable point of departure, and (4) an appropriate method of extrapolation from the range of observation to the range of prediction when a chemical's mode of action remains uncertain.

A physiological toxicokinetic model for exogenous and endogenous ethylene and ethylene oxide in rat, mouse, and human: formation of 2-hydroxyethyl adducts with hemoglobin and DNA.

Ethylene (ET) is a gaseous olefin of considerable industrial importance. It is also ubiquitous in the environment and is produced in plants, mammals, and humans. Uptake of exogenous ET occurs via inhalation. ET is biotransformed to ethylene oxide (EO), which is also an important volatile industrial chemical. This epoxide forms hydroxyethyl adducts with macromolecules such as hemoglobin and DNA and is mutagenic in vivo and in vitro and carcinogenic in experimental animals. It is metabolically eliminated by epoxide hydrolase and glutathione S-transferase and a small fraction is exhaled unchanged. To estimate the body burden of EO in rodents and human resulting from exposures to EO and ET, we developed a physiological toxicokinetic model. It describes uptake of ET and EO following inhalation and intraperitoneal administration, endogenous production of ET, enzyme-mediated oxidation of ET to EO, bioavailability of EO, EO metabolism, and formation of 2-hydroxyethyl adducts of hemoglobin and DNA. The model includes compartments representing arterial, venous, and pulmonary blood, liver, muscle, fat, and richly perfused tissues. Partition coefficients and metabolic parameters were derived from experimental data or published values. Model simulations were compared with a series of data collected in rodents or humans. The model describes well the uptake, elimination, and endogenous production of ET in all three species. Simulations of EO concentrations in blood and exhaled air of rodents and humans exposed to EO or ET were in good agreement with measured data. Using published rate constants for the formation of 2-hydroxyethyl adducts with hemoglobin and DNA, adduct levels were predicted and compared with values reported. In humans, predicted hemoglobin adducts resulting from exposure to EO or ET are in agreement with measured values. In rodents, simulated and measured DNA adduct levels agreed generally well, but hemoglobin adducts were underpredicted by a factor of 2 to 3. Obviously, there are inconsistencies between measured DNA and hemoglobin adduct levels.

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).

Development of a physiologically based pharmacokinetic model of 2-methoxyethanol and 2-methoxyacetic acid disposition in pregnant rats.

An accurate description of developing embryos' exposure to a xenobiotic is a desirable component of mechanism-based risk assessments for humans exposed to potential developmental toxicants during pregnancy. 2-Methoxyethanol (2-ME), a solvent used in the manufacture of semiconductors, is embryotoxic and teratogenic in all species tested including nonhuman primates. 2-Methoxyacetic acid (2-MAA) is the primary metabolite of 2-ME and the proximate embryotoxic agent. The objective of the work described here was to adapt an existing physiologically based pharmacokinetic (PBPK) model for 2-ME and 2-MAA kinetics during midorganogenesis in mice to rats on gestation days (GD) 13 and 15. Blood and tissue data were analyzed using the extrapolated PBPK model that was modified to simulate 2-ME and 2-MAA kinetics in maternal plasma and total embryo tissues in pregnant rats. The original mouse model was simplified by combining the embryos and placenta with the richly perfused tissue compartment. The model includes a description of the growth of the developing embryo and changes in the physiology of the dam during pregnancy. Biotransformation pathways of 2-ME to either ethylene glycol (EG) or to 2-MAA were described as first-order processes based on the data collected from rats by Green et al., (Occup. Hyg. 2, 67-75, 1996). Tissue partition coefficients (PCs) for 2-ME and 2-MAA were determined for a variety of maternal tissues and the embryos. Model simulations closely reflected the biological measurement of 2-ME and 2-MAA concentrations in blood and embryo tissue following gavage or iv administration of 2-ME or 2-MAA. The PBPK model for rats as described here is well suited for extrapolation to pregnant women and for assessment of 2-MAA dosimetry under various conditions of possible human exposure to 2-ME.

Kinetic characterization of CYP2E1 inhibition in vivo and in vitro by the chloroethylenes.

Trans- and cis-1,2-dichloroethylene (DCE) isomers inhibit their own metabolism in vivo by inactivation of the metabolizing enzyme, presumably the cytochrome P450 isoform, CYP2E1. In this study, we examined cytochrome P450 isoform-specific inhibition by three chloroethylenes, cis-DCE, trans-DCE, and trichloroethylene (TCE), and evaluated several kinetic mechanisms of enzyme inhibition with physiological models of inhibition. Trans-DCE was more potent than cis-DCE, and both were much more effective than TCE in inhibiting CYP2E1. The kinetics of in vitro loss of p-nitrophenol hydroxylase (pNP-OH) activity (a marker of CYP2E1) in microsomal incubations and of the in vivo gas uptake results were most consistent with a mechanism in which inhibition of the metabolizing enzyme (CYP2E1) was presumed to be related to interaction of a reactive DCE metabolite with remaining substrate-bound, active CYP2E1. The kinetics of inhibition by TCE, a weak inhibitor in vitro, were very different from that of the dichloroethylenes. With TCE, parent compound concentrations influenced enzyme loss. Trans-DCE was a more potent inhibitor of CYP2E1 than cis-DCE based on both in vivo and in vitro studies. Quantitative differences in the inhibitory properties of the 1,2-DCE isomers may be due to the different stability of epoxides formed from bioactivation by CYP2E1. Epoxide intermediates of DCE metabolism, reacting by water addition, would yield dialdehyde, a potent cross-linking reagent.

Human ingestion of chromium (VI) in drinking water: pharmacokinetics following repeated exposure.

Regulatory agencies have established safe drinking water concentrations for hexavalent chromium [Cr(VI)] based in part on the presumed capability of human gastric juices to rapidly reduce Cr(VI) to nontoxic trivalent chromium [Cr(III)] prior to systemic absorption. This study examines dose-related pharmacokinetics in humans following repeated oral exposure to Cr(VI) in drinking water. In particular, we sought to examine whether plausible drinking water exposures to Cr(VI) caused a sustained increase in red blood cell chromium levels, a specific marker for systemic uptake of Cr(VI). Adult male volunteers ingested a liter (in three volumes of 333 ml, at approximate 6-hr intervals) of deionized water containing Cr(VI) concentrations ranging from 0.1 to 10.0 mg/liter. Samples of urine, plasma, and red blood cells were collected and analyzed for chromium. A dose-related increase in urinary chromium excretion was observed in all volunteers. Red blood cell and plasma chromium concentrations became elevated in certain individuals at the highest doses. The RBC chromium profiles suggest that the ingested Cr(VI) was reduced to Cr(III) before entering the bloodstream, since the chromium concentration in the RBCs dropped rapidly postexposure. These findings suggest that the human gastrointestinal tract has the capacity to reduce ingested Cr(VI) following ingestion of up to 1 liter of water containing 10.0 mg/liter of Cr(VI), which is consistent with USEPA's position that the Cr(VI) drinking water standard of 0.10 mg Cr(VI)/liter is below the reductive capacity of the stomach.

Refinement and verification of the physiologically based dosimetry description for acrylonitrile in rats.

The physiologically based dosimetry description for acrylonitrile (ACN) and its mutagenic epoxide metabolite 2-cyanoethylene oxide (CEO) in F-344 rats (M. L. Gargas, M. E. Anderson, S.K.O. Teo, R. Batra, T. R. Fennell, and G. L. Kedderis, 1995, Toxicol. Appl. Pharmacol. 134, 185-194) has been refined to include a physiological stomach compartment and the reactions of ACN with tissue glutathione (GSH). The second-order rate constant for reaction of ACN and GSH at pH 7.3 was measured and included in the dosimetry description. Metabolic parameters for ACN and CEO were estimated from oral bolus pharmacokinetic studies and previously obtained iv bolus data (3.4, 47, 55, or 84 mg ACN/kg). Rats were given bolus oral doses of 3, 10, or 30 mg ACN/kg in water, and blood samples were collected at selected time points. ACN and CEO blood concentrations were determined by gas chromatography. The brain and liver concentrations of ACN and CEO were also measured after 10 mg ACN/kg po. ACN elimination from blood was described by saturable P450 epoxidation (Vmax of 5.0 mg/hr/kg and K(M) of 1.5 mg/liter) and first-order GSH conjugation (73 hr(-1)/kg). CEO elimination was described by first-order GSH conjugation (500 hr(-1)/kg). The pharmacokinetic data were well simulated, although CEO blood concentrations after bolus oral dosing were somewhat overestimated. Sensitivity analysis of the dosimetry description indicated that the inhalation exposure route was much more sensitive to changes in metabolic and physiological parameters than either the iv or oral bolus routes. Therefore, inhalation pharmacokinetic data were obtained and compared to simulations of the dosimetry description. Rats were exposed to 186, 254, or 291 ppm ACN for 3 hr. ACN and CEO concentrations were measured in blood, brain, and liver at selected postexposure time points. The dosimetry description accurately simulated the ACN inhalation pharmacokinetic data, providing verification of the parameter estimates. The verified rat dosimetry description for ACN and CEO will be used as the basis for development of a dosimetry description for ACN in people.

PBPK modeling/Monte Carlo simulation of methylene chloride kinetic changes in mice in relation to age and acute, subchronic, and chronic inhalation exposure.

During a 2-year chronic inhalation study on methylene chloride (2000 or 0 ppm; 6 hr/day, 5 days/week), gas-uptake pharmacokinetic studies and tissue partition coefficient determinations were conducted on female B6C3F1, mice after 1 day, 1 month, 1 year, and 2 years of exposure. Using physiologically based pharmacokinetic (PBPK) modeling coupled with Monte Carlo simulation and bootstrap resampling for data analyses, a significant induction in the mixed function oxidase (MFO) rate constant (Vmaxc) was observed at the 1-day and 1-month exposure points when compared to concurrent control mice while decreases in glutathione S-transferase (GST) rate constant (Kfc) were observed in the 1-day and 1-month exposed mice. Within exposure groups, the apparent Vmaxc maintained significant increases in the 1-month and 2-year control groups. Although the same initial increase exists in the exposed group, the 2-year Vmaxc is significantly smaller than the 1-month group (p < 0.001). Within group differences in median Kfc values show a significant decrease in both 1-month and 2-year groups among control and exposed mice (p < 0.001). Although no changes in methylene chloride solubility as a result of prior exposure were observed in blood, muscle, liver, or lung, a marginal decrease in the fat:air partition coefficient was found in the exposed mice at p = 0.053. Age related solubility differences were found in muscle:air, liver:air, lung:air, and fat:air partition coefficients at p < 0.001, while the solubility of methylene chloride in blood was not affected by age (p = 0.461). As a result of this study, we conclude that age and prior exposure to methylene chloride can produce notable changes in disposition and metabolism and may represent important factors in the interpretation for toxicologic data and its application to risk assessment.