Encoding protein-ligand interaction patterns in fingerprints and graphs

Structural analysis of neomycin B and kanamycin A binding Aminoglycosides Modifying Enzymes (AME) and bacterial ribosomal RNA

Aminoglycosides are crucial antibiotics facing challenges from bacterial resistance. This study addresses the importance of aminoglycoside modifying enzymes in the context of escalating resistance. Drawing upon over two decades of structural data in the Protein Data Bank, we focused on two key antibiotics, neomycin B and kanamycin A, to explore how the aminoglycoside structure is exploited by this family of enzymes. A systematic comparison across diverse enzymes and the RNA A-site target identified common characteristics in the recognition mode, while assessing the adaptability of neomycin B and kanamycin A in various environments.

SiteMine: Large-scale binding site similarity searching in protein structure databases

Drug discovery and design challenges, such as drug repurposing, analyzing protein-ligand and protein-protein complexes, ligand promiscuity studies, or function prediction, can be addressed by protein binding site similarity analysis. Although numerous tools exist, they all have individual strengths and drawbacks with regard to run time, provision of structure superpositions, and applicability to diverse application domains. Here, we introduce SiteMine, an all-in-one database-driven, alignment-providing binding site similarity search tool to tackle the most pressing challenges of binding site comparison. The performance of SiteMine is evaluated on the ProSPECCTs benchmark, showing a promising performance on most of the data sets. The method performs convincingly regarding all quality criteria for reliable binding site comparison, offering a novel state-of-the-art approach for structure-based molecular design based on binding site comparisons. In a SiteMine showcase, we discuss the high structural similarity between cathepsin L and calpain 1 binding sites and give an outlook on the impact of this finding on structure-based drug design. SiteMine is available at https://uhh.de/naomi.

VirtuousPocketome: a computational tool for screening protein-ligand complexes to identify similar binding sites

Protein residues within binding pockets play a critical role in determining the range of ligands that can interact with a protein, influencing its structure and function. Identifying structural similarities in proteins offers valuable insights into their function and activation mechanisms, aiding in predicting protein-ligand interactions, anticipating off-target effects, and facilitating the development of therapeutic agents. Numerous computational methods assessing global or local similarity in protein cavities have emerged, but their utilization is impeded by complexity, impractical automation for amino acid pattern searches, and an inability to evaluate the dynamics of scrutinized protein-ligand systems. Here, we present a general, automatic and unbiased computational pipeline, named VirtuousPocketome, aimed at screening huge databases of proteins for similar binding pockets starting from an interested protein-ligand complex. We demonstrate the pipeline's potential by exploring a recently-solved human bitter taste receptor, i.e. the TAS2R46, complexed with strychnine. We pinpointed 145 proteins sharing similar binding sites compared to the analysed bitter taste receptor and the enrichment analysis highlighted the related biological processes, molecular functions and cellular components. This work represents the foundation for future studies aimed at understanding the effective role of tastants outside the gustatory system: this could pave the way towards the rationalization of the diet as a supplement to standard pharmacological treatments and the design of novel tastants-inspired compounds to target other proteins involved in specific diseases or disorders. The proposed pipeline is publicly accessible, can be applied to any protein-ligand complex, and could be expanded to screen any database of protein structures.

Systematic analysis, aggregation and visualisation of interaction fingerprints for molecular dynamics simulation data

Computational methods such as molecular docking or molecular dynamics (MD) simulations have been developed to simulate and explore the interactions between biomolecules. However, the interactions obtained using these methods are difficult to analyse and evaluate. Interaction fingerprints (IFPs) have been proposed to derive interactions from static 3D coordinates and transform them into 1D bit vectors. More recently, the concept has been applied to derive IFPs from MD simulations, which adds a layer of complexity by adding the temporal motion and dynamics of a system. As a result, many IFPs are obtained from one MD simulation, resulting in a large number of individual IFPs that are difficult to analyse compared to IFPs derived from static 3D structures. Scientific contribution: We introduce a new method to systematically aggregate IFPs derived from MD simulation data. In addition, we propose visualisations to effectively analyse and compare IFPs derived from MD simulation data to account for the temporal evolution of interactions and to compare IFPs across different MD simulations. This has been implemented as a freely available Python library and can therefore be easily adopted by other researchers and to different MD simulation datasets.

Fingerprinting Interactions between Proteins and Ligands for Facilitating Machine Learning in Drug Discovery

Molecular recognition is fundamental in biology, underpinning intricate processes through specific protein-ligand interactions. This understanding is pivotal in drug discovery, yet traditional experimental methods face limitations in exploring the vast chemical space. Computational approaches, notably quantitative structure-activity/property relationship analysis, have gained prominence. Molecular fingerprints encode molecular structures and serve as property profiles, which are essential in drug discovery. While two-dimensional (2D) fingerprints are commonly used, three-dimensional (3D) structural interaction fingerprints offer enhanced structural features specific to target proteins. Machine learning models trained on interaction fingerprints enable precise binding prediction. Recent focus has shifted to structure-based predictive modeling, with machine-learning scoring functions excelling due to feature engineering guided by key interactions. Notably, 3D interaction fingerprints are gaining ground due to their robustness. Various structural interaction fingerprints have been developed and used in drug discovery, each with unique capabilities. This review recapitulates the developed structural interaction fingerprints and provides two case studies to illustrate the power of interaction fingerprint-driven machine learning. The first elucidates structure-activity relationships in β2 adrenoceptor ligands, demonstrating the ability to differentiate agonists and antagonists. The second employs a retrosynthesis-based pre-trained molecular representation to predict protein-ligand dissociation rates, offering insights into binding kinetics. Despite remarkable progress, challenges persist in interpreting complex machine learning models built on 3D fingerprints, emphasizing the need for strategies to make predictions interpretable. Binding site plasticity and induced fit effects pose additional complexities. Interaction fingerprints are promising but require continued research to harness their full potential.

Structure-based, deep-learning models for protein-ligand binding affinity prediction

The launch of AlphaFold series has brought deep-learning techniques into the molecular structural science. As another crucial problem, structure-based prediction of protein-ligand binding affinity urgently calls for advanced computational techniques. Is deep learning ready to decode this problem? Here we review mainstream structure-based, deep-learning approaches for this problem, focusing on molecular representations, learning architectures and model interpretability. A model taxonomy has been generated. To compensate for the lack of valid comparisons among those models, we realized and evaluated representatives from a uniform basis, with the advantages and shortcomings discussed. This review will potentially benefit structure-based drug discovery and related areas.

Predicting the duration of action of β2-adrenergic receptor agonists: Ligand and structure-based approaches

Agonists of the β2 adrenergic receptor (ADRB2) are an important class of medications used for the treatment of respiratory diseases. They can be classified as short acting (SABA) or long acting (LABA), with each class playing a different role in patient management. In this work we explored both ligand-based and structure-based high-throughput approaches to classify β2-agonists based on their duration of action. A completely in-silico prediction pipeline using an AlphaFold generated structure was used for structure-based modelling. Our analysis identified the ligands' 3D structure and lipophilicity as the most relevant features for the prediction of the duration of action. Interaction-based methods were also able to select ligands with the desired duration of action, incorporating the bias directly in the structure-based drug discovery pipeline without the need for further processing.

Rigid Scaffolds Are Promising for Designing Macrocyclic Kinase Inhibitors

Macrocyclic kinase inhibitors (MKIs) are gaining attention due to their favorable selectivity and potential to overcome drug resistance, yet they remain challenging to design because of their novel structures. To facilitate the design and discovery of MKIs, we investigate MKI rational design starting from initial acyclic compounds by performing microsecond-scale atomistic simulations for multiple MKIs, constructing an MKI database, and analyzing MKIs using hierarchical cluster analysis. Our studies demonstrate that the binding modes of MKIs are like those of their corresponding acyclic counterparts against the same kinase targets. Importantly, within the respective binding sites, the MKI scaffolds retain the same conformations as their corresponding acyclic counterparts, demonstrating the rigidity of scaffolds before and after molecular cyclization. The MKI database includes 641 nanomole-level MKIs from 56 human kinases elucidating the features of rigid scaffolds and the core structures of MKIs. Collectively these results and resources can facilitate MKI development.

AQDnet: Deep Neural Network for Protein-Ligand Docking Simulation

We have developed an innovative system, AI QM Docking Net (AQDnet), which utilizes the three-dimensional structure of protein-ligand complexes to predict binding affinity. This system is novel in two respects: first, it significantly expands the training dataset by generating thousands of diverse ligand configurations for each protein-ligand complex and subsequently determining the binding energy of each configuration through quantum computation. Second, we have devised a method that incorporates the atom-centered symmetry function (ACSF), highly effective in describing molecular energies, for the prediction of protein-ligand interactions. These advancements have enabled us to effectively train a neural network to learn the protein-ligand quantum energy landscape (P-L QEL). Consequently, we have achieved a 92.6% top 1 success rate in the CASF-2016 docking power, placing first among all models assessed in the CASF-2016, thus demonstrating the exceptional docking performance of our model.

Roles of Virtual Screening and Molecular Dynamics Simulations in Discovering and Understanding Antimalarial Drugs

Malaria continues to be a global health threat, with approximately 247 million cases worldwide. Despite therapeutic interventions being available, patient compliance is a problem due to the length of treatment. Moreover, drug-resistant strains have emerged over the years, necessitating urgent identification of novel and more potent treatments. Given that traditional drug discovery often requires a great deal of time and resources, most drug discovery efforts now use computational methods. In silico techniques such as quantitative structure-activity relationship (QSAR), docking, and molecular dynamics (MD) can be used to study protein-ligand interactions and determine the potency and safety profile of a set of candidate compounds to help prioritize those tested using assays and animal models. This paper provides an overview of antimalarial drug discovery and the application of computational methods in identifying candidate inhibitors and elucidating their potential mechanisms of action. We conclude with the continued challenges and future perspectives in the field of antimalarial drug discovery.

Mining the Protein Data Bank to inspire fragment library design

The fragment approach has emerged as a method of choice for drug design, as it allows difficult therapeutic targets to be addressed. Success lies in the choice of the screened chemical library and the biophysical screening method, and also in the quality of the selected fragment and structural information used to develop a drug-like ligand. It has recently been proposed that promiscuous compounds, i.e., those that bind to several proteins, present an advantage for the fragment approach because they are likely to give frequent hits in screening. In this study, we searched the Protein Data Bank for fragments with multiple binding modes and targeting different sites. We identified 203 fragments represented by 90 scaffolds, some of which are not or hardly present in commercial fragment libraries. By contrast to other available fragment libraries, the studied set is enriched in fragments with a marked three-dimensional character (download at 10.5281/zenodo.7554649).

AI-accelerated protein-ligand docking for SARS-CoV-2 is 100-fold faster with no significant change in detection

Protein-ligand docking is a computational method for identifying drug leads. The method is capable of narrowing a vast library of compounds down to a tractable size for downstream simulation or experimental testing and is widely used in drug discovery. While there has been progress in accelerating scoring of compounds with artificial intelligence, few works have bridged these successes back to the virtual screening community in terms of utility and forward-looking development. We demonstrate the power of high-speed ML models by scoring 1 billion molecules in under a day (50 k predictions per GPU seconds). We showcase a workflow for docking utilizing surrogate AI-based models as a pre-filter to a standard docking workflow. Our workflow is ten times faster at screening a library of compounds than the standard technique, with an error rate less than 0.01% of detecting the underlying best scoring 0.1% of compounds. Our analysis of the speedup explains that another order of magnitude speedup must come from model accuracy rather than computing speed. In order to drive another order of magnitude of acceleration, we share a benchmark dataset consisting of 200 million 3D complex structures and 2D structure scores across a consistent set of 13 million "in-stock" molecules over 15 receptors, or binding sites, across the SARS-CoV-2 proteome. We believe this is strong evidence for the community to begin focusing on improving the accuracy of surrogate models to improve the ability to screen massive compound libraries 100 × or even 1000 × faster than current techniques and reduce missing top hits. The technique outlined aims to be a fast drop-in replacement for docking for screening billion-scale molecular libraries.

A multilayer dynamic perturbation analysis method for predicting ligand-protein interactions

BACKGROUND

Ligand-protein interactions play a key role in defining protein function, and detecting natural ligands for a given protein is thus a very important bioengineering task. In particular, with the rapid development of AI-based structure prediction algorithms, batch structural models with high reliability and accuracy can be obtained at low cost, giving rise to the urgent requirement for the prediction of natural ligands based on protein structures. In recent years, although several structure-based methods have been developed to predict ligand-binding pockets and ligand-binding sites, accurate and rapid methods are still lacking, especially for the prediction of ligand-binding regions and the spatial extension of ligands in the pockets.

RESULTS

In this paper, we proposed a multilayer dynamics perturbation analysis (MDPA) method for predicting ligand-binding regions based solely on protein structure, which is an extended version of our previously developed fast dynamic perturbation analysis (FDPA) method. In MDPA/FDPA, ligand binding tends to occur in regions that cause large changes in protein conformational dynamics. MDPA, examined using a standard validation dataset of ligand-protein complexes, yielded an averaged ligand-binding site prediction Matthews coefficient of 0.40, with a prediction precision of at least 50% for 71% of the cases. In particular, for 80% of the cases, the predicted ligand-binding region overlaps the natural ligand by at least 50%. The method was also compared with other state-of-the-art structure-based methods.

CONCLUSIONS

MDPA is a structure-based method to detect ligand-binding regions on protein surface. Our calculations suggested that a range of spaces inside the protein pockets has subtle interactions with the protein, which can significantly impact on the overall dynamics of the protein. This work provides a valuable tool as a starting point upon which further docking and analysis methods can be used for natural ligand detection in protein functional annotation. The source code of MDPA method is freely available at: https://github.com/mingdengming/mdpa .

One class classification for the detection of β2 adrenergic receptor agonists using single-ligand dynamic interaction data

G protein-coupled receptors are involved in many biological processes, relaying the extracellular signal inside the cell. Signaling is regulated by the interactions between receptors and their ligands, it can be stimulated by agonists, or inhibited by antagonists or inverse agonists. The development of a new drug targeting a member of this family requires to take into account the pharmacological profile of the designed ligands in order to elicit the desired response. The structure-based virtual screening of chemical libraries may prioritize a specific class of ligands by combining docking results and ligand binding information provided by crystallographic structures. The performance of the method depends on the relevance of the structural data, in particular the conformation of the targeted site, the binding mode of the reference ligand, and the approach used to compare the interactions formed by the docked ligand with those formed by the reference ligand in the crystallographic structure. Here, we propose a new method based on the conformational dynamics of a single protein–ligand reference complex to improve the biased selection of ligands with specific pharmacological properties in a structure-based virtual screening exercise. Interactions patterns between a reference agonist and the receptor, here exemplified on the β2 adrenergic receptor, were extracted from molecular dynamics simulations of the agonist/receptor complex and encoded in graphs used to train a one-class machine learning classifier. Different conditions were tested: low to high affinity agonists, varying simulation duration, considering or ignoring hydrophobic contacts, and tuning of the classifier parametrization. The best models applied to post-process raw data from retrospective virtual screening obtained by docking of test libraries effectively filtered out irrelevant poses, discarding inactive and non-agonist ligands while identifying agonists. Taken together, our results suggest that consistency of the binding mode during the simulation is a key to the success of the method.

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.

Harnessing systematic protein-ligand interaction fingerprints for drug discovery

Determining protein-ligand interaction characteristics and mechanisms is crucial to the drug discovery process. Here, we review recent progress and successful applications of a systematic protein-ligand interaction fingerprint (IFP) approach for investigating proteome-wide protein-ligand interactions for drug development. Specifically, we review the use of this IFP approach for revealing polypharmacology across the kinome, predicting promising targets from which to design allosteric inhibitors and covalent kinase inhibitors, uncovering the binding mechanisms of drugs of interest, and demonstrating resistant mechanisms of specific drugs. Together, we demonstrate that the IFP strategy is efficient and practical for drug design research for protein kinases as targets and is extensible to other protein families.

Fast Local Alignment of Protein Pockets (FLAPP): A System-Compiled Program for Large-Scale Binding Site Alignment

Protein function is a direct consequence of its sequence, structure, and the arrangement at the binding site. Bioinformatics using sequence analysis is typically used to gain a first insight into protein function. Protein structures, on the other hand, provide a higher resolution platform into understanding functions. As the protein structural information is increasingly becoming available through experimental structure determination and through advances in computational methods for structure prediction, the opportunity to utilize these data is also increasing. Structural analysis of small molecule ligand binding sites in particular provides a direct and more accurate window to infer protein function. However, it remains a poorly utilized resource due to the huge computational cost of existing methods that make large-scale structural comparisons of binding sites prohibitive. Here, we present an algorithm called FLAPP that produces very rapid atomic level alignments. By combining clique matching in graphs and the power of modern CPU architectures, FLAPP aligns a typical pair of binding sites at ∼12.5 ms using a single CPU core, ∼1 ms using 12 cores on a standard desktop machine, and performs a PDB-wide scan in 1-2 min. We perform rigorous validation of the algorithm at multiple levels of complexity and show that FLAPP provides accurate alignments. We also present a case study involving vitamin B12 binding sites to showcase the usefulness of FLAPP for performing an exhaustive alignment-based PDB-wide scan. We expect that this tool will be invaluable to the scientific community to quickly align millions of site pairs on a normal desktop machine to gain insights into protein function and drug discovery for drug target and off-target identification and polypharmacology.

fingeRNAt—A novel tool for high-throughput analysis of nucleic acid-ligand interactions

Computational methods play a pivotal role in drug discovery and are widely applied in virtual screening, structure optimization, and compound activity profiling. Over the last decades, almost all the attention in medicinal chemistry has been directed to protein-ligand binding, and computational tools have been created with this target in mind. With novel discoveries of functional RNAs and their possible applications, RNAs have gained considerable attention as potential drug targets. However, the availability of bioinformatics tools for nucleic acids is limited. Here, we introduce fingeRNAt—a software tool for detecting non-covalent interactions formed in complexes of nucleic acids with ligands. The program detects nine types of interactions: (i) hydrogen and (ii) halogen bonds, (iii) cation-anion, (iv) pi-cation, (v) pi-anion, (vi) pi-stacking, (vii) inorganic ion-mediated, (viii) water-mediated, and (ix) lipophilic interactions. However, the scope of detected interactions can be easily expanded using a simple plugin system. In addition, detected interactions can be visualized using the associated PyMOL plugin, which facilitates the analysis of medium-throughput molecular complexes. Interactions are also encoded and stored as a bioinformatics-friendly Structural Interaction Fingerprint (SIFt)—a binary string where the respective bit in the fingerprint is set to 1 if a particular interaction is present and to 0 otherwise. This output format, in turn, enables high-throughput analysis of interaction data using data analysis techniques. We present applications of fingeRNAt-generated interaction fingerprints for visual and computational analysis of RNA-ligand complexes, including analysis of interactions formed in experimentally determined RNA-small molecule ligand complexes deposited in the Protein Data Bank. We propose interaction fingerprint-based similarity as an alternative measure to RMSD to recapitulate complexes with similar interactions but different folding. We present an application of interaction fingerprints for the clustering of molecular complexes. This approach can be used to group ligands that form similar binding networks and thus have similar biological properties. The fingeRNAt software is freely available at https://github.com/n-szulc/fingeRNAt.

Delta Machine Learning to Improve Scoring-Ranking-Screening Performances of Protein-Ligand Scoring Functions

Protein-ligand scoring functions are widely used in structure-based drug design for fast evaluation of protein-ligand interactions, and it is of strong interest to develop scoring functions with machine-learning approaches. In this work, by expanding the training set, developing physically meaningful features, employing our recently developed linear empirical scoring function Lin_F9 (Yang, C. J. Chem. Inf. Model. 2021, 61, 4630-4644) as the baseline, and applying extreme gradient boosting (XGBoost) with Δ-machine learning, we have further improved the robustness and applicability of machine-learning scoring functions. Besides the top performances for scoring-ranking-screening power tests of the CASF-2016 benchmark, the new scoring function Δ_{Lin_F9}XGB also achieves superior scoring and ranking performances in different structure types that mimic real docking applications. The scoring powers of Δ_{Lin_F9}XGB for locally optimized poses, flexible redocked poses, and ensemble docked poses of the CASF-2016 core set achieve Pearson's correlation coefficient (R) values of 0.853, 0.839, and 0.813, respectively. In addition, the large-scale docking-based virtual screening test on the LIT-PCBA data set demonstrates the reliability and robustness of Δ_{Lin_F9}XGB in virtual screening application. The Δ_{Lin_F9}XGB scoring function and its code are freely available on the web at (https://yzhang.hpc.nyu.edu/Delta_LinF9_XGB).

PharmRF: A machine-learning scoring function to identify the best protein-ligand complexes for structure-based pharmacophore screening with high enrichments

Structure-based pharmacophore models are often developed by selecting a single protein-ligand complex with good resolution and better binding affinity data which prevents the analysis of other structures having a similar potential to act as better templates. PharmRF is a pharmacophore-based scoring function for selecting the best crystal structures with the potential to attain high enrichment rates in pharmacophore-based virtual screening prospectively. The PharmRF scoring function is trained and tested on the PDBbind v2018 protein-ligand complex dataset and employs a random forest regressor to correlate protein pocket descriptors and ligand pharmacophoric elements with binding affinity. PharmRF score represents the calculated binding affinity which identifies high-affinity ligands by thorough pruning of all the PDB entries available for a particular protein of interest with a high PharmRF score. Ligands with high PharmRF scores can provide a better basis for structure-based pharmacophore enumerations with a better enrichment rate. Evaluated on 10 protein-ligand systems of the DUD-E dataset, PharmRF achieved superior performance (average success rate: 77.61%, median success rate: 87.16%) than Vina docking score (75.47%, 79.39%). PharmRF was further evaluated using the CASF-2016 benchmark set yielding a moderate correlation of 0.591 with experimental binding affinity, similar in performance to 25 scoring functions tested on this dataset. Independent assessment of PharmRF on 8 protein-ligand systems of LIT-PCBA dataset exhibited average and median success rates of 57.55% and 74.72% with 4 targets attaining success rate > 90%. The PharmRF scoring model, scripts, and related resources can be accessed at https://github.com/Prasanth-Kumar87/PharmRF.

Polypharmacology: The science of multi-targeting molecules

Polypharmacology is a concept where a molecule can interact with two or more targets simultaneously. It offers many advantages as compared to the conventional single-targeting molecules. A multi-targeting drug is much more efficacious due to its cumulative efficacy at all of its individual targets making it much more effective in complex and multifactorial diseases like cancer, where multiple proteins and pathways are involved in the onset and development of the disease. For a molecule to be polypharmacologic in nature, it needs to possess promiscuity which is the ability to interact with multiple targets; and at the same time avoid binding to antitargets which would otherwise result in off-target adverse effects. There are certain structural features and physicochemical properties which when present would help researchers to predict if the designed molecule would possess promiscuity or not. Promiscuity can also be identified via advanced state-of-the-art computational methods. In this review, we also elaborate on the methods by which one can intentionally incorporate promiscuity in their molecules and make them polypharmacologic. The polypharmacology paradigm of "one drug-multiple targets" has numerous applications especially in drug repurposing where an already established drug is redeveloped for a new indication. Though designing a polypharmacological drug is much more difficult than designing a single-targeting drug, with the current technologies and information regarding different diseases and chemical functional groups, it is plausible for researchers to intentionally design a polypharmacological drug and unlock its advantages.

Protein-ligand binding affinity prediction based on profiles of intermolecular contacts

As a key element in structure-based drug design, binding affinity prediction (BAP) for putative protein-ligand complexes can be efficiently achieved by the incorporation of structural descriptors and machine-learning models. However, developing concise descriptors that will lead to accurate and interpretable BAP remains a difficult problem in this field. Herein, we introduce the profiles of intermolecular contacts (IMCPs) as descriptors for machine-learning-based BAP. IMCPs describe each group of protein-ligand contacts by the count and average distance of the group members, and collaborate closely with classical machine-learning models. Performed on multiple validation sets, IMCP-based models often result in better BAP accuracy than those originating from other similar descriptors. Additionally, IMCPs are simple and concise, and easy to interpret in model training. These descriptors highly conclude the structural information of protein-ligand complexes and can be easily updated with personalized profile features. IMCPs have been implemented in the BAP Toolkit on github ( https://github.com/debbydanwang/BAP).

Structure-based protein-ligand interaction fingerprints for binding affinity prediction

Binding affinity prediction (BAP) using protein–ligand complex structures is crucial to computer-aided drug design, but remains a challenging problem. To achieve efficient and accurate BAP, machine-learning scoring functions (SFs) based on a wide range of descriptors have been developed. Among those descriptors, protein–ligand interaction fingerprints (IFPs) are competitive due to their simple representations, elaborate profiles of key interactions and easy collaborations with machine-learning algorithms. In this paper, we have adopted a building-block-based taxonomy to review a broad range of IFP models, and compared representative IFP-based SFs in target-specific and generic scoring tasks. Atom-pair-counts-based and substructure-based IFPs show great potential in these tasks.

Navigating Chemical Space by Interfacing Generative Artificial Intelligence and Molecular Docking

Here, we report the implementation and application of a simple, structure-aware framework to generate target-specific screening libraries. Our approach combines advances in generative artificial intelligence (AI) with conventional molecular docking to explore chemical space conditioned on the unique physicochemical properties of the active site of a biomolecular target. As a demonstration, we used our framework, which we refer to as sample-and-dock, to construct focused libraries for cyclin-dependent kinase type-2 (CDK2) and the active site of the main protease (M^pro) of the SARS-CoV-2 virus. We envision that the sample-and-dock framework could be used to generate theoretical maps of the chemical space specific to a given target and so provide information about its molecular recognition characteristics.

The Structural Determinants for α1-Adrenergic/Serotonin Receptors Activity among Phenylpiperazine-Hydantoin Derivatives

Several studies confirmed the reciprocal interactions between adrenergic and serotoninergic systems and the influence of these phenomena on the pathogenesis of anxiety. Hence, searching for chemical agents with a multifunctional pharmacodynamic profile may bring highly effective therapy for CNS disorders. This study presents a deep structural insight into the hydantoin-arylpiperazine group and their serotonin/α-adrenergic activity. The newly synthesized compounds were tested in the radioligand binding assay and the intrinsic activity was evaluated for the selected derivatives. The computer-aided SAR analysis enabled us to answer questions about the influence of particular structural fragments on selective vs. multifunctional activity. As a result of the performed investigations, there were two leading structures: (a) compound 12 with multifunctional adrenergic-serotonin activity, which is a promising candidate to be an effective anxiolytic agent; (b) compound 14 with high α_1A/α_1D affinity and selectivity towards α_1B, which is recommended due to the elimination of probable cardiotoxic effect. The structural conclusions of this work provide significant support for future lead optimization in order to achieve the desired pharmacodynamic profile in searching for new CNS-modulating agents.

2D SIFt: a matrix of ligand-receptor interactions

Depicting a ligand-receptor complex via Interaction Fingerprints has been shown to be both a viable data visualization and an analysis tool. The spectrum of its applications ranges from simple visualization of the binding site through analysis of molecular dynamics runs, to the evaluation of the homology models and virtual screening. Here we present a novel tool derived from the Structural Interaction Fingerprints providing a detailed and unique insight into the interactions between receptor and specific regions of the ligand (grouped into pharmacophore features) in the form of a matrix, a 2D-SIFt descriptor. The provided implementation is easy to use and extends the python library, allowing the generation of interaction matrices and their manipulation (reading and writing as well as producing the average 2D-SIFt). The library for handling the interaction matrices is available via repository http://bitbucket.org/zchl/sift2d.

Subtype-selective mechanisms of negative allosteric modulators binding to group I metabotropic glutamate receptors

Group I metabotropic glutamate receptors (mGlu₁ and mGlu₅) are promising targets for multiple psychiatric and neurodegenerative disorders. Understanding the subtype selectivity of mGlu₁ and mGlu₅ allosteric sites is essential for the rational design of novel modulators with single- or dual-target mechanism of action. In this study, starting from the deposited mGlu₁ and mGlu₅ crystal structures, we utilized computational modeling approaches integrating docking, molecular dynamics simulation, and efficient post-trajectory analysis to reveal the subtype-selective mechanism of mGlu₁ and mGlu₅ to 10 diverse drug scaffolds representing known negative allosteric modulators (NAMs) in the literature. The results of modeling identified six pairs of non-conserved residues and four pairs of conserved ones as critical features to distinguish the selective NAMs binding to the corresponding receptors. In addition, nine pairs of residues are beneficial to the development of novel dual-target NAMs of group I metabotropic glutamate receptors. Furthermore, the binding modes of a reported dual-target NAM (VU0467558) in mGlu₁ and mGlu₅ were predicted to verify the identified residues that play key roles in the receptor selectivity and the dual-target binding. The results of this study can guide rational structure-based design of novel NAMs, and the approach can be generally applicable to characterize the features of selectivity for other G-protein-coupled receptors.

MolADI: A Web Server for Automatic Analysis of Protein-Small Molecule Dynamic Interactions

Protein-ligand interaction analysis is important for drug discovery and rational protein design. The existing online tools adopt only a single conformation of the complex structure for calculating and displaying the interactions, whereas both protein residues and ligand molecules are flexible to some extent. The interactions evolved with time in the trajectories are of greater interest. MolADI is a user-friendly online tool which analyzes the protein-ligand interactions in detail for either a single structure or a trajectory. Interactions can be viewed easily with both 2D graphs and 3D representations. MolADI is available as a web application.

Ligand Binding Site Comparison - LiBiSCo - a web-based tool for analyzing interactions between proteins and ligands to explore amino acid specificity within active sites

Interaction between protein and ligands are ubiquitous in a biological cell, and understanding these interactions at the atom level in protein-ligand complexes is crucial for structural bioinformatics and drug discovery. Here, we present a web-based protein-ligand interaction application named Ligand Binding Site Comparison (LiBiSCo) for comparing the amino acid residues interacting with atoms of a ligand molecule between different protein-ligand complexes available in the Protein Data Bank (PDB) database. The comparison is performed at the ligand atom level irrespectively of having binding site similarity or not between the protein structures of interest. The input used in LiBiSCo is one or several PDB IDs of protein-ligand complex(es) and the tool returns a list of identified interactions at ligand atom level including both bonded and non-bonded interactions. A sequence profile for the interaction for each ligand atoms is provided as a WebLogo. The LiBiSco is useful in understanding ligand binding specificity and structural promiscuity among families that are structurally unrelated. The LiBiSCo tool can be accessed through https://albiorix.bioenv.gu.se/LiBiSCo/HomePage.py.

Structure-Activity Relationship Analysis of Cocrystallized Gliptin-like Pyrrolidine, Trifluorophenyl, and Pyrimidine-2,4-Dione Dipeptidyl Peptidase-4 Inhibitors

Approved and potent reported dipeptidyl peptidase-4 (DPP-4) inhibitors with gliptin-like structures are classified here according to their structures and mechanisms of the inhibition in three groups: (i) those with pyrrolidine or analogs as P1 fragment with α-aminoacyl linker, (ii) structures with trifluorophenyl moiety or analogs as P1 fragment with β-aminobutanoyl linker, and (iii) DPP-4 inhibitors with pyrimidine-2,4-dione or analogs as P1' fragment. The structure-activity relationship analysis was performed for those whose cocrystallized structures with the enzyme were published. While inhibitors with pyrrolidine and trifluorophenyl moiety or analogs as P1 fragment bind in a similar way in S1, S2 and S2 extensive domains of the enzyme, the binding mode of pyrimidine-2,4-dione derivatives/analogs differs with additional interactions in S1' and S2' pockets. Three general schemes of fragmented gliptins and gliptin-like structures with the enzyme and protein-ligand interaction fingerprints were made, which might be useful in the creation of DPP-4 inhibitor's design strategies.

Exploring the stability of inhibitor binding to SIK2 using molecular dynamics simulation and binding free energy calculation

Salt inducible kinase 2 (SIK2) is a calcium/calmodulin-dependent protein kinase-like kinase that is implicated in a variety of biological phenomena, including cellular metabolism, growth, and apoptosis. SIK2 is the key target for various cancers, including ovarian, breast, prostate, and lung cancers. Although potent inhibitors of SIK2 are being developed, their binding stability and functional role are not presently known. In this work, we studied the detailed interactions between SIK2 and four of its inhibitors, HG-9-91-01, KIN112, MRT67307, and MRT199665, using molecular docking, molecular dynamics simulation, binding free energy calculation, and interaction fingerprint analysis. Intermolecular interactions revealed that HG-9-91-01 and KIN112 have stronger interactions with SIK2 than those of MRT199665 and MRT67307. The key residues involved in binding with SIK2 are conserved among all four inhibitors. Our results explain the detailed interaction of SIK2 with its inhibitors at the molecular level, thus paving the way for the development of targeted efficient anti-cancer drugs.

True Accuracy of Fast Scoring Functions to Predict High-Throughput Screening Data from Docking Poses: The Simpler the Better

Hundreds of fast scoring functions have been developed over the last 20 years to predict binding free energies from three-dimensional structures of protein-ligand complexes. Despite numerous statistical promises, we believe that none of them has been properly validated for daily prospective high-throughput virtual screening studies, mostly because in silico screening challenges usually employ artificially built and biased datasets. We here carry out a fully unbiased evaluation of four scoring functions (Pafnucy, Δ_vinaRF₂₀, IFP, and GRIM) on an in-house developed data collection of experimental high-confidence screening data (LIT-PCBA) covering about 3 million data points on 15 diverse pharmaceutical targets. All four scoring functions were applied to rescore the docking poses of LIT-PCBA compounds in conditions mimicking exactly standard drug discovery scenarios and were compared in terms of propensity to enrich true binders in the top 1%-ranked hit lists. Interestingly, rescoring based on simple interaction fingerprints or interaction graphs outperforms state-of-the-art machine learning and deep learning scoring functions in most of the cases. The current study notably highlights the strong tendency of deep learning methods to predict affinity values within a very narrow range centered on the mean value of samples used for training. Moreover, it suggests that knowledge of pre-existing binding modes is the key to detecting the most potent binders.

Computational representations of protein-ligand interfaces for structure-based virtual screening

Introduction: Structure-based virtual screening (SBVS) is an essential strategy for hit identification. SBVS primarily uses molecular docking, which exploits the protein-ligand binding mode and associated affinity score for compound ranking. Previous studies have shown that computational representation of protein-ligand interfaces and the later establishment of machine learning models are efficacious in improving the accuracy of SBVS.Areas covered: The authors review the computational methods for representing protein-ligand interfaces, which include the traditional ones that use deliberately designed fingerprints and descriptors and the more recent methods that automatically extract features with deep learning. The effects of these methods on the performance of machine learning models are briefly discussed. Additionally, case studies that applied various computational representations to machine learning are cited with remarks.Expert opinion: It has become a trend to extract binding features automatically by deep learning, which uses a completely end-to-end representation. However, there is still plenty of scope for improvement . The interpretability of deep-learning models, the organization of data management, the quantity and quality of available data, and the optimization of hyperparameters could impact the accuracy of feature extraction. In addition, other important structural factors such as water molecules and protein flexibility should be considered.

Applications of Virtual Screening in Bioprospecting: Facts, Shifts, and Perspectives to Explore the Chemo-Structural Diversity of Natural Products

Natural products are continually explored in the development of new bioactive compounds with industrial applications, attracting the attention of scientific research efforts due to their pharmacophore-like structures, pharmacokinetic properties, and unique chemical space. The systematic search for natural sources to obtain valuable molecules to develop products with commercial value and industrial purposes remains the most challenging task in bioprospecting. Virtual screening strategies have innovated the discovery of novel bioactive molecules assessing in silico large compound libraries, favoring the analysis of their chemical space, pharmacodynamics, and their pharmacokinetic properties, thus leading to the reduction of financial efforts, infrastructure, and time involved in the process of discovering new chemical entities. Herein, we discuss the computational approaches and methods developed to explore the chemo-structural diversity of natural products, focusing on the main paradigms involved in the discovery and screening of bioactive compounds from natural sources, placing particular emphasis on artificial intelligence, cheminformatics methods, and big data analyses.

Deep Learning in Virtual Screening: Recent Applications and Developments

Drug discovery is a cost and time-intensive process that is often assisted by computational methods, such as virtual screening, to speed up and guide the design of new compounds. For many years, machine learning methods have been successfully applied in the context of computer-aided drug discovery. Recently, thanks to the rise of novel technologies as well as the increasing amount of available chemical and bioactivity data, deep learning has gained a tremendous impact in rational active compound discovery. Herein, recent applications and developments of machine learning, with a focus on deep learning, in virtual screening for active compound design are reviewed. This includes introducing different compound and protein encodings, deep learning techniques as well as frequently used bioactivity and benchmark data sets for model training and testing. Finally, the present state-of-the-art, including the current challenges and emerging problems, are examined and discussed.

Site2Vec: a reference frame invariant algorithm for vector embedding of protein–ligand binding sites

Protein–ligand interactions are one of the fundamental types of molecular interactions in living systems. Ligands are small molecules that interact with protein molecules at specific regions on their surfaces called binding sites. Binding sites would also determine ADMET properties of a drug molecule. Tasks such as assessment of protein functional similarity and detection of side effects of drugs need identification of similar binding sites of disparate proteins across diverse pathways. To this end, methods for computing similarities between binding sites are still evolving and is an active area of research even today. Machine learning methods for similarity assessment require feature descriptors of binding sites. Traditional methods based on hand engineered motifs and atomic configurations are not scalable across several thousands of sites. In this regard, deep neural network algorithms are now deployed which can capture very complex input feature space. However, one fundamental challenge in applying deep learning to structures of binding sites is the input representation and the reference frame. We report here a novel algorithm, Site2Vec, that derives reference frame invariant vector embedding of a protein–ligand binding site. The method is based on pairwise distances between representative points and chemical compositions in terms of constituent amino acids of a site. The vector embedding serves as a locality sensitive hash function for proximity queries and determining similar sites. The method has been the top performer with more than 95% quality scores in extensive benchmarking studies carried over 10 data sets and against 23 other site comparison methods in the field. The algorithm serves for high throughput processing and has been evaluated for stability with respect to reference frame shifts, coordinate perturbations and residue mutations. We also provide the method as a standalone executable and a web service hosted at (http://services.iittp.ac.in/bioinfo/home).

Potential repurposing of four FDA approved compounds with antiplasmodial activity identified through proteome scale computational drug discovery and in vitro assay

Malaria elimination can benefit from time and cost-efficient approaches for antimalarials such as drug repurposing. In this work, 796 DrugBank compounds were screened against 36 Plasmodium falciparum targets using QuickVina-W. Hits were selected after rescoring using GRaph Interaction Matching (GRIM) and ligand efficiency metrics: surface efficiency index (SEI), binding efficiency index (BEI) and lipophilic efficiency (LipE). They were further evaluated in Molecular dynamics (MD). Twenty-five protein-ligand complexes were finally retained from the 28,656 (36 × 796) dockings. Hit GRIM scores (0.58 to 0.78) showed their molecular interaction similarity to co-crystallized ligands. Minimum LipE (3), SEI (23) and BEI (7) were in at least acceptable thresholds for hits. Binding energies ranged from -6 to -11 kcal/mol. Ligands showed stability in MD simulation with good hydrogen bonding and favorable protein-ligand interactions energy (the poorest being -140.12 kcal/mol). In vitro testing showed 4 active compounds with two having IC50 values in the single-digit μM range.

Merging Ligand-Based and Structure-Based Methods in Drug Discovery: An Overview of Combined Virtual Screening Approaches

Virtual screening (VS) is an outstanding cornerstone in the drug discovery pipeline. A variety of computational approaches, which are generally classified as ligand-based (LB) and structure-based (SB) techniques, exploit key structural and physicochemical properties of ligands and targets to enable the screening of virtual libraries in the search of active compounds. Though LB and SB methods have found widespread application in the discovery of novel drug-like candidates, their complementary natures have stimulated continued efforts toward the development of hybrid strategies that combine LB and SB techniques, integrating them in a holistic computational framework that exploits the available information of both ligand and target to enhance the success of drug discovery projects. In this review, we analyze the main strategies and concepts that have emerged in the last years for defining hybrid LB + SB computational schemes in VS studies. Particularly, attention is focused on the combination of molecular similarity and docking, illustrating them with selected applications taken from the literature.

Benchmarking Data Sets from PubChem BioAssay Data: Current Scenario and Room for Improvement

Developing realistic data sets for evaluating virtual screening methods is a task that has been tackled by the cheminformatics community for many years. Numerous artificially constructed data collections were developed, such as DUD, DUD-E, or DEKOIS. However, they all suffer from multiple drawbacks, one of which is the absence of experimental results confirming the impotence of presumably inactive molecules, leading to possible false negatives in the ligand sets. In light of this problem, the PubChem BioAssay database, an open-access repository providing the bioactivity information of compounds that were already tested on a biological target, is now a recommended source for data set construction. Nevertheless, there exist several issues with the use of such data that need to be properly addressed. In this article, an overview of benchmarking data collections built upon experimental PubChem BioAssay input is provided, along with a thorough discussion of noteworthy issues that one must consider during the design of new ligand sets from this database. The points raised in this review are expected to guide future developments in this regard, in hopes of offering better evaluation tools for novel in silico screening procedures.

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.

Revealing Acquired Resistance Mechanisms of Kinase-Targeted Drugs Using an on-the-Fly, Function-Site Interaction Fingerprint Approach

Although kinase-targeted drugs have achieved significant clinical success, they are frequently subject to the limitations of drug resistance, which has become a primary vulnerability to targeted drug therapy. Therefore, deciphering resistance mechanisms is an important step in designing more efficacious, antiresistant drugs. Here we studied two FDA-approved kinase drugs: Crizotinib and Ceritinib, which are first- and second-generation anaplastic lymphoma kinase (ALK) targeted inhibitors, to unravel drug-resistance mechanisms. We used an on-the-fly, function-site interaction fingerprint (on-the-fly Fs-IFP) approach, combining binding free-energy surface calculations with the Fs-IFPs. Establishing the potentials of mean force and monitoring the atomic-scale protein-ligand interactions, before and after L1196M-induced drug resistance, revealed insights into drug-resistance/antiresistant mechanisms. Crizotinib prefers to bind the wild-type ALK kinase domain, whereas Ceritinib binds more favorably to the mutated ALK kinase domain, in agreement with experimental results. We determined that ALK kinase-drug interactions in the region of the front pocket are associated with drug resistance. Additionally, we find that the L1196M mutation does not simply alter the binding modes of inhibitors but also affects the flexibility of the entire ALK kinase domain. Our work provides an understanding of the mechanisms of ALK drug resistance, confirms the usefulness of the on-the-fly Fs-IFP approach, and provides a practical paradigm to study drug-resistance mechanisms in prospective drug discovery.

Drug repositioning or target repositioning: A structural perspective of drug-target-indication relationship for available repurposed drugs

Drug repositioning aims to find new indications for existing drugs in order to reduce drug development cost and time. Currently,there are numerous stories of successful drug repositioning that have been reported and many repurposed drugs are already available on the market. Although drug repositioning is often a product of serendipity, repositioning opportunities can be uncovered systematically. There are three systematic approaches to drug repositioning: disease-centric approach, target-centric and drug-centric. Disease-centric approaches identify close relationships between an old and a new indication. A target-centric approach links a known target and its established drug to a new indication. Lastly, a drug-centric approach connects a known drug to a new target and its associated indication. These three approaches differ in their potential and their limitations, but above all else, in the required start information and computing power. This raises the question of which approach prevails in current drug discovery and what that implies for future developments. To address this question, we systematically evaluated over 100 drugs, 200 target structures and over 300 indications from the Drug Repositioning Database. Each analyzed case was classified as one of the three repositioning approaches. For the majority of cases (more than 60%) the disease-centric definition was assigned. Almost 30% of the cases were classified as target-centric and less than 10% as drug-centric approaches. We concluded that, despite the use of umbrella term “drug” repositioning, disease- and target-centric approaches have dominated the field until now. We propose the use of drug-centric approaches while discussing reasons, such as structure-based repositioning techniques, to exploit the full potential of drug-target-disease connections.

AlphaSpace 2.0: Representing Concave Biomolecular Surfaces Using β-Clusters

Modern rational modulator design and structure-function characterization often concentrate on concave regions of biomolecular surfaces, ranging from well-defined small-molecule binding sites to large protein-protein interaction interfaces. Here, we introduce a β-cluster as a pseudomolecular representation of fragment-centric pockets detected by AlphaSpace [J. Chem. Inf. Model. 2015, 55, 1585], a recently developed computational analysis tool for topographical mapping of biomolecular concavities. By mimicking the shape as well as atomic details of potential molecular binders, this new β-cluster representation allows direct pocket-to-ligand shape comparison and can be used to guide ligand optimization. Furthermore, we defined the β-score, the optimal Vina score of the β-cluster, as an indicator of pocket ligandability and developed an ensemble β-cluster approach, which allows one-to-one pocket mapping and comparison among aligned protein structures. We demonstrated the utility of β-cluster representation by applying the approach to a wide variety of problems including binding site detection and comparison, characterization of protein-protein interactions, and fragment-based ligand optimization. These new β-cluster functionalities have been implemented in AlphaSpace 2.0, which is freely available on the web at http://www.nyu.edu/projects/yzhang/AlphaSpace2.

visGReMLIN: graph mining-based detection and visualization of conserved motifs at 3D protein-ligand interface at the atomic level

Background

Interactions between proteins and non-proteic small molecule ligands play important roles in the biological processes of living systems. Thus, the development of computational methods to support our understanding of the ligand-receptor recognition process is of fundamental importance since these methods are a major step towards ligand prediction, target identification, lead discovery, and more. This article presents visGReMLIN, a web server that couples a graph mining-based strategy to detect motifs at the protein-ligand interface with an interactive platform to visually explore and interpret these motifs in the context of protein-ligand interfaces.

Results

To illustrate the potential of visGReMLIN, we conducted two cases in which our strategy was compared with previous experimentally and computationally determined results. visGReMLIN allowed us to detect patterns previously documented in the literature in a totally visual manner. In addition, we found some motifs that we believe are relevant to protein-ligand interactions in the analyzed datasets.

Conclusions

We aimed to build a visual analytics-oriented web server to detect and visualize common motifs at the protein-ligand interface. visGReMLIN motifs can support users in gaining insights on the key atoms/residues responsible for protein-ligand interactions in a dataset of complexes.

The binding mode of vilazodone in the human serotonin transporter elucidated by ligand docking and molecular dynamics simulations

Vilazodone is a novel antidepressant used for the treatment of major depressive disorder (MDD) with a primary action mechanism of inhibiting the human serotonin reuptake transporter (hSERT) and acting as a 5-HT_1A receptor partial agonist. The interaction between vilazodone and the 5-HT_1A receptor has been reported, however, the binding mode of vilazodone in the hSERT remains elusive. In the current study, to elucidate the molecular mechanism of vilazodone binding in the hSERT, the drug and its five analogs were docked into the hSERT crystal structure as initial conformations and were sampled by 400 ns molecular dynamics (MD) simulations. Through the analysis of the profiles of protein-ligand binding free energies, interaction fingerprints, and conformational rearrangements, the binding mode of vilazodone in the hSERT was revealed. As a result, unlike the classical antidepressants located in the S1 site of the hSERT, vilazodone adopted a linear pose in the binding pocket. Its arylpiperazine fragment occupies the central site (S1) and interacts with Y95, D98, I172, Y176, F335, F341, S438, and T439, while the indole fragment extends to the allosteric site (S2) via interacting with the ionic switch (R104/E403) between the two sites. The new insights obtained are not only helpful in understanding the binding mode of vilazodone in the hSERT, but also provide valuable guidance to the discovery of novel antidepressant drugs.

Structure-based drug repositioning: Potential and limits

Drug repositioning, the assignment of new therapeutic purposes to known drugs, is an established strategy with many repurposed drugs on the market and many more at experimental stage. We review three use cases, a herpes drug with benefits in cancer, a cancer drug with potential in autoimmune disease, and a selective and an unspecific drug binding the same target (GPCR). We explore these use cases from a structural point of view focusing on a deep understanding of the underlying drug-target interactions. We review tools and data needed for such a drug-centric structural repositioning approach. Finally, we show that the availability of data on targets is an important limiting factor to realize the full potential of structural drug-repositioning.

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.

Targeting the Side-Chain Convergence of Hydrophobic α-Helical Hot Spots To Design Small-Molecule Mimetics: Key Binding Features for i, i + 3, and i + 7

The conformational convergence of hydrophobic α-helical hot spots was revealed by analyzing α-helix-mediated protein-protein interaction (PPI) complex structures. The pharmacophore models were derived for hydrophobic α-helical hot spots at positions i, i + 3, and i + 7. These provide the foundation for designing generalizable scaffolds that can directly mimic the binding mode of the side chains of α-helical hot spots, offering a new class of small-molecule α-helix mimetics. For the first time, the protocol was developed to identify the PPI targets that have similar binding pockets, allowing evaluation of inhibitor selectivities between α-helix-mediated PPIs. The mimicry efficiency of the previously designed scaffold 1 was disclosed. The close positioning of this small molecule to the additional α-helical hot spots suggests that the decoration of this series of generalizable scaffolds can conveniently reach the binding pockets of additional α-helical hot spots to produce potent small-molecule inhibitors for α-helix-mediated PPIs.