Search

Chemically induced steatosis is characterized by lipid accumulation associated with mitochondrial dysfunction, oxidative stress and nucleus distortion. New approach methods integrating in vitro and in silico models are needed to identify chemicals that may induce these cellular events as potential risk factors for steatosis and associated hepatotoxicity. In this study we used high-content imaging for the simultaneous quantification of four cellular markers as sentinels for hepatotoxicity and steatosis in chemically exposed human liver cells in vitro. Furthermore, we evaluated the results with a computational model for the extrapolation of human oral equivalent doses (OED). First, we tested 16 reference chemicals with known capacities to induce cellular alterations in nuclear morphology, lipid accumulation, mitochondrial membrane potential and oxidative stress. Then, using physiologically based pharmacokinetic modeling and reverse dosimetry, OEDs were extrapolated from data of any stimulated individual sentinel response. The extrapolated OEDs were confirmed to be within biologically relevant exposure ranges for the reference chemicals. Next, we tested 14 chemicals found in food, selected from thousands of putative chemicals on the basis of structure-based prediction for nuclear receptor activation. Amongst these, orotic acid had an extrapolated OED overlapping with realistic exposure ranges. Thus, we were able to characterize known steatosis-inducing chemicals as well as data-scarce food-related chemicals, amongst which we confirmed orotic acid to induce hepatotoxicity. This strategy addresses needs of next generation risk assessment and can be used as a first chemical prioritization hazard screening step in a tiered approach to identify chemical risk factors for steatosis and hepatotoxicity-associated events., Archives of Toxicology, 97 (6), ISSN:0340-5761, ISSN:1432-0738

In silico methods to estimate and/or quantify skin absorption of chemicals as a function of chemistry are needed to realistically predict pharmacological, occupational, and environmental exposures. The Potts-Guy equation is a well- established approach, using multi-linear regression analysis describing skin permeability (Kp) in terms of the octanol/water partition coefficient (logP) and molecular weight (MW). In this work, we obtained regression equations for different human datasets relevant to environmental and cosmetic chemicals. Since the Potts-Guy equation was published in 1992, we explored recent datasets that include different skin layers, such as dermatomed (including dermis to a defined thickness) and full-skin. Our work was consistent with others who have observed that fits to the Potts-Guy equation are stronger for experiments focused on the epidermis. Permeability estimates for dermatomed skin and full skin resulted in low regression coefficients when compared to epidermis datasets. An updated regression equation uses a combination of fitted permeability values obtained with a published 2D compartmental model previously evaluated. The resulting regression equation was: logKp = = -2.55 +0.65*logP - 0.0085*MW, R

Recycled materials are found in many consumer products as part of a circular economy; however, the chemical content of recycled products is generally uncharacterized. A suspect screening analysis using two-dimensional gas chromatography time-of-flight mass spectrometry (GC × GC-TOFMS) was applied to 210 products (154 recycled, 56 virgin) across seven categories. Chemicals in products were tentatively identified using a standard spectral library or confirmed using chemical standards. A total of 918 probable chemical structures identified (112 of which were confirmed) in recycled materials versus 587 (110 confirmed) in virgin materials. Identified chemicals were characterized in terms of their functional use and structural class. Recycled paper products and construction materials contained greater numbers of chemicals than virgin products; 733 identified chemicals had greater occurrence in recycled compared to virgin materials. Products made from recycled materials contained greater numbers of fragrances, flame retardants, solvents, biocides, and dyes. The results were clustered to identify groups of chemicals potentially associated with unique chemical sources, and identified chemicals were prioritized for further study using high-throughput hazard and exposure information. While occurrence is not necessarily indicative of risk, these results can be used to inform the expansion of existing models or identify exposure pathways currently neglected in exposure assessments.

INTRODUCTION: Toxicity data are unavailable for many thousands of chemicals in commerce and the environment. Therefore, risk assessors need to rapidly screen these chemicals for potential risk to public health. High-throughput screening (HTS) for in vitro bioactivity, when used with high-throughput toxicokinetic (HTTK) data and models, allows characterization of these thousands of chemicals. AREAS COVERED: This review covers generic physiologically-based toxicokinetic (PBTK) models and high-throughput PBTK modeling for in vitro-in vivo extrapolation (IVIVE) of HTS data. We focus on “httk”, a public, open-source set of computational modeling tools and in vitro toxicokinetic (TK) data. EXPERT OPINION: HTTK benefits chemical risk assessors with its ability to support rapid chemical screening/prioritization, perform IVIVE, and provide provisional TK modeling for large numbers of chemicals using only limited chemical-specific data. Although generic TK model design can increase prediction uncertainty, these models provide offsetting benefits by increasing model implementation accuracy. Also, public distribution of the models and data enhances reproducibility. For the httk package, the modular and open-source design can enable the tool to be used and continuously improved by a broad user community in support of the critical need for high-throughput chemical prioritization and rapid dose estimation to facilitate rapid hazard assessments.

Toxicokinetic (TK) data needed for chemical risk assessment are not available for most chemicals. To support a greater number of chemicals, the U.S. Environmental Protection Agency (EPA) created the open-source R package "httk" (High Throughput ToxicoKinetics). The "httk" package provides functions and data tables for simulation and statistical analysis of chemical TK, including a population variability simulator that uses biometrics data from the National Health and Nutrition Examination Survey (NHANES).Here we modernize the "HTTK-Pop" population variability simulator based on the currently available data and literature. We provide explanations of the algorithms used by "httk" for variability simulation and uncertainty propagation.We updated and revised the population variability simulator in the "httk" package with the most recent NHANES biometrics (up to the 2017-18 NHANES cohort). Model equations describing glomerular filtration rate (GFR) were revised to more accurately represent physiology and population variability. The model output from the updated "httk" package was compared with the current version.The revised population variability simulator in the "httk" package now provides refined, more relevant, and better justified estimations.Fulfilling the U.S. EPA's mission to provide open-source data and models for evaluations and applications by the broader scientific community, and continuously improving the accuracy of the "httk" package based on the currently available data and literature.

Currently it is difficult to prospectively estimate human toxicokinetics (particularly for novel chemicals) in a high-throughput manner. The R software package httk has been developed, in part, to address this deficiency, and the aim of this investigation was to develop a generalized inhalation model for httk. The structure of the inhalation model was developed from two previously published physiologically based models from Jongeneelen and Berge (Ann Occup Hyg 55:841–864, 2011) and Clewell et al. (Toxicol Sci 63:160–172, 2001), while calculated physicochemical data was obtained from EPA’s CompTox Chemicals Dashboard. In total, 142 exposure scenarios across 41 volatile organic chemicals were modeled and compared to published data. The slope of the regression line of best fit between log-transformed simulated and observed blood and exhaled breath concentrations was 0.46 with an r2 = 0.45 and a root mean square error (RMSE) of direct comparison between the log-transformed simulated and observed values of 1.11. Approximately 5.1% (n = 108) of the data points analyzed were >2 orders of magnitude different than expected. The volatile organic chemicals examined in this investigation represent small, generally lipophilic molecules. Ultimately this paper details a generalized inhalation component that integrates with the httk physiologically based toxicokinetic model to provide high-throughput estimates of inhalation chemical exposures.

Estimates of exposure are critical to prioritize and assess chemicals based on risk posed to public health and the environment. The U.S. Environmental Protection Agency (EPA) is responsible for regulating thousands of chemicals in commerce and the environment for which exposure data are limited. Since 2009 the EPA's ExpoCast ("Exposure Forecasting") project has sought to develop the data, tools, and evaluation approaches required to generate rapid and scientifically defensible exposure predictions for the full universe of existing and proposed commercial chemicals. This review article aims to summarize issues in exposure science that have been addressed through initiatives affiliated with ExpoCast. ExpoCast research has generally focused on chemical exposure as a statistical systems problem intended to inform thousands of chemicals. The project exists as a companion to EPA's ToxCast ("Toxicity Forecasting") project which has used in vitro high-throughput screening technologies to characterize potential hazard posed by thousands of chemicals for which there are limited toxicity data. Rapid prediction of chemical exposures and in vitro-in vivo extrapolation (IVIVE) of ToxCast data allow for prioritization based upon risk of adverse outcomes due to environmental chemical exposure. ExpoCast has developed (1) integrated modeling approaches to reliably predict exposure and IVIVE dose, (2) highly efficient screening tools for chemical prioritization, (3) efficient and affordable tools for generating new exposure and dose data, and (4) easily accessible exposure databases. The development of new exposure models and databases along with the application of technologies like non-targeted analysis and machine learning have transformed exposure science for data-poor chemicals. By developing high-throughput tools for chemical exposure analytics and translating those tools into public health decisions ExpoCast research has served as a crucible for identifying and addressing exposure science knowledge gaps.

1AbstractNonalcoholic fatty liver disease (NAFLD) is an increasingly prevalent human disease with accumulating evidence linking its pathophysiology and co-morbidities to chemical exposures. The complex pathophysiology of NAFLD has limited the elucidation of potential chemical etiologies. In this study we generated a high-content imaging analysis method for the simultaneous quantification of sentinel steatosis cellular markers in chemically exposed human liver cells in vitro combined with a computational model for the extrapolation of human oral equivalent doses (OED). First, the in vitro test method was generated using 14 reference chemicals with known capacities to induce cellular alterations in nuclear morphology, lipid accumulation, mitochondrial membrane potential and oxidative stress. These effects were quantified on a single cell- and population-level, and then, using physiologically based pharmacokinetic modelling and reverse dosimetry, OEDs were extrapolated from these in vitro data. The extrapolated OEDs were confirmed to be within biologically relevant exposure ranges for the reference chemicals. Next, we tested 14 chemicals found in food, selected from thousands of putative chemicals on the basis of structure-based prediction for nuclear receptor activation. Amongst these, orotic acid had an extrapolated OED overlapping with realistic exposure ranges. By the strategy developed in this study, we were able to characterize known NAFLD-inducing chemicals and translate to data scarce food-related chemicals, amongst which we identified orotic acid to induce steatosis. This strategy addresses needs of next generation risk assessment, and can be used as a first chemical prioritization hazard screening step in a tiered approach to identify chemical risk factors for NAFLD.

Chemical risk assessment considers potentially susceptible populations including pregnant women and developing fetuses. Humans encounter thousands of chemicals in their environments, few of which have been fully characterized. Toxicokinetic (TK) information is needed to relate chemical exposure to potentially bioactive tissue concentrations. Observational data describing human gestational exposures are unavailable for most chemicals, but physiologically based TK (PBTK) models estimate such exposures. Development of chemical-specific PBTK models requires considerable time and resources. As an alternative, generic PBTK approaches describe a standardized physiology and characterize chemicals with a set of standard physical and TK descriptors - primarily plasma protein binding and hepatic clearance. Here we report and evaluate a generic PBTK model of a human mother and developing fetus. We used a published set of formulas describing the major anatomical and physiological changes that occur during pregnancy to augment the High-Throughput Toxicokinetics (httk) software package. We simulated the ratio of concentrations in maternal and fetal plasma and compared to literature in vivo measurements. We evaluated the model with literature in vivo time-course measurements of maternal plasma concentrations in pregnant and non-pregnant women. Finally, we prioritized chemicals measured in maternal serum based on predicted fetal brain concentrations. This new model can be used for TK simulations of 859 chemicals with existing human-specific in vitro TK data as well as any new chemicals for which such data become available. This gestational model may allow for in vitro to in vivo extrapolation of point of departure doses relevant to reproductive and developmental toxicity.

In the real world, individuals are exposed to chemicals from sources that vary over space and time. However, traditional risk assessments based on in vivo animal studies typically use a chemical-by-chemical approach and apical disease endpoints. New approach methodologies (NAMs) in toxicology, such as in vitro high-throughput (HTS) assays generated in Tox21 and ToxCast, can more readily provide mechanistic chemical hazard information for chemicals with no existing data than in vivo methods. In this paper, we establish a workflow to assess the joint action of 41 modeled ambient chemical exposures in the air from the USA-wide National Air Toxics Assessment by integrating human exposures with hazard data from curated HTS (cHTS) assays to identify counties where exposure to the local chemical mixture may perturb a common biological target. We exemplify this proof-of-concept using CYP1A1 mRNA up-regulation. We first estimate internal exposure and then convert the inhaled concentration to a steady state plasma concentration using physiologically based toxicokinetic modeling parameterized with county-specific information on ages and body weights. We then use the estimated blood plasma concentration and the concentration-response curve from the in vitro cHTS assay to determine the chemical-specific effects of the mixture components. Three mixture modeling methods were used to estimate the joint effect from exposure to the chemical mixture on the activity levels, which were geospatially mapped. Finally, a Monte Carlo uncertainty analysis was performed to quantify the influence of each parameter on the combined effects. This workflow demonstrates how NAMs can be used to predict early-stage biological perturbations that can lead to adverse health outcomes that result from exposure to chemical mixtures. As a result, this work will advance mixture risk assessment and other early events in the effects of chemicals.

With the increasing availability of high-resolution mass spectrometers, suspect screening and non-targeted analysis are becoming popular compound identification tools for environmental researchers. Samples of interest often contain a large (unknown) number of chemicals spanning the detectable mass range of the instrument. In an effort to separate these chemicals prior to injection into the mass spectrometer, a chromatography method is often utilized. There are numerous types of gas and liquid chromatographs that can be coupled to commercially available mass spectrometers. Depending on the type of instrument used for analysis, the researcher is likely to observe a different subset of compounds based on the amenability of those chemicals to the selected experimental techniques and equipment. It would be advantageous if this subset of chemicals could be predicted prior to conducting the experiment, in order to minimize potential false-positive and false-negative identifications. In this work, we utilize experimental datasets to predict the amenability of chemical compounds to detection with liquid chromatography-electrospray ionization-mass spectrometry (LC-ESI-MS). The assembled dataset totals 5517 unique chemicals either explicitly detected or not detected with LC-ESI-MS. The resulting detected/not-detected matrix has been modeled using specific molecular descriptors to predict which chemicals are amenable to LC-ESI-MS, and to which form(s) of ionization. Random forest models, including a measure of the applicability domain of the model for both positive and negative modes of the electrospray ionization source, were successfully developed. The outcome of this work will help to inform future suspect screening and non-targeted analyses of chemicals by better defining the potential LC-ESI-MS detectable chemical landscape of interest.

Advances in the biological sciences have led to an ongoing paradigm shift in toxicity testing based on expanded application of high-throughput in vitro screening and in silico methods to assess potential health risks of environmental agents. This review examines progress on the vision for toxicity testing elaborated by the US National Research Council (NRC) during the decade that has passed since the 2007 NRC report on Toxicity Testing in the 21st Century (TT21C). Concomitant advances in exposure assessment, including computational approaches and high-throughput exposomics, are also documented. A vision for the next generation of risk science, incorporating risk assessment methodologies suitable for the analysis of new toxicological and exposure data, resulting in human exposure guidelines is described. Case study prototypes indicating how these new approaches to toxicity testing, exposure measurement, and risk assessment are beginning to be applied in practice are presented. Overall, progress on the 20-year transition plan laid out by the US NRC in 2007 has been substantial. Importantly, government agencies within the United States and internationally are beginning to incorporate the new approach methodologies envisaged in the original TT21C vision into regulatory practice. Future perspectives on the continued evolution of toxicity testing to strengthen regulatory risk assessment are provided.

High(er) throughput toxicokinetics (HTTK) encompasses in vitro measures of key determinants of chemical toxicokinetics and reverse dosimetry approaches for in vitro-in vivo extrapolation (IVIVE). With HTTK, the bioactivity identified by any in vitro assay can be converted to human equivalent doses and compared with chemical intake estimates. Biological variability in HTTK has been previously considered, but the relative impact of measurement uncertainty has not. Bayesian methods were developed to provide chemical-specific uncertainty estimates for 2 in vitro toxicokinetic parameters: unbound fraction in plasma (fup) and intrinsic hepatic clearance (Clint). New experimental measurements of fup and Clint are reported for 418 and 467 chemicals, respectively. These data raise the HTTK chemical coverage of the ToxCast Phase I and II libraries to 57%. Although the standard protocol for Clint was followed, a revised protocol for fup measured unbound chemical at 10%, 30%, and 100% of physiologic plasma protein concentrations, allowing estimation of protein binding affinity. This protocol reduced the occurrence of chemicals with fup too low to measure from 44% to 9.1%. Uncertainty in fup was also reduced, with the median coefficient of variation dropping from 0.4 to 0.1. Monte Carlo simulation was used to propagate both measurement uncertainty and biological variability into IVIVE. The uncertainty propagation techniques used here also allow incorporation of other sources of uncertainty such as in silico predictors of HTTK parameters. These methods have the potential to inform risk-based prioritization based on the relationship between in vitro bioactivities and exposures.

Acceptance of new approach methodologies (NAMs) for use in characterizing chemical hazard and risk requires informed expectations regarding the minimum precision and maximum accuracy of their results. Uncertainty in NAMs derived from variability in the traditional reference data used to train or validate performance of the NAM, and uncertainty in the modeling procedures themselves, limits NAM performance. Herein, we review current approaches to characterizing uncertainty in NAMs. We discuss variability in in vivo data used as a reference for NAM development and validation; the quantitative uncertainty in concentration–response modeling for high-throughput in vitro bioactivity screening; the uncertainties associated with in vitro to in vivo extrapolation using toxicokinetic information; and the quantitative uncertainty in the experimental inputs and modeled outputs from quantitative structure activity relationship models for prediction of point of departure doses. Communication of the amount of uncertainty, both in the input and output for NAMs, often involves a confidence or prediction interval around a given potency estimate, derived from an understanding of the variability in the data modeled. Tuning expectations of NAM performance to an understanding of the reproducibility and variability, both of traditional approaches and NAM approaches, provides a path for the adoption of NAMs as alternatives in screening chemicals for risk.

Chemical risk assessment relies on knowledge of hazard, the dose–response relationship, and exposure to characterize potential risks to public health and the environment. A chemical with minimal toxicity might pose a risk if exposures are extensive, repeated, and/or occurring during critical windows across the human life span. Exposure assessment involves understanding human activity, and this activity is confounded by interindividual variability that is both biological and behavioral. Exposures further vary between the general population and susceptible or occupationally exposed populations. Recent computational exposure efforts have tackled these problems through the creation of new tools and predictive models. These tools include machine learning to draw inferences from existing data and computer-enhanced screening analyses to generate new data. Mathematical models provide frameworks describing chemical exposure processes. These models can be statistically evaluated to establish rigorous confidence in their predictions. The computational exposure tools reviewed here are oriented toward ‘high-throughput’ application, that is, they are suitable for dealing with the thousands of chemicals in commerce with limited sources of chemical exposure information. These new tools and models are moving chemical exposure and risk assessment forward in the 21st century.

The intrinsic metabolic clearance rate (Cl(int)) and fraction of chemical unbound in plasma (f(up)) serve as important parameters for high throughput toxicokinetic models, but experimental data are limited for many chemicals. Open-source quantitative structure-activity relationship (QSAR) models for both parameters were developed to offer reliable in silico predictions for a diverse set of chemicals regulated under U.S. law, including pharmaceuticals, pesticides, and industrial chemicals. As a case study to demonstrate their utility, model predictions served as inputs to the TK component of a risk-based prioritization approach based on Bioactivity: Exposure Ratios (BER), in which a BER < 1 indicates exposures are predicted to exceed a biological activity threshold. When applied to a subset of the Tox21 screening library (6484 chemicals) we found that the proportion of chemicals with BER 1 using either in silico or in vitro parameters (768/848, 90.5%). Thus, the presented QSARs may be suitable for prioritizing the risk posed by many chemicals for which measured in vitro TK data are lacking.

The intrinsic metabolic clearance rate (Clint) and fraction of chemical unbound in plasma (fup) serve as important parameters for high throughput toxicokinetic models, but experimental data are limited for many chemicals. Open-source quantitative structure-activity relationship (QSAR) models for both parameters were developed to offer reliable in silico predictions for a diverse set of chemicals regulated under U.S. law, including pharmaceuticals, pesticides, and industrial chemicals. As a case study to demonstrate their utility, model predictions served as inputs to the TK component of a risk-based prioritization approach based on Bioactivity: Exposure Ratios (BER), in which a BER 1 using either in silico or in vitro parameters (776/850, 91.30%). Thus, the presented QSARs may be suitable for prioritizing the risk posed by many chemicals for which measured in vitro TK data are lacking.

The toxicokinetic (TK) parameters fraction of the chemical unbound to plasma proteins and metabolic clearance are critical for relating exposure and internal dose when building in vitro-based risk assessment models. However, experimental toxicokinetic studies have only been carried out on limited chemicals of environmental interest (~1000 chemicals with TK data relative to tens of thousands of chemicals of interest). This work evaluated the utility of chemical structure information to predict TK parameters in silico; development of cluster-based read-across and quantitative structure-activity relationship models of fraction unbound or fub (regression) and intrinsic clearance or Cl(int) (classification and regression) using a dataset of 1487 chemicals; utilization of predicted TK parameters to estimate uncertainty in steady-state plasma concentration (C(ss)); and subsequent in vitro–in vivo extrapolation analyses to derive bioactivity-exposure ratio (BER) plot to compare human oral equivalent doses and exposure predictions using androgen and estrogen receptor activity data for 233 chemicals as an example dataset. The results demonstrate that fub is structurally more predictable than Cl(int). The model with the highest observed performance for fub had an external test set RMSE/σ=0.62 and R(2)=0.61, for Cl(int) classification had an external test set accuracy = 65.9%, and for intrinsic clearance regression had an external test set RMSE/σ=0.90 and R(2)=0.20. This relatively low performance is in part due to the large uncertainty in the underlying Cl(int) data. We show that C(ss) is relatively insensitive to uncertainty in Cl(int). The models were benchmarked against the ADMET Predictor software. Finally, the BER analysis allowed identification of 14 out of 136 chemicals for further risk assessment demonstrating the utility of these models in aiding risk-based chemical prioritization.

In vitro chemical safety testing methods offer the potential for efficient and economical tools to provide relevant assessments of human health risk. To realize this potential, methods are needed to relate in vitro effects to in vivo responses, i.e., in vitro to in vivo extrapolation (IVIVE). Currently available IVIVE approaches need to be refined before they can be utilized for regulatory decision-making. To explore the capabilities and limitations of IVIVE within this context, the U.S. Environmental Protection Agency Office of Research and Development and the National Toxicology Program Interagency Center for the Evaluation of Alternative Toxicological Methods co-organized a workshop and webinar series. Here, we integrate content from the webinars and workshop to discuss activities and resources that would promote inclusion of IVIVE in regulatory decision-making. We discuss properties of models that successfully generate predictions of in vivo doses from effective in vitro concentration, including the experimental systems that provide input parameters for these models, areas of success, and areas for improvement to reduce model uncertainty. Finally, we provide case studies on the uses of IVIVE in safety assessments, which highlight the respective differences, information requirements, and outcomes across various approaches when applied for decision-making.

Chemical ionization plays an important role in many aspects of pharmacokinetic (PK) processes such as protein binding, tissue partitioning, and apparent volume of distribution at steady state (Vdss). Here, estimates of ionization equilibrium constants (i.e., pKa) were analyzed for 8,132 pharmaceuticals and 24,281 other compounds to which humans might be exposed in the environment. Results revealed broad differences in the ionization of pharmaceutical chemicals and chemicals with either near-field (in the home) or far-field sources. The utility of these high-throughput ionization predictions was evaluated via a case-study of predicted PK Vdss for 22 compounds monitored in the blood and serum of the U.S. population by the U.S. Centers for Disease Control and Prevention National Health and Nutrition Examination Survey (NHANES). The chemical distribution ratio between water and tissue was estimated using predicted ionization states characterized by pKa. Probability distributions corresponding to ionizable atom types (IATs) were then used to analyze the sensitivity of predicted Vdss on predicted pKa using Monte Carlo methods. 8 of the 22 compounds were predicted to be ionizable. For 5 of the 8 the predictions based upon ionization are significantly different from what would be predicted for a neutral compound. For all but one (foramsulfuron), the probability distribution of predicted Vdss generated by IAT sensitivity analysis spans both the neutral prediction and the prediction using ionization. As new data sets of chemical-specific information on metabolism and excretion for hundreds of chemicals are being made available (e.g., Wetmore et al., 2015), high-throughput methods for calculating Vdss and tissue-specific PK distribution coefficients will allow the rapid construction of PK models to provide context for both biomonitoring data and high-throughput toxicity screening studies such as Tox21 and ToxCast.

Human health risks from chronic exposures to environmental chemicals are typically estimated from potential human exposure estimates and dose-response data obtained from repeated-dose animal toxicity studies. Various criteria are available for selecting the top (highest) dose used in these animal studies. For example, toxicokinetic (TK) and toxicological data provided by shorter-term or dose range finding studies can be evaluated in a weight of evidence approach to provide insight into the dose range that would provide dose-response data that are relevant to human exposures. However, there are concerns that a top dose resulting from the consideration of TK data may be too low compared to other criteria, such as the limit dose or the maximum tolerated dose. In this paper, we address several concerns related to human exposures by discussing 1) the resources and methods available to predict human exposure levels and the associated uncertainty and variability, and 2) the margin between predicted human exposure levels and the dose levels used in repeated-dose animal studies. A series of case studies, ranging from data-rich to data-poor chemicals, are presented to demonstrate that expected human exposures to environmental chemicals are typically orders of magnitude lower than no-observed-adverse-effect levels/lowest-observed-adverse-effect levels (NOAELs/LOAELs) when available (used as conservative surrogates for top doses). The results of these case studies support that a top dose based, in part, on TK data is typically orders of magnitude higher than expected human exposure levels.

Top dose selection for repeated dose animal studies has generally focused on identification of apical endpoints, use of the limit dose, or determination of a maximum tolerated dose (MTD). The intention is to optimize the ability of toxicity tests performed in a small number of animals to detect effects for hazard identification. An alternative approach, the kinetically derived maximum dose (KMD), has been proposed as a mechanism to integrate toxicokinetic (TK) data into the dose selection process. The approach refers to the dose above which the systemic exposures depart from being proportional to external doses. This non-linear external-internal dose relationship arises from saturation or limitation of TK process(es), such as absorption or metabolism. The importance of TK information is widely acknowledged when assessing human health risks arising from exposures to environmental chemicals, as TK determines the amount of chemical at potential sites of toxicological responses. However, there have been differing opinions and interpretations within the scientific and regulatory communities related to the validity and application of the KMD concept. A multi-stakeholder working group, led by the Health and Environmental Sciences Institute (HESI), was formed to provide an opportunity for impacted stakeholders to address commonly raised scientific and technical issues related to this topic and, more specifically, a weight of evidence approach is recommended to inform design and dose selection for repeated dose animal studies. Commonly raised challenges related to the use of TK data for dose selection are discussed, recommendations are provided, and illustrative case examples are provided to address these challenges or refute misconceptions.

Tens-of-thousands of chemicals are registered in the U.S. for use in countless processes and products. Recent evidence suggests that many of these chemicals are measureable in environmental and/or biological systems, indicating the potential for widespread exposures. Traditional public health research tools, including in vivo studies and targeted analytical chemistry methods, have been unable to meet the needs of screening programs designed to evaluate chemical safety. As such, new tools have been developed to enable rapid assessment of potentially harmful chemical exposures and their attendant biological responses. One group of tools, known as "non-targeted analysis" (NTA) methods, allows the rapid characterization of thousands of never-before-studied compounds in a wide variety of environmental, residential, and biological media. This article discusses current applications of NTA methods, challenges to their effective use in chemical screening studies, and ways in which shared resources (e.g., chemical standards, databases, model predictions, and media measurements) can advance their use in risk-based chemical prioritization. A brief review is provided of resources and projects within EPA's Office of Research and Development (ORD) that provide benefit to, and receive benefits from, NTA research endeavors. A summary of EPA's Non-Targeted Analysis Collaborative Trial (ENTACT) is also given, which makes direct use of ORD resources to benefit the global NTA research community. Finally, a research framework is described that shows how NTA methods will bridge chemical prioritization efforts within ORD. This framework exists as a guide for institutions seeking to understand the complexity of chemical exposures, and the impact of these exposures on living systems.

Under the ExpoCast program, United States Environmental Protection Agency (EPA) researchers have developed a high-throughput (HT) framework for estimating aggregate exposures to chemicals from multiple pathways to support rapid prioritization of chemicals. Here, we present methods to estimate HT exposures to chemicals migrating into food from food contact substances (FCS). These methods consisted of combining an empirical model of chemical migration with estimates of daily population food intakes derived from food diaries from the National Health and Nutrition Examination Survey (NHANES). A linear regression model for migration at equilibrium was developed by fitting available migration measurements as a function of temperature, food type (i.e., fatty, aqueous, acidic, alcoholic), initial chemical concentration in the FCS (C0) and chemical properties. The most predictive variables in the resulting model were C0, molecular weight, log Kow, and food type (R2=0.71, p

Toxicokinetics (TK) provides critical information for integrating chemical toxicity and exposure assessments in order to determine potential chemical risk (i.e., the margin between toxic doses and plausible exposures). For thousands of chemicals that are present in our environment, in vivo TK data are lacking. The publicly available R package “httk” (version 1.8, named for “high throughput TK”) draws from a database of in vitro data and physico-chemical properties in order to run physiologically-based TK (PBTK) models for 553 compounds. The PBTK model parameters include tissue:plasma partition coefficients (Kp) which the httk software predicts using the model of Schmitt (Toxicol In Vitro 22 (2):457–467, 2008). In this paper we evaluated and modified httk predictions, and quantified confidence using in vivo literature data. We used 964 rat Kp measured by in vivo experiments for 143 compounds. Initially, predicted Kp were significantly larger than measured Kp for many lipophilic compounds (log10 octanol:water partition coefficient > 3). Hence the approach for predicting Kp was revised to account for possible deficiencies in the in vitro protein binding assay, and the method for predicting membrane affinity was revised. These changes yielded improvements ranging from a factor of 10 to nearly a factor of 10,000 for 83 Kp across 23 compounds with only 3 Kp worsening by more than a factor of 10. The vast majority (92%) of Kp were predicted within a factor of 10 of the measured value (overall root mean squared error of 0.59 on log10-transformed scale). After applying the adjustments, regressions were performed to calibrate and evaluate the predictions for 12 tissues. Predictions for some tissues (e.g., spleen, bone, gut, lung) were observed to be better than predictions for other tissues (e.g., skin, brain, fat), indicating that confidence in the application of in silico tools to predict chemical partitioning varies depending upon the tissues involved. Our calibrated model was then evaluated using a second data set of human in vivo measurements of volume of distribution (Vss) for 498 compounds reviewed by Obach et al. (Drug Metab Dispos 36(7):1385–1405, 2008). We found that calibration of the model improved performance: a regression of the measured values as a function of the predictions has a slope of 1.03, intercept of − 0.04, and R2 of 0.43. Through careful evaluation of predictive methods for chemical partitioning into tissues, we have improved and calibrated these methods and quantified confidence for TK predictions in humans and rats.

The thousands of chemicals present in the environment (USGAO, 2013) must be triaged to identify priority chemicals for human health risk research. Most chemicals have little of the toxicokinetic (TK) data that are necessary for relating exposures to tissue concentrations that are believed to be toxic. Ongoing efforts have collected limited, in vitro TK data for a few hundred chemicals. These data have been combined with biomonitoring data to estimate an approximate margin between potential hazard and exposure. The most “at risk” 95th percentile of adults have been identified from simulated populations that are generated either using standard “average” adult human parameters or very specific cohorts such as Northern Europeans. To better reflect the modern U.S. population, we developed a population simulation using physiologies based on distributions of demographic and anthropometric quantities from the most recent U.S. Centers for Disease Control and Prevention National Health and Nutrition Examination Survey (NHANES) data. This allowed incorporation of inter-individual variability, including variability across relevant demographic subgroups. Variability was analyzed with a Monte Carlo approach that accounted for the correlation structure in physiological parameters. To identify portions of the U.S. population that are more at risk for specific chemicals, physiologic variability was incorporated within an open-source high-throughput (HT) TK modeling framework. We prioritized 50 chemicals based on estimates of both potential hazard and exposure. Potential hazard was estimated from in vitro HT screening assays (i.e., the Tox21 and ToxCast programs). Bioactive in vitro concentrations were extrapolated to doses that produce equivalent concentrations in body tissues using a reverse dosimetry approach in which generic TK models are parameterized with: 1) chemical-specific parameters derived from in vitro measurements and predicted from chemical structure; and 2) with physiological parameters for a virtual population. For risk-based prioritization of chemicals, predicted bioactive equivalent doses were compared to demographic-specific inferences of exposure rates that were based on NHANES urinary analyte biomonitoring data. The inclusion of NHANES-derived inter-individual variability decreased predicted bioactive equivalent doses by 12% on average for the total population when compared to previous methods. However, for some combinations of chemical and demographic groups the margin was reduced by as much as three quarters. This TK modeling framework allows targeted risk prioritization of chemicals for demographic groups of interest, including potentially sensitive life stages and subpopulations. Keywords: Risk assessment, Toxicokinetics, High throughput, Environmental chemicals, IVIVE