Reconstructing population exposures to environmental chemicals from biomarkers: challenges and opportunities.

A conceptual/computational framework for exposure reconstruction from biomarker data combined with auxiliary exposure-related data is presented, evaluated with example applications, and examined in the context of future needs and opportunities. This framework employs physiologically based toxicokinetic (PBTK) modeling in conjunction with numerical "inversion" techniques. To quantify the value of different types of exposure data "accompanying" biomarker data, a study was conducted focusing on reconstructing exposures to chlorpyrifos, from measurements of its metabolite levels in urine. The study employed biomarker data as well as supporting exposure-related information from the National Human Exposure Assessment Survey (NHEXAS), Maryland, while the MENTOR-3P system (Modeling ENvironment for TOtal Risk with Physiologically based Pharmacokinetic modeling for Populations) was used for PBTK modeling. Recently proposed, simple numerical reconstruction methods were applied in this study, in conjunction with PBTK models. Two types of reconstructions were studied using (a) just the available biomarker and supporting exposure data and (b) synthetic data developed via augmenting available observations. Reconstruction using only available data resulted in a wide range of variation in estimated exposures. Reconstruction using synthetic data facilitated evaluation of numerical inversion methods and characterization of the value of additional information, such as study-specific data that can be collected in conjunction with the biomarker data. Although the NHEXAS data set provides a significant amount of supporting exposure-related information, especially when compared to national studies such as the National Health and Nutrition Examination Survey (NHANES), this information is still not adequate for detailed reconstruction of exposures under several conditions, as demonstrated here. The analysis presented here provides a starting point for introducing improved designs for future biomonitoring studies, from the perspective of exposure reconstruction; identifies specific limitations in existing exposure reconstruction methods that can be applied to population biomarker data; and suggests potential approaches for addressing exposure reconstruction from such data.

Parameters for Carbamate Pesticide QSAR and PBPK/PD Models for Human Risk Assessment.

Our interest in providing parameters for the development of quantitative structure physiologically based pharmacokinetic/pharmacodynamic (QSPBPK/PD) models for assessing health risks to carbamates (USEPA 2005) comes from earlier work with organophosphorus (OP) insecticides (Knaak et al. 2004). Parameters specific to each carbamate are needed in the construction of PBPK/PD models along with their metabolic pathways. Parameters may be obtained by (1) development of QSAR models, (2) collecting pharmacokinetic data, and (3) determining pharmacokinetic parameters by fitting to experimental data. The biological parameters are given in Table 1 (Blancato et al. 2000). Table 1 Biological Parameters Required for Carbamate Pesticide Physiologically Based Pharmacokinetic/Pharmacodynamic (PBPK/PD) Models.(a).

A physiologically based pharmacokinetic/pharmacodynamic model for carbofuran in Sprague-Dawley rats using the exposure-related dose estimating model.

Carbofuran (2,3-dihydro-2,2-dimethyl-7-benzofuranyl-N-methylcarbamate), a broad spectrum N-methyl carbamate insecticide, and its metabolite, 3-hydroxycarbofuran, exert their toxicity by reversibly inhibiting acetylcholinesterase (AChE). To characterize AChE inhibition from carbofuran exposure, a physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model was developed in the Exposure-Related Dose Estimating Model (ERDEM) platform for the Sprague-Dawley (SD) rat. Experimental estimates of physiological, biochemical, and physicochemical model parameters were obtained or based on data from the open literature. The PBPK/PD model structure included carbofuran metabolism in the liver to 16 known metabolites, enterohepatic circulation of glucuronic acid conjugates, and excretion in urine and feces. Bolus doses by ingestion of 50 microg/kg and 0.5 mg/kg carbofuran were simulated for the blood and brain AChE activity. The carbofuran ERDEM simulated a half-life of 5.2 h for urinary clearance, and the experimental AChE activity data were reproduced for the blood and brain. Thirty model parameters were found influential to the model outputs and were chosen for perturbation in Monte Carlo simulations to evaluate the impact of their variability on the model predictions. Results of the simulation runs indicated that the minimum AChE activity in the blood ranged from 29.3 to 79.0% (as 5th and 95th percentiles) of the control level with a mean of 55.9% (standard deviation = 15.1%) compared to an experimental value of 63%. The constructed PBPK/PD model for carbofuran in the SD rat provides a foundation for extrapolating to a human model that can be used for future risk assessment.

In vitro metabolism of the fungicide and environmental contaminant trans-bromuconazole and implications for risk assessment.

trans-Bromuconazole is a chiral chemical representative of a class of triazole derivatives known to inhibit specific fungal cytochrome P-450 (CYP) reactions. Kinetic measurements and delineation of metabolic pathways for triazole chemicals within in vitro hepatic microsomes are needed for accurate risk assessment and predictive in vivo physiological modeling. The studies described here were conducted with rat liver microsomes to determine Michaelis-Menten saturation kinetic parameters (Vmax and KM) for trans-bromuconazole using both substrate depletion and product formation reaction velocities. Kinetic parameters determined for trans-bromuconazole depletion at varying protein levels incubated at physiological temperature 37 degrees C resulted in a KM value of 1.69 microM and a Vmax value of 1398 pmol/min/mg protein. The concomitant linear formation of two metabolites identified using liquid chromatography/time-of-flight mass spectrometry (LC/MS-TOF) and LC-MS/MS indicated hydroxylation of the trans-bromuconazole dichlorophenyl ring moiety. KM values determined for the hydroxylated metabolites were 0.87 and 1.03 microM, with Vmax values of 449 and 694 pmol/min/mg protein, respectively. Chemical inhibition assays and studies conducted with individual purified human recombinant enzymes indicated the CYP3A subfamily was primarily responsible for biotransformation of the parent substrate. Additionally, trans-bromuconazole was found to undergo stereoselective metabolism as evidenced by a change in the enantiomeric ratio (trans-/trans+) with respect to time.