Impact of Binding Site Comparisons on Medicinal Chemistry and Rational Molecular Design

Subpocket Similarity-Based Hit Identification for Challenging Targets: Application to the WDR Domain of LRRK2

We herewith applied a priori a generic hit identification method (POEM) for difficult targets of known three-dimensional structure, relying on the simple knowledge of physicochemical and topological properties of a user-selected cavity. Searching for local similarity to a set of fragment-bound protein microenvironments of known structure, a point cloud registration algorithm is first applied to align known subpockets to the target cavity. The resulting alignment then permits us to directly pose the corresponding seed fragments in a target cavity space not typically amenable to classical docking approaches. Last, linking potentially connectable atoms by a deep generative linker enables full ligand enumeration. When applied to the WD40 repeat (WDR) central cavity of leucine-rich repeat kinase 2 (LRRK2), an unprecedented binding site, POEM was able to quickly propose 94 potential hits, five of which were subsequently confirmed to bind in vitro to LRRK2-WDR.

User-centric design of a 3D search interface for protein-ligand complexes

In this work, we present the frontend of GeoMine and showcase its application, focusing on the new features of its latest version. GeoMine is a search engine for ligand-bound and predicted empty binding sites in the Protein Data Bank. In addition to its basic text-based search functionalities, GeoMine offers a geometric query type for searching binding sites with a specific relative spatial arrangement of chemical features such as heavy atoms and intermolecular interactions. In contrast to a text search that requires simple and easy-to-formulate user input, a 3D input is more complex, and its specification can be challenging for users. GeoMine's new version aims to address this issue from the graphical user interface perspective by introducing an additional visualization concept and a new query template type. In its latest version, GeoMine extends its query-building capabilities primarily through input formulation in 2D. The 2D editor is fully synchronized with GeoMine's 3D editor and provides the same functionality. It enables template-free query generation and template-based query selection directly in 2D pose diagrams. In addition, the query generation with the 3D editor now supports predicted empty binding sites for AlphaFold structures as query templates. GeoMine is freely accessible on the ProteinsPlus web server ( https://proteins.plus ).

An optimized relational database for querying structural patterns in proteins

A database is an essential component in almost any software system, and its creation involves more than just data modeling and schema design. It also includes query optimization and tuning. This paper focuses on a web system called GSP4PDB, which is used for searching structural patterns in proteins. The system utilizes a normalized relational database, which has proven to be inefficient even for simple queries. This article discusses the optimization of the GSP4PDB database by implementing two techniques: denormalization and indexing. The empirical evaluation described in the article shows that combining these techniques enhances the efficiency of the database when querying both real and artificial graph-based structural patterns.

Integrated Molecular Modeling and Machine Learning for Drug Design

Modern therapeutic development often involves several stages that are interconnected, and multiple iterations are usually required to bring a new drug to the market. Computational approaches have increasingly become an indispensable part of helping reduce the time and cost of the research and development of new drugs. In this Perspective, we summarize our recent efforts on integrating molecular modeling and machine learning to develop computational tools for modulator design, including a pocket-guided rational design approach based on AlphaSpace to target protein-protein interactions, delta machine learning scoring functions for protein-ligand docking as well as virtual screening, and state-of-the-art deep learning models to predict calculated and experimental molecular properties based on molecular mechanics optimized geometries. Meanwhile, we discuss remaining challenges and promising directions for further development and use a retrospective example of FDA approved kinase inhibitor Erlotinib to demonstrate the use of these newly developed computational tools.

Searching similar local 3D micro-environments in protein structure databases with MicroMiner

The available protein structure data are rapidly increasing. Within these structures, numerous local structural sites depict the details characterizing structure and function. However, searching and analyzing these sites extensively and at scale poses a challenge. We present a new method to search local sites in protein structure databases using residue-defined local 3D micro-environments. We implemented the method in a new tool called MicroMiner and demonstrate the capabilities of residue micro-environment search on the example of structural mutation analysis. Usually, experimental structures for both the wild-type and the mutant are unavailable for comparison. With MicroMiner, we extracted $>255 \times 10^{6}$ amino acid pairs in protein structures from the PDB, exemplifying single mutations' local structural changes for single chains and $>45 \times 10^{6}$ pairs for protein-protein interfaces. We further annotate existing data sets of experimentally measured mutation effects, like $\Delta \Delta G$ measurements, with the extracted structure pairs to combine the mutation effect measurement with the structural change upon mutation. In addition, we show how MicroMiner can bridge the gap between mutation analysis and structure-based drug design tools. MicroMiner is available as a command line tool and interactively on the https://proteins.plus/ webserver.

TWN-RENCOD: A novel method for protein binding site comparison

Several diverse proteins possess similar binding sites. Protein binding site comparison provides valuable insights for the drug discovery and development. Binding site similarities are useful in understanding polypharmacology, identifying potential off-targets and repurposing of known drugs. Many binding site analysis and comparison methods are available today, however, these methods may not be adequate to explain variation in the activity of a drug or a small molecule against a number of similar proteins. Water molecules surrounding the protein surface contribute to structure and function of proteins. Water molecules form diverse types of hydrogen-bonded cyclic water-ring networks known as topological water networks (TWNs). Analysis of TWNs in binding site of proteins may improve understanding of the characteristics of binding sites. We propose TWN-based residue encoding (TWN-RENCOD), a novel binding site comparison method which compares the aqueous environment in binding sites of similar proteins. As compared to other existing methods, results obtained using our method correlated better with differences in wide range of activity of a known drug (Sunitinib) against nine different protein kinases (KIT, PDGFRA, VEGFR2, PHKG2, ITK, HPK1, MST3, PAK6 and CDK2).

A high-throughput effector screen identifies a novel small molecule scaffold for inhibition of ten-eleven translocation dioxygenase 2

Ten-eleven translocation dioxygenases (TETs) are the erasers of 5-methylcytosine (mC), the central epigenetic regulator of mammalian DNA. TETs convert mC to three oxidized derivatives with unique physicochemical properties and inherent regulatory potential, and it initializes active demethylation by the base excision repair pathway. Potent small molecule inhibitors would be useful tools to study TET functions by conditional control. To facilitate the discovery of such tools, we here report a high-throughput screening pipeline and its application to screen and validate 31.5k compounds for inhibition of TET2. Using a homogenous fluorescence assay, we discover a novel quinoline-based scaffold that we further validate with an orthogonal semi-high throughput MALDI-MS assay for direct monitoring of substrate turnover. Structure-activity relationship (SAR) studies involving >20 derivatives of this scaffold led to the identification of optimized inhibitors, and together with computational studies suggested a plausible model for its mode of action.

Protein binding sites for drug design

Drug development is a lengthy and challenging process that can be accelerated at early stages by new mathematical approaches and modern computers. To address this important issue, we are developing new mathematical solutions for the detection and characterization of protein binding sites that are important for new drug development. In this review, we present algorithms based on graph theory combined with molecular dynamics simulations that we have developed for studying biological target proteins to provide important data for optimizing the early stages of new drug development. A particular focus is the development of new protein binding site prediction algorithms (ProBiS) and new web tools for modeling pharmaceutically interesting molecules-ProBiS Tools (algorithm, database, web server), which have evolved into a full-fledged graphical tool for studying proteins in the proteome. ProBiS differs from other structural algorithms in that it can align proteins with different folds without prior knowledge of the binding sites. It allows detection of similar binding sites and can predict molecular ligands of various types of pharmaceutical interest that could be advanced to drugs to treat a disease, based on the entire Protein Data Bank (PDB) and AlphaFold database, including proteins not yet in the PDB. All ProBiS Tools are freely available to the academic community at http://insilab.org and https://probis.nih.gov.

Classification of Protein-Binding Sites Using a Spherical Convolutional Neural Network

The analysis and comparison of protein-binding sites aid various applications in the drug discovery process, e.g., hit finding, drug repurposing, and polypharmacology. Classification of binding sites has been a hot topic for the past 30 years, and many different methods have been published. The rapid development of machine learning computational algorithms, coupled with the large volume of publicly available protein-ligand 3D structures, makes it possible to apply deep learning techniques in binding site comparison. Our method uses a cutting-edge spherical convolutional neural network based on the DeepSphere architecture to learn global representations of protein-binding sites. The model was trained on TOUGH-C1 and TOUGH-M1 data and validated with the ProSPECCTs datasets. Our results show that our model can (1) perform well in protein-binding site similarity and classification tasks and (2) learn and separate the physicochemical properties of binding sites. Lastly, we tested the model on a set of kinases, where the results show that it is able to cluster the different kinase subfamilies effectively. This example demonstrates the method's promise for lead hopping within or outside a protein target, directly based on binding site information.

ProBiS-Fold Approach for Annotation of Human Structures from the AlphaFold Database with No Corresponding Structure in the PDB to Discover New Druggable Binding Sites

ProBiS (Protein Binding Sites), a local structure-based comparison algorithm, is used in the new ProBiS-Fold web server to annotate human structures from the AlphaFold database without a corresponding structure in the Protein Data Bank (PDB) to discover new druggable binding sites. The ProBiS algorithm is used to compare each query protein structure predicted by the AlphaFold approach with the protein structures from the PDB to identify similarities between known binding sites found in the PDB and yet unknown binding sites in the AlphaFold database. Ligands bound in these identified similar PDB sites are then transferred to each query protein from the AlphaFold database, and binding sites are identified as ligand clusters on an AlphaFold protein. Small molecule binding sites are assigned druggability scores based on the similarity of their ligands to known drugs, allowing them to be ranked according to their perceived and actual potential for drug development. ProBiS-Fold provides interactive and downloadable binding sites for the entire human structural proteome, including more than 3000 new druggable binding sites that have no corresponding structure in the PDB, taking into account AlphaFold's model quality, to enable protein structure-function relationship studies and pharmaceutical drug discovery research. The web server is freely accessible to academic users at http://probis-fold.insilab.org.

Structural fractal analysis of the active sites of acetylcholinesterase from various organisms

Acetylcholinesterase (AChE) is the object of many studies due to the fact that it plays an important role in the vital activity of organisms. In particular, when new AChE inhibitors are developed, much attention is paid to the specificity of their action. One of the approaches used to study the specificity is to compare AChE taken from various organisms. In this work, crystallographic data are used to investigate the active sites of AChE (ASAs) in the free (uncomplexed) state for the following five organisms: Homo sapiens (HS), Mus musculus (MM), Torpedo californica (TC), Electrophorus electricus (EE), and Drosophila melanogaster (DM). The structural fractal analysis (SFA) proposed by us earlier is used as a research method. This method is based on the calculation and comparison of the fractal dimensions of molecular structures. SFA demonstrates that there are no significant structural differences between the active sites of human AChE and other AChEs. However, differences are found for the MM/EE pair. Further analysis of individual AARs has revealed two different areas of active sites. Ser203, Trp236, Phe338, and Tyr341 are found to belong to a variable region, and the remaining AARs belong to a conservative region of the ASAs. The fraction of "variability" is low, 0.8%.

Estimating the Similarity between Protein Pockets

With the exponential increase in publicly available protein structures, the comparison of protein binding sites naturally emerged as a scientific topic to explain observations or generate hypotheses for ligand design, notably to predict ligand selectivity for on- and off-targets, explain polypharmacology, and design target-focused libraries. The current review summarizes the state-of-the-art computational methods applied to pocket detection and comparison as well as structural druggability estimates. The major strengths and weaknesses of current pocket descriptors, alignment methods, and similarity search algorithms are presented. Lastly, an exhaustive survey of both retrospective and prospective applications in diverse medicinal chemistry scenarios illustrates the capability of the existing methods and the hurdle that still needs to be overcome for more accurate predictions.

GraphSite: Ligand Binding Site Classification with Deep Graph Learning

The binding of small organic molecules to protein targets is fundamental to a wide array of cellular functions. It is also routinely exploited to develop new therapeutic strategies against a variety of diseases. On that account, the ability to effectively detect and classify ligand binding sites in proteins is of paramount importance to modern structure-based drug discovery. These complex and non-trivial tasks require sophisticated algorithms from the field of artificial intelligence to achieve a high prediction accuracy. In this communication, we describe GraphSite, a deep learning-based method utilizing a graph representation of local protein structures and a state-of-the-art graph neural network to classify ligand binding sites. Using neural weighted message passing layers to effectively capture the structural, physicochemical, and evolutionary characteristics of binding pockets mitigates model overfitting and improves the classification accuracy. Indeed, comprehensive cross-validation benchmarks against a large dataset of binding pockets belonging to 14 diverse functional classes demonstrate that GraphSite yields the class-weighted F1-score of 81.7%, outperforming other approaches such as molecular docking and binding site matching. Further, it also generalizes well to unseen data with the F1-score of 70.7%, which is the expected performance in real-world applications. We also discuss new directions to improve and extend GraphSite in the future.

Natural products for infectious microbes and diseases: an overview of sources, compounds, and chemical diversities

As coronavirus disease 2019 (COVID-19) threatens human health globally, infectious disorders have become one of the most challenging problem for the medical community. Natural products (NP) have been a prolific source of antimicrobial agents with widely divergent structures and a range vast biological activities. A dataset comprising 618 articles, including 646 NP-based compounds from 672 species of natural sources with biological activities against 21 infectious pathogens from five categories, was assembled through manual selection of published articles. These data were used to identify 268 NP-based compounds classified into ten groups, which were used for network pharmacology analysis to capture the most promising lead-compounds such as agelasine D, dicumarol, dihydroartemisinin and pyridomycin. The distribution of maximum Tanimoto scores indicated that compounds which inhibited parasites exhibited low diversity, whereas the chemistries inhibiting bacteria, fungi, and viruses showed more structural diversity. A total of 331 species of medicinal plants with compounds exhibiting antimicrobial activities were selected to classify the family sources. The family Asteraceae possesses various compounds against C. neoformans, the family Anacardiaceae has compounds against Salmonella typhi, the family Cucurbitacea against the human immunodeficiency virus (HIV), and the family Ancistrocladaceae against Plasmodium. This review summarizes currently available data on NP-based antimicrobials against refractory infections to provide information for further discovery of drugs and synthetic strategies for anti-infectious agents.

A New Strategy for Multitarget Drug Discovery/Repositioning Through the Identification of Similar 3D Amino Acid Patterns Among Proteins Structures: The Case of Tafluprost and its Effects on Cardiac Ion Channels

The identification of similar three-dimensional (3D) amino acid patterns among different proteins might be helpful to explain the polypharmacological profile of many currently used drugs. Also, it would be a reasonable first step for the design of novel multitarget compounds. Most of the current computational tools employed for this aim are limited to the comparisons among known binding sites, and do not consider several additional important 3D patterns such as allosteric sites or other conserved motifs. In the present work, we introduce Geomfinder2.0, which is a new and improved version of our previously described algorithm for the deep exploration and discovery of similar and druggable 3D patterns. As compared with the original version, substantial improvements that have been incorporated to our software allow: (i) to compare quaternary structures, (ii) to deal with a list of pairs of structures, (iii) to know how druggable is the zone where similar 3D patterns are detected and (iv) to significantly reduce the execution time. Thus, the new algorithm achieves up to 353x speedup as compared to the previous sequential version, allowing the exploration of a significant number of quaternary structures in a reasonable time. In order to illustrate the potential of the updated Geomfinder version, we show a case of use in which similar 3D patterns were detected in the cardiac ions channels NaV1.5 and TASK-1. These channels are quite different in terms of structure, sequence and function and both have been regarded as important targets for drugs aimed at treating atrial fibrillation. Finally, we describe the in vitro effects of tafluprost (a drug currently used to treat glaucoma, which was identified as a novel putative ligand of NaV1.5 and TASK-1) upon both ion channels’ activity and discuss its possible repositioning as a novel antiarrhythmic drug.

ProBiS-Dock: A Hybrid Multitemplate Homology Flexible Docking Algorithm Enabled by Protein Binding Site Comparison

The protein data bank (PDB) is a rich source of protein ligand structures, but ligands are not explicitly used in current docking algorithms. We have developed ProBiS-Dock, a docking algorithm complementary to the ProBiS-Dock Database (J. Chem. Inf. Model. 2021, 61, 4097-4107) that treats small molecules and proteins as fully flexible entities and allows conformational changes in both after ligand binding. A new scoring function is described that consists of a binding site-specific scoring function (ProBiS-Score) and a general statistical scoring function. ProBiS-Dock enables rapid docking of small molecules to proteins and has been successfully validated in silico against standard benchmarks. It enables rapid search for new active ligands by leveraging existing knowledge in the PDB. The potential of the software for drug development has been confirmed in vitro by the discovery of new inhibitors of human indoleamine 2,3-dioxygenase 1, an enzyme that is an attractive target for cancer therapy and catalyzes the first rate-determining step of l-tryptophan metabolism via the kynurenine pathway. The software is freely available to academic users at http://insilab.org/probisdock.

MeCOM: A Method for Comparing Three-Dimensional Metalloenzyme Active Sites

Since metalloenzymes are a large collection of metal ion(s) dependent enzymes, comparison analyses of metalloenzyme active sites are critical for metalloenzyme de novo design, function investigation, and inhibitor development. Here, we report a method named MeCOM for comparing metalloenzyme active sites. It is characterized by metal ion(s) centric active site recognition and three-dimensional superimposition using α-carbon or pharmacophore features. The test results revealed that for the given metalloenzymes, MeCOM could effectively recognize the active sites, extract active site features, and superimpose the active sites; it also could correctly identify similar active sites, differentiate dissimilar active sites, and evaluate the similarity degree. Moreover, MeCOM showed potential to establish new associations between structurally distinct metalloenzymes by active site comparison. MeCOM is freely available at https://mecom.ddtmlab.org.

Fragment Hotspot Mapping to Identify Selectivity-Determining Regions between Related Proteins

Selectivity is a crucial property in small molecule development. Binding site comparisons within a protein family are a key piece of information when aiming to modulate the selectivity profile of a compound. Binding site differences can be exploited to confer selectivity for a specific target, while shared areas can provide insights into polypharmacology. As the quantity of structural data grows, automated methods are needed to process, summarize, and present these data to users. We present a computational method that provides quantitative and data-driven summaries of the available binding site information from an ensemble of structures of the same protein. The resulting ensemble maps identify the key interactions important for ligand binding in the ensemble. The comparison of ensemble maps of related proteins enables the identification of selectivity-determining regions within a protein family. We applied the method to three examples from the well-researched human bromodomain and kinase families, demonstrating that the method is able to identify selectivity-determining regions that have been used to introduce selectivity in past drug discovery campaigns. We then illustrate how the resulting maps can be used to automate comparisons across a target protein family.

Role of Fyn Kinase Inhibitors in Switching Neuroinflammatory Pathways

Fyn kinase is a member of the Src non-receptor tyrosine kinase family. Fyn is involved in multiple signaling pathways extending from cell proliferation and differentiation to cell adhesion and cell motility, and it has been found to be overexpressed in various types of cancers. In the central nervous system, Fyn exerts several different functions such as axon-glial signal transduction, oligodendrocyte maturation and myelination, and it is implicated in neuroinflammatory processes. Based on these premises, Fyn emerges as an attractive target in cancer and neurodegenerative disease therapy, particularly Alzheimer disease (AD), based on its activation by Aβ via cellular prion protein and its interaction with tau protein. However, Fyn is also a challenging target since the Fyn inhibitors discovered so far, due to the relevant homology of Fyn with other kinases, suffer from off-target effects. This review covers the efforts performed in the last decade to identify and optimize small molecules that effectively inhibit Fyn, both in enzymatic and in cell assays, including drug repositioning practices, as an opportunity of therapeutic intervention in neurodegeneration.

CRNNTL: Convolutional Recurrent Neural Network and Transfer Learning for QSAR Modeling in Organic Drug and Material Discovery

Molecular latent representations, derived from autoencoders (AEs), have been widely used for drug or material discovery over the past couple of years. In particular, a variety of machine learning methods based on latent representations have shown excellent performance on quantitative structure–activity relationship (QSAR) modeling. However, the sequence feature of them has not been considered in most cases. In addition, data scarcity is still the main obstacle for deep learning strategies, especially for bioactivity datasets. In this study, we propose the convolutional recurrent neural network and transfer learning (CRNNTL) method inspired by the applications of polyphonic sound detection and electrocardiogram classification. Our model takes advantage of both convolutional and recurrent neural networks for feature extraction, as well as the data augmentation method. According to QSAR modeling on 27 datasets, CRNNTL can outperform or compete with state-of-art methods in both drug and material properties. In addition, the performances on one isomers-based dataset indicate that its excellent performance results from the improved ability in global feature extraction when the ability of the local one is maintained. Then, the transfer learning results show that CRNNTL can overcome data scarcity when choosing relative source datasets. Finally, the high versatility of our model is shown by using different latent representations as inputs from other types of AEs.

Applications of Artificial Intelligence in Drug Design: Opportunities and Challenges

Artificial intelligence (AI) has undergone rapid development in recent years and has been successfully applied to real-world problems such as drug design. In this chapter, we review recent applications of AI to problems in drug design including virtual screening, computer-aided synthesis planning, and de novo molecule generation, with a focus on the limitations of the application of AI therein and opportunities for improvement. Furthermore, we discuss the broader challenges imposed by AI in translating theoretical practice to real-world drug design; including quantifying prediction uncertainty and explaining model behavior.

Enhancing preclinical drug discovery with artificial intelligence

Artificial intelligence (AI) is becoming an integral part of drug discovery. It has the potential to deliver across the drug discovery and development value chain, starting from target identification and reaching through clinical development. In this review, we provide an overview of current AI technologies and a glimpse of how AI is reimagining preclinical drug discovery by highlighting examples where AI has made a real impact. Considering the excitement and hyperbole surrounding AI in drug discovery, we aim to present a realistic view by discussing both opportunities and challenges in adopting AI in drug discovery.

Unexpected similarity between HIV-1 reverse transcriptase and tumor necrosis factor binding sites revealed by computer vision

Rationalizing the identification of hidden similarities across the repertoire of druggable protein cavities remains a major hurdle to a true proteome-wide structure-based discovery of novel drug candidates. We recently described a new computational approach (ProCare), inspired by numerical image processing, to identify local similarities in fragment-based subpockets. During the validation of the method, we unexpectedly identified a possible similarity in the binding pockets of two unrelated targets, human tumor necrosis factor alpha (TNF-α) and HIV-1 reverse transcriptase (HIV-1 RT). Microscale thermophoresis experiments confirmed the ProCare prediction as two of the three tested and FDA-approved HIV-1 RT inhibitors indeed bind to soluble human TNF-α trimer. Interestingly, the herein disclosed similarity could be revealed neither by state-of-the-art binding sites comparison methods nor by ligand-based pairwise similarity searches, suggesting that the point cloud registration approach implemented in ProCare, is uniquely suited to identify local and unobvious similarities among totally unrelated targets.

Atomistic insights into the selective therapeutic activity of 6-(2,4-difluorophenoxy)-5-((ethylmethyl)pyridine-3-yl)-8-methylpyrrolo[1,2-a]pyrazin-1(2H)-one towards bromodomain-containing proteins

Cross-target effect has been one of the major mechanisms of drug toxicity, this has necessitated the design of inhibitors that are specifically tailored to target particular biomolecules. 6-(2,4-difluorophenoxy)-5-((ethylmethyl)pyridine-3-yl)-8-methylpyrrolo[1,2-a] pyrazin-1(2H)-one (Cpd38) is an inhibitor possessing high inhibition rate and tailored specificity towards bromodomain-containing protein 4 (BRD4). In this research, we used an array of computational techniques to provide insight at the atomistic level the specific targeting of BRD4 by Cpd38 relative to the binding of Cpd38 with E1A binding protein P300 (EP300); another bromodomain-containing protein (BCP). Comparatively, binding of Cpd38 improved the conformational stability and compactness of BRD4 protein when compared to the Cpd38 bound EP300. Also, Cpd38 induced a conformational change in the active site of BRD4 that facilitated a complementary pose between Cpd38 and BRD4 suitable for effective atomistic interactions. Expectedly, thermodynamic calculations revealed that the Cpd38-BRD4 system had higher binding energy (-36.11 Kcal/mol) than the Cpd38-EP300 system with a free binding energy of -15.86 Kcal/mol. Noteworthy is the opposing role Trp81 (acting as hydrogen bond acceptor) and Pro1074 (acting as hydrogen bond donor) found on the WPF and LPF loops respectively play in maintaining Cpd38 stability. Furthermore, the hydrogen bond acceptor/donator ratio was approximately 4:1 in Cpd38-BRD4 system compared with 2:1 in Cpd38-EP300 system. Taken together, atomistic insights and structural perspectives detailed in this report supplements the experimental report supporting the improved selectivity of Cpd38 for BRD4 ahead of other BCPs while providing leeway for the future design of BET selective agents with better pharmacological profile.

A fast protein binding site comparison algorithm for proteome-wide protein function prediction and drug repurposing

The expansion of three-dimensional protein structures and enhanced computing power have significantly facilitated our understanding of protein sequence/structure/function relationships. A challenge in structural genomics is to predict the function of uncharacterized proteins. Protein function deconvolution based on global sequence or structural homology is impracticable when a protein relates to no other proteins with known function, and in such cases, functional relationships can be established by detecting their local ligand binding site similarity. Here, we introduce a sequence order-independent comparison algorithm, PocketShape, for structural proteome-wide exploration of protein functional site by fully considering the geometry of the backbones, orientation of the sidechains, and physiochemical properties of the pocket-lining residues. PocketShape is efficient in distinguishing similar from dissimilar ligand binding site pairs by retrieving 99.3% of the similar pairs while rejecting 100% of the dissimilar pairs on a dataset containing 1538 binding site pairs. This method successfully classifies 83 enzyme structures with diverse functions into 12 clusters, which is highly in accordance with the actual structural classification of proteins classification. PocketShape also achieves superior performances than other methods in protein profiling based on experimental data. Potential new applications for representative SARS-CoV-2 drugs Remdesivir and 11a are predicted. The high accuracy and time-efficient characteristics of PocketShape will undoubtedly make it a promising complementary tool for proteome-wide protein function inference and drug repurposing study.

CAVIAR: a method for automatic cavity detection, description and decomposition into subcavities

The accurate description of protein binding sites is essential to the determination of similarity and the application of machine learning methods to relate the binding sites to observed functions. This work describes CAVIAR, a new open source tool for generating descriptors for binding sites, using protein structures in PDB and mmCIF format as well as trajectory frames from molecular dynamics simulations as input. The applicability of CAVIAR descriptors is showcased by computing machine learning predictions of binding site ligandability. The method can also automatically assign subcavities, even in the absence of a bound ligand. The defined subpockets mimic the empirical definitions used in medicinal chemistry projects. It is shown that the experimental binding affinity scales relatively well with the number of subcavities filled by the ligand, with compounds binding to more than three subcavities having nanomolar or better affinities to the target. The CAVIAR descriptors and methods can be used in any machine learning-based investigations of problems involving binding sites, from protein engineering to hit identification. The full software code is available on GitHub and a conda package is hosted on Anaconda cloud.

Exploring protein hotspots by optimized fragment pharmacophores

Fragment-based drug design has introduced a bottom-up process for drug development, with improved sampling of chemical space and increased effectiveness in early drug discovery. Here, we combine the use of pharmacophores, the most general concept of representing drug-target interactions with the theory of protein hotspots, to develop a design protocol for fragment libraries. The SpotXplorer approach compiles small fragment libraries that maximize the coverage of experimentally confirmed binding pharmacophores at the most preferred hotspots. The efficiency of this approach is demonstrated with a pilot library of 96 fragment-sized compounds (SpotXplorer0) that is validated on popular target classes and emerging drug targets. Biochemical screening against a set of GPCRs and proteases retrieves compounds containing an average of 70% of known pharmacophores for these targets. More importantly, SpotXplorer0 screening identifies confirmed hits against recently established challenging targets such as the histone methyltransferase SETD2, the main protease (3CLPro) and the NSP3 macrodomain of SARS-CoV-2.

Discovery of novel diarylamides as orally active diuretics targeting urea transporters

Urea transporters (UT) play a vital role in the mechanism of urine concentration and are recognized as novel targets for the development of salt-sparing diuretics. Thus, UT inhibitors are promising for development as novel diuretics. In the present study, a novel UT inhibitor with a diarylamide scaffold was discovered by high-throughput screening. Optimization of the inhibitor led to the identification of a promising preclinical candidate, N-[4-(acetylamino)phenyl]-5-nitrofuran-2-carboxamide (1H), with excellent in vitro UT inhibitory activity at the submicromolar level. The half maximal inhibitory concentrations of 1H against UT-B in mouse, rat, and human erythrocyte were 1.60, 0.64, and 0.13 μmol/L, respectively. Further investigation suggested that 8 μmol/L 1H more powerfully inhibited UT-A1 at a rate of 86.8% than UT-B at a rate of 73.9% in MDCK cell models. Most interestingly, we found for the first time that oral administration of 1H at a dose of 100 mg/kg showed superior diuretic effect in vivo without causing electrolyte imbalance in rats. Additionally, 1H did not exhibit apparent toxicity in vivo and in vitro, and possessed favorable pharmacokinetic characteristics. 1H shows promise as a novel diuretic to treat hyponatremia accompanied with volume expansion and may cause few side effects.

Chemogenomics and bioinformatics approaches for prioritizing kinases as drug targets for neglected tropical diseases

Neglected tropical diseases (NTDs) are a group of twenty-one diseases classified by the World Health Organization that prevail in regions with tropical and subtropical climate and affect more than one billion people. There is an urgent need to develop new and safer drugs for these diseases. Protein kinases are a potential class of targets for developing new drugs against NTDs, since they play crucial role in many biological processes, such as signaling pathways, regulating cellular communication, division, metabolism and death. Bioinformatics is a field that aims to organize large amounts of biological data as well as develop and use tools for understanding and analyze them in order to produce meaningful information in a biological manner. In combination with chemogenomics, which analyzes chemical-biological interactions to screen ligands against selected targets families, these approaches can be used to stablish a rational strategy for prioritizing new drug targets for NTDs. Here, we describe how bioinformatics and chemogenomics tools can help to identify protein kinases and their potential inhibitors for the development of new drugs for NTDs. We present a review of bioinformatics tools and techniques that can be used to define an organisms kinome for drug prioritization, drug and target repurposing, multi-quinase inhibition approachs and selectivity profiling. We also present some successful examples of the application of such approaches in recent case studies.

Applications of Genome-Wide Screening and Systems Biology Approaches in Drug Repositioning

Modern drug discovery through de novo drug discovery entails high financial costs, low success rates, and lengthy trial periods. Drug repositioning presents a suitable approach for overcoming these issues by re-evaluating biological targets and modes of action of approved drugs. Coupling high-throughput technologies with genome-wide essentiality screens, network analysis, genome-scale metabolic modeling, and machine learning techniques enables the proposal of new drug-target signatures and uncovers unanticipated modes of action for available drugs. Here, we discuss the current issues associated with drug repositioning in light of curated high-throughput multi-omic databases, genome-wide screening technologies, and their application in systems biology/medicine approaches.

In Silico Drug Repositioning for Chagas Disease

Chagas disease is an infectious tropical disease included within the group of neglected tropical diseases. Though historically endemic to Latin America, it has lately spread to high-income countries due to human migration. At present, there are only two available drugs, nifurtimox and benznidazole, approved for this treatment, both with considerable side-effects (which often result in treatment interruption) and limited efficacy in the chronic stage of the disease in adults. Drug repositioning involves finding novel therapeutic indications for known drugs, including approved, withdrawn, abandoned and investigational drugs. It is today a broadly applied approach to develop innovative medications, since indication shifts are built on existing safety, ADME and manufacturing information, thus greatly shortening development timeframes. Drug repositioning has been signaled as a particularly interesting strategy to search for new therapeutic solutions for neglected and rare conditions, which traditionally present limited commercial interest and are mostly covered by the public sector and not-for-profit initiatives and organizations. Here, we review the applications of computer-aided technologies as systematic approaches to drug repositioning in the field of Chagas disease. In silico screening represents the most explored approach, whereas other rational methods such as network-based and signature-based approximations have still not been applied.

Virtual screening web servers: designing chemical probes and drug candidates in the cyberspace

The interplay between life sciences and advancing technology drives a continuous cycle of chemical data growth; these data are most often stored in open or partially open databases. In parallel, many different types of algorithms are being developed to manipulate these chemical objects and associated bioactivity data. Virtual screening methods are among the most popular computational approaches in pharmaceutical research. Today, user-friendly web-based tools are available to help scientists perform virtual screening experiments. This article provides an overview of internet resources enabling and supporting chemical biology and early drug discovery with a main emphasis on web servers dedicated to virtual ligand screening and small-molecule docking. This survey first introduces some key concepts and then presents recent and easily accessible virtual screening and related target-fishing tools as well as briefly discusses case studies enabled by some of these web services. Notwithstanding further improvements, already available web-based tools not only contribute to the design of bioactive molecules and assist drug repositioning but also help to generate new ideas and explore different hypotheses in a timely fashion while contributing to teaching in the field of drug development.

GSP4PDB: a web tool to visualize, search and explore protein-ligand structural patterns

BACKGROUND

In the field of protein engineering and biotechnology, the discovery and characterization of structural patterns is highly relevant as these patterns can give fundamental insights into protein-ligand interaction and protein function. This paper presents GSP4PDB, a bioinformatics web tool that enables the user to visualize, search and explore protein-ligand structural patterns within the entire Protein Data Bank.

RESULTS

We introduce the notion of graph-based structural pattern (GSP) as an abstract model for representing protein-ligand interactions. A GSP is a graph where the nodes represent entities of the protein-ligand complex (amino acids and ligands) and the edges represent structural relationships (e.g. distances ligand - amino acid). The novel feature of GSP4PDB is a simple and intuitive graphical interface where the user can "draw" a GSP and execute its search in a relational database containing the structural data of each PDB entry. The results of the search are displayed using the same graph-based representation of the pattern. The user can further explore and analyse the results using a wide range of filters, or download their related information for external post-processing and analysis.

CONCLUSIONS

GSP4PDB is a user-friendly and efficient application to search and discover new patterns of protein-ligand interaction.

The Thiol-polyamine Metabolism of Trypanosoma cruzi: Molecular Targets and Drug Repurposing Strategies

Chagas´ disease continues to be a challenging and neglected public health problem in many American countries. The etiologic agent, Trypanosoma cruzi, develops intracellularly in the mammalian host, which hinders treatment efficacy. Progress in the knowledge of parasite biology and host-pathogen interaction has not been paralleled by the development of novel, safe and effective therapeutic options. It is then urgent to seek for novel therapeutic candidates and to implement drug discovery strategies that may accelerate the discovery process. The most appealing targets for pharmacological intervention are those essential for the pathogen and, whenever possible, absent or significantly different from the host homolog. The thiol-polyamine metabolism of T. cruzi offers interesting candidates for a rational design of selective drugs. In this respect, here we critically review the state of the art of the thiolpolyamine metabolism of T. cruzi and the pharmacological potential of its components. On the other hand, drug repurposing emerged as a valid strategy to identify new biological activities for drugs in clinical use, while significantly shortening the long time and high cost associated with de novo drug discovery approaches. Thus, we also discuss the different drug repurposing strategies available with a special emphasis in their applications to the identification of drug candidates targeting essential components of the thiol-polyamine metabolism of T. cruzi.

Binding site matching in rational drug design: algorithms and applications

Interactions between proteins and small molecules are critical for biological functions. These interactions often occur in small cavities within protein structures, known as ligand-binding pockets. Understanding the physicochemical qualities of binding pockets is essential to improve not only our basic knowledge of biological systems, but also drug development procedures. In order to quantify similarities among pockets in terms of their geometries and chemical properties, either bound ligands can be compared to one another or binding sites can be matched directly. Both perspectives routinely take advantage of computational methods including various techniques to represent and compare small molecules as well as local protein structures. In this review, we survey 12 tools widely used to match pockets. These methods are divided into five categories based on the algorithm implemented to construct binding-site alignments. In addition to the comprehensive analysis of their algorithms, test sets and the performance of each method are described. We also discuss general pharmacological applications of computational pocket matching in drug repurposing, polypharmacology and side effects. Reflecting on the importance of these techniques in drug discovery, in the end, we elaborate on the development of more accurate meta-predictors, the incorporation of protein flexibility and the integration of powerful artificial intelligence technologies such as deep learning.

Exhaustive Repertoire of Druggable Cavities at Protein-Protein Interfaces of Known Three-Dimensional Structure

Protein-protein interactions (PPIs) offer the unique opportunity to tailor ligands aimed at specifically stabilizing or disrupting the corresponding interfaces and providing a safer alternative to conventional ligands targeting monomeric macromolecules. Selecting biologically relevant protein-protein interfaces for either stabilization or disruption by small molecules is usually biology-driven on a case-by-case basis and does not follow a structural rationale that could be applied to an entire interactome. We herewith provide a first step to the latter goal by using a fully automated and structure-based workflow, applicable to any PPI of known three-dimensional (3D) structure, to identify and prioritize druggable cavities at and nearby PPIs of pharmacological interest. When applied to the entire Protein Data Bank, 164 514 druggable cavities were identified and classified in four groups (interfacial, rim, allosteric, orthosteric) according to their properties and spatial locations. Systematic comparison of PPI cavities with pockets deduced from druggable protein-ligand complexes shows almost no overlap in property space, suggesting that even the most druggable PPI cavities are unlikely to be addressed with conventional drug-like compound libraries. The archive is freely accessible at http://drugdesign.unistra.fr/ppiome .

An Integrated System Biology Approach Yields Drug Repositioning Candidates for the Treatment of Heart Failure

Identifying effective pharmacological treatments for heart failure (HF) patients remains critically important. Given that the development of drugs de novo is expensive and time consuming, drug repositioning has become an increasingly important branch. In the present study, we propose a two-step drug repositioning pipeline and investigate the novel therapeutic effects of existing drugs approved by the US Food and Drug Administration to discover potential therapeutic drugs for HF. In the first step, we compared the gene expression pattern of HF patients with drug-induced gene expression profiles to obtain preliminary candidates. In the second step, we performed a systems biology approach based on the known protein–protein interaction information and targets of drugs to narrow down preliminary candidates to obtain final candidates. Drug set enrichment analysis and literature search were applied to assess the performance of our repositioning approach. We also constructed a mode of action network for each candidate and performed pathway analysis for each candidate using genes contained in their mode of action network to uncover pathways that potentially reflect the mechanisms of candidates’ therapeutic efficacy to HF. We discovered numerous preliminary candidates, some of which are used in clinical practice and supported by the literature. The final candidates contained nearly all of the preliminary candidates supported by previous studies. Drug set enrichment analysis and literature search support the validity of our repositioning approach. The mode of action network for each candidate not only displayed the underlying mechanisms of drug efficacy but also uncovered potential biomarkers and therapeutic targets for HF. Our two-step drug repositioning approach is efficient to find candidates with potential therapeutic efficiency to HF supported by the literature and might be of particular use in the discovery of novel effective pharmacological therapies for HF.

Insights into Computational Drug Repurposing for Neurodegenerative Disease

Computational drug repurposing has the ability to remarkably reduce drug development time and cost in an era where these factors are prohibitively high. Several examples of successful repurposed drugs exist in fields such as oncology, diabetes, leprosy, inflammatory bowel disease, among others, however computational drug repurposing in neurodegenerative disease has presented several unique challenges stemming from the lack of validation methods and difficulty in studying heterogenous diseases of aging. Here, we examine existing approaches to computational drug repurposing, including molecular, clinical, and biophysical methods, and propose data sources and methods to advance computational drug repurposing in neurodegenerative disease using Alzheimer's disease as an example.

3D-PP: A Tool for Discovering Conserved Three-Dimensional Protein Patterns

Discovering conserved three-dimensional (3D) patterns among protein structures may provide valuable insights into protein classification, functional annotations or the rational design of multi-target drugs. Thus, several computational tools have been developed to discover and compare protein 3D-patterns. However, most of them only consider previously known 3D-patterns such as orthosteric binding sites or structural motifs. This fact makes necessary the development of new methods for the identification of all possible 3D-patterns that exist in protein structures (allosteric sites, enzyme-cofactor interaction motifs, among others). In this work, we present 3D-PP, a new free access web server for the discovery and recognition all similar 3D amino acid patterns among a set of proteins structures (independent of their sequence similarity). This new tool does not require any previous structural knowledge about ligands, and all data are organized in a high-performance graph database. The input can be a text file with the PDB access codes or a zip file of PDB coordinates regardless of the origin of the structural data: X-ray crystallographic experiments or in silico homology modeling. The results are presented as lists of sequence patterns that can be further analyzed within the web page. We tested the accuracy and suitability of 3D-PP using two sets of proteins coming from the Protein Data Bank: (a) Zinc finger containing and (b) Serotonin target proteins. We also evaluated its usefulness for the discovering of new 3D-patterns, using a set of protein structures coming from in silico homology modeling methodologies, all of which are overexpressed in different types of cancer. Results indicate that 3D-PP is a reliable, flexible and friendly-user tool to identify conserved structural motifs, which could be relevant to improve the knowledge about protein function or classification. The web server can be freely utilized at https://appsbio.utalca.cl/3d-pp/.

Binding site characterization - similarity, promiscuity, and druggability

The elucidation of non-obvious binding site similarities has provided useful indications for the establishment of polypharmacology, the identification of potential off-targets, or the repurposing of known drugs. The concept underlying all of these approaches is promiscuous binding which can be analyzed from a ligand-based or a binding site-based perspective. Herein, we applied methods for the automated analysis and comparison of protein binding sites to study promiscuous binding on a novel dataset of sites in complex with ligands sharing common shape and physicochemical properties. We show the suitability of this dataset for the benchmarking of novel binding site comparison methods. Our investigations also reveal promising directions for further in-depth analyses of promiscuity and druggability in a pocket-centered manner. Drawbacks concerning binding site similarity assessment and druggability prediction are outlined, enabling researchers to avoid the typical pitfalls of binding site analyses.

PocketPipe: A computational pipeline for integratedPocketome prediction and comparison

Functional characterisation of proteins often depends on specific interactions with other molecules. In the drug discovery scenario, the ability of a protein to bind with drug-like molecule with a high affinity is referred as druggability. Deciphering such druggable binding pockets on proteins plays an important role in structure-based drug designing studies. Moreover, availability of plethora of structural data poses a need automated pipelines which can efficiently integrate robust algorithms towards large-scale pocket identification and comparison. These pipelines have direct applicability on off-target analysis, drug repurposing and structural prioritization of drug targets in pathogenic microbes. However, currently there is a paucity of such efficient pipelines. Hence, by this study a highly optimized shell script based pipeline (PocketPipe) has been developed with seamless integration of robust algorithms namely, P2Rank (predicts binding sites based on machine learning) and PocketMatch-v2.1 (compares binding pockets by residue-based method), for pocketome generation and comparison, respectively. The process of input workflow and various steps carried out by PocketPipe and the output results are well documented in the operating manual. On execution, the pipeline features seamless operability, high scalability, dynamic file handling and results parsing. PocketPipe is distributed freely under GNU GPL license and can be downloaded at https://github.com/inpacdb/PocketPipe

TeachOpenCADD: a teaching platform for computer-aided drug design using open source packages and data

Owing to the increase in freely available software and data for cheminformatics and structural bioinformatics, research for computer-aided drug design (CADD) is more and more built on modular, reproducible, and easy-to-share pipelines. While documentation for such tools is available, there are only a few freely accessible examples that teach the underlying concepts focused on CADD, especially addressing users new to the field. Here, we present TeachOpenCADD, a teaching platform developed by students for students, using open source compound and protein data as well as basic and CADD-related Python packages. We provide interactive Jupyter notebooks for central CADD topics, integrating theoretical background and practical code. TeachOpenCADD is freely available on GitHub: https://github.com/volkamerlab/TeachOpenCADD .

Binding site comparisons for target-centered drug discovery

INTRODUCTION

The success of binding site comparisons in drug discovery is based on the recognized fact that many different proteins have similar binding sites. Indeed, binding site comparisons have found many uses in drug development and have the potential to dramatically cut the cost and shorten the time necessary for the development of new drugs. Areas covered: The authors review recent methods for comparing protein binding sites and their use in drug repurposing and polypharmacology. They examine emerging fields including the use of binding site comparisons in precision medicine, the prediction of structured water molecules, the search for targets of natural compounds, and their application in the development of protein-based drugs by loop modeling and for comparison of RNA binding sites. Expert opinion: Binding site comparisons have produced many interesting results in drug development, but relatively little work has been done on protein-protein interaction sites, which are particularly relevant in view of the success of biological drugs. Growth of protein loop modeling for modulating biological drugs is anticipated. The fusion of currently distinct methods for the comparison of RNA and protein binding sites into a single comprehensive approach could allow the search for new selective ribosomal antibiotics and initiate pharmaceutical research into other nucleoproteins.