Structure-Based Prediction of G-Protein-Coupled Receptor Ligand Function: A β-Adrenoceptor Case Study

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.

GPCR-IPL score: multilevel featurization of GPCR-ligand interaction patterns and prediction of ligand functions from selectivity to biased activation

G-protein-coupled receptors (GPCRs) mediate diverse cell signaling cascades after recognizing extracellular ligands. Despite the successful history of known GPCR drugs, a lack of mechanistic insight into GPCR challenges both the deorphanization of some GPCRs and optimization of the structure-activity relationship of their ligands. Notably, replacing a small substituent on a GPCR ligand can significantly alter extracellular GPCR-ligand interaction patterns and motion of transmembrane helices in turn to occur post-binding events of the ligand. In this study, we designed 3D multilevel features to describe the extracellular interaction patterns. Subsequently, these 3D features were utilized to predict the post-binding events that result from conformational dynamics from the extracellular to intracellular areas. To understand the adaptability of GPCR ligands, we collected the conformational information of flexible residues during binding and performed molecular featurization on a broad range of GPCR-ligand complexes. As a result, we developed GPCR-ligand interaction patterns, binding pockets, and ligand features as score (GPCR-IPL score) for predicting the functional selectivity of GPCR ligands (agonism versus antagonism), using the multilevel features of (1) zoomed-out 'residue level' (for flexible transmembrane helices of GPCRs), (2) zoomed-in 'pocket level' (for sophisticated mode of action) and (3) 'atom level' (for the conformational adaptability of GPCR ligands). GPCR-IPL score demonstrated reliable performance, achieving area under the receiver operating characteristic of 0.938 and area under the precision-recall curve of 0.907 (available in gpcr-ipl-score.onrender.com). Furthermore, we used the molecular features to predict the biased activation of downstream signaling (Gi/o, Gq/11, Gs and β-arrestin) as well as the functional selectivity. The resulting models are interpreted and applied to out-of-set validation with three scenarios including the identification of a new MRGPRX antagonist.

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.

Rational design, synthesis, and pharmacological evaluation of a cohort of novel beta-adrenergic receptors ligands enables an assessment of structure-activity relationships

Biomedical applications of molecules that are able to modulate β-adrenergic signaling have become increasingly attractive over the last decade, revealing that β-adrenergic receptors (β-ARs) are key targets for a plethora of therapeutic interventions, including cancer. Despite successes in β-AR drug discovery, identification of β-AR ligands that are useful as selective chemical tools in pharmacological studies of the three β-AR subtypes, or lead compounds for drug development is still a highly challenging task. This is mainly due to the intrinsic plasticity of β-ARs as G protein-coupled receptors in conjunction with the requirement for functional receptor subtype selectivity, tissue specificity and minimal off-target effects. With the aim to provide insight into structure-activity relationships for the three β-AR subtypes, we have synthesized and obtained the pharmacological profile of a series of structurally diverse compounds (named MC) that were designed based on the aryloxy-propanolamine scaffold of SR59230A. Comparative analysis of their predicted binding mode within the active and inactive states of the receptors in combination with their pharmacological profile revealed key structural elements that control their activity as agonists or antagonists, in addition to clues about substituents that mediate selectivity for one receptor subtype over the others. We anticipate that these results will facilitate selective β-AR drug development efforts.

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.

Accelerating GPCR Drug Discovery With Conformation-Stabilizing VHHs

The human genome encodes 850 G protein-coupled receptors (GPCRs), half of which are considered potential drug targets. GPCRs transduce extracellular stimuli into a plethora of vital physiological processes. Consequently, GPCRs are an attractive drug target class. This is underlined by the fact that approximately 40% of marketed drugs modulate GPCRs. Intriguingly 60% of non-olfactory GPCRs have no drugs or candidates in clinical development, highlighting the continued potential of GPCRs as drug targets. The discovery of small molecules targeting these GPCRs by conventional high throughput screening (HTS) campaigns is challenging. Although the definition of success varies per company, the success rate of HTS for GPCRs is low compared to other target families (Fujioka and Omori, 2012; Dragovich et al., 2022). Beyond this, GPCR structure determination can be difficult, which often precludes the application of structure-based drug design approaches to arising HTS hits. GPCR structural studies entail the resource-demanding purification of native receptors, which can be challenging as they are inherently unstable when extracted from the lipid matrix. Moreover, GPCRs are flexible molecules that adopt distinct conformations, some of which need to be stabilized if they are to be structurally resolved. The complexity of targeting distinct therapeutically relevant GPCR conformations during the early discovery stages contributes to the high attrition rates for GPCR drug discovery programs. Multiple strategies have been explored in an attempt to stabilize GPCRs in distinct conformations to better understand their pharmacology. This review will focus on the use of camelid-derived immunoglobulin single variable domains (VHHs) that stabilize disease-relevant pharmacological states (termed ConfoBodies by the authors) of GPCRs, as well as GPCR:signal transducer complexes, to accelerate drug discovery. These VHHs are powerful tools for supporting in vitro screening, deconvolution of complex GPCR pharmacology, and structural biology purposes. In order to demonstrate the potential impact of ConfoBodies on translational research, examples are presented of their role in active state screening campaigns and structure-informed rational design to identify de novo chemical space and, subsequently, how such matter can be elaborated into more potent and selective drug candidates with intended pharmacology.

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.

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.

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

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

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.

α1-adrenoceptor activity of β-adrenoceptor ligands - An expected drug property with limited clinical relevance

Many β-adrenoceptor agonists and antagonists including several clinically used drugs have been reported to also exhibit binding to α₁-adrenoceptors. Such promiscuity within the adrenoceptor family appears to occur more often than off-target effects of drugs in general. It should not be considered surprising based on the amino acid homology among the nine adrenoceptor subtypes including the counter-ions for binding the endogenous catecholamines. When β-adrenoceptor ligands also bind to α₁-adrenoceptors, they almost always act as antagonists, regardless of being agonists or antagonists at the β-adrenoceptor. The α₁-adrenoceptor affinity of β-adrenoceptor ligands in most cases is at least one, and often more log units lower than at their cognate receptor. Consistent evidence from multiple investigators indicates that β-adrenoceptor ligands relatively have the highest affinity for α_1A- and lowest for α_1B-adrenoceptors. While promiscuity among adrenoceptor subtypes causes misleading interpretation of experimental in vitro data, it is proposed based on the law of mass action that α₁-adrenoceptor binding of β-adrenoceptor ligands rarely contributes to the clinical profile of such drugs, particularly if they are agonists at the β-adrenoceptor.

A Systematic Review of Inverse Agonism at Adrenoceptor Subtypes

As many, if not most, ligands at G protein-coupled receptor antagonists are inverse agonists, we systematically reviewed inverse agonism at the nine adrenoceptor subtypes. Except for β₃-adrenoceptors, inverse agonism has been reported for each of the adrenoceptor subtypes, most often for β₂-adrenoceptors, including endogenously expressed receptors in human tissues. As with other receptors, the detection and degree of inverse agonism depend on the cells and tissues under investigation, i.e., they are greatest when the model has a high intrinsic tone/constitutive activity for the response being studied. Accordingly, they may differ between parts of a tissue, for instance, atria vs. ventricles of the heart, and within a cell type, between cellular responses. The basal tone of endogenously expressed receptors is often low, leading to less consistent detection and a lesser extent of observed inverse agonism. Extent inverse agonism depends on specific molecular properties of a compound, but inverse agonism appears to be more common in certain chemical classes. While inverse agonism is a fascinating facet in attempts to mechanistically understand observed drug effects, we are skeptical whether an a priori definition of the extent of inverse agonism in the target product profile of a developmental candidate is a meaningful option in drug discovery and development.

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.

Chemokine Receptor Crystal Structures: What Can Be Learned from Them?

Chemokine receptors belong to the class A of G protein-coupled receptors (GPCRs) and are implicated in a wide variety of physiologic functions, mostly related to the homeostasis of the immune system. Chemokine receptors are also involved in multiple pathologic processes, including immune and autoimmune diseases, as well as cancer. Hence, several members of this GPCR subfamily are considered to be very relevant therapeutic targets. Since drug discovery efforts can be significantly reinforced by the availability of crystal structures, substantial efforts in the area of chemokine receptor structural biology could dramatically increase the outcome of drug discovery campaigns. This short review summarizes the available data on chemokine receptor crystal structures, discusses the numerous applications from chemokine receptor structures that can enhance the daily work of molecular pharmacologists, and describes the challenges and pitfalls to consider when relying on crystal structures for further research applications. SIGNIFICANCE STATEMENT: This short review summarizes the available data on chemokine receptor crystal structures, discusses the numerous applications from chemokine receptor structures that can enhance the daily work of molecular pharmacologists, and describes the challenges and pitfalls to consider when relying on crystal structures for further research applications.

Ligand-binding affinity of alternative conformers of human β2 -adrenergic receptor in the presence of intracellular loop 3 (ICL3) and their potential use in virtual screening studies

This study investigates the structural distinctiveness of orthosteric ligand-binding sites of several human β₂ adrenergic receptor (β₂ -AR) conformations that have been obtained from a set of independent molecular dynamics (MD) simulations in the presence of intracellular loop 3 (ICL3). A docking protocol was established in order to classify each receptor conformation via its binding affinity to selected ligands with known efficacy. This work's main goal was to reveal many subtle features of the ligand-binding site, presenting alternative conformations, which might be considered as either active- or inactive-like but mostly specific for that ligand. Agonists, inverse agonists, and antagonists were docked to each MD conformer with distinct binding pockets, using different docking tools and scoring functions. Mostly favored receptor conformation persistently observed in all docking/scoring evaluations was classified as active or inactive based on the type of ligand's biological effect. Classified MD conformers were further tested for their ability to discriminate agonists from inverse agonists/antagonists, and several conformers were proposed as important targets to be used in virtual screening experiments that were often limited to a single X-ray structure.

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.

A benchmark driven guide to binding site comparison: An exhaustive evaluation using tailor-made data sets (ProSPECCTs)

The automated comparison of protein-ligand binding sites provides useful insights into yet unexplored site similarities. Various stages of computational and chemical biology research can benefit from this knowledge. The search for putative off-targets and the establishment of polypharmacological effects by comparing binding sites led to promising results for numerous projects. Although many cavity comparison methods are available, a comprehensive analysis to guide the choice of a tool for a specific application is wanting. Moreover, the broad variety of binding site modeling approaches, comparison algorithms, and scoring metrics impedes this choice. Herein, we aim to elucidate strengths and weaknesses of binding site comparison methodologies. A detailed benchmark study is the only possibility to rationalize the selection of appropriate tools for different scenarios. Specific evaluation data sets were developed to shed light on multiple aspects of binding site comparison. An assembly of all applied benchmark sets (ProSPECCTs–Protein Site Pairs for the Evaluation of Cavity Comparison Tools) is made available for the evaluation and optimization of further and still emerging methods. The results indicate the importance of such analyses to facilitate the choice of a methodology that complies with the requirements of a specific scientific challenge.

Binding site similarities are useful in the context of promiscuity prediction, drug repurposing, the analysis of protein-ligand and protein-protein complexes, function prediction, and further fields of general interest in chemical biology and biochemistry. Many years of research have led to the development of a multitude of methods for binding site analysis and comparison. On the one hand, their availability supports research. On the other hand, the huge number of methods hampers the efficient selection of a specific tool. Our research is dedicated to the analysis of different cavity comparison tools. We use several binding site data sets to establish guidelines which can be applied to ensure a successful application of comparison methods by circumventing potential pitfalls.

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.

A Structural Framework for GPCR Chemogenomics: What's In a Residue Number?

The recent surge of crystal structures of G protein-coupled receptors (GPCRs), as well as comprehensive collections of sequence, structural, ligand bioactivity, and mutation data, has enabled the development of integrated chemogenomics workflows for this important target family. This chapter will focus on cross-family and cross-class studies of GPCRs that have pinpointed the need for, and the implementation of, a generic numbering scheme for referring to specific structural elements of GPCRs. Sequence- and structure-based numbering schemes for different receptor classes will be introduced and the remaining caveats will be discussed. The use of these numbering schemes has facilitated many chemogenomics studies such as consensus binding site definition, binding site comparison, ligand repurposing (e.g. for orphan receptors), sequence-based pharmacophore generation for homology modeling or virtual screening, and class-wide chemogenomics studies of GPCRs.

Structural Features and Ligand Selectivity for 10 Intermediates in the Activation Process of β2-Adrenergic Receptor

It has already been suggested by researchers that there should be multiple intermediate states in the activation process for G-protein-coupled receptors (GPCRs). However, the intermediate states are very short-lived and hardly captured by the experiments, leading to very limited understanding of their structural features and drug efficacies. In this work, a novel joint strategy of targeted molecular dynamics simulation, conventional molecular dynamics simulation, and virtual screening is developed to address the problems. The results from 10 intermediate conformations obtained from the work reveal that the ligand pocket is very unstable and fluctuates between the inactive state and the active one in the case of ligand-free, in particular for ECL2 as a gate-keeper of the ligand-binding. The ligand-binding site could be stable in the active state with a small volume and a completely closed ECL2, only when the G-protein-binding region is fully activated. In addition, the activations of the ligand-binding pocket and G-protein-binding site are relatively independent and exhibit a loose allosteric coupling, which contributes to the existence of multiple intermediate conformations. Interestingly, the screening performance of the agonists does not increase on increasing the overall activity of the intermediate state, but is dependent on the activated extent of the ligand pocket. The receptor is prone to bind the agonist when closing ECL2 and reducing the ligand-binding pocket volume, whereas it is more favorable for binding the antagonist when opening ECL2 and increasing the pocket volume. These observations added to previous studies could help us better understand the activation mechanism of GPCRs and provide valuable information for drug design.

Methods for Virtual Screening of GPCR Targets: Approaches and Challenges

Virtual screening (VS) has become an integral part of the drug discovery process and is a valuable tool for finding novel chemical starting points for GPCR targets. Ligand-based VS makes use of biochemical data for known, active compounds and has been applied successfully to many diverse GPCRs. Recent progress in GPCR X-ray crystallography has made it possible to incorporate detailed structural information into the VS process. This chapter outlines the latest VS techniques along with examples that highlight successful applications of these methods. Best practices for increasing the likelihood of VS success, as well as ongoing challenges, are also discussed.

Dynamical Binding Modes Determine Agonistic and Antagonistic Ligand Effects in the Prostate-Specific G-Protein Coupled Receptor (PSGR)

We analysed the ligand-based activation mechanism of the prostate-specific G-protein coupled receptor (PSGR), which is an olfactory receptor that mediates cellular growth in prostate cancer cells. Furthermore, it is an olfactory receptor with a known chemically near identic antagonist/agonist pair, α- and β-ionone. Using a combined theoretical and experimental approach, we propose that this receptor is activated by a ligand-induced rearrangement of a protein-internal hydrogen bond network. Surprisingly, this rearrangement is not induced by interaction of the ligand with the network, but by dynamic van der Waals contacts of the ligand with the involved amino acid side chains, altering their conformations and intraprotein connectivity. Ligand recognition in this GPCR is therefore highly stereo selective, but seemingly lacks any ligand recognition via polar contacts. A putative olfactory receptor-based drug design scheme will have to take this unique mode of protein/ligand action into account.

Improving virtual screening of G protein-coupled receptors via ligand-directed modeling

G protein-coupled receptors (GPCRs) play crucial roles in cell physiology and pathophysiology. There is increasing interest in using structural information for virtual screening (VS) of libraries and for structure-based drug design to identify novel agonist or antagonist leads. However, the sparse availability of experimentally determined GPCR/ligand complex structures with diverse ligands impedes the application of structure-based drug design (SBDD) programs directed to identifying new molecules with a select pharmacology. In this study, we apply ligand-directed modeling (LDM) to available GPCR X-ray structures to improve VS performance and selectivity towards molecules of specific pharmacological profile. The described method refines a GPCR binding pocket conformation using a single known ligand for that GPCR. The LDM method is a computationally efficient, iterative workflow consisting of protein sampling and ligand docking. We developed an extensive benchmark comparing LDM-refined binding pockets to GPCR X-ray crystal structures across seven different GPCRs bound to a range of ligands of different chemotypes and pharmacological profiles. LDM-refined models showed improvement in VS performance over origin X-ray crystal structures in 21 out of 24 cases. In all cases, the LDM-refined models had superior performance in enriching for the chemotype of the refinement ligand. This likely contributes to the LDM success in all cases of inhibitor-bound to agonist-bound binding pocket refinement, a key task for GPCR SBDD programs. Indeed, agonist ligands are required for a plethora of GPCRs for therapeutic intervention, however GPCR X-ray structures are mostly restricted to their inactive inhibitor-bound state.

G protein-coupled receptors (GPCRs) are a major target for drug discovery. These receptors are highly dynamic membrane proteins, and have had limited tractability using with biophysical screens that are widely adopted for globular protein targets. Thus, structure-based virtual screening (SBVS) holds great promise as a complement to physical screening for rational design of novel drugs. Indeed, the increasing number of atomic-detail GPCR X-ray crystal structures has coincided with an increase in prospective SBVS studies that have identified novel compounds. However, experimentally solved GPCR structures do not meet the full demand for SBVS, as the GPCR structural landscape is incomplete, lacking both in coverage of available GPCRs, and diversity in both receptor conformations and the chemistry of co-crystalised ligands. Here we present a novel computational GPCR binding pocket refinement method that can generate predictive GPCR/ligand complexes with improved SBVS performance. This ligand-directed modeling workflow uses parallel processing and efficient algorithms to search the GPCR/ligand conformational space faster and more efficiently than the widely used protein refinement method molecular dynamics. In this study, the resulting models are evaluated both structurally, and in retrospective SBVS. We demonstrate improved performance of refined models over their starting structures in the majority of our test cases.

Structure-Guided Screening for Functionally Selective D2 Dopamine Receptor Ligands from a Virtual Chemical Library

Functionally selective ligands stabilize conformations of G protein-coupled receptors (GPCRs) that induce a preference for signaling via a subset of the intracellular pathways activated by the endogenous agonists. The possibility to fine-tune the functional activity of a receptor provides opportunities to develop drugs that selectively signal via pathways associated with a therapeutic effect and avoid those causing side effects. Animal studies have indicated that ligands displaying functional selectivity at the D₂ dopamine receptor (D₂R) could be safer and more efficacious drugs against neuropsychiatric diseases. In this work, computational design of functionally selective D₂R ligands was explored using structure-based virtual screening. Molecular docking of known functionally selective ligands to a D₂R homology model indicated that such compounds were anchored by interactions with the orthosteric site and extended into a common secondary pocket. A tailored virtual library with close to 13 000 compounds bearing 2,3-dichlorophenylpiperazine, a privileged orthosteric scaffold, connected to diverse chemical moieties via a linker was docked to the D₂R model. Eighteen top-ranked compounds that occupied both the orthosteric and allosteric site were synthesized, leading to the discovery of 16 partial agonists. A majority of the ligands had comparable maximum effects in the G protein and β-arrestin recruitment assays, but a subset displayed preference for a single pathway. In particular, compound 4 stimulated β-arrestin recruitment (EC₅₀ = 320 nM, E_max = 16%) but had no detectable G protein signaling. The use of structure-based screening and virtual libraries to discover GPCR ligands with tailored functional properties will be discussed.

Structural features embedded in G protein-coupled receptor co-crystal structures are key to their success in virtual screening

Structure based drug discovery on GPCRs harness atomic detail X-ray binding pockets and large libraries of potential drug lead candidates in virtual screening (VS) to identify novel lead candidates. Relatively small conformational differences between such binding pockets can be critical to the success of VS. Retrospective VS on GPCR/ligand co-crystal structures revealed stark differences in the ability of different structures to identify known ligands, despite being co-crystallized with the same ligand. When using the OpenEye toolkit and the ICM modeling package, we identify criteria associated with the predictive power of binding pockets in VS that consists of a combination of ligand/receptor interaction pattern and predicted ligand/receptor interaction strength. These findings can guide the selection and refinement of GPCR binding pockets for use in SBDD programs and may also provide a potential framework for evaluating the ability of computational GPCR binding pocket refinement tools in improving the predictive power of binding pockets.

Structural Analysis of Chemokine Receptor-Ligand Interactions

This review focuses on the construction and application of structural chemokine receptor models for the elucidation of molecular determinants of chemokine receptor modulation and the structure-based discovery and design of chemokine receptor ligands. A comparative analysis of ligand binding pockets in chemokine receptors is presented, including a detailed description of the CXCR4, CCR2, CCR5, CCR9, and US28 X-ray structures, and their implication for modeling molecular interactions of chemokine receptors with small-molecule ligands, peptide ligands, and large antibodies and chemokines. These studies demonstrate how the integration of new structural information on chemokine receptors with extensive structure–activity relationship and site-directed mutagenesis data facilitates the prediction of the structure of chemokine receptor–ligand complexes that have not been crystallized. Finally, a review of structure-based ligand discovery and design studies based on chemokine receptor crystal structures and homology models illustrates the possibilities and challenges to find novel ligands for chemokine receptors.

Interacting with GPCRs: Using Interaction Fingerprints for Virtual Screening

The expanding number of crystal structures of G protein-coupled receptors (GPCRs) has increased the knowledge on receptor function and their ability to recognize ligands. Although structure-based virtual screening has been quite successful on GPCRs, scores obtained by docking are typically not indicative for ligand affinity. Methods capturing interactions between protein and ligand in a more explicit manner, such as interaction fingerprints (IFPs), have been applied as an addition or alternative to docking. Originally IFPs captured the interactions of amino acid residues with ligands with specific definitions for the various interaction types. More complex IFPs now capture atom-atom interactions, such as in SYBYL, or fragment-fragment co-occurrences such as in SPLIF. Overall, most of the IFPs have been studied in comparison with docking in retrospective studies. For GPCRs it remains unclear which IFP should be used, if at all, and in what manner. Thus, the performance between five different IFPs was compared on five different representative GPCRs, including several extensions of the original implementations,. Results show that the more detailed IFPs, SYBYL and SPLIF, perform better than the other IFPs (Deng, Credo, and Elements). SPLIF was further tuned based on the number of poses, fingerprint similarity coefficient, and using an ensemble of structures. Enrichments were obtained that were significantly higher than initial enrichments and those obtained by 2D-similarity. With the increase in available crystal structures for GPCRs, and given that IFPs such as SPLIF enhance enrichment in virtual screens, it is anticipated that IFPs will be used in conjunction with docking, especially for GPCRs with a large binding pocket.

Molecular interaction fingerprint approaches for GPCR drug discovery

Protein-ligand interaction fingerprints (IFPs) are binary 1D representations of the 3D structure of protein-ligand complexes encoding the presence or absence of specific interactions between the binding pocket amino acids and the ligand. Various implementations of IFPs have been developed and successfully applied for post-processing molecular docking results for G Protein-Coupled Receptor (GPCR) ligand binding mode prediction and virtual ligand screening. Novel interaction fingerprint methods enable structural chemogenomics and polypharmacology predictions by complementing the increasing amount of GPCR structural data. Machine learning methods are increasingly used to derive relationships between bioactivity data and fingerprint descriptors of chemical and structural information of binding sites, ligands, and protein-ligand interactions. Factors that influence the application of IFPs include structure preparation, binding site definition, fingerprint similarity assessment, and data processing and these factors pose challenges as well possibilities to optimize interaction fingerprint methods for GPCR drug discovery.

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

The ability of scoring functions to correctly select and rank docking poses of small molecules in protein binding sites is highly target dependent, which presents a challenge for structure-based drug discovery. Here we describe a virtual screening method that combines an energy-based docking scoring function with a molecular interaction fingerprint (IFP) to identify new ligands based on G protein-coupled receptor (GPCR) crystal structures. The consensus scoring method is prospectively evaluated by: 1) the discovery of chemically novel, fragment-like, high affinity histamine H₁ receptor (H₁R) antagonists/inverse agonists, 2) the selective structure-based identification of ß₂-adrenoceptor (ß₂R) agonists, and 3) the experimental validation and comparison of the combined and individual scoring approaches. Systematic retrospective virtual screening simulations allowed the definition of scoring cut-offs for the identification of H₁R and ß₂R ligands and the selection of an optimal ß-adrenoceptor crystal structure for the discrimination between ß₂R agonists and antagonists. The consensus approach resulted in the experimental validation of 53% of the ß₂R and 73% of the H₁R virtual screening hits with up to nanomolar affinities and potencies. The selective identification of ß₂R agonists shows the possibilities of structure-based prediction of GPCR ligand function by integrating protein-ligand binding mode information.

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

Modern rational drug design not only deals with the search for ligands binding to interesting and promising validated targets but also aims to identify the function and ligands of yet uncharacterized proteins having impact on different diseases. Additionally, it contributes to the design of inhibitors with distinct selectivity patterns and the prediction of possible off-target effects. The identification of similarities between binding sites of various proteins is a useful approach to cope with those challenges. The main scope of this perspective is to describe applications of different protein binding site comparison approaches to outline their applicability and impact on molecular design. The article deals with various substantial application domains and provides some outstanding examples to show how various binding site comparison methods can be applied to promote in silico drug design workflows. In addition, we will also briefly introduce the fundamental principles of different protein binding site comparison methods.