Identification of a dual FLT3 and MNK2 inhibitor for acute myeloid leukemia treatment using a structure-based virtual screening approach.

Fms-like tyrosine kinase 3 (FLT3) is considered a promising therapeutic target for acute myeloid leukemia (AML) in the clinical. However, monotherapy with FLT3 inhibitor is usually accompanied by drug resistance. Dual inhibitors might be therapeutically beneficial to patients with AML due to their ability to overcome drug resistance. Mitogen-activated protein kinase (MAPK)-interacting kinases (MNKs) phosphorylate eukaryotic translation initiation factor 4E (eIF4E), which brings together the RAS/RAF/ERK and PI3K/AKT/mTOR oncogenic pathways. Therefore, dual inhibition of FLT3 and MNK2 might have an additive effect against AML. Herein, a structure-based virtual screening approach was performed to identify dual inhibitors of FLT3 and MNK2 from the ChemDiv database. Compound K783-0308 was identified as a dual inhibitor of FLT3 and MNK2 with IC50 values of 680 and 406 nM, respectively. In addition, the compound showed selectivity for both FLT3 and MNK2 in a panel of 82 kinases. The structure-activity relationship analysis and common interactions revealed interactions between K783-0308 analogs and FLT3 and MNK2. Furthermore, K783-0308 inhibited MV-4-11 and MOLM-13 AML cell growth and induced G0/G1 cell cycle arrest. Taken together, the dual inhibitor K783-0308 showed promising results and can be potentially optimized as a lead compound for AML treatment.

Structure-based virtual screening and biological evaluation of novel small-molecule BTK inhibitors.

Bruton tyrosine kinase (BTK) is linked to multiple signalling pathways that regulate cellular survival, activation, and proliferation. A covalent BTK inhibitor has shown favourable outcomes for treating B cell malignant leukaemia. However, covalent inhibitors require a high reactive warhead that may contribute to unexpected toxicity, poor selectivity, or reduced effectiveness in solid tumours. Herein, we report the identification of a novel noncovalent BTK inhibitor. The binding interactions (i.e. interactions from known BTK inhibitors) for the BTK binding site were identified and incorporated into a structure-based virtual screening (SBVS). Top-rank compounds were selected and testing revealed a BTK inhibitor with >50% inhibition at 10 µM concentration. Examining analogues revealed further BTK inhibitors. When tested across solid tumour cell lines, one inhibitor showed favourable inhibitory activity, suggesting its potential for targeting BTK malignant tumours. This inhibitor could serve as a basis for developing an effective BTK inhibitor targeting solid cancers.

Harnessing In Silico, In Vitro, and In Vivo Data to Understand the Toxicity Landscape of Polycyclic Aromatic Compounds (PACs).

Polycyclic aromatic compounds (PACs) are compounds with a minimum of two six-atom aromatic fused rings. PACs arise from incomplete combustion or thermal decomposition of organic matter and are ubiquitous in the environment. Within PACs, carcinogenicity is generally regarded to be the most important public health concern. However, toxicity in other systems (reproductive and developmental toxicity, immunotoxicity) has also been reported. Despite the large number of PACs identified in the environment, research attention to understand exposure and health effects of PACs has focused on a relatively limited subset, namely polycyclic aromatic hydrocarbons (PAHs), the PACs with only carbon and hydrogen atoms. To triage the rest of the vast number of PACs for more resource-intensive testing, we developed a data-driven approach to contextualize hazard characterization of PACs, by leveraging the available data from various data streams (in silico toxicity, in vitro activity, structural fingerprints, and in vivo data availability). The PACs were clustered on the basis of their in silico toxicity profiles containing predictions from 8 different categories (carcinogenicity, cardiotoxicity, developmental toxicity, genotoxicity, hepatotoxicity, neurotoxicity, reproductive toxicity, and urinary toxicity). We found that PACs with the same parent structure (e.g., fluorene) could have diverse in silico toxicity profiles. In contrast, PACs with similar substituted groups (e.g., alkylated-PAHs) or heterocyclics (e.g., N-PACs) with varying ring sizes could have similar in silico toxicity profiles, suggesting that these groups are better candidates for toxicity read-across analysis. The clusters/regions associated with certain in silico toxicity, in vitro activity, and structural fingerprints were identified. We found that genotoxicity/carcinogenicity (in silico toxicity) and xenobiotic homeostasis and stress response (in vitro activity), respectively, dominate the toxicity/activity variation seen in the PACs. The "hot spots" with enriched toxicity/activity in conjunction with availability of in vivo carcinogenicity data revealed regions of either data-poor (hydroxylated-PAHs) or data-rich (unsubstituted, parent PAHs) PACs. These regions offer potential targets for prioritization of further in vivo assessment and for chemical read-across efforts. The analysis results are searchable through an interactive web application (https://ntp.niehs.nih.gov/go/pacs_tableau), allowing for alternative hypothesis generation.

A Gene Expression Biomarker Identifies Chemical Modulators of Estrogen Receptor α in an MCF-7 Microarray Compendium.

Identification of chemicals that affect hormone-regulated systems will help to predict endocrine disruption. In our previous study, a 46 gene biomarker was found to be an accurate predictor of estrogen receptor (ER) α modulation in chemically treated MCF-7 cells. Here, potential ERα modulators were identified using the biomarker by screening a microarray compendium consisting of ∼1600 gene expression comparisons representing exposure to ∼1200 chemicals. A total of ∼170 chemicals were identified as potential ERα modulators. In the Connectivity Map 2.0 collection, 75 and 39 chemicals were predicted to activate or suppress ERα, and they included 12 and six known ERα agonists and antagonists/selective ERα modulators, respectively. Nineteen and eight of the total number were also identified as active in an ERα transactivation assay carried out in an MCF-7-derived cell line used to screen the Tox21 10K chemical library in agonist or antagonist modes, respectively. Chemicals predicted to modulate ERα in MCF-7 cells were examined further using global and targeted gene expression in wild-type and ERα-null cells, transactivation assays, and cell-free ERα coregulator interaction assays. Environmental chemicals classified as weak and very weak agonists were confirmed to activate ERα including apigenin, kaempferol, and oxybenzone. Novel activators included digoxin, nabumetone, ivermectin, and six progestins. Novel suppressors included emetine, mifepristone, niclosamide, and proscillaridin. Our strategy will be useful to identify environmentally relevant ERα modulators in future high-throughput transcriptomic screens.

Pose Filter-Based Ensemble Learning Enables Discovery of Orally Active, Nonsteroidal Farnesoid X Receptor Agonists.

Farnesoid X receptor (FXR) agonists can reverse dysregulated bile acid metabolism, and thus, they are potential therapeutics to prevent and treat nonalcoholic fatty liver disease. The low success rate of FXR agonists' R&D and the side effects of clinical candidates such as obeticholic acid make it urgent to discover new chemotypes. Unfortunately, structure-based virtual screening (SBVS) that can speed up drug discovery has rarely been reported with success for FXR, which was likely hindered by the failure in addressing protein flexibility. To address this issue, we devised human FXR (hFXR)-specific ensemble learning models based on pose filters from 24 agonist-bound hFXR crystal structures and coupled them to traditional SBVS approaches of the FRED docking plus Chemgauss4 scoring function. It turned out that the hFXR-specific pose filter ensemble (PFE) was able to improve ligand enrichment significantly, which rendered 3RUT-based SBVS with its PFE the ideal approach for FXR agonist discovery. By screening of the Specs chemical library and in vitro FXR transactivation bioassay, we identified a new class of FXR agonists with compound XJ034 as the representative, which would have been missed if the PFE was not coupled. Following that, we performed in-depth biological studies which demonstrated that XJ034 resulted in a downtrend of intracellular triglyceride in vitro, significantly decreased the serum/liver TG in high fat diet-induced C57BL/6J obese mice, and more importantly, showed metabolic stabilities in both plasma and liver microsomes. To provide insight into further structure-based lead optimization, we solved the crystal structure of hFXR complexed with compound XJ034, uncovering a unique hydrogen bond between compound XJ034 and residue Y375. The current work highlights the power of our pose filter-based ensemble learning approach in terms of scaffold hopping and provides a promising lead compound for further development.

Identifying Compounds with Genotoxicity Potential Using Tox21 High-Throughput Screening Assays.

Genotoxicity is a critical component of a comprehensive toxicological profile. The Tox21 Program used five quantitative high-throughput screening (qHTS) assays measuring some aspect of DNA damage/repair to provide information on the genotoxic potential of over 10 000 compounds. Included were assays detecting activation of p53, increases in the DNA repair protein ATAD5, phosphorylation of H2AX, and enhanced cytotoxicity in DT40 cells deficient in DNA-repair proteins REV3 or KU70/RAD54. Each assay measures a distinct component of the DNA damage response signaling network; >70% of active compounds were detected in only one of the five assays. When qHTS results were compared with results from three standard genotoxicity assays (bacterial mutation, in vitro chromosomal aberration, and in vivo micronucleus), a maximum of 40% of known, direct-acting genotoxicants were active in one or more of the qHTS genotoxicity assays, indicating low sensitivity. This suggests that these qHTS assays cannot in their current form be used to replace traditional genotoxicity assays. However, despite the low sensitivity, ranking chemicals by potency of response in the qHTS assays revealed an enrichment for genotoxicants up to 12-fold compared with random selection, when allowing a 1% false positive rate. This finding indicates these qHTS assays can be used to prioritize chemicals for further investigation, allowing resources to focus on compounds most likely to induce genotoxic effects. To refine this prioritization process, models for predicting the genotoxicity potential of chemicals that were active in Tox21 genotoxicity assays were constructed using all Tox21 assay data, yielding a prediction accuracy up to 0.83. Data from qHTS assays related to stress-response pathway signaling (including genotoxicity) were the most informative for model construction. By using the results from qHTS genotoxicity assays, predictions from models based on qHTS data, and predictions from commercial bacterial mutagenicity QSAR models, we prioritized Tox21 chemicals for genotoxicity characterization.

Genetic toxicology in silico protocol.

In silico toxicology (IST) approaches to rapidly assess chemical hazard, and usage of such methods is increasing in all applications but especially for regulatory submissions, such as for assessing chemicals under REACH as well as the ICH M7 guideline for drug impurities. There are a number of obstacles to performing an IST assessment, including uncertainty in how such an assessment and associated expert review should be performed or what is fit for purpose, as well as a lack of confidence that the results will be accepted by colleagues, collaborators and regulatory authorities. To address this, a project to develop a series of IST protocols for different hazard endpoints has been initiated and this paper describes the genetic toxicity in silico (GIST) protocol. The protocol outlines a hazard assessment framework including key effects/mechanisms and their relationships to endpoints such as gene mutation and clastogenicity. IST models and data are reviewed that support the assessment of these effects/mechanisms along with defined approaches for combining the information and evaluating the confidence in the assessment. This protocol has been developed through a consortium of toxicologists, computational scientists, and regulatory scientists across several industries to support the implementation and acceptance of in silico approaches.

In silico toxicology protocols.

The present publication surveys several applications of in silico (i.e., computational) toxicology approaches across different industries and institutions. It highlights the need to develop standardized protocols when conducting toxicity-related predictions. This contribution articulates the information needed for protocols to support in silico predictions for major toxicological endpoints of concern (e.g., genetic toxicity, carcinogenicity, acute toxicity, reproductive toxicity, developmental toxicity) across several industries and regulatory bodies. Such novel in silico toxicology (IST) protocols, when fully developed and implemented, will ensure in silico toxicological assessments are performed and evaluated in a consistent, reproducible, and well-documented manner across industries and regulatory bodies to support wider uptake and acceptance of the approaches. The development of IST protocols is an initiative developed through a collaboration among an international consortium to reflect the state-of-the-art in in silico toxicology for hazard identification and characterization. A general outline for describing the development of such protocols is included and it is based on in silico predictions and/or available experimental data for a defined series of relevant toxicological effects or mechanisms. The publication presents a novel approach for determining the reliability of in silico predictions alongside experimental data. In addition, we discuss how to determine the level of confidence in the assessment based on the relevance and reliability of the information.

The Development of Target-Specific Pose Filter Ensembles To Boost Ligand Enrichment for Structure-Based Virtual Screening.

Structure-based virtual screening (SBVS) has become an indispensable technique for hit identification at the early stage of drug discovery. However, the accuracy of current scoring functions is not high enough to confer success to every target and thus remains to be improved. Previously, we had developed binary pose filters (PFs) using knowledge derived from the protein-ligand interface of a single X-ray structure of a specific target. This novel approach had been validated as an effective way to improve ligand enrichment. Continuing from it, in the present work we attempted to incorporate knowledge collected from diverse protein-ligand interfaces of multiple crystal structures of the same target to build PF ensembles (PFEs). Toward this end, we first constructed a comprehensive data set to meet the requirements of ensemble modeling and validation. This set contains 10 diverse targets, 118 well-prepared X-ray structures of protein-ligand complexes, and large benchmarking actives/decoys sets. Notably, we designed a unique workflow of two-layer classifiers based on the concept of ensemble learning and applied it to the construction of PFEs for all of the targets. Through extensive benchmarking studies, we demonstrated that (1) coupling PFE with Chemgauss4 significantly improves the early enrichment of Chemgauss4 itself and (2) PFEs show greater consistency in boosting early enrichment and larger overall enrichment than our prior PFs. In addition, we analyzed the pairwise topological similarities among cognate ligands used to construct PFEs and found that it is the higher chemical diversity of the cognate ligands that leads to the improved performance of PFEs. Taken together, the results so far prove that the incorporation of knowledge from diverse protein-ligand interfaces by ensemble modeling is able to enhance the screening competence of SBVS scoring functions.

An Intuitive Approach for Predicting Potential Human Health Risk with the Tox21 10k Library.

In vitro-in vivo extrapolation (IVIVE) analyses translating high-throughput screening (HTS) data to human relevance have been limited. This study represents the first report applying IVIVE approaches and exposure comparisons using the entirety of the Tox21 federal collaboration chemical screening data, incorporating assay response efficacy and quality of concentration-response fits, and providing quantitative anchoring to first address the likelihood of human in vivo interactions with Tox21 compounds. This likelihood was assessed using a maximum blood concentration to in vitro response ratio approach (Cmax/AC50), analogous to decision-making methods for clinical drug-drug interactions. Fraction unbound in plasma (fup) and intrinsic hepatic clearance (CLint) parameters were estimated in silico and incorporated in a three-compartment toxicokinetic (TK) model to first predict Cmax for in vivo corroboration using therapeutic scenarios. Toward lower exposure scenarios, 36 compounds of 3925 unique chemicals with curated activity in the HTS data using high-quality dose-response model fits and ≥40% efficacy gave "possible" human in vivo interaction likelihoods lower than median human exposures predicted in the United States Environmental Protection Agency's ExpoCast program. A publicly available web application has been designed to provide all Tox21-ToxCast dose-likelihood predictions. Overall, this approach provides an intuitive framework to relate in vitro toxicology data rapidly and quantitatively to exposures using either in vitro or in silico derived TK parameters and can be thought of as an important step toward estimating plausible biological interactions in a high-throughput risk-assessment framework.

Development of an in Vitro-Based Risk Assessment Framework for Predicting Ambient Particulate Matter-Bound Polycyclic Aromatic Hydrocarbon-Activated Toxicity Pathways.

Polycyclic aromatic hydrocarbons (PAHs) are widely distributed throughout the atmosphere as mixtures attached to ambient particulate matter (PM). PAHs usually elicit similar toxicological pathways but do so with varying levels of efficacy. In this study, we utilized high-throughput screening (HTS) in vitro data of PAHs to predict health risks associated with coarse and fine PM. PM samples with 22 PAH compounds obtained from residential areas close to industrial parks in central Taiwan were analyzed. On the basis of the PM-bound PAH concentrations and their activities reported in HTS assays, we developed a probabilistic model for estimating cumulative exposure of humans to PAHs. Activity-to-exposure ratio (AER) values were calculated to compare relative risks of activating the aryl hydrocarbon receptor (AhR), nuclear factor erythroid 2-related factor 2 (Nrf2), and tumor suppressor gene (p53) when children or adults were exposed to fine or coarse PM in different seasons. On the basis of AER values, the risk of fine PM exposure was relatively higher than the risk of exposure to coarse PM in pathway activation. Children as a susceptible population had a risk of the activating AhR pathway greater than that of adults. Particularly higher risks were observed in winter than in summer. Among three pathways, AhR was the most sensitive one activated by exposure to PAHs. In addition, the activation of the AhR, Nrf2, and p53 pathways was compared by in vitro reporter assays with and without the pre-extraction of PAHs from PM. Our proposed novel approach accounts for mixture toxicities in characterizing in vitro pathway-based risks via inhalation exposure to ambient PAHs.

Gender differences in cancer susceptibility: role of oxidative stress.

Cancer is a leading cause of death worldwide and environmental factors, including chemicals, have been suggested as major etiological incitements. Cancer statistics indicates that men get more cancer than women. However, differences in the known risk factors including life style or occupational exposure only offer partial explanation. Using a text mining tool, we have investigated the scientific literature concerning male- and female-specific rat carcinogens that induced tumors only in one gender in NTP 2-year cancer bioassay. Our evaluation shows that oxidative stress, although frequently reported for both male- and female-specific rat carcinogens, was mentioned significantly more in literature concerning male-specific rat carcinogens. Literature analysis of testosterone and estradiol showed the same pattern. Tox21 high-throughput assay results, although showing only weak association of oxidative stress-related processes for male- and female-specific rat carcinogens, provide additional support. We also analyzed the literature concerning 26 established human carcinogens (IARC group 1). Oxidative stress was more frequently reported for the majority of these carcinogens, and the Tox21 data resembled that of male-specific rat carcinogens. Thus, our data, based on about 600000 scientific abstracts and Tox21 screening assays, suggest a link between male-specific carcinogens, testosterone and oxidative stress. This implies that a different cellular response to oxidative stress in men and women may be a critical factor in explaining the greater cancer susceptibility observed in men. Although the IARC carcinogens are classified as human carcinogens, their classification largely based on epidemiological evidence from male cohorts, which raises the question whether carcinogen classifications should be gender specific.

Differential modulation of FXR activity by chlorophacinone and ivermectin analogs.

Chemicals that alter normal function of farnesoid X receptor (FXR) have been shown to affect the homeostasis of bile acids, glucose, and lipids. Several structural classes of environmental chemicals and drugs that modulated FXR transactivation were previously identified by quantitative high-throughput screening (qHTS) of the Tox21 10K chemical collection. In the present study, we validated the FXR antagonist activity of selected structural classes, including avermectin anthelmintics, dihydropyridine calcium channel blockers, 1,3-indandione rodenticides, and pyrethroid pesticides, using in vitro assay and quantitative structural-activity relationship (QSAR) analysis approaches. (Z)-Guggulsterone, chlorophacinone, ivermectin, and their analogs were profiled for their ability to alter CDCA-mediated FXR binding using a panel of 154 coregulator motifs and to induce or inhibit transactivation and coactivator recruitment activities of constitutive androstane receptor (CAR), liver X receptor alpha (LXRα), or pregnane X receptor (PXR). Our results showed that chlorophacinone and ivermectin had distinct modes of action (MOA) in modulating FXR-coregulator interactions and compound selectivity against the four aforementioned functionally-relevant nuclear receptors. These findings collectively provide mechanistic insights regarding compound activities against FXR and possible explanations for in vivo toxicological observations of chlorophacinone, ivermectin, and their analogs.

Integrated Approach for Testing and Assessment for Developmental Neurotoxicity (DNT) to Prioritize Aromatic Organophosphorus Flame Retardants.

Organophosphorus flame retardants (OPFRs) are abundant and persistent in the environment but have limited toxicity information. Their similarity in structure to organophosphate pesticides presents great concern for developmental neurotoxicity (DNT). However, current in vivo testing is not suitable to provide DNT information on the amount of OPFRs that lack data. Over the past decade, an in vitro battery was developed to enhance DNT assessment, consisting of assays that evaluate cellular processes in neurodevelopment and function. In this study, behavioral data of small model organisms were also included. To assess if these assays provide sufficient mechanistic coverage to prioritize chemicals for further testing and/or identify hazards, an integrated approach to testing and assessment (IATA) was developed with additional information from the Integrated Chemical Environment (ICE) and the literature. Human biomonitoring and exposure data were identified and physiologically-based toxicokinetic models were applied to relate in vitro toxicity data to human exposure based on maximum plasma concentration. Eight OPFRs were evaluated, including aromatic OPFRs (triphenyl phosphate (TPHP), isopropylated phenyl phosphate (IPP), 2-ethylhexyl diphenyl phosphate (EHDP), tricresyl phosphate (TMPP), isodecyl diphenyl phosphate (IDDP), tert-butylphenyl diphenyl phosphate (BPDP)) and halogenated FRs ((Tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tris(2-chloroethyl) phosphate (TCEP)). Two representative brominated flame retardants (BFRs) (2,2'4,4'-tetrabromodiphenyl ether (BDE-47) and 3,3',5,5'-tetrabromobisphenol A (TBBPA)) with known DNT potential were selected for toxicity benchmarking. Data from the DNT battery indicate that the aromatic OPFRs have activity at similar concentrations as the BFRs and should therefore be evaluated further. However, these assays provide limited information on the mechanism of the compounds. By integrating information from ICE and the literature, endocrine disruption was identified as a potential mechanism. This IATA case study indicates that human exposure to some OPFRs could lead to a plasma concentration similar to those exerting in vitro activities, indicating potential concern for human health.

An ensemble machine learning model generates a focused screening library for the identification of CDK8 inhibitors.

The identification of an effective inhibitor is an important starting step in drug development. Unfortunately, many issues such as the characterization of protein binding sites, the screening library, materials for assays, etc., make drug screening a difficult proposition. As the size of screening libraries increases, more resources will be inefficiently consumed. Thus, new strategies are needed to preprocess and focus a screening library towards a targeted protein. Herein, we report an ensemble machine learning (ML) model to generate a CDK8-focused screening library. The ensemble model consists of six different algorithms optimized for CDK8 inhibitor classification. The models were trained using a CDK8-specific fragment library along with molecules containing CDK8 activity. The optimized ensemble model processed a commercial library containing 1.6 million molecules. This resulted in a CDK8-focused screening library containing 1,672 molecules, a reduction of more than 99.90%. The CDK8-focused library was then subjected to molecular docking, and 25 candidate compounds were selected. Enzymatic assays confirmed six CDK8 inhibitors, with one compound producing an IC50 value of ≤100 nM. Analysis of the ensemble ML model reveals the role of the CDK8 fragment library during training. Structural analysis of molecules reveals the hit compounds to be structurally novel CDK8 inhibitors. Together, the results highlight a pipeline for curating a focused library for a specific protein target, such as CDK8.

Interlaboratory Study on Zebrafish in Toxicology: Systematic Evaluation of the Application of Zebrafish in Toxicology's (SEAZIT's) Evaluation of Developmental Toxicity.

Embryonic zebrafish represent a useful test system to screen substances for their ability to perturb development. The exposure scenarios, endpoints captured, and data analysis vary among the laboratories who conduct screening. A lack of harmonization impedes the comparison of the substance potency and toxicity outcomes across laboratories and may hinder the broader adoption of this model for regulatory use. The Systematic Evaluation of the Application of Zebrafish in Toxicology (SEAZIT) initiative was developed to investigate the sources of variability in toxicity testing. This initiative involved an interlaboratory study to determine whether experimental parameters altered the developmental toxicity of a set of 42 substances (3 tested in duplicate) in three diverse laboratories. An initial dose-range-finding study using in-house protocols was followed by a definitive study using four experimental conditions: chorion-on and chorion-off using both static and static renewal exposures. We observed reasonable agreement across the three laboratories as 33 of 42 test substances (78.6%) had the same activity call. However, the differences in potency seen using variable in-house protocols emphasizes the importance of harmonization of the exposure variables under evaluation in the second phase of this study. The outcome of the Def will facilitate future practical discussions on harmonization within the zebrafish research community.

Establishment of extensive artificial intelligence models for kinase inhibitor prediction: Identification of novel PDGFRB inhibitors.

Identifying hit compounds is an important step in drug development. Unfortunately, this process continues to be a challenging task. Several machine learning models have been generated to aid in simplifying and improving the prediction of candidate compounds. Models tuned for predicting kinase inhibitors have been established. However, an effective model can be limited by the size of the chosen training dataset. In this study, we tested several machine learning models to predict potential kinase inhibitors. A dataset was curated from a number of publicly available repositories. This resulted in a comprehensive dataset covering more than half of the human kinome. More than 2,000 kinase models were established using different model approaches. The performances of the models were compared, and the Keras-MLP model was determined to be the best performing model. The model was then used to screen a chemical library for potential inhibitors targeting platelet-derived growth factor receptor-ß (PDGFRB). Several PDGFRB candidates were selected, and in vitro assays confirmed four compounds with PDGFRB inhibitory activity and IC50 values in the nanomolar range. These results show the effectiveness of machine learning models trained on the reported dataset. This report would aid in the establishment of machine learning models as well as in the discovery of novel kinase inhibitors.

Systematic Evaluation of the Application of Zebrafish in Toxicology (SEAZIT): Developing a Data Analysis Pipeline for the Assessment of Developmental Toxicity with an Interlaboratory Study.

The embryonic zebrafish is a useful vertebrate model for assessing the effects of substances on growth and development. However, cross-laboratory developmental toxicity outcomes can vary and reported developmental defects in zebrafish may not be directly comparable between laboratories. To address these limitations for gaining broader adoption of the zebrafish model for toxicological screening, we established the Systematic Evaluation of the Application of Zebrafish in Toxicology (SEAZIT) program to investigate how experimental protocol differences can influence chemical-mediated effects on developmental toxicity (i.e., mortality and the incidence of altered phenotypes). As part of SEAZIT, three laboratories were provided a common and blinded dataset (42 substances) to evaluate substance-mediated effects on developmental toxicity in the embryonic zebrafish model. To facilitate cross-laboratory comparisons, all the raw experimental data were collected, stored in a relational database, and analyzed with a uniform data analysis pipeline. Due to variances in laboratory-specific terminology for altered phenotypes, we utilized ontology terms available from the Ontology Lookup Service (OLS) for Zebrafish Phenotype to enable additional cross-laboratory comparisons. In this manuscript, we utilized data from the first phase of screening (dose range finding, DRF) to highlight the methodology associated with the development of the database and data analysis pipeline, as well as zebrafish phenotype ontology mapping.

Tox21Enricher-Shiny: an R Shiny application for toxicity functional annotation analysis.

Inference of toxicological and mechanistic properties of untested chemicals through structural or biological similarity is a commonly employed approach for initial chemical characterization and hypothesis generation. We previously developed a web-based application, Tox21Enricher-Grails, on the Grails framework that identifies enriched biological/toxicological properties of chemical sets for the purpose of inferring properties of untested chemicals within the set. It was able to detect significantly overrepresented biological (e.g., receptor binding), toxicological (e.g., carcinogenicity), and chemical (e.g., toxicologically relevant chemical substructures) annotations within sets of chemicals screened in the Tox21 platform. Here, we present an R Shiny application version of Tox21Enricher-Grails, Tox21Enricher-Shiny, with more robust features and updated annotations. Tox21Enricher-Shiny allows users to interact with the web application component (available at http://hurlab.med.und.edu/Tox21Enricher/) through a user-friendly graphical user interface or to directly access the application's functions through an application programming interface. This version now supports InChI strings as input in addition to CASRN and SMILES identifiers. Input chemicals that contain certain reactive functional groups (nitrile, aldehyde, epoxide, and isocyanate groups) may react with proteins in cell-based Tox21 assays: this could cause Tox21Enricher-Shiny to produce spurious enrichment analysis results. Therefore, this version of the application can now automatically detect and ignore such problematic chemicals in a user's input. The application also offers new data visualizations, and the architecture has been greatly simplified to allow for simple deployment, version control, and porting. The application may be deployed onto a Posit Connect or Shiny server, and it uses Postgres for database management. As other Tox21-related tools are being migrated to the R Shiny platform, the development of Tox21Enricher-Shiny is a logical transition to use R's strong data analysis and visualization capacities and to provide aesthetic and developmental consistency with other Tox21 applications developed by the Division of Translational Toxicology (DTT) at the National Institute of Environmental Health Sciences (NIEHS).

Corrigendum: Tox21Enricher-Shiny: an R Shiny application for toxicity functional annotation analysis.

[This corrects the article DOI: 10.3389/ftox.2023.1147608.].