AI-Powered Virtual Screening of Large Compound Libraries Leads to the Discovery of Novel Inhibitors of Sirtuin-1.

The discovery of new scaffolds and chemotypes via high-throughput screening is tedious and resource intensive. Yet, there are millions of small molecules commercially available, rendering comprehensive in vitro tests intractable. We show how smart algorithms reduce large screening collections to target-specific sets of just a few hundred small molecules, allowing for a much faster and more cost-effective hit discovery process. We showcase the application of this virtual screening strategy by preselecting 434 compounds for Sirtuin-1 inhibition from a library of 2.6 million compounds, corresponding to 0.02% of the original library. Multistage in vitro validation ultimately confirmed nine chemically novel inhibitors. When compared to a competitive benchmark study for Sirtuin-1, our method shows a 12-fold higher hit rate. The results demonstrate how AI-driven preselection from large screening libraries allows for a massive reduction in the number of small molecules to be tested in vitro while still retaining a large number of hits.

An interaction-based drug discovery screen explains known SARS-CoV-2 inhibitors and predicts new compound scaffolds.

The recent outbreak of the COVID-19 pandemic caused by severe acute respiratory syndrome-Coronavirus-2 (SARS-CoV-2) has shown the necessity for fast and broad drug discovery methods to enable us to react quickly to novel and highly infectious diseases. A well-known SARS-CoV-2 target is the viral main 3-chymotrypsin-like cysteine protease (Mpro), known to control coronavirus replication, which is essential for the viral life cycle. Here, we applied an interaction-based drug repositioning algorithm on all protein-compound complexes available in the protein database (PDB) to identify Mpro inhibitors and potential novel compound scaffolds against SARS-CoV-2. The screen revealed a heterogeneous set of 692 potential Mpro inhibitors containing known ones such as Dasatinib, Amodiaquine, and Flavin mononucleotide, as well as so far untested chemical scaffolds. In a follow-up evaluation, we used publicly available data published almost two years after the screen to validate our results. In total, we are able to validate 17% of the top 100 predictions with publicly available data and can furthermore show that predicted compounds do cover scaffolds that are yet not associated with Mpro. Finally, we detected a potentially important binding pattern consisting of 3 hydrogen bonds with hydrogen donors of an oxyanion hole within the active side of Mpro. Overall, these results give hope that we will be better prepared for future pandemics and that drug development will become more efficient in the upcoming years.

Structural binding site comparisons reveal Crizotinib as a novel LRRK2 inhibitor.

Mutations in leucine-rich repeat kinase 2 (LRRK2) are a frequent cause of autosomal dominant Parkinson's disease (PD) and have been associated with familial and sporadic PD. Reducing the kinase activity of LRRK2 is a promising therapeutic strategy since pathogenic mutations increase the kinase activity. Several small-molecule LRRK2 inhibitors are currently under investigation for the treatment of PD. However, drug discovery and development are always accompanied by high costs and a risk of late failure. The use of already approved drugs for a new indication, which is known as drug repositioning, can reduce the cost and risk. In this study, we applied a structure-based drug repositioning approach to identify new LRRK2 inhibitors that are already approved for a different indication. In a large-scale structure-based screening, we compared the protein-ligand interaction patterns of known LRRK2 inhibitors with protein-ligand complexes in the PDB. The screening yielded 6 drug repositioning candidates. Two of these candidates, Sunitinib and Crizotinib, demonstrated an inhibition potency (IC50) and binding affinity (Kd) in the nanomolar to micromolar range. While Sunitinib has already been known to inhibit LRRK2, Crizotinib is a novel LRRK2 binder. Our results underscore the potential of structure-based methods for drug discovery and development. In light of the recent breakthroughs in cryo-electron microscopy and structure prediction, we believe that structure-based approaches like ours will grow in importance.

PLIP 2021: expanding the scope of the protein-ligand interaction profiler to DNA and RNA.

With the growth of protein structure data, the analysis of molecular interactions between ligands and their target molecules is gaining importance. PLIP, the protein-ligand interaction profiler, detects and visualises these interactions and provides data in formats suitable for further processing. PLIP has proven very successful in applications ranging from the characterisation of docking experiments to the assessment of novel ligand-protein complexes. Besides ligand-protein interactions, interactions with DNA and RNA play a vital role in many applications, such as drugs targeting DNA or RNA-binding proteins. To date, over 7% of all 3D structures in the Protein Data Bank include DNA or RNA. Therefore, we extended PLIP to encompass these important molecules. We demonstrate the power of this extension with examples of a cancer drug binding to a DNA target, and an RNA-protein complex central to a neurological disease. PLIP is available online at https://plip-tool.biotec.tu-dresden.de and as open source code. So far, the engine has served over a million queries and the source code has been downloaded several thousand times.

Repositioned Drugs for Chagas Disease Unveiled via Structure-Based Drug Repositioning.

Chagas disease, caused by the parasite Trypanosoma cruzi, affects millions of people in South America. The current treatments are limited, have severe side effects, and are only partially effective. Drug repositioning, defined as finding new indications for already approved drugs, has the potential to provide new therapeutic options for Chagas. In this work, we conducted a structure-based drug repositioning approach with over 130,000 3D protein structures to identify drugs that bind therapeutic Chagas targets and thus represent potential new Chagas treatments. The screening yielded over 500 molecules as hits, out of which 38 drugs were prioritized following a rigorous filtering process. About half of the latter were already known to have trypanocidal activity, while the others are novel to Chagas disease. Three of the new drug candidates-ciprofloxacin, naproxen, and folic acid-showed a growth inhibitory activity in the micromolar range when tested ex vivo on T. cruzi trypomastigotes, validating the prediction. We show that our drug repositioning approach is able to pinpoint relevant drug candidates at a fraction of the time and cost of a conventional screening. Furthermore, our results demonstrate the power and potential of structure-based drug repositioning in the context of neglected tropical diseases where the pharmaceutical industry has little financial interest in the development of new drugs.

The structural basis of the genetic code: amino acid recognition by aminoacyl-tRNA synthetases.

Storage and directed transfer of information is the key requirement for the development of life. Yet any information stored on our genes is useless without its correct interpretation. The genetic code defines the rule set to decode this information. Aminoacyl-tRNA synthetases are at the heart of this process. We extensively characterize how these enzymes distinguish all natural amino acids based on the computational analysis of crystallographic structure data. The results of this meta-analysis show that the correct read-out of genetic information is a delicate interplay between the composition of the binding site, non-covalent interactions, error correction mechanisms, and steric effects.

Computational Drug Repositioning for Chagas Disease Using Protein-Ligand Interaction Profiling.

Chagas disease, caused by Trypanosoma cruzi (T. cruzi), affects nearly eight million people worldwide. There are currently only limited treatment options, which cause several side effects and have drug resistance. Thus, there is a great need for a novel, improved Chagas treatment. Bifunctional enzyme dihydrofolate reductase-thymidylate synthase (DHFR-TS) has emerged as a promising pharmacological target. Moreover, some human dihydrofolate reductase (HsDHFR) inhibitors such as trimetrexate also inhibit T. cruzi DHFR-TS (TcDHFR-TS). These compounds serve as a starting point and a reference in a screening campaign to search for new TcDHFR-TS inhibitors. In this paper, a novel virtual screening approach was developed that combines classical docking with protein-ligand interaction profiling to identify drug repositioning opportunities against T. cruzi infection. In this approach, some food and drug administration (FDA)-approved drugs that were predicted to bind with high affinity to TcDHFR-TS and whose predicted molecular interactions are conserved among known inhibitors were selected. Overall, ten putative TcDHFR-TS inhibitors were identified. These exhibited a similar interaction profile and a higher computed binding affinity, compared to trimetrexate. Nilotinib, glipizide, glyburide and gliquidone were tested on T. cruzi epimastigotes and showed growth inhibitory activity in the micromolar range. Therefore, these compounds could lead to the development of new treatment options for Chagas disease.

Structure-based drug repositioning explains ibrutinib as VEGFR2 inhibitor.

Many drugs are promiscuous and bind to multiple targets. On the one hand, these targets may be linked to unwanted side effects, but on the other, they may achieve a combined desired effect (polypharmacology) or represent multiple diseases (drug repositioning). With the growth of 3D structures of drug-target complexes, it is today possible to study drug promiscuity at the structural level and to screen vast amounts of drug-target interactions to predict side effects, polypharmacological potential, and repositioning opportunities. Here, we pursue such an approach to identify drugs inactivating B-cells, whose dysregulation can function as a driver of autoimmune diseases. Screening over 500 kinases, we identified 22 candidate targets, whose knock out impeded the activation of B-cells. Among these 22 is the gene KDR, whose gene product VEGFR2 is a prominent cancer target with anti-VEGFR2 drugs on the market for over a decade. The main result of this paper is that structure-based drug repositioning for the identified kinase targets identified the cancer drug ibrutinib as micromolar VEGFR2 inhibitor with a very high therapeutic index in B-cell inactivation. These findings prove that ibrutinib is not only acting on the Bruton's tyrosine kinase BTK, against which it was designed. Instead, it may be a polypharmacological drug, which additionally targets angiogenesis via inhibition of VEGFR2. Therefore ibrutinib carries potential to treat other VEGFR2 associated disease. Structure-based drug repositioning explains ibrutinib's anti VEGFR2 action through the conservation of a specific pattern of interactions of the drug with BTK and VEGFR2. Overall, structure-based drug repositioning was able to predict these findings at a fraction of the time and cost of a conventional screen.

Backbone Brackets and Arginine Tweezers delineate Class I and Class II aminoacyl tRNA synthetases.

The origin of the machinery that realizes protein biosynthesis in all organisms is still unclear. One key component of this machinery are aminoacyl tRNA synthetases (aaRS), which ligate tRNAs to amino acids while consuming ATP. Sequence analyses revealed that these enzymes can be divided into two complementary classes. Both classes differ significantly on a sequence and structural level, feature different reaction mechanisms, and occur in diverse oligomerization states. The one unifying aspect of both classes is their function of binding ATP. We identified Backbone Brackets and Arginine Tweezers as most compact ATP binding motifs characteristic for each Class. Geometric analysis shows a structural rearrangement of the Backbone Brackets upon ATP binding, indicating a general mechanism of all Class I structures. Regarding the origin of aaRS, the Rodin-Ohno hypothesis states that the peculiar nature of the two aaRS classes is the result of their primordial forms, called Protozymes, being encoded on opposite strands of the same gene. Backbone Brackets and Arginine Tweezers were traced back to the proposed Protozymes and their more efficient successors, the Urzymes. Both structural motifs can be observed as pairs of residues in contemporary structures and it seems that the time of their addition, indicated by their placement in the ancient aaRS, coincides with the evolutionary trace of Proto- and Urzymes.

Pioneering topological methods for network-based drug-target prediction by exploiting a brain-network self-organization theory.

The bipartite network representation of the drug-target interactions (DTIs) in a biosystem enhances understanding of the drugs' multifaceted action modes, suggests therapeutic switching for approved drugs and unveils possible side effects. As experimental testing of DTIs is costly and time-consuming, computational predictors are of great aid. Here, for the first time, state-of-the-art DTI supervised predictors custom-made in network biology were compared-using standard and innovative validation frameworks-with unsupervised pure topological-based models designed for general-purpose link prediction in bipartite networks. Surprisingly, our results show that the bipartite topology alone, if adequately exploited by means of the recently proposed local-community-paradigm (LCP) theory-initially detected in brain-network topological self-organization and afterwards generalized to any complex network-is able to suggest highly reliable predictions, with comparable performance with the state-of-the-art-supervised methods that exploit additional (non-topological, for instance biochemical) DTI knowledge. Furthermore, a detailed analysis of the novel predictions revealed that each class of methods prioritizes distinct true interactions; hence, combining methodologies based on diverse principles represents a promising strategy to improve drug-target discovery. To conclude, this study promotes the power of bio-inspired computing, demonstrating that simple unsupervised rules inspired by principles of topological self-organization and adaptiveness arising during learning in living intelligent systems (like the brain) can efficiently equal perform complicated algorithms based on advanced, supervised and knowledge-based engineering.

From malaria to cancer: Computational drug repositioning of amodiaquine using PLIP interaction patterns.

Drug repositioning identifies new indications for known drugs. Here we report repositioning of the malaria drug amodiaquine as a potential anti-cancer agent. While most repositioning efforts emerge through serendipity, we have devised a computational approach, which exploits interaction patterns shared between compounds. As a test case, we took the anti-viral drug brivudine (BVDU), which also has anti-cancer activity, and defined ten interaction patterns using our tool PLIP. These patterns characterise BVDU's interaction with its target s. Using PLIP we performed an in silico screen of all structural data currently available and identified the FDA approved malaria drug amodiaquine as a promising repositioning candidate. We validated our prediction by showing that amodiaquine suppresses chemoresistance in a multiple myeloma cancer cell line by inhibiting the chaperone function of the cancer target Hsp27. This work proves that PLIP interaction patterns are viable tools for computational repositioning and can provide search query information from a given drug and its target to identify structurally unrelated candidates, including drugs approved by the FDA, with a known safety and pharmacology profile. This approach has the potential to reduce costs and risks in drug development by predicting novel indications for known drugs and drug candidates.

Discovery of Mycobacterium tuberculosis InhA Inhibitors by Binding Sites Comparison and Ligands Prediction.

Drug discovery is usually focused on a single protein target; in this process, existing compounds that bind to related proteins are often ignored. We describe ProBiS plugin, extension of our earlier ProBiS-ligands approach, which for a given protein structure allows prediction of its binding sites and, for each binding site, the ligands from similar binding sites in the Protein Data Bank. We developed a new database of precalculated binding site comparisons of about 290000 proteins to allow fast prediction of binding sites in existing proteins. The plugin enables advanced viewing of predicted binding sites, ligands' poses, and their interactions in three-dimensional graphics. Using the InhA query protein, an enoyl reductase enzyme in the Mycobacterium tuberculosis fatty acid biosynthesis pathway, we predicted its possible ligands and assessed their inhibitory activity experimentally. This resulted in three previously unrecognized inhibitors with novel scaffolds, demonstrating the plugin's utility in the early drug discovery process.

New HSP27 inhibitors efficiently suppress drug resistance development in cancer cells.

Drug resistance is an important open problem in cancer treatment. In recent years, the heat shock protein HSP27 (HSPB1) was identified as a key player driving resistance development. HSP27 is overexpressed in many cancer types and influences cellular processes such as apoptosis, DNA repair, recombination, and formation of metastases. As a result cancer cells are able to suppress apoptosis and develop resistance to cytostatic drugs. To identify HSP27 inhibitors we follow a novel computational drug repositioning approach. We exploit a similarity between a predicted HSP27 binding site to a viral thymidine kinase to generate lead inhibitors for HSP27. Six of these leads were verified experimentally. They bind HSP27 and down-regulate its chaperone activity. Most importantly, all six compounds inhibit development of drug resistance in cellular assays. One of the leads - chlorpromazine - is an antipsychotic, which has a positive effect on survival time in human breast cancer. In summary, we make two important contributions: First, we put forward six novel leads, which inhibit HSP27 and tackle drug resistance. Second, we demonstrate the power of computational drug repositioning.

Computational Drug Repositioning by Target Hopping: A Use Case in Chagas Disease.

BACKGROUND: Drug repositioning aims to identify novel indications for existing drugs. One approach to repositioning exploits shared binding sites between the drug targets and other proteins. Here, we review the principle and algorithms of such target hopping and illustrate them in Chagas disease, an in Latin America widely spread, but neglected disease. CONCLUSION: We demonstrate how target hopping recovers known treatments for Chagas disease and predicts novel drugs, such as the antiviral foscarnet, which we predict to target Farnesyl Pyrophosphate Synthase in Trypanosoma cruzi, the causative agent of Chagas disease.

PLIP: fully automated protein-ligand interaction profiler.

The characterization of interactions in protein-ligand complexes is essential for research in structural bioinformatics, drug discovery and biology. However, comprehensive tools are not freely available to the research community. Here, we present the protein-ligand interaction profiler (PLIP), a novel web service for fully automated detection and visualization of relevant non-covalent protein-ligand contacts in 3D structures, freely available at projects.biotec.tu-dresden.de/plip-web. The input is either a Protein Data Bank structure, a protein or ligand name, or a custom protein-ligand complex (e.g. from docking). In contrast to other tools, the rule-based PLIP algorithm does not require any structure preparation. It returns a list of detected interactions on single atom level, covering seven interaction types (hydrogen bonds, hydrophobic contacts, pi-stacking, pi-cation interactions, salt bridges, water bridges and halogen bonds). PLIP stands out by offering publication-ready images, PyMOL session files to generate custom images and parsable result files to facilitate successive data processing. The full python source code is available for download on the website. PLIP's command-line mode allows for high-throughput interaction profiling.

Polypharmacology rescored: protein-ligand interaction profiles for remote binding site similarity assessment.

Detection of remote binding site similarity in proteins plays an important role for drug repositioning and off-target effect prediction. Various non-covalent interactions such as hydrogen bonds and van-der-Waals forces drive ligands' molecular recognition by binding sites in proteins. The increasing amount of available structures of protein-small molecule complexes enabled the development of comparative approaches. Several methods have been developed to characterize and compare protein-ligand interaction patterns. Usually implemented as fingerprints, these are mainly used for post processing docking scores and (off-)target prediction. In the latter application, interaction profiles detect similarities in the bound interactions of different ligands and thus identify essential interactions between a protein and its small molecule ligands. Interaction pattern similarity correlates with binding site similarity and is thus contributing to a higher precision in binding site similarity assessment of proteins with distinct global structure. This renders it valuable for existing drug repositioning approaches in structural bioinformatics. Current methods to characterize and compare structure-based interaction patterns - both for protein-small-molecule and protein-protein interactions - as well as their potential in target prediction will be reviewed in this article. The question of how the set of interaction types, flexibility or water-mediated interactions, influence the comparison of interaction patterns will be discussed. Due to the wealth of protein-ligand structures available today, predicted targets can be ranked by comparing their ligand interaction pattern to patterns of the known target. Such knowledge-based methods offer high precision in comparison to methods comparing whole binding sites based on shape and amino acid physicochemical similarity.

Drug Promiscuity in PDB: Protein Binding Site Similarity Is Key.

Drug repositioning applies established drugs to new disease indications with increasing success. A pre-requisite for drug repurposing is drug promiscuity (polypharmacology) - a drug's ability to bind to several targets. There is a long standing debate on the reasons for drug promiscuity. Based on large compound screens, hydrophobicity and molecular weight have been suggested as key reasons. However, the results are sometimes contradictory and leave space for further analysis. Protein structures offer a structural dimension to explain promiscuity: Can a drug bind multiple targets because the drug is flexible or because the targets are structurally similar or even share similar binding sites? We present a systematic study of drug promiscuity based on structural data of PDB target proteins with a set of 164 promiscuous drugs. We show that there is no correlation between the degree of promiscuity and ligand properties such as hydrophobicity or molecular weight but a weak correlation to conformational flexibility. However, we do find a correlation between promiscuity and structural similarity as well as binding site similarity of protein targets. In particular, 71% of the drugs have at least two targets with similar binding sites. In order to overcome issues in detection of remotely similar binding sites, we employed a score for binding site similarity: LigandRMSD measures the similarity of the aligned ligands and uncovers remote local similarities in proteins. It can be applied to arbitrary structural binding site alignments. Three representative examples, namely the anti-cancer drug methotrexate, the natural product quercetin and the anti-diabetic drug acarbose are discussed in detail. Our findings suggest that global structural and binding site similarity play a more important role to explain the observed drug promiscuity in the PDB than physicochemical drug properties like hydrophobicity or molecular weight. Additionally, we find ligand flexibility to have a minor influence.

Drug repositioning through incomplete bi-cliques in an integrated drug-target-disease network.

Recently, there has been much interest in gene-disease networks and polypharmacology as a basis for drug repositioning. Here, we integrate data from structural and chemical databases to create a drug-target-disease network for 147 promiscuous drugs, their 553 protein targets, and 44 disease indications. Visualizing and analyzing such complex networks is still an open problem. We approach it by mining the network for network motifs of bi-cliques. In our case, a bi-clique is a subnetwork in which every drug is linked to every target and disease. Since the data are incomplete, we identify incomplete bi-cliques, whose completion introduces novel, predicted links from drugs to targets and diseases. We demonstrate the power of this approach by repositioning cardiovascular drugs to parasitic diseases, by predicting the cancer-related kinase PIK3CG as a novel target of resveratrol, and by identifying for five drugs a shared binding site in four serine proteases and novel links to cancer, cardiovascular, and parasitic diseases.

Old friends in new guise: repositioning of known drugs with structural bioinformatics.

Developing a drug de novo is a laborious and costly endeavor. Thus, the repositioning of already approved drugs for the treatment of new diseases is promising and valuable. One computational approach to repositioning exploits the structural similarity of binding sites of known and new targets. Here, we review computational methods to represent and align binding sites. We review available tools, present success stories and discuss limits of the approach.

Chemoinformatic analysis of biologically active macrocycles.

In this review, macrocycles are defined as organic molecules containing at least one non-bridged cycle of 12 or more covalently connected atoms. Common statistical aspects as well as structure activity relationships (SAR) of macrocycles based on chemoinformatics methods are discussed. Can macrocyclic structural features be linked to activities against cancer cells, viruses like HIV, or bacteria? General challenges and problems in using chemoinformatics for more detailed analyses of macrocycles based on large compound databases are outlined.