Considerations for Improving Metabolism Predictions for In Vitro to In Vivo Extrapolation.

High-throughput (HT) in vitro to in vivo extrapolation (IVIVE) is an integral component in new approach method (NAM)-based risk assessment paradigms, for rapidly translating in vitro toxicity assay results into the context of in vivo exposure. When coupled with rapid exposure predictions, HT-IVIVE supports the use of HT in vitro assays for risk-based chemical prioritization. However, the reliability of prioritization based on HT bioactivity data and HT-IVIVE can be limited as the domain of applicability of current HT-IVIVE is generally restricted to intrinsic clearance measured primarily in pharmaceutical compounds. Further, current approaches only consider parent chemical toxicity. These limitations occur because current state-of-the-art HT prediction tools for clearance and metabolite kinetics do not provide reliable data to support HT-IVIVE. This paper discusses current challenges in implementation of IVIVE for prioritization and risk assessment and recommends a path forward for addressing the most pressing needs and expanding the utility of IVIVE.

Application of a combined aggregate exposure pathway and adverse outcome pathway (AEP-AOP) approach to inform a cumulative risk assessment: A case study with phthalates.

Advancements in measurement and modeling capabilities are providing unprecedented access to estimates of chemical exposure and bioactivity. With this influx of new data, there is a need for frameworks that help organize and disseminate information on chemical hazard and exposure in a manner that is accessible and transparent. A case study approach was used to demonstrate integration of the Adverse Outcome Pathway (AOP) and Aggregate Exposure Pathway (AEP) frameworks to support cumulative risk assessment of co-exposure to two phthalate esters that are ubiquitous in the environment and that are associated with disruption of male sexual development in the rat: di(2-ethylhexyl) phthalate (DEHP) and di-n-butyl phthalate (DnBP). A putative AOP was developed to guide selection of an in vitro assay for derivation of bioactivity values for DEHP and DnBP and their metabolites. AEPs for DEHP and DnBP were used to extract key exposure data as inputs for a physiologically based pharmacokinetic (PBPK) model to predict internal metabolite concentrations. These metabolite concentrations were then combined using in vitro-based relative potency factors for comparison with an internal dose metric, resulting in an estimated margin of safety of ~13,000. This case study provides an adaptable workflow for integrating exposure and toxicity data by coupling AEP and AOP frameworks and using in vitro and in silico methodologies for cumulative risk assessment.

A weight of evidence approach to investigate potential common mechanisms in pesticide groups to support cumulative risk assessment: A case study with dinitroaniline pesticides.

In 2016, the United States Environmental Protection Agency's (EPA) Office of Pesticide Programs published guidelines for establishing candidate common mechanism groups (CMGs) for cumulative risk assessment (CRA) weight-of-evidence-based screenings. A candidate CMG is a group of chemicals that may share similar structure, apical endpoints, and/or mechanistic data that suggest the potential for a common mechanism of toxicity among them. Here, a weight-of-evidence approach is presented to establish candidacy of a CMG for a group of nine dinitroaniline pesticides. This approach involves review of available in vivo toxicity information and literature to determine mode of action, along with analyses of in vitro toxicity data and chemical structure. Despite structural similarity among some dinitroanilines and some shared target organs identified through toxicity observed in in vivo studies, there were no consistencies among groups, suggesting lack of a common mechanism when all analyses are considered together. For example, two structurally similar compounds with thyroid/liver in vivo effects were not found active in any Toxicity Forecaster (ToxCast) in vitro assays. The weight-of-evidence is insufficient to support the testable hypothesis that dinitroanilines could form a CMG, and highlights the importance of establishing a consensus among multiple lines of evidence prior to CRA.

A Workflow for Identifying Metabolically Active Chemicals to Complement in vitro Toxicity Screening.

The new paradigm of toxicity testing approaches involves rapid screening of thousands of chemicals across hundreds of biological targets through use of in vitro assays. Such assays may lead to false negatives when the complex metabolic processes that render a chemical bioactive in a living system are unable to be replicated in an in vitro environment. In the current study, a workflow is presented for complementing in vitro testing results with in silico and in vitro techniques to identify inactive parents that may produce active metabolites. A case study applying this workflow involved investigating the influence of metabolism for over 1,400 chemicals considered inactive across18 in vitro assays related to the estrogen receptor (ER) pathway. Over 7,500 first-generation and second-generation metabolites were generated for these in vitro inactive chemicals using an in silico software program. Next, a consensus model comprised of four individual quantitative structure activity relationship (QSAR) models was used to predict ER-binding activity for each of the metabolites. Binding activity was predicted for ~8-10% of metabolites in each generation, with these metabolites linked to 259 in vitro inactive parent chemicals. Metabolites were enriched in substructures consisting of alcohol, aromatic, and phenol bonds relative to their inactive parent chemicals, suggesting these features are potentially favorable for ER-binding. The workflow presented here can be used to identify parent chemicals that can be potentially bioactive, to aid confidence in high throughput risk screening.

Aggregate Exposure Pathways in Support of Risk Assessment.

Over time, risk assessment has shifted from establishing relationships between exposure to a single chemical and a resulting adverse health outcome, to evaluation of multiple chemicals and disease outcomes simultaneously. As a result, there is an increasing need to better understand the complex mechanisms that influence risk of chemical and non-chemical stressors, beginning at their source and ending at a biological endpoint relevant to human or ecosystem health risk assessment. Just as the Adverse Outcome Pathway (AOP) framework has emerged as a means of providing insight into mechanism-based toxicity, the exposure science community has seen the recent introduction of the Aggregate Exposure Pathway (AEP) framework. AEPs aid in making exposure data applicable to the FAIR (i.e., findable, accessible, interoperable, and reusable) principle, especially by (1) organizing continuous flow of disjointed exposure information;(2) identifying data gaps, to focus resources on acquiring the most relevant data; (3) optimizing use and repurposing of existing exposure data; and (4) facilitating interoperability among predictive models. Herein, we discuss integration of the AOP and AEP frameworks and how such integration can improve confidence in both traditional and cumulative risk assessment approaches.

Challenges Associated With Applying Physiologically Based Pharmacokinetic Modeling for Public Health Decision-Making.

The development and application of physiologically based pharmacokinetic (PBPK) models in chemical toxicology have grown steadily since their emergence in the 1980s. However, critical evaluation of PBPK models to support public health decision-making across federal agencies has thus far occurred for only a few environmental chemicals. In order to encourage decision-makers to embrace the critical role of PBPK modeling in risk assessment, several important challenges require immediate attention from the modeling community. The objective of this contemporary review is to highlight 3 of these challenges, including: (1) difficulties in recruiting peer reviewers with appropriate modeling expertise and experience; (2) lack of confidence in PBPK models for which no tissue/plasma concentration data exist for model evaluation; and (3) lack of transferability across modeling platforms. Several recommendations for addressing these 3 issues are provided to initiate dialog among members of the PBPK modeling community, as these issues must be overcome for the field of PBPK modeling to advance and for PBPK models to be more routinely applied in support of public health decision-making.

Refining the aggregate exposure pathway.

Advancements in measurement technologies and modeling capabilities continue to result in an abundance of exposure information, adding to that currently in existence. However, fragmentation within the exposure science community acts as an obstacle for realizing the vision set forth in the National Research Council's report on Exposure Science in the 21st century to consider exposures from source to dose, on multiple levels of integration, and to multiple stressors. The concept of an Aggregate Exposure Pathway (AEP) was proposed as a framework for organizing and integrating diverse exposure information that exists across numerous repositories and among multiple scientific fields. A workshop held in May 2016 followed introduction of the AEP concept, allowing members of the exposure science community to provide extensive evaluation and feedback regarding the framework's structure, key components, and applications. The current work briefly introduces topics discussed at the workshop and attempts to address key challenges involved in refining this framework. The resulting evolution in the AEP framework's features allows for facilitating acquisition, integration, organization, and transparent application and communication of exposure knowledge in a manner that is independent of its ultimate use, thereby enabling reuse of such information in many applications.

Investigating the state of physiologically based kinetic modelling practices and challenges associated with gaining regulatory acceptance of model applications.

Physiologically based kinetic (PBK) models are used widely throughout a number of working sectors, including academia and industry, to provide insight into the dosimetry related to observed adverse health effects in humans and other species. Use of these models has increased over the last several decades, especially in conjunction with emerging alternative methods to animal testing, such as in vitro studies and data-driven in silico quantitative-structure-activity-relationship (QSAR) predictions. Experimental information derived from these new approach methods can be used as input for model parameters and allows for increased confidence in models for chemicals that did not have in vivo data for model calibration. Despite significant advancements in good modelling practice (GMP) for model development and evaluation, there remains some reluctance among regulatory agencies to use such models during the risk assessment process. Here, the results of a survey disseminated to the modelling community are presented in order to inform the frequency of use and applications of PBK models in science and regulatory submission. Additionally, the survey was designed to identify a network of investigators involved in PBK modelling and knowledgeable of GMP so that they might be contacted in the future for peer review of PBK models, especially in regards to vetting the models to such a degree as to gain a greater acceptance for regulatory purposes.

Assessing toxicity of contaminants in riverine suspended sediments to freshwater mussels.

The Clinch River in Virginia and Tennessee, USA, is well known for its diverse native freshwater mussel assemblages; however, notable declines in mussel populations in recent decades have prompted much concern and subsequent research. The authors examined the toxicity of recently deposited sediments on juveniles of the freshwater mussel Epioblasma brevidens by collecting time-integrated sediment samples from the water column with sediment traps from 11 sites in the Clinch River basin, including 6 sites within an 88-km reach deemed a "mussel zone of decline." Mussels were exposed to the riverine sediments and to 3 control sediments for 28 d; survival, shell length, and biomass were then assessed. Sediment treatment (i.e., river location) had a significant effect on mussel survival (p < 0.01) and biomass (p = 0.02) but did not affect length (p = 0.37), and sediments from 2 of the tributaries were the most toxic. Inorganic and organic analyses of sediments indicated the presence of metals and polycyclic aromatic hydrocarbons at all sites. Manganese was negatively correlated with mussel survival and biomass, as was ammonia with survival and total organic carbon with biomass. Current land uses in the watershed indicate that fossil fuel mining and agriculture may be associated with elevated manganese and ammonia, respectively. The authors found that sediments collected with sediment traps over relatively short deployment durations can help elucidate recent contaminant influx and its potential for inducing toxicity in benthic organisms. Environ Toxicol Chem 2017;36:395-407. © 2016 SETAC.

Extending the toxicity-testing paradigm for freshwater mussels: Assessing chronic reproductive effects of the synthetic estrogen 17α-ethinylestradiol on the unionid mussel Elliptio complanata.

Surface water concentrations of the synthetic estrogen 17α-ethinylestradiol (EE2) as low as 1ng/L can cause adverse reproductive effects in fish under acute and chronic exposure conditions, whereas higher concentrations (> 5ng/L) in acute studies are necessary to elicit adverse effects in freshwater mussels. Prolonged chronic exposures of freshwater mussels to EE2 remain un-evaluated. An extended duration testing paradigm was used to examine reproductive and biochemical (carbohydrate, lipid, protein) effects of EE2 on the unionid mussel, Elliptio complanata, throughout its reproductive cycle. Mussels were exposed to a control and EE2 concentrations (5 and 50ng/L) in six discrete and sequential 28 d tests, and in one discrete and simultaneous 180 d test, from February through August. Foot protrusion and siphoning behavior were recorded daily, along with conglutinate releases and larval (glochidia) mortality. Gonad, hemolymph, and gonad fluid samples were taken for biochemical and vitellogenin-like protein (Vtg) analysis post-exposure. Female mussels released eggs and conglutinates during the months of April to June, indicating sexual maturation during this time. Conglutinates released in the 5ng/L treatment in 28 d exposures contained fewer glochidia and more eggs, and increased concentrations of Vtg in hemolymph were observed from April to August in the 5ng/L treatment during the 180 d exposure. Results indicate that the 180 d test approach, concurrent with the sequence of 28 d tests, enabled a more robust evaluation of mussel behavior and physiology than would have been possible with a single short-term (28 d) test.

Evaluating the Impact of Uncertainties in Clearance and Exposure When Prioritizing Chemicals Screened in High-Throughput Assays.

The toxicity-testing paradigm has evolved to include high-throughput (HT) methods for addressing the increasing need to screen hundreds to thousands of chemicals rapidly. Approaches that involve in vitro screening assays, in silico predictions of exposure concentrations, and pharmacokinetic (PK) characteristics provide the foundation for HT risk prioritization. Underlying uncertainties in predicted exposure concentrations or PK behaviors can significantly influence the prioritization of chemicals, though the impact of such influences is unclear. In the current study, a framework was developed to incorporate absorbed doses, PK properties, and in vitro dose-response data into a PK/pharmacodynamic (PD) model to allow for placement of chemicals into discrete priority bins. Literature-reported or predicted values for clearance rates and absorbed doses were used in the PK/PD model to evaluate the impact of their uncertainties on chemical prioritization. Scenarios using predicted absorbed doses resulted in a larger number of bin misassignments than those scenarios using predicted clearance rates, when comparing to bin placement using literature-reported values. Sensitivity of parameters on the model output of toxicological activity was examined across possible ranges for those parameters to provide insight into how uncertainty in their predicted values might impact uncertainty in activity.

Estimating Margin of Exposure to Thyroid Peroxidase Inhibitors Using High-Throughput in vitro Data, High-Throughput Exposure Modeling, and Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling.

Some pharmaceuticals and environmental chemicals bind the thyroid peroxidase (TPO) enzyme and disrupt thyroid hormone production. The potential for TPO inhibition is a function of both the binding affinity and concentration of the chemical within the thyroid gland. The former can be determined through in vitro assays, and the latter is influenced by pharmacokinetic properties, along with environmental exposure levels. In this study, a physiologically based pharmacokinetic (PBPK) model was integrated with a pharmacodynamic (PD) model to establish internal doses capable of inhibiting TPO in relation to external exposure levels predicted through exposure modeling. The PBPK/PD model was evaluated using published serum or thyroid gland chemical concentrations or circulating thyroxine (T4) and triiodothyronine (T3) hormone levels measured in rats and humans. After evaluation, the model was used to estimate human equivalent intake doses resulting in reduction of T4 and T3 levels by 10% (ED10) for 6 chemicals of varying TPO-inhibiting potencies. These chemicals were methimazole, 6-propylthiouracil, resorcinol, benzophenone-2, 2-mercaptobenzothiazole, and triclosan. Margin of exposure values were estimated for these chemicals using the ED10 and predicted population exposure levels for females of child-bearing age. The modeling approach presented here revealed that examining hazard or exposure alone when prioritizing chemicals for risk assessment may be insufficient, and that consideration of pharmacokinetic properties is warranted. This approach also provides a mechanism for integrating in vitro data, pharmacokinetic properties, and exposure levels predicted through high-throughput means when interpreting adverse outcome pathways based on biological responses.

Developing a Physiologically-Based Pharmacokinetic Model Knowledgebase in Support of Provisional Model Construction.

Developing physiologically-based pharmacokinetic (PBPK) models for chemicals can be resource-intensive, as neither chemical-specific parameters nor in vivo pharmacokinetic data are easily available for model construction. Previously developed, well-parameterized, and thoroughly-vetted models can be a great resource for the construction of models pertaining to new chemicals. A PBPK knowledgebase was compiled and developed from existing PBPK-related articles and used to develop new models. From 2,039 PBPK-related articles published between 1977 and 2013, 307 unique chemicals were identified for use as the basis of our knowledgebase. Keywords related to species, gender, developmental stages, and organs were analyzed from the articles within the PBPK knowledgebase. A correlation matrix of the 307 chemicals in the PBPK knowledgebase was calculated based on pharmacokinetic-relevant molecular descriptors. Chemicals in the PBPK knowledgebase were ranked based on their correlation toward ethylbenzene and gefitinib. Next, multiple chemicals were selected to represent exact matches, close analogues, or non-analogues of the target case study chemicals. Parameters, equations, or experimental data relevant to existing models for these chemicals and their analogues were used to construct new models, and model predictions were compared to observed values. This compiled knowledgebase provides a chemical structure-based approach for identifying PBPK models relevant to other chemical entities. Using suitable correlation metrics, we demonstrated that models of chemical analogues in the PBPK knowledgebase can guide the construction of PBPK models for other chemicals.

Completing the Link between Exposure Science and Toxicology for Improved Environmental Health Decision Making: The Aggregate Exposure Pathway Framework.

Driven by major scientific advances in analytical methods, biomonitoring, computation, and a newly articulated vision for a greater impact in public health, the field of exposure science is undergoing a rapid transition from a field of observation to a field of prediction. Deployment of an organizational and predictive framework for exposure science analogous to the "systems approaches" used in the biological sciences is a necessary step in this evolution. Here we propose the aggregate exposure pathway (AEP) concept as the natural and complementary companion in the exposure sciences to the adverse outcome pathway (AOP) concept in the toxicological sciences. Aggregate exposure pathways offer an intuitive framework to organize exposure data within individual units of prediction common to the field, setting the stage for exposure forecasting. Looking farther ahead, we envision direct linkages between aggregate exposure pathways and adverse outcome pathways, completing the source to outcome continuum for more meaningful integration of exposure assessment and hazard identification. Together, the two frameworks form and inform a decision-making framework with the flexibility for risk-based, hazard-based, or exposure-based decision making.

A Workflow to Investigate Exposure and Pharmacokinetic Influences on High-Throughput in Vitro Chemical Screening Based on Adverse Outcome Pathways.

BACKGROUND: Adverse outcome pathways (AOPs) link adverse effects in individuals or populations to a molecular initiating event (MIE) that can be quantified using in vitro methods. Practical application of AOPs in chemical-specific risk assessment requires incorporation of knowledge on exposure, along with absorption, distribution, metabolism, and excretion (ADME) properties of chemicals. OBJECTIVES: We developed a conceptual workflow to examine exposure and ADME properties in relation to an MIE. The utility of this workflow was evaluated using a previously established AOP, acetylcholinesterase (AChE) inhibition. METHODS: Thirty chemicals found to inhibit human AChE in the ToxCast™ assay were examined with respect to their exposure, absorption potential, and ability to cross the blood-brain barrier (BBB). Structures of active chemicals were compared against structures of 1,029 inactive chemicals to detect possible parent compounds that might have active metabolites. RESULTS: Application of the workflow screened 10 "low-priority" chemicals of 30 active chemicals. Fifty-two of the 1,029 inactive chemicals exhibited a similarity threshold of ≥ 75% with their nearest active neighbors. Of these 52 compounds, 30 were excluded due to poor absorption or distribution. The remaining 22 compounds may inhibit AChE in vivo either directly or as a result of metabolic activation. CONCLUSIONS: The incorporation of exposure and ADME properties into the conceptual workflow eliminated 10 "low-priority" chemicals that may otherwise have undergone additional, resource-consuming analyses. Our workflow also increased confidence in interpretation of in vitro results by identifying possible "false negatives." CITATION: Phillips MB, Leonard JA, Grulke CM, Chang DT, Edwards SW, Brooks R, Goldsmith MR, El-Masri H, Tan YM. 2016. A workflow to investigate exposure and pharmacokinetic influences on high-throughput in vitro chemical screening based on adverse outcome pathways. Environ Health Perspect 124:53-60; http://dx.doi.org/10.1289/ehp.1409450.

Metabolomic, behavioral, and reproductive effects of the aromatase inhibitor fadrozole hydrochloride on the unionid mussel Lampsilis fasciola.

Androgen-induced masculinization of female aquatic biota poses concerns for natural population stability. This research evaluated the effects of a twelve day exposure of fadrozole hydrochloride on the metabolism and reproductive status of the unionid mussel Lampsilis fasciola. Although this compound is not considered to be widespread in the aquatic environment, it was selected as a model aromatase (enzyme that converts testosterone to estradiol) inhibitor. Adult mussels were exposed to a control and 3 concentrations of fadrozole (2µg/L, 20µg/L, and 50µg/L), and samples of gill tissue were taken on days 4 and 12 for metabolomics analysis. Gills were used because of the variety of critical processes they mediate, such as feeding, ion exchange, and siphoning. Daily observed mussel behavior included female mantle display, foot protrusion, siphoning, and larval (glochidia) releases. Glochidia mortality was significantly higher in the 20µg/L treatment. Fewer conglutinate (packets of glochidia) releases were observed in the 50µg/L treatment, and mortality was highly correlated to release numbers. Foot protrusion was significantly higher in females in nearly all treatments, including the control, during the first 4days of observations. However, this sex difference was observed only in the 50µg/L treatment during the last 8days. Generally, metabolites were significantly altered in female gill tissue in the 2µg/L treatment whereas males were mostly affected only at the highest (50µg/L) treatment. Both sexes also revealed significant reductions in fadrozole-induced metabolic effects in gill tissue sampled after 12days compared to tissue sampled after 4days, indicating time-dependent mechanisms of disruptions in metabolic pathways and homeostatic processes to compensate for such disruptions.

Metabolomic, behavioral, and reproductive effects of the synthetic estrogen 17 α-ethinylestradiol on the unionid mussel Lampsilis fasciola.

The endocrine disrupting effects of estrogenic compounds in surface waters on fish, such as feminization of males and altered sex ratios, may also occur in aquatic invertebrates. However, the underlying mechanisms of action and toxicity, especially in native freshwater mussels (Order Unionoida), remain undefined. This study evaluated the effects of a 12-day exposure of 17 α-ethinylestradiol (EE2), a synthetic estrogen in oral contraceptives commonly found in surface waters, on the behavior, condition, metabolism, and reproductive status of Lampsilis fasciola. Adult mussels of both sexes were exposed to a control and two concentrations of EE2 (0 ng/L, 5 ng/L considered to be environmentally relevant, and 1,000 ng/L designed to provide a positive metabolic response), and samples of gill tissue were taken on days 4 and 12; gills were used because of the variety of critical processes they mediate, such as feeding, ion exchange, and siphoning. Observations of mussel behavior (mantle display, siphoning, and foot movement) were made daily, and condition of conglutinates (packets of eggs and/or glochidia) released by females was examined. No significant effects of EE2 on glochidia mortality, conglutinate condition, female marsupial gill condition, or mussel foot extension were observed. However, exposure to both concentrations of EE2 significantly reduced male siphoning and mantle display behavior of females. Metabolomics analyses identified 207 known biochemicals in mussel gill tissue and showed that environmentally relevant EE2 concentrations led to decreases in glycogen metabolism end products, glucose, and several essential fatty acids in females after 12 days, indicating reductions in energy reserves that could otherwise be used for growth or reproduction. Moreover, males and females showed significant alterations in metabolites involved in signal transduction, immune response, and neuromodulation. Most of these changes were apparent at 1,000 ng/L EE2, but similar metabolites and pathways were also affected at 5 ng/L EE2. Components of the extracellular matrix of gill tissue were also altered. These results demonstrate the utility of metabolomics when used in conjunction with traditional physiological and behavioral toxicity test endpoints and establish the usefulness of this approach in determining possible underlying toxicological mechanisms of EE2 in exposed freshwater mussels.