Function-specific virtual screening for GPCR ligands using a combined scoring method

AlphaFold accelerated discovery of psychotropic agonists targeting the trace amine-associated receptor 1

Artificial intelligence is revolutionizing protein structure prediction, providing unprecedented opportunities for drug design. To assess the potential impact on ligand discovery, we compared virtual screens using protein structures generated by the AlphaFold machine learning method and traditional homology modeling. More than 16 million compounds were docked to models of the trace amine-associated receptor 1 (TAAR1), a G protein-coupled receptor of unknown structure and target for treating neuropsychiatric disorders. Sets of 30 and 32 highly ranked compounds from the AlphaFold and homology model screens, respectively, were experimentally evaluated. Of these, 25 were TAAR1 agonists with potencies ranging from 12 to 0.03 μM. The AlphaFold screen yielded a more than twofold higher hit rate (60%) than the homology model and discovered the most potent agonists. A TAAR1 agonist with a promising selectivity profile and drug-like properties showed physiological and antipsychotic-like effects in wild-type but not in TAAR1 knockout mice. These results demonstrate that AlphaFold structures can accelerate drug discovery.

Know your molecule: pharmacological characterization of drug candidates to enhance efficacy and reduce late-stage attrition

Despite advances in chemical, computational and biological sciences, the rate of attrition of drug candidates in clinical development is still high. A key point in the small-molecule discovery process that could provide opportunities to help address this challenge is the pharmacological characterization of hit and lead compounds, culminating in the selection of a drug candidate. Deeper characterization is increasingly important, because the 'quality' of drug efficacy, at least for G protein-coupled receptors (GPCRs), is now understood to be much more than activation of commonly evaluated pathways such as cAMP signalling, with many more 'efficacies' of ligands that could be harnessed therapeutically. Such characterization is being enabled by novel assays to characterize the complex behaviour of GPCRs, such as biased signalling and allosteric modulation, as well as advances in structural biology, such as cryo-electron microscopy. This article discusses key factors in the assessments of the pharmacology of hit and lead compounds in the context of GPCRs as a target class, highlighting opportunities to identify drug candidates with the potential to address limitations of current therapies and to improve the probability of them succeeding in clinical development.

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.

Computational drug development for membrane protein targets

The application of computational biology in drug development for membrane protein targets has experienced a boost from recent developments in deep learning-driven structure prediction, increased speed and resolution of structure elucidation, machine learning structure-based design and the evaluation of big data. Recent protein structure predictions based on machine learning tools have delivered surprisingly reliable results for water-soluble and membrane proteins but have limitations for development of drugs that target membrane proteins. Structural transitions of membrane proteins have a central role during transmembrane signaling and are often influenced by therapeutic compounds. Resolving the structural and functional basis of dynamic transmembrane signaling networks, especially within the native membrane or cellular environment, remains a central challenge for drug development. Tackling this challenge will require an interplay between experimental and computational tools, such as super-resolution optical microscopy for quantification of the molecular interactions of cellular signaling networks and their modulation by potential drugs, cryo-electron microscopy for determination of the structural transitions of proteins in native cell membranes and entire cells, and computational tools for data analysis and prediction of the structure and function of cellular signaling networks, as well as generation of promising drug candidates.

Pharmacological Characterization of the Zebrafish (Danio Rerio) Histamine H1 Receptor Reveals the Involvement of the Second Extracellular Loop in the Binding of Histamine

The zebrafish (Danio rerio) histamine H1 receptor gene (zfH1R) was cloned in 2007 and reported to be involved in fish locomotion. Yet, no detailed characterization of its pharmacology and signaling properties have so far been reported. In this study, we pharmacologically characterized the zfH1R expressed in HEK-293T cells by means of [3H]-mepyramine binding and G protein-signaling assays. The zfH1R [dissociation constant (KD), 0.7 nM] displayed similar affinity for the antagonist [3H]-mepyramine as the human histamine H1 receptor (hH1R) (KD, 1.5 nM), whereas the affinity for histamine is 100-fold higher than for the human H1R. The zfH1R couples to Gαq/11 proteins and activates several reporter genes, i.e., NFAT, NFϰB, CRE, VEGF, COX-2, SRE, and AP-1, and zfH1R-mediated signaling is prevented by the Gαq/11 inhibitor YM-254890 and the antagonist mepyramine. Molecular modeling of the zfH1R and human H1R shows that the binding pockets are identical, implying that variations along the ligand binding pathway could underly the differences in histamine affinity instead. Targeting differentially charged residues in extracellular loop 2 (ECL2) using site-directed mutagenesis revealed that Arg21045x55 is most likely involved in the binding process of histamine in zfH1R. This study aids the understanding of the pharmacological differences between H1R orthologs and the role of ECL2 in histamine binding and provides fundamental information for the understanding of the histaminergic system in the zebrafish. SIGNIFICANCE STATEMENT: The use of the zebrafish as in vivo models in neuroscience is growing exponentially, which asks for detailed characterization of the aminergic neurotransmitter systems in this model. This study is the first to pharmacologically characterize the zebrafish histamine H1 receptor after expression in HEK-293T cells. The results show a high pharmacological and functional resemblance with the human ortholog but also reveal interesting structural differences and unveils an important role of the second extracellular loop in histamine binding.

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.

Combined docking and machine learning identify key molecular determinants of ligand pharmacological activity on β2 adrenoceptor

G protein‐coupled receptors (GPCRs) are valuable therapeutic targets for many diseases. A central question of GPCR drug discovery is to understand what determines the agonism or antagonism of ligands that bind them. Ligands exert their action via the interactions in the ligand binding pocket. We hypothesized that there is a common set of receptor interactions made by ligands of diverse structures that mediate their action and that among a large dataset of different ligands, the functionally important interactions will be over‐represented. We computationally docked ~2700 known β2AR ligands to multiple β2AR structures, generating ca 75 000 docking poses and predicted all atomic interactions between the receptor and the ligand. We used machine learning (ML) techniques to identify specific interactions that correlate with the agonist or antagonist activity of these ligands. We demonstrate with the application of ML methods that it is possible to identify the key interactions associated with agonism or antagonism of ligands. The most representative interactions for agonist ligands involve K97^2.68×67, F194^ECL2, S203^5.42×43, S204^5.43×44, S207^5.46×641, H296^6.58×58, and K305^7.32×31. Meanwhile, the antagonist ligands made interactions with W286^6.48×48 and Y316^7.43×42, both residues considered to be important in GPCR activation. The interpretation of ML analysis in human understandable form allowed us to construct an exquisitely detailed structure‐activity relationship that identifies small changes to the ligands that invert their pharmacological activity and thus helps to guide the drug discovery process. This approach can be readily applied to any drug target.

Optimization of Peptide Linker-Based Fluorescent Ligands for the Histamine H1 Receptor

The histamine H₁ receptor (H₁R) has recently been implicated in mediating cell proliferation and cancer progression; therefore, high-affinity H₁R-selective fluorescent ligands are desirable tools for further investigation of this behavior in vitro and in vivo. We previously reported a H₁R fluorescent ligand, bearing a peptide-linker, based on antagonist VUF13816 and sought to further explore structure-activity relationships (SARs) around the linker, orthostere, and fluorescent moieties. Here, we report a series of high-affinity H₁R fluorescent ligands varying in peptide linker composition, orthosteric targeting moiety, and fluorophore. Incorporation of a boron-dipyrromethene (BODIPY) 630/650-based fluorophore conferred high binding affinity to our H₁R fluorescent ligands, remarkably overriding the linker SAR observed in corresponding unlabeled congeners. Compound 31a, both potent and subtype-selective, enabled H₁R visualization using confocal microscopy at a concentration of 10 nM. Molecular docking of 31a with the human H₁R predicts that the optimized peptide linker makes interactions with key residues in the receptor.

Consensus scoring evaluated using the GPCR-Bench dataset: Reconsidering the role of MM/GBSA

The recent availability of large numbers of GPCR crystal structures has provided an unprecedented opportunity to evaluate their performance in virtual screening protocols using established benchmarking datasets. In this study, we evaluated the ability of MM/GBSA in consensus scoring-based virtual screening enrichment together with nine classical scoring functions, using the GPCR-Bench dataset consisting of 24 GPCR crystal structures and 254,646 actives and decoys. While the performance of consensus scoring was modest overall, combinations which included MM/GBSA performed relatively well compared to combinations of classical scoring functions. Combinations of MM/GBSA and good-performing scoring functions provided the highest proportion of improvements, with improvements observed in 32% and 19% of all combinations across all targets at the EF_1% and EF_5% levels respectively. Combinations of MM/GBSA and poor-performing scoring functions still outperformed classical scoring functions, with improvements observed in 26% and 17% of all combinations at the EF_1% and EF_5% levels. In comparison, only 14-22% and 6-11% of combinations of classical scoring functions produced improvements at EF_1% and EF_5% respectively. Efforts to improve performance by increasing the number of scoring functions in consensus scoring to three were mostly ineffective. We also observed that consensus scoring performed better for individual scoring functions possessing initially low enrichment factors, potentially implying their benefits are more relevant in such scenarios. Overall, this study demonstrated the first implementation of MM/GBSA in consensus scoring using the GPCR-Bench dataset and could provide a valuable benchmark of the performance of MM/GBSA in comparison to classical scoring functions in consensus scoring for GPCRs.

Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

G protein-coupled receptors (GPCRs) constitute the largest family of membrane proteins in the human genome and are important therapeutic targets. During the last decade, the number of atomic-resolution structures of GPCRs has increased rapidly, providing insights into drug binding at the molecular level. These breakthroughs have created excitement regarding the potential of using structural information in ligand design and initiated a new era of rational drug discovery for GPCRs. The molecular docking method is now widely applied to model the three-dimensional structures of GPCR-ligand complexes and screen for chemical probes in large compound libraries. In this review article, we first summarize the current structural coverage of the GPCR superfamily and the understanding of receptor-ligand interactions at atomic resolution. We then present the general workflow of structure-based virtual screening and strategies to discover GPCR ligands in chemical libraries. We assess the state of the art of this research field by summarizing prospective applications of virtual screening based on experimental structures. Strategies to identify compounds with specific efficacy and selectivity profiles are discussed, illustrating the opportunities and limitations of the molecular docking method. Our overview shows that structure-based virtual screening can discover novel leads and will be essential in pursuing the next generation of GPCR drugs. SIGNIFICANCE STATEMENT: Extraordinary advances in the structural biology of G protein-coupled receptors have revealed the molecular details of ligand recognition by this large family of therapeutic targets, providing novel avenues for rational drug design. Structure-based docking is an efficient computational approach to identify novel chemical probes from large compound libraries, which has the potential to accelerate the development of drug candidates.

From structure to clinic: Design of a muscarinic M1 receptor agonist with potential to treatment of Alzheimer's disease

Current therapies for Alzheimer's disease seek to correct for defective cholinergic transmission by preventing the breakdown of acetylcholine through inhibition of acetylcholinesterase, these however have limited clinical efficacy. An alternative approach is to directly activate cholinergic receptors responsible for learning and memory. The M1-muscarinic acetylcholine (M1) receptor is the target of choice but has been hampered by adverse effects. Here we aimed to design the drug properties needed for a well-tolerated M1-agonist with the potential to alleviate cognitive loss by taking a stepwise translational approach from atomic structure, cell/tissue-based assays, evaluation in preclinical species, clinical safety testing, and finally establishing activity in memory centers in humans. Through this approach, we rationally designed the optimal properties, including selectivity and partial agonism, into HTL9936-a potential candidate for the treatment of memory loss in Alzheimer's disease. More broadly, this demonstrates a strategy for targeting difficult GPCR targets from structure to clinic.

Nonantimicrobial Actions of Macrolides: Overview and Perspectives for Future Development

Macrolides are among the most widely prescribed broad spectrum antibacterials, particularly for respiratory infections. It is now recognized that these drugs, in particular azithromycin, also exert time-dependent immunomodulatory actions that contribute to their therapeutic benefit in both infectious and other chronic inflammatory diseases. Their increased chronic use in airway inflammation and, more recently, of azithromycin in COVID-19, however, has led to a rise in bacterial resistance. An additional crucial aspect of chronic airway inflammation, such as chronic obstructive pulmonary disease, as well as other inflammatory disorders, is the loss of epithelial barrier protection against pathogens and pollutants. In recent years, azithromycin has been shown with time to enhance the barrier properties of airway epithelial cells, an action that makes an important contribution to its therapeutic efficacy. In this article, we review the background and evidence for various immunomodulatory and time-dependent actions of macrolides on inflammatory processes and on the epithelium and highlight novel nonantibacterial macrolides that are being studied for immunomodulatory and barrier-strengthening properties to circumvent the risk of bacterial resistance that occurs with macrolide antibacterials. We also briefly review the clinical effects of macrolides in respiratory and other inflammatory diseases associated with epithelial injury and propose that the beneficial epithelial effects of nonantibacterial azithromycin derivatives in chronic inflammation, even given prophylactically, are likely to gain increasing attention in the future. SIGNIFICANCE STATEMENT: Based on its immunomodulatory properties and ability to enhance the protective role of the lung epithelium against pathogens, azithromycin has proven superior to other macrolides in treating chronic respiratory inflammation. A nonantibiotic azithromycin derivative is likely to offer prophylactic benefits against inflammation and epithelial damage of differing causes while preserving the use of macrolides as antibiotics.

RealVS: Toward Enhancing the Precision of Top Hits in Ligand-Based Virtual Screening of Drug Leads from Large Compound Databases

Accurate modeling of compound bioactivities is essential for the virtual screening of drug leads. In real-world scenarios, pharmacists tend to choose from the top-k hit compounds ranked by predicted bioactivities from a large database with interest to continue wet experiments for drug discovery. Significant improvement of the precision of the top hits in ligand-based virtual screening of drug leads is more valuable than conventional schemes for accurately predicting the bioactivities of all compounds from a large database. Here, we proposed a new method, RealVS, to significantly improve the top hits' precision and learn interpretable key substructures associated with compound bioactivities. The features of RealVS involve the following points. (1) Abundant transferable information from the source domain was introduced for alleviating the insufficiency of inactive ligands associated with drug targets. (2) The adversarial domain alignment was adopted to fit the distribution of generated features of compounds from the training data set and that from the screening database for greater model generalization ability. (3) A novel objective function was proposed to simultaneously optimize the classification loss, regression loss, and adversarial loss, where most inactive ligands tend to be screened out before activity regression prediction. (4) Graph attention networks were adopted for learning key substructures associated with ligand bioactivities for better model interpretability. The results on a large number of benchmark data sets show that our method has significantly improved the precision of top hits under various k values in ligand-based virtual screening of drug leads from large compound databases, which is of great value in real-world scenarios. The web server of RealVS is freely available at noveldelta.com/RealVS for academic purposes, where virtual screening of hits from large compound databases is accessible.

Main protease inhibitors and drug surface hotspots for the treatment of COVID-19: A drug repurposing and molecular docking approach

Here, drug repurposing and molecular docking were employed to screen approved MPP inhibitors and their derivatives to suggest a specific therapeutic agent for the treatment of COVID-19. The approved MPP inhibitors against HIV and HCV were prioritized, while RNA dependent RNA Polymerase (RdRp) inhibitor remdesivir including Favipiravir, alpha-ketoamide were studied as control groups. The target drug surface hotspot was also investigated through the molecular docking technique. Molecular dynamics was performed to determine the binding stability of docked complexes. Absorption, distribution, metabolism, and excretion analysis was conducted to understand the pharmacokinetics and drug-likeness of the screened MPP inhibitors. The results of the study revealed that Paritaprevir (-10.9 kcal/mol) and its analog (CID 131982844) (-16.3 kcal/mol) showed better binding affinity than the approved MPP inhibitors compared in this study, including remdesivir, Favipiravir, and alpha-ketoamide. A comparative study among the screened putative MPP inhibitors revealed that the amino acids T25, T26, H41, M49, L141, N142, G143, C145, H164, M165, E166, D187, R188, and Q189 are at potentially critical positions for being surface hotspots in the MPP of SARS-CoV-2. The top 5 predicted drugs (Paritaprevir, Glecaprevir, Nelfinavir, and Lopinavir) and the topmost analog showed conformational stability in the active site of the SARS-CoV-2 MP protein. The study also suggested that Paritaprevir and its analog (CID 131982844) might be effective against SARS-CoV-2. The current findings are limited to in silico analysis and lack in vivo efficacy testing; thus, we strongly recommend a quick assessment of Paritaprevir and its analog (CID 131982844) in a clinical trial.

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

Ligands of Adrenergic Receptors: A Structural Point of View

Adrenergic receptors are G protein-coupled receptors for epinephrine and norepinephrine. They are targets of many drugs for various conditions, including treatment of hypertension, hypotension, and asthma. Adrenergic receptors are intensively studied in structural biology, displayed for binding poses of different types of ligands. Here, we summarized molecular mechanisms of ligand recognition and receptor activation exhibited by structure. We also reviewed recent advances in structure-based ligand discovery against adrenergic receptors.

Discovery of Novel Allosteric Modulators Targeting an Extra-Helical Binding Site of GLP-1R Using Structure- and Ligand-Based Virtual Screening

Allosteric modulators have emerged with many potential pharmacological advantages as they do not compete the binding of agonist or antagonist to the orthosteric sites but ultimately affect downstream signaling. To identify allosteric modulators targeting an extra-helical binding site of the glucagon-like peptide-1 receptor (GLP-1R) within the membrane environment, the following two computational approaches were applied: structure-based virtual screening with consideration of lipid contacts and ligand-based virtual screening with the maintenance of specific allosteric pocket residue interactions. Verified by radiolabeled ligand binding and cAMP accumulation experiments, two negative allosteric modulators and seven positive allosteric modulators were discovered using structure-based and ligand-based virtual screening methods, respectively. The computational approach presented here could possibly be used to discover allosteric modulators of other G protein-coupled receptors.

Comparison of structure- and ligand-based scoring functions for deep generative models: a GPCR case study

Deep generative models have shown the ability to devise both valid and novel chemistry, which could significantly accelerate the identification of bioactive compounds. Many current models, however, use molecular descriptors or ligand-based predictive methods to guide molecule generation towards a desirable property space. This restricts their application to relatively data-rich targets, neglecting those where little data is available to sufficiently train a predictor. Moreover, ligand-based approaches often bias molecule generation towards previously established chemical space, thereby limiting their ability to identify truly novel chemotypes. In this work, we assess the ability of using molecular docking via Glide-a structure-based approach-as a scoring function to guide the deep generative model REINVENT and compare model performance and behaviour to a ligand-based scoring function. Additionally, we modify the previously published MOSES benchmarking dataset to remove any induced bias towards non-protonatable groups. We also propose a new metric to measure dataset diversity, which is less confounded by the distribution of heavy atom count than the commonly used internal diversity metric. With respect to the main findings, we found that when optimizing the docking score against DRD2, the model improves predicted ligand affinity beyond that of known DRD2 active molecules. In addition, generated molecules occupy complementary chemical and physicochemical space compared to the ligand-based approach, and novel physicochemical space compared to known DRD2 active molecules. Furthermore, the structure-based approach learns to generate molecules that satisfy crucial residue interactions, which is information only available when taking protein structure into account. Overall, this work demonstrates the advantage of using molecular docking to guide de novo molecule generation over ligand-based predictors with respect to predicted affinity, novelty, and the ability to identify key interactions between ligand and protein target. Practically, this approach has applications in early hit generation campaigns to enrich a virtual library towards a particular target, and also in novelty-focused projects, where de novo molecule generation either has no prior ligand knowledge available or should not be biased by it.

Molecular Modeling of Histamine Receptors-Recent Advances in Drug Discovery

The recent developments of fast reliable docking, virtual screening and other algorithms gave rise to discovery of many novel ligands of histamine receptors that could be used for treatment of allergic inflammatory disorders, central nervous system pathologies, pain, cancer and obesity. Furthermore, the pharmacological profiles of ligands clearly indicate that these receptors may be considered as targets not only for selective but also for multi-target drugs that could be used for treatment of complex disorders such as Alzheimer's disease. Therefore, analysis of protein-ligand recognition in the binding site of histamine receptors and also other molecular targets has become a valuable tool in drug design toolkit. This review covers the period 2014-2020 in the field of theoretical investigations of histamine receptors mostly based on molecular modeling as well as the experimental characterization of novel ligands of these receptors.

New machine learning and physics-based scoring functions for drug discovery

Scoring functions are essential for modern in silico drug discovery. However, the accurate prediction of binding affinity by scoring functions remains a challenging task. The performance of scoring functions is very heterogeneous across different target classes. Scoring functions based on precise physics-based descriptors better representing protein–ligand recognition process are strongly needed. We developed a set of new empirical scoring functions, named DockTScore, by explicitly accounting for physics-based terms combined with machine learning. Target-specific scoring functions were developed for two important drug targets, proteases and protein–protein interactions, representing an original class of molecules for drug discovery. Multiple linear regression (MLR), support vector machine and random forest algorithms were employed to derive general and target-specific scoring functions involving optimized MMFF94S force-field terms, solvation and lipophilic interactions terms, and an improved term accounting for ligand torsional entropy contribution to ligand binding. DockTScore scoring functions demonstrated to be competitive with the current best-evaluated scoring functions in terms of binding energy prediction and ranking on four DUD-E datasets and will be useful for in silico drug design for diverse proteins as well as for specific targets such as proteases and protein–protein interactions. Currently, the MLR DockTScore is available at www.dockthor.lncc.br.

G protein-coupled receptors: structure- and function-based drug discovery

As one of the most successful therapeutic target families, G protein-coupled receptors (GPCRs) have experienced a transformation from random ligand screening to knowledge-driven drug design. We are eye-witnessing tremendous progresses made recently in the understanding of their structure-function relationships that facilitated drug development at an unprecedented pace. This article intends to provide a comprehensive overview of this important field to a broader readership that shares some common interests in drug discovery.

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.

Consensus Docking in Drug Discovery

Background:

Molecular docking is probably the most popular and profitable approach in computer-aided drug design, being the staple technique for predicting the binding mode of bioactive compounds and for performing receptor-based virtual screening studies. The growing attention received by docking, as well as the need for improving its reliability in pose prediction and virtual screening performance, has led to the development of a wide plethora of new docking algorithms and scoring functions. Nevertheless, it is unlikely to identify a single procedure outperforming the other ones in terms of reliability and accuracy or demonstrating to be generally suitable for all kinds of protein targets.

Methods:

In this context, consensus docking approaches are taking hold in computer-aided drug design. These computational protocols consist in docking ligands using multiple docking methods and then comparing the binding poses predicted for the same ligand by the different methods. This analysis is usually carried out calculating the root-mean-square deviation among the different docking results obtained for each ligand, in order to identify the number of docking methods producing the same binding pose.

Results:

The consensus docking approaches demonstrated to improve the quality of docking and virtual screening results compared to the single docking methods. From a qualitative point of view, the improvement in pose prediction accuracy was obtained by prioritizing ligand binding poses produced by a high number of docking methods, whereas with regards to virtual screening studies, high hit rates were obtained by prioritizing the compounds showing a high level of pose consensus.

Conclusion:

In this review, we provide an overview of the results obtained from the performance assessment of various consensus docking protocols and we illustrate successful case studies where consensus docking has been applied in virtual screening studies.

A Focus on Unusual ECL2 Interactions Yields β2 -Adrenergic Receptor Antagonists with Unprecedented Scaffolds

The binding pockets of aminergic G protein-coupled receptors are often targeted by drugs and virtual screening campaigns. In order to find ligands with unprecedented scaffolds for one of the best-investigated receptors of this subfamily, the β2 -adrenergic receptor, we conducted a docking-based screen insisting that molecules would address previously untargeted residues in extracellular loop 2. We here report the discovery of ligands with a previously undescribed coumaran-based scaffold. Furthermore, we provide an analysis of the added value that X-ray structures in different conformations deliver for such docking screens.

Application of MM-PBSA Methods in Virtual Screening

Computer-aided drug design techniques are today largely applied in medicinal chemistry. In particular, receptor-based virtual screening (VS) studies, in which molecular docking represents the gold standard in silico approach, constitute a powerful strategy for identifying novel hit compounds active against the desired target receptor. Nevertheless, the need for improving the ability of docking in discriminating true active ligands from inactive compounds, thus boosting VS hit rates, is still pressing. In this context, the use of binding free energy evaluation approaches can represent a profitable tool for rescoring ligand-protein complexes predicted by docking based on more reliable estimations of ligand-protein binding affinities than those obtained with simple scoring functions. In the present review, we focused our attention on the Molecular Mechanics-Poisson Boltzman Surface Area (MM-PBSA) method for the calculation of binding free energies and its application in VS studies. We provided examples of successful applications of this method in VS campaigns and evaluation studies in which the reliability of this approach has been assessed, thus providing useful guidelines for employing this approach in VS.

Developing an effective polarizable bond method for small molecules with application to optimized molecular docking

Electrostatic interaction plays an essential role in protein-ligand binding. Due to the polarization effect, electrostatic interactions are largely impacted by their local environments. However, traditional force fields use fixed point charge-charge interactions to describe electrostatic interactions but is unable to include the polarization effect. The lack of the polarization effect in the force field representation can result in substantial error in biomolecular studies, such as molecular dynamics and molecular docking. Docking programs usually employ traditional force fields to estimate the binding energy between a ligand and a protein for pose selection or scoring. The intermolecular interaction energy mainly consists of van der Waals and electrostatic interaction in the force field representation. In the current study, we developed an Effective Polarizable Bond (EPB) method for small organic molecules and applied this EPB method to optimize protein-ligand docking in computational tests for a variety of protein-ligand systems. We tested the method on a set of 38 cocrystallized structures taken from the Protein Data Bank (PDB) and found that the maximum error was reduced from 7.98 Å to 2.03 Å when using EPB Dock, providing strong evidence that the use of EPB charges is important. We found that our optimized docking approach with EPB charges could improve the docking performance, sometimes dramatically, and the maximum error was reduced from 12.88 Å to 1.57 Å in Optimized Docking (in the case of 1fqx). The average RMSD decreased from 2.83 Å to 1.85 Å. Further investigations showed that the use of the EBP method could enhance intermolecular hydrogen bonding, which is a major contributing factor to improved docking performance. Developed tools for the calculation of the polarized ligand charge from a protein-ligand complex structure with the EPB method are freely available on GitHub (https://github.com/Xundrug/EPB).

Exploring the molecular structures that confer ligand selectivity for galanin type II and III receptors

Galanin receptors (GALRs) belong to the superfamily of G-protein coupled receptors. The three GALR subtypes (GALR1, GALR2, and GALR3) are activated by their endogenous ligands: spexin (SPX) and galanin (GAL). The synthetic SPX-based GALR2-specific agonist, SG2A, plays a dual role in the regulation of appetite and depression-like behaviors. Little is known, however, about the molecular interaction between GALR2 and SG2A. Using site-directed mutagenesis and domain swapping between GALR2 and GALR3, we identified residues in GALR2 that promote interaction with SG2A and residues in GALR3 that inhibit interaction with SG2A. In particular, Phe¹⁰³, Phe¹⁰⁶, and His¹¹⁰ in the transmembrane helix 3 (TM3) domain; Val¹⁹³, Phe¹⁹⁴, and Ser¹⁹⁵ in the TM5 domain; and Leu²⁷³ in the extracellular loop 3 (ECL3) domain of GALR2 provide favorable interactions with the Asn⁵, Ala⁷, Phe¹¹, and Pro¹³ residues of SG2A. Our results explain how SG2A achieves selective interaction with GALR2 and inhibits interaction with GALR3. The results described here can be used broadly for in silico virtual screening of small molecules for the development of GALR subtype-specific agonists and/or antagonists.

Consensus Analyses in Molecular Docking Studies Applied to Medicinal Chemistry

The increasing number of computational studies in medicinal chemistry involving molecular docking has put the technique forward as promising in Computer-Aided Drug Design. Considering the main method in the virtual screening based on the structure, consensus analysis of docking has been applied in several studies to overcome limitations of algorithms of different programs and mainly to increase the reliability of the results and reduce the number of false positives. However, some consensus scoring strategies are difficult to apply and, in some cases, are not reliable due to the small number of datasets tested. Thus, for such a methodology to be successful, it is necessary to understand why, when and how to use consensus docking. Therefore, the present study aims to present different approaches to docking consensus, applications, and several scoring strategies that have been successful and can be applied in future studies.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

Fragment-Based Computational Method for Designing GPCR Ligands

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors, which is arguably the most important family of drug target. With the technology breakthroughs in X-ray crystallography and cryo-electron microscopy, more than 300 GPCR-ligand complex structures have been publicly reported since 2007, covering about 60 unique GPCRs. Such abundant structural information certainly will facilitate the structure-based drug design by targeting GPCRs. In this study, we have developed a fragment-based computational method for designing novel GPCR ligands. We first extracted the characteristic interaction patterns (CIPs) on the binding interfaces between GPCRs and their ligands. The CIPs were used as queries to search the chemical fragments derived from GPCR ligands, which were required to form similar interaction patterns with GPCR. Then, the selected chemical fragments were assembled into complete molecules by using the AutoT&T2 software. In this work, we chose β-adrenergic receptor (β-AR) and muscarinic acetylcholine receptor (mAChR) as the targets to validate this method. Based on the designs suggested by our method, samples of 63 compounds were purchased and tested in a cell-based functional assay. A total of 15 and 22 compounds were identified as active antagonists for β₂-AR and mAChR M1, respectively. Molecular dynamics simulations and binding free energy analysis were performed to explore the key interactions (e.g., hydrogen bonds and π-π interactions) between those active compounds and their target GPCRs. In summary, our work presents a useful approach to the de novo design of GPCR ligands based on the relevant 3D structural information.

Virtual screening identifies a PIN1 inhibitor with possible antiovarian cancer effects

Peptidyl-prolyl cis-trans isomerase, NIMA-interacting 1 (PIN1) is a peptidyl-prolyl isomerase that binds phospho-Ser/Thr-Pro motifs in proteins and catalyzes the cis-trans isomerization of proline peptide bonds. PIN1 is overexpressed in several cancers including high-grade serous ovarian cancer. Since few therapies are effective against this cancer, PIN1 could be a therapeutic target but effective PIN1 inhibitors are lacking. To identify molecules with in vivo inhibitory effects on PIN1, we used consensus docking to model existing PIN1-ligand X-ray structures and to screen a chemical database for candidate inhibitors. Ten molecules were selected and tested in cellular assays, leading to the identification of VS10 that bound and inhibited PIN1. VS10 treatment reduced the viability of ovarian cancer cell lines by inducing proteasomal PIN1 degradation, without effects on PIN1 transcription, and also reduced the levels of downstream targets β-catenin, cyclin D1, and pSer473-Akt. VS10 is a selective PIN1 inhibitor that may offer new opportunities for treating PIN1-overexpressing tumors.

Improving the Virtual Screening Ability of Target-Specific Scoring Functions Using Deep Learning Methods

Scoring functions play an important role in structure-based virtual screening. It has been widely accepted that target-specific scoring functions (TSSFs) may achieve better performance compared with universal scoring functions in actual drug research and development processes. A method that can effectively construct TSSFs will be of great value to drug design and discovery. In this work, we proposed a deep learning–based model named DeepScore to achieve this goal. DeepScore adopted the form of PMF scoring function to calculate protein–ligand binding affinity. However, different from PMF scoring function, in DeepScore, the score for each protein–ligand atom pair was calculated using a feedforward neural network. Our model significantly outperformed Glide Gscore on validation data set DUD-E. The average ROC-AUC on 102 targets was 0.98. We also combined Gscore and DeepScore together using a consensus method and put forward a consensus model named DeepScoreCS. The comparison results showed that DeepScore outperformed other machine learning–based TSSFs building methods. Furthermore, we presented a strategy to visualize the prediction of DeepScore. All of these results clearly demonstrated that DeepScore would be a useful model in constructing TSSFs and represented a novel way incorporating deep learning and drug design.

Importance of protein dynamics in the structure-based drug discovery of class A G protein-coupled receptors (GPCRs)

Demand for novel GPCR modulators is increasing as the association between the GPCR signaling pathway and numerous diseases such as cancers, psychological and metabolic disorders continues to be established. In silico structure-based drug design (SBDD) offers an outlet where researchers could exploit the accumulating structural information of GPCR to expedite the process of drug discovery. The coupling of structure-based approaches such as virtual screening and molecular docking with molecular dynamics and/or Monte Carlo simulation aids in reflecting the dynamics of proteins in nature into previously static docking studies, thus enhancing the accuracy of rationally designed ligands. This review will highlight recent computational strategies that incorporate protein flexibility into SBDD of GPCR-targeted ligands.

ImmtorLig_DB: repertoire of virtually screened small molecules against immune receptors to bolster host immunity

Host directed therapies to boost immunity against infection are gaining considerable impetus following the observation that use of antibiotics has become a continuous source for the emergence of drug resistant strains of pathogens. Receptors expressed by the cells of immune system play a cardinal role in initiating sequence of events necessary to ameliorate many morbid conditions. Although, ligands for the immune receptors are available; but their use is limited due to complex structure, synthesis and cost-effectiveness. Virtual screening (VS) is an integral part of chemoinformatics and computer-aided drug design (CADD) and aims to streamline the process of drug discovery. ImmtorLig_DB is a repertoire of 5000 novel small molecules, screened from ZINC database and ranked using structure based virtual screening (SBVS) against 25 immune receptors which play a pivotal role in defending and initiating the activation of immune system. Consequently, in the current study, small molecules were screened by docking on the essential domains present on the receptors expressed by cells of immune system. The screened molecules exhibited efficacious binding to immune receptors, and indicated a possibility of discovering novel small molecules. Other features of ImmtorLig_DB include information about availability, clustering analysis, and estimation of absorption, distribution, metabolism, and excretion (ADME) properties of the screened small molecules. Structural comparisons indicate that predicted small molecules may be considered novel. Further, this repertoire is available via a searchable graphical user interface (GUI) through http://bioinfo.imtech.res.in/bvs/immtor/ .

Ultra-large library docking for discovering new chemotypes

Despite intense interest in expanding chemical space, libraries containing hundreds-of-millions to billions of diverse molecules have remained inaccessible. Here we investigate structure-based docking of 170 million make-on-demand compounds from 130 well-characterized reactions. The resulting library is diverse, representing over 10.7 million scaffolds that are otherwise unavailable. For each compound in the library, docking against AmpC β-lactamase (AmpC) and the D4 dopamine receptor were simulated. From the top-ranking molecules, 44 and 549 compounds were synthesized and tested for interactions with AmpC and the D4 dopamine receptor, respectively. We found a phenolate inhibitor of AmpC, which revealed a group of inhibitors without known precedent. This molecule was optimized to 77 nM, which places it among the most potent non-covalent AmpC inhibitors known. Crystal structures of this and other AmpC inhibitors confirmed the docking predictions. Against the D4 dopamine receptor, hit rates fell almost monotonically with docking score, and a hit-rate versus score curve predicted that the library contained 453,000 ligands for the D4 dopamine receptor. Of 81 new chemotypes discovered, 30 showed submicromolar activity, including a 180-pM subtype-selective agonist of the D4 dopamine receptor.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

Approaches for Differentiation and Interconverting GPCR Agonists and Antagonists

Predicting the functional preferences of the ligands was always a highly demanding task, much harder that predicting whether a ligand can bind to the receptor. This is because of significant similarities of agonists, antagonists and inverse agonists which are binding usually in the same binding site of the receptor and only small structural changes can push receptor toward a particular activation state. For G protein-coupled receptors, due to a large progress in crystallization techniques and also in receptor thermal stabilization, it was possible to obtain a large number of high-quality structures of complexes of these receptors with agonists and non-agonists. Additionally, the long-time-scale molecular dynamics simulations revealed how the activation processes of GPCRs can take place. Using both theoretical and experimental knowledge it was possible to employ many clever and sophisticated methods which can help to differentiate agonists and non-agonists, so one can interconvert them in search of the optimal drug.

Breakthrough in GPCR Crystallography and Its Impact on Computer-Aided Drug Design

Recent crystallographic structures of G protein-coupled receptors (GPCRs) have greatly advanced our understanding of the recognition of their diverse agonist and antagonist ligands. We illustrate here how this applies to A_2A adenosine receptors (ARs) and to P2Y₁ and P2Y₁₂ receptors (P2YRs) for ADP. These X-ray structures have impacted the medicinal chemistry aimed at discovering new ligands for these two receptor families, including receptors that have not yet been crystallized but are closely related to the known structures. In this Chapter, we discuss recent structure-based drug design projects that led to the discovery of: (a) novel A₃AR agonists based on a highly rigidified (N)-methanocarba scaffold for the treatment of chronic neuropathic pain and other conditions, (b) fluorescent probes of the ARs and P2Y₁₄R, as chemical tools for structural probing of these GPCRs and for improving assay capabilities, and (c) new more drug-like antagonists of the inflammation-related P2Y₁₄R. We also describe the computationally enabled molecular recognition of positive (for A₃AR) and negative (P2Y₁R) allosteric modulators that in some cases are shown to be consistent with structure-activity relationship (SAR) data. Thus, computational modeling has become an essential tool for the design of purine receptor ligands.

3D-e-Chem: Structural Cheminformatics Workflows for Computer-Aided Drug Discovery

eScience technologies are needed to process the information available in many heterogeneous types of protein-ligand interaction data and to capture these data into models that enable the design of efficacious and safe medicines. Here we present scientific KNIME tools and workflows that enable the integration of chemical, pharmacological, and structural information for: i) structure-based bioactivity data mapping, ii) structure-based identification of scaffold replacement strategies for ligand design, iii) ligand-based target prediction, iv) protein sequence-based binding site identification and ligand repurposing, and v) structure-based pharmacophore comparison for ligand repurposing across protein families. The modular setup of the workflows and the use of well-established standards allows the re-use of these protocols and facilitates the design of customized computer-aided drug discovery workflows.

Exploring G Protein-Coupled Receptors (GPCRs) Ligand Space via Cheminformatics Approaches: Impact on Rational Drug Design

The primary goal of rational drug discovery is the identification of selective ligands which act on single or multiple drug targets to achieve the desired clinical outcome through the exploration of total chemical space. To identify such desired compounds, computational approaches are necessary in predicting their drug-like properties. G Protein-Coupled Receptors (GPCRs) represent one of the largest and most important integral membrane protein families. These receptors serve as increasingly attractive drug targets due to their relevance in the treatment of various diseases, such as inflammatory disorders, metabolic imbalances, cardiac disorders, cancer, monogenic disorders, etc. In the last decade, multitudes of three-dimensional (3D) structures were solved for diverse GPCRs, thus referring to this period as the "golden age for GPCR structural biology." Moreover, accumulation of data about the chemical properties of GPCR ligands has garnered much interest toward the exploration of GPCR chemical space. Due to the steady increase in the structural, ligand, and functional data of GPCRs, several cheminformatics approaches have been implemented in its drug discovery pipeline. In this review, we mainly focus on the cheminformatics-based paradigms in GPCR drug discovery. We provide a comprehensive view on the ligand- and structure-based cheminformatics approaches which are best illustrated via GPCR case studies. Furthermore, an appropriate combination of ligand-based knowledge with structure-based ones, i.e., integrated approach, which is emerging as a promising strategy for cheminformatics-based GPCR drug design is also discussed.

Benchmark Study Based on 2P2IDB to Gain Insights into the Discovery of Small-Molecule PPI Inhibitors

Protein-protein interactions (PPIs) have been regarded as novel and highly promising drug targets in drug discovery. Numerous new experimental techniques and computational approaches have been developed to assist the design of PPI modulators during the past two decades. However, identification and optimization of small-molecule inhibitors targeting PPIs is still a particularly challenging task due to the "undruggable" profiles of PPI interfaces. Nowadays, in silico screening, especially docking-based virtual screening, has emerged as an effective method to complement experimental high-throughput screening in identifying novel and potent small-molecule PPI inhibitors. Here, on the basis of the 2P2I_DB database, we explored the structural features of the known small-molecule PPI inhibitors and analyzed the characteristics of the PPI binding pockets. More importantly, we evaluated the sampling power and screening power of six popular docking programs for PPI targets. Our results indicate that the chlorinated conjugate group and amidelike linkage are two types of privileged fragments of PPI inhibitors; the average druggability of the binding sites of the PPI targets in 2P2IDB is slightly worse than that of traditional ones; both academic and commercial docking programs exhibit an acceptable accuracy on pose prediction for PPI inhibitors, but their screening powers for identifying PPI inhibitors are still not satisfactory. It is expected that our work can provide valuable guidance on the construction of PPI-focused library, the determination of druggable PPI binding pocket, and the selection of docking program for the screening of small-molecule PPI inhibitors.