Affinity Selection and Mass Spectrometry-Based Strategies to Identify Lead Compounds in Combinatorial Libraries

The screening of diverse libraries of small molecules created by combinatorial synthetic methods is a recent development which has the potential to accelerate the identification of lead compounds in drug discovery. We have developed a direct and rapid method to identify lead compounds in libraries involving affinity selection and mass spectrometry. In our strategy, the receptor or target molecule of interest is used to isolate the active components from the library physically, followed by direct structural identification of the active compounds bound to the target molecule by mass spectrometry. In a drug design strategy, structurally diverse libraries can be used for the initial identification of lead compounds. Once lead compounds have been identified, libraries containing compounds chemically similar to the lead compound can be generated and used to optimize the binding characteristics. These strategies have also been adopted for more detailed studies of protein–ligand interactions.     

Protein Interactions in Genome Maintenance as Novel Antibacterial Targets

Antibacterial compounds typically act by directly inhibiting essential bacterial enzyme activities. Although this general mechanism of action has fueled traditional antibiotic discovery efforts for decades, new antibiotic development has not kept pace with the emergence of drug resistant bacterial strains. These limitations have severely restricted the therapeutic tools available for treating bacterial infections. Here we test an alternative antibacterial lead-compound identification strategy in which essential protein-protein interactions are targeted rather than enzymatic activities. Bacterial single-stranded DNA-binding proteins (SSBs) form conserved protein interaction “hubs” that are essential for recruiting many DNA replication, recombination, and repair proteins to SSB/DNA nucleoprotein substrates. Three small molecules that block SSB/protein interactions are shown to have antibacterial activity against diverse bacterial species. Consistent with a model in which the compounds target multiple SSB/protein interactions, treatment of Bacillus subtilis cultures with the compounds leads to rapid inhibition of DNA replication and recombination, and ultimately to cell death. The compounds also have unanticipated effects on protein synthesis that could be due to a previously unknown role for SSB/protein interactions in translation or to off-target effects. Our results highlight the potential of targeting protein-protein interactions, particularly those that mediate genome maintenance, as a powerful approach for identifying new antibacterial compounds.     

Predicting drug targets based on protein domains

The identification of interactions between drugs and proteins plays key roles in understanding mechanisms underlying drug actions and can lead to new drug design strategies. Here, we present a novel statistical approach, namely PDTD (Predicting Drug Targets with Domains), to predict potential target proteins of new drugs based on derived interactions between drugs and protein domains. The known target proteins of those drugs that have similar therapeutic effects allow us to infer interactions between drugs and protein domains which in turn leads to identification of potential drug-protein interactions. Benchmarking with known drug-protein interactions shows that our proposed methodology outperforms previous methods that exploit either protein sequences or compound structures to predict drug targets, which demonstrates the predictive power of our proposed PDTD method.     

Protein domain-based prediction of drug/compound–target interactions and experimental validation on LIM kinases

Predictive approaches such as virtual screening have been used in drug discovery with the objective of reducing developmental time and costs. Current machine learning and network-based approaches have issues related to generalization, usability, or model interpretability, especially due to the complexity of target proteins’ structure/function, and bias in system training datasets. Here, we propose a new method “DRUIDom” (DRUg Interacting Domain prediction) to identify bio-interactions between drug candidate compounds and targets by utilizing the domain modularity of proteins, to overcome problems associated with current approaches. DRUIDom is composed of two methodological steps. First, ligands/compounds are statistically mapped to structural domains of their target proteins, with the aim of identifying their interactions. As such, other proteins containing the same mapped domain or domain pair become new candidate targets for the corresponding compounds. Next, a million-scale dataset of small molecule compounds, including those mapped to domains in the previous step, are clustered based on their molecular similarities, and their domain associations are propagated to other compounds within the same clusters. Experimentally verified bioactivity data points, obtained from public databases, are meticulously filtered to construct datasets of active/interacting and inactive/non-interacting drug/compound–target pairs (~2.9M data points), and used as training data for calculating parameters of compound–domain mappings, which led to 27,032 high-confidence associations between 250 domains and 8,165 compounds, and a finalized output of ~5 million new compound–protein interactions. DRUIDom is experimentally validated by syntheses and bioactivity analyses of compounds predicted to target LIM-kinase proteins, which play critical roles in the regulation of cell motility, cell cycle progression, and differentiation through actin filament dynamics. We showed that LIMK-inhibitor-2 and its derivatives significantly block the cancer cell migration through inhibition of LIMK phosphorylation and the downstream protein cofilin. One of the derivative compounds (LIMKi-2d) was identified as a promising candidate due to its action on resistant Mahlavu liver cancer cells. The results demonstrated that DRUIDom can be exploited to identify drug candidate compounds for intended targets and to predict new target proteins based on the defined compound–domain relationships. Datasets, results, and the source code of DRUIDom are fully-available at: https://github.com/cansyl/DRUIDom.     

The Overlap of Small Molecule and Protein Binding Sites within Families of Protein Structures

Protein–protein interactions are challenging targets for modulation by small molecules. Here, we propose an approach that harnesses the increasing structural coverage of protein complexes to identify small molecules that may target protein interactions. Specifically, we identify ligand and protein binding sites that overlap upon alignment of homologous proteins. Of the 2,619 protein structure families observed to bind proteins, 1,028 also bind small molecules (250–1000 Da), and 197 exhibit a statistically significant (p<0.01) overlap between ligand and protein binding positions. These “bi-functional positions”, which bind both ligands and proteins, are particularly enriched in tyrosine and tryptophan residues, similar to “energetic hotspots” described previously, and are significantly less conserved than mono-functional and solvent exposed positions. Homology transfer identifies ligands whose binding sites overlap at least 20% of the protein interface for 35% of domain–domain and 45% of domain–peptide mediated interactions. The analysis recovered known small-molecule modulators of protein interactions as well as predicted new interaction targets based on the sequence similarity of ligand binding sites. We illustrate the predictive utility of the method by suggesting structural mechanisms for the effects of sanglifehrin A on HIV virion production, bepridil on the cellular entry of anthrax edema factor, and fusicoccin on vertebrate developmental pathways. The results, available at http://pibase.janelia.org, represent a comprehensive collection of structurally characterized modulators of protein interactions, and suggest that homologous structures are a useful resource for the rational design of interaction modulators.     

The Role of Protein Interactions in Mediating Essentiality and Synthetic Lethality

Genes are characterized as essential if their knockout is associated with a lethal phenotype, and these “essential genes” play a central role in biological function. In addition, some genes are only essential when deleted in pairs, a phenomenon known as synthetic lethality. Here we consider genes displaying synthetic lethality as “essential pairs” of genes, and analyze the properties of yeast essential genes and synthetic lethal pairs together. As gene duplication initially produces an identical pair or sets of genes, it is often invoked as an explanation for synthetic lethality. However, we find that duplication explains only a minority of cases of synthetic lethality. Similarly, disruption of metabolic pathways leads to relatively few examples of synthetic lethality. By contrast, the vast majority of synthetic lethal gene pairs code for proteins with related functions that share interaction partners. We also find that essential genes and synthetic lethal pairs cluster in the protein-protein interaction network. These results suggest that synthetic lethality is strongly dependent on the formation of protein-protein interactions. Compensation by duplicates does not usually occur mainly because the genes involved are recent duplicates, but is more commonly due to functional similarity that permits preservation of essential protein complexes. This unified view, combining genes that are individually essential with those that form essential pairs, suggests that essentiality is a feature of physical interactions between proteins protein-protein interactions, rather than being inherent in gene and protein products themselves.     

A Multidimensional Strategy to Detect Polypharmacological Targets in the Absence of Structural and Sequence Homology

Conventional drug design embraces the “one gene, one drug, one disease” philosophy. Polypharmacology, which focuses on multi-target drugs, has emerged as a new paradigm in drug discovery. The rational design of drugs that act via polypharmacological mechanisms can produce compounds that exhibit increased therapeutic potency and against which resistance is less likely to develop. Additionally, identifying multiple protein targets is also critical for side-effect prediction. One third of potential therapeutic compounds fail in clinical trials or are later removed from the market due to unacceptable side effects often caused by off-target binding. In the current work, we introduce a multidimensional strategy for the identification of secondary targets of known small-molecule inhibitors in the absence of global structural and sequence homology with the primary target protein. To demonstrate the utility of the strategy, we identify several targets of 4,5-dihydroxy-3-(1-naphthyldiazenyl)-2,7-naphthalenedisulfonic acid, a known micromolar inhibitor of Trypanosoma brucei RNA editing ligase 1. As it is capable of identifying potential secondary targets, the strategy described here may play a useful role in future efforts to reduce drug side effects and/or to increase polypharmacology.     

Identification of Distant Drug Off-Targets by Direct Superposition of Binding Pocket Surfaces

Correctly predicting off-targets for a given molecular structure, which would have the ability to bind a large range of ligands, is both particularly difficult and important if they share no significant sequence or fold similarity with the respective molecular target (“distant off-targets”). A novel approach for identification of off-targets by direct superposition of protein binding pocket surfaces is presented and applied to a set of well-studied and highly relevant drug targets, including representative kinases and nuclear hormone receptors. The entire Protein Data Bank is searched for similar binding pockets and convincing distant off-target candidates were identified that share no significant sequence or fold similarity with the respective target structure. These putative target off-target pairs are further supported by the existence of compounds that bind strongly to both with high topological similarity, and in some cases, literature examples of individual compounds that bind to both. Also, our results clearly show that it is possible for binding pockets to exhibit a striking surface similarity, while the respective off-target shares neither significant sequence nor significant fold similarity with the respective molecular target (“distant off-target”).     

DeepCPI:A Deep Learning-based Framework for Large-scale in silico Drug Screening

Accurate identification of compound–protein interactions(CPIs) in silico may deepen our understanding of the underlying mechanisms of drug action and thus remarkably facilitate drug discovery and development.Conventional similarity-or docking-based computational methods for predicting CPIs rarely exploit latent features from currently available large-scale unlabeled compound and protein data and often limit their usage to relatively small-scale datasets.In the present study,we propose Deep CPI,a novel general and scalable computational framework that combines effective feature embedding(a technique of representation learning) with powerful deep learning methods to accurately predict CPIs at a large scale.Deep CPI automatically learns the implicit yet expressive low-dimensional features of compounds and proteins from a massive amount of unlabeled data.Evaluations of the measured CPIs in large-scale databases,such as Ch EMBL and Binding DB,as well as of the known drug–target interactions from Drug Bank,demonstrated the superior predictive performance of Deep CPI.Furthermore,several interactions among smallmolecule compounds and three G protein-coupled receptor targets(glucagon-like peptide-1 receptor,glucagon receptor,and vasoactive intestinal peptide receptor) predicted using Deep CPI were experimentally validated.The present study suggests that Deep CPI is a useful and powerful tool for drug discovery and repositioning.The source code of Deep CPI can be downloaded from https://github.com/Fangping Wan/Deep CPI.     

Knowledge-based Fragment Binding Prediction

Target-based drug discovery must assess many drug-like compounds for potential activity. Focusing on low-molecular-weight compounds (fragments) can dramatically reduce the chemical search space. However, approaches for determining protein-fragment interactions have limitations. Experimental assays are time-consuming, expensive, and not always applicable. At the same time, computational approaches using physics-based methods have limited accuracy. With increasing high-resolution structural data for protein-ligand complexes, there is now an opportunity for data-driven approaches to fragment binding prediction. We present FragFEATURE, a machine learning approach to predict small molecule fragments preferred by a target protein structure. We first create a knowledge base of protein structural environments annotated with the small molecule substructures they bind. These substructures have low-molecular weight and serve as a proxy for fragments. FragFEATURE then compares the structural environments within a target protein to those in the knowledge base to retrieve statistically preferred fragments. It merges information across diverse ligands with shared substructures to generate predictions. Our results demonstrate FragFEATURE''s ability to rediscover fragments corresponding to the ligand bound with 74% precision and 82% recall on average. For many protein targets, it identifies high scoring fragments that are substructures of known inhibitors. FragFEATURE thus predicts fragments that can serve as inputs to fragment-based drug design or serve as refinement criteria for creating target-specific compound libraries for experimental or computational screening.     

Target identification by chromatographic co-elution: monitoring of drug-protein interactions without immobilization or chemical derivatization

Bioactive molecules typically mediate their biological effects through direct physical association with one or more cellular proteins. The detection of drug-target interactions is therefore essential for the characterization of compound mechanism of action and off-target effects, but generic label-free approaches for detecting binding events in biological mixtures have remained elusive. Here, we report a method termed target identification by chromatographic co-elution (TICC) for routinely monitoring the interaction of drugs with cellular proteins under nearly physiological conditions in vitro based on simple liquid chromatographic separations of cell-free lysates. Correlative proteomic analysis of drug-bound protein fractions by shotgun sequencing is then performed to identify candidate target(s). The method is highly reproducible, does not require immobilization or derivatization of drug or protein, and is applicable to diverse natural products and synthetic compounds. The capability of TICC to detect known drug-protein target physical interactions (K(d) range: micromolar to nanomolar) is demonstrated both qualitatively and quantitatively. We subsequently used TICC to uncover the sterol biosynthetic enzyme Erg6p as a novel putative anti-fungal target. Furthermore, TICC identified Asc1 and Dak1, a core 40 S ribosomal protein that represses gene expression, and dihydroxyacetone kinase involved in stress adaptation, respectively, as novel yeast targets of a dopamine receptor agonist.     

Enhancements to the Rosetta Energy Function Enable Improved Identification of Small Molecules that Inhibit Protein-Protein Interactions

Protein-protein interactions are among today’s most exciting and promising targets for therapeutic intervention. To date, identifying small-molecules that selectively disrupt these interactions has proven particularly challenging for virtual screening tools, since these have typically been optimized to perform well on more “traditional” drug discovery targets. Here, we test the performance of the Rosetta energy function for identifying compounds that inhibit protein interactions, when these active compounds have been hidden amongst pools of “decoys.” Through this virtual screening benchmark, we gauge the effect of two recent enhancements to the functional form of the Rosetta energy function: the new “Talaris” update and the “pwSHO” solvation model. Finally, we conclude by developing and validating a new weight set that maximizes Rosetta’s ability to pick out the active compounds in this test set. Looking collectively over the course of these enhancements, we find a marked improvement in Rosetta’s ability to identify small-molecule inhibitors of protein-protein interactions.     

Protein structure similarity clustering (PSSC) and natural product structure as inspiration sources for drug development and chemical genomics

Finding small molecules that modulate protein function is of primary importance in drug development and in the emerging field of chemical genomics. To facilitate the identification of such molecules, we developed a novel strategy making use of structural conservatism found in protein domain architecture and natural product inspired compound library design. Domains and proteins identified as being structurally similar in their ligand-sensing cores are grouped in a protein structure similarity cluster (PSSC). Natural products can be considered as evolutionary pre-validated ligands for multiple proteins and therefore natural products that are known to interact with one of the PSSC member proteins are selected as guiding structures for compound library synthesis. Application of this novel strategy for compound library design provided enhanced hit rates in small compound libraries for structurally similar proteins.     

Discovery of Anthelmintic Drug Targets and Drugs Using Chokepoints in Nematode Metabolic Pathways

Parasitic roundworm infections plague more than 2 billion people (1/3 of humanity) and cause drastic losses in crops and livestock. New anthelmintic drugs are urgently needed as new drug resistance and environmental concerns arise. A “chokepoint reaction” is defined as a reaction that either consumes a unique substrate or produces a unique product. A chokepoint analysis provides a systematic method of identifying novel potential drug targets. Chokepoint enzymes were identified in the genomes of 10 nematode species, and the intersection and union of all chokepoint enzymes were found. By studying and experimentally testing available compounds known to target proteins orthologous to nematode chokepoint proteins in public databases, this study uncovers features of chokepoints that make them successful drug targets. Chemogenomic screening was performed on drug-like compounds from public drug databases to find existing compounds that target homologs of nematode chokepoints. The compounds were prioritized based on chemical properties frequently found in successful drugs and were experimentally tested using Caenorhabditis elegans. Several drugs that are already known anthelmintic drugs and novel candidate targets were identified. Seven of the compounds were tested in Caenorhabditis elegans and three yielded a detrimental phenotype. One of these three drug-like compounds, Perhexiline, also yielded a deleterious effect in Haemonchus contortus and Onchocerca lienalis, two nematodes with divergent forms of parasitism. Perhexiline, known to affect the fatty acid oxidation pathway in mammals, caused a reduction in oxygen consumption rates in C. elegans and genome-wide gene expression profiles provided an additional confirmation of its mode of action. Computational modeling of Perhexiline and its target provided structural insights regarding its binding mode and specificity. Our lists of prioritized drug targets and drug-like compounds have potential to expedite the discovery of new anthelmintic drugs with broad-spectrum efficacy.     

Combinatorial computational approaches to identify tetracycline derivatives as flavivirus inhibitors

Limited structural information of drug targets, cellular toxicity possessed by lead compounds, and large amounts of potential leads are the major issues facing the design-oriented approach of discovering new leads. In an attempt to tackle these issues, we have developed a process of virtual screening based on the observation that conformational rearrangements of the dengue virus envelope protein are essential for the mediation of viral entry into host cells via membrane fusion. Screening was based solely on the structural information of the Dengue virus envelope protein and was focused on a target site that is presumably important for the conformational rearrangements necessary for viral entry. To circumvent the issue of lead compound toxicity, we performed screening based on molecular docking using structural databases of medical compounds. To enhance the identification of hits, we further categorized and selected candidates according to their novel structural characteristics. Finally, the selected candidates were subjected to a biological validation assay to assess inhibition of Dengue virus propagation in mammalian host cells using a plaque formation assay. Among the 10 compounds examined, rolitetracycline and doxycycline significantly inhibited plaque formation, demonstrating their inhibitory effect on dengue virus propagation. Both compounds were tetracycline derivatives with IC(50)s estimated to be 67.1 microM and 55.6 microM, respectively. Their docked conformations displayed common hydrophobic interactions with critical residues that affected membrane fusion during viral entry. These interactions will therefore position the tetracyclic ring moieties of both inhibitors to bind firmly to the target and, subsequently, disrupt conformational rearrangement and block viral entry. This process can be applied to other drug targets in which conformational rearrangement is critical to function.     

Identification of Small-Molecule Inhibitors against Meso-2, 6-Diaminopimelate Dehydrogenase from Porphyromonas gingivalis

Species-specific antimicrobial therapy has the potential to combat the increasing threat of antibiotic resistance and alteration of the human microbiome. We therefore set out to demonstrate the beginning of a pathogen-selective drug discovery method using the periodontal pathogen Porphyromonas gingivalis as a model. Through our knowledge of metabolic networks and essential genes we identified a “druggable” essential target, meso-diaminopimelate dehydrogenase, which is found in a limited number of species. We adopted a high-throughput virtual screen method on the ZINC chemical library to select a group of potential small-molecule inhibitors. Meso-diaminopimelate dehydrogenase from P. gingivalis was first expressed and purified in Escherichia coli then characterized for enzymatic inhibitor screening studies. Several inhibitors with similar structural scaffolds containing a sulfonamide core and aromatic substituents showed dose-dependent inhibition. These compounds were further assayed showing reasonable whole-cell activity and the inhibition mechanism was determined. We conclude that the establishment of this target and screening strategy provides a model for the future development of new antimicrobials.     

Modeling the CRL4A ligase complex to predict target protein ubiquitination induced by cereblon-recruiting PROTACs

PROteolysis TArgeting Chimeras (PROTACs) are hetero-bifunctional small molecules that can simultaneously recruit target proteins and E3 ligases to form a ternary complex, promoting target protein ubiquitination and degradation via the Ubiquitin-Proteasome System (UPS). PROTACs have gained increasing attention in recent years due to certain advantages over traditional therapeutic modalities and enabling targeting of previously “undruggable” proteins. To better understand the mechanism of PROTAC-induced Target Protein Degradation (TPD), several computational approaches have recently been developed to study and predict ternary complex formation. However, mounting evidence suggests that ubiquitination can also be a rate-limiting step in PROTAC-induced TPD. Here, we propose a structure-based computational approach to predict target protein ubiquitination induced by cereblon (CRBN)-based PROTACs by leveraging available structural information of the CRL4A ligase complex (CRBN/DDB1/CUL4A/Rbx1/NEDD8/E2/Ub). We generated ternary complex ensembles with Rosetta, modeled multiple CRL4A ligase complex conformations, and predicted ubiquitination efficiency by separating the ternary ensemble into productive and unproductive complexes based on the proximity of the ubiquitin to accessible lysines on the target protein. We validated our CRL4A ligase complex models with published ternary complex structures and additionally employed our modeling workflow to predict ubiquitination efficiencies and sites of a series of cyclin-dependent kinases (CDKs) after treatment with TL12–186, a pan-kinase PROTAC. Our predictions are consistent with CDK ubiquitination and site-directed mutagenesis of specific CDK lysine residues as measured using a NanoBRET ubiquitination assay in HEK293 cells. This work structurally links PROTAC-induced ternary formation and ubiquitination, representing an important step toward prediction of target “degradability.”