Nisha S. Sipes

Zebrafish developmental screening of the ToxCast™ Phase I chemical library

Zebrafish (Danio rerio) is an emerging toxicity screening model for both human health and ecology. As part of the Computational Toxicology Research Program of the U.S. EPA, the toxicity of the 309 ToxCast™ Phase I chemicals was assessed using a zebrafish screen for developmental toxicity. All exposures were by immersion from 6-8 h post fertilization (hpf) to 5 days post fertilization (dpf); nominal concentration range of 1 nM-80 μM. On 6 dpf larvae were assessed for death and overt structural defects. Results revealed that the majority (62%) of chemicals were toxic to the developing zebrafish; both toxicity incidence and potency was correlated with chemical class and hydrophobicity (logP); and inter-and intra-plate replicates showed good agreement. The zebrafish embryo screen, by providing an integrated model of the developing vertebrate, compliments the ToxCast assay portfolio and has the potential to provide information relative to overt and organismal toxicity.

Zebrafish: as an integrative model for twenty-first century toxicity testing.

The zebrafish embryo is a useful small model for investigating vertebrate development because of its transparency, low cost, transgenic and morpholino capabilities, conservation of cell signaling, and concordance with mammalian developmental phenotypes. From these advantages, the zebrafish embryo has been considered as an alternative model for traditional in vivo developmental toxicity screening. The use of this organism in conjunction with traditional in vivo developmental toxicity testing has the potential to reduce cost and increase throughput of testing the chemical universe, prioritize chemicals for targeted toxicity testing, generate predictive models of developmental toxicants, and elucidate mechanisms and adverse outcome pathways for abnormal development. This review gives an overview of the zebrafish embryo for pre dictive toxicology and 21st century toxicity testing. Developmental eye defects were selected as an example to evaluate data from the U.S. Environmental Protection Agencys ToxCast program comparing responses in zebrafish embryos with those from pregnant rats and rabbits for a subset of 24 environmental chemicals across >600 in vitro assay targets. Cross-species comparisons implied a common basis for biological pathways associated with neuronal defects, extracellular matrix remodeling, and mitotic arrest.

Predictive Models of Prenatal Developmental Toxicity from ToxCast High-Throughput Screening Data

Environmental Protection Agencys ToxCast project is profiling the in vitro bioactivity of chemicals to assess pathway-level and cell-based signatures that correlate with observed in vivo toxicity. We hypothesized that developmental toxicity in guideline animal studies captured in the ToxRefDB database would correlate with cell-based and cell-free in vitro high-throughput screening (HTS) data to reveal meaningful mechanistic relationships and provide models identifying chemicals with the potential to cause developmental toxicity. To test this hypothesis, we built statistical associations based on HTS and in vivo developmental toxicity data from ToxRefDB. Univariate associations were used to filter HTS assays based on statistical correlation with distinct in vivo endpoint. This revealed 423 total associations with distinctly different patterns for rat (301 associations) and rabbit (122 associations) across multiple HTS assay platforms. From these associations, linear discriminant analysis with cross-validation was used to build the models. Species-specific models of predicted developmental toxicity revealed strong balanced accuracy (> 70%) and unique correlations between assay targets such as transforming growth factor beta, retinoic acid receptor, and G-protein-coupled receptor signaling in the rat and inflammatory signals, such as interleukins (IL) (IL1a and IL8) and chemokines (CCL2), in the rabbit. Species-specific toxicity endpoints were associated with one another through common Gene Ontology biological processes, such as cleft palate to urogenital defects through placenta and embryonic development. This work indicates the utility of HTS assays for developing pathway-level models predictive of developmental toxicity.

Activity profiles of 309 ToxCast™ chemicals evaluated across 292 biochemical targets☆

Understanding the potential health risks posed by environmental chemicals is a significant challenge elevated by the large number of diverse chemicals with generally uncharacterized exposures, mechanisms, and toxicities. The present study is a performance evaluation and critical analysis of assay results for an array of 292 high-throughput cell-free assays aimed at preliminary toxicity evaluation of 320 environmental chemicals in EPAs ToxCast™ project (Phase I). The chemicals (309 unique, 11 replicates) were mainly precursors or the active agent of commercial pesticides, for which a wealth of in vivo toxicity data is available. Biochemical HTS (high-throughput screening) profiled cell and tissue extracts using semi-automated biochemical and pharmacological methodologies to evaluate a subset of G-protein coupled receptors (GPCRs), CYP450 enzymes (CYPs), kinases, phosphatases, proteases, HDACs, nuclear receptors, ion channels, and transporters. The primary screen tested all chemicals at a relatively high concentration 25 μM concentration (or 10 μM for CYP assays), and a secondary screen re-tested 9132 chemical-assay pairs in 8-point concentration series from 0.023 to 50 μM (or 0.009-20 μM for CYPs). Mapping relationships across 93,440 chemical-assay pairs based on half-maximal activity concentration (AC50) revealed both known and novel targets in signaling and metabolic pathways. The primary dataset, summary data and details on quality control checks are available for download at http://www.epa.gov/ncct/toxcast/.

Profiling 976 ToxCast chemicals across 331 enzymatic and receptor signaling assays.

Understanding potential health risks is a significant challenge due to the large numbers of diverse chemicals with poorly characterized exposures and mechanisms of toxicities. The present study analyzes 976 chemicals (including failed pharmaceuticals, alternative plasticizers, food additives, and pesticides) in Phases I and II of the U.S. EPA’s ToxCast project across 331 cell-free enzymatic and ligand-binding high-throughput screening (HTS) assays. Half-maximal activity concentrations (AC50) were identified for 729 chemicals in 256 assays (7,135 chemical–assay pairs). Some of the most commonly affected assays were CYPs (CYP2C9 and CYP2C19), transporters (mitochondrial TSPO, norepinephrine, and dopaminergic), and GPCRs (aminergic). Heavy metals, surfactants, and dithiocarbamate fungicides showed promiscuous but distinctly different patterns of activity, whereas many of the pharmaceutical compounds showed promiscuous activity across GPCRs. Literature analysis confirmed >50% of the activities for the most potent chemical–assay pairs (54) but also revealed 10 missed interactions. Twenty-two chemicals with known estrogenic activity were correctly identified for the majority (77%), missing only the weaker interactions. In many cases, novel findings for previously unreported chemical–target combinations clustered with known chemical–target interactions. Results from this large inventory of chemical–biological interactions can inform read-across methods as well as link potential targets to molecular initiating events in adverse outcome pathways for diverse toxicities.

Environmental Impact on Vascular Development Predicted by High-Throughput Screening

Background: Understanding health risks to embryonic development from exposure to environmental chemicals is a significant challenge given the diverse chemical landscape and paucity of data for most of these compounds. High-throughput screening (HTS) in the U.S. Environmental Protection Agency (EPA) ToxCast™ project provides vast data on an expanding chemical library currently consisting of > 1,000 unique compounds across > 500 in vitro assays in phase I (complete) and Phase II (under way). This public data set can be used to evaluate concentration-dependent effects on many diverse biological targets and build predictive models of prototypical toxicity pathways that can aid decision making for assessments of human developmental health and disease. Objective: We mined the ToxCast phase I data set to identify signatures for potential chemical disruption of blood vessel formation and remodeling. Methods: ToxCast phase I screened 309 chemicals using 467 HTS assays across nine assay technology platforms. The assays measured direct interactions between chemicals and molecular targets (receptors, enzymes), as well as downstream effects on reporter gene activity or cellular consequences. We ranked the chemicals according to individual vascular bioactivity score and visualized the ranking using ToxPi (Toxicological Priority Index) profiles. Results: Targets in inflammatory chemokine signaling, the vascular endothelial growth factor pathway, and the plasminogen-activating system were strongly perturbed by some chemicals, and we found positive correlations with developmental effects from the U.S. EPA ToxRefDB (Toxicological Reference Database) in vivo database containing prenatal rat and rabbit guideline studies. We observed distinctly different correlative patterns for chemicals with effects in rabbits versus rats, despite derivation of in vitro signatures based on human cells and cell-free biochemical targets, implying conservation but potentially differential contributions of developmental pathways among species. Follow-up analysis with antiangiogenic thalidomide analogs and additional in vitro vascular targets showed in vitro activity consistent with the most active environmental chemicals tested here. Conclusions: We predicted that blood vessel development is a target for environmental chemicals acting as putative vascular disruptor compounds (pVDCs) and identified potential species differences in sensitive vascular developmental pathways.

A Computational Model Predicting Disruption of Blood Vessel Development

Vascular development is a complex process regulated by dynamic biological networks that vary in topology and state across different tissues and developmental stages. Signals regulating de novo blood vessel formation (vasculogenesis) and remodeling (angiogenesis) come from a variety of biological pathways linked to endothelial cell (EC) behavior, extracellular matrix (ECM) remodeling and the local generation of chemokines and growth factors. Simulating these interactions at a systems level requires sufficient biological detail about the relevant molecular pathways and associated cellular behaviors, and tractable computational models that offset mathematical and biological complexity. Here, we describe a novel multicellular agent-based model of vasculogenesis using the CompuCell3D (http://www.compucell3d.org/) modeling environment supplemented with semi-automatic knowledgebase creation. The model incorporates vascular endothelial growth factor signals, pro- and anti-angiogenic inflammatory chemokine signals, and the plasminogen activating system of enzymes and proteases linked to ECM interactions, to simulate nascent EC organization, growth and remodeling. The model was shown to recapitulate stereotypical capillary plexus formation and structural emergence of non-coded cellular behaviors, such as a heterologous bridging phenomenon linking endothelial tip cells together during formation of polygonal endothelial cords. Molecular targets in the computational model were mapped to signatures of vascular disruption derived from in vitro chemical profiling using the EPAs ToxCast high-throughput screening (HTS) dataset. Simulating the HTS data with the cell-agent based model of vascular development predicted adverse effects of a reference anti-angiogenic thalidomide analog, 5HPP-33, on in vitro angiogenesis with respect to both concentration-response and morphological consequences. These findings support the utility of cell agent-based models for simulating a morphogenetic series of events and for the first time demonstrate the applicability of these models for predictive toxicology.

Evaluation of 309 Environmental Chemicals Using a Mouse Embryonic Stem Cell Adherent Cell Differentiation and Cytotoxicity Assay

The vast landscape of environmental chemicals has motivated the need for alternative methods to traditional whole-animal bioassays in toxicity testing. Embryonic stem (ES) cells provide an in vitro model of embryonic development and an alternative method for assessing developmental toxicity. Here, we evaluated 309 environmental chemicals, mostly food-use pesticides, from the ToxCast™ chemical library using a mouse ES cell platform. ES cells were cultured in the absence of pluripotency factors to promote spontaneous differentiation and in the presence of DMSO-solubilized chemicals at different concentrations to test the effects of exposure on differentiation and cytotoxicity. Cardiomyocyte differentiation (α,β myosin heavy chain; MYH6/MYH7) and cytotoxicity (DRAQ5™/Sapphire700™) were measured by In-Cell Western™ analysis. Half-maximal activity concentration (AC50) values for differentiation and cytotoxicity endpoints were determined, with 18% of the chemical library showing significant activity on either endpoint. Mining these effects against the ToxCast Phase I assays (∼500) revealed significant associations for a subset of chemicals (26) that perturbed transcription-based activities and impaired ES cell differentiation. Increased transcriptional activity of several critical developmental genes including BMPR2, PAX6 and OCT1 were strongly associated with decreased ES cell differentiation. Multiple genes involved in reactive oxygen species signaling pathways (NRF2, ABCG2, GSTA2, HIF1A) were strongly associated with decreased ES cell differentiation as well. A multivariate model built from these data revealed alterations in ABCG2 transporter was a strong predictor of impaired ES cell differentiation. Taken together, these results provide an initial characterization of metabolic and regulatory pathways by which some environmental chemicals may act to disrupt ES cell growth and differentiation.

Predictive Endocrine Testing in the 21st Century Using in Vitro Assays of Estrogen Receptor Signaling Responses

Thousands of environmental chemicals are subject to regulatory review for their potential to be endocrine disruptors (ED). In vitro high-throughput screening (HTS) assays have emerged as a potential tool for prioritizing chemicals for ED-related whole-animal tests. In this study, 1814 chemicals including pesticide active and inert ingredients, industrial chemicals, food additives, and pharmaceuticals were evaluated in a panel of 13 in vitro HTS assays. The panel of in vitro assays interrogated multiple end points related to estrogen receptor (ER) signaling, namely binding, agonist, antagonist, and cell growth responses. The results from the in vitro assays were used to create an ER Interaction Score. For 36 reference chemicals, an ER Interaction Score >0 showed 100% sensitivity and 87.5% specificity for classifying potential ER activity. The magnitude of the ER Interaction Score was significantly related to the potency classification of the reference chemicals (p < 0.0001). ERα/ERβ selectivity was also evaluated, but relatively few chemicals showed significant selectivity for a specific isoform. When applied to a broader set of chemicals with in vivo uterotrophic data, the ER Interaction Scores showed 91% sensitivity and 65% specificity. Overall, this study provides a novel method for combining in vitro concentration response data from multiple assays and, when applied to a large set of ER data, accurately predicted estrogenic responses and demonstrated its utility for chemical prioritization.

Editor's Highlight: Analysis of the Effects of Cell Stress and Cytotoxicity on In Vitro Assay Activity Across a Diverse Chemical and Assay Space.

Chemical toxicity can arise from disruption of specific biomolecular functions or through more generalized cell stress and cytotoxicity-mediated processes. Here, responses of 1060 chemicals including pharmaceuticals, natural products, pesticidals, consumer, and industrial chemicals across a battery of 815 in vitro assay endpoints from 7 high-throughput assay technology platforms were analyzed in order to distinguish between these types of activities. Both cell-based and cell-free assays showed a rapid increase in the frequency of responses at concentrations where cell stress/cytotoxicity responses were observed in cell-based assays. Chemicals that were positive on at least 2 viability/cytotoxicity assays within the concentration range tested (typically up to 100 μM) activated a median of 12% of assay endpoints whereas those that were not cytotoxic in this concentration range activated 1.3% of the assays endpoints. The results suggest that activity can be broadly divided into: (1) specific biomolecular interactions against one or more targets (eg, receptors or enzymes) at concentrations below which overt cytotoxicity-associated activity is observed; and (2) activity associated with cell stress or cytotoxicity, which may result from triggering specific cell stress pathways, chemical reactivity, physico-chemical disruption of proteins or membranes, or broad low-affinity non-covalent interactions. Chemicals showing a greater number of specific biomolecular interactions are generally designed to be bioactive (pharmaceuticals or pesticidal active ingredients), whereas intentional food-use chemicals tended to show the fewest specific interactions. The analyses presented here provide context for use of these data in ongoing studies to predict in vivo toxicity from chemicals lacking extensive hazard assessment.