Target identification remains a major challenge for modern drug discovery programs aimed at understanding the molecular mechanisms of drugs. Computational target prediction approaches like 2D chemical similarity searches have been widely used but are limited to structures sharing high chemical similarity. Here, we present a new computational approach called chemical similarity network analysis pull-down 3D (CSNAP3D) that combines 3D chemical similarity metrics and network algorithms for structure-based drug target profiling, ligand deorphanization, and automated identification of scaffold hopping compounds. In conjunction with 2D chemical similarity fingerprints, CSNAP3D achieved a >95% success rate in correctly predicting the drug targets of 206 known drugs. Significant improvement in target prediction was observed for HIV reverse transcriptase (HIVRT) compounds, which consist of diverse scaffold hopping compounds targeting the nucleotidyltransferase binding site. CSNAP3D was further applied to a set of antimitotic compounds identified in a cell-based chemical screen and identified novel small molecules that share a pharmacophore with Taxol and display a Taxol-like mechanism of action, which were validated experimentally using in vitro microtubule polymerization assays and cell-based assays.

Leucine carboxyl methyltransferase-1 (LCMT1) and protein phosphatase methylesterase-1 (PME-1) are essential enzymes that regulate the methylation of the protein phosphatase 2A catalytic subunit (PP2AC). LCMT1 and PME-1 have been linked to the regulation of cell growth and proliferation, but the underlying mechanisms have remained elusive. We show here an important role for an LCMT1-PME-1 methylation equilibrium in controlling mitotic spindle size. Depletion of LCMT1 or overexpression of PME-1 led to long spindles. In contrast, depletion of PME-1, pharmacological inhibition of PME-1 or overexpression of LCMT1 led to short spindles. Furthermore, perturbation of the LCMT1-PME-1 methylation equilibrium led to mitotic arrest, spindle assembly checkpoint activation, defective cell divisions, induction of apoptosis and reduced cell viability. Thus, we propose that the LCMT1-PME-1 methylation equilibrium is critical for regulating mitotic spindle size and thereby proper cell division.

Discovery of first-in-class medicines for treating cancer is limited by concerns with their toxicity and safety profiles, while repurposing known drugs for new anticancer indications has become a viable alternative. Here, we have developed a new approach that utilizes cell cycle arresting patterns as unique molecular signatures for prioritizing FDA-approved drugs with repurposing potential. As proof-of-principle, we conducted large-scale cell cycle profiling of 884 FDA-approved drugs. Using cell cycle indexes that measure changes in cell cycle profile patterns upon chemical perturbation, we identified 36 compounds that inhibited cancer cell viability including 6 compounds that were previously undescribed. Further cell cycle fingerprint analysis and 3D chemical structural similarity clustering identified unexpected FDA-approved drugs that induced DNA damage, including clinically relevant microtubule destabilizers, which was confirmed experimentally via cell-based assays. Our study shows that computational cell cycle profiling can be used as an approach for prioritizing FDA-approved drugs with repurposing potential, which could aid the development of cancer therapeutics.

Small molecules that bind to tubulin exert powerful effects on cell division and apoptosis (programmed cell death). Cell-based high-throughput screening combined with chemo/bioinformatic and biochemical analyses recently revealed a novel compound MI-181 as a potent mitotic inhibitor with heightened activity towards melanomas. MI-181 causes tubulin depolymerization, activates the spindle assembly checkpoint arresting cells in mitosis, and induces apoptotic cell death. C2 is an unrelated compound previously shown to have lethal effects on microtubules in tumorigenic cell lines. We report 2.60 Å and 3.75 Å resolution structures of MI-181 and C2, respectively, bound to a ternary complex of αβ-tubulin, the tubulin-binding protein stathmin, and tubulin tyrosine ligase. In the first of these structures, our crystallographic results reveal a unique binding mode for MI-181 extending unusually deep into the well-studied colchicine-binding site on β-tubulin. In the second structure the C2 compound occupies the colchicine-binding site on β-tubulin with two chemical moieties recapitulating contacts made by colchicine, in combination with another system of atomic contacts. These insights reveal the source of the observed effects of MI-181 and C2 on microtubules, mitosis, and cultured cancer cell lines. The structural details of the interaction between tubulin and the described compounds may guide the development of improved derivative compounds as therapeutic candidates or molecular probes to study cancer cell division.

The polo-like kinase 1 (Plk1) is a critical regulator of cell division that is overexpressed in many types of tumors. Thus, a strategy in the treatment of cancer has been to target the kinase activity (ATPase domain) or substrate-binding domain (Polo-box Domain, PBD) of Plk1. However, only few synthetic small molecules have been identified that target the Plk1-PBD. Here, we have applied an integrative approach that combines pharmacophore modeling, molecular docking, virtual screening, and in vitro testing to discover novel Plk1-PBD inhibitors. Nine Plk1-PBD crystal structures were used to generate structure-based hypotheses. A common pharmacophore model (Hypo1) composed of five chemical features was selected from the 9 structure-based hypotheses and used for virtual screening of a drug-like database consisting of 159,757 compounds to identify novel Plk1-PBD inhibitors. The virtual screening technique revealed 9,327 compounds with a maximum fit value of 3 or greater, which were selected and subjected to molecular docking analyses. This approach yielded 93 compounds that made good interactions with critical residues within the Plk1-PBD active site. The testing of these 93 compounds in vitro for their ability to inhibit the Plk1-PBD, showed that many of these compounds had Plk1-PBD inhibitory activity and that compound Chemistry_28272 was the most potent Plk1-PBD inhibitor. Thus Chemistry_28272 and the other top compounds are novel Plk1-PBD inhibitors and could be used for the development of cancer therapeutics.

Target identification is one of the most critical steps following cell-based phenotypic chemical screens aimed at identifying compounds with potential uses in cell biology and for developing novel disease therapies. Current in silico target identification methods, including chemical similarity database searches, are limited to single or sequential ligand analysis that have limited capabilities for accurate deconvolution of a large number of compounds with diverse chemical structures. Here, we present CSNAP (Chemical Similarity Network Analysis Pulldown), a new computational target identification method that utilizes chemical similarity networks for large-scale chemotype (consensus chemical pattern) recognition and drug target profiling. Our benchmark study showed that CSNAP can achieve an overall higher accuracy (>80%) of target prediction with respect to representative chemotypes in large (>200) compound sets, in comparison to the SEA approach (60-70%). Additionally, CSNAP is capable of integrating with biological knowledge-based databases (Uniprot, GO) and high-throughput biology platforms (proteomic, genetic, etc) for system-wise drug target validation. To demonstrate the utility of the CSNAP approach, we combined CSNAP's target prediction with experimental ligand evaluation to identify the major mitotic targets of hit compounds from a cell-based chemical screen and we highlight novel compounds targeting microtubules, an important cancer therapeutic target. The CSNAP method is freely available and can be accessed from the CSNAP web server (http://services.mbi.ucla.edu/CSNAP/).

Targeting the leukemia proliferation cycle has been a successful approach to developing antileukemic therapies. However, drug screening efforts to identify novel antileukemic agents have been hampered by the lack of a suitable high-throughput screening platform for suspension cells that does not rely on flow-cytometry analyses. We report the development of a novel leukemia cell-based high-throughput chemical screening platform for the discovery of cell cycle phase specific inhibitors that utilizes chemical cell cycle profiling. We have used this approach to analyze the cell cycle response of acute lymphoblastic leukemia CCRF-CEM cells to each of 181420 druglike compounds. This approach yielded cell cycle phase specific inhibitors of leukemia cell proliferation. Further analyses of the top G2-phase and M-phase inhibitors identified the leukemia specific inhibitor 1 (Leusin-1). Leusin-1 arrests cells in G2 phase and triggers an apoptotic cell death. Most importantly, Leusin-1 was more active in acute lymphoblastic leukemia cells than other types of leukemias, non-blood cancers, or normal cells and represents a lead molecule for developing antileukemic drugs.

STARD9 is a largely uncharacterized mitotic kinesin and putative cancer target that is critical for regulating pericentriolar material cohesion during bipolar spindle assembly. To begin to understand the mechanisms regulating STARD9 function and their importance to cell division, we took a multidisciplinary approach to define the cis and trans factors that regulate the stability of the STARD9 motor domain. We show that, unlike the other ∼50 mammalian kinesins, STARD9 contains an insertion in loop 12 of its motor domain (MD). Working with the STARD9-MD, we show that it is phosphorylated in mitosis by mitotic kinases that include Plk1. These phosphorylation events are important for targeting a pool of STARD9-MD for ubiquitination by the SCFβ-TrCP ubiquitin ligase and proteasome-dependent degradation. Of interest, overexpression of nonphosphorylatable/nondegradable STARD9-MD mutants leads to spindle assembly defects. Our results with STARD9-MD imply that in vivo the protein levels of full-length STARD9 could be regulated by Plk1 and SCFβ-TrCP to promote proper mitotic spindle assembly.

Short-rib polydactyly syndromes (SRPS) arise from mutations in genes involved in retrograde intraflagellar transport (IFT) and basal body homeostasis, which are critical for cilia assembly and function. Recently, mutations in WDR34 or WDR60 (candidate dynein intermediate chains) were identified in SRPS. We have identified and characterized Tctex1d2, which associates with Wdr34, Wdr60 and other dynein complex 1 and 2 subunits. Tctex1d2 and Wdr60 localize to the base of the cilium and their depletion causes defects in ciliogenesis. We propose that Tctex1d2 is a novel dynein light chain important for trafficking to the cilium and potentially retrograde IFT and is a new molecular link to understanding SRPS pathology.

The Katanin family of microtubule-severing enzymes is critical for remodeling microtubule-based structures that influence cell division, motility, morphogenesis and signaling. Katanin is composed of a catalytic p60 subunit (A subunit, KATNA1) and a regulatory p80 subunit (B subunit, KATNB1). The mammalian genome also encodes two additional A-like subunits (KATNAL1 and KATNAL2) and one additional B-like subunit (KATNBL1) that have remained poorly characterized. To better understand the factors and mechanisms controlling mammalian microtubule-severing, we have taken a mass proteomic approach to define the protein interaction module for each mammalian Katanin subunit and to generate the mammalian Katanin family interaction network (Katan-ome). Further, we have analyzed the function of the KATNBL1 subunit and determined that it associates with KATNA1 and KATNAL1, it localizes to the spindle poles only during mitosis and it regulates Katanin A subunit microtubule-severing activity in vitro Interestingly, during interphase, KATNBL1 is sequestered in the nucleus through an N-terminal nuclear localization signal. Finally KATNB1 was able to compete the interaction of KATNBL1 with KATNA1 and KATNAL1. These data indicate that KATNBL1 functions as a regulator of Katanin A subunit microtubule-severing activity during mitosis and that it likely coordinates with KATNB1 to perform this function.