Two distinct conformations of helix 6 observed in antagonist-bound structures of a beta1-adrenergic receptor

The role of ligands on the equilibria between functional states of a G protein-coupled receptor

G protein-coupled receptors exhibit a wide variety of signaling behaviors in response to different ligands. When a small label was incorporated on the cytosolic interface of transmembrane helix 6 (Cys-265), (19)F NMR spectra of the β2 adrenergic receptor (β2AR) reconstituted in maltose/neopentyl glycol detergent micelles revealed two distinct inactive states, an activation intermediate state en route to activation, and, in the presence of a G protein mimic, a predominant active state. Analysis of the spectra as a function of temperature revealed that for all ligands, the activation intermediate is entropically favored and enthalpically disfavored. β2AR enthalpy changes toward activation are notably lower than those observed with rhodopsin, a likely consequence of basal activity and the fact that the ionic lock and other interactions stabilizing the inactive state of β2AR are weaker. Positive entropy changes toward activation likely reflect greater mobility (configurational entropy) in the cytoplasmic domain, as confirmed through an order parameter analysis. Ligands greatly influence the overall changes in enthalpy and entropy of the system and the corresponding changes in population and amplitude of motion of given states, suggesting a complex landscape of states and substates.

Molecular modeling of vasopressin receptor and in silico screening of V1b receptor antagonists

INTRODUCTION

G protein-coupled receptors (GPCRs) are integral membrane proteins which contain seven-transmembrane-spanning alpha-helices. GPCR-mediated signaling has been associated with various human diseases, positioning GPCRs as attractive targets in the drug discovery field. Recently, through advances in protein engineering and crystallography, the number of resolved GPCR structures has increased dramatically. This growing availability of GPCR structures has greatly accelerated structure-based drug design (SBDD) and in silico screening for GPCR-targeted drug discovery.

AREAS COVERED

The authors introduce the current status of X-ray crystallography of GPCRs and what has been revealed from the resolved crystal structures. They also review the recent advances in SBDD and in silico screening for GPCR-targeted drug discovery and discuss a docking study, using homology modeling, with the discovery of potent antagonists of the vasopressin 1b receptor.

EXPERT OPINION

Several innovative protein engineering techniques and crystallographic methods have greatly accelerated SBDD, not only for already-resolved GPCRs but also for those structures which remain unclear. These technological advances are expected to enable the determination of GPCR-fragment complexes, making it practical to perform fragment-based drug discovery. This paves the way for a new era of GPCR-targeted drug discovery.

The muscarinic M3 acetylcholine receptor exists as two differently sized complexes at the plasma membrane

The literature on GPCR (G-protein-coupled receptor) homo-oligomerization encompasses conflicting views that range from interpretations that GPCRs must be monomeric, through comparatively newer proposals that they exist as dimers or higher-order oligomers, to suggestions that such quaternary structures are rather ephemeral or merely accidental and may serve no functional purpose. In the present study we use a novel method of FRET (Förster resonance energy transfer) spectrometry and controlled expression of energy donor-tagged species to show that M3Rs (muscarinic M3 acetylcholine receptors) at the plasma membrane exist as stable dimeric complexes, a large fraction of which interact dynamically to form tetramers without the presence of trimers, pentamers, hexamers etc. That M3R dimeric units interact dynamically was also supported by co-immunoprecipitation of receptors synthesized at distinct times. On the basis of all these findings, we propose a conceptual framework that may reconcile the conflicting views on the quaternary structure of GPCRs.

Stabilizing membrane proteins through protein engineering

Integral membrane proteins (IMPs) are crucial components of all cells but are difficult to study in vitro because they are generally unstable when removed from their native membranes using detergents. Despite the major biomedical relevance of IMPs, less than 1% of Protein Data Bank (PDB) entries are IMP structures, reflecting the technical gap between studies of soluble proteins compared to IMPs. Stability can be engineered into IMPs by inserting stabilizing mutations, thereby generating proteins that can be successfully applied to biochemical and structural studies when solubilized in detergent micelles. The identification of stabilizing mutations is not trivial, and this review will focus on the methods that have been used to identify stabilized membrane proteins, including alanine scanning and screening, directed evolution and computational design.

Study of structural dynamics of ligand-activated membrane receptors by means of principal component analysis

The structural dynamics of three different ligand-activated G-protein coupled receptors (GPCRs) and the photoreactive receptor rhodopsin from mammals were comparatively studied. As a result, diagrams demonstrating the main structural differences between the studied membrane receptors were obtained. These diagrams represent the projection of the crystal structures of rhodopsin photointermediates and ligand-activated receptors onto the plane defined by the principal components. Thus, we were able to associate the activation process of the receptors with large-scale movements of their individual transmembrane (TM) domains. In addition, the dynamics of extracellular loops of ligand-activated receptors responsible for recognition and initial binding of ligands was studied. Based on these results, two parameters of functionally significant structural dynamics of membrane receptors can be thoroughly analyzed simultaneously - movements of individual TM helices and of extracellular loops.

Biophysical fragment screening of the β1-adrenergic receptor: identification of high affinity arylpiperazine leads using structure-based drug design

Biophysical fragment screening of a thermostabilized β1-adrenergic receptor (β1AR) using surface plasmon resonance (SPR) enabled the identification of moderate affinity, high ligand efficiency (LE) arylpiperazine hits 7 and 8. Subsequent hit to lead follow-up confirmed the activity of the chemotype, and a structure-based design approach using protein-ligand crystal structures of the β1AR resulted in the identification of several fragments that bound with higher affinity, including indole 19 and quinoline 20. In the first example of GPCR crystallography with ligands derived from fragment screening, structures of the stabilized β1AR complexed with 19 and 20 were determined at resolutions of 2.8 and 2.7 Å, respectively.

X-ray structural information of GPCRs in drug design: what are the limitations and where do we go?

INTRODUCTION

In 2007, the X-ray structural determination of non-rhodopsin G-Protein coupled receptors (GPCRs), considered the most extensively targeted protein class for marketed drugs, commenced. With the relatively rapid availability of additional structures, an assessment of the progression made is needed in addition to the assessment of the understandings gleaned, deployment successes and forthcoming prospects.

AREAS COVERED

The author reviews the approaches and tools that have made it possible to determine the three dimensional structures of GPCRs using X-ray crystallography. Furthermore, the author describes the methods suited for crystallization of membrane bound GPCR proteins including the lipidic cubic phase and various protein modification approaches. The author also provides highlights, from the literature, of the structures determined to date including targets solved, the nature of the content provided (such as selectivity, activating vs. inactivating determinants) and how these structural features relate to drug design strategies.

EXPERT OPINION

The GPCR X-ray structures that have been so far determined have yielded significant information. This has presented dramatic evidence concerning their ability to impact the discovery of compounds through their action as traditional, orthosteric modulators. It is, however, noted that more challenging design strategies, such as identifying biased agonists and the use of sites remote from the orthosteric site for allosteric modulation, are still in their infancy.

Biomolecular simulation: a computational microscope for molecular biology

Molecular dynamics simulations capture the behavior of biological macromolecules in full atomic detail, but their computational demands, combined with the challenge of appropriately modeling the relevant physics, have historically restricted their length and accuracy. Dramatic recent improvements in achievable simulation speed and the underlying physical models have enabled atomic-level simulations on timescales as long as milliseconds that capture key biochemical processes such as protein folding, drug binding, membrane transport, and the conformational changes critical to protein function. Such simulation may serve as a computational microscope, revealing biomolecular mechanisms at spatial and temporal scales that are difficult to observe experimentally. We describe the rapidly evolving state of the art for atomic-level biomolecular simulation, illustrate the types of biological discoveries that can now be made through simulation, and discuss challenges motivating continued innovation in this field.

Full and partial agonists of thromboxane prostanoid receptor unveil fine tuning of receptor superactive conformation and G protein activation

The intrahelical salt bridge between E/D^3.49 and R^3.50 within the E/DRY motif on helix 3 (H3) and the interhelical hydrogen bonding between the E/DRY and residues on H6 are thought to be critical in stabilizing the class A G protein-coupled receptors in their inactive state. Removal of these interactions is expected to generate constitutively active receptors. This study examines how neutralization of E^3.49/6.30 in the thromboxane prostanoid (TP) receptor alters ligand binding, basal, and agonist-induced activity and investigates the molecular mechanisms of G protein activation. We demonstrate here that a panel of full and partial agonists showed an increase in affinity and potency for E129V and E240V mutants. Yet, even augmenting the sensitivity to detect constitutive activity (CA) with overexpression of the receptor or the G protein revealed resistance to an increase in basal activity, while retaining fully the ability to cause agonist-induced signaling. However, direct G protein activation measured through bioluminescence resonance energy transfer (BRET) indicates that these mutants more efficiently communicate and/or activate their cognate G proteins. These results suggest the existence of additional constrains governing the shift of TP receptor to its active state, together with an increase propensity of these mutants to agonist-induced signaling, corroborating their definition as superactive mutants. The particular nature of the TP receptor as somehow “resistant” to CA should be examined in the context of its pathophysiological role in the cardiovascular system. Evolutionary forces may have favored regulation mechanisms leading to low basal activity and selected against more highly active phenotypes.

Molecular signatures of G-protein-coupled receptors

G-protein-coupled receptors (GPCRs) are physiologically important membrane proteins that sense signalling molecules such as hormones and neurotransmitters, and are the targets of several prescribed drugs. Recent exciting developments are providing unprecedented insights into the structure and function of several medically important GPCRs. Here, through a systematic analysis of high-resolution GPCR structures, we uncover a conserved network of non-covalent contacts that defines the GPCR fold. Furthermore, our comparative analysis reveals characteristic features of ligand binding and conformational changes during receptor activation. A holistic understanding that integrates molecular and systems biology of GPCRs holds promise for new therapeutics and personalized medicine.

Endogenous lipid activated G protein-coupled receptors: emerging structural features from crystallography and molecular dynamics simulations

Class A G-protein coupled receptors (GPCRs) are thought to have a common topology that includes seven transmembrane alpha helices (TMHs) that are arranged to form a closed bundle. This bundle forms the ligand binding pocket into which ligands are commonly thought to enter via the extracellular milieu. This ligand approach direction makes sense for GPCRs that have small positively charged ligands, such as the beta-2-adrenergic or the dopamine D2 receptor. However, there is a growing sub-group of Class A GPCRs that bind lipid-derived endogenous ligands, such as the cannabinoid CB1 and CB2 receptors (with endogenous ligands, N-arachidonoylethanolamine (anandamide) and sn-2-arachidonylglycerol (2-AG)) and the S1P1-5 receptors (with endogenous ligand, sphingosine-1-phosphate). Even the widely studied Class A GPCR, rhodopsin, binds a highly lipophillic chromophore (11-cis-retinal). For these receptors, ligand approach from the extracellular milieu has seemed unlikely given that the ligands of these receptors readily partition into lipid or are actually synthesized in the lipid bilayer. The recent X-ray-crystal structure of the sub-type 1 sphingosine-1-phosphate receptor (S1P1) provides important information on the key structural variations that may be the hallmarks for a Class A GPCR that binds lipid-derived ligands. These include an extracellular domain that is closed off to the extracellular milieu and the existence of an opening between transmembrane helices that may serve as a portal for ligand entry via the lipid bilayer. This review examines structural aspects that the cannabinoid receptors may share with the S1P1 receptor based upon sequence homology. This review also examines experimental and simulation results that suggest ligand entry via a lipid portal is quite likely for this emerging sub-group.

Crystal structure of oligomeric β1-adrenergic G protein-coupled receptors in ligand-free basal state

G protein-coupled receptors (GPCRs) mediate transmembrane signaling. Before ligand binding, GPCRs exist in a basal state. Crystal structures of several GPCRs bound with antagonists or agonists have been solved. However, the crystal structure of the ligand-free basal state of a GPCR, the starting point of GPCR activation and function, has not been determined. Here we report the X-ray crystal structure of the first ligand-free basal state of a GPCR in a lipid membrane-like environment. Oligomeric turkey β₁-adrenergic receptors display two alternating dimer interfaces. One interface involves the transmembrane domain (TM) 1, TM2, the C-terminal H8, and the extracellular loop 1. The other interface engages residues from TM4, TM5, the intracellular loop 2 and the extracellular loop 2. Structural comparisons show that this ligand-free state is in an inactive conformation. This provides the structural information regarding GPCR dimerization and oligomerization.

Structure of β-adrenergic receptors

β-Adrenergic receptors (βARs) control key physiological functions by transducing signals encoded in catecholamine hormones and neurotransmitters to activate intracellular signaling pathways. As members of the large family of G protein-coupled receptors (GPCRs), βARs have a seven-transmembrane helix topology and signal via G protein- and arrestin-dependent pathways. Until 2007, three-dimensional structural information of GPCRs activated by diffusible ligands, including βARs, was limited to homology models that used the related photoreceptor rhodopsin as a template. Over many years, several labs have developed strategies that have finally allowed the structures of the turkey β(1)AR and the human β(2)AR to be determined experimentally. The challenges to overcome included heterologous receptor overexpression, design of stabilized and crystallizable modified receptor constructs, ligand-affinity purification of active receptor and the development of novel techniques in crystallization and microcrystallography. The structures of βARs in complex with inverse agonists, antagonists, and agonists have revealed the binding mode of ligands with different efficacies, have allowed to obtain insights into ligand selectivity, and have provided better templates for drug design. Also, the structures of β(2)AR in complex with a G protein and a G protein-mimicking nanobody have provided important insights into the mechanism of receptor activation and G protein coupling. This chapter summarizes the strategies and methods that have been successfully applied to the structural studies of βARs. These are exemplified with detailed protocols toward the structure determination of stabilized turkey β(1)AR-ligand complexes. We also discuss the spectacular insights into adrenergic receptor function that were obtained from the structures.

Structural determinants of arrestin functions

Arrestins are a small protein family with only four members in mammals. Arrestins demonstrate an amazing versatility, interacting with hundreds of different G protein-coupled receptor (GPCR) subtypes, numerous nonreceptor signaling proteins, and components of the internalization machinery, as well as cytoskeletal elements, including regular microtubules and centrosomes. Here, we focus on the structural determinants that mediate various arrestin functions. The receptor-binding elements in arrestins were mapped fairly comprehensively, which set the stage for the construction of mutants targeting particular GPCRs. The elements engaged by other binding partners are only now being elucidated and in most cases we have more questions than answers. Interestingly, even very limited and imprecise identification of structural requirements for the interaction with very few other proteins has enabled the development of signaling-biased arrestin mutants. More comprehensive understanding of the structural underpinning of different arrestin functions will pave the way for the construction of arrestins that can link the receptor we want to the signaling pathway of our choosing.

GPCR activation: a mutagenic spotlight on crystal structures

The crystal structures of antagonist and agonist complexes of isolated β(2) and β(1) adrenoceptors have recently been supplemented by antagonist structures of M(2) and M(3) muscarinic acetylcholine receptors. Importantly, a structure of an agonist-ligated β(2) adrenoceptor complexed with its cognate G protein has provided the first view of a ternary complex representing the transition state in agonist-mediated G protein activation. This review interprets these G-protein-coupled receptor (GPCR) structures through the focus provided by extensive mutagenesis studies on muscarinic receptors, revealing an activation mechanism that is both modular and dynamic. Specific motifs, based around highly conserved residues, functionalise the seven-transmembrane architecture of these receptors. While exploiting conserved motifs, the ligand binding and signal transduction pathways work around and through water-containing cavities, an emerging feature of GPCR structures. These cavities may have undergone evolutionary selection to adapt GPCRs to particular signalling niches, and may provide targeting opportunities to enhance drug selectivity.

Restructuring G-protein- coupled receptor activation

G-protein-coupled receptors serve as key signal transduction conduits, linking extracellular inputs with diverse cellular responses. These receptors eluded structural characterization for decades following their identification. A landmark structure of rhodopsin provided a basis for structure-function studies and homology modeling, but advances in receptor biology suffered from a lack of receptor-specific structural insights. The recent explosion in GPCR structures confirms some features predicted by rhodopsin-based models, and more importantly, it reveals unexpected ligand-binding modes and critical aspects of the receptor activation process. The new structures also promise to foster studies testing emerging models for GPCR function such as receptor dimerization and ligand-biased signaling.

Synthesis and structure-affinity relationships of selective high-affinity 5-HT(4) receptor antagonists: application to the design of new potential single photon emission computed tomography tracers

The work described herein aims at finding new potential ligands for the brain imaging of 5-HT(4) receptors (5-HT(4)Rs) using single-photon emission computed tomography (SPECT). Starting from the nonsubstituted phenanthridine compound 4a, exhibiting a K(i) value of 51 nM on the 5-HT(4)R, we explored the structure-affinity in this series. We found that substitution in position 4 of the tricycle with a fluorine atom gave the best result. Introduction of an additional nitrogen atom inside the tricyclic framework led to an increase of both the affinity and selectivity for 5-HT(4)R, suggesting the design of the antagonist 4v, exhibiting a high affinity of 0.04 nM. Several iodinated analogues were then synthesized as potential SPECT tracers. The iodinated compound 11d was able to displace the reference radioiodinated 5-HT(4)R antagonist (1-butylpiperidin-4-yl)methyl-8-amino-7-iodo[(123)I]-2,3-dihydrobenzo[b][1,4]dioxine-5-carboxylate {[(123)I]1, [(123)I]SB 207710} both in vitro and in vivo in brain. Compound 11d was radiolabeled with [(125)I]iodine, providing a potential SPECT candidate for brain imaging of 5-HT(4)R.

Structure-function of the G protein-coupled receptor superfamily

During the past few years, crystallography of G protein-coupled receptors (GPCRs) has experienced exponential growth, resulting in the determination of the structures of 16 distinct receptors-9 of them in 2012 alone. Including closely related subtype homology models, this coverage amounts to approximately 12% of the human GPCR superfamily. The adrenergic, rhodopsin, and adenosine receptor systems are also described by agonist-bound active-state structures, including a structure of the receptor-G protein complex for the β(2)-adrenergic receptor. Biochemical and biophysical techniques, such as nuclear magnetic resonance and hydrogen-deuterium exchange coupled with mass spectrometry, are providing complementary insights into ligand-dependent dynamic equilibrium between different functional states. Additional details revealed by high-resolution structures illustrate the receptors as allosteric machines that are controlled not only by ligands but also by ions, lipids, cholesterol, and water. This wealth of data is helping redefine our knowledge of how GPCRs recognize such a diverse array of ligands and how they transmit signals 30 angstroms across the cell membrane; it also is shedding light on a structural basis of GPCR allosteric modulation and biased signaling.

[G protein-coupled receptors in the spotlight]

G protein-coupled receptors (GPCR) are the largest family of integral membrane proteins found in the plasma membrane of mammalian cells. GPCR respond to a large variety of ligands such as amines, lipids, hormones and amino-acids, which are involved in inter-cellular signalling events in a multitude of physiological and pathological processes. GPCRs are therefore key regulators of signal transduction by which cells respond to variations in their environment. During the last five years, striking progress has been made to solve high-resolution structure of GPCR. The most recent successes are the structures of the β(1) and β(2) adrenoreceptors and the adenosine A(2A) receptor bound to a variety of agonists. Most importantly, the structure of the β(2) adrenoreceptor in complex with a trimeric G protein, Gs, was recently reported. This review will present an overview of the X-ray structure determination of the GPCR and of their activation mechanism.

β-Adrenergic agonists mediate enhancement of β1-adrenergic receptor N-terminal cleavage and stabilization in vivo and in vitro

The β(1)-adrenergic receptor (β(1)AR) is the predominant βAR in the heart and is the main target for β-adrenergic antagonists, widely used in the treatment of cardiovascular diseases. Previously, we have shown that the human (h) β(1)AR is cleaved in its N terminus by a metalloproteinase, both constitutively and in a receptor activation-dependent manner. In this study, we investigated the specific events involved in β(1)AR regulation, focusing on the effects of long-term treatment with β-adrenergic ligands on receptor processing in stably transfected human embryonic kidney 293(i) cells. The key findings were verified using the transiently transfected hβ(1)AR and the endogenously expressed receptor in neonatal rat cardiomyocytes. By using flow cytometry and Western blotting, we demonstrated that isoproterenol, S-propranolol, CGP-12177 [4-[3-[(1,1-dimethylethyl)amino]2-hydroxypropoxy]-1,3-dihydro-2H-benzimidazol-2-one], pindolol, and timolol, which displayed agonistic properties toward the β(1)AR in either the adenylyl cyclase or the mitogen-activated protein kinase signaling pathways, induced cleavage of the mature cell-surface receptor. In contrast, metoprolol, bisoprolol, and CGP-20712 [1-[2-((3-carbamoyl-4-hydroxy)phenoxy)ethylamino]-3-[4-(1-methyl-4-trifluoromethyl-2-imidazolyl)phenoxy]-2-propanol], which showed no agonistic activity, had only a marginal or no effect. Importantly, the agonists also stabilized intracellular receptor precursors, possibly via their pharmacological chaperone action, and they stabilized the receptor in vitro. The opposing effects on the two receptor forms thus led to an increase in the amount of cleaved receptor fragments at the plasma membrane. The results underscore the pluridimensionality of β-adrenergic ligands and extend this property from receptor activation and signaling to the regulation of β(1)AR levels. This phenomenon may contribute to the exceptional resistance of β(1)ARs to downregulation and tendency toward upregulation following long-term ligand treatments.

Crystal structures of a stabilized β1-adrenoceptor bound to the biased agonists bucindolol and carvedilol

The β₁-adrenoceptor (β₁AR) is the site of action of beta blockers used in the treatment of cardiac-related illnesses. Two beta blockers, carvedilol and bucindolol, show distinctive activities compared to other beta blockers and have been proposed as treatments tailored to the Arg/Gly389^8.56 polymorphism of the human β₁AR. Both carvedilol and bucindolol are classified as biased agonists, because they stimulate G protein-independent signaling, while acting as either inverse or partial agonists of the G protein pathway. We have determined the crystal structures of a thermostabilized avian β₁AR mutant bound to bucindolol and to carvedilol at 3.2 and 2.3 Å resolution, respectively. In comparison to other beta blockers, bucindolol and carvedilol interact with additional residues, in extracellular loop 2 and transmembrane helix 7, which may promote G protein-independent signaling. The structures also suggest that there may be a structural explanation for the pharmacological differences arising from the Arg/Gly389^8.56 polymorphism.

Action of molecular switches in GPCRs--theoretical and experimental studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

Modeling of mammalian olfactory receptors and docking of odorants

Olfactory receptors (ORs) belong to the superfamily of G protein-coupled receptors (GPCRs), the second largest class of genes after those related to immunity, and account for about 3 % of mammalian genomes. ORs are present in all multicellular organisms and represent more than half the GPCRs in mammalian species (e.g., the mouse OR repertoire contains >1,000 functional genes). ORs are mainly expressed in the olfactory epithelium where they detect odorant molecules, but they are also expressed in a number of other cells, such as sperm cells, although their functions in these cells remain mostly unknown. It has recently been reported that ORs are present in tumoral tissues where they are expressed at different levels than in healthy tissues. A specific OR is over-expressed in prostate cancer cells, and activation of this OR has been shown to inhibit the proliferation of these cells. Odorant stimulation of some of these receptors results in inhibition of cell proliferation. Even though their biological role has not yet been elucidated, these receptors might constitute new targets for diagnosis and therapeutics. It is important to understand the activation mechanism of these receptors at the molecular level, in particular to be able to predict which ligands are likely to activate a particular receptor ('deorphanization') or to design antagonists for a given receptor. In this review, we describe the in silico methodologies used to model the three-dimensional (3D) structure of ORs (in the more general framework of GPCR modeling) and to dock ligands into these 3D structures.

Leu128(3.43) (l128) and Val247(6.40) (V247) of CXCR1 are critical amino acid residues for g protein coupling and receptor activation

CXCR1, a classic GPCR that binds IL-8, plays a key role in neutrophil activation and migration by activating phospholipase C (PLC)β through Gα₁₅ and Gα_i which generates diacylglycerol and inositol phosphates (IPs). In this study, two conserved amino acid residues of CXCR1 on the transmembrane domain (TM) 3 and TM6, Leu128^3.43 (L128) and Val247^6.40 (V247), respectively, were selectively substituted with other amino acids to investigate the role of these conserved residues in CXCR1 activation. Although two selective mutants on Leu128, Leu128Ala (L128A) and Leu128Arg (L128R), demonstrated high binding affinity to IL-8, they were not capable of coupling to G proteins and consequently lost the functional response of the receptors. By contrast, among the four mutants at residue Val247 (TM6.40), replacing Val247 with Ala (V247A) and Asn (V247N) led to constitutive activation of mutant receptors when cotransfected with Gα₁₅. The V247N mutant also constitutively activated the Gα_i protein. These results indicate that L128 on TM3.43 is involved in G protein coupling and receptor activation but is unimportant for ligand binding. On the other hand, V247 on TM6.40 plays a critical role in maintaining the receptor in the inactive state, and the substitution of V247 impaired the receptor constraint and stabilized an active conformation. Functionally, there was an increase in chemotaxis in response to IL-8 in cells expressing V247A and V247N. Our findings indicate that Leu128^3.43 and Val247^6.40 are critical for G protein coupling and activation of signaling effectors, providing a valuable insight into the mechanism of CXCR1 activation.

Induction of cardiac fibrosis by β-blocker in G protein-independent and G protein-coupled receptor kinase 5/β-arrestin2-dependent Signaling pathways

G-protein coupled receptors (GPCRs) have long been known as receptors that activate G protein-dependent cellular signaling pathways. In addition to the G protein-dependent pathways, recent reports have revealed that several ligands called "biased ligands" elicit G protein-independent and β-arrestin-dependent signaling through GPCRs (biased agonism). Several β-blockers are known as biased ligands. All β-blockers inhibit the binding of agonists to the β-adrenergic receptors. In addition to β-blocking action, some β-blockers are reported to induce cellular responses through G protein-independent and β-arrestin-dependent signaling pathways. However, the physiological significance induced by the β-arrestin-dependent pathway remains much to be clarified in vivo. Here, we demonstrate that metoprolol, a β(1)-adrenergic receptor-selective blocker, could induce cardiac fibrosis through a G protein-independent and β-arrestin2-dependent pathway. Metoprolol, a β-blocker, increased the expression of fibrotic genes responsible for cardiac fibrosis in cardiomyocytes. Furthermore, metoprolol induced the interaction between β(1)-adrenergic receptor and β-arrestin2, but not β-arrestin1. The interaction between β(1)-adrenergic receptor and β-arrestin2 by metoprolol was impaired in the G protein-coupled receptor kinase 5 (GRK5)-knockdown cells. Metoprolol-induced cardiac fibrosis led to cardiac dysfunction. However, the metoprolol-induced fibrosis and cardiac dysfunction were not evoked in β-arrestin2- or GRK5-knock-out mice. Thus, metoprolol is a biased ligand that selectively activates a G protein-independent and GRK5/β-arrestin2-dependent pathway, and induces cardiac fibrosis. This study demonstrates the physiological importance of biased agonism, and suggests that G protein-independent and β-arrestin-dependent signaling is a reason for the diversity of the effectiveness of β-blockers.

A crystal clear solution for determining G-protein-coupled receptor structures

G-protein-coupled receptors (GPCRs) are medically important membrane proteins that are targeted by over 30% of small molecule drugs. At the time of writing, 15 unique GPCR structures have been determined, with 77 structures deposited in the PDB database, which offers new opportunities for drug development and for understanding the molecular mechanisms of GPCR activation. Many different factors have contributed to this success, but if there is one single factor that can be singled out as the foundation for producing well-diffracting GPCR crystals, it is the stabilisation of the detergent-solubilised receptor-ligand complex. This review will focus predominantly on one of the successful strategies for the stabilisation of GPCRs, namely the thermostabilisation of GPCRs using systematic mutagenesis coupled with thermostability assays. Structures of thermostabilised GPCRs bound to a wide variety of ligands have been determined, which has led to an understanding of ligand specificity; why some ligands act as agonists as opposed to partial or inverse agonists; and the structural basis for receptor activation.

Lifting the lid on GPCRs: the role of extracellular loops

GPCRs exhibit a common architecture of seven transmembrane helices (TMs) linked by intracellular loops and extracellular loops (ECLs). Given their peripheral location to the site of G-protein interaction, it might be assumed that ECL segments merely link the important TMs within the helical bundle of the receptor. However, compelling evidence has emerged in recent years revealing a critical role for ECLs in many fundamental aspects of GPCR function, which supported by recent GPCR crystal structures has provided mechanistic insights. This review will present current understanding of the key roles of ECLs in ligand binding, activation and regulation of both family A and family B GPCRs.

Small molecules with similar structures exhibit agonist, neutral antagonist or inverse agonist activity toward angiotensin II type 1 receptor

Small differences in the chemical structures of ligands can be responsible for agonism, neutral antagonism or inverse agonism toward a G-protein-coupled receptor (GPCR). Although each ligand may stabilize the receptor conformation in a different way, little is known about the precise conformational differences. We synthesized the angiotensin II type 1 receptor blocker (ARB) olmesartan, R239470 and R794847, which induced inverse agonism, antagonism and agonism, respectively, and then investigated the ligand-specific changes in the receptor conformation with respect to stabilization around transmembrane (TM)3. The results of substituted cysteine accessibility mapping studies support the novel concept that ligand-induced changes in the conformation of TM3 play a role in stabilizing GPCR. Although the agonist-, neutral antagonist and inverse agonist-binding sites in the AT(1) receptor are similar, each ligand induced specific conformational changes in TM3. In addition, all of the experimental data were obtained with functional receptors in a native membrane environment (in situ).

New insights for drug design from the X-ray crystallographic structures of G-protein-coupled receptors

Methodological advances in X-ray crystallography have made possible the recent solution of X-ray structures of pharmaceutically important G protein-coupled receptors (GPCRs), including receptors for biogenic amines, peptides, a nucleoside, and a sphingolipid. These high-resolution structures have greatly increased our understanding of ligand recognition and receptor activation. Conformational changes associated with activation common to several receptors entail outward movements of the intracellular side of transmembrane helix 6 (TM6) and movements of TM5 toward TM6. Movements associated with specific agonists or receptors have also been described [e.g., extracellular loop (EL) 3 in the A(2A) adenosine receptor]. The binding sites of different receptors partly overlap but differ significantly in ligand orientation, depth, and breadth of contact areas in TM regions and the involvement of the ELs. A current challenge is how to use this structural information for the rational design of novel potent and selective ligands. For example, new chemotypes were discovered as antagonists of various GPCRs by subjecting chemical libraries to in silico docking in the X-ray structures. The vast majority of GPCR structures and their ligand complexes are still unsolved, and no structures are known outside of family A GPCRs. Molecular modeling, informed by supporting information from site-directed mutagenesis and structure-activity relationships, has been validated as a useful tool to extend structural insights to related GPCRs and to analyze docking of other ligands in already crystallized GPCRs.

Structural transitions of transmembrane helix 6 in the formation of metarhodopsin I

Absorption of light by the visual pigment rhodopsin triggers a rapid cis-trans photoisomerization of its retinal chromophore and a series of conformational changes in both the retinal and protein. The largest structural change is an outward tilt of transmembrane helix H6 that increases the separation of the intracellular ends of H6 and H3 and opens up the G-protein binding site. In the dark state of rhodopsin, Glu247 at the intracellular end of H6 forms a salt bridge with Arg135 on H3 to tether H6 in an inactive conformation. The Arg135-Glu247 interaction is broken in the active state of the receptor, and Arg135 is then stabilized by interactions with Tyr223, Met257, and Tyr306 on helices H5, H6, and H7, respectively. To address the mechanism of H6 motion, solid-state NMR measurements are undertaken of Metarhodopsin I (Meta I), the intermediate preceding the active Metarhodopsin II (Meta II) state of the receptor. (13)C NMR dipolar recoupling measurements reveal an interhelical contact of (13)Cζ-Arg135 with (13)Cε-Met257 in Meta I but not with (13)Cζ-Tyr223 or (13)Cζ-Tyr306. These observations suggest that helix H6 has rotated in the formation of Meta I but that structural changes involving helices H5 and H7 have not yet occurred. Together, our results provide insights into the sequence of events leading up to the outward motion of H6, a hallmark of G protein-coupled receptor activation.

Comparison of dynamics of extracellular accesses to the β(1) and β(2) adrenoceptors binding sites uncovers the potential of kinetic basis of antagonist selectivity

From the molecular mechanism of antagonist unbinding in the β(1) and β(2) adrenoceptors investigated by steered molecular dynamics, we attempt to provide further possibilities of ligand subtype and subspecies selectivity. We have simulated unbinding of β(1)-selective Esmolol and β(2)-selective ICI-118551 from both receptors to the extracellular environment and found distinct molecular features of unbinding. By calculating work profiles, we show different preference in antagonist unbinding pathways between the receptors, in particular, perpendicular to the membrane pathway is favourable in the β(1) adrenoceptor, whereas the lateral pathway involving helices 5, 6 and 7 is preferable in the β(2) adrenoceptor. The estimated free energy change of unbinding based on the preferable pathway correlates with the experimental ligand selectivity. We then show that the non-conserved K347 (6.58) appears to facilitate in guiding Esmolol to the extracellular surface via hydrogen bonds in the β(1) adrenoceptor. In contrast, hydrophobic and aromatic interactions dominate in driving ICI-118551 through the easiest pathway in the β(2) adrenoceptor. We show how our study can stimulate design of selective antagonists and discuss other possible molecular reasons of ligand selectivity, involving sequential binding of agonists and glycosylation of the receptor extracellular surface.

Comparative analysis of the heptahelical transmembrane bundles of G protein-coupled receptors

Background

G protein-coupled receptors represent a large family of eukaryotic membrane proteins, and are involved in almost all physiological processes in humans. Recent advances in the crystallographic study of these receptors enable a detailed comparative analysis of the commonly shared heptahelical transmembrane bundle. Systematic comparison of the bundles from a variety of receptors is indispensable for understanding not only of the structural diversification optimized for the binding of respective ligands but also of the structural conservation required for the common mechanism of activation accompanying the interaction changes among the seven helices.

Methodology/Principal Findings

We have examined the bundles of 94 polypeptide chains from almost all available structures of 11 receptors, which we classified into either inactivated chain or activated chain, based on the type of bound ligand. For the inactivated chains, superposition of 200 residue bundles by secondary structure matching demonstrated that the bound ligands share a laterally limited cavity in the extracellular section of the bundle. Furthermore, a distinct feature was found for helix III of bovine rhodopsin, which might have evolved to lower its activity in the presence of 11-cis-retinal, to a level that other receptors could hardly achieve with any currently available ligands.

Conclusions/Significance

Systematic analysis described here would be valuable for understanding of the rearrangement of seven helices which depends on the ligand specificity and activation state of the receptors.

Ligand-dependent conformations and dynamics of the serotonin 5-HT(2A) receptor determine its activation and membrane-driven oligomerization properties

From computational simulations of a serotonin 2A receptor (5-HT_2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT_2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization.

The 5-HT_2A receptor for the neurotransmitter serotonin (5-HT) belongs to family A (rhodopsin-like) G-protein coupled receptors (GPCRs), one of the most important classes of membrane proteins that are targeted by an extensive and diverse collection of external stimuli. Recently we learned that different ligands targeting the same GPCR can elicit different biological responses, but the mechanisms remain unknown. We address this fundamental question for the serotonin 5-HT_2A receptor, because it is known to respond to the binding of structurally diverse ligands by producing similar stimuli in the cell, and to the binding of quite similar ligands with dramatically different responses. Molecular dynamics simulations of molecular models of the serotonin 5-HT_2A receptor in complex with pharmacologically distinct ligands show the dynamic rearrangements of the receptor molecule to be different for these ligands, and the nature and extents of the rearrangements reflect the known pharmacological properties of the ligands as full, partial or inverse activators of the receptor. The different rearrangements of the receptor molecule are shown to produce different rearrangements of the surrounding membrane, a remodeling of the environment that can have differential ligand-determined effects on receptor function and association in the cell's membrane.

Agonist-bound structures of G protein-coupled receptors

G protein-coupled receptors (GPCRs) play a major role in intercellular communication by binding small diffusible ligands (agonists) at the extracellular surface. Agonist-binding induces a conformational change in the receptor, which results in the binding and activation of heterotrimeric G proteins within the cell. Ten agonist-bound structures of non-rhodopsin GPCRs published last year defined for the first time the molecular details of receptor activated states and how inverse agonists, partial agonists and full agonists bind to produce different effects on the receptor. In addition, the structure of the β(2)-adrenoceptor coupled to a heterotrimeric G protein showed how the opening of a cleft in the cytoplasmic face of the receptor as a consequence of agonist binding results in G protein coupling and activation of the G protein.

Revisiting the homology modeling of G-protein coupled receptors: β1-adrenoceptor as an example

G-protein coupled receptors (GPCRs) are recognized to constitute the largest family of membrane proteins. Due to the disproportion in the quantity of crystal structures and their amino acid sequences, homology modeling contributes a reasonable and feasible approach to GPCR theoretical coordinates. With the brand new crystal structures resolved recently, herein we deliberated how to designate them as templates to carry out homology modeling in four aspects: (1) various sequence alignment methods; (2) protein weight matrix; (3) different sets of multiple templates; (4) active and inactive state of templates. The accuracy of models was evaluated by comparing the similarity of stereo conformation and molecular docking results between models and the experimental structure of Meleagris gallopavo β(1)-adrenergic receptor (Mg_Adrb1) that we desired to develop as an example. Our results proposed that: (1) Cobalt and MAFFT, two algorithms of sequence alignment, were suitable for single- and multiple-template modeling, respectively; (2) Blosum30 is applicable to align sequences in the case of low sequence identity; (3) multiple-template modeling is not always better than single-template one; (4) the state of template is an influential factor in simulating the GPCR structures as well.

Fluorescence/bioluminescence resonance energy transfer techniques to study G-protein-coupled receptor activation and signaling

Fluorescence and bioluminescence resonance energy transfer (FRET and BRET) techniques allow the sensitive monitoring of distances between two labels at the nanometer scale. Depending on the placement of the labels, this permits the analysis of conformational changes within a single protein (for example of a receptor) or the monitoring of protein-protein interactions (for example, between receptors and G-protein subunits). Over the past decade, numerous such techniques have been developed to monitor the activation and signaling of G-protein-coupled receptors (GPCRs) in both the purified, reconstituted state and in intact cells. These techniques span the entire spectrum from ligand binding to the receptors down to intracellular second messengers. They allow the determination and the visualization of signaling processes with high temporal and spatial resolution. With these techniques, it has been demonstrated that GPCR signals may show spatial and temporal patterning. In particular, evidence has been provided for spatial compartmentalization of GPCRs and their signals in intact cells and for distinct physiological consequences of such spatial patterning. We review here the FRET and BRET technologies that have been developed for G-protein-coupled receptors and their signaling proteins (G-proteins, effectors) and the concepts that result from such experiments.

Cell-based and in-silico studies on the high intrinsic activity of two boron-containing salbutamol derivatives at the human β₂-adrenoceptor

Salbutamol is a well-known β(2) adrenoceptor (β(2)AR) partial agonist. We synthesized two boron-containing salbutamol derivatives (BCSDs) with greater potency and efficacy, compared to salbutamol, for inducing β(2)AR-mediated smooth-muscle relaxation in guinea-pig tracheal rings. However, the mechanism involved in this pharmacological effect remains unclear. In order to gain insight, we carried out binding and functional assays for BCSDs in HEK-293T cells transfected with the human β(2)AR (hβ(2)AR). The transfected hβ(2)AR showed similar affinity for BCSDs and salbutamol, but adenosine 3',5'-cyclic phosphate (cAMP) accumulation induced by both BCSDs was similar to that elicited by isoproterenol and greater than that induced by salbutamol. The boron-containing precursors (boric and phenylboronic acids, 100 μM) had no significant effect on salbutamol binding or salbutamol-induced cAMP accumulation. These experimental results are in agreement with theoretical docking simulations on lipid bilayer membrane-embedded hβ(2)AR structures. These receptors showed slightly higher affinity for BCSDs than for salbutamol. An essential change between putative active and inactive conformational states depended on the interaction of the tested ligands with the fifth, sixth and seventh transmembrane domains. Overall, these data suggest that BCSDs induce and stabilize conformational states of the hβ(2)AR that are highly capable of stimulating cAMP production.

ER-bound steps in the biosynthesis of G protein-coupled receptors

The polypeptide of a G protein-coupled receptor is inserted into the membrane of the endoplasmic reticulum while being translated and this process by itself may be sufficient to establish the proper receptor fold. X-ray structures reveal a common polypeptide topology with little variation in the alignment and orientation of the seven transmembrane segments, the proximal carboxyl terminus (C-tail) and parts of the extracellular loops. These define a structural core the stability of which probably represents a major criterion for the receptor to pass endoplasmic reticulum (ER) quality control; point mutations affecting the structure of the core have an extraordinary chance of causing receptor retention. In contrast, cytoplasmic loops 2 and 3 and the distal C-tail are poorly ordered at least in the absence of an interaction partner. Similarly, the amino terminal tail of rhodopsin-related receptors (but not of receptor subtypes where ligand binding requires a stable fold of the N-tail) is unlikely to establish a stable fold. These segments can cause ER retention when mutated to inappropriately expose hydrophobic peptide patches; to prevent protein aggregation chaperone molecules attach to them thus initiating selection for ER-associated degradation. It is less clear however if there are additional mechanisms to specifically survey the transmembrane core at the level of the lipid bilayer or if insufficient packing is detected due to misalignment of the cytoplasmic or extracellular face of the receptor.

Ligand and decoy sets for docking to G protein-coupled receptors

We compiled a G protein-coupled receptor (GPCR) ligand library (GLL) for 147 targets, selecting for each ligand 39 decoy molecules, collected in the GPCR Decoy Database (GDD). Decoys were chosen ensuring a ligand-decoy similarity of six physical properties, while enforcing ligand-decoy chemical dissimilarity. The performance in docking of the GDD was evaluated on 19 GPCRs, showing a marked decrease in enrichment compared to bias-uncorrected decoy sets. Both the GLL and GDD are freely available for the scientific community.

Diversity and modularity of G protein-coupled receptor structures

G protein-coupled receptors (GPCRs) comprise the most 'prolific' family of cell membrane proteins, which share a general mechanism of signal transduction, but greatly vary in ligand recognition and function. Crystal structures are now available for rhodopsin, adrenergic, and adenosine receptors in both inactive and activated forms, as well as for chemokine, dopamine, and histamine receptors in inactive conformations. Here we review common structural features, outline the scope of structural diversity of GPCRs at different levels of homology, and briefly discuss the impact of the structures on drug discovery. Given the current set of GPCR crystal structures, a distinct modularity is now being observed between the extracellular (ligand-binding) and intracellular (signaling) regions. The rapidly expanding repertoire of GPCR structures provides a solid framework for experimental and molecular modeling studies, and helps to chart a roadmap for comprehensive structural coverage of the whole superfamily and an understanding of GPCR biological and therapeutic mechanisms.

Engineering an ultra-thermostable β(1)-adrenoceptor

Conformational thermostabilisation of G-protein-coupled receptors is a successful strategy for their structure determination. The thermostable mutants tolerate short-chain detergents, such as octylglucoside and nonylglucoside, which are ideal for crystallography, and in addition, the receptors are preferentially in a single conformational state. The first thermostabilised receptor to have its structure determined was the β(1)-adrenoceptor mutant β(1)AR-m23 bound to the antagonist cyanopindolol, and recently, additional structures have been determined with agonist bound. Here, we describe further stabilisation of β(1)AR-m23 by the addition of three thermostabilising mutations (I129V, D322K, and Y343L) to make a mutant receptor that is 31 °C more thermostable than the wild-type receptor in dodecylmaltoside and is 13 °C more thermostable than β(1)AR-m23 in nonylglucoside. Although a number of thermostabilisation methods were tried, including rational design of disulfide bonds and engineered zinc bridges, the two most successful strategies to improve the thermostability of β(1)AR-m23 were an engineered salt bridge and leucine scanning mutagenesis. The three additional thermostabilising mutations did not significantly affect the pharmacological properties of β(1)AR-m23, but the new mutant receptor was significantly more stable in short-chain detergents such as heptylthioglucoside and denaturing detergents such as SDS.

Structural insights into agonist-induced activation of G-protein-coupled receptors

Recent years have seen tremendous breakthroughs in structure determination of G-protein-coupled receptors (GPCRs). In 2011, two agonist-bound active-state structures of rhodopsin have been published. Together with structures of several rhodopsin activation intermediates and a wealth of biochemical and spectroscopic information, they provide a unique structural framework on which to understand GPCR activation. Here we use this framework to compare the recent crystal structures of the agonist-bound active states of the β(2) adrenergic receptor (β(2)AR) and the A(2A) adenosine receptor (A(2A)AR). While activation of these three GPCRs results in rearrangements of TM5 and TM6, the extent of this conformational change varies considerably. Displacements of the cytoplasmic side of TM6 ranges between 3 and 8Å depending on whether selective stabilizers of the active conformation are used (i.e. a G-protein peptide in the case of rhodopsin or a conformationally selective nanobody in the case of the β(2)AR) or not (A(2A)AR). The agonist-induced conformational changes in the ligand-binding pocket are largely receptor specific due to the different chemical nature of the agonists. However, several similarities can be observed, including a relocation of conserved residues W6.48 and F6.44 towards L5.51 and P5.50, and of I/L3.40 away from P5.50. This transmission switch links agonist binding to the movement of TM5 and TM6 through the rearrangement of the TM3-TM5-TM6 interface, and possibly constitutes a common theme of GPCR activation.