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.

Impact of GPCR Structures on Drug Discovery.

Structures of 70 unique G protein-coupled receptors (GPCRs) have been determined, with over 370 structures in total bound to different ligands and the receptors in various conformational states. Structure-based drug design has been applied to an increasing number of GPCR targets over the past decade and now a few of these drug candidates have entered clinical trials. Given the length of time required for a drug to reach the market, there are no documented examples of licensed drugs being developed with the aid of a structure, but this is likely to change as current efforts come to fruition.

Conformational dynamics underlying atypical chemokine receptor 3 activation.

Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon ß-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of ß-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins.

Pharmacological Characterization and Radiolabeling of VUF15485, a High-Affinity Small-Molecule Agonist for the Atypical Chemokine Receptor ACKR3.

Atypical chemokine receptor 3 (ACKR3), formerly referred to as CXCR7, is considered to be an interesting drug target. In this study, we report on the synthesis, pharmacological characterization and radiolabeling of VUF15485, a new ACKR3 small-molecule agonist, that will serve as an important new tool to study this ß-arrestin-biased chemokine receptor. VUF15485 binds with nanomolar affinity (pIC50 = 8.3) to human ACKR3, as measured in [125I]CXCL12 competition binding experiments. Moreover, in a bioluminescence resonance energy transfer-based ß-arrestin2 recruitment assay VUF15485 acts as a potent ACKR3 agonist (pEC50 = 7.6) and shows a similar extent of receptor activation compared with CXCL12 when using a newly developed, fluorescence resonance energy transfer-based ACKR3 conformational sensor. Moreover, the ACKR3 agonist VUF15485, tested against a (atypical) chemokine receptor panel (agonist and antagonist mode), proves to be selective for ACKR3. VUF15485 labeled with tritium at one of its methoxy groups ([3H]VUF15485), binds ACKR3 saturably and with high affinity (K d = 8.2 nM). Additionally, [3H]VUF15485 shows rapid binding kinetics and consequently a short residence time (<2 minutes) for binding to ACKR3. The selectivity of [3H]VUF15485 for ACKR3, was confirmed by binding studies, whereupon CXCR3, CXCR4, and ACKR3 small-molecule ligands were competed for binding against the radiolabeled agonist. Interestingly, the chemokine ligands CXCL11 and CXCL12 are not able to displace the binding of [3H]VUF15485 to ACKR3. The radiolabeled VUF15485 was subsequently used to evaluate its binding pocket. Site-directed mutagenesis and docking studies using a recently solved cryo-EM structure propose that VUF15485 binds in the major and the minor binding pocket of ACKR3. SIGNIFICANCE STATEMENT: The atypical chemokine receptor atypical chemokine receptor 3 (ACKR3) is considered an interesting drug target in relation to cancer and multiple sclerosis. The study reports on new chemical biology tools for ACKR3, i.e., a new agonist that can also be radiolabeled and a new ACKR3 conformational sensor, that both can be used to directly study the interaction of ACKR3 ligands with the G protein-coupled receptor.

Predicting the target landscape of kinase inhibitors using 3D convolutional neural networks.

Many therapies in clinical trials are based on single drug-single target relationships. To further extend this concept to multi-target approaches using multi-targeted drugs, we developed a machine learning pipeline to unravel the target landscape of kinase inhibitors. This pipeline, which we call 3D-KINEssence, uses a new type of protein fingerprints (3D FP) based on the structure of kinases generated through a 3D convolutional neural network (3D-CNN). These 3D-CNN kinase fingerprints were matched to molecular Morgan fingerprints to predict the targets of each respective kinase inhibitor based on available bioactivity data. The performance of the pipeline was evaluated on two test sets: a sparse drug-target set where each drug is matched in most cases to a single target and also on a densely-covered drug-target set where each drug is matched to most if not all targets. This latter set is more challenging to train, given its non-exclusive character. Our model's root-mean-square error (RMSE) based on the two datasets was 0.68 and 0.8, respectively. These results indicate that 3D FP can predict the target landscape of kinase inhibitors at around 0.8 log units of bioactivity. Our strategy can be utilized in proteochemometric or chemogenomic workflows by consolidating the target landscape of kinase inhibitors.

Comparative Study of Allosteric GPCR Binding Sites and Their Ligandability Potential.

The steadily growing number of experimental G-protein-coupled receptor (GPCR) structures has revealed diverse locations of allosteric modulation, and yet few drugs target them. This gap highlights the need for a deeper understanding of allosteric modulation in GPCR drug discovery. The current work introduces a systematic annotation scheme to structurally classify GPCR binding sites based on receptor class, transmembrane helix contacts, and, for membrane-facing sites, membrane sublocation. This GPCR specific annotation scheme was applied to 107 GPCR structures bound by small molecules contributing to 24 distinct allosteric binding sites for comparative evaluation of three binding site detection methods (BioGPS, SiteMap, and FTMap). BioGPS identified the most in 22 of 24 sites. In addition, our property analysis showed that extrahelical allosteric ligands and binding sites represent a distinct chemical space characterized by shallow pockets with low volume, and the corresponding allosteric ligands showed an enrichment of halogens. Furthermore, we demonstrated that combining receptor and ligand similarity can be a viable method for ligandability assessment. One challenge regarding site prediction is the ligand shaping effect on the observed binding site, especially for extrahelical sites where the ligand-induced effect was most pronounced. To our knowledge, this is the first study presenting a binding site annotation scheme standardized for GPCRs, and it allows a comparison of allosteric binding sites across different receptors in an objective way. The insight from this study provides a framework for future GPCR binding site studies and highlights the potential of targeting allosteric sites for drug development.

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.

Human GLP-1 receptor transmembrane domain structure in complex with allosteric modulators.

The glucagon-like peptide-1 receptor (GLP-1R) and the glucagon receptor (GCGR) are members of the secretin-like class B family of G-protein-coupled receptors (GPCRs) and have opposing physiological roles in insulin release and glucose homeostasis. The treatment of type 2 diabetes requires positive modulation of GLP-1R to inhibit glucagon secretion and stimulate insulin secretion in a glucose-dependent manner. Here we report crystal structures of the human GLP-1R transmembrane domain in complex with two different negative allosteric modulators, PF-06372222 and NNC0640, at 2.7 and 3.0 Å resolution, respectively. The structures reveal a common binding pocket for negative allosteric modulators, present in both GLP-1R and GCGR and located outside helices V-VII near the intracellular half of the receptor. The receptor is in an inactive conformation with compounds that restrict movement of the intracellular tip of helix VI, a movement that is generally associated with activation mechanisms in class A GPCRs. Molecular modelling and mutagenesis studies indicate that agonist positive allosteric modulators target the same general region, but in a distinct sub-pocket at the interface between helices V and VI, which may facilitate the formation of an intracellular binding site that enhances G-protein coupling.

Structure of the full-length glucagon class B G-protein-coupled receptor.

The human glucagon receptor, GCGR, belongs to the class B G-protein-coupled receptor family and plays a key role in glucose homeostasis and the pathophysiology of type 2 diabetes. Here we report the 3.0 Å crystal structure of full-length GCGR containing both the extracellular domain and transmembrane domain in an inactive conformation. The two domains are connected by a 12-residue segment termed the stalk, which adopts a ß-strand conformation, instead of forming an α-helix as observed in the previously solved structure of the GCGR transmembrane domain. The first extracellular loop exhibits a ß-hairpin conformation and interacts with the stalk to form a compact ß-sheet structure. Hydrogen-deuterium exchange, disulfide crosslinking and molecular dynamics studies suggest that the stalk and the first extracellular loop have critical roles in modulating peptide ligand binding and receptor activation. These insights into the full-length GCGR structure deepen our understanding of the signalling mechanisms of class B G-protein-coupled receptors.

KLIFS: an overhaul after the first 5 years of supporting kinase research.

Kinases are a prime target of drug development efforts with >60 drug approvals in the past two decades. Due to the research into this protein family, a wealth of data has been accumulated that keeps on growing. KLIFS-Kinase-Ligand Interaction Fingerprints and Structures-is a structural database focusing on how kinase inhibitors interact with their targets. The aim of KLIFS is to support (structure-based) kinase research through the systematic collection, annotation, and processing of kinase structures. Now, 5 years after releasing the initial KLIFS website, the database has undergone a complete overhaul with a new website, new logo, and new functionalities. In this article, we start by looking back at how KLIFS has been used by the research community, followed by a description of the renewed KLIFS, and conclude with showcasing the functionalities of KLIFS. Major changes include the integration of approved drugs and inhibitors in clinical trials, extension of the coverage to atypical kinases, and a RESTful API for programmatic access. KLIFS is available at the new domain https://klifs.net.

Advanced fluorescence microscopy reveals disruption of dynamic CXCR4 dimerization by subpocket-specific inverse agonists.

Although class A G protein-coupled receptors (GPCRs) can function as monomers, many of them form dimers and oligomers, but the mechanisms and functional relevance of such oligomerization is ill understood. Here, we investigate this problem for the CXC chemokine receptor 4 (CXCR4), a GPCR that regulates immune and hematopoietic cell trafficking, and a major drug target in cancer therapy. We combine single-molecule microscopy and fluorescence fluctuation spectroscopy to investigate CXCR4 membrane organization in living cells at densities ranging from a few molecules to hundreds of molecules per square micrometer of the plasma membrane. We observe that CXCR4 forms dynamic, transient homodimers, and that the monomer-dimer equilibrium is governed by receptor density. CXCR4 inverse agonists that bind to the receptor minor pocket inhibit CXCR4 constitutive activity and abolish receptor dimerization. A mutation in the minor binding pocket reduced the dimer-disrupting ability of these ligands. In addition, mutating critical residues in the sixth transmembrane helix of CXCR4 markedly diminished both basal activity and dimerization, supporting the notion that CXCR4 basal activity is required for dimer formation. Together, these results link CXCR4 dimerization to its density and to its activity. They further suggest that inverse agonists binding to the minor pocket suppress both dimerization and constitutive activity and may represent a specific strategy to target CXCR4.

SAR exploration of the non-imidazole histamine H3 receptor ligand ZEL-H16 reveals potent inverse agonism.

Histamine H3 receptor (H3 R) agonists without an imidazole moiety remain very scarce. Of these, ZEL-H16 (1) has been reported previously as a high-affinity non-imidazole H3 R (partial) agonist. Our structure-activity relationship analysis using derivatives of 1 identified both basic moieties as key interaction motifs and the distance of these from the central core as a determinant for H3 R affinity. However, in spite of the reported H3 R (partial) agonism, in our hands, 1 acts as an inverse agonist for Gαi signaling in a CRE-luciferase reporter gene assay and using an H3 R conformational sensor. Inverse agonism was also observed for all of the synthesized derivatives of 1. Docking studies and molecular dynamics simulations suggest ionic interactions/hydrogen bonds to H3 R residues D1143.32 and E2065.46 as essential interaction points.

Classification Model for the Second Extracellular Loop of Class A GPCRs.

The extracellular loop 2 (ECL2) is the longest and the most diverse loop among class A G protein-coupled receptors (GPCRs). It connects the transmembrane (TM) helices 4 and 5 and contains a highly conserved cysteine through which it is bridged with TM3. In this paper, experimental ECL2 structures were analyzed based on their sequences, shapes, and intramolecular contacts. To take into account the flexibility, we incorporated into our analyses information from the molecular dynamics trajectories available on the GPCRmd website. Despite the high sequence variability, shapes of the analyzed structures, defined by the backbone volume overlaps, can be clustered into seven main groups. Conformational differences within the clusters can be then identified by intramolecular interactions with other GPCR structural domains. Overall, our work provides a reorganization of the structural information of the ECL2 of class A GPCR subfamilies, highlighting differences and similarities on sequence and conformation levels.

Extending the Structural View of Class B GPCRs.

The secretin-like class B family of G protein-coupled receptors (GPCRs) are key players in hormonal homeostasis. Recent structures of various receptors in complex with a variety of orthosteric and allosteric ligands provide fundamental new insights into the function and mechanism of class B GPCRs, including: (i) ligand-induced changes in the relative orientation of the extracellular and transmembrane receptor domains; (ii) intramolecular interaction networks that stabilize conformational changes to accommodate intracellular G protein binding; and (iii) allosteric modulation of receptor activation. This review provides a comprehensive analysis of the structural, biochemical, and pharmacological data on class B GPCRs for understanding ligand-receptor interaction and modulation mechanisms and assessing the potential implications for drug discovery for the secretin-like GPCR family.

Controlling the selectivity of aminergic GPCR ligands from the extracellular vestibule.

In addition to the orthosteric binding pocket (OBP) of GPCRs, recent structural studies have revealed that there are several allosteric sites available for pharmacological intervention. The secondary binding pocket (SBP) of aminergic GPCRs is located in the extracellular vestibule of these receptors, and it has been suggested to be a potential selectivity pocket for bitopic ligands. Here, we applied a virtual screening protocol based on fragment docking to the SBP of the orthosteric ligand-receptor complex. This strategy was employed for a number of aminergic receptors. First, we designed dopamine D3 preferring bitopic compounds from a D2 selective orthosteric ligand. Next, we designed 5-HT2B selective bitopic compounds starting from the 5-HT1B preferring ergoline core of LSD. Comparing the serotonergic profiles of the new derivatives to that of LSD, we found that these derivatives became significantly biased towards the desired 5-HT2B receptor target. Finally, addressing the known limitations of H1 antihistamines, our protocol was successfully used to eliminate the well-known side effects related to the muscarinic M1 activity of amitriptyline while preserving H1 potency in some of the designed bitopic compounds. These applications highlight the usefulness of our new virtual screening protocol and offer a powerful strategy towards bitopic GPCR ligands with designed receptor profiles.

X-Ray Crystallography and Free Energy Calculations Reveal the Binding Mechanism of A2A Adenosine Receptor Antagonists.

We present a robust protocol based on iterations of free energy perturbation (FEP) calculations, chemical synthesis, biophysical mapping and X-ray crystallography to reveal the binding mode of an antagonist series to the A2A adenosine receptor (AR). Eight A2A AR binding site mutations from biophysical mapping experiments were initially analyzed with sidechain FEP simulations, performed on alternate binding modes. The results distinctively supported one binding mode, which was subsequently used to design new chromone derivatives. Their affinities for the A2A AR were experimentally determined and investigated through a cycle of ligand-FEP calculations, validating the binding orientation of the different chemical substituents proposed. Subsequent X-ray crystallography of the A2A AR with a low and a high affinity chromone derivative confirmed the predicted binding orientation. The new molecules and structures here reported were driven by free energy calculations, and provide new insights on antagonist binding to the A2A AR, an emerging target in immuno-oncology.

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.

Modulators of CXCR4 and CXCR7/ACKR3 Function.

The two G protein-coupled receptors (GPCRs) C-X-C chemokine receptor type 4 (CXCR4) and atypical chemokine receptor 3 (ACKR3) are part of the class A chemokine GPCR family and represent important drug targets for human immunodeficiency virus (HIV) infection, cancer, and inflammation diseases. CXCR4 is one of only three chemokine receptors with a US Food and Drug Administration approved therapeutic agent, the small-molecule modulator AMD3100. In this review, known modulators of the two receptors are discussed in detail. Initially, the structural relationship between receptors and ligands is reviewed on the basis of common structural motifs and available crystal structures. To date, no atypical chemokine receptor has been crystallized, which makes ligand design and predictions for these receptors more difficult. Next, the selectivity, receptor activation, and the resulting ligand-induced signaling output of chemokines and other peptide ligands are reviewed. Binding of pepducins, a class of lipid-peptides whose basis is the internal loop of a GPCR, to CXCR4 is also discussed. Finally, small-molecule modulators of CXCR4 and ACKR3 are reviewed. These modulators have led to the development of radio- and fluorescently labeled tool compounds, enabling the visualization of ligand binding and receptor characterization both in vitro and in vivo. SIGNIFICANCE STATEMENT: To investigate the pharmacological modulation of CXCR4 and ACKR3, significant effort has been focused on the discovery and development of a range of ligands, including small-molecule modulators, pepducins, and synthetic peptides. Imaging tools, such as fluorescent probes, also play a pivotal role in the field of drug discovery. This review aims to provide an overview of the aforementioned modulators that facilitate the study of CXCR4 and ACKR3 receptors.

Identification of a novel allosteric GLP-1R antagonist HTL26119 using structure- based drug design.

A series of novel allosteric antagonists of the GLP-1 receptor (GLP-1R), exemplified by HTL26119, are described. SBDD approaches were employed to identify HTL26119, exploiting structural understanding of the allosteric binding site of the closely related Glucagon receptor (GCGR) (Jazayeri et al., 2016) and the homology relationships between GCGR and GLP-1R. The region around residue C3476.36b of the GLP-1R receptor represents a key difference from GCGR and was targeted for selectivity for GLP-1R.

Understanding Ligand Binding Selectivity in a Prototypical GPCR Family.

Adenosine receptors are involved in many pathological conditions and are thus promising drug targets. However, developing drugs that target this GPCR subfamily is a challenging task. A number of drug candidates fail due to lack of selectivity which results in unwanted side effects. The extensive structural similarity of adenosine receptors complicates the design of selective ligands. The problem of selective targeting is a general concern in GPCRs, and in this respect adenosine receptors are a prototypical example. Here we use enhanced sampling simulations to decipher the determinants of selectivity of ligands in A2a and A1 adenosine receptors. Our model shows how small differences in the binding pocket and in the water network around the ligand can be leveraged to achieve selectivity.