New insights from structural biology into the druggability of G protein-coupled receptors

Structural-Functional Features of the Thyrotropin Receptor: A Class A G-Protein-Coupled Receptor at Work

The thyroid-stimulating hormone receptor (TSHR) is a member of the glycoprotein hormone receptors, a sub-group of class A G-protein-coupled receptors (GPCRs). TSHR and its endogenous ligand thyrotropin (TSH) are of essential importance for growth and function of the thyroid gland and proper function of the TSH/TSHR system is pivotal for production and release of thyroid hormones. This receptor is also important with respect to pathophysiology, such as autoimmune (including ophthalmopathy) or non-autoimmune thyroid dysfunctions and cancer development. Pharmacological interventions directly targeting the TSHR should provide benefits to disease treatment compared to currently available therapies of dysfunctions associated with the TSHR or the thyroid gland. Upon TSHR activation, the molecular events conveying conformational changes from the extra- to the intracellular side of the cell across the membrane comprise reception, conversion, and amplification of the signal. These steps are highly dependent on structural features of this receptor and its intermolecular interaction partners, e.g., TSH, antibodies, small molecules, G-proteins, or arrestin. For better understanding of signal transduction, pathogenic mechanisms such as autoantibody action and mutational modifications or for developing new pharmacological strategies, it is essential to combine available structural data with functional information to generate homology models of the entire receptor. Although so far these insights are fragmental, in the past few decades essential contributions have been made to investigate in-depth the involved determinants, such as by structure determination via X-ray crystallography. This review summarizes available knowledge (as of December 2016) concerning the TSHR protein structure, associated functional aspects, and based on these insights we suggest several receptor complex models. Moreover, distinct TSHR properties will be highlighted in comparison to other class A GPCRs to understand the molecular activation mechanisms of this receptor comprehensively. Finally, limitations of current knowledge and lack of information are discussed highlighting the need for intensified efforts toward TSHR structure elucidation.

Construction of Structural Mimetics of the Thyrotropin Receptor Intracellular Domain

The signaling of a G protein-coupled receptor (GPCR) is dictated by the complementary responsiveness of interacting intracellular effectors such as G proteins. Many GPCRs are known to couple to more than one G protein subtype and induce a multitude of signaling pathways, although the in vivo relevance of particular pathways is mostly unrecognized. Dissecting GPCR signaling in terms of the pathways that are activated will boost our understanding of the molecular fundamentals of hormone action. The structural determinants governing the selectivity of GPCR/G protein coupling, however, remain obscure. Here, we describe the design of soluble GPCR mimetics to study the details of the interplay between G-proteins and activators. We constructed functional mimetics of the intracellular domain of a model GPCR, the thyrotropin receptor. We based the construction on a unique scaffold, 6-Helix, an artificial protein that was derived from the elements of the trimer-of-hairpins structure of HIV gp41 and represents a bundle of six α-helices. The 6-Helix scaffold, which endowed the substituted thyrotropin receptor intracellular domain elements with spatial constraints analogous to those found in native receptors, enabled the reconstitution of a microdomain that consists of intracellular loops 2 and 3, and is capable of binding and activating Gα-(s). The 6-Helix-based mimetics could be used as a platform to study the molecular basis of GPCR/G protein recognition. Such knowledge could help investigators develop novel therapeutic strategies for GPCR-related disorders by targeting the GPCR/G protein interfaces and counteracting cellular dysfunctions via focused tuning of GPCR signaling.

Deorphaning the Macromolecular Targets of the Natural Anticancer Compound Doliculide

The cyclodepsipeptide doliculide is a marine natural product with strong actin-polymerizing and anticancer activities. Evidence for doliculide acting as a potent and subtype-selective antagonist of prostanoid E receptor 3 (EP3) is presented. Computational target prediction suggested that this membrane receptor is a likely macromolecular target and enabled immediate in vitro validation. This proof-of-concept study demonstrates the in silico deorphanization of phenotypic screening hits as a viable concept for future natural-product-inspired chemical biology and drug discovery efforts.

Native phasing of x-ray free-electron laser data for a G protein-coupled receptor

Serial femtosecond crystallography (SFX) takes advantage of extremely bright and ultrashort pulses produced by x-ray free-electron lasers (XFELs), allowing for the collection of high-resolution diffraction intensities from micrometer-sized crystals at room temperature with minimal radiation damage, using the principle of "diffraction-before-destruction." However, de novo structure factor phase determination using XFELs has been difficult so far. We demonstrate the ability to solve the crystallographic phase problem for SFX data collected with an XFEL using the anomalous signal from native sulfur atoms, leading to a bias-free room temperature structure of the human A_2A adenosine receptor at 1.9 Å resolution. The advancement was made possible by recent improvements in SFX data analysis and the design of injectors and delivery media for streaming hydrated microcrystals. This general method should accelerate structural studies of novel difficult-to-crystallize macromolecules and their complexes.

The use of conformationally thermostabilised GPCRs in drug discovery: application to fragment, structure and biophysical techniques

Recent developments in receptor stabilisation have facilitated major advances in G protein-coupled receptor (GPCR) research, notably structural biology, over the past eight years. Here we review the application of fragment, structure and biophysical techniques using stabilised GPCRs (StaR proteins), and their impact in the drug discovery process. These techniques have, most recently, been utilised in the discovery of the non-alkyne mGlu₅ negative allosteric modulator HTL14242, in addition to the dual orexin receptor antagonist HTL6641, with differentiated residence time kinetics.

Oligomerization of GPCRs involved in endocrine regulation

More than 800 different human membrane-spanning G-protein-coupled receptors (GPCRs) serve as signal transducers at biological barriers. These receptors are activated by a wide variety of ligands such as peptides, ions and hormones, and are able to activate a diverse set of intracellular signaling pathways. GPCRs are of central importance in endocrine regulation, which underpins the significance of comprehensively studying these receptors and interrelated systems. During the last decade, the capacity for multimerization of GPCRs was found to be a common and functionally relevant property. The interaction between GPCR monomers results in higher order complexes such as homomers (identical receptor subtype) or heteromers (different receptor subtypes), which may be present in a specific and dynamic monomer/oligomer equilibrium. It is widely accepted that the oligomerization of GPCRs is a mechanism for determining the fine-tuning and expansion of cellular processes by modification of ligand action, expression levels, and related signaling outcome. Accordingly, oligomerization provides exciting opportunities to optimize pharmacological treatment with respect to receptor target and tissue selectivity or for the development of diagnostic tools. On the other hand, GPCR heteromerization may be a potential reason for the undesired side effects of pharmacological interventions, faced with numerous and common mutual signaling modifications in heteromeric constellations. Finally, detailed deciphering of the physiological occurrence and relevance of specific GPCR/GPCR-ligand interactions poses a future challenge. This review will tackle the aspects of GPCR oligomerization with specific emphasis on family A GPCRs involved in endocrine regulation, whereby only a subset of these receptors will be discussed in detail.

Structurally Enabled Discovery of Adenosine A2A Receptor Antagonists

Over the past decade there has been a revolution in the field of G protein-coupled receptor (GPCR) structural biology. Many years of innovative research from different areas have come together to fuel this significant change in the fortunes of this field, which for many years was characterized by the paucity of high-resolution structures. The determination to succeed has been in part due to the recognized importance of these proteins as drug targets, and although the pharmaceutical industry has been focusing on these receptors, it can be justifiably argued and demonstrated that many of the approved and commercially successful GPCR drugs can be significantly improved to increase efficacy and/or reduce undesired side effects. In addition, many validated targets in this class remain to be drugged. It is widely recognized that application of structure-based drug design approaches can help medicinal chemists a long way toward discovering better drugs. The achievement of structural biologists in providing high-resolution insight is beginning to transform drug discovery efforts, and there are a number of GPCR drugs that have been discovered by use of structural information that are in clinical development. This review aims to highlight the key developments that have brought success to GPCR structure resolution efforts and exemplify the practical application of structural information for the discovery of adenosine A_2A receptor antagonists that have potential to treat multiple conditions.

ALIX Regulates the Ubiquitin-Independent Lysosomal Sorting of the P2Y1 Purinergic Receptor via a YPX3L Motif

Endocytic sorting and lysosomal degradation are integral to the regulation of G protein-coupled receptor (GPCR) function. Upon ligand binding, classical GPCRs are activated, internalized and recycled or sorted to lysosomes for degradation, a process that requires receptor ubiquitination. However, recent studies have demonstrated that numerous GPCRs are sorted to lysosomes independent of receptor ubiquitination. Here, we describe an ubiquitin-independent lysosomal sorting pathway for the purinergic GPCR P2Y₁. After activation, P2Y₁ sorts to lysosomes for degradation independent of direct ubiquitination that is mediated by a YPX₃L motif within the second intracellular loop that serves as a binding site for the adaptor protein ALIX. Depletion of ALIX or site-directed mutation of the YPX₃L motif inhibits P2Y₁ sorting into the lumen of multivesicular endosomes/lysosomes and degradation. These findings confirm the function of YPX₃L motifs as lysosomal targeting sequences for GPCRs and demonstrate that ALIX mediates the ubiquitin-independent degradation of certain GPCRs.

G-Protein/β-Arrestin-Linked Fluctuating Network of G-Protein-Coupled Receptors for Predicting Drug Efficacy and Bias Using Short-Term Molecular Dynamics Simulation

The efficacy and bias of signal transduction induced by a drug at a target protein are closely associated with the benefits and side effects of the drug. In particular, partial agonist activity and G-protein/β-arrestin-biased agonist activity for the G-protein-coupled receptor (GPCR) family, the family with the most target proteins of launched drugs, are key issues in drug discovery. However, designing GPCR drugs with appropriate efficacy and bias is challenging because the dynamic mechanism of signal transduction induced by ligand—receptor interactions is complicated. Here, we identified the G-protein/β-arrestin-linked fluctuating network, which initiates large-scale conformational changes, using sub-microsecond molecular dynamics (MD) simulations of the β₂-adrenergic receptor (β₂AR) with a diverse collection of ligands and correlation analysis of their G protein/β-arrestin efficacy. The G-protein-linked fluctuating network extends from the ligand-binding site to the G-protein-binding site through the connector region, and the β-arrestin-linked fluctuating network consists of the NPxxY motif and adjacent regions. We confirmed that the averaged values of fluctuation in the fluctuating network detected are good quantitative indexes for explaining G protein/β-arrestin efficacy. These results indicate that short-term MD simulation is a practical method to predict the efficacy and bias of any compound for GPCRs.

Application of advanced X-ray methods in life sciences

BACKGROUND

Synchrotron radiation (SR) sources provide diverse X-ray methods for the investigation of structure-function relationships in biological macromolecules.

SCOPE OF REVIEW

Recent developments in SR sources and in the X-ray tools they offer for life sciences are reviewed. Specifically, advances in macromolecular crystallography, small angle X-ray solution scattering, X-ray absorption and fluorescence spectroscopy, and imaging are discussed with examples.

MAJOR CONCLUSIONS

SR sources offer a range of X-ray techniques that can be used in a complementary fashion in studies of biological systems at a wide range of resolutions from atomic to cellular scale. Emerging applications of X-ray techniques include the characterization of disordered proteins, noncrystalline and nonequilibrium systems, elemental imaging of tissues, cells and organs, and detection of time-resolved changes in molecular structures.

GENERAL SIGNIFICANCE

X-ray techniques are in the center of hybrid approaches that are used to gain insight into complex problems relating to biomolecular mechanisms, disease and possible therapeutic solutions. This article is part of a Special Issue entitled "Science for Life". Guest Editors: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Guided reconstitution of membrane protein fragments

Structural analysis by NMR of G protein-coupled receptors (GPCRs) has proven to be extremely challenging. To reduce the number of peaks in the NMR spectra by segmentally labeling a GPCR, we have developed a Guided Reconstitution method that includes the use of charged residues and Cys activation to drive heterodimeric disulfide bond formation. Three different cysteine-activating reagents: 5-5'-dithiobis(2-nitrobenzoic acid) [DTNB], 2,2'-dithiobis(5-nitropyridine) [DTNP], and 4,4'-dipyridyl disulfide [4-PDS] were analyzed to determine their efficiency in heterodimer formation at different pHs. Short peptides representing the N-terminal (NT) and C-terminal (CT) regions of the first extracellular loop (EL1) of Ste2p, the Saccharomyces cerevisiae alpha-factor mating receptor, were activated using these reagents and the efficiencies of activation and rates of heterodimerization were analyzed. Activation of NT peptides with DTNP and 4-PDS resulted in about 60% yield, but heterodimerization was rapid and nearly quantitative. Double transmembrane domain protein fragments were biosynthesized and used in Guided Reconstitution reactions. A 102-residue fragment, 2TM-tail [Ste2p(G31-I120C)], was heterodimerized with CT-EL1-tail(DTNP) at pH 4.6 with a yield of ∼75%. A 132-residue fragment, 2TMlong-tail [Ste2p(M1-I120C)], was expressed in both unlabeled and (15)N-labeled forms and used with a peptide comprising the third transmembrane domain, to generate a 180-residue segmentally labeled 3TM protein that was found to be segmentally labeled using [(15)N,(1)H]-HSQC analysis. Our data indicate that the Guided Reconstitution method would be applicable to the segmental labeling of a membrane protein with 3 transmembrane domains and may prove useful in the preparation of an intact reconstituted GPCR for use in biophysical analysis and structure determination.

GPCR-Bench: A Benchmarking Set and Practitioners' Guide for G Protein-Coupled Receptor Docking

Virtual screening is routinely used to discover new ligands and in particular new ligand chemotypes for G protein-coupled receptors (GPCRs). To prepare for a virtual screen, we often tailor a docking protocol that will enable us to select the best candidates for further screening. To aid this, we created GPCR-Bench, a publically available docking benchmarking set in the spirit of the DUD and DUD-E reference data sets for validation studies, containing 25 nonredundant high-resolution GPCR costructures with an accompanying set of diverse ligands and computational decoy molecules for each target. Benchmarking sets are often used to compare docking protocols; however, it is important to evaluate docking methods not by "retrospective" hit rates but by the actual likelihood that they will produce novel prospective hits. Therefore, docking protocols must not only rank active molecules highly but also produce good poses that a chemist will select for purchase and screening. Currently, no simple objective machine-scriptable function exists that can do this; instead, docking hit lists must be subjectively examined in a consistent way to compare between docking methods. We present here a case study highlighting considerations we feel are of importance when evaluating a method, intended to be useful as a practitioners' guide.

Complete Reversible Refolding of a G-Protein Coupled Receptor on a Solid Support

The factors defining the correct folding and stability of integral membrane proteins are poorly understood. Folding of only a few select membrane proteins has been scrutinised, leaving considerable deficiencies in knowledge for large protein families, such as G protein coupled receptors (GPCRs). Complete reversible folding, which is problematic for any membrane protein, has eluded this dominant receptor family. Moreover, attempts to recover receptors from denatured states are inefficient, yielding at best 40–70% functional protein. We present a method for the reversible unfolding of an archetypal family member, the β₁-adrenergic receptor, and attain 100% recovery of the folded, functional state, in terms of ligand binding, compared to receptor which has not been subject to any unfolding and retains its original, folded structure. We exploit refolding on a solid support, which could avoid unwanted interactions and aggregation that occur in bulk solution. We determine the changes in structure and function upon unfolding and refolding. Additionally, we employ a method that is relatively new to membrane protein folding; pulse proteolysis. Complete refolding of β₁-adrenergic receptor occurs in n-decyl-β-D-maltoside (DM) micelles from a urea-denatured state, as shown by regain of its original helical structure, ligand binding and protein fluorescence. The successful refolding strategy on a solid support offers a defined method for the controlled refolding and recovery of functional GPCRs and other membrane proteins that suffer from instability and irreversible denaturation once isolated from their native membranes.

7TM Domain Structure of Adhesion GPCRs

Schematic presentation of the overall adhesion G Protein-Coupled Receptor (aGPCR) structure and functional domains, covering an extracellular N-terminal fragment (NTF), a membrane-spanning C-terminal fragment (CTF) and a GPCR proteolysis site (GPS). (Left side) aGPCR model constructed based on the seven-transmembrane (7TM) structure (blue) of secretin family glucagon receptor (GCGR) (PDB, 4L6R) [11] and the GPCR autoproteolysis inducing (GAIN) domain (magenta) structure of latrophilin 1 (PDB, 4DLQ) [9]. The β-13 strand residues are depicted in green. (Right side) The experimentally validated full-length secretin family GCGR structure combining structural and experimental information from the GCGR 7TM crystal structure (PDB, 4L6R) (blue), the GCGR extracellular domain (ECD) structure (PDB, 4ERS) (magenta) and the ECD structure of glucagon-like peptide-1 (GLP-1)-bound glucagon-like peptide-1 receptor (GLP-1R) (PDB, 3IOL) (green), complemented by site-directed mutagenesis, electron microscopy (EM), hydrogen-deuterium exchange (HDX) and cross-linking studies [11-13]) Despite the recent breakthroughs in the elucidation of the three-dimensional structures of the seven transmembrane (7TM) domain of the G protein-coupled receptor (GPCR) superfamily, a corresponding structure of a member of the adhesion GPCR (aGPCR) family has not yet been solved. In this chapter, we give an overview of the current knowledge of the 7TM domain of aGPCRs by comparative structure-based sequence similarity analyses between aGPCRs and GPCRs with known crystal structure. Of the GPCR superfamily, only the secretin family shares some sequence similarity with aGPCRs. This chapter will therefore emphasize on the comparison of these two GPCR families. Two 7TM domain structures of secretin family GPCRs are known that provide insight into the structure-function relationships of conserved sequence motifs that play important roles and are also present in most aGPCRs. This suggests that the 7TM domains of aGPCRs and secretin family GPCRs share a similar structural fold and that the conserved residues in both families may be involved in similar intermolecular interaction networks and facilitate similar conformational changes. Comparison of the residues that line the large peptide hormone binding pocket in the 7TM domain of secretin family GPCRs with corresponding residues in aGPCRs indicates that in the latter, the corresponding pocket in the 7TM domain is relatively hydrophobic and may be even larger. Improved knowledge on these conserved sequence motifs will help to understand the interactions of the aGPCR 7TM domain with ligands and gain insight into the activation mechanism of aGPCRs.

Application of GPCR Structures for Modelling of Free Fatty Acid Receptors

Five G protein-coupled receptors (GPCRs) have been identified to be activated by free fatty acids (FFA). Among them, FFA1 (GPR40) and FFA4 (GPR120) bind long-chain fatty acids, FFA2 (GPR43) and FFA3 (GPR41) bind short-chain fatty acids and GPR84 binds medium-chain fatty acids. Free fatty acid receptors have now emerged as potential targets for the treatment of diabetes, obesity and immune diseases. The recent progress in crystallography of GPCRs has now enabled the elucidation of the structure of FFA1 and provided reliable templates for homology modelling of other FFA receptors. Analysis of the crystal structure and improved homology models, along with mutagenesis data and structure activity, highlighted an unusual arginine charge-pairing interaction in FFA1-3 for receptor modulation, distinct structural features for ligand binding to FFA1 and FFA4 and an arginine of the second extracellular loop as a possible anchoring point for FFA at GPR84. Structural data will be helpful for searching novel small-molecule modulators at the FFA receptors.

Beyond Membrane Protein Structure: Drug Discovery, Dynamics and Difficulties

Most of the previous content of this book has focused on obtaining the structures of membrane proteins. In this chapter we explore how those structures can be further used in two key ways. The first is their use in structure based drug design (SBDD) and the second is how they can be used to extend our understanding of their functional activity via the use of molecular dynamics. Both aspects now heavily rely on computations. This area is vast, and alas, too large to consider in depth in a single book chapter. Thus where appropriate we have referred the reader to recent reviews for deeper assessment of the field. We discuss progress via the use of examples from two main drug target areas; G-protein coupled receptors (GPCRs) and ion channels. We end with a discussion of some of the main challenges in the area.

The Fragment Molecular Orbital Method Reveals New Insight into the Chemical Nature of GPCR–Ligand Interactions

Our interpretation of ligand-protein interactions is often informed by high-resolution structures, which represent the cornerstone of structure-based drug design. However, visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum Mechanics approaches are often too computationally expensive, but one method, Fragment Molecular Orbital (FMO), offers an excellent compromise and has the potential to reveal key interactions that would otherwise be hard to detect. To illustrate this, we have applied the FMO method to 18 Class A GPCR-ligand crystal structures, representing different branches of the GPCR genome. Our work reveals key interactions that are often omitted from structure-based descriptions, including hydrophobic interactions, nonclassical hydrogen bonds, and the involvement of backbone atoms. This approach provides a more comprehensive picture of receptor-ligand interactions than is currently used and should prove useful for evaluation of the chemical nature of ligand binding and to support structure-based drug design.

Fragment-Based Discovery of Subtype-Selective Adenosine Receptor Ligands from Homology Models

Fragment-based lead discovery (FBLD) holds great promise for drug discovery, but applications to G protein-coupled receptors (GPCRs) have been limited by a lack of sensitive screening techniques and scarce structural information. If virtual screening against homology models of GPCRs could be used to identify fragment ligands, FBLD could be extended to numerous important drug targets and contribute to efficient lead generation. Access to models of multiple receptors may further enable the discovery of fragments that bind specifically to the desired target. To investigate these questions, we used molecular docking to screen >500 000 fragments against homology models of the A3 and A1 adenosine receptors (ARs) with the goal to discover A3AR-selective ligands. Twenty-one fragments with predicted A3AR-specific binding were evaluated in live-cell fluorescence-based assays; of eight verified ligands, six displayed A3/A1 selectivity, and three of these had high affinities ranging from 0.1 to 1.3 μM. Subsequently, structure-guided fragment-to-lead optimization led to the identification of a >100-fold-selective antagonist with nanomolar affinity from commercial libraries. These results highlight that molecular docking screening can guide fragment-based discovery of selective ligands even if the structures of both the target and antitarget receptors are unknown. The same approach can be readily extended to a large number of pharmaceutically important targets.

How Can Mutations Thermostabilize G-Protein-Coupled Receptors?

Structures of over 30 different G-protein-coupled receptors (GPCRs) have advanced our understanding of cell signaling and have provided a foundation for structure-guided drug design. This exciting progress has required the development of three complementary methods to facilitate GPCR crystallization, one of which is the thermostabilization of receptors by systematic mutagenesis. However, the reason why a particular mutation, or combination of mutations, stabilizes the receptor is not always evident from a static crystal structure. Molecular dynamics (MD) simulations have been used to identify and estimate the energetic factors that affect thermostability through comparing the dynamics of the thermostabilized receptors with structure-based models of the wild-type receptor. The data indicate that receptors are stabilized through a combination of factors, including an increase in receptor rigidity, a decrease in collective motion, reduced stress at specific residues, and the presence of ordered water molecules. Predicting thermostabilizing mutations computationally represents a major challenge for the field.

The α-arrestin ARRDC3 mediates ALIX ubiquitination and G protein-coupled receptor lysosomal sorting

The novel ALIX-dependent GPCR sorting pathway is regulated by the a-arrestin ARRDC3. A critical role is also shown for the E3 ubiquitin ligase WWP2 in regulation of ALIX ubiquitination and lysosomal sorting of GPCRs.

The sorting of G protein–coupled receptors (GPCRs) to lysosomes is critical for proper signaling and cellular responses. We previously showed that the adaptor protein ALIX regulates lysosomal degradation of protease-activated receptor-1 (PAR1), a GPCR for thrombin, independent of ubiquitin-binding ESCRTs and receptor ubiquitination. However, the mechanisms that regulate ALIX function during PAR1 lysosomal sorting are not known. Here we show that the mammalian α-arrestin arrestin domain–containing protein-3 (ARRDC3) regulates ALIX function in GPCR sorting via ubiquitination. ARRDC3 colocalizes with ALIX and is required for PAR1 sorting at late endosomes and degradation. Depletion of ARRDC3 by small interfering RNA disrupts ALIX interaction with activated PAR1 and the CHMP4B ESCRT-III subunit, suggesting that ARRDC3 regulates ALIX activity. We found that ARRDC3 is required for ALIX ubiquitination induced by activation of PAR1. A screen of nine mammalian NEDD4-family E3 ubiquitin ligases revealed a critical role for WWP2. WWP2 interacts with ARRDC3 and not ALIX. Depletion of WWP2 inhibited ALIX ubiquitination and blocked ALIX interaction with activated PAR1 and CHMP4B. These findings demonstrate a new role for the α-arrestin ARRDC3 and the E3 ubiquitin ligase WWP2 in regulation of ALIX ubiquitination and lysosomal sorting of GPCRs.

Biased signalling: the instinctive skill of the cell in the selection of appropriate signalling pathways

GPCRs (G-protein-coupled receptors) are members of a family of proteins which are generally regarded as the largest group of therapeutic drug targets. Ligands of GPCRs do not usually activate all cellular signalling pathways linked to a particular seven-transmembrane receptor in a uniform manner. The fundamental idea behind this concept is that each ligand has its own ability, while interacting with the receptor, to activate different signalling pathways (or a particular set of signalling pathways) and it is this concept which is known as biased signalling. The importance of biased signalling is that it may selectively activate biological responses to favour therapeutically beneficial signalling pathways and to avoid adverse effects. There are two levels of biased signalling. First, bias can arise from the ability of GPCRs to couple to a subset of the available G-protein subtypes: Gαs, Gαq/11, Gαi/o or Gα12/13. These subtypes produce the diverse effects of GPCRs by targeting different effectors. Secondly, biased GPCRs may differentially activate G-proteins or β-arrestins. β-Arrestins are ubiquitously expressed and function to terminate or inhibit classic G-protein signalling and initiate distinct β-arrestin-mediated signalling processes. The interplay of G-protein and β-arrestin signalling largely determines the cellular consequences of the administration of GPCR-targeted drugs. In the present review, we highlight the particular functionalities of biased signalling and discuss its biological effects subsequent to GPCR activation. We consider that biased signalling is potentially allowing a choice between signalling through ‘beneficial’ pathways and the avoidance of ‘harmful’ ones.

The impact of binding thermodynamics on medicinal chemistry optimizations

Ligand binding thermodynamics has been attracted considerable interest in the past decade owing to the recognized relation between binding thermodynamic profile and the physicochemical and druglike properties of compounds. In this review, the relation between optimization strategies and ligand properties is presented based on the structural and thermodynamic analysis of ligand–protein complex formation. The control of the binding thermodynamic profile is beneficial for the balanced affinity and physicochemical properties of drug candidates, and early phase optimization gives more opportunity to this control.

From G Protein-coupled Receptor Structure Resolution to Rational Drug Design

A number of recent technical solutions have led to significant advances in G protein-coupled receptor (GPCR) structural biology. Apart from a detailed mechanistic view of receptor activation, the new structures have revealed novel ligand binding sites. Together, these insights provide avenues for rational drug design to modulate the activities of these important drug targets. The application of structural data to GPCR drug discovery ushers in an exciting era with the potential to improve existing drugs and discover new ones. In this review, we focus on technical solutions that have accelerated GPCR crystallography as well as some of the salient findings from structures that are relevant to drug discovery. Finally, we outline some of the approaches used in GPCR structure based drug design.

Systems Pharmacology Links GPCRs with Retinal Degenerative Disorders

In most biological systems, second messengers and their key regulatory and effector proteins form links between multiple cellular signaling pathways. Such signaling nodes can integrate the deleterious effects of genetic aberrations, environmental stressors, or both in complex diseases, leading to cell death by various mechanisms. Here we present a systems (network) pharmacology approach that, together with transcriptomics analyses, was used to identify different G protein-coupled receptors that experimentally protected against cellular stress and death caused by linked signaling mechanisms. We describe the application of this concept to degenerative and diabetic retinopathies in appropriate mouse models as an example. Systems pharmacology also provides an attractive framework for devising strategies to combat complex diseases by using (repurposing) US Food and Drug Administration-approved pharmacological agents.

Clustered distribution of natural product leads of drugs in the chemical space as influenced by the privileged target-sites

Some natural product leads of drugs (NPLDs) have been found to congregate in the chemical space. The extent, detailed patterns, and mechanisms of this congregation phenomenon have not been fully investigated and their usefulness for NPLD discovery needs to be more extensively tested. In this work, we generated and evaluated the distribution patterns of 442 NPLDs of 749 pre-2013 approved and 263 clinical trial small molecule drugs in the chemical space represented by the molecular scaffold and fingerprint trees of 137,836 non-redundant natural products. In the molecular scaffold trees, 62.7% approved and 37.4% clinical trial NPLDs congregate in 62 drug-productive scaffolds/scaffold-branches. In the molecular fingerprint tree, 82.5% approved and 63.0% clinical trial NPLDs are clustered in 60 drug-productive clusters (DCs) partly due to their preferential binding to 45 privileged target-site classes. The distribution patterns of the NPLDs are distinguished from those of the bioactive natural products. 11.7% of the NPLDs in these DCs have remote-similarity relationship with the nearest NPLD in their own DC. The majority of the new NPLDs emerge from preexisting DCs. The usefulness of the derived knowledge for NPLD discovery was demonstrated by the recognition of the new NPLDs of 2013-2014 approved drugs.

Comprehensive, structurally-informed alignment and phylogeny of vertebrate biogenic amine receptors

Biogenic amine receptors play critical roles in regulating behavior and physiology in both vertebrates and invertebrates, particularly within the central nervous system. Members of the G-protein coupled receptor (GPCR) family, these receptors interact with endogenous bioamine ligands such as dopamine, serotonin, and epinephrine, and are targeted by a wide array of pharmaceuticals. Despite the clear clinical and biological importance of these receptors, their evolutionary history remains poorly characterized. In particular, the relationships among biogenic amine receptors and any specific evolutionary constraints acting within distinct receptor subtypes are largely unknown. To advance and facilitate studies in this receptor family, we have constructed a comprehensive, high-quality sequence alignment of vertebrate biogenic amine receptors. In particular, we have integrated a traditional multiple sequence approach with robust structural domain predictions to ensure that alignment columns accurately capture the highly-conserved GPCR structural domains, and we demonstrate how ignoring structural information produces spurious inferences of homology. Using this alignment, we have constructed a structurally-partitioned maximum-likelihood phylogeny from which we deduce novel biogenic amine receptor relationships and uncover previously unrecognized lineage-specific receptor clades. Moreover, we find that roughly 1% of the 3039 sequences in our final alignment are either misannotated or unclassified, and we propose updated classifications for these receptors. We release our comprehensive alignment and its corresponding phylogeny as a resource for future research into the evolution and diversification of biogenic amine receptors.

Structure of Class B GPCRs: new horizons for drug discovery

Class B GPCRs of the secretin family are important drug targets in many human diseases including diabetes, neurodegeneration, cardiovascular disease and psychiatric disorders. X-ray crystal structures for the glucagon receptor and corticotropin-releasing factor receptor 1 have now been published. In this review, we analyse the new structures and how they compare with each other and with Class A and F receptors. We also consider the differences in druggability and possible similarity in the activation mechanisms. Finally, we discuss the potential for the design of small-molecule modulators for these important targets in drug discovery. This new structural insight allows, for the first time, structure-based drug design methods to be applied to Class B GPCRs.

Docking and Virtual Screening Strategies for GPCR Drug Discovery

Progress in structure determination of G protein-coupled receptors (GPCRs) has made it possible to apply structure-based drug design (SBDD) methods to this pharmaceutically important target class. The quality of GPCR structures available for SBDD projects fall on a spectrum ranging from high resolution crystal structures (<2 Å), where all water molecules in the binding pocket are resolved, to lower resolution (>3 Å) where some protein residues are not resolved, and finally to homology models that are built using distantly related templates. Each GPCR project involves a distinct set of opportunities and challenges, and requires different approaches to model the interaction between the receptor and the ligands. In this review we will discuss docking and virtual screening to GPCRs, and highlight several refinement and post-processing steps that can be used to improve the accuracy of these calculations. Several examples are discussed that illustrate specific steps that can be taken to improve upon the docking and virtual screening accuracy. While GPCRs are a unique target class, many of the methods and strategies outlined in this review are general and therefore applicable to other protein families.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.

Docking and free energy perturbation studies of ligand binding in the kappa opioid receptor

The kappa opioid receptor (KOR) is an important target for pain and depression therapeutics that lack harmful and addictive qualities of existing medications. We present a model for the binding of morphinan ligands and JDTic to the JDTic/KOR crystal structure based on an atomic level description of the water structure within its active site. The model contains two key interaction motifs that are supported by experimental evidence. The first is the formation of a salt bridge between the ligand and Asp 138(3.32) in transmembrane domain (TM) 3. The second is the stabilization by the ligand of two high energy, isolated, and ice-like waters near TM5 and TM6. This model is incorporated via energetic terms into a new empirical scoring function, WScore, designed to assess interactions between ligands and localized water in a binding site. Pairing WScore with the docking program Glide discriminates known active KOR ligands from large sets of decoy molecules much better than Glide's older generation scoring functions, SP and XP. We also use rigorous free energy perturbation calculations to provide evidence for the proposed mechanism of interaction between ligands and KOR. The molecular description of ligand binding in KOR should provide a good starting point for future drug discovery efforts for this receptor.

W246(6.48) opens a gate for a continuous intrinsic water pathway during activation of the adenosine A2A receptor

The question how G-protein-coupled receptors transduce an extracellular signal by a sequence of transmembrane conformational transitions into an intracellular response remains to be solved at molecular detail. Herein, we use molecular dynamics simulations to reveal distinct conformational transitions of the adenosine A2A receptor, and we found that the conserved W246(6.48) residue in transmembrane helix TM6 performs a key rotamer toggle switch. Agonist binding induces the sidechain of W246(6.48) to fluctuate between two distinct conformations enabling the diffusion of water molecules from the bulk into the center of the receptor. After passing the W246(6.48) gate, the internal water molecules induce another conserved residue, Y288(7.53), to switch to a distinct rotamer conformation establishing a continuous transmembrane water pathway. Further, structural changes of TM6 and TM7 induce local structural changes of the adjacent lipid bilayer.

Endothelial C-type natriuretic peptide maintains vascular homeostasis

The endothelium plays a fundamental role in maintaining vascular homeostasis by releasing factors that regulate local blood flow, systemic blood pressure, and the reactivity of leukocytes and platelets. Accordingly, endothelial dysfunction underpins many cardiovascular diseases, including hypertension, myocardial infarction, and stroke. Herein, we evaluated mice with endothelial-specific deletion of Nppc, which encodes C-type natriuretic peptide (CNP), and determined that this mediator is essential for multiple aspects of vascular regulation. Specifically, disruption of CNP leads to endothelial dysfunction, hypertension, atherogenesis, and aneurysm. Moreover, we identified natriuretic peptide receptor-C (NPR-C) as the cognate receptor that primarily underlies CNP-dependent vasoprotective functions and developed small-molecule NPR-C agonists to target this pathway. Administration of NPR-C agonists promotes a vasorelaxation of isolated resistance arteries and a reduction in blood pressure in wild-type animals that is diminished in mice lacking NPR-C. This work provides a mechanistic explanation for genome-wide association studies that have linked the NPR-C (Npr3) locus with hypertension by demonstrating the importance of CNP/NPR-C signaling in preserving vascular homoeostasis. Furthermore, these results suggest that the CNP/NPR-C pathway has potential as a disease-modifying therapeutic target for cardiovascular disorders.

Structure-based drug design for G protein-coupled receptors

Our understanding of the structural biology of G protein-coupled receptors has undergone a transformation over the past 5 years. New protein-ligand complexes are described almost monthly in high profile journals. Appreciation of how small molecules and natural ligands bind to their receptors has the potential to impact enormously how medicinal chemists approach this major class of receptor targets. An outline of the key topics in this field and some recent examples of structure- and fragment-based drug design are described. A table is presented with example views of each G protein-coupled receptor for which there is a published X-ray structure, including interactions with small molecule antagonists, partial and full agonists. The possible implications of these new data for drug design are discussed.

Improving the apo-state detergent stability of NTS1 with CHESS for pharmacological and structural studies

The largest single class of drug targets is the G protein-coupled receptor (GPCR) family. Modern high-throughput methods for drug discovery require working with pure protein, but this has been a challenge for GPCRs, and thus the success of screening campaigns targeting soluble, catalytic protein domains has not yet been realized for GPCRs. Therefore, most GPCR drug screening has been cell-based, whereas the strategy of choice for drug discovery against soluble proteins is HTS using purified proteins coupled to structure-based drug design. While recent developments are increasing the chances of obtaining GPCR crystal structures, the feasibility of screening directly against purified GPCRs in the unbound state (apo-state) remains low. GPCRs exhibit low stability in detergent micelles, especially in the apo-state, over the time periods required for performing large screens. Recent methods for generating detergent-stable GPCRs, however, offer the potential for researchers to manipulate GPCRs almost like soluble enzymes, opening up new avenues for drug discovery. Here we apply cellular high-throughput encapsulation, solubilization and screening (CHESS) to the neurotensin receptor 1 (NTS1) to generate a variant that is stable in the apo-state when solubilized in detergents. This high stability facilitated the crystal structure determination of this receptor and also allowed us to probe the pharmacology of detergent-solubilized, apo-state NTS1 using robotic ligand binding assays. NTS1 is a target for the development of novel antipsychotics, and thus CHESS-stabilized receptors represent exciting tools for drug discovery.

Alpha-bulges in G protein-coupled receptors

Agonist binding is related to a series of motions in G protein-coupled receptors (GPCRs) that result in the separation of transmembrane helices III and VI at their cytosolic ends and subsequent G protein binding. A large number of smaller motions also seem to be associated with activation. Most helices in GPCRs are highly irregular and often contain kinks, with extensive literature already available about the role of prolines in kink formation and the precise function of these kinks. GPCR transmembrane helices also contain many α-bulges. In this article we aim to draw attention to the role of these α-bulges in ligand and G-protein binding, as well as their role in several aspects of the mobility associated with GPCR activation. This mobility includes regularization and translation of helix III in the extracellular direction, a rotation of the entire helix VI, an inward movement of the helices near the extracellular side, and a concerted motion of the cytosolic ends of the helices that makes their orientation appear more circular and that opens up space for the G protein to bind. In several cases, α-bulges either appear or disappear as part of the activation process.

An unbiased method to build benchmarking sets for ligand-based virtual screening and its application to GPCRs

Benchmarking data sets have become common in recent years for the purpose of virtual screening, though the main focus had been placed on the structure-based virtual screening (SBVS) approaches. Due to the lack of crystal structures, there is great need for unbiased benchmarking sets to evaluate various ligand-based virtual screening (LBVS) methods for important drug targets such as G protein-coupled receptors (GPCRs). To date these ready-to-apply data sets for LBVS are fairly limited, and the direct usage of benchmarking sets designed for SBVS could bring the biases to the evaluation of LBVS. Herein, we propose an unbiased method to build benchmarking sets for LBVS and validate it on a multitude of GPCRs targets. To be more specific, our methods can (1) ensure chemical diversity of ligands, (2) maintain the physicochemical similarity between ligands and decoys, (3) make the decoys dissimilar in chemical topology to all ligands to avoid false negatives, and (4) maximize spatial random distribution of ligands and decoys. We evaluated the quality of our Unbiased Ligand Set (ULS) and Unbiased Decoy Set (UDS) using three common LBVS approaches, with Leave-One-Out (LOO) Cross-Validation (CV) and a metric of average AUC of the ROC curves. Our method has greatly reduced the “artificial enrichment” and “analogue bias” of a published GPCRs benchmarking set, i.e., GPCR Ligand Library (GLL)/GPCR Decoy Database (GDD). In addition, we addressed an important issue about the ratio of decoys per ligand and found that for a range of 30 to 100 it does not affect the quality of the benchmarking set, so we kept the original ratio of 39 from the GLL/GDD.