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.

Current and Future Roles of Artificial Intelligence in Medicinal Chemistry Synthesis.

Artificial intelligence and machine learning have demonstrated their potential role in predictive chemistry and synthetic planning of small molecules; there are at least a few reports of companies employing in silico synthetic planning into their overall approach to accessing target molecules. A data-driven synthesis planning program is one component being developed and evaluated by the Machine Learning for Pharmaceutical Discovery and Synthesis (MLPDS) consortium, comprising MIT and 13 chemical and pharmaceutical company members. Together, we wrote this perspective to share how we think predictive models can be integrated into medicinal chemistry synthesis workflows, how they are currently used within MLPDS member companies, and the outlook for this field.

UniLectin3D, a database of carbohydrate binding proteins with curated information on 3D structures and interacting ligands.

Lectins, and related receptors such as adhesins and toxins, are glycan-binding proteins from all origins that decipher the glycocode, i.e. the structural information encoded in the conformation of complex carbohydrates present on the surface of all cells. Lectins are still poorly classified and annotated, but since their functions are based on ligand recognition, their 3D-structures provide a solid foundation for characterization. UniLectin3D is a curated database that classifies lectins on origin and fold, with cross-links to literature, other databases in glycosciences and functional data such as known specificity. The database provides detailed information on lectins, their bound glycan ligands, and features their interactions using the Protein-Ligand Interaction Profiler (PLIP) server. Special care was devoted to the description of the bound glycan ligands with the use of simple graphical representation and numerical format for cross-linking to other databases in glycoscience. We conceived the design of the database architecture and the navigation tools to account for all organisms, as well as to search for oligosaccharide epitopes complexed within specified binding sites. UniLectin3D is accessible at https://www.unilectin.eu/unilectin3D.

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.

MAGPIE: Simplifying access and execution of computational models in the life sciences.

Over the past decades, quantitative methods linking theory and observation became increasingly important in many areas of life science. Subsequently, a large number of mathematical and computational models has been developed. The BioModels database alone lists more than 140,000 Systems Biology Markup Language (SBML) models. However, while the exchange within specific model classes has been supported by standardisation and database efforts, the generic application and especially the re-use of models is still limited by practical issues such as easy and straight forward model execution. MAGPIE, a Modeling and Analysis Generic Platform with Integrated Evaluation, closes this gap by providing a software platform for both, publishing and executing computational models without restrictions on the programming language, thereby combining a maximum on flexibility for programmers with easy handling for non-technical users. MAGPIE goes beyond classical SBML platforms by including all models, independent of the underlying programming language, ranging from simple script models to complex data integration and computations. We demonstrate the versatility of MAGPIE using four prototypic example cases. We also outline the potential of MAGPIE to improve transparency and reproducibility of computational models in life sciences. A demo server is available at magpie.imb.medizin.tu-dresden.de.

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.

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.