Agonist-induced conformational changes in thyrotropin-releasing hormone receptor type I: disulfide cross-linking and molecular modeling approaches

The Antagonist pGlu-βGlu-Pro-NH2 Binds to an Allosteric Site of the Thyrotropin-Releasing Hormone Receptor

After we identified pGlu-βGlu-Pro-NH2 as the first functional antagonist of the cholinergic central actions of the thyrotropin-releasing hormone (TRH, pGlu-His-Pro-NH2), we became interested in finding the receptor-associated mechanism responsible for this antagonism. By utilizing a human TRH receptor (hTRH-R) homology model, we first refined the active binding site within the transmembrane bundle of this receptor to enhance TRH's binding affinity. However, this binding site did not accommodate the TRH antagonist. This directed us to consider a potential allosteric binding site in the extracellular domain (ECD). Searches for ECD binding pockets prompted the remodeling of the extracellular loops and the N-terminus. We found that different trajectories of ECDs produced novel binding cavities that were then systematically probed with TRH, as well as its antagonist. This led us to establish not only a surface-recognition binding site for TRH, but also an allosteric site that exhibited a selective and high-affinity binding for pGlu-βGlu-Pro-NH2. The allosteric binding of this TRH antagonist is more robust than TRH's binding to its own active site. The findings reported here may shed light on the mechanisms and the multimodal roles by which the ECD of a TRH receptor is involved in agonist and/or antagonist actions.

Structural features of the G-protein/GPCR interactions

BACKGROUND

The details of the functional interaction between G proteins and the G protein coupled receptors (GPCRs) have long been subjected to extensive investigations with structural and functional assays and a large number of computational studies.

SCOPE OF REVIEW

The nature and sites of interaction in the G-protein/GPCR complexes, and the specificities of these interactions selecting coupling partners among the large number of families of GPCRs and G protein forms, are still poorly defined.

MAJOR CONCLUSIONS

Many of the contact sites between the two proteins in specific complexes have been identified, but the three dimensional molecular architecture of a receptor-Gα interface is only known for one pair. Consequently, many fundamental questions regarding this macromolecular assembly and its mechanism remain unanswered.

GENERAL SIGNIFICANCE

In the context of current structural data we review the structural details of the interfaces and recognition sites in complexes of sub-family A GPCRs with cognate G-proteins, with special emphasis on the consequences of activation on GPCR structure, the prevalence of preassembled GPCR/G-protein complexes, the key structural determinants for selective coupling and the possible involvement of GPCR oligomerization in this process.

PTH and PTH antagonist induce different conformational changes in the PTHR1 receptor

Interaction of ligands with their specific receptors is accompanied by conformational shifts culminating in receptor activation and expression of hormonal activity. Using an engineered disulfide bond formation strategy, we characterized the relative conformational changes taking place within the PTH type 1 receptor (PTHR1) at the interface of transmembrane (TM)5 and TM6 on binding the PTH agonist, PTH(1-34), compared with the antagonist PTH(7-34). Cysteines were singly incorporated into a portion of the extracellular-facing region of TM5 (365-370), while simultaneously a second cysteine was introduced at position 420, 423, or 425 at the extracellular end of TM6, leading to a total of 18 double cysteine-containing PTHR1 mutants. All mutants, except P366C/V423C and P366C/M425C, were expressed in the cell membrane preparations. In the presence of agonist, H420C and M425C in TM6 formed disulfide bonds with all and with most, respectively, of the substituted cysteines incorporated in TM5. In contrast to the conformational shift induced (or stabilized) by agonist in activating the receptor, antagonist binding produced no detectable change from the basal (inactive) conformation of PTHR1. Our studies provide physicochemical evidence that the extracellular-facing ligand binding regions of receptor, TM5 and TM6, are dynamic and move relative to each other on ligand binding. The distinct differences in receptor conformation induced (or stabilized) by agonist PTH(1-34) compared with antagonist PTH(7-34) begin to provide insight into the early events in and mechanism of PTHR1 activation.

Two arginine-glutamate ionic locks near the extracellular surface of FFAR1 gate receptor activation

Activation of a number of class A G protein-coupled receptors (GPCRs) is thought to involve two molecular switches, a rotamer toggle switch within the transmembrane domain and an ionic lock at the cytoplasmic surface of the receptor; however, the mechanism by which agonist binding changes these molecular interactions is not understood. Importantly, 80% of GPCRs including free fatty acid receptor 1 (FFAR1) lack the complement of amino acid residues implicated in either or both of these two switches; the mechanism of activation of these GPCRs is therefore less clear. By homology modeling, we identified two Glu residues (Glu-145 and Glu-172) in the second extracellular loop of FFAR1 that form putative interactions individually with two transmembrane Arg residues (Arg-183(5.39) and Arg-258(7.35)) to create two ionic locks. Molecular dynamics simulations showed that binding of agonists to FFAR1 leads to breakage of these Glu-Arg interactions. In mutagenesis experiments, breakage of these two putative interactions by substituting Ala for Glu-145 and Glu-172 caused constitutive receptor activation. Our results therefore reveal a molecular switch for receptor activation present on the extracellular surface of FFAR1 that is broken by agonist binding. Similar ionic locks between the transmembrane domains and the extracellular loops may constitute a mechanism common to other class A GPCRs also.

Conformational changes involved in G-protein-coupled-receptor activation

Little is known about the nature of the conformational changes that convert G-protein-coupled receptors (GPCRs), which bind diffusible ligands, from their resting into their active states. To gain structural insight into this process, various laboratories have used disulfide cross-linking strategies involving cysteine-substituted mutant GPCRs. Several recent disulfide cross-linking studies using the M(3) muscarinic acetylcholine receptor as a model system have led to novel insights into the conformational changes associated with the activation of this prototypical class I GPCR. These structural changes are predicted to involve multiple receptor regions, primarily distinct segments of transmembrane helices III, VI and VII and helix 8. Given the high degree of structural homology found among most GPCRs, it is likely that these findings will be of considerable general relevance. A better understanding of the molecular mechanisms underlying GPCR activation might lead to novel strategies aimed at modulating GPCR function for therapeutic purposes.

On the applicability of GPCR homology models to computer-aided drug discovery: a comparison between in silico and crystal structures of the beta2-adrenergic receptor

The publication of the crystal structure of the beta2-adrenergic receptor (beta2-AR) proved that G protein-coupled receptors (GPCRs) share a structurally conserved rhodopsin-like 7TM core. Here, to probe to which extent realistic GPCR structures can be recreated through modeling, carazolol was docked at two rhodopsin-based homology models of the human beta 2-AR. The first featured a rhodopsin-like second extracellular loop, which interfered with ligand docking and with the orientation of several residues in the binding pocket. The second featured a second extracellular loop built completely de novo, which afforded a more accurate model of the binding pocket and a better docking of the ligand. Furthermore, incorporating available biochemical and computational data to the model by correcting the conformation of a single residue lining the binding pocket --Phe290(6.52)--, resulted in significantly improved docking poses. These results support the applicability of GPCR modeling to the design of site-directed mutagenesis experiments and to drug discovery.

Understanding the structural and functional differences between mouse thyrotropin-releasing hormone receptors 1 and 2

Multiple computational methods have been employed in a comparative study of thyrotropin-releasing hormone receptors 1 and 2 (TRH-R1 and TRH-R2) to explore the structural bases for the different functional properties of these G protein-coupled receptors. Three-dimensional models of both murine TRH receptors have been built and optimized by means of homology modeling based on the crystal structure of bovine rhodopsin, molecular dynamics simulations, and energy minimizations in a membrane-aqueous environment. The comparison between the two models showed a correlation between the higher flexibility and higher basal activity of TRH-R2 versus the lesser flexibility and lower basal activity of TRH-R1 and supported the involvement of the highly conserved W6.48 in the signaling process. A correlation between the level of basal activity and conformational changes of TM5 was detected also. Comparison between models of the wild type receptors and their W6.48A mutants, which have reversed basal activities compared with their respective wild types, further supported these correlations. A flexible molecular docking procedure revealed that TRH establishes a direct interaction with W6.48 in TRH-R2 but not in TRH-R1. We designed and performed new mutagenesis experiments that strongly supported these observations.

Conformational change of the transmembrane helices II and IV of metabotropic glutamate receptor involved in G protein activation

G protein-coupled receptors are classified into several families on the basis of their amino acid sequences and the members of the same family exhibit sequence similarity but those of different families do not. In family 1 GPCRs such as rhodopsin and adrenergic receptor, extensive studies have revealed the stimulus-dependent conformational change of the receptor: the rearrangement of transmembrane helices III and VI is essential for G protein activation. In contrast, in family 3 GPCRs such as metabotropic glutamate receptor (mGluR), the inter-protomer relocation upon ligand binding has been observed but there is much less information about the structural changes of the transmsmbrane helices and the cytoplasmic domains. Here we identified constitutively active mutation sites at the cytoplasmic borders of helices II and IV of mGluR8 and successfully inhibited the G protein activation ability by engineering disulfide cross-linking between these cytoplasmic regions. The analysis of all possible single substitution mutants of these residues revealed that some steric interactions around these sites would be important to keep the receptor protein inactive. These results provided the model that the conformational changes at the cytoplasmic ends of helices II and IV of mGluR are involved in the efficient G protein coupling.

Activation of G protein-coupled receptors: beyond two-state models and tertiary conformational changes

Transformation of G protein-coupled receptors (GPCRs) from a quiescent to an active state initiates signal transduction. All GPCRs share a common architecture comprising seven transmembrane-spanning alpha-helices, which accommodates signal propagation from a diverse repertoire of external stimuli across biological membranes to a heterotrimeric G protein. Signal propagation through the transmembrane helices likely involves mechanistic features common to all GPCRs. The structure of the light receptor rhodopsin may serve as a prototype for the transmembrane architecture of GPCRs. Early biochemical, biophysical, and pharmacological studies led to the conceptualization of receptor activation based on the context of two-state equilibrium models and conformational changes in protein structure. More recent studies indicate a need to move beyond these classical paradigms and to consider additional aspects of the molecular character of GPCRs, such as the oligomerization and dynamics of the receptor.

Conformational changes in the parathyroid hormone receptor associated with activation by agonist

Binding of hormones to their cognate G protein-coupled receptors (GPCRs) induces conformational shifts within the receptor based on evidence from a few hormone-receptor systems. Employing an engineered disulfide bond formation strategy and guided by a previously established model of the PTH-PTH receptor (PTHR)1 bimolecular complex, we set out to document and characterize the nature of agonist-induced changes in this family B GPCR. A mutant PTHR1 was generated which incorporates a Factor Xa cleavage site in the third intracellular loop. Treatment with Factor Xa fragments the receptor. However, if a new disulfide bond was formed before exposure to the enzyme, the fragments remain held together. A set of double cysteine-containing mutants were designed to probe the internal relative movements of transmembrane (TM) helices 2 and TM7. PTH enhanced formation of disulfide bonds in the K240C/F447C and A242C/F447C mutants. For the F238C/F447C mutant, a disulfide bond is formed in the basal state, but is disrupted by interaction with PTH. For the D241C/F447C PTHR1 construct, no disulfide bond formation was observed in either the basal or hormone-bound state. These findings demonstrate that the conformation of PTHR1 is altered from the basal state when PTH is bound. Novel information regarding spatial proximities between TM2 and TM7 of PTHR1 and the nature of relative movements between the two transmembrane regions was revealed. The data confirm and extend the experimentally derived model of the PTH-PTHR1 complex and provide insights at a new level of detail into the early events in PTHR1 activation by PTH.

Heterotrimeric G protein activation by G-protein-coupled receptors

Heterotrimeric G proteins have a crucial role as molecular switches in signal transduction pathways mediated by G-protein-coupled receptors. Extracellular stimuli activate these receptors, which then catalyse GTP-GDP exchange on the G protein alpha-subunit. The complex series of interactions and conformational changes that connect agonist binding to G protein activation raise various interesting questions about the structure, biomechanics, kinetics and specificity of signal transduction across the plasma membrane.

Thyrotropin-releasing hormone and its receptors — A hypothesis for binding and receptor activation

Thyrotropin-releasing hormone (TRH), a tripeptide, exerts its biological effects through stimulation of cell-surface receptors, TRH-R, belonging to the superfamily of G protein-coupled receptors (GPCR). Because of the intermediate size of TRH, it is smaller than polypeptide ligands that interact at GPCR ectodomains and larger than biogenic amines, which interact within GPCR transmembrane domains (TMD), the TRH/TRH-R complex probably shares properties of these 2 extremes, representing a unique system to study GPCR/ligand interactions. In this review, we summarize the current knowledge of the structure-activity relationships in the TRH/TRH-R system. Based on experimental data and the structural information acquired from computer simulations, we formulate a working hypothesis to describe the molecular events underlying the processes of TRH binding and TRH-R activation. This hypothesis represents a starting point for understanding the biology of the TRH/TRH-R system on a molecular level and provides a basis for potential design of new potent and selective modulators of TRH-R's activity.

Contribution of presenilin transmembrane domains 6 and 7 to a water-containing cavity in the gamma-secretase complex

Gamma-secretase is a multiprotein complex responsible for the intramembranous cleavage of the amyloid precursor protein and other type I transmembrane proteins. Mutations in Presenilin, the catalytic core of this complex, cause Alzheimer disease. Little is known about the structure of the protein and even less about the catalytic mechanism, which involves proteolytic cleavage in the hydrophobic environment of the cell membrane. It is basically unclear how water, needed to perform hydrolysis, is provided to this reaction. Presenilin transmembrane domains 6 and 7 seem critical in this regard, as each bears a critical aspartate contributing to catalytic activity. Current models imply that both aspartyl groups should closely oppose each other and have access to water. This is, however, still to be experimentally verified. Here, we have performed cysteine-scanning mutagenesis of both domains and have demonstrated that several of the introduced residues are exposed to water, providing experimental evidence for the existence of a water-filled cavity in the catalytic core of Presenilin. In addition, we have demonstrated that the two aspartates reside within this cavity and are opposed to each other in the native complex. We have also identified the conserved tyrosine 389 as a critical partner in the catalytic mechanism. Several additional amino acid substitutions affect differentially the processing of gamma-secretase substrates, implying that they contribute to enzyme specificity. Our data suggest the possibility that more selective gamma-secretase inhibitors could be designed.

RADIOLYTIC PROTEIN FOOTPRINTING WITH MASS SPECTROMETRY TO PROBE THE STRUCTURE OF MACROMOLECULAR COMPLEXES

Structural proteomics approaches using mass spectrometry are increasingly used in biology to examine the composition and structure of macromolecules. Hydroxyl radical-mediated protein footprinting using mass spectrometry has recently been developed to define structure, assembly, and conformational changes of macromolecules in solution based on measurements of reactivity of amino acid side chain groups with covalent modification reagents. Accurate measurements of side chain reactivity are achieved using quantitative liquid-chromatography-coupled mass spectrometry, whereas the side chain modification sites are identified using tandem mass spectrometry. In addition, the use of footprinting data in conjunction with computational modeling approaches is a powerful new method for testing and refining structural models of macromolecules and their complexes. In this review, we discuss the basic chemistry of hydroxyl radical reactions with peptides and proteins, highlight various approaches to map protein structure using radical oxidation methods, and describe state-of-the-art approaches to combine computational and footprinting data.

Structure modeling of all identified G protein-coupled receptors in the human genome

G protein–coupled receptors (GPCRs), encoded by about 5% of human genes, comprise the largest family of integral membrane proteins and act as cell surface receptors responsible for the transduction of endogenous signal into a cellular response. Although tertiary structural information is crucial for function annotation and drug design, there are few experimentally determined GPCR structures. To address this issue, we employ the recently developed threading assembly refinement (TASSER) method to generate structure predictions for all 907 putative GPCRs in the human genome. Unlike traditional homology modeling approaches, TASSER modeling does not require solved homologous template structures; moreover, it often refines the structures closer to native. These features are essential for the comprehensive modeling of all human GPCRs when close homologous templates are absent. Based on a benchmarked confidence score, approximately 820 predicted models should have the correct folds. The majority of GPCR models share the characteristic seven-transmembrane helix topology, but 45 ORFs are predicted to have different structures. This is due to GPCR fragments that are predominantly from extracellular or intracellular domains as well as database annotation errors. Our preliminary validation includes the automated modeling of bovine rhodopsin, the only solved GPCR in the Protein Data Bank. With homologous templates excluded, the final model built by TASSER has a global C_α root-mean-squared deviation from native of 4.6 Å, with a root-mean-squared deviation in the transmembrane helix region of 2.1 Å. Models of several representative GPCRs are compared with mutagenesis and affinity labeling data, and consistent agreement is demonstrated. Structure clustering of the predicted models shows that GPCRs with similar structures tend to belong to a similar functional class even when their sequences are diverse. These results demonstrate the usefulness and robustness of the in silico models for GPCR functional analysis. All predicted GPCR models are freely available for noncommercial users on our Web site (http://www.bioinformatics.buffalo.edu/GPCR).

G protein–coupled receptors (GPCRs) are a large superfamily of integral membrane proteins that transduce signals across the cell membrane. Because of the breadth and importance of the physiological roles undertaken by the GPCR family, many of its members are important pharmacological targets. Although the knowledge of a protein's native structure can provide important insight into understanding its function and for the design of new drugs, the experimental determination of the three-dimensional structure of GPCR membrane proteins has proved to be very difficult. This is demonstrated by the fact that there is only one solved GPCR structure (from bovine rhodopsin) deposited in the Protein Data Bank library. In contrast, there are no human GPCR structures in the Protein Data Bank. To address the need for the tertiary structures of human GPCRs, using just sequence information, the authors use a newly developed threading-assembly-refinement method to generate models for all 907 registered GPCRs in the human genome. About 820 GPCRs are anticipated to have correct topology and transmembrane helix arrangement. A subset of the resulting models is validated by comparison with mutagenesis experimental data, and consistent agreement is demonstrated.

Carboxyl Tail Cysteine Mutants of the Thyrotropin-Releasing Hormone Receptor Type 1 Exhibit Constitutive Signaling: Role of Palmitoylation

We studied the role of carboxyl tail cysteine residues and their palmitoylation in constitutive signaling by the thyrotropin-releasing hormone (TRH) receptor type 1 (TRH-R1) in transfected mammalian cells and in Xenopus laevis oocytes. To study palmitoylation, we inserted a factor Xa cleavage site within the third extracellular loop of TRH-R1, added a carboxyl-terminal C9 immunotag and expressed the mutant receptor in Chinese hamster ovary cells. We identified TRH-R1-specific palmitoylation in the transmembrane helix-7/carboxyl-tail receptor fragment mainly at Cys-335 and Cys-337. In contrast to a mutant truncated at Cys-335 that was reported previously to be constitutively active, a receptor truncated at Lys-338 (K338Stop), which preserves Cys-335 and Cys-337, and C337Stop and N336Stop, which preserve Cys-335, did not exhibit increased constitutive signaling. TRH-R1 mutants substituted singly by Gly or Ser at Cys-335 or Cys-337 did not exhibit constitutive signaling. By contrast, substitution of both cysteines (C335G/C337G or C335S/C337S) yielded TRH-R1 mutants that exhibited marked constitutive signaling in mammalian cells. In the oocyte, constitutive signaling by C335G/C337G resulted in homologous (of C335G/C337G) and heterologous (of M1 muscarinic receptor) desensitization. Because both Cys-335 and Cys-337 have to be substituted or deleted for constitutive signaling, we propose that a single palmitoylation site in the proximal carboxyl tail is sufficient to constrain TRH-R1 in an inactive conformation.