GLUT5, a fructose-transporting member of the facilitative glucose transporter (GLUT, SLC2) family, is a therapeutic target for diabetes and cancer but has no potent inhibitors. We virtually screened a library of 6 million chemicals onto a GLUT5 model and identified N-[4-(methylsulfonyl)-2-nitrophenyl]-1,3-benzodioxol-5-amine (MSNBA) as an inhibitor of GLUT5 fructose transport in proteoliposomes. MSNBA inhibition was specific to GLUT5; this inhibitor did not affect the fructose transport of human GLUT2 or the glucose transport of human GLUT1-4 or bacterial GlcPSe. In MCF7 cells, a human breast cancer cell line, MSNBA competitively inhibited GLUT5 fructose uptake with a KI of 3.2 ± 0.4 μM. Ligand docking, mutagenesis and functional studies indicate that MSNBA binds near the active site and inhibitor discrimination involves H387 of GLUT5. Thus, MSNBA is a selective and potent inhibitor of fructose transport via GLUT5, and the first chemical probe for this transporter. Our data indicate that active site differences in GLUT members could be exploited to further enhance ligand specificity.

The 'druggable genome' encompasses several protein families, but only a subset of targets within them have attracted significant research attention and thus have information about them publicly available. The Illuminating the Druggable Genome (IDG) program was initiated in 2014, has the goal of developing experimental techniques and a Knowledge Management Center (KMC) that would collect and organize information about protein targets from four families, representing the most common druggable targets with an emphasis on understudied proteins. Here, we describe two resources developed by the KMC: the Target Central Resource Database (TCRD) which collates many heterogeneous gene/protein datasets and Pharos (https://pharos.nih.gov), a multimodal web interface that presents the data from TCRD. We briefly describe the types and sources of data considered by the KMC and then highlight features of the Pharos interface designed to enable intuitive access to the IDG knowledgebase. The aim of Pharos is to encourage 'serendipitous browsing', whereby related, relevant information is made easily discoverable. We conclude by describing two use cases that highlight the utility of Pharos and TCRD.

Rho family GTPases (including Rac, Rho and Cdc42) collectively control cell proliferation, adhesion and migration and are of interest as functional therapeutic targets in numerous epithelial cancers. Based on high throughput screening of the Prestwick Chemical Library® and cheminformatics we identified the R-enantiomers of two approved drugs (naproxen and ketorolac) as inhibitors of Rac1 and Cdc42. The corresponding S-enantiomers are considered the active component in racemic drug formulations, acting as non-steroidal anti-inflammatory drugs (NSAIDs) with selective activity against cyclooxygenases. Here, we show that the S-enantiomers of naproxen and ketorolac are inactive against the GTPases. Additionally, more than twenty other NSAIDs lacked inhibitory action against the GTPases, establishing the selectivity of the two identified NSAIDs. R-naproxen was first identified as a lead compound and tested in parallel with its S-enantiomer and the non-chiral 6-methoxy-naphthalene acetic acid (active metabolite of nabumetone, another NSAID) as a structural series. Cheminformatics-based substructure analyses-using the rotationally constrained carboxylate in R-naproxen-led to identification of racemic [R/S] ketorolac as a suitable FDA-approved candidate. Cell based measurement of GTPase activity (in animal and human cell lines) demonstrated that the R-enantiomers specifically inhibit epidermal growth factor stimulated Rac1 and Cdc42 activation. The GTPase inhibitory effects of the R-enantiomers in cells largely mimic those of established Rac1 (NSC23766) and Cdc42 (CID2950007/ML141) specific inhibitors. Docking predicts that rotational constraints position the carboxylate moieties of the R-enantiomers to preferentially coordinate the magnesium ion, thereby destabilizing nucleotide binding to Rac1 and Cdc42. The S-enantiomers can be docked but are less favorably positioned in proximity to the magnesium. R-naproxen and R-ketorolac have potential for rapid translation and efficacy in the treatment of several epithelial cancer types on account of established human toxicity profiles and novel activities against Rho-family GTPases.

Kaposi's sarcoma-associated herpesvirus (KSHV) is linked to human malignancies. The majority of tumor cells harbor latent virus, and a small percentage undergo spontaneous lytic replication. Both latency and lytic replication are important for viral pathogenesis and spread, but the cellular players involved in the switch between the two viral life-cycle phases are not clearly understood. We conducted a small interfering RNA (siRNA) screen targeting the cellular kinome and identified Tousled-like kinases (TLKs) as cellular kinases that control KSHV reactivation from latency. Upon treatment of latent KSHV-infected cells with siRNAs targeting TLKs, we saw robust viral reactivation. Knockdown of TLKs in latent KSHV-infected cells induced expression of viral lytic proteins and production of infectious virus. TLKs were also found to play a role in regulating reactivation from latency of another related oncogenic gammaherpesvirus, Epstein-Barr virus. Our results establish the TLKs as cellular repressors of gammaherpesvirus reactivation.

Multiplexed small molecule inhibitors covalently bound to Sepharose beads (MIBs) were used to capture functional kinases in luminal, HER2-enriched and triple negative (basal-like and claudin-low) breast cancer cell lines and tumors. Kinase MIB-binding profiles at baseline without perturbation proteomically distinguished the four breast cancer subtypes. Understudied kinases, whose disease associations and pharmacology are generally unexplored, were highly represented in MIB-binding taxonomies and are integrated into signaling subnetworks with kinases that have been previously well characterized in breast cancer. Computationally it was possible to define subtypes using profiles of less than 50 of the more than 300 kinases bound to MIBs that included understudied as well as metabolic and lipid kinases. Furthermore, analysis of MIB-binding profiles established potential functional annotations for these understudied kinases. Thus, comprehensive MIBs-based capture of kinases provides a unique proteomics-based method for integration of poorly characterized kinases of the understudied kinome into functional subnetworks in breast cancer cells and tumors that is not possible using genomic strategies. The MIB-binding profiles readily defined subtype-selective differential adaptive kinome reprogramming in response to targeted kinase inhibition, demonstrating how MIB profiles can be used in determining dynamic kinome changes that result in subtype selective phenotypic state changes.

The human kinome consists of roughly 500 kinases, including 150 that have been proposed as therapeutic targets. Protein kinases regulate an array of signalling pathways that control metabolism, cell cycle progression, cell death, differentiation and survival. It is not surprising, then, that new kinase inhibitors developed to treat cancer, including sorafenib, also exhibit cardiotoxicity. We hypothesized that sorafenib cardiotoxicity is related to its deleterious effects on specific cardiac metabolic pathways given the critical roles of protein kinases in cardiac metabolism.

Rett syndrome (RTT) is a neurodevelopmental disorder caused by mutations in the X-linked methyl-CpG binding protein 2 (MeCP2) gene. While MeCP2 mutations are lethal in most males, females survive birth but show severe neurological defects. Because X-chromosome inactivation (XCI) is a random process, approximately 50% of the cells silence the wild-type (WT) copy of the MeCP2 gene. Thus, reactivating the silent WT copy of MeCP2 could provide therapeutic intervention for RTT.

DrugCentral is a public resource (http://drugcentral.org) that serves the scientific community by providing up-to-date drug information, as described in previous papers. The current release includes 109 newly approved (October 2018 through March 2020) active pharmaceutical ingredients in the US, Europe, Japan and other countries; and two molecular entities (e.g. mefuparib) of interest for COVID19. New additions include a set of pharmacokinetic properties for ∼1000 drugs, and a sex-based separation of side effects, processed from FAERS (FDA Adverse Event Reporting System); as well as a drug repositioning prioritization scheme based on the market availability and intellectual property rights forFDA approved drugs. In the context of the COVID19 pandemic, we also incorporated REDIAL-2020, a machine learning platform that estimates anti-SARS-CoV-2 activities, as well as the 'drugs in news' feature offers a brief enumeration of the most interesting drugs at the present moment. The full database dump and data files are available for download from the DrugCentral web portal.

This Formal Comment responds to a recent Meta-Research Article by identifying initiatives that are already in place for funding risky exploratory research that illuminate mysteries of the dark genome.

Inhibition of the HER2/ERBB2 receptor is a keystone to treating HER2-positive malignancies, particularly breast cancer, but a significant fraction of HER2-positive (HER2+) breast cancers recur or fail to respond. Anti-HER2 monoclonal antibodies, like trastuzumab or pertuzumab, and ATP active site inhibitors like lapatinib, commonly lack durability because of adaptive changes in the tumor leading to resistance. HER2+ cell line responses to inhibition with lapatinib were analyzed by RNAseq and ChIPseq to characterize transcriptional and epigenetic changes. Motif analysis of lapatinib-responsive genomic regions implicated the pioneer transcription factor FOXA1 as a mediator of adaptive responses. Lapatinib in combination with FOXA1 depletion led to dysregulation of enhancers, impaired adaptive upregulation of HER3, and decreased proliferation. HER2-directed therapy using clinically relevant drugs (trastuzumab with or without lapatinib or pertuzumab) in a 7-day clinical trial designed to examine early pharmacodynamic response to antibody-based anti-HER2 therapy showed reduced FOXA1 expression was coincident with decreased HER2 and HER3 levels, decreased proliferation gene signatures, and increased immune gene signatures. This highlights the importance of the immune response to anti-HER2 antibodies and suggests that inhibiting FOXA1-mediated adaptive responses in combination with HER2 targeting is a potential therapeutic strategy.

Dysregulation of PI3K/Akt signaling is a dominant feature in basal-like or triple-negative breast cancers (TNBC). However, the mechanisms regulating this pathway are largely unknown in this subset of aggressive tumors. Here we demonstrate that the transcription factor SOX4 is a key regulator of PI3K signaling in TNBC. Genomic and proteomic analyses coupled with mechanistic studies identified TGFBR2 as a direct transcriptional target of SOX4 and demonstrated that TGFBR2 is required to mediate SOX4-dependent PI3K signaling. We further report that SOX4 and the SWI/SNF ATPase SMARCA4, which are uniformly overexpressed in basal-like tumors, form a previously unreported complex that is required to maintain an open chromatin conformation at the TGFBR2 regulatory regions in order to mediate TGFBR2 expression and PI3K signaling. Collectively, our findings delineate the mechanism by which SOX4 and SMARCA4 cooperatively regulate PI3K/Akt signaling and suggest that this complex may play an essential role in TNBC genesis and/or progression.

Despite decades of intensive search for compounds that modulate the activity of particular protein targets, a large proportion of the human kinome remains as yet undrugged. Effective approaches are therefore required to map the massive space of unexplored compound-kinase interactions for novel and potent activities. Here, we carry out a crowdsourced benchmarking of predictive algorithms for kinase inhibitor potencies across multiple kinase families tested on unpublished bioactivity data. We find the top-performing predictions are based on various models, including kernel learning, gradient boosting and deep learning, and their ensemble leads to a predictive accuracy exceeding that of single-dose kinase activity assays. We design experiments based on the model predictions and identify unexpected activities even for under-studied kinases, thereby accelerating experimental mapping efforts. The open-source prediction algorithms together with the bioactivities between 95 compounds and 295 kinases provide a resource for benchmarking prediction algorithms and for extending the druggable kinome.

Birth defects are functional and structural abnormalities that impact about 1 in 33 births in the United States. They have been attributed to genetic and other factors such as drugs, cosmetics, food, and environmental pollutants during pregnancy, but for most birth defects there are no known causes.

DrugCentral, accessible at https://drugcentral.org , is an open-access online drug information repository. It covers over 4950 drugs, incorporating structural, physicochemical, and pharmacological details to support drug discovery, development, and repositioning. With around 20,000 bioactivity data points, manual curation enhances information from several major digital sources. Approximately 724 mechanism-of-action (MoA) targets offer updated drug target insights. The platform captures clinical data: over 14,300 on- and off-label uses, 27,000 contraindications, and around 340,000 adverse drug events from pharmacovigilance reports. DrugCentral encompasses information from molecular structures to marketed formulations, providing a comprehensive pharmaceutical reference. Users can easily navigate basic drug information and key features, making DrugCentral a versatile, unique resource. Furthermore, we present a use-case example where we utilize experimentally determined data from DrugCentral to support drug repurposing. A minimum activity threshold t should be considered against novel targets to repurpose a drug. Analyzing 1156 bioactivities for human MoA targets suggests a general threshold of 1 µM: t = 6 when expressed as - log[Activity(M)]). This applies to 87% of the drugs. Moreover, t can be refined empirically based on water solubility (S): t = 3 - logS, for logS < - 3. Alongside the drug repurposing classification scheme, which considers intellectual property rights, market exclusivity protections, and market accessibility, DrugCentral provides valuable data to prioritize candidates for drug repurposing programs efficiently.

ABCB6 is a member of the adenosine triphosphate (ATP)-binding cassette family of transporter proteins that is increasingly recognized as a relevant physiological and therapeutic target. Evaluation of modulators of ABCB6 activity would pave the way toward a more complete understanding of the significance of this transport process in tumor cell growth, proliferation and therapy-related drug resistance. In addition, this effort would improve our understanding of the function of ABCB6 in normal physiology with respect to heme biosynthesis, and cellular adaptation to metabolic demand and stress responses. To search for modulators of ABCB6, we developed a novel cell-based approach that, in combination with flow cytometric high-throughput screening (HTS), can be used to identify functional modulators of ABCB6. Accumulation of protoporphyrin, a fluorescent molecule, in wild-type ABCB6 expressing K562 cells, forms the basis of the HTS assay. Screening the Prestwick Chemical Library employing the HTS assay identified four compounds, benzethonium chloride, verteporfin, tomatine hydrochloride and piperlongumine, that reduced ABCB6 mediated cellular porphyrin levels. Validation of the identified compounds employing the hemin-agarose affinity chromatography and mitochondrial transport assays demonstrated that three out of the four compounds were capable of inhibiting ABCB6 mediated hemin transport into isolated mitochondria. However, only verteporfin and tomatine hydrochloride inhibited ABCB6's ability to compete with hemin as an ABCB6 substrate. This assay is therefore sensitive, robust, and suitable for automation in a high-throughput environment as demonstrated by our identification of selective functional modulators of ABCB6. Application of this assay to other libraries of synthetic compounds and natural products is expected to identify novel modulators of ABCB6 activity.

Protein kinases play key roles in oncogenic signaling and are a major focus in the development of targeted cancer therapies. Imatinib, a BCR-Abl tyrosine kinase inhibitor, is a successful front-line treatment for chronic myelogenous leukemia (CML). However, resistance to imatinib may be acquired by BCR-Abl mutations or hyperactivation of Src family kinases such as Lyn. We have used multiplexed kinase inhibitor beads (MIBs) and quantitative mass spectrometry (MS) to compare kinase expression and activity in an imatinib-resistant (MYL-R) and -sensitive (MYL) cell model of CML. Using MIB/MS, expression and activity changes of over 150 kinases were quantitatively measured from various protein kinase families. Statistical analysis of experimental replicates assigned significance to 35 of these kinases, referred to as the MYL-R kinome profile. MIB/MS and immunoblotting confirmed the over-expression and activation of Lyn in MYL-R cells and identified additional kinases with increased (MEK, ERK, IKKα, PKCβ, NEK9) or decreased (Abl, Kit, JNK, ATM, Yes) abundance or activity. Inhibiting Lyn with dasatinib or by shRNA-mediated knockdown reduced the phosphorylation of MEK and IKKα. Because MYL-R cells showed elevated NF-κB signaling relative to MYL cells, as demonstrated by increased IκBα and IL-6 mRNA expression, we tested the effects of an IKK inhibitor (BAY 65-1942). MIB/MS and immunoblotting revealed that BAY 65-1942 increased MEK/ERK signaling and that this increase was prevented by co-treatment with a MEK inhibitor (AZD6244). Furthermore, the combined inhibition of MEK and IKKα resulted in reduced IL-6 mRNA expression, synergistic loss of cell viability and increased apoptosis. Thus, MIB/MS analysis identified MEK and IKKα as important downstream targets of Lyn, suggesting that co-targeting these kinases may provide a unique strategy to inhibit Lyn-dependent imatinib-resistant CML. These results demonstrate the utility of MIB/MS as a tool to identify dysregulated kinases and to interrogate kinome dynamics as cells respond to targeted kinase inhibition.

Many bioactivity databases offer information regarding the biological activity of small molecules on protein targets. Information in these databases is often hard to resolve with certainty because of subsetting different data in a variety of formats; use of different bioactivity metrics; use of different identifiers for chemicals and proteins; and having to access different query interfaces, respectively. Given the multitude of data sources, interfaces and standards, it is challenging to gather relevant facts and make appropriate connections and decisions regarding chemical-protein associations. The CARLSBAD database has been developed as an integrated resource, focused on high-quality subsets from several bioactivity databases, which are aggregated and presented in a uniform manner, suitable for the study of the relationships between small molecules and targets. In contrast to data collection resources, CARLSBAD provides a single normalized activity value of a given type for each unique chemical-protein target pair. Two types of scaffold perception methods have been implemented and are available for datamining: HierS (hierarchical scaffolds) and MCES (maximum common edge subgraph). The 2012 release of CARLSBAD contains 439 985 unique chemical structures, mapped onto 1,420 889 unique bioactivities, and annotated with 277 140 HierS scaffolds and 54 135 MCES chemical patterns, respectively. Of the 890 323 unique structure-target pairs curated in CARLSBAD, 13.95% are aggregated from multiple structure-target values: 94 975 are aggregated from two bioactivities, 14 544 from three, 7 930 from four and 2214 have five bioactivities, respectively. CARLSBAD captures bioactivities and tags for 1435 unique chemical structures of active pharmaceutical ingredients (i.e. 'drugs'). CARLSBAD processing resulted in a net 17.3% data reduction for chemicals, 34.3% reduction for bioactivities, 23% reduction for HierS and 25% reduction for MCES, respectively. The CARLSBAD database supports a knowledge mining system that provides non-specialists with novel integrative ways of exploring chemical biology space to facilitate knowledge mining in drug discovery and repurposing. Database URL: http://carlsbad.health.unm.edu/carlsbad/.

Cerebral cavernous malformations (CCM) are vascular abnormalities of the central nervous system predisposing blood vessels to leakage, leading to hemorrhagic stroke. Three genes, Krit1 (CCM1), OSM (CCM2), and PDCD10 (CCM3) are involved in CCM development. PDCD10 binds specifically to PtdIns(3,4,5)P3 and OSM. Using threading analysis and multi-template modeling, we constructed a three-dimensional model of PDCD10. PDCD10 appears to be a six-helical-bundle protein formed by two heptad-repeat-hairpin structures (alpha1-3 and alpha4-6) sharing the closest 3D homology with the bacterial phosphate transporter, PhoU. We identified a stretch of five lysines forming an amphipathic helix, a potential PtdIns(3,4,5)P3 binding site, in the alpha5 helix. We generated a recombinant wild-type (WT) and three PDCD10 mutants that have two (Delta2KA), three (Delta3KA), and five (Delta5KA) K to A mutations. Delta2KA and Delta3KA mutants hypothetically lack binding residues to PtdIns(3,4,5)P3 at the beginning and the end of predicted helix, while Delta5KA completely lacks all predicted binding residues. The WT, Delta2KA, and Delta3KA mutants maintain their binding to PtdIns(3,4,5)P3. Only the Delta5KA abolishes binding to PtdIns(3,4,5)P3. Both Delta5KA and WT show similar secondary and tertiary structures; however, Delta5KA does not bind to OSM. When WT and Delta5KA are co-expressed with membrane-bound constitutively-active PI3 kinase (p110-CAAX), the majority of the WT is co-localized with p110-CAAX at the plasma membrane where PtdIns(3,4,5)P3 is presumably abundant. In contrast, the Delta5KA remains in the cytoplasm and is not present in the plasma membrane. Combining computational modeling and biological data, we propose that the CCM protein complex functions in the PI3K signaling pathway through the interaction between PDCD10 and PtdIns(3,4,5)P3.

Neurofibromatosis type 2 (NF2) is a dominantly inherited autosomal disease characterized by schwannomas of the 8th cranial nerve. The NF2 tumor suppressor gene encodes for Merlin, a protein implicated as a suppressor of multiple cellular signaling pathways. To identify potential drug targets in NF2-associated malignancies we assessed the consequences of inhibiting the tyrosine kinase receptor MET. We identified crizotinib, a MET and ALK inhibitor, as a potent inhibitor of NF2-null Schwann cell proliferation in vitro and tumor growth in vivo. To identify the target/s of crizotnib we employed activity-based protein profiling (ABPP), leading to identification of FAK1 (PTK2) as the relevant target of crizotinib inhibition in NF2-null schwannoma cells. Subsequent studies confirm that inhibition of FAK1 is sufficient to suppress tumorigenesis in animal models of NF2 and that crizotinib-resistant forms of FAK1 can rescue the effects of treatment. These studies identify a FDA approved drug as a potential treatment for NF2 and delineate the mechanism of action in NF2-null Schwann cells.

Bioassay data analysis continues to be an essential, routine, yet challenging task in modern drug discovery and chemical biology research. The challenge is to infer reliable knowledge from big and noisy data. Some aspects of this problem are general with solutions informed by existing and emerging data science best practices. Some aspects are domain specific, and rely on expertise in bioassay methodology and chemical biology. Testing compounds for biological activity requires complex and innovative methodology, producing results varying widely in accuracy, precision, and information content. Hit selection criteria involve optimizing such that the overall probability of success in a project is maximized, and resource-wasteful "false trails" are avoided. This "fail-early" approach is embraced both in pharmaceutical and academic drug discovery, since follow-up capacity is resource-limited. Thus, early identification of likely promiscuous compounds has practical value.