Bioactivity Signatures of Drugs vs. Environmental Chemicals Revealed by Tox21 High-Throughput Screening Assays

Screening and characterization of 133 physiologically-relevant environmental chemicals for reproductive toxicity

Reproduction is a functional outcome that relies on complex cellular, tissue, and organ interactions that span the developmental period to adulthood. Thus, the assessment of its disruption by environmental chemicals would benefit significantly from scalable and innovative approaches to testing using functionally comparable reproductive models such as the nematode C. elegans. We adapted a previously described low-throughput in vivo chromosome segregation assay using C. elegans predictive of reproductive toxicity and leveraged available public data sources (ToxCast, ICE) to screen and characterize 133 physiologically-relevant chemicals in a high-throughput manner. The screening outcome was further validated in a second, independent in vivo assay assessing embryonic viability. In total, 13 chemicals were classified as reproductive toxicants with the two most active chemicals belonging to the large family of Quaternary Ammonium Compounds (QACs) commonly used as disinfectants but with limited available reproductive toxicity data. We compared the results from the C. elegans assay with ToxCast in vitro data compiled from 700+ cell response assays and 300+ signaling pathways-based assays. We did not observe a difference in the bioactivity or in the average potency (AC50) between the top and bottom chemicals. However, the intended target categories were significantly different between the classified chemicals with, in particular, an over-representation of steroid hormone targets for the high Z-score chemicals. Taken together, these results point to the value of in vivo models that scale to high-throughput level for reproductive toxicity assessment and to the need to prioritize the assessment of QACs impacts on reproduction.

Screening and characterization of 133 physiologically-relevant environmental chemicals for reproductive toxicity

Reproduction is a functional outcome that relies on complex cellular, tissue, and organ interactions that span the developmental period to adulthood. Thus, the assessment of its disruption by environmental chemicals is remarkably painstaking in conventional toxicological animal models and does not scale up to the number of chemicals present in our environment and requiring testing. We adapted a previously described low-throughput in vivo chromosome segregation assay using C. elegans predictive of reproductive toxicity and leveraged available public data sources (ToxCast, ICE) to screen and characterize 133 physiologically-relevant chemicals in a high-throughput manner. The screening outcome was further validated in a second, independent in vivo assay assessing embryonic viability. In total, 13 chemicals were classified as reproductive toxicants with the two most active chemicals belonging to the large family of Quaternary Ammonium Compounds (QACs) commonly used as disinfectants but with limited available reproductive toxicity data. We compared the results from the C. elegans assay with ToxCast in vitro data compiled from 700+ cell response assays and 300+ signaling pathways-based assays. We did not observe a difference in the bioactivity or in average potency (AC50) between the top and bottom chemicals. However, the intended target categories were significantly different between the classified chemicals with, in particular, an over-representation of steroid hormone targets for the high Z-score chemicals. Taken together, these results point to the value of in vivo models that scale to high-throughput level for reproductive toxicity assessment and to the need to prioritize the assessment of QACs impacts on reproduction.

In vitro profiling of pesticides within the Tox21 10K compound library for bioactivity and potential toxicity

Pesticides include a diverse class of toxic chemicals, often having numerous modes of actions when used in agriculture against targeted organisms to control insect infestation, halt unwanted vegetation, and prevent the spread of disease. In this study, the in vitro assay activity of pesticides within the Tox21 10K compound library were examined. The assays in which pesticides showed significantly more activities than non-pesticide chemicals revealed potential targets and mechanisms of action for pesticides. Furthermore, pesticides that showed promiscuous activity against many targets and cytotoxicity were identified, which warrant further toxicological evaluation. Several pesticides were shown to require metabolic activation, demonstrating the importance of introducing metabolic capacity to in vitro assays. Overall, the activity profiles of pesticides highlighted in this study can contribute to the knowledge gaps surrounding pesticide mechanisms and to the better understanding of the on- and off-target organismal effects of pesticides.

In vivo quantitative high-throughput screening for drug discovery and comparative toxicology

Quantitative high-throughput screening (qHTS) pharmacologically evaluates chemical libraries for therapeutic uses, toxicological risk and, increasingly, for academic probe discovery. Phenotypic high-throughput screening assays interrogate molecular pathways, often relying on cell culture systems, historically less focused on multicellular organisms. Caenorhabditis elegans has served as a eukaryotic model organism for human biology by virtue of genetic conservation and experimental tractability. Here, a paradigm enabling C. elegans qHTS using 384-well microtiter plate laser-scanning cytometry is described, in which GFP-expressing organisms revealing phenotype-modifying structure-activity relationships guide subsequent life-stage and proteomic analyses, and Escherichia coli bacterial ghosts, a non-replicating nutrient source, allow compound exposures over two life cycles, mitigating bacterial overgrowth complications. We demonstrate the method with libraries of anti-infective agents, or substances of toxicological concern. Each was tested in seven-point titration to assess the feasibility of nematode-based in vivo qHTS, and examples of follow-up strategies were provided to study organism-based chemotype selectivity and subsequent network perturbations with a physiological impact. We anticipate that this qHTS approach will enable analysis of C. elegans orthologous phenotypes of human pathologies to facilitate drug library profiling for a range of therapeutic indications.

Repurposing drugs as COVID-19 therapies: a toxicity evaluation

Drug repurposing is an appealing method to address the Coronavirus 2019 (COVID-19) pandemic because of the low cost and efficiency. We analyzed our in-house database of approved drug screens and compared their activity profiles with results from a severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) cytopathic effect (CPE) assay. The activity profiles of the human ether-à-go-go-related gene (hERG), phospholipidosis (PLD), and many cytotoxicity screens were found significantly correlated with anti-SARS-CoV-2 activity. hERG inhibition is a nonspecific off-target effect that has contributed to promiscuous drug interactions, whereas drug-induced PLD is an undesirable effect linked to hERG blockers. Thus, this study identifies preferred drug candidates as well as chemical structures that should be avoided because of their potential to induce toxicity. Lastly, we highlight the hERG liability of anti-SARS-CoV-2 drugs currently enrolled in clinical trials.

Prediction Models for Agonists and Antagonists of Molecular Initiation Events for Toxicity Pathways Using an Improved Deep-Learning-Based Quantitative Structure-Activity Relationship System

In silico approaches have been studied intensively to assess the toxicological risk of various chemical compounds as alternatives to traditional in vivo animal tests. Among these approaches, quantitative structure-activity relationship (QSAR) analysis has the advantages that it is able to construct models to predict the biological properties of chemicals based on structural information. Previously, we reported a deep learning (DL) algorithm-based QSAR approach called DeepSnap-DL for high-performance prediction modeling of the agonist and antagonist activity of key molecules in molecular initiating events in toxicological pathways using optimized hyperparameters. In the present study, to achieve high throughput in the DeepSnap-DL system-which consists of the preparation of three-dimensional molecular structures of chemical compounds, the generation of snapshot images from the three-dimensional chemical structures, DL, and statistical calculations-we propose an improved DeepSnap-DL approach. Using this improved system, we constructed 59 prediction models for the agonist and antagonist activity of key molecules in the Tox21 10K library. The results indicate that modeling of the agonist and antagonist activity with high prediction performance and high throughput can be achieved by optimizing suitable parameters in the improved DeepSnap-DL system.

The Tox21 10K Compound Library: Collaborative Chemistry Advancing Toxicology

Since 2009, the Tox21 project has screened ∼8500 chemicals in more than 70 high-throughput assays, generating upward of 100 million data points, with all data publicly available through partner websites at the United States Environmental Protection Agency (EPA), National Center for Advancing Translational Sciences (NCATS), and National Toxicology Program (NTP). Underpinning this public effort is the largest compound library ever constructed specifically for improving understanding of the chemical basis of toxicity across research and regulatory domains. Each Tox21 federal partner brought specialized resources and capabilities to the partnership, including three approximately equal-sized compound libraries. All Tox21 data generated to date have resulted from a confluence of ideas, technologies, and expertise used to design, screen, and analyze the Tox21 10K library. The different programmatic objectives of the partners led to three distinct, overlapping compound libraries that, when combined, not only covered a diversity of chemical structures, use-categories, and properties but also incorporated many types of compound replicates. The history of development of the Tox21 "10K" chemical library and data workflows implemented to ensure quality chemical annotations and allow for various reproducibility assessments are described. Cheminformatics profiling demonstrates how the three partner libraries complement one another to expand the reach of each individual library, as reflected in coverage of regulatory lists, predicted toxicity end points, and physicochemical properties. ToxPrint chemotypes (CTs) and enrichment approaches further demonstrate how the combined partner libraries amplify structure-activity patterns that would otherwise not be detected. Finally, CT enrichments are used to probe global patterns of activity in combined ToxCast and Tox21 activity data sets relative to test-set size and chemical versus biological end point diversity, illustrating the power of CT approaches to discern patterns in chemical-activity data sets. These results support a central premise of the Tox21 program: A collaborative merging of programmatically distinct compound libraries would yield greater rewards than could be achieved separately.

The Tox21 10K Compound Library: Collaborative Chemistry Advancing Toxicology

Since 2009, the Tox21 project has screened ∼8500 chemicals in more than 70 high-throughput assays, generating upward of 100 million data points, with all data publicly available through partner websites at the United States Environmental Protection Agency (EPA), National Center for Advancing Translational Sciences (NCATS), and National Toxicology Program (NTP). Underpinning this public effort is the largest compound library ever constructed specifically for improving understanding of the chemical basis of toxicity across research and regulatory domains. Each Tox21 federal partner brought specialized resources and capabilities to the partnership, including three approximately equal-sized compound libraries. All Tox21 data generated to date have resulted from a confluence of ideas, technologies, and expertise used to design, screen, and analyze the Tox21 10K library. The different programmatic objectives of the partners led to three distinct, overlapping compound libraries that, when combined, not only covered a diversity of chemical structures, use-categories, and properties but also incorporated many types of compound replicates. The history of development of the Tox21 "10K" chemical library and data workflows implemented to ensure quality chemical annotations and allow for various reproducibility assessments are described. Cheminformatics profiling demonstrates how the three partner libraries complement one another to expand the reach of each individual library, as reflected in coverage of regulatory lists, predicted toxicity end points, and physicochemical properties. ToxPrint chemotypes (CTs) and enrichment approaches further demonstrate how the combined partner libraries amplify structure-activity patterns that would otherwise not be detected. Finally, CT enrichments are used to probe global patterns of activity in combined ToxCast and Tox21 activity data sets relative to test-set size and chemical versus biological end point diversity, illustrating the power of CT approaches to discern patterns in chemical-activity data sets. These results support a central premise of the Tox21 program: A collaborative merging of programmatically distinct compound libraries would yield greater rewards than could be achieved separately.

Tissue distribution, metabolism, and elimination of Victoria Pure Blue BO in rainbow trout: Main metabolite as an appropriate residue marker

Victoria Pure Blue BO is a dye that bears some therapeutic activity and that can be retrieved in effluent or may be used in aquaculture as a prohibited drug. In this study, the metabolism and tissue distribution during uptake and depuration of VPBO were investigated in order to propose a residue marker of illegal treatment in fish. The dye was administered to rainbow trout (oncorhynchus mykiss) for one day by water bath at a dose of 0.1 mg.L^-1. The concentrations of VPBO in all tissues increased rapidly during the treatment period, reaching a C_max of 567 ± 301 μg.L^-1 in plasma and 1846 μg kg^-1 ±517 for liver after 2 h. After placing the rainbow trout in a clean water bath for a 64 day-period of depuration, the concentrations in the tissues and plasma decreased to reach comparable levels for muscle and for skin after 33 days. The concentrations measured were still above the LOQ at 2.26 ± 0.48 μg kg^-1 for muscle and 2.85 ± 1.99 μg kg^-1 for skin at the end of the depuration period. The results indicated the existence of 14 phase I metabolites and one glucuronide conjugated metabolite. Non-compartmental analysis was applied to assess the pharmacokinetic parameters. The half-life in edible muscle of the main metabolite detected, deethyl-leuco-VPBO, was found to be 22.5 days compared to a half-life of 19.7 days for the parent VPBO. This study provides new information to predict a VPBO drug treatment of aquacultured species via a proposed new residue marker.