Ligand-specific receptor states: implications for opiate receptor signalling and regulation

Opiate drugs produce their effects by acting upon G protein coupled receptors (GPCRs) and although they are among the most effective analgesics available, their clinical use is restricted by unwanted side effects such as tolerance, physical dependence, respiratory depression, nausea and constipation. As a class, opiates share a common profile of unwanted effects but there are also significant differences in ligand liability for producing these actions. A growing number of studies show that GPCRs may exist in multiple active states that differ in their signalling and regulatory properties and which may distinctively bind different agonists. In this review we summarize evidence supporting the existence of multiple active conformations for MORs and DORs, analyze information favouring the existence of ligand-specific receptor states and assess how ligand-selective efficacy may contribute to the production of longer lasting, better tolerated opiate analgesics.     

Nanobody stabilization of G protein-coupled receptor conformational states

Remarkable progress has been made in the field of G protein-coupled receptor (GPCR) structural biology during the past four years. Several obstacles to generating diffraction quality crystals of GPCRs have been overcome by combining innovative methods ranging from protein engineering to lipid-based screens and microdiffraction technology. The initial GPCR structures represent energetically stable inactive-state conformations. However, GPCRs signal through different G protein isoforms or G protein-independent effectors upon ligand binding suggesting the existence of multiple ligand-specific active states. These active-state conformations are unstable in the absence of specific cytosolic signaling partners representing new challenges for structural biology. Camelid single chain antibody fragments (nanobodies) show promise for stabilizing active GPCR conformations and as chaperones for crystallogenesis.     

Molecular modelling of human 5-hydroxytryptamine receptor (5-HT2A) and virtual screening studies towards the identification of agonist and antagonist molecules

The serotonin receptors, also known as 5-hydroxytryptamine (5-HT) receptors, are a group of G protein-coupled receptors (GPCRs) and ligand-gated ion channels found in the central and peripheral nervous systems. GPCRs have a characteristic feature of activating different signalling pathways upon ligand binding and these ligands display several efficacy levels to differentially activate the receptor. GPCRs are primary drug targets due to their central role in several signal transduction pathways. Drug design for GPCRs is also most challenging due to their inherent promiscuity in ligand recognition, which gives rise to several side effects of existing drugs. Here, we have performed the ligand interaction study using the two prominent states of GPCR, namely the active and inactive state of the 5-HT_2A receptor. Active state of 5-HT_2A receptor model enhances the understanding of conformational difference which influences the ligand-binding site. A 5-HT_2A receptor active state model was constructed by homology modelling using active state β₂-adrenergic receptor (β₂-AR). In addition, virtual screening and docking studies with both active and inactive state models reveal potential small molecule hits which could be considered as agonist-like and antagonist-like molecules. The results from the all-atom molecular dynamics simulations further confirmed that agonists and antagonists interact in different modes with the receptor.     

Toward the active conformations of rhodopsin and the beta2-adrenergic receptor

Using sets of experimental distance restraints, which characterize active or inactive receptor conformations, and the X-ray crystal structure of the inactive form of bovine rhodopsin as a starting point, we have constructed models of both the active and inactive forms of rhodopsin and the beta2-adrenergic G-protein coupled receptors (GPCRs). The distance restraints were obtained from published data for site-directed crosslinking, engineered zinc binding, site-directed spin-labeling, IR spectroscopy, and cysteine accessibility studies conducted on class A GPCRs. Molecular dynamics simulations in the presence of either "active" or "inactive" restraints were used to generate two distinguishable receptor models. The process for generating the inactive and active models was validated by the hit rates, yields, and enrichment factors determined for the selection of antagonists in the inactive model and for the selection of agonists in the active model from a set of nonadrenergic GPCR drug-like ligands in a virtual screen using ligand docking software. The simulation results provide new insights into the relationships observed between selected biochemical data, the crystal structure of rhodopsin, and the structural rearrangements that occur during activation.     

Biased agonism of the calcium-sensing receptor

After the discovery of molecules modulating G protein-coupled receptors (GPCRs) that are able to selectively affect one signaling pathway over others for a specific GPCR, thereby "biasing" the signaling, it has become obvious that the original model of GPCRs existing in either an "on" or "off" conformation is too simple. The current explanation for this biased agonism is that GPCRs can adopt multiple active conformations stabilized by different molecules, and that each conformation affects intracellular signaling in a different way. In the present study we sought to investigate biased agonism of the calcium-sensing receptor (CaSR), by looking at 12 well-known orthosteric CaSR agonists in 3 different CaSR signaling pathways: G(q/11) protein, G(i/o) protein, and extracellular signal-regulated kinases 1 and 2 (ERK1/2). Here we show that apart from G(q/11) and G(i/o) signaling, ERK1/2 is activated through recruitment of β-arrestins. Next, by measuring activity of all three signaling pathways we found that barium, spermine, neomycin, and tobramycin act as biased agonist in terms of efficacy and/or potency. Finally, polyamines and aminoglycosides in general were biased in their potencies toward ERK1/2 signaling. In conclusion, the results of this study indicate that several active conformations of CaSR, stabilized by different molecules, exist, which affect intracellular signaling distinctly.     

Putative Active States of a Prototypic G-Protein-Coupled Receptor from Biased Molecular Dynamics

A major current focus of structural work on G-protein-coupled receptors (GPCRs) pertains to the investigation of their active states. However, for virtually all GPCRs, active agonist-bound intermediate states have been difficult to characterize experimentally owing to their higher conformational flexibility, and thus intrinsic instability, as compared to inactive inverse agonist-bound states. In this work, we explored possible activation pathways of the prototypic GPCR bovine rhodopsin by means of biased molecular dynamics simulations. Specifically, we used an explicit atomistic representation of the receptor and its environment, and sampled the conformational transition from the crystal structure of a photoactivated deprotonated state of rhodopsin to the low pH crystal structure of opsin in the presence of 11-trans-retinal, using adiabatic biased molecular dynamics simulations. We then reconstructed the system free-energy landscape along the predetermined transition trajectories using a path collective variable approach based on metadynamics. Our results suggest that the two experimental endpoints of rhodopsin/opsin are connected by at least two different pathways, and that the conformational transition is populated by at least four metastable states of the receptor, characterized by a different amplitude of the outward movement of transmembrane helix 6.     

Conserved activation pathways in G-protein-coupled receptors

GPCRs (G-protein-coupled receptors) are seven-transmembrane helix proteins that transduce exogenous and endogenous signals to modulate the activity of downstream effectors inside the cell. Despite the relevance of these proteins in human physiology and pharmaceutical research, we only recently started to understand the structural basis of their activation mechanism. In the period 2008-2011, nine active-like structures of GPCRs were solved. Among them, we have determined the structure of light-activated rhodopsin with all the features of the active metarhodopsin-II, which represents so far the most native-like model of an active GPCR. This structure, together with the structures of other inactive, intermediate and active states of rhodopsin constitutes a unique structural framework on which to understand the conserved aspects of the activation mechanism of GPCRs. This mechanism can be summarized as follows: retinal isomerization triggers a series of local structural changes in the binding site that are amplified into three intramolecular activation pathways through TM (transmembrane helix) 5/TM3, TM6 and TM7/TM2. Sequence analysis strongly suggests that these pathways are conserved in other GPCRs. Differential activation of these pathways by ligands could be translated into the stabilization of different active states of the receptor with specific signalling properties.     

Elucidating structural and molecular mechanisms of β-arrestin-biased agonism at GPCRs via MS-based proteomics

The discovery of β-arrestin-dependent GPCR signaling has led to an exciting new field in GPCR pharmacology: to develop “biased agonists” that can selectively target a specific downstream signaling pathway that elicits beneficial therapeutic effects without activating other pathways that elicit negative side effects. This new trend in GPCR drug discovery requires us to understand the structural and molecular mechanisms of β-arrestin-biased agonism, which largely remain unclear. We have used cutting-edge mass spectrometry (MS)-based proteomics, combined with systems, chemical and structural biology to study protein function, macromolecular interaction, protein expression and posttranslational modifications in the β-arrestin-dependent GPCR signaling. These high-throughput proteomic studies have provided a systems view of β-arrestin-biased agonism from several perspectives: distinct receptor phosphorylation barcode, multiple receptor conformations, distinct β-arrestin conformations, and ligand-specific signaling. The information obtained from these studies offers new insights into the molecular basis of GPCR regulation by β-arrestin and provides a potential platform for developing novel therapeutic interventions through GPCRs.     

G protein-coupled receptors of class A harness the energy of membrane potential to increase their sensitivity and selectivity

The human genome contains about 700 genes of G protein-coupled receptors (GPCRs) of class A; these seven-helical membrane proteins are the targets of almost half of all known drugs. In the middle of the helix bundle, crystal structures reveal a highly conserved sodium-binding site, which is connected with the extracellular side by a water-filled tunnel. This binding site contains a sodium ion in those GPCRs that are crystallized in their inactive conformations but does not in those GPCRs that are trapped in agonist-bound active conformations. The escape route of the sodium ion upon the inactive-to-active transition and its very direction have until now remained obscure. Here, by modeling the available experimental data, we show that the sodium gradient over the cell membrane increases the sensitivity of GPCRs if their activation is thermodynamically coupled to the sodium ion translocation into the cytoplasm but decreases it if the sodium ion retreats into the extracellular space upon receptor activation. The model quantitatively describes the available data on both activation and suppression of distinct GPCRs by membrane voltage. The model also predicts selective amplification of the signal from (endogenous) agonists if only they, but not their (partial) analogs, induce sodium translocation. Comparative structure and sequence analyses of sodium-binding GPCRs indicate a key role for the conserved leucine residue in the second transmembrane helix (Leu2.46) in coupling sodium translocation to receptor activation. Hence, class A GPCRs appear to harness the energy of the transmembrane sodium potential to increase their sensitivity and selectivity.     

Quantifying conformational changes in GPCRs: glimpse of a common functional mechanism

Background

G-protein-coupled receptors (GPCRs) are important drug targets and a better understanding of their molecular mechanisms would be desirable. The crystallization rate of GPCRs has accelerated in recent years as techniques have become more sophisticated, particularly with respect to Class A GPCRs interacting with G-proteins. These developments have made it possible for a quantitative analysis of GPCR geometrical features and binding-site conformations, including a statistical comparison between Class A GPCRs in active (agonist-bound) and inactive (antagonist-bound) states.

Results

Here we implement algorithms for the analysis of interhelical angles, distances, interactions and binding-site volumes in the transmembrane domains of 25 Class A GPCRs (7 active and 18 inactive). Two interhelical angles change in a statistically significant way between average inactive and active states: TM3-TM6 (by -9°) and TM6-TM7 (by +12°). A third interhelical angle: TM5-TM6 shows a trend, changing by -9°. In the transition from inactive to active states, average van der Waals interactions between TM3 and TM7 significantly increase as the average distance between them decreases by >2 Å. Average H-bonding between TM3 and TM6 decreases but is seemingly compensated by an increase in H-bonding between TM5 and TM6. In five Class A GPCRs, crystallized in both active and inactive states, increased H-bonding of agonists to TM6 and TM7, relative to antagonists, is observed. These protein-agonist interactions likely favour a change in the TM6-TM7 angle, which creates a narrowing in the binding pocket of activated receptors and an average ~200 ų reduction in volume.

Conclusions

In terms of similar conformational changes and agonist binding pattern, Class A GPCRs appear to share a common mechanism of activation, which can be exploited in future drug development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0567-3) contains supplementary material, which is available to authorized users.     

Structural Features of β2 Adrenergic Receptor: Crystal Structures and Beyond

The beta2-adrenergic receptor (β2AR) family, which is the largest family of cell surface receptors in humans. Extra attention has been focused on the human GPCRs because they have been studied as important protein targets for pharmaceutical drug development. In fact, approximately 40% of marketed drugs directly work on GPCRs. GPCRs respond to various extracellular stimuli, such as sensory signals, neurotransmitters, chemokines, and hormones, to induce structural changes at the cytoplasmic surface, activating downstream signaling pathways, primarily through interactions with heterotrimeric G proteins or through G-protein independent pathways, such as arrestin. Most GPCRs, except for rhodhopsin, which contains covalently linked 11 cis-retinal, bind to diffusible ligands, having various conformational states between inactive and active structures. The first human GPCR structure was determined using an inverse agonist bound β2AR in 2007 and since then, more than 20 distinct GPCR structures have been solved. However, most GPCR structures were solved as inactive forms, and an agonist bound fully active structure is still hard to obtain. In a structural point of view, β2AR is relatively well studied since its fully active structure as a complex with G protein as well as several inactive structures are available. The structural comparison of inactive and active states gives an important clue in understanding the activation mechanism of β2AR. In this review, structural features of inactive and active states of β2AR, the interaction of β2AR with heterotrimeric G protein, and the comparison with β1AR will be discussed.     

Bioluminescence resonance energy transfer assays reveal ligand-specific conformational changes within preformed signaling complexes containing delta-opioid receptors and heterotrimeric G proteins

Heptahelical receptors communicate extracellular information to the cytosolic compartment by binding an extensive variety of ligands. They do so through conformational changes that propagate to intracellular signaling partners as the receptor switches from a resting to an active conformation. This active state has been classically considered unique and responsible for regulation of all signaling pathways controlled by a receptor. However, recent functional studies have challenged this notion and called for a paradigm where receptors would exist in more than one signaling conformation. This study used bioluminescence resonance energy transfer assays in combination with ligands of different functional profiles to provide in vivo physical evidence of conformational diversity of delta-opioid receptors (DORs). DORs and alpha(i1)beta(1)gamma(2) G protein subunits were tagged with Luc or green fluorescent protein to produce bioluminescence resonance energy transfer pairs that allowed monitoring DOR-G protein interactions from different vantage points. Results showed that DORs and heterotrimeric G proteins formed a constitutive complex that underwent structural reorganization upon ligand binding. Conformational rearrangements could not be explained by a two-state model, supporting the idea that DORs adopt ligand-specific conformations. In addition, conformational diversity encoded by the receptor was conveyed to the interaction among heterotrimeric subunits. The existence of multiple active receptor states has implications for the way we conceive specificity of signal transduction.     

Active state-like conformational elements in the beta2-AR and a photoactivated intermediate of rhodopsin identified by dynamic properties of GPCRs

G-Protein-coupled receptors (GPCRs) adopt various functionally relevant conformational states in cell signaling processes. Recently determined crystal structures of rhodopsin and the beta 2-adrenergic receptor (beta 2-AR) offer insight into previously uncharacterized active conformations, but the molecular states of these GPCRs are likely to contain both inactive and active-like conformational elements. We have identified conformational rearrangements in the dynamics of the TM7-HX8 segment that relate to the properties of the conserved NPxxY(x)5,6F motif and show that they can be used to identify active state-like conformational elements in the corresponding regions of the new structures of rhodopsin and the beta 2-AR.     

Allosteric Effects of Sodium Ion Binding on Activation of the M3 Muscarinic G-Protein-Coupled Receptor

G-protein-coupled receptors (GPCRs) are important membrane proteins that mediate cellular signaling and represent primary targets for about one-third of currently marketed drugs. Recent x-ray crystallographic studies identified distinct conformations of GPCRs in the active and inactive states. An allosteric sodium ion was found bound to a highly conserved D2.50 residue in inactive GPCRs, whereas the D2.50 allosteric pocket became collapsed in active GPCR structures. However, the dynamic mechanisms underlying these observations remain elusive. In this study, we aimed to understand the mechanistic effects of sodium ion binding on dynamic activation of the M3 muscarinic GPCR through long-timescale accelerated molecular dynamics (aMD) simulations. Results showed that with the D2.50 residue deprotonated, the M3 receptor is bound by an allosteric sodium ion and confined mostly in the inactive state with remarkably reduced flexibility. In contrast, the D2.50-protonated receptor does not exhibit sodium ion binding to the D2.50 allosteric site and samples a significantly larger conformational space. The receptor activation is captured and characterized by large-scale structural rearrangements of the transmembrane helices via dynamic hydrogen bond and salt bridge interactions. The residue motions are highly correlated during receptor activation. Further network analysis revealed that the allosteric signaling between residue D2.50 and key residues in the intracellular, extracellular, and orthosteric pockets is significantly weakened upon sodium ion binding.     

The CholecystokininB Receptor Is Coupled to Two Effector Pathways Through Pertussis Toxin-Sensitive and -Insensitive G Proteins

Previous binding studies have suggested the existence of two affinity states for type B cholecystokinin receptors (CCK(B)R), which could correspond to different coupling states of the receptor to G proteins. To test this hypothesis, we have further investigated signal transduction pathways coupled to rat CCK(B)R stably transfected in Chinese hamster ovary cells. We show that CCK(B)R are coupled to two distinct transduction pathways involving two different G proteins, a pertussis toxin-insensitive/phospholipase C pathway leading to the production of inositol phosphate and arachidonic acid, and a pertussis toxin-sensitive/phospholipase A2 pathway leading to the release of arachidonic acid. We further demonstrate that the relative degree of activation of each effector pathway by different specific CCK(B)R agonists is the same, and that a specific CCK(B)R antagonist, RB213, can differentially antagonize the two signal transduction pathways elicited by these agonists. Taken all together, these data could be explained by the recently proposed theory assuming that the receptor can exist in a three-state model in which two active conformations corresponding to the complex formed by the receptor with two different G proteins coexist. According to this model, agonists or antagonists could recognize preferentially either conformation of the activated receptor, leading to variable behavior in a system containing a single receptor type.     

Ligand based conformational space studies of the μ-opioid receptor

BackgroundG protein-coupled receptors (GPCRs) comprise a family of membrane proteins that can be activated by a variety of external factors. The μ-opioid receptor (MOR), a class A GPCR, is the main target of morphine. Recently, enhanced sampling molecular dynamics simulations of a constitutively active mutant of MOR in its apo form allowed us to capture the novel intermediate states of activation, as well as the active state. This prompted us to apply the same techniques to wild type MOR in complex with ligands, in order to explore their contributions to the receptor conformational changes in the activation process.MethodsMOR was modeled in complex with agonists (morphine, BU72), a partial agonist (naloxone benzoylhydrazone) and an antagonist (naloxone). Replica exchange with solute tempering (REST2) molecular dynamics simulations were carried out for all systems. Trajectory frames were clustered, and the activation state of each cluster was assessed by two different methods.ResultsCluster sizes and activation indices show that while agonists stabilized structures in a higher activation state, the antagonist behaved oppositely. Morphine tends to drive the receptor towards increasing R165-T279 distances, while naloxone tends to increase the NPxxYA motif conformational change.ConclusionsDespite not observing a full transition between inactive and active states, an important conformational change of transmembrane helix 5 was observed and associated with a ligand-driven step of the process.General significanceThe activation process of GPCRs is widely studied but still not fully understood. Here we carried out a step forward in the direction of gaining more details of this process.     

Generation of an agonistic binding site for blockers of the M(3) muscarinic acetylcholine receptor

GPCRs (G-protein-coupled receptors) exist in a spontaneous equilibrium between active and inactive conformations that are stabilized by agonists and inverse agonists respectively. Because ligand binding of agonists and inverse agonists often occurs in a competitive manner, one can assume an overlap between both binding sites. Only a few studies report mutations in GPCRs that convert receptor blockers into agonists by unknown mechanisms. Taking advantage of a genetically modified yeast strain, we screened libraries of mutant M(3)Rs {M(3) mAChRs [muscarinic ACh (acetylcholine) receptors)]} and identified 13 mutants which could be activated by atropine (EC50 0.3-10 microM), an inverse agonist on wild-type M(3)R. Many of the mutations sensitizing M(3)R to atropine activation were located at the junction of intracellular loop 3 and helix 6, a region known to be involved in G-protein coupling. In addition to atropine, the pharmacological switch was found for other M(3)R blockers such as scopolamine, pirenzepine and oxybutynine. However, atropine functions as an agonist on the mutant M(3)R only when expressed in yeast, but not in mammalian COS-7 cells, although high-affinity ligand binding was comparable in both expression systems. Interestingly, we found that atropine still blocks carbachol-induced activation of the M(3)R mutants in the yeast expression system by binding at the high-affinity-binding site (Ki approximately 10 nM). Our results indicate that blocker-to-agonist converting mutations enable atropine to function as both agonist and antagonist by interaction with two functionally distinct binding sites.     

Modeling activation and desensitization of G-protein coupled receptors provides insight into ligand efficacy.

Signaling through G-protein coupled receptors is one of the most prevalent and important methods of transmitting information to the inside of cells. Many mathematical models have been proposed to describe this type of signal transduction, and the ternary complex (ligand/receptor/G-protein) model and its derivatives are among the most widely accepted. Current versions of these equilibrium models include both active (i.e. signaling) and inactive conformations of the receptor, but do not include the dynamics of G-protein activation or receptor desensitization. Yet understanding how these dynamic events effect response behavior is crucial to determining ligand efficacy. We developed a mathematical model for G-protein coupled receptor signaling that includes G-protein activation and receptor desensitization, and used it to predict how activation and desensitization would change if either the conformational selectivity (the effect of ligand binding on the distribution of active and inactive receptor states) or the desensitization rate constant was ligand-specific. In addition, the model was used to explore the implications of measuring responses far downstream from G-protein activation. By comparing the experimental data from the beta(2)-adrenergic, micro-opioid, D(1)dopamine, and neutrophil N -formyl peptide receptors with the predictions of our model, we found that the conformational selectivity is the predominant factor in determining the amounts of activation and desensitization caused by a particular ligand.     

Deciphering conformational selectivity in the A2A adenosine G protein-coupled receptor by free energy simulations

Transmembranal G Protein-Coupled Receptors (GPCRs) transduce extracellular chemical signals to the cell, via conformational change from a resting (inactive) to an active (canonically bound to a G-protein) conformation. Receptor activation is normally modulated by extracellular ligand binding, but mutations in the receptor can also shift this equilibrium by stabilizing different conformational states. In this work, we built structure-energetic relationships of receptor activation based on original thermodynamic cycles that represent the conformational equilibrium of the prototypical A_2A adenosine receptor (AR). These cycles were solved with efficient free energy perturbation (FEP) protocols, allowing to distinguish the pharmacological profile of different series of A_2AAR agonists with different efficacies. The modulatory effects of point mutations on the basal activity of the receptor or on ligand efficacies could also be detected. This methodology can guide GPCR ligand design with tailored pharmacological properties, or allow the identification of mutations that modulate receptor activation with potential clinical implications.     

Direct assessment of CXCR4 mutant conformations reveals complex link between receptor structure and G(alpha)(i) activation

Ligand binding to G protein-coupled receptors (GPCRs) is thought to induce changes in receptor conformation that translate into activation of downstream effectors. The link between receptor conformation and activity is still insufficiently understood, as current models of GPCR activation fail to take an increasing amount of experimental data into account. To elucidate structure-function relationships in GPCR activation, we used bioluminescence resonance energy transfer to directly assess the conformation of mutants of the chemokine receptor CXCR4. We analyzed substitutions in the arginine cage DRY motif and in the conserved asparagine N(3.35)119, which are pivotal molecular switches for receptor conformation and activation. G(alpha)(i) activation of the mutants was either similar to wild-type CXCR4 (D133N, Y135A, and N119D) or resulted in loss of activity (R134A and N119K). Mutant N119S was constitutively active but further activated by agonist. Bioluminescence resonance energy transfer analysis suggested no simple correlation between conformational changes in response to ligand binding and activation of G(alpha)(i) by the mutants. Different conformations of active receptors were detected (for wild-type CXCR4, D133N, and N119S), suggesting that different receptor conformations are able to trigger G(alpha)(i) activity. Several conformations were also found for inactive mutants. These data provide biophysical evidence for different receptor conformations being active with respect to a single readout. They support models of GPCR structure-activity relationships that take this conformational flexibility of active receptors into account.