The role of helices 5 and 6 on the human β1-adrenoceptor activation mechanism

The human β₁-adrenoceptor (β₁AR) is a G-protein-coupled receptor (GPCR) involved in sympathetic system regulation through agonist-induced activation. The conserved CWXP-motif in helix 6 (rotamer toggle switch) is one of the most important activation switches in Class A GPCRs. In order to investigate how the agonist binding disturbs this switch, we carried out molecular dynamics simulations of a hβ₁AR model in the apo and R-noradrenaline-bound forms. The results show that the agonist binding changes the β₁-angle distribution of Cys336, Trp337 and Phe341 residues and increases the helix 6 bending. Overall, we provide a functional hβ₁AR model, showing how the rotamer toggle switch mechanism works at atomic level.     

Linking agonist binding to histamine H1 receptor activation

G protein-coupled receptors (GPCRs) constitute a large and functionally diverse family of transmembrane proteins. They are fundamental in the transfer of extracellular stimuli to intracellular signaling pathways and are among the most targeted proteins in drug discovery. The detailed molecular mechanism for agonist-induced activation of rhodopsin-like GPCRs has not yet been described. Using a combination of site-directed mutagenesis and molecular modeling, we characterized important steps in the activation of the human histamine H1 receptor. Both Ser3.36 and Asn7.45 are important links between histamine binding and previously proposed conformational changes in helices 6 and 7. Ser3.36 acts as a rotamer toggle switch that, upon agonist binding, initiates the activation of the receptor through Asn7.45. The proposed transduction involves specific residues that are conserved among rhodopsin-like GPCRs.     

Ligand binding-dependent limited proteolysis of the atrial natriuretic peptide receptor: juxtamembrane hinge structure essential for transmembrane signal transduction

The atrial natriuretic peptide (ANP) receptor is a 130-kDa transmembrane protein containing an extracellular ANP-binding domain, a single transmembrane sequence, an intracellular kinase-homologous domain, and a guanylate cyclase (GCase) domain. We observed that the receptor, when bound with ANP, was rapidly cleaved by endogenous or exogenously added protease to yield a 65-kDa ANP-binding fragment. No cleavage occurred without bound ANP. This ligand-induced cleavage abolished GCase activation by ANP. Cleavage occurred in an extracellular, juxtamembrane region containing six closely spaced Pro residues and a disulfide bond. Such structural features are shared among the A-type and B-type ANP receptors but not by ANP clearance receptors. The potential role of the hinge structure was examined by mutagenesis experiments. Mutation of Pro(417), but not other Pro residues, to Ala abolished GCase activation by ANP. Elimination of the disulfide bond by Cys to Ser mutations yielded a constitutively active receptor. Pro(417), and Cys(423) and Cys(432) forming the disulfide bond are strictly conserved among GCase-coupled receptors, while other residues are largely variable. The conserved Pro(417) and the disulfide bond may represent a consensus signaling motif in the juxtamembrane hinge structure that undergoes a marked conformational change upon ligand binding and apparently mediates transmembrane signal transduction.     

A Ser residue influences the structure and stability of a Pro-kinked transmembrane helix dimer

When localized adjacent to a Pro-kink, Thr and Ser residues can form hydrogen bonds between their polar hydroxyl group and a backbone carbonyl oxygen and thereby modulate the actual bending angle of a distorted transmembrane α-helix. We have used the homo-dimeric transmembrane cytochrome b(559)' to analyze the potential role of a highly conserved Ser residue for assembly and stabilization of transmembrane proteins. Mutation of the conserved Ser residue to Ala resulted in altered heme binding properties and in increased stability of the holo-protein, most likely by tolerating subtle structural rearrangements upon heme binding. The results suggest a crucial impact of an intrahelical Ser hydrogen bond in defining the structure of a Pro-kinked transmembrane helix dimer.     

Activation of the beta 2-adrenergic receptor involves disruption of an ionic lock between the cytoplasmic ends of transmembrane segments 3 and 6

The movements of transmembrane segments (TMs) 3 and 6 at the cytoplasmic side of the membrane play an important role in the activation of G-protein-coupled receptors. Here we provide evidence for the existence of an ionic lock that constrains the relative mobility of the cytoplasmic ends of TM3 and TM6 in the inactive state of the beta(2)-adrenergic receptor. We propose that the highly conserved Arg-131(3.50) at the cytoplasmic end of TM3 interacts both with the adjacent Asp-130(3.49) and with Glu-268(6.30) at the cytoplasmic end of TM6. Such a network of ionic interactions has now been directly supported by the high-resolution structure of the inactive state of rhodopsin. We hypothesized that the network of interactions would serve to constrain the receptor in the inactive state, and the release of this ionic lock could be a key step in receptor activation. To test this hypothesis, we made charge-neutralizing mutations of Glu-268(6.30) and of Asp-130(3.49) in the beta(2)-adrenergic receptor. Alone and in combination, we observed a significant increase in basal and pindolol-stimulated cAMP accumulation in COS-7 cells transiently transfected with the mutant receptors. Moreover, based on the increased accessibility of Cys-285(6.47) in TM6, we provide evidence for a conformational rearrangement of TM6 that is highly correlated with the extent of constitutive activity of the different mutants. The present experimental data together with the recent high-resolution structure of rhodopsin suggest that ionic interactions between Asp/Glu(3.49), Arg(3.50), and Glu(6.30) may constitute a common switch governing the activation of many rhodopsin-like G-protein-coupled receptors.     

Serine and threonine residues bend alpha-helices in the chi(1) = g(-) conformation

The relationship between the Ser, Thr, and Cys side-chain conformation (chi(1) = g(-), t, g(+)) and the main-chain conformation (phi and psi angles) has been studied in a selection of protein structures that contain alpha-helices. The statistical results show that the g(-) conformation of both Ser and Thr residues decreases their phi angles and increases their psi angles relative to Ala, used as a control. The additional hydrogen bond formed between the O(gamma) atom of Ser and Thr and the i-3 or i-4 peptide carbonyl oxygen induces or stabilizes a bending angle in the helix 3-4 degrees larger than for Ala. This is of particular significance for membrane proteins. Incorporation of this small bending angle in the transmembrane alpha-helix at one side of the cell membrane results in a significant displacement of the residues located at the other side of the membrane. We hypothesize that local alterations of the rotamer configurations of these Ser and Thr residues may result in significant conformational changes across transmembrane helices, and thus participate in the molecular mechanisms underlying transmembrane signaling. This finding has provided the structural basis to understand the experimentally observed influence of Ser residues on the conformational equilibrium between inactive and active states of the receptor, in the neurotransmitter subfamily of G protein-coupled receptors.     

Metal ion site engineering indicates a global toggle switch model for seven-transmembrane receptor activation

Much evidence indicates that, during activation of seven-transmembrane (7TM) receptors, the intracellular segments of the transmembrane helices (TMs) move apart with large amplitude, rigid body movements of especially TM-VI and TM-VII. In this study, AspIII:08 (Asp113), the anchor point for monoamine binding in TM-III, was used as the starting point to engineer activating metal ion sites between the extracellular segments of the beta2-adrenergic receptor. Cu(II) and Zn(II) alone and in complex with aromatic chelators acted as potent (EC50 decreased to 0.5 microm) and efficacious agonists in sites constructed between positions III:08 (Asp or His), VI:16 (preferentially Cys), and/or VII:06 (preferentially Cys). In molecular models built over the backbone conformation of the inactive rhodopsin structure, the heavy atoms that coordinate the metal ion were located too far away from each other to form high affinity metal ion sites in both the bidentate and potential tridentate settings. This indicates that the residues involved in the main ligand-binding pocket will have to move closer to each other during receptor activation. On the basis of the distance constraints from these activating metal ion sites, we propose a global toggle switch mechanism for 7TM receptor activation in which inward movement of the extracellular segments of especially TM-VI and, to some extent, TM-VII is coupled to the well established outward movement of the intracellular segments of these helices. We suggest that the pivots for these vertical seesaw movements are the highly conserved proline bends of the involved helices.     

Dopamine Receptor Microdomains Involved in Molecular Recognition and the Regulation of Drug Affinity and Function

A cationic protonatable amine moiety on dopaminergic ligands forms a high affinity reinforced ionic bond with an anionic aspartic acid at position 3.32 of dopamine receptors. When present, catechol hydroxyls of the ligands form hydrogen bonds with serines at position 5.42, 5.43, and 5.46, and this network of hydrogen bonds serves to orient ligands in the binding-site crevice and increase their binding affinity. A steric clash between aromatic moieties of the ligands and aromatic amino acids of the receptor (e.g., H6.55, F6.52 or F6.51 and W6.48) is likely to be propagated in domino-like fashion along the length of TM6, which is believed to trigger activation of the receptor. Specifically, it is the change in the conformation of W6.48 from an orientation perpendicular to the plane of the lipid membrane to one that is parallel that is believed to result in activation. Molecular determinants that mediate the D4/D2-selectivity of many extremely D4-selective 1,4-DAP ligands, include a nonconserved cluster of bulky amino acids at the TM2/TM3 interface (positions 2.61, 3.28 and 3.29).     

Dopamine receptor microdomains involved in molecular recognition and the regulation of drug affinity and function

A cationic protonatable amine moiety on dopaminergic ligands forms a high affinity reinforced ionic bond with an anionic aspartic acid at position 3.32 of dopamine receptors. When present, catechol hydroxyls of the ligands form hydrogen bonds with serines at position 5.42, 5.43, and 5.46, and this network of hydrogen bonds serves to orient ligands in the binding-site crevice and increase their binding affinity. A steric clash between aromatic moieties of the ligands and aromatic amino acids of the receptor (e.g., H6.55, F6.52 or F6.51 and W6.48) is likely to be propagated in domino-like fashion along the length of TM6, which is believed to trigger activation of the receptor. Specifically, it is the change in the conformation of W6.48 from an orientation perpendicular to the plane of the lipid membrane to one that is parallel that is believed to result in activation. Molecular determinants that mediate the D4/D2-selectivity of many extremely D4-selective 1,4-DAP ligands, include a nonconserved cluster of bulky amino acids at the TM2/TM3 interface (positions 2.61, 3.28 and 3.29).     

Ser and Thr Residues Modulate the Conformation of Pro-Kinked Transmembrane α-Helices

Functionally required conformational plasticity of transmembrane proteins implies that specific structural motifs have been integrated in transmembrane helices. Surveying a database of transmembrane helices and the large family of G-protein coupled receptors we identified a series of overrepresented motifs associating Pro with either Ser or Thr. Thus, we have studied the conformation of Pro-kinked transmembrane helices containing Ser or Thr residues, in both g+ and g− rotamers, by molecular dynamics simulations in a hydrophobic environment. Analysis of the simulations shows that Ser or Thr can significantly modulate the deformation of the Pro. A series of motifs, such as (S/T)P and (S/T)AP in the g+ rotamer and the TAP and PAA(S/T) motifs in the g− rotamer, induce an increase in bending angle of the helix compared to a standard Pro-kink, apparently due to the additional hydrogen bond formed between the side chain of Ser/Thr and the backbone carbonyl oxygen. In contrast, (S/T)AAP and PA(S/T) motifs, in both g+ and g−, and PAA(S/T) in g+ rotamers decrease the bending angle of the helix by either reducing the steric clash between the pyrrolidine ring of Pro and the helical backbone, or by adding a constrain in the form of a hydrogen bond in the curved-in face of the helix. Together with a number of available experimental data, our results strongly suggest that association of Ser and Thr with Pro is commonly used in transmembrane helices to accommodate the structural needs of specific functions.     

The role of transmembrane segment TM3 in the xanthine permease XanQ of Escherichia coli

The xanthine permease XanQ of Escherichia coli is used as a study prototype for function-structure analysis of the ubiquitous nucleobase-ascorbate transporter (NAT/NCS2) family. Our previous mutagenesis study of polar residues of XanQ has shown that Asn-93 at the middle of putative TM3 is a determinant of substrate affinity and specificity. To study the role of TM3 in detail we employed Cys-scanning mutagenesis. Using a functional mutant devoid of Cys residues (C-less), each amino acid residue in sequence 79-107 (YGIVGSGLLSIQSVNFSFVTVMIALGSSM) including TM3 (underlined) and flanking sequences was replaced individually with Cys. Of 29 single-Cys mutants, 20 accumulate xanthine to 40-110% of the steady state observed with C-less, six (S88C, F94C, A102C, G104C, S106C) accumulate to low levels (10-30%) and three (G83C, G85C, N93C) are inactive. Extensive mutagenesis reveals that Gly-83 and, to a lesser extent, Gly-85, are crucial for expression in the membrane. Replacements of Asn-93 disrupt affinity (Thr) or permit recognition of 8-methylxanthine which is not a wild-type ligand (Ala, Ser, Asp) and utilization of uric acid which is not a wild-type substrate (Ala, Ser). Replacements of Phe-94 impair affinity for 2-thio and 6-thioxanthine (Tyr) or 3-methylxanthine (Ile). Single-Cys mutants S84C, L86C, L87C, and S95C are highly sensitive to inactivation by N-ethylmaleimide. Our data reveal that key residues of TM3 cluster in two conserved sequence motifs, (83)GSGLL(87) and (93)NFS(95), and highlight the importance of Asn-93 and Phe-94 in substrate recognition and specificity; these findings are supported by structural modeling on the recently described x-ray structure of the uracil-transporting homolog UraA.     

Identification of an efficacy switch region in the ghrelin receptor responsible for interchange between agonism and inverse agonism

The carboxyamidated wFwLL peptide was used as a core ligand to probe the structural basis for agonism versus inverse agonism in the constitutively active ghrelin receptor. In the ligand, an efficacy switch could be built at the N terminus, as exemplified by AwFwLL, which functioned as a high potency agonist, whereas KwFwLL was an equally high potency inverse agonist. The wFw-containing peptides, agonists as well as inverse agonists, were affected by receptor mutations covering the whole main ligand-binding pocket with key interaction sites being an aromatic cluster in transmembrane (TM)-VI and -VII and residues on the opposing face of TM-III. Gain-of-function in respect of either increased agonist or inverse agonist potency or swap between high potency versions of these properties was obtained by substitutions at a number of positions covering a broad area of the binding pocket on TM-III, -IV, and -V. However, in particular, space-generating substitutions at position III:04 shifted the efficacy of the ligands from inverse agonism toward agonism, whereas similar substitutions at position III: 08, one helical turn below, shifted the efficacy from agonism toward inverse agonism. It is suggested that the relative position of the ligand in the binding pocket between this "efficacy shift region" on TM-III and the opposing aromatic cluster on TM-VI and TM-VII leads either to agonism, i.e. in a superficial binding mode, or it leads to inverse agonism, i.e. in a more profound binding mode. This relationship between different binding modes and opposite efficacy is in accordance with the Global Toggle Switch model for 7TM receptor activation.     

Structural mechanisms of constitutive activation in the C5a receptors with mutations in the extracellular loops: molecular modeling study

Previously we demonstrated by random saturation mutagenesis a set of mutations in the extracellular (EC) loops that constitutively activate the C5a receptor (C5aR) (Klco et al., Nat Struct Mol Biol 2005;12:320-326; Klco et al., J Biol Chem 2006;281:12010-12019). In this study, molecular modeling revealed possible conformations for the extracellular loops of the C5a receptors with mutations in the EC2 loop or in the EC3 loop. Comparison of low-energy conformations of the EC loops defined two distinct clusters of conformations typical either for strongly constitutively active mutants of C5aR (the CAM cluster) or for nonconstitutively active mutants (the non-CAM cluster). In the CAM cluster, the EC3 loop was turned towards the transmembrane (TM) helical bundle and more closely interacted with EC2 than in the non-CAM cluster. This suggested a structural mechanism of constitutive activity where EC3 contacts EC2 leading to EC2 interactions with helix TM3, thus triggering movement of TM7 towards TM2 and TM3. The movement initiates rearrangement of the system of hydrogen bonds between TM2, TM3 and TM7 including formation of the hydrogen bond between the side chains of D82(2.50) in TM2 and N296(7.49) in TM7, which is crucial for formation of the activated states of the C5a receptors (Nikiforovich et al., Proteins: Struct Funct Gene 2011;79:787-802). Since the relative large length of EC3 in C5aR (13 residues) is comparable with those in many other members of rhodopsin family of GPCRs (13-19 residues), our findings might reflect general mechanisms of receptor constitutive activation. The very recent X-ray structure of the agonist-induced constitutively active mutant of rhodopsin (Standfuss et al., Nature 2011;471:656-660) is discussed in view of our modeling results.     

Rho-kinase phosphorylates eNOS at threonine 495 in endothelial cells

Endothelial nitric oxide synthase (eNOS) produces nitric oxide (NO), which is involved in various physiological functions of the cardiovascular system. eNOS is activated by dephosphorylation at Thr495 and phosphorylation at Ser1177. Inhibition of Rho-kinase, an effector of the small GTPase RhoA, leads to activation of Akt/PKB, which phosphorylates eNOS at Ser1177 and thereby promotes NO production. However, little is known about the effects of Rho-kinase on phosphorylation of Thr495. We here found that the constitutively active form of Rho-kinase phosphorylated eNOS at Thr495 in vitro. Expression of the constitutively active form of RhoA or Rho-kinase increased this phosphorylation in COS-7 cells. Addition of thrombin to cultured human umbilical vein endothelial cells induced phosphorylation of eNOS at Thr495. Treatment with Y27632, a Rho-kinase inhibitor, suppressed thrombin-induced phosphorylation at Thr495. These results indicate that Rho-kinase can directly phosphorylate eNOS at Thr495 to suppress NO production in endothelium.     

Identification of novel phosphorylation sites required for activation of MAPKAP kinase-2.

MAP kinase-activated protein (MAPKAP) kinase-2 is activated in vivo by reactivating kinase (RK), a MAP kinase (MAPK) homologue stimulated by cytokines and cellular stresses. Here we show that in vitro RK phosphorylates human GST-MAPKAP kinase-2 at Thr25 in the proline-rich N-terminal region Thr222 and Ser272 in the catalytic domain and Thr334 in the C-terminal domain. Using novel methodology we demonstrate that activation of MAPKAP kinase-2 requires the phosphorylation of any two of the three residues Thr222, Ser272 and Thr334. Ser9, Thr25, Thr222, Ser272, Thr334 and Thr338 became 32P-labelled in stressed KB cells and labelling was prevented by the specific RK inhibitor SB 203580, establishing that RK phosphorylates Thr25, Thr222, Ser272 and Thr334 in vivo. The 32P-labelling of Thr338 is likely to result from autophosphorylation. GST-MAPKAP kinase-2 lacking the N-terminal domain was inactive, but activated fully when phosphorylated at Thr222, Ser272 and Thr334 by p42 MAPK or RK. In contrast, full-length GST-MAPKAP kinase-2 was phosphorylated at Thr25 (and not activated) by p42 MAPK, suggesting a role for the N-terminal domain in specifying activation by RK in vivo. The mutant Asp222/Asp334 was 20% as active as phosphorylated MAPKAP kinase-2, and this constitutively active form may be useful for studying its physiological roles.     

Structural determinants for ligand-receptor conformational selection in a peptide G protein-coupled receptor

G protein coupled receptors (GPCRs) modulate the majority of physiological processes through specific intermolecular interactions with structurally diverse ligands and activation of differential intracellular signaling. A key issue yet to be resolved is how GPCRs developed selectivity and diversity of ligand binding and intracellular signaling during evolution. We have explored the structural basis of selectivity of naturally occurring gonadotropin-releasing hormones (GnRHs) from different species in the single functional human GnRH receptor. We found that the highly variable amino acids in position 8 of the naturally occurring isoforms of GnRH play a discriminating role in selecting receptor conformational states. The human GnRH receptor has a higher affinity for the cognate GnRH I but a lower affinity for GnRH II and GnRHs from other species possessing substitutions for Arg(8). The latter were partial agonists in the human GnRH receptor. Mutation of Asn(7.45) in transmembrane domain (TM) 7 had no effect on GnRH I affinity but specifically increased affinity for other GnRHs and converted them to full agonists. Using molecular modeling and site-directed mutagenesis, we demonstrated that the highly conserved Asn(7.45) makes intramolecular interactions with a highly conserved Cys(6.47) in TM 6, suggesting that disruption of this intramolecular interaction induces a receptor conformational change which allosterically alters ligand specific binding sites and changes ligand selectivity and signaling efficacy. These results reveal GnRH ligand and receptor structural elements for conformational selection, and support co-evolution of GnRH ligand and receptor conformations.     

Correlations in palmitoylation and multiple phosphorylation of rat bradykinin B2 receptor in Chinese hamster ovary cells

Rat bradykinin B2 receptor from unstimulated Chinese hamster ovary cells transfected with the corresponding cDNA has been isolated, and subsequent mass spectrometric analysis of multiple phosphorylated species and of the palmitoylation attachment site is described. Bradykinin B2 receptor was isolated on oligo(dT)-cellulose using N-(epsilon-maleimidocaproyloxy)succinimide-Met-Lys-bradykinin coupled to a protected (dA)30-mer. This allowed a one-step isolation of the receptor on an oligo(dT)-cellulose column via variation solely of salt concentration. After enzymatic in-gel digestion, matrix-assisted laser desorption ionization and electrospray ion trap mass spectrometric analysis of the isolated rat bradykinin B2 receptor showed phosphorylation at Ser365, Ser371, Ser378, Ser380, and Thr374. Further phosphorylation at Tyr352 and Tyr161 was observed. Rat bradykinin receptor B2 receptor is also palmitoylated at Cys356. All of the phosphorylation sites except for Tyr161 cluster at the carboxyl-terminal domain of the receptor located on the cytoplasmic face of the cell membrane. Surprisingly, many of the post-translational modifications were shown by MSn mass spectroscopic analysis to be correlated pairwise, e.g. diphosphorylation at Ser365 and Ser371, at Ser378 and Ser380, and at Thr374 and Ser380 as well as mutually exclusive phosphorylation at Tyr352 and palmitoylation at Cys356. The last correlation may be involved in a receptor internalization motif. Pairwise correlations and mutual exclusion of phosphorylation and palmitoylation suggest critical roles of multiple post-translational modifications for the regulation of activity, coupling to intracelluar signaling pathways, and/or sequestration of the bradykinin receptor.     

Structure-function relationships of the human bitter taste receptor hTAS2R1: insights from molecular modeling studies

Bitter taste receptors (T2Rs) belong to G-protein-coupled receptors (GPCRs). Despite extensive studies, the precise mechanisms of GPCR activation are still poorly understood. In this study, the models of the human bitter taste receptor hTAS2R1 alone and in complex with various ligands were constructed on the basis of template-based modeling and molecular docking. Then these models were subjected to all-atom molecular dynamics (MD) simulations in explicit lipid bilayers. The binding pocket of hTAS2R1 is mainly formed by transmembrane helix (TM) III, TM V, TM VI, and TM VII. Most of the residues contributing to ligand binding are positionally conserved comparing with other hTAS2Rs. By comparing the final conformations obtained by extensive MD simulations, we identified the changes in the transmembrane helices and the intra- and extracellular loops, which were expected to initiate the activation of the receptor. The intracellular loop II (ICL2) and TM III were found to play prominent roles in the process of activation. We proposed that a set of interactions between the aromatic Phe115 in the middle of ICL2 and three residues (Tyr103, Lys106, and Val107) at the cytoplasmic end of TM III may serve as a conformational switch of hTAS2R1 activation. All of the residues involved in the switch are highly conserved among T2Rs. This indicates that the control switch we proposed may be universal in T2Rs. Besides, our results also suggest that the formation of a short helical segment in ICL2 may be necessary for the activation of hTAS2R1.