Development and application of a human PBPK model for bromodichloromethane to investigate the impacts of multi-route exposure.

As a result of its presence in water as a volatile disinfection byproduct, bromodichloromethane (BDCM), which is mutagenic, poses a potential health risk from exposure via oral, dermal and inhalation routes. We developed a refined human physiologically based pharmacokinetic (PBPK) model for BDCM (including new chemical-specific human parameters) to evaluate the impact of BDCM exposure during showering and bathing on important measures of internal dose compared with oral exposure. The refined model adequately predicted data from the published literature for oral, dermal and bathing/showering exposures. A liter equivalency approach (L-eq) was used to estimate BDCM concentration in a liter of water consumed by the oral route that would be required to produce the same internal dose of BDCM resulting from a 20-min bath or a 10-min shower in water containing 10 µg l(-1) BDCM. The oral liter equivalent concentrations for the bathing scenario were 605, 803 and 5 µg l(-1) BDCM for maximum venous blood concentration (Cmax), the area under the curve (AUCv) and the amount metabolized in the liver per hour (MBDCM), respectively. For a 10-min showering exposure, the oral L-eq concentrations were 282, 312 and 2.1 µg l(-1) for Cmax, AUC and MBDCM, respectively. These results demonstrate large contributions of dermal and inhalation exposure routes to the internal dose of parent chemical reaching the systemic circulation, which could be transformed to mutagenic metabolites in extrahepatic target tissues. Thus, consideration of the contribution of multiple routes of exposure when evaluating risks from water-borne BDCM is needed, and this refined human model will facilitate improved assessment of internal doses from real-world exposures. Published 2015. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

The impact of variation in scaling factors on the estimation of internal dose metrics: a case study using bromodichloromethane (BDCM).

A rate for hepatic metabolism (Vmax) determined in vitro must be scaled for in vivo use in a physiologically based pharmacokinetic (PBPK) model. This requires the use of scaling factors such as mg of microsomal protein per gram of liver (MPPGL) and liver mass (FVL). Variation in MPPGL and FVL impacts variation in Vmax, and hence PBPK model-derived estimates of internal dose used in dose response analysis. The impacts of adult human variation in MPPGL and FVL on estimates of internal dose were assessed using a human PBPK model for bromodichloromethane (BDCM), a water disinfection byproduct, for multiple internal dose metrics for two exposure scenarios (single 0.25 liter drink of water or 10 min shower) under plausible (5 µg/L) and high level (20 µg/L) water concentrations. For both concentrations, all internal dose metrics were changed less than 5% for the showering scenario (combined inhalation and dermal exposure). In contrast, a 27-fold variation in area under the curve (AUC) for BDCM in venous blood was observed at both oral exposure concentrations, whereas total amount of BDCM metabolized in liver was relatively unchanged. This analysis demonstrates that variability in the scaling factors used for in vitro to in vivo extrapolation (IVIVE) for metabolic rate parameters can have a significant route-dependent impact on estimates of internal dose under environmentally relevant exposure scenarios. This indicates the need to evaluate both uncertainty and variability for scaling factors used for IVIVE.

Identifiability of PBPK models with applications to dimethylarsinic acid exposure.

Any statistical model should be identifiable in order for estimates and tests using it to be meaningful. We consider statistical analysis of physiologically-based pharmacokinetic (PBPK) models in which parameters cannot be estimated precisely from available data, and discuss different types of identifiability that occur in PBPK models and give reasons why they occur. We particularly focus on how the mathematical structure of a PBPK model and lack of appropriate data can lead to statistical models in which it is impossible to estimate at least some parameters precisely. Methods are reviewed which can determine whether a purely linear PBPK model is globally identifiable. We propose a theorem which determines when identifiability at a set of finite and specific values of the mathematical PBPK model (global discete identifiability) implies identifiability of the statistical model. However, we are unable to establish conditions that imply global discrete identifiability, and conclude that the only safe approach to analysis of PBPK models involves Bayesian analysis with truncated priors. Finally, computational issues regarding posterior simulations of PBPK models are discussed. The methodology is very general and can be applied to numerous PBPK models which can be expressed as linear time-invariant systems. A real data set of a PBPK model for exposure to dimethyl arsinic acid (DMA(V)) is presented to illustrate the proposed methodology.

Genomic comparisons between hepatocarcinogenic and non-hepatocarcinogenic organophosphate insecticides in the mouse liver.

Short-term biomarkers of toxicity have an increasingly important role in the screening and prioritization of new chemicals. In this study, we examined early indicators of liver toxicity for three reference organophosphate (OP) chemicals, which are among the most widely used insecticides in the world. The OP methidathion was previously shown to increase the incidence of liver toxicity, including hepatocellular tumors, in male mice. To provide insights into the adverse outcome pathway (AOP) that underlies these tumors, effects of methidathion in the male mouse liver were examined after 7 and 28 day exposures and compared to those of two other OPs that either do not increase (fenthion) or possibly suppress liver cancer (parathion) in mice. None of the chemicals caused increases in liver weight/body weight or histopathological changes in the liver. Parathion decreased liver cell proliferation after 7 and 28 days while the other chemicals had no effects. There was no evidence for hepatotoxicity in any of the treatment groups. Full-genome microarray analysis of the livers from the 7 and 28 day treatments demonstrated that methidathion and fenthion regulated a large number of overlapping genes, while parathion regulated a unique set of genes. Examination of cytochrome P450 enzyme activities and use of predictive gene expression biomarkers found no consistent evidence for activation of AhR, CAR, PXR, or PPARα. Parathion suppressed the male-specific gene expression pattern through STAT5b, similar to genetic and dietary conditions that decrease liver tumor incidence in mice. Overall, these findings indicate that methidathion causes liver cancer by a mechanism that does not involve common mechanisms of liver cancer induction.

Arsenic toxicokinetic modeling and risk analysis: Progress, needs and applications.

Arsenic (As) poses unique challenges in PBTK model development and risk analysis applications. Arsenic metabolism is complex, adequate information to attribute specific metabolites to particular adverse effects in humans is sparse, and measurement of relevant metabolites in biological media can be difficult. Multiple As PBTK models have been published and used or adapted for use in various exposure and risk analysis applications. These applications illustrate the broad utility of PBTK models for exposure and dose-response analysis, particularly for arsenic where multi-pathway, multi-route exposures and multiple toxic effects are of concern. Arsenic PBTK models have been used together with exposure reconstruction and dose-response functions to estimate risk of specific adverse health effects due to drinking water exposure and consumption of specific foodstuffs (e.g. rice, seafood), as well as to derive safe exposure levels and develop consumption advisories. Future refinements to arsenic PBTK models can enhance the confidence in such analyses. Improved estimates for methylation biotransformation parameters based on in vitro to in vivo extrapolation (IVIVE) methods and estimation of interindividual variability in key model parameters for specific toxicologically relevant metabolites are two important areas for consideration.

Aging and susceptibility to toluene in rats: a pharmacokinetic, biomarker, and physiological approach.

Aging adults are a growing segment of the U.S. population and are likely to exhibit increased susceptibility to many environmental toxicants. However, there is little information on the susceptibility of the aged to toxicants. The toxicity of toluene has been well characterized in young adult rodents but there is little information in the aged. Three approaches were used: (1) pharmacokinetic (PK), (2) cardiac biomarkers, and (3) whole-animal physiology to assess whether aging increases susceptibility to toluene in the Brown Norway (BN) rat. Three life stages, young adult, middle aged, and aged (4, 12, and 24 mo, respectively), were administered toluene orally at doses of 0, 0.3, 0.65, or 1 g/kg and subjected to the following: terminated at 45 min or 4 h post dosing, and blood and brain toluene concentration were measured; terminated at 4 h post dosing, and biomarkers of cardiac function were measured; or monitor heart rate (HR), core temperature (Tc), and motor activity (MA) by radiotelemetry before and after dosing. Brain toluene concentration was significantly elevated in aged rats at 4 h after dosing with either 0.3 or 1 g/kg. Blood toluene concentrations were unaffected by age. There were various interactions between aging and toluene-induced effects on cardiac biomarkers. Most notably, toluene exposure led to reductions in mRNA markers for oxidative stress in aged but not younger animals. Toluene also produced a reduction in cardiac endothelin-1 in aged rats. Higher doses of toluene led to tachycardia, hypothermia, and a transient elevation in MA. Aged rats were less sensitive to the tachycardic effects of toluene but showed a prolonged hypothermic response. Elevated brain levels of toluene in aged rats may be attributed to their suppressed cardiovascular and respiratory responses. The expression of several cardiac biochemical markers of toluene exposure in the aged may also reflect differential susceptibility to this toxicant.

Comparison of in vivo derived and scaled in vitro metabolic rate constants for several volatile organic compounds (VOCs).

Metabolic rate parameters estimation using in vitro data is necessary due to numbers of chemicals for which data are needed, trend towards minimizing laboratory animal use, and limited opportunity to collect data in human subjects. We evaluated how well metabolic rate parameters derived from in vitro data predict overall in vivo metabolism for a set of environmental chemicals for which well validated and established methods exist. We compared values of VmaxC derived from in vivo vapor uptake studies with estimates of VmaxC scaled up from in vitro hepatic microsomal metabolism studies for VOCs for which data were available in male F344 rats. For 6 of 7 VOCs, differences between the in vivo and scaled up in vitro VmaxC estimates were less than 2.6-fold. For bromodichloromethane (BDCM), the in vivo derived VmaxC was approximately 4.4-fold higher than the in vitro derived and scaled up VmaxC. The more rapid rate of BDCM metabolism estimated based in vivo studies suggests other factors such as extrahepatic metabolism, binding or other non-specific losses making a significant contribution to overall clearance. Systematic and reliable utilization of scaled up in vitro biotransformation rate parameters in PBPK models will require development of methods to predict cases in which extrahepatic metabolism and binding as well as other factors are likely to be significant contributors.

The Impact of Scaling Factor Variability on Risk-Relevant Pharmacokinetic Outcomes in Children: A Case Study Using Bromodichloromethane (BDCM).

Biotransformation rates extrapolated from in vitro data are used increasingly in human physiologically based pharmacokinetic (PBPK) models. This practice requires use of scaling factors, including microsomal content (mg of microsomal protein/g liver, MPPGL), enzyme specific content, and liver mass as a fraction of body weight (FVL). Previous analyses indicated that scaling factor variability impacts pharmacokinetic (PK) outcomes used in adult population dose-response studies. This analysis was extended to pediatric populations because large inter-individual differences in enzyme ontogeny likely would further contribute to scaling factor variability. An adult bromodichloromethane (BDCM) model (Kenyon, E. M., Eklund, C., Leavens, T. L., and Pegram, R. A. (2016a). Development and application of a human PBPK model for bromodichloromethane (BDCM) to investigate impacts of multi-route exposure. J. Appl. Toxicol. 36, 1095-1111) was re-parameterized for neonates, infants, and toddlers. Monte Carlo analysis was used to assess the impact of pediatric scaling factor variation on model-derived PK outcomes compared with adult findings. BDCM dose metrics were estimated following a single 0.05-liter drink of water or a 20-min bath, under typical (5 µg/l) and plausible higher (20 µg/l) BDCM concentrations. MPPGL, CYP2E1, and FVL values reflected the distribution of reported pediatric population values. The impact of scaling factor variability on PK outcome variation was different for each exposure scenario, but similar for each BDCM water concentration. The higher CYP2E1 expression variability during early childhood was reflected in greater variability in predicted PK outcomes in younger age groups, particularly for the oral exposure route. Sensitivity analysis confirmed the most influential parameter for this variability was CYP2E1, particularly in neonates. These findings demonstrate the importance of age-dependent scaling factor variation used for in vitro to in vivo extrapolation of biotransformation rates.

How can biologically-based modeling of arsenic kinetics and dynamics inform the risk assessment process? - A workshop review.

Quantitative biologically-based models describing key events in the continuum from arsenic exposure to the development of adverse health effects provide a framework to integrate information obtained across diverse research areas. For example, genetic polymorphisms in arsenic metabolizing enzymes can lead to differences in target tissue dosimetry for key metabolites causative in toxic and carcinogenic response. This type of variation can be quantitatively incorporated into pharmacokinetic (PK) models and used together with population-based modeling approaches to evaluate the impact of genetic variation in methylation capacity on dose of key metabolites to target tissue. The PK model is an essential bridge to the pharmacodynamic (PD) models. A particular benefit of PD modeling for arsenic is that alternative models can be constructed for multiple proposed modes of action for arsenicals. Genomics data will prove useful for identifying the key pathways involved in particular responses and aid in determining other types of data needed for quantitative modeling. These models, when linked with PK models, can be used to better understand and explain dose- and time-response behaviors. This in turn assists in prioritizing modes of action with respect to their risk assessment relevance and future research. This type of integrated modeling approach can form the basis for a highly informative mode-of-action directed risk assessment for inorganic arsenic (iAs). This paper will address both practical and theoretical aspects of integrating PK and PD data in a modeling framework, including practical barriers to its application.

Modeling the toxicokinetics of inhaled toluene in rats: influence of physical activity and feeding status.

Toluene is found in petroleum-based fuels and used as a solvent in consumer products and industrial applications. The critical effects following inhalation exposure involve the brain and nervous system in both humans and experimental animals, whether exposure duration is acute or chronic. The goals of this physiologically based pharmacokinetic (PBPK) model development effort were twofold: (1) to evaluate and explain the influence of feeding status and activity level on toluene pharmacokinetics utilizing our own data from toluene-exposed Long Evans (LE) rats, and (2) to evaluate the ability of the model to simulate data from the published literature and explain differing toluene kinetics. Compartments in the model were lung, slowly and rapidly perfused tissue groups, fat, liver, gut, and brain; tissue transport was blood-flow limited and metabolism occurred in the liver. Chemical-specific parameters and initial organ volumes and blood flow rates were obtained from the literature. Sensitivity analysis revealed that the single most influential parameter for our experimental conditions was alveolar ventilation; other moderately influential parameters (depending upon concentration) included cardiac output, rate of metabolism, and blood flow to fat. Based on both literature review and sensitivity analysis, other parameters (e.g., partition coefficients and metabolic rate parameters) were either well defined (multiple consistent experimental results with low variability) or relatively noninfluential (e.g. organ volumes). Rats that were weight-maintained compared to free-fed rats in our studies could be modeled with a single set of parameters because feeding status did not have a significant impact on toluene pharmacokinetics. Heart rate (HR) measurements in rats performing a lever-pressing task indicated that the HR increased in proportion to task intensity. For rats acclimated to eating in the lab during the day, both sedentary rats and rats performing the lever-pressing task required different alveolar ventilation rates to successfully predict the data. Model evaluation using data from diverse sources together with statistical evaluation of the resulting fits revealed that the model appropriately predicted blood and brain toluene concentrations with some minor exceptions. These results (1) emphasize the importance of experimental conditions and physiological status in explaining differing kinetic data, and (2) demonstrate the need to consider simulation conditions when estimating internal dose metrics for toxicity studies in which kinetic data were not collected.

Development of a human physiologically based pharmacokinetic (PBPK) model for inorganic arsenic and its mono- and di-methylated metabolites.

A physiologically-based pharmacokinetic (PBPK) model was developed to estimate levels of arsenic and its metabolites in human tissues and urine after oral exposure to arsenate (As(V)), arsenite (As(III)) or organoarsenical pesticides. The model consists of interconnected individual PBPK models for inorganic arsenic (As(V) and As(III)), monomethylarsenic acid (MMA(V)), and, dimethylarsenic acid (DMA(V)). Reduction of MMA(V) and DMA(V) to their respective trivalent forms also occurs in the lung, liver, and kidney including excretion in urine. Each submodel was constructed using flow limited compartments describing the mass balance of the chemicals in GI tract (lumen and tissue), lung, liver, kidney, muscle, skin, heart, and brain. The choice of tissues was based on physiochemical properties of the arsenicals (solubility), exposure routes, target tissues, and sites for metabolism. Metabolism of inorganic arsenic in liver was described as a series of reduction and oxidative methylation steps incorporating the inhibitory influence of metabolites on methylation. The inhibitory effects of As(III) on the methylation of MMA(III) to DMA, and MMA(III) on the methylation of As(III) to MMA were modeled as noncompetitive. To avoid the uncertainty inherent in estimation of many parameters from limited human data, a priori independent parameter estimates were derived using data from diverse experimental systems with priority given to data derived using human cells and tissues. This allowed the limited data for human excretion of arsenicals in urine to be used to estimate only parameters that were most sensitive to this type of data. Recently published urinary excretion data, not previously used in model development, are also used to evaluate model predictions.

Quantitative comparisons of the acute neurotoxicity of toluene in rats and humans.

The behavioral and neurophysiological effects of acute exposure to toluene are the most thoroughly explored of all the hydrocarbon solvents. Behavioral effects have been experimentally studied in humans and other species, for example, rats. The existence of both rat and human dosimetric data offers the opportunity to quantitatively compare the relative sensitivity to acute toluene exposure. The purpose of this study was to fit dose-effect curves to existing data and to estimate the dose-equivalence equation (DEE) between rats and humans. The DEE gives the doses that produce the same magnitude of effect in the two species. Doses were brain concentrations of toluene estimated from physiologically based pharmacokinetic models. Human experiments measuring toluene effects on choice reaction time (CRT) were meta-analyzed. Rat studies employed various dependent variables: amplitude of visual-evoked potentials (VEPs), signal detection (SIGDET) accuracy (ACCU) and reaction time (RT), and escape-avoidance (ES-AV) behaviors. Comparison of dose-effect functions showed that human and rat sensitivity was practically the same for those two task regimens that exerted the least control over the behaviors being measured (VEP in rats and CRT in humans) and the sensitivity was progressively lower for SIGDET RT, SIGDET ACCU, and ES-AV behaviors in rats. These results suggested that the sensitivity to impairment by toluene depends on the strength of control over the measured behavior rather than on the species being tested. This interpretation suggests that (1) sensitivity to toluene would be equivalent in humans and rats if both species performed behaviors that were controlled to the same extent, (2) the most sensitive tests of neurobehavioral effects would be those in which least control is exerted on the behavior being measured, and (3) effects of toluene in humans may be estimated using the DEEs from rat studies despite differences in the amount of control exerted by the experimental regimen or differences in the behaviors under investigation.

A dosimetric analysis of the acute behavioral effects of inhaled toluene in rats.

Knowledge of the appropriate metric of dose for a toxic chemical facilitates quantitative extrapolation of toxicity observed in the laboratory to the risk of adverse effects in the human population. Here, we utilize a physiologically based toxicokinetic (PBTK) model for toluene, a common volatile organic compound (VOC), to illustrate that its acute behavioral effects in rats can be quantitatively predicted on the basis of its concentration in the brain. Rats previously trained to perform a visual signal detection task for food reward performed the task while inhaling toluene (0, 1200, 1600, 2000, and 2400 ppm in different test sessions). Accuracy and speed of responding were both decreased by toluene; the magnitude of these effects increased with increasing concentration of the vapor and with increasing duration of exposure. Converting the exposure conditions to brain toluene concentration using the PBTK model yielded a family of overlapping curves for each end point, illustrating that the effects of toluene can be described quantitatively by its internal dose at the time of behavioral assessment. No other dose metric, including inhaled toluene concentration, duration of exposure, the area under the curve of either exposure (ppm h), or modeled brain toluene concentration (mg-h/kg), provided unambiguous predictions of effect. Thus, the acute behavioral effects of toluene (and of other VOCs with a similar mode of action) can be predicted for complex exposure scenarios by simulations that estimate the concentration of the VOC in the brain from the exposure scenario.

Acute toluene exposure and rat visual function in proportion to momentary brain concentration.

Acute exposure to toluene was assessed in two experiments to determine the relationship between brain toluene concentration and changes in neurophysiological function. The concentration of toluene in brain tissue at the time of assessment was estimated using a physiologically based pharmacokinetic model. Brain neurophysiological function was measured using pattern-elicited visual evoked potentials (VEP) recorded from electrodes located over visual cortex of adult male Long-Evans rats. In the first experiment, VEPs were recorded before and during exposure to control air or toluene at 1000 ppm for 4 h, 2000 ppm for 2 h, 3000 ppm for 1.3 h, or 4000 ppm for 1 h. In the second experiment, VEPs were recorded during and after exposure to clean air or 3000 or 4000 ppm toluene. In both experiments, the response amplitude of the major spectral component of the VEP (F2 at twice the stimulus rate in steady-state responses) was reduced by toluene. A logistic function was fit to baseline-adjusted F2 amplitudes from the first experiment that described a significant relationship between brain toluene concentration and VEP amplitude deficits. In the second experiment, 3000 ppm caused equivalent VEP deficits during or after exposure as a function of estimated brain concentration, but 4000 ppm showed a rapid partial adaptation to the acute effects of toluene after exposure. In general, however, the neurophysiological deficits caused by acute toluene exposure could be described by estimates of the momentary concentration of toluene in the brain at the time of VEP evaluation.

From the Cover: Genomic Effects of Androstenedione and Sex-Specific Liver Cancer Susceptibility in Mice.

Current strategies for predicting carcinogenic mode of action for nongenotoxic chemicals are based on identification of early key events in toxicity pathways. The goal of this study was to evaluate short-term key event indicators resulting from exposure to androstenedione (A4), an androgen receptor agonist and known liver carcinogen in mice. Liver cancer is more prevalent in men compared with women, but androgen-related pathways underlying this sex difference have not been clearly identified. Short-term hepatic effects of A4 were compared with reference agonists of the estrogen receptor (ethinyl estradiol, EE) and glucocorticoid receptor (prednisone, PRED). Male B6C3F1 mice were exposed for 7 or 28 days to A4, EE, or PRED. EE increased and PRED suppressed hepatocyte proliferation, while A4 had no detectable effects. In a microarray analysis, EE and PRED altered >3000 and >670 genes, respectively, in a dose-dependent manner, whereas A4 did not significantly alter any genes. Gene expression was subsequently examined in archival liver samples from male and female B6C3F1 mice exposed to A4 for 90 days. A4 altered more genes in females than males and did not alter expression of genes linked to activation of the mitogenic xenobiotic receptors AhR, CAR, and PPARα in either sex. A gene expression biomarker was used to show that in female mice, the high dose of A4 activated the growth hormone-regulated transcription factor STAT5b, which controls sexually dimorphic gene expression in the liver. These findings suggest that A4 induces subtle age-related effects on STAT5b signaling that may contribute to the higher risk of liver cancer in males compared with females.

A General Physiological and Toxicokinetic (GPAT) Model for Simulating Complex Toluene Exposure Scenarios in Humans.

Realistic simulation of environmental exposure scenarios requires dynamic methods in which exposures and human activities vary continuously as a function of time. Simulation of such complex scenarios is, with conventional physiologically based methods, a complex and programming-intensive task. The goal of the present effort was to simplify this task by combining a commercially available general whole-body human physiological model (QCP2004) with a slightly extended physiologically based toxicokinetic (PBTK) model from the literature. The QCP2004 model is a differential equation-based model similar to PBTK models except that normal organ function is simulated and the body organs are appropriately interlinked. Here QCP2004 provided estimates of physiological parameters required by the PBTK model. These were updated as the model was iteratively executed appropriate to the varying activity of the human subject. The combined general physiological model and the PBTK model was called a general physiological and toxicokinetic (GPAT) model. The GPAT model was tested and (within the constraints of available toluene exposure experiments in the literature) found to predict toluene blood concentrations, even in dynamic situations. A model of the structure used in the present work is capable of expansion as new knowledge is developed and greater detail is desired. Similarly, multiple toxicant PBTK models can be developed and incorporated for applications to mixtures risk assessment. Additionally, toxicant effects on organ systems can be achieved by altering organ function during a simulation as a function of the internal dose of toxicants. By cumulatively adding detail to the model as new physiological and chemical-specific information becomes available, the model can become a repository of knowledge for increasingly sophisticated risk-assessment applications.

Dose and Effect Thresholds for Early Key Events in a PPARα-Mediated Mode of Action.

Current strategies for predicting adverse health outcomes of environmental chemicals are centered on early key events in toxicity pathways. However, quantitative relationships between early molecular changes in a given pathway and later health effects are often poorly defined. The goal of this study was to evaluate short-term key event indicators using qualitative and quantitative methods in an established pathway of mouse liver tumorigenesis mediated by peroxisome proliferator-activated receptor alpha (PPARα). Male B6C3F1 mice were exposed for 7 days to di (2-ethylhexyl) phthalate (DEHP), di-n-octyl phthalate (DNOP), and n-butyl benzyl phthalate (BBP), which vary in PPARα activity and liver tumorigenicity. Each phthalate increased expression of select PPARα target genes at 7 days, while only DEHP significantly increased liver cell proliferation labeling index (LI). Transcriptional benchmark dose (BMDT) estimates for dose-related genomic markers stratified phthalates according to hypothetical tumorigenic potencies, unlike BMDs for non-genomic endpoints (relative liver weights or proliferation). The 7-day BMDT values for Acot1 as a surrogate measure for PPARα activation were 29, 370, and 676 mg/kg/day for DEHP, DNOP, and BBP, respectively, distinguishing DEHP (liver tumor BMD of 35 mg/kg/day) from non-tumorigenic DNOP and BBP. Effect thresholds were generated using linear regression of DEHP effects at 7 days and 2-year tumor incidence values to anchor early response molecular indicators and a later phenotypic outcome. Thresholds varied widely by marker, from 2-fold (Pdk4 and proliferation LI) to 30-fold (Acot1) induction to reach hypothetical tumorigenic expression levels. These findings highlight key issues in defining thresholds for biological adversity based on molecular changes.

An integrated pharmacokinetic and pharmacodynamic study of arsenite action 2. Heme oxygenase induction in mice.

Heme oxygenase (HO) is the rate-limiting enzyme in heme degradation and its activity has a significant impact on intracellular heme pools. Rat studies indicate that HO induction is a sensitive, dose-dependent response to arsenite (As(III)) exposure in both liver and kidney. The objective of this study was to evaluate the relationship of HO induction to administered As(III) dose, and concentrations of inorganic arsenic (iAs) in tissues and urine. Levels of iAs, mono- (MMA) and dimethylated arsenic (DMA) as well as HO activity were determined in liver, lung and kidney over time in female B6C3F1 mice given a single oral dose of 0, 1, 10, 30 or 100 micromol/kg As(III). Increased HO activity was a time and dose-dependent response in liver and kidney, but not in lung. Activity peaked in the 4-6 h time range in liver and kidney with the responsiveness in liver being approximately 2- to 3-fold greater than kidney. The lowest observed effect levels (LOELs) in this study for HO induction are 30 and 100 micromol/kg, respectively, in liver and kidney. The predominant form of arsenic (As) was iAs in liver at all doses, whereas DMA was the predominant form of As in kidney at all doses. Three- to four-fold higher levels of iAs were achieved in liver compared to kidney. MMA was the least abundant form of As in liver and kidney, never exceeding more than 20% of the total As present. The concentration of iAs in tissue or urine demonstrated the strongest correlation with HO activity in both liver and kidney. Results of this study suggest that HO induction is a biomarker of effect that is specific for tissue iAs because a high, but nontoxic, acute dose of DMA (5220 micromol/kg) did not induce HO in mice. Thus, HO induction has potential for use as a biomarker of effect for inorganic arsenic exposure and may be used as an indicator response to further the development of a biologically-based dose response model for As.

Kinetic modeling of beta-chloroprene metabolism: II. The application of physiologically based modeling for cancer dose response analysis.

beta-Chloroprene (2-chloro-1,3-butadiene; CD), which is used in the synthesis of polychloroprene, caused significant incidences of several tumor types in B6C3F1 mice and Fischer rats, but not in Wistar rats or Syrian hamsters. This project investigates the relevance of the bioassay lung tumor findings to human health risk by developing a physiologically based toxicokinetic (PBTK) model and exploring a tissue specific exposure-dose-response relationship. Key steps included identification of the plausible genotoxic mode of action, experimental quantification of tissue-to-air partition coefficients, scaling of in vitro parameters of CD metabolism for input into the PBTK model, comparing the model with in vivo experimental gas uptake data, selecting an appropriate tissue dosimetric, and predicting a corresponding human exposure concentration. The total daily milligram amount of CD metabolized per gram of lung was compared with the animal bioassay response data, specifically combined bronchiolar adenoma/carcinoma. The faster rate of metabolism in mouse lung agreed with the markedly greater incidence of lung tumors compared with the other rodent species. A lung tissue dose was predicted for the combined rodent lung tumor bioassay data at a 10% benchmark response. A human version of the PBTK model predicted that the lung tissue dose in humans would be equivalent to continuous lifetime daily exposure of 23 ppm CD. PBTK model sensitivity analysis indicated greater dependence of model predictions of dosimetry on physiological than biochemical parameters. The combined analysis of lung tumor response across species using the PBTK-derived internal dose provides an improved alternative to default pharmacokinetic interspecies adjustments for application to human health risk assessment.