Origin of Groundwater Arsenic in a Rural Pleistocene Aquifer in Bangladesh Depressurized by Distal Municipal Pumping.

Across South Asia, millions of villagers have reduced their exposure to high-arsenic (As) groundwater by switching to low-As wells. Isotopic tracers and flow modeling are used in this study to understand the groundwater flow system of a semi-confined aquifer of Pleistocene (>10 kyr) age in Bangladesh that is generally low in As but has been perturbed by massive pumping at a distance of about 25 km for the municipal water supply of Dhaka. A 10- to 15-m-thick clay aquitard caps much of the intermediate aquifer (>40- to 90-m depth) in the 3-km2 study area, with some interruptions by younger channel sand deposits indicative of river scouring. Hydraulic heads in the intermediate aquifer below the clay-capped areas are 1-2 m lower than in the high-As shallow aquifer above the clay layer. In contrast, similar heads in the shallow and intermediate aquifer are observed where the clay layer is missing. The head distribution suggests a pattern of downward flow through interruptions in the aquitard and lateral advection from the sandy areas to the confined portion of the aquifer. The interpreted flow system is consistent with 3H-3He ages, stable isotope data, and groundwater flow modeling. Lateral flow could explain an association of elevated As with high methane concentrations within layers of gray sand below certain clay-capped portions of the Pleistocene aquifer. An influx of dissolved organic carbon from the clay layer itself leading to a reduction of initially orange sands has also likely contributed to the rise of As.

Early diastolic septal movement in patients with myocarditis.

AIM: To evaluate early diastolic septal relaxation as a parameter in the diagnostic workup via cardiovascular magnetic resonance imaging (CMRI) in patients with myocarditis. MATERIALS AND METHODS: Early diastolic septal movement was evaluated (EDS) prospectively via frame-by-frame analysis in 255 consecutive patients with presenting signs of myocarditis and in 64 controls matched 4:1 for gender and age. ECG-triggered, T2-weighted, fast spin echo triple inversion recovery sequences and late gadolinium enhancement were obtained, as well as left ventricular (LV) function and dimensions in patients and controls. RESULTS: EDS was detected in 66.7% of the patients and 18.7% of the controls (p<0.001). Sensitivity was 69.4% and specificity 79.7%. Patients with EDS had a significant lower LV ejection fraction (LV-EF) of 61.1±0.6% and significant higher end-diastolic volume (EDV) of 158.5±2.7 ml than in patients without EDS (LV-EF 65.3±0.9%, p=0.0001; EDV 148.4±3.9 ml, p=0.04). A significant negative correlation was observed between LV-EF and EDS in patients, and a lower LV-EF correlated with a more frequent occurrence of EDS (r=-0.24, p=0.0001). Scar tissue was also more frequent in patients than controls (63.1% and 7.8%, p=0.007). CONCLUSIONS: EDS is a parameter obtained non-invasively by CMRI and is present in a high percentage of patients with myocarditis. Cardiac functional parameters are significantly altered in patients with EDS. EDS is a feasible parameter that can play an important role in the diagnosis of myocarditis.

Using physiologically based pharmacokinetic modeling and benchmark dose methods to derive an occupational exposure limit for N-methylpyrrolidone.

The developmental effects of NMP are well studied in Sprague-Dawley rats following oral, inhalation, and dermal routes of exposure. Short-term and chronic occupational exposure limit (OEL) values were derived using an updated physiologically based pharmacokinetic (PBPK) model for NMP, along with benchmark dose modeling. Two suitable developmental endpoints were evaluated for human health risk assessment: (1) for acute exposures, the increased incidence of skeletal malformations, an effect noted only at oral doses that were toxic to the dam and fetus; and (2) for repeated exposures to NMP, changes in fetal/pup body weight. Where possible, data from multiple studies were pooled to increase the predictive power of the dose-response data sets. For the purposes of internal dose estimation, the window of susceptibility was estimated for each endpoint, and was used in the dose-response modeling. A point of departure value of 390 mg/L (in terms of peak NMP in blood) was calculated for skeletal malformations based on pooled data from oral and inhalation studies. Acceptable dose-response model fits were not obtained using the pooled data for fetal/pup body weight changes. These data sets were also assessed individually, from which the geometric mean value obtained from the inhalation studies (470 mg*hr/L), was used to derive the chronic OEL. A PBPK model for NMP in humans was used to calculate human equivalent concentrations corresponding to the internal dose point of departure values. Application of a net uncertainty factor of 20-21, which incorporates data-derived extrapolation factors, to the point of departure values yields short-term and chronic occupational exposure limit values of 86 and 24 ppm, respectively.

Recharge of low-arsenic aquifers tapped by community wells in Araihazar, Bangladesh, inferred from environmental isotopes.

More than 100,000 community wells have been installed in the 150-300 m depth range throughout Bangladesh over the past decade to provide low-arsenic drinking water (<10 µg/L As), but little is known about how aquifers tapped by these wells are recharged. Within a 25 km2 area of Bangladesh east of Dhaka, groundwater from 65 low-As wells in the 35-240 m depth range was sampled for tritium (3H), oxygen and hydrogen isotopes of water (18O/16O and 2H/1H), carbon isotope ratios in dissolved inorganic carbon (DIC, 14C/12C and 13C/12C), noble gases, and a suite of dissolved constituents, including major cations, anions, and trace elements. At shallow depths (<90 m), 24 out of 42 wells contain detectable 3H of up to 6 TU, indicating the presence of groundwater recharged within 60 years. Radiocarbon (14C) ages in DIC range from modern to 10 kyr. In the 90-240 m depth range, however, only 5 wells shallower than 150 m contain detectable 3H (<0.3 TU) and 14C ages of DIC cluster around 10 kyr. The radiogenic helium (4He) content in groundwater increases linearly across the entire range of 14C ages at a rate of 2.5×10-12 ccSTP 4He g-1 yr-1. Within the samples from depths >90 m, systematic relationships between 18O/16O, 2H/1H, 13C/12C and 14C/12C, and variations in noble gas temperatures, suggest that changes in monsoon intensity and vegetation cover occurred at the onset of the Holocene, when the sampled water was recharged. Thus, the deeper low-As aquifers remain relatively isolated from the shallow, high-As aquifer.

An evaluation of in vivo models for toxicokinetics of hexavalent chromium in the stomach.

Hexavalent chromium (Cr6) is a drinking water contaminant that has been detected in most of the water systems throughout the United States. In 2-year drinking water bioassays, the National Toxicology Program (NTP) found clear evidence of carcinogenic activity in male and female rats and mice. Because reduction of Cr6 to trivalent chromium (Cr3) is an important detoxifying step in the gastrointestinal (GI) tract prior to systemic absorption, models have been developed to estimate the extent of reduction in humans and animals. The objective of this work was to use a revised model of ex vivo Cr6 reduction kinetics in gastric juice to analyze the potential reduction kinetics under in vivo conditions for mice, rats and humans. A published physiologically-based pharmacokinetic (PBPK) model was adapted to incorporate the new reduction model. This paper focuses on the toxicokinetics of Cr6 in the stomach compartment, where most of the extracellular Cr6 reduction is believed to occur in humans. Within the range of doses administered by the NTP bioassays, neither the original nor revised models predict saturation of stomach reducing capacity to occur in vivo if applying default parameters. However, both models still indicate that mice exhibit the lowest extent of reduction in the stomach, meaning that a higher percentage of the Cr6 dose may escape stomach reduction in that species. Similarly, both models predict that humans exhibit the highest extent of reduction at low doses.

Physiologically based pharmacokinetic model use in risk assessment--Why being published is not enough.

A panel of experts in physiologically based pharmacokinetic (PBPK) modeling and relevant quantitative methods was convened to describe and discuss model evaluation criteria, issues, and choices that arise in model application and computational tools for improving model quality for use in human health risk assessments (HHRAs). Although publication of a PBPK model in a peer-reviewed journal is a mark of good science, subsequent evaluation of published models and the supporting computer code is necessary for their consideration for use in HHRAs. Standardized model evaluation criteria and a thorough and efficient review process can reduce the number of review and revision iterations and hence the time needed to prepare a model for application. Efficient and consistent review also allows for rapid identification of needed model modifications to address HHRA-specific issues. This manuscript reports on the workshop where a process and criteria that were created for PBPK model review were discussed along with other issues related to model review and application in HHRA. Other issues include (1) model code availability, portability, and validity; (2) probabilistic (e.g., population-based) PBPK models and critical choices in parameter values to fully characterize population variability; and (3) approaches to integrating PBPK model outputs with other HHRA tools, including benchmark dose modeling. Two specific case study examples are provided to illustrate challenges that were encountered during the review and application process. By considering the frequent challenges encountered in the review and application of PBPK models during the model development phase, scientists may be better able to prepare their models for use in HHRAs.

Arsenic migration to deep groundwater in Bangladesh influenced by adsorption and water demand.

Drinking shallow groundwater with naturally elevated concentrations of arsenic is causing widespread disease in many parts of South and Southeast Asia. In the Bengal Basin, growing reliance on deep (>150 m) groundwater has lowered exposure. In the most affected districts of Bangladesh, shallow groundwater concentrations average 100 to 370 µg L(-1), while deep groundwater is typically < 10 µg L(-1). Groundwater flow simulations have suggested that, even when deep pumping is restricted to domestic use, deep groundwater in some areas of the Bengal Basin is at risk of contamination. However, these simulations have neglected the impedance of As migration by adsorption to aquifer sediments. Here we quantify for the first time As sorption on deeper sediments in situ by replicating the intrusion of shallow groundwater through injection of 1,000 L of deep groundwater modified with 200 µg L(-1) of As into a deeper aquifer. Arsenic concentrations in the injected water were reduced by 70% due to adsorption within a single day. Basin-scale modelling indicates that while As adsorption extends the sustainable use of deep groundwater, some areas remain vulnerable; these areas can be prioritized for management and monitoring.

Degradation rates of CFC-11, CFC-12 and CFC-113 in anoxic shallow aquifers of Araihazar, Bangladesh.

Chlorofluorocarbons CFC-11 (CCl(3)F), CFC-12 (CCl(2)F(2)), and CFC-113 (CCl(2)F-CClF(2)) are used in hydrology as transient tracers under the assumption of conservative behavior in the unsaturated and saturated soil zones. However, laboratory and field studies have shown that these compounds are not stable under anaerobic conditions. To determine the degradation rates of CFCs in a tropical environment, atmospheric air, unsaturated zone soil gas, and anoxic groundwater samples were collected in Araihazar upazila, Bangladesh. Observed CFC concentrations in both soil gas and groundwater were significantly below those expected from atmospheric levels. The CFC deficits in the unsaturated zone can be explained by gas exchange with groundwater undersaturated in CFCs. The CFC deficits observed in (3)H/(3)He dated groundwater were used to estimate degradation rates in the saturated zone. The results show that CFCs are degraded to the point where practically no (<5%) CFC-11, CFC-12, or CFC-113 remains in groundwater with (3)H/(3)He ages above 10 yr. In groundwater sampled at our site CFC-11 and CFC-12 appear to degrade at similar rates with estimated degradation rates ranging from approximately 0.25 yr(-1) to approximately 6 yr(-1). Degradation rates increased as a function of reducing conditions. This indicates that CFC dating of groundwater in regions of humid tropical climate has to be carried out with great caution.

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

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

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

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

A review of quantitative studies of benzene metabolism.

Benzene is a ubiquitous, highly flammable, colorless liquid that is a known hematotoxin, myelotoxin, and human leukemogen. Benzene-induced toxicity in animals is clearly mediated by its metabolism. The mechanisms of acute hemato- and myelotoxicity in humans are almost certainly the same as in animals, and there is compelling evidence that metabolism is requisite for the induction of leukemia in humans. A very large number of experimental investigations of benzene metabolism have been conducted with animals, both in vivo and in vitro. There have also been many investigations of benzene metabolism in humans and with human tissues, Although the blood or tissue concentrations of benzene metabolites in humans resulting from benzene exposure have never been measured. Further, a number of mathematical models of benzene metabolism and dosimetry have been developed. In this article, we consider results from both experimental and mathematical modeling research, with particular emphasis on the last decade, and discuss the factors that are likely to be most influential in the metabolism of benzene.

Quantitative exposure assessment: application of physiologically-based pharmacokinetic (PBPK) modeling of low-dose, long-term exposures of organic acid toxicant in the brain.

Our objective was to construct a physiologically-based pharmacokinetic (PBPK) model describing the kinetic behavior of 2,4-dichlorophenoxyacetic acid (2,4-D) on rats after long-term exposures to low doses. Our study demonstrated the model's ability to simulate uptake of 2,4-D in discrete areas of the rat brain. The model was derived from the generic PBPK model that was first developed for high-dose, single exposures of 2,4-D to rats or rabbits (Kim, C.S., Gargas, M.L., Andersen, M.E., 1994. Pharmacokinetic modeling of 2,4-dichlorophenoxyacetic acid (2,4-D) in rats and rabbits brain following single dose administration. Toxicol. Lett. 74, 189-201; Kim, C.S., Slikker, W., Jr., Binienda, Z., Gargas, M.L., Andersen, M.E., 1995. Development of a physiologically based pharmacokinetic (PBPK) model for 2,4-dichlorophenoxyacetic acid (2,4-D) dosimetry in discrete areas of the brain following a single intraperitoneal or intravenous dose. Neurotox. Teratol. 17, 111-120.), to which a subcutaneous (hypodermal) compartment was incorporated for low-dose, long-term infusion. It consisted of two body compartments, along with compartments for venous and arterial blood, cerebrospinal fluid, brain plasma and six brain regions. Uptake of the toxin was membrane-limited by the blood-brain barrier with clearance from the brain provided by cerebrospinal fluid 'sink' mechanisms. This model predicted profiles of 2,4-D levels in brain and blood over a 28-day period that compared well with concentrations measured in vivo with rats that had been given 2,4-D (1 or 10 mg/kg per day) with [14C]-2,4-D subcutaneously (s.c.) for 7, 14, or 28 days, respectively. This PBPK model should be an effective tool for evaluating the target tissue doses that may produce the neurotoxicity of organic acid toxicants after low-dose, long-term exposures to contaminated foods or the environment.

Physiologically based pharmacokinetic modeling of benzene metabolism in mice through extrapolation from in vitro to in vivo.

Benzene (C6H6) is a highly flammable, colorless liquid. Ubiquitous exposures result from its presence in gasoline vapors, cigarette smoke, and industrial processes. Benzene increases the incidence of leukemia in humans when they are exposed to high doses for extended periods; however, leukemia risks in humans at low exposures are uncertain. The exposure-dose-response relationship of benzene in humans is expected to be nonlinear because benzene undergoes a series of metabolic transformations, detoxifying and activating, in the liver, resulting in multiple metabolites that exert toxic effects on the bone marrow. We developed a physiologically based pharmacokinetic model for the uptake and elimination of benzene in mice to relate the concentration of inhaled and orally administered benzene to the tissue doses of benzene and its key metabolites, benzene oxide, phe nol, and hydroquinone. As many parameter values as possible were taken from the literature; in particular, metabolic parameters obtained from in vitro studies with mouse liver were used since comparable parameters are also available for humans. Parameters estimated by fitting the model to published data were first-order rate constants for pathways lacking in vitro data and the concentrations of microsomal and cytosolic protein, which effectively alter overall enzyme activity. The model was constrained by using the in vitro metabolic parameters (maximum velocities, first-order rate constants, and saturation parameters), and data from multiple laboratories and experiments were used. Despite these constraints and sources of variability, the model simulations matched the data reasonably well in most cases, showing that in vitro metabolic constants can be successfully extrapolated to predict in vivo data for benzene metabolism and dosimetry. Therefore in vitro metabolic constants for humans can subsequently be extrapolated to predict the dosimetry of benzene and its metabolites in humans. This will allow us to better estimate the risks of adverse effects from low-level benzene exposures.

A model of gonadotropin regulation during the menstrual cycle in women: qualitative features.

Increasing concerns that environmental contaminants may disrupt the endocrine system require development of mathematical tools to predict the potential for such compounds to significantly alter human endocrine function. The endocrine system is largely self-regulating, compensating for moderate changes in dietary phytoestrogens (e.g., in soy products) and normal variations in physiology. However, severe changes in dietary or oral exposures or in health status (e.g., anorexia), can completely disrupt the menstrual cycle in women. Thus, risk assessment tools should account for normal regulation and its limits. We present a mathematical model for the synthesis and release of luteinizing hormone (LH) and follicle-stimulating hormone (FSH) in women as a function of estrogen, progesterone, and inhibin blood levels. The model reproduces the time courses of LH and FSH during the menstrual cycle and correctly predicts observed effects of administered estrogen and progesterone on LH and FSH during clinical studies. The model should be useful for predicting effects of hormonally active substances, both in the pharmaceutical sciences and in toxicology and risk assessment.

Haber's rule: a special case in a family of curves relating concentration and duration of exposure to a fixed level of response for a given endpoint.

The concept that the product of the concentration (C) of a substance and the length of time (t) it is administered produces a fixed level of effect for a given endpoint has been ascribed to Fritz Haber, who was a German scientist in the early 1900s. He contended that the acute lethality of war gases could be assessed by the amount of the gas in a cubic meter of air (i.e. the concentration) multiplied by the time in min that the animal had to breathe the air before death ensued (i.e. C x t=k). While Haber recognized that C x t=k was applicable only under certain conditions, many toxicologists have used his rule to analyze experimental data whether or not their chemicals, biological endpoints, and exposure scenarios were suitable candidates for the rule. The fact that the relationship between C and t is linear on a log-log scale and could easily be solved by hand, led to early acceptance among toxicologists, particularly in the field of entomology. In 1940, a statistician named Bliss provided an elegant treatment on the relationships among exposure time, concentration, and the toxicity of insecticides. He proposed solutions for when the log-log plot of C and t was composed of two or more rectilinear segments, for when the log-log plot was curvilinear, and for when the slope of the dosage-mortality curve was a function of C. Despite the fact that Haber's rule can underestimate or overestimate effects (and consequently risks), it has been used in various settings by regulatory bodies. Examples are presented from the literature of data sets that follow Haber's rule as well as those that do not. Haber's rule is put into perspective by showing that it is simply a special case in a family of power law curves relating concentration and duration of exposure to a fixed level of response for a given endpoint. Also shown is how this power law family can be used to examine the three-dimensional surface relating C, t, and varying levels of response. The time has come to move beyond the limited view of C and t relationships inferred by Haber's rule to the use of the broader family of curves of which this rule is a special case.

Inhibition of cytochrome P450 2E1 decreases, but does not eliminate, genotoxicity mediated by 1,3-butadiene.

1,3-Butadiene (BD), a rodent carcinogen, is metabolized to mutagenic and DNA-reactive epoxides. In vitro data suggest that this oxidation is mediated by cytochrome P450 2E1 (CYP2E1). In this study, we tested the hypothesis that oxidation of BD by CYP2E1 is required for genotoxicity to occur. Inhalation exposures were conducted with B6C3F1 mice using a closed-chamber technique, and the maximum rate of butadiene oxidation was estimated. The total amount of butadiene metabolized was then correlated with the frequency of micronuclei (MN). Three treatment groups were used: (1) mice with no pretreatment; (2) mice pretreated with 1,2-trans-dichloroethylene (DCE), a specific CYP2E1 inhibitor; and (3) mice pretreated with 1-aminobenzotriazole (ABT), an irreversible inhibitor of cytochromes P450. Mice in all 3 groups were exposed to an initial BD concentration of 1100 ppm, and the decline in concentration of BD in the inhalation chamber with time, due to uptake and metabolism of BD, was monitored using gas chromatography. A physiologically based pharmacokinetic model was used to analyze the gas uptake data, estimate V(max) for BD oxidation, and compute the total amount of BD metabolized. Model simulations of the gas uptake data predicted the maximum rate of BD oxidation would be reduced by 60% and 100% for the DCE- and ABT-pretreated groups, respectively. Bone marrow was harvested 24 h after the onset of the inhalation exposure and analyzed for frequency of micronuclei in polychromatic erythrocytes (MN-PCE). The frequency of MN-PCE per 1000 PCE in mice exposed to BD was 28.2 +/- 3.1, 19.8 +/- 2.5, and 12.3 +/- 1.9, for the mice with no pretreatment, DCE-pretreated mice and ABT-pretreated mice, respectively. Although inhibition of CYP2E1 decreased BD-mediated genotoxicity, it did not completely eliminate genotoxic effects. These data suggest that other P450 isoforms may contribute significantly to the metabolic activation of BD and resultant genotoxicity.

Relative roles of convection and chemical reaction for the disposition of formaldehyde and ozone in nasal mucus.

The mucociliary apparatus is an important respiratory-tract defense system that may provide significant protection of the underlying epithelium from gases and vapors. Limiting-case calculations were performed to determine the significance of convective mucus transport and chemical reaction for formaldehyde (HCHO) and ozone (O(3)) in rat nasal respiratory epithelial mucus. Less than 4.6% of absorbed HCHO can be bound to amino groups (serum albumin) after 20 min of exposure. Thus, at the slowest measured mucus flow rates in rats, approximately 1 mm/min, a fluid element of mucus could travel more than 2 cm before binding 5% of absorbed HCHO, by which time the element would probably leave the nose (the site of toxic responses). In other calculations, HCHO removed by chemical reaction from a volume of mucus exposed for longer times was determined to be less than 0.54% of that removed by mucus flow (convection). Given the solubility of HCHO in mucus (water) and estimates of total mucus flow, however, as much as 22-42% of inhaled HCHO may be removed by total mucus flow. Alternately, O(3) dissolved in mucus would react completely with unsaturated fatty acids in 8.3 x 10(-4) s, in which time the mucus could flow no more than approximately 0.42 microm at the maximum reported flow rate of 30 mm/min. Even if a volume of mucus is flushed by net flow in 1 s, the amount of O(3) removed by flow would only be 0.12% of that removed by chemical reaction. Finally, based on the solubility of ozone, less than 8.0 x 10(-5)% of inhaled material could be removed from the nose by mucus flow. These results indicate which mucociliary processes are significant in site-specific dosimetry modeling.

Use of a mathematical model of rodent in vitro benzene metabolism to predict human in vitro metabolism data.

Benzene, a ubiquitous environmental pollutant, is known to cause leukemia and aplastic anemia in humans and hematotoxicity and myelotoxicity in rodents. Toxicity is thought to be exerted through oxidative metabolites formed in the liver, primarily via pathways mediated by cytochrome P450 2E1 (CYP2E1). Phenol, hydroquinone and trans-trans-muconaldehyde have all been hypothesized to be involved in benzene-induced toxicity. Recent reports indicate that benzene oxide is produced in vitro and in vivo and may be sufficiently stable to reach the bone marrow. Our goal was to improve existing mathematical models of microsomal benzene metabolism by including time course data for benzene oxide, by obtaining better parameter estimates and by determining if enzymes other than CYP2E1 are involved. Microsomes from male B6C3F1 mice and F344 rats were incubated with [(14)C]benzene (14 microM), [(14)C]phenol (303 microM) and [(14)C]hydroquinone (8 microM). Benzene and phenol were also incubated with mouse microsomes in the presence of trans-dichloroethylene, a CYP2E1 inhibitor, and benzene was incubated with trichloropropene oxide, an epoxide hydrolase inhibitor. These experiments did not indicate significant contributions of enzymes other than CYP2E1. Mathematical model parameters were fitted to rodent data and the model was validated by predicting human data. Model simulations predicted the qualitative behavior of three human time course data sets and explained up to 81% of the total variation in data from incubations of benzene for 16 min with microsomes from nine human individuals. While model predictions did deviate systematically from the data for benzene oxide and trihydroxybenzene, overall model performance in predicting the human data was good. The model should be useful in quantifying human risk due to benzene exposure and explicitly accounts for interindividual variation in CYP2E1 activity.