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.

Reduction of deepwater formation in the greenland sea during the 1980s: evidence from tracer data.

Hydrographic observations and measurements of the concentrations of chlorofluorocarbons (CFCs) have suggested that the formation of Greenland Sea Deep Water (GSDW) slowed down considerably during the 1980s. Such a decrease is related to weakened convection in the Greenland Sea and thus could have significant impact on the properties of the waters flowing over the Scotland-Iceland-Greenland ridge system into the deep Atlantic. Study of the variability of GSDW formation is relevant for understanding the impact of the circulation in the European Polar seas on regional and global deep water characteristics. New long-term multitracer observations from the Greenland Sea show that GSDW formation indeed was greatly reduced during the 1980s. A box model of deepwater formation and exchange in the European Polar seas tuned by the tracer data indicates that the reduction rate of GSDW formation was about 80 percent and that the start date of the reduction was between 1978 and 1982.

Paleotemperatures in the southwestern United States derived from noble gases in ground water.

A paleotemperature record based on measurements of atmospheric noble gases dissolved in ground water of the Carrizo aquifer (Texas) shows that the annual mean temperature in the southwestern United States during the last glacial maximum was about 5 degrees C lower than the present-day value. In combination with evidence for fluctuations in mountain snow lines, this cooling indicates that the glacial lapse rate was approximately the same as it is today. In contrast, measurements on deep-sea sediments indicate that surface temperatures in the ocean basins adjacent to our study area decreased by only about 2 degrees C. This difference between continental and oceanic records poses questions concerning our current understanding of paleoclimate and climate-controlling processes.

Cooling of Tropical Brazil (5{degrees}C) During the Last Glacial Maximum.

A 30,000-year paleotemperature record derived from noble gases dissolved in carbon-14-dated ground water indicates that the climate in lowland Brazil (Piaui Province, 7 degrees S, 41.5 degrees W; altitude, 400 meters) was 5.4 degrees +/- 0.6 degrees C cooler during the last glacial maximum than today. This result suggests a rather uniform cooling of the Americas between 40 degrees S and 40 degrees N. A 5.4 degrees C cooling of tropical South America is consistent with pollen records, snow line reconstructions, and strontium/calcium ratios and delta(18)O coral records but is inconsistent with the sea-surface temperature reconstruction of CLIMAP (Climate: Long-Range Investigation, Mapping and Prediction). On the basis of these results, it appears that the tropical Americas are characterized by a temperature sensitivity comparable to that found in higher latitudes.

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.

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.

Critical issues in benzene toxicity and metabolism: the effect of interactions with other organic chemicals on risk assessment.

Benzene, an important industrial solvent, is also present in unleaded gasoline and cigarette smoke. The hematotoxic effects of benzene are well documented and include aplastic anemia and pancytopenia. Some individuals exposed repeatedly to cytotoxic concentrations of benzene develop acute myeloblastic anemia. It has been hypothesized that metabolism of benzene is required for its toxicity, although administration of no single benzene metabolite duplicates the toxicity of benzene. Several investigators have demonstrated that a combination of metabolites (hydroquinone and phenol, for example) is necessary to duplicate the hematotoxic effect of benzene. Enzymes implicated in the metabolic activation of benzene and its metabolites include the cytochrome P450 monooxygenases and myeloperoxidase. Since benzene and its hydroxylated metabolites (phenol, hydroquinone, and catechol) are substrates for the same cytochrome P450 enzymes, competitive interactions among the metabolites are possible. In vivo data on metabolite formation by mice exposed to various benzene concentrations are consistent with competitive inhibition of phenol oxidation by benzene. Other organic molecules that are substrates for cytochrome P450 can inhibit the metabolism of benzene. For example, toluene has been shown to inhibit the oxidation of benzene in a noncompetitive manner. Enzyme inducers, such as ethanol, can alter the target tissue dosimetry of benzene metabolites by inducing enzymes responsible for oxidation reactions involved in benzene metabolism. The dosimetry of benzene and its metabolites in the target tissue, bone marrow, depends on the balance of activation processes, such as enzymatic oxidation, and deactivation processes, like conjugation and excretion.(ABSTRACT TRUNCATED AT 250 WORDS)

Mechanistic considerations in benzene physiological model development.

Benzene, an important industrial solvent, is also present in unleaded gasoline and cigarette smoke. The hematotoxic effects of benzene in humans are well documented and include aplastic anemia, pancytopenia, and acute myelogenous leukemia. However, the risks of leukemia at low exposure concentrations have not been established. A combination of metabolites (hydroquinone and phenol, for example) may be necessary to duplicate the hematotoxic effect of benzene, perhaps due in part to the synergistic effect of phenol on myeloperoxidase-mediated oxidation of hydroquinone to the reactive metabolite benzoquinone. Because benzene and its hydroxylated metabolites (phenol, hydroquinone, and catechol) are substrates for the same cytochrome P450 enzymes, competitive interactions among the metabolites are possible. In vivo data on metabolite formation by mice exposed to various benzene concentrations are consistent with competitive inhibition of phenol oxidation by benzene. In vitro studies of the metabolic oxidation of benzene, phenol, and hydroquinone are consistent with the mechanism of competitive interaction among the metabolites. The dosimetry of benzene and its metabolites in the target tissue, bone marrow, depends on the balance of activation processes such as enzymatic oxidation and deactivation processes such as conjugation and excretion. Phenol, the primary benzene metabolite, can undergo both oxidation and conjugation. Thus the potential exists for competition among various enzymes for phenol. Zonal localization of phase I and phase II enzymes in various regions of the liver acinus also impacts this competition. Biologically based dosimetry models that incorporate the important determinants of benzene flux, including interactions with other chemicals, will enable prediction of target tissue doses of benzene and metabolites at low exposure concentrations relevant for humans.

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.

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.

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.

Benzene: a case study in parent chemical and metabolite interactions.

Benzene, an important industrial solvent, is also present in unleaded gasoline and cigarette smoke. The hematotoxic effects of benzene in humans are well documented and include aplastic anemia and pancytopenia, and acute myelogenous leukemia. A combination of metabolites (hydroquinone and phenol for example) is apparently necessary to duplicate the hematotoxic effect of benzene, perhaps due in part to the synergistic effect of phenol on myeloperoxidase-mediated oxidation of hydroquinone to the reactive metabolite benzoquinone. Since benzene and its hydroxylated metabolites (phenol, hydroquinone and catechol) are substrates for the same cytochrome P450 enzymes, competitive interactions among the metabolites are possible. In vivo data on metabolite formation by mice exposed to various benzene concentrations are consistent with competitive inhibition of phenol oxidation by benzene. In vitro studies of the metabolic oxidation of benzene, phenol and hydroquinone are consistent with the mechanism of competitive interaction among the metabolites. The dosimetry of benzene and its metabolites in the target tissue, bone marrow, depends on the balance of activation processes such as enzymatic oxidation and deactivation processes such as conjugation and excretion. Phenol, the primary benzene metabolite, can undergo both oxidation and conjugation. Thus, the potential exists for competition among various enzymes for phenol. However, zonal localization of Phase I and Phase II enzymes in various regions of the liver acinus regulates this competition. Biologically-based dosimetry models that incorporate the important determinants of benzene flux, including interactions with other chemicals, will enable prediction of target tissue doses of benzene and metabolites at low exposure concentrations relevant for humans.

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.

Physiologically based pharmacokinetic modeling of blood and tissue epoxide measurements for butadiene.

In vitro and in vivo butadiene (BD) metabolism data from laboratory animals were integrated into a rodent physiologically based pharmacokinetic (PBPK) model with flow- and diffusion-limited compartments. The resulting model describes experimental data from multiple sources under scenarios such as closed chamber inhalation and nose-only flow-through inhalation exposures. Incorporation of diurnal glutathione (GSH) variation allows accurate simulation of GSH changes observed in air control nose-only exposures and BD exposures. An isolated tissue model based on rate parameters determined in vitro predicts the decrease in epoxide concentrations in intact animals during the time lag between exsanguination and tissue removal for tissues capable of epoxide biotransformation, providing a better indication of in vivo dosimetry. Further refinements of the model are required relative to model predictions of an important BD metabolite, diepoxybutane.

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.

An integrated modeling-experimental strategy for the analysis of metabolic pathways.

An inherent problem in studying the behavior of a metabolic pathway is the impossibility of developing a complete, detailed model that includes all the cellular processes that have an impact on the set of fluxes in such a pathway. Lacking this, one requires some means of modeling the interactions between a metabolic pathway and other cellular processes for the purpose of analyzing pathway characteristics within the cell (e.g., determining sensitivity coefficients for various steps in the pathway) with a minimal amount of time and effort. A general framework is developed for studying these issues in a rigorous manner. Using this framework, detailed knowledge about a metabolic pathway (i.e., a set of rate expressions for steps in the pathway) can be combined with the results from a relatively simple set of experiments in order to obtain estimates for the sensitivity of the pathway to enzyme activities, inhibition constants, and other parameters that determine the pathway's behavior, while accounting for the pathway's interaction with the rest of the cellular metabolism. A model system representing amino acid production is used to illustrate the problem and to provide results based on computational experiments. The modeling strategy described here should be useful in genetic design to improve pathway fluxes and metabolic network selectivity.