Theoretical study of the electrostatically driven step of receptor-G protein recognition

Constitutive activity and inverse agonism at the α(₁a) and α(₁b) adrenergic receptor subtypes

The α(1b)-adrenergic receptor (AR) was, after rhodopsin, the first G protein-coupled receptor (GPCR) in which point mutations were shown to trigger constitutive (agonist-independent) activity. Constitutively activating mutations have been found in other AR subtypes as well as in several GPCRs. This chapter briefly summarizes the main findings on constitutively active mutants of the α(1a)- and α(1b)-AR subtypes and the methods used to predict activating mutations, to measure constitutive activity of Gq-coupled receptors and to investigate inverse agonism. In addition, it highlights the implications of studies on constitutively active AR mutants on elucidating the molecular mechanisms of receptor activation and drug action.

Coupling of beta2-adrenoceptors to XLalphas and Galphas: a new insight into ligand-induced G protein activation

Galpha(s) and extra-large Galpha(s) (XLalpha(s)) can both transduce receptor activation into intracellular cAMP generation. It is unknown, however, whether these two GNAS-locus products display distinct properties with respect to receptor coupling. Here, we show that XLalpha(s) couples to the beta2-adrenoceptor more efficiently than Galpha(s). In transfected human embryonic kidney 293 cells and mouse embryonic fibroblasts null for both Galpha(s) and XLalpha(s) (2B2 cells), basal cAMP accumulation mediated by XLalpha(s) was higher than that mediated by Galpha(s). Inverse agonist treatment reduced Galpha(s)-mediated basal activity, whereas its effect was markedly blunted on XLalpha(s)-mediated basal activity. Rank order of ligand efficacies regarding cAMP accumulation was the same when the receptor was coupled to XLalpha(s) or Galpha(s). However, ligand-induced and XLalpha(s)-mediated cAMP generation was higher than that mediated by Galpha(s). The relatively high efficiency of XLalpha(s)-mediated cAMP generation was conditional, disappearing with increased level of receptor expression or increased efficacy of ligand. Full-agonist responses in XLalpha(s)- and Galpha(s)-expressing cells were comparable even at low receptor levels, whereas partial agonist responses became comparable only when the receptor expression was increased (>3 pmol/mg). Radioligand binding studies showed that the high-affinity component in agonist binding to beta2-adrenoceptor was more pronounced in cells expressing XLalpha(s) than those expressing Galpha(s). We discuss these findings in the framework of current receptor-G protein activation models and offer an extended ternary complex model that can fully explain our observations.

Exploring allosteric coupling in the alpha-subunit of heterotrimeric G proteins using evolutionary and ensemble-based approaches

BACKGROUND

Allosteric coupling, which can be defined as propagation of a perturbation at one region of the protein molecule (such as ligand binding) to distant sites in the same molecule, constitutes the most general mechanism of regulation of protein function. However, unlike molecular details of ligand binding, structural elements involved in allosteric effects are difficult to diagnose. Here, we identified allosteric linkages in the alpha-subunits of heterotrimeric G proteins, which were evolved to transmit membrane receptor signals by allosteric mechanisms, by using two different approaches that utilize fundamentally different and independent information.

RESULTS

We analyzed: 1) correlated mutations in the family of G protein alpha-subunits, and 2) cooperativity of the native state ensemble of the Galphai1 or transducin. The combination of these approaches not only recovered already-known details such as the switch regions that change conformation upon nucleotide exchange, and those regions that are involved in receptor, effector or Gbetagamma interactions (indicating that the predictions of the analyses can be viewed with a measure of confidence), but also predicted new sites that are potentially involved in allosteric communication in the Galpha protein. A summary of the new sites found in the present analysis, which were not apparent in crystallographic data, is given along with known functional and structural information. Implications of the results are discussed.

CONCLUSION

A set of residues and/or structural elements that are potentially involved in allosteric communication in Galpha is presented. This information can be used as a guide to structural, spectroscopic, mutational, and theoretical studies on the allosteric network in Galpha proteins, which will provide a better understanding of G protein-mediated signal transduction.

Dark and photoactivated rhodopsin share common binding modes to transducin

The structure of the photoactivated deprotonated rhodopsin intermediate was compared with two different structures of dark rhodopsin. Structure comparisons relied on the computation of molecular indices and on docking simulations with heterotrimeric transducin (Gt). The results of this study provide the first evidence that dark and photoactivated rhodopsins share a common recognition mode to Gt, characterized by the docking of the Gt alpha C-tail in the proximity to the E/DRY motif of rhodopsin.

Modeling the human oxytocin receptor for drug discovery efforts

The oxytocin receptor belongs to class A receptors within the great family of G protein-coupled receptors. The endogenous ligand oxytocin is a nonapeptide hormone that induces the uterine contractions at parturition and is used to induce the labor. The peptide oxytocin and, even more, its non-peptide antagonist, could be valuable tools in tocolysis. The knowledge of the three-dimensional structure of the oxytocin receptor and the determination of the main interaction points between the receptor and the ligands may help to develop selective oxytocin agonists and antagonist. This review summarizes the knowledge about the mapping of the binding domain of the oxytocin receptor and the efforts in the field of molecular modeling studies related to oxytocin receptor-ligand interactions.

Contributions of intracellular loops 2 and 3 of the lutropin receptor in Gs coupling

A number of amino acids essential for Gs coupling, i.e. hot spots, were identified after in vitro Ala-scanning mutagenesis of the cytosolic extensions of helices 3, 5, and 6 and of intracellular loops 2 and 3 (IL2 and IL3) of the human LH receptor (LHR). Consistent with the results of in vitro experiments involving ligand binding and ligand-mediated signaling in transiently transfected human embryonic kidney 293 cells, computational modeling of the isolated receptor and of the receptor-G protein complexes suggests an important role of the cytosolic extension of helix 3 and the N-terminal portion of the IL2 in Gs(alpha) interaction, whereas the contribution of IL3 is marginal. Mapping the hot spots into the computational models of LHR and the LHR-Gs complexes allowed for a distinction between receptor sites required for intramolecular structural changes (i.e. I460, T461, H466, and I549) and receptor sites more likely involved in G protein recognition (i.e. R464, T467, I468, Y470, Y550, and D564). The latter sites include the highly conserved arginine of the (E/D)R(Y/W) motif, which is therefore likely to be a receptor recognition point for Gs rather than a switch of receptor activation. The results of in vitro and in silico experiments carried out in this study represent the first comprehensive delineation of functionality of the individual residues in the intracellular domains of LHR and establish potential switches of receptor activation as well as a map of the primary receptor recognition sites for Gs. A novel way to consider constitutively active mutants was inferred from this study, i.e. receptor states with improved complementarity for the G protein compared to the wild-type receptor.

Constitutive activity and inverse agonism at the α1adrenoceptors

Mutations of G protein-coupled receptors (GPCR) can increase their constitutive (agonist-independent) activity. Some of these mutations have been artificially introduced by site-directed mutagenesis, others occur spontaneously in human diseases. The alpha(1B)adrenoceptor was the first GPCR in which point mutations were shown to trigger receptor activation. This article briefly summarizes some of the findings reported in the last several years on constitutive activity of the alpha(1)adrenoceptor subtypes, the location where mutations have been found in the receptors, the spontaneous activity of native receptors in recombinant as well as physiological systems. In addition, it will highlight how the analysis of the pharmacological and molecular properties of the constitutively active adrenoceptor mutants provided an important contribution to our understanding of the molecular mechanisms underlying the mechanism of receptor activation and inverse agonism.

Common structural requirements for heptahelical domain function in class A and class C G protein-coupled receptors

G protein-coupled receptors (GPCRs) are key players in cell communication. Several classes of such receptors have been identified. Although all GPCRs possess a heptahelical domain directly activating G proteins, important structural and sequence differences within receptors from different classes suggested distinct activation mechanisms. Here we show that highly conserved charged residues likely involved in an interaction network between transmembrane domains (TM) 3 and 6 at the cytoplasmic side of class C GPCRs are critical for activation of the gamma-aminobutyric acid type B receptor. Indeed, the loss of function resulting from the mutation of the conserved lysine residue into aspartate or glutamate in the TM3 of gamma-aminobutyric acid type B(2) can be partly rescued by mutating the conserved acidic residue of TM6 into either lysine or arginine. In addition, mutation of the conserved lysine into an acidic residue leads to a nonfunctional receptor that displays a high agonist affinity. This is reminiscent of a similar ionic network that constitutes a lock stabilizing the inactive state of many class A rhodopsin-like GPCRs. These data reveal that despite their original structure, class C GPCRs share with class A receptors at least some common structural feature controlling G protein activation.

Inactive and active states and supramolecular organization of GPCRs: insights from computational modeling

Herein we make an overview of the results of our computational experiments aimed at gaining insight into the molecular mechanisms of GPCR functioning either in their normal conditions or when hit by gain-of-function or loss-of-function mutations. Molecular simulations of a number of GPCRs in their wild type and mutated as well as free and ligand-bound forms were instrumental in inferring the structural features, which differentiate the mutation- and ligand-induced active from the inactive states. These features essentially reside in the interaction pattern of the E/DRY arginine and in the degree of solvent exposure of selected cytosolic domains. Indeed, the active states differ from the inactive ones in the weakening of the interactions made by the highly conserved arginine and in the increase in solvent accessibility of the cytosolic interface between helices 3 and 6. Where possible, the structural hallmarks of the active and inactive receptor states are translated into molecular descriptors useful for in silico functional screening of novel receptor mutants or ligands. Computational modeling of the supramolecular organization of GPCRs and their intracellular partners is the current challenge toward a deep understanding of their functioning mechanisms.

Structural features of the inactive and active states of the melanin-concentrating hormone receptors: Insights from molecular simulations

Comparative molecular dynamics simulations of both subtypes 1 and 2 of the melanin-concentrating hormone receptor (MCHR1 and MCHR2, respectively) in their free and hormone-bound forms have been carried out. The hormone has been used in its full-length and truncated forms, as well as in 16 mutated forms. Moreover, MCHR1 has been simulated in complex with T-226296, a novel orally active and selective antagonist. The comparative analysis of an extended number of receptor configurations suggests that the differences between inactive (i.e., free and antagonist-bound) and active (i.e., agonist-bound) states of MCHRs involve the receptor portions close to the E/DRY and NPxxY motifs, with prominence to the cytosolic extensions of helices 2, 3, 6, and 7. In fact, the active forms of these receptors share the release of selected intramolecular interactions found in the inactive forms, such as that between R3.50 of the E/DRY motif and D2.40, and that between Y7.53 of the NPxxY motif and F7.60. Another feature of the active forms of both MCHRs is the approach of "helix 8" to the cytosolic extension of helix 3. These features of the active forms are concurrent with the opening of a cleft at the cytosolic end of the helix bundle. For both MCHRs, the agonist-induced chemical information transfer from the extracellular to the cytosolic domains is mediated by a cluster of aromatic amino acids in helix 6, following the ligand interaction with selected amino acids in the extracellular half of the receptor.

Sequence analysis reveals how G protein-coupled receptors transduce the signal to the G protein

Sequence entropy-variability plots based on alignments of very large numbers of sequences-can indicate the location in proteins of the main active site and modulator sites. In the previous article in this issue, we applied this observation to a series of well-studied proteins and concluded that it was possible to detect most of the residues with a known functional role. Here, we apply the method to rhodopsin-like G protein-coupled receptors. Our conclusion is that G protein binding is the main evolutionary constraint on these receptors, and that other ligands, such as agonists, act as modulators. The activation of the receptors can be described as a simple, two-step process, and the residues involved in signal transduction can be identified.

Molecular Dynamics Simulations of the Ligand-Induced Chemical Information Transfer in the 5-HT1A Receptor

Comparative molecular dynamics simulations of the 5-HT(1A) receptor in its empty as well as agonist- (i.e. active) and antagonist-bound (i.e. nonactive) forms have been carried out. The agonists 5-HT and (R)-8-OH-DPAT as well as the antagonist WAY100635 have been employed. The results of this study strengthen the hypothesis that the receptor portions close to the E/DRY/W motif, with prominence to the cytosolic extensions of helices 3 and 6, are particularly susceptible to undergo structural modification in response to agonist binding. Despite the differences in the structural/dynamics behavior of the two agonists when docked into the 5-HT(1A) receptor, they both exert a destabilization of the intrahelical and interhelical interactions found in the empty and antagonist-bound receptor forms between the arginine of the E/DRY sequence and both D133(3.49) and E340(6.30). For both agonists, the chemical information transfer from the extracellular to the cytosolic domains is mediated by a cluster of aromatic amino acids in helix 6, following the ligand interaction with selected amino acids in the extracellular half of the receptor, such as D116(3.32), S199(5.42), Y195(5.38), and F361(6.51). A significant reduction in the bend at P360(6.50), as compared to the empty and the antagonist-bound receptor forms, is one of the features of the agonist-bound forms that is related to the breakage of the interhelical salt bridge between the E/DRY arginine and E340(6.30). Another structural feature, shared by the agonist-bound receptor forms and not by the empty and antagonist-bound forms, is the detachment of helices 2 and 4, as marked by the movement of W161(4.50) away from helix 2, toward the membrane space.

Opioid receptor random mutagenesis reveals a mechanism for G protein-coupled receptor activation

The high resolution structure of rhodopsin has greatly enhanced current understanding of G protein-coupled receptor (GPCR) structure in the off-state, but the activation process remains to be clarified. We investigated molecular mechanisms of delta-opioid receptor activation without a preconceived structural hypothesis. Using random mutagenesis of the entire receptor, we identified 30 activating point mutations. Three-dimensional modeling revealed an activation path originating from the third extracellular loop and propagating through tightly packed helices III, VI and VII down to a VI-VII cytoplasmic switch. N- and C-terminal determinants also influence receptor activity. Findings for this therapeutically important receptor may apply to other GPCRs that respond to diffusible ligands.

The role of a conserved region of the second intracellular loop in AT1 angiotensin receptor activation and signaling

The pleiotropic actions of angiotensin II are mediated by the primarily G(q) protein-coupled type 1 angiotensin (AT(1)) receptor. In this study a mutational analysis of the function of the conserved DRYXXV/IXXPL domain in the second intracellular loop of the rat AT(1A) receptor was performed in COS7 cells. Alanine substitution studies showed that single replacement of the highly conserved Asp(125) and Arg(126), but not Tyr(127), moderately impaired angiotensin II-induced inositol phosphate signaling. However, concomitant substitution of both Asp(125) and Arg(126) caused marked reduction of both inositol phosphate signaling and receptor internalization. Alanine scanning of the adjacent residues showed that substitution of Ile(130), His(132), and Pro(133) reduced agonist-induced inositol phosphate signal generation, whereas mutations of Met(134) also impaired receptor internalization. Expression of the D125A mutant AT(1A) receptor in COS7 cells endowed the receptor with moderate constitutive activity, as indicated by its enhanced basal Elk1 promoter activity and inositol phosphate response to partial agonists. Angiotensin II-induced stimulation of the Elk1 promoter showed parallel impairment with inositol phosphate signal generation in receptors containing mutations in this region of the AT(1A) receptor. These data confirm that Ca(2+) signal generation is required for the nuclear effects of angiotensin II-induced ERK activation. They are also consistent with the role of the conserved DRY sequence of the AT(1A) receptor in receptor activation, and of Asp(125) in constraining the receptor in its inactive conformation. Furthermore, in the cytoplasmic helical extension of the third helix, an apolar surface that includes Ile(130) and Met(134) appears to have a direct role in G protein coupling.

Promiscuous coupling at receptor-Galpha fusion proteins. The receptor of one covalent complex interacts with the alpha-subunit of another

Fusion proteins between heptahelical receptors (GPCR) and G protein alpha-subunits show enhanced signaling efficiency in transfected cells. This is believed to be the result of molecular proximity, because the interaction between linked modules of one protein chain, if not constrained by structure, should be strongly favored compared with the same in which partners react as free species. To test this assumption we made a series of fusion proteins (type 1 and 4 opioid receptors with G(o) and beta(2) adrenergic and dopamine 1 receptors with G(sL)) and some mutated analogs carrying different tags and defective GPCR or Galpha subunits. Using cotransfection experiments with readout protocols able to distinguish activation at fused and non-fused alpha-subunits, we found that both the GPCR and the Galpha limb of one fusion protein can freely interact with non-fused proteins and the tethered partners of a neighboring fusion complex. Moreover, a bulky polyanionic inhibitor can suppress with identical potency receptor-Galpha interaction, either when occurring between latched domains of a fused system or separate elements of distinct molecules, indicating that the binding surfaces are equally accessible in both cases. These data demonstrate that there is no entropy drive from the linked condition of fusion proteins and suggest that their signaling may result from the GPCR of one complex interacting with the alpha-subunit of another. Moreover, the enhanced coupling efficiency commonly observed for fusion proteins is not due to the receptor tether, but to the transmembrane helix that anchors Galpha to the membrane.

Mutagenesis and modelling of the alpha(1b)-adrenergic receptor highlight the role of the helix 3/helix 6 interface in receptor activation

Computer simulations on a new model of the alpha1b-adrenergic receptor based on the crystal structure of rhodopsin have been combined with experimental mutagenesis to investigate the role of residues in the cytosolic half of helix 6 in receptor activation. Our results support the hypothesis that a salt bridge between the highly conserved arginine (R143(3.50)) of the E/DRY motif of helix 3 and a conserved glutamate (E289(6.30)) on helix 6 constrains the alpha1b-AR in the inactive state. In fact, mutations of E289(6.30) that weakened the R143(3.50)-E289(6.30) interaction constitutively activated the receptor. The functional effect of mutating other amino acids on helix 6 (F286(6.27), A292(6.33), L296(6.37), V299(6.40,) V300(6.41), and F303(6.44)) correlates with the extent of their interaction with helix 3 and in particular with R143(3.50) of the E/DRY sequence.

The lutropin/choriogonadotropin receptor, a 2002 perspective

Reproduction cannot take place without the proper functioning of the lutropin/choriogonadotropin receptor (LHR). When the LHR does not work properly, ovulation does not occur in females and Leydig cells do not develop normally in the male. Also, because the LHR is essential for sustaining the elevated levels of progesterone needed to maintain pregnancy during the first trimester, disruptions in the functions of the LHR during pregnancy have catastrophic consequences. As such, a full understanding of the biology of the LHR is essential to the survival of our species. In this review we summarize our current knowledge of the structure, functions, and regulation of this important receptor.

Suramin affects coupling of rhodopsin to transducin

Suramin, a polysulfonated naphthylurea, is under investigation for the treatment of several cancers. It interferes with signal transduction through G(s), G(i), and G(o), but structural and kinetic aspects of the molecular mechanism are not well understood. Here, we have investigated the influence of suramin on coupling of bovine rhodopsin to G(t), where G-protein activation and receptor structure can be monitored by spectroscopic in vitro assays. G(t) fluorescence changes in response to rhodopsin-catalyzed nucleotide exchange reveal that suramin inhibits G(t) activation by slowing down the rate of complex formation between metarhodopsin-II and G(t). The metarhodopsin-I/-II photoproduct equilibrium, GTPase activity, and nucleotide uptake by G(t) are unaffected. Attenuated total reflection Fourier transform infrared spectroscopy shows that the structure of rhodopsin, metarhodopsin-II, and the metarhodopsin-II G(t) complex is also not altered. Instead, suramin dissociates G(t) from disk membranes in the dark, whereas metarhodopsin-II G(t) complexes are stable. Förster resonance energy transfer suggests a suramin-binding site near Trp(207) on the G(t alpha) subunit (K(d) approximately 0.5 microM). The kinetic analyses and the structural data are consistent with a specific perturbation by suramin of the membrane attachment site on G(t alpha). Disruption of membrane anchoring may contribute to some of the effects of suramin exerted on other G-proteins.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

Quasi-irreversible binding of agonist to beta-adrenoceptors and formation of non-dissociating receptor-G(s) complex in the absence of guanine nucleotides

Here, we tested the hypothesis that receptor-G protein and agonist may form an irreversible complex in the absence of guanine nucleotides. We used the beta-adrenoceptor-G(s) system of guinea pig lung parenchymal membranes as a model. Two groups of membranes were used in the experiments: (1) washed with nucleotide-free buffer in the presence of isoproterenol (isoproterenol-treated), and (2) washed with buffer alone or with agonist+GDP (both were treated as control). Results were as follows: (1) the iodopindolol binding capacity of isoproterenol-treated membranes was reduced by about 30%. (2) No such reduction was observed in control membranes. (3) Addition of GDP to the isoproterenol-treated membranes completely restored the pindolol binding capacity. We interpreted this result as indicating irreversible agonist-receptor complex is formed when the receptor interacts with nucleotide-free G(salpha). (4) We observed a single peak of beta(2)-adrenoceptor activity in the control group by size-exclusion chromatography of the solubilized membranes. Inclusion of isoproterenol in the washing buffer led to an additional (heavier) peak of beta(2)-adrenoceptor activity. This peak disappeared when GDP was added to the detergent extract before high-pressure liquid chromatography (HPLC) analysis. Western blot analysis of these HPLC fractions showed that the agonist-induced heavier peak contained significantly more G(salpha) protein than did the other fractions. We interpreted this result as indicating that a practically irreversible complex of receptor and G protein is formed in the absence of GDP. We suggest that the tightly bound (nucleotide-free) receptor-G protein complex also contains the agonist, and that this complex can be reversed only by the addition of nucleotides. The implications of these results are also discussed.

Differences in conformational properties of the second intracellular loop (IL2) in 5HT(2C) receptors modified by RNA editing can account for G protein coupling efficiency

Adenosine-to-inosine RNA editing events that have been demonstrated for 5HT (2C) receptors resulted in alterations of the amino acid sequence at positions 156, 158 and 160 in the intracellular loop 2 (IL2) region. The edited receptor isoforms were shown to have reduced basal activity, but similar maximum responses to agonist binding. To identify the molecular mechanism of these pharmacological effects of editing we explored the conformational properties of the edited IL2 in comparison with the wild type. The results from conformational studies of the IL2 isoforms, using biased Monte Carlo simulations with an implicit solvent model based on a screened Coulomb potential, show that the compared loops differ in their preferred spatial orientations as a result of differences in the conformational space that is accessible to them by energy criteria. For the IL2 of the unedited (5HT (2C-INI) ) receptor, the preference for structures oriented towards the 7TM bundle is larger than for the 5HT (2C-VGV) edited receptor. This difference in preferred orientation can affect the association of IL2 with other intracellular loop domains involved in G protein coupling and hence the coupling efficiency. The results illustrate the high sensitivity of the system to small changes in the interaction surface presented to other intracellular loops, and/or the G protein.

Theoretical study on mutation-induced activation of the luteinizing hormone receptor

Here, three-dimensional model building and molecular dynamics simulations of the luteinizing hormone receptor have been employed to generate hypotheses about the molecular mechanisms underlying the activation of the receptor induced by naturally occurring activating mutations. The comparative analysis of the wild-type receptor and of 16 constitutively active or inactive mutants has been instrumental in inferring the structural/dynamic features which could characterize the inactive and the active forms of the receptor. These features have been also employed for predicting the functional behavior of new receptor mutants. The results of this study might provide a structural framework to interpret the pathological effects induced by mutations of the luteinizing hormone receptor. In addition, the proposed theoretical model could be useful for engineering new mutations or ligands able to modulate receptor function.