A hydrophobic cluster between transmembrane helices 5 and 6 constrains the thyrotropin-releasing hormone receptor in an inactive conformation

Central Hypothyroidism Due to a TRHR Mutation Causing Impaired Ligand Affinity and Transactivation of Gq

CONTEXT

Central congenital hypothyroidism (CCH) is an underdiagnosed disorder characterized by deficient production and bioactivity of thyroid-stimulating hormone (TSH) leading to low thyroid hormone synthesis. Thyrotropin-releasing hormone (TRH) receptor (TRHR) defects are rare recessive disorders usually associated with incidentally identified CCH and short stature in childhood.

OBJECTIVES

Clinical and genetic characterization of a consanguineous family of Roma origin with central hypothyroidism and identification of underlying molecular mechanisms.

DESIGN

All family members were phenotyped with thyroid hormone profiles, pituitary magnetic resonance imaging, TRH tests, and dynamic tests for other pituitary hormones. Candidate TRH, TRHR, TSHB, and IGSF1 genes were screened for mutations. A mutant TRHR was characterized in vitro and by molecular modeling.

RESULTS

A homozygous missense mutation in TRHR (c.392T > C; p.I131T) was identified in an 8-year-old boy with moderate hypothyroidism (TSH: 2.61 mIU/L, Normal: 0.27 to 4.2; free thyroxine: 9.52 pmol/L, Normal: 10.9 to 25.7) who was overweight (body mass index: 20.4 kg/m2, p91) but had normal stature (122 cm; -0.58 standard deviation). His mother, two brothers, and grandmother were heterozygous for the mutation with isolated hyperthyrotropinemia (TSH: 4.3 to 8 mIU/L). The I131T mutation, in TRHR intracellular loop 2, decreases TRH affinity and increases the half-maximal effective concentration for signaling. Modeling of TRHR-Gq complexes predicts that the mutation disrupts the interaction between receptor and a hydrophobic pocket formed by Gq.

CONCLUSIONS

A unique missense TRHR defect identified in a consanguineous family is associated with central hypothyroidism in homozygotes and hyperthyrotropinemia in heterozygotes, suggesting compensatory elevation of TSH with reduced biopotency. The I131T mutation decreases TRH binding and TRHR-Gq coupling and signaling.

Interactions of vasopressin and oxytocin receptors with vasopressin analogues substituted in position 2 with 3,3'-diphenylalanine--a molecular docking study

Vasopressin and oxytocin receptors belong to the superfamily of G protein-coupled receptors and play an important role in many physiological functions. They are also involved in a number of pathological conditions being important drug targets. In this work, four vasopressin analogues substituted at position 2 with 3,3'-diphenylalanine have been docked into partially flexible vasopressin and oxytocin receptors. The bulky residue at position 2 acts as a structural restraint much stronger in the oxytocin receptor (OTR) than in the vasopressin V2 receptor (V2R), resulting in a different location of the analogues in these receptors. This explains the different, either agonistic or antagonistic, activities of the analogues in V2R and OTR, respectively. In all complexes, the conserved polar residues serve as anchor points for the ligand both in OTR and V2R. Strong interactions of the C-terminus of analogue II ([Mpa(1) ,d-Dpa(2) ,Val(4) ,d-Arg(8) ]VP) with extracellular loop 3 may be responsible for its highest activity at V2R. It also appears that V2R adapts more readily to the docking analogues by conformational changes in the aromatic side chains triggering receptor activation. A weak activity at V1a vasopressin receptor appears to be caused by weak receptor-ligand interactions. Results of this study may facilitate a rational design of new analogues with the highest activity/selectivity at vasopressin and OTRs.

Critical hydrogen bond formation for activation of the angiotensin II type 1 receptor

G protein-coupled receptors contain selectively important residues that play central roles in the conformational changes that occur during receptor activation. Asparagine 111 (N111(3.35)) is such a residue within the angiotensin II type 1 (AT(1)) receptor. Substitution of N111(3.35) for glycine leads to a constitutively active receptor, whereas substitution for tryptophan leads to an inactivable receptor. Here, we analyzed the AT(1) receptor and two mutants (N111G and N111W) by molecular dynamics simulations, which revealed a novel molecular switch involving the strictly conserved residue D74(2.50). Indeed, D74(2.50) forms a stable hydrogen bond (H-bond) with the residue in position 111(3.35) in the wild-type and the inactivable receptor. However, in the constitutively active mutant N111G-AT(1) receptor, residue D74 is reoriented to form a new H-bond with another strictly conserved residue, N46(1.50). When expressed in HEK293 cells, the mutant N46G-AT(1) receptor was poorly activable, although it retained a high binding affinity. Interestingly, the mutant N46G/N111G-AT(1) receptor was also inactivable. Molecular dynamics simulations also revealed the presence of a cluster of hydrophobic residues from transmembrane domains 2, 3, and 7 that appears to stabilize the inactive form of the receptor. Whereas this hydrophobic cluster and the H-bond between D74(2.50) and W111(3.35) are more stable in the inactivable N111W-AT(1) receptor, the mutant N111W/F77A-AT(1) receptor, designed to weaken the hydrophobic core, showed significant agonist-induced signaling. These results support the potential for the formation of an H-bond between residues D74(2.50) and N46(1.50) in the activation of the AT(1) receptor.

Involvement of tryptophan W276 and of two surrounding amino acid residues in the high constitutive activity of the ghrelin receptor GHS-R1a

The human ghrelin receptor (GHS-R1a) is known to display a high level of signaling in the absence of ligand. The Trp276, located in the fully conserved CWXP motif of G protein-coupled receptors, is believed to function as a rotameric switch in these receptors. A comparative modelling of GHS-R1a with the motilin receptor, the most related G protein-coupled receptor to GHS-R1a known to date, but characterized by a very low ligand-independent signaling level, revealed that only two surrounding residues of Trp276, that are Val131 and Ile134, were different from these receptors. We mutated them at once in GHS-R1a to create a "motilin receptor-like" environment of Trp276 in order to study the consequences on GHS-R1a activation. We studied the pharmacological properties of the W276A, V131L-I134M GHS-R1a mutants. Basal as well as maximal ghrelin-induced signaling was assessed both by inositol-phosphate accumulation and SRE pathways. As compared to the wild type receptor, the SRE-luciferase assay displayed a markedly impaired basal activity for W276A whereas that of V131L-I134M was, strikingly, two fold increased. Nevertheless, the efficacy of ghrelin to bind or to stimulate mutant receptors remained unchanged. It is concluded that Trp276, Val131 and Ile134 have a significant impact on constitutive signaling of GHS-R1a, V131L-I134M being the first example of a GHS-R1a mutant with a higher basal activity than the wild type receptor.

A conserved aromatic lock for the tryptophan rotameric switch in TM-VI of seven-transmembrane receptors

The conserved tryptophan in position 13 of TM-VI (Trp-VI:13 or Trp-6.48) of the CWXP motif located at the bottom of the main ligand-binding pocket in TM-VI is believed to function as a rotameric microswitch in the activation process of seven-transmembrane (7TM) receptors. Molecular dynamics simulations in rhodopsin demonstrated that rotation around the chi1 torsion angle of Trp-VI:13 brings its side chain close to the equally highly conserved Phe-V:13 (Phe-5.47) in TM-V. In the ghrelin receptor, engineering of high affinity metal-ion sites between these positions confirmed their close spatial proximity. Mutational analysis was performed in the ghrelin receptor with multiple substitutions and with Ala substitutions in GPR119, GPR39, and the beta(2)-adrenergic receptor as well as the NK1 receptor. In all of these cases, it was found that mutation of the Trp-VI:13 rotameric switch itself eliminated the constitutive signaling and strongly impaired agonist-induced signaling without affecting agonist affinity and potency. Ala substitution of Phe-V:13, the presumed interaction partner for Trp-VI:13, also in all cases impaired both the constitutive and the agonist-induced receptor signaling, but not to the same degree as observed in the constructs where Trp-VI:13 itself was mutated, but again without affecting agonist potency. In a proposed active receptor conformation generated by molecular simulations, where the extracellular segment of TM-VI is tilted inwards in the main ligand-binding pocket, Trp-VI:13 could rotate into a position where it obtained an ideal aromatic-aromatic interaction with Phe-V:13. It is concluded that Phe-V:13 can serve as an aromatic lock for the proposed active conformation of the Trp-VI:13 rotameric switch, being involved in the global movement of TM-V and TM-VI in 7TM receptor activation.

Highlights on the development of A(2A) adenosine receptor agonists and antagonists

Although significant progress has been made in the past few decades demonstrating that adenosine modulates a variety of physiological and pathophysiological processes through the interaction with four subtypes of a family of cell-surface G-protein-coupled receptors, clinical evaluation of some adenosine receptor ligands has been discontinued. Major problems include side effects due to the wide distribution of adenosine receptors, low brain penetration (which is important for the targeting of CNS diseases), short half-life of compounds, or a lack of effects, in some cases perhaps due to receptor desensitization or to low receptor density in the targeted tissue. Currently, three A(2A) adenosine receptor agonists have begun phase III studies. Two of them are therapeutically evaluated as pharmacologic stress agents and the third proved to be effective in the treatment of acute spinal cord injury (SCI), while avoiding the adverse effects of steroid agents. On the other hand, the great interest in the field of A(2A) adenosine receptor antagonists is related to their application in neurodegenerative disorders, in particular, Parkinson's disease, and some of them are currently in various stages of evaluation. This review presents an update of medicinal chemistry and molecular recognition of A(2A) adenosine receptor agonists and antagonists, and stresses the strong need for more selective ligands at the A(2A) human subtype.

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.

Analysis of the activation mechanism of the guinea-pig Histamine H1-receptor

The Histamine H(1)-receptor (H1R), belonging to the amine receptor-class of family A of the G-protein coupled receptors (GPCRs) gets activated by agonists. The consequence is a conformational change of the receptor, which may involve the binding-pocket. So, for a good prediction of the binding-mode of an agonist, it is necessary to have knowledge about these conformational changes. Meanwhile some experimental data about the structural changes of GPCRs during activation exist. Based on homology modeling of the guinea-pig H1R (gpH1R), using the crystal structure of bovine rhodopsin as template, we performed several MD simulations with distance restraints in order to get an inactive and an active structure of the gpH1R. The calculations led to a Phe6.44/Trp6.48/Phe6.52-switch and linearization of the proline kinked transmembrane helix VI during receptor activation. Our calculations showed that the Trp6.48/Phe6.52-switch induces a conformational change in Phe6.44, which slides between transmembrane helices III and VI. Additionally we observed a hydrogen bond interaction of Ser3.39 with Asn7.45 in the inactive gpH1R, but because of a counterclockwise rotation of transmembrane helix III Ser3.39 establishes a water-mediated hydrogen bond to Asp2.50 in the active gpH1R. Additionally we simulated a possible mechanism for receptor activation with a modified LigPath-algorithm.

Bidirectional, iterative approach to the structural delineation of the functional "chemoprint" in GPR40 for agonist recognition

GPR40, free fatty acid receptor 1 (FFAR1), is a member of the GPCR superfamily and a possible target for the treatment of type 2 diabetes. In this work, we conducted a bidirectional iterative investigation, including computational modeling and site-directed mutagenesis, aimed at delineating amino acid residues forming the functional "chemoprint" of GPR40 for agonist recognition. The computational and experimental studies revolved around the recognition of the potent synthetic agonist GW9508. Our experimentally supported model suggested that H137(4.56), R183(5.39), N244(6.55), and R258(7.35) are directly involved in interactions with the ligand. We have proposed a polarized NH-pi interaction between H137(4.56) and GW9508 as one of the contributing forces leading to the high potency of GW9508. The modeling approach presented in this work provides a general strategy for the exploration of receptor-ligand interactions in G-protein coupled receptors beginning prior to acquisition of experimental data.

Suppression of thyrotropin receptor constitutive activity by a monoclonal antibody with inverse agonist activity

TSH binding to the TSH receptor (TSHR) induces thyrocyte growth and proliferation primarily by activating the adenylyl cyclase signaling pathway. Relative to the other glycoprotein hormone receptors, the TSHR has considerable ligand-independent (constitutive) activity. We describe a TSHR monoclonal antibody (CS-17) with the previously unrecognized property of being an inverse agonist for TSHR constitutive activity. This property is retained, even when constitutive activity is extremely high consequent to diverse TSHR extracellular region mutations. A similar effect on an activating mutation at the base of the sixth transmembrane helix (not accessible to direct CS-17 contact) indicates that CS-17 is acting allosterically. Administered to mice in vivo, CS-17 reduces serum T(4) levels. The CS-17 epitope is conformational and a significant portion lies in the C-terminal region of the TSHR leucine-rich domain (residues 260-289). By interacting with the large TSHR extracellular domain, CS-17 is, to our knowledge, the first antibody reported to be an inverse agonist for a member of the G protein receptor superfamily. After humanization of its murine constant region, CS-17 has the potential to be an adjunctive therapeutic agent in athyreotic patients with residual well-differentiated thyroid carcinoma as well as pending definitive treatment in some selected hyperthyroidism states.

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.

Analysis of interactions responsible for vasopressin binding to human neurohypophyseal hormone receptors-molecular dynamics study of the activated receptor-vasopressin-G(alpha) systems

Vasopressin (CYFQNCPRG-NH(2), AVP) is a semicyclic endogenous peptide, which exerts a variety of biological effects in mammals. The main physiological roles of AVP are the regulation of water balance and the control of blood pressure and adrenocorticotropin hormone (ACTH) secretion, mediated via three different subtypes of vasopressin receptors: V1a, V1b and V2 receptors (V1aR, V1bR and V2R, respectively). They are the members of the class A, G-protein-coupled receptors (GPCRs). AVP also modulates several behavioral and social functions. In this study, the interactions responsible for AVP binding to vasopressin V1a and V2 receptors versus the closely related oxytocin ([I3,L8]AVP, OT) receptor (OTR) have been investigated. Three-dimensional models of the activated receptors were constructed using multiple sequence alignment, followed by homology modeling using the complex of activated rhodopsin with Gt(alpha) C-terminal peptide of transducin MII-Gt(338-350) prototype as a template. AVP was docked into the receptor-G(alpha) systems. The three lowest-energy pairs of receptor-AVP-G(alpha) (two complexes per each receptor) were selected. The 1-ns unconstrained molecular dynamics (MD) of complexes embedded into the fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) lipid bilayer was conducted in the AMBER 7.0 force field. Six relaxed receptor-AVP-G(alpha) models were obtained. The residues responsible for AVP binding to vasopressin receptors have been identified and a different mechanism of AVP binding to V2R than to V1aR has been proposed.

A chemogenomic analysis of the transmembrane binding cavity of human G-protein-coupled receptors

The amino acid sequences of 369 human nonolfactory G-protein-coupled receptors (GPCRs) have been aligned at the seven transmembrane domain (TM) and used to extract the nature of 30 critical residues supposed--from the X-ray structure of bovine rhodopsin bound to retinal--to line the TM binding cavity of ground-state receptors. Interestingly, the clustering of human GPCRs from these 30 residues mirrors the recently described phylogenetic tree of full-sequence human GPCRs (Fredriksson et al., Mol Pharmacol 2003;63:1256-1272) with few exceptions. A TM cavity could be found for all investigated GPCRs with physicochemical properties matching that of their cognate ligands. The current approach allows a very fast comparison of most human GPCRs from the focused perspective of the predicted TM cavity and permits to easily detect key residues that drive ligand selectivity or promiscuity.

Investigation of mechanism of desmopressin binding in vasopressin V2 receptor versus vasopressin V1a and oxytocin receptors: Molecular dynamics simulation of the agonist-bound state in the membrane–aqueous system

The vasopressin V2 receptor (V2R) belongs to the Class A G protein-coupled receptors (GPCRs). V2R is expressed in the renal collecting duct (CD), where it mediates the antidiuretic action of the neurohypophyseal hormone arginine vasopressin (CYFQNCPRG-NH2, AVP). Desmopressin ([1-deamino, 8-D]AVP, dDAVP) is strong selective V2R agonist with negligible pressor and uterotonic activity. In this paper, the interactions responsible for binding of dDAVP to vasopressin V2 receptor versus vasopressin V1a and oxytocin receptors has been examined. Three-dimensional activated models of the receptors were constructed using the multiple sequence alignment and the complex of activated rhodopsin with Gt(alpha) C-terminal peptide of transducin MII-Gt(alpha) (338-350) prototype (Slusarz, R.; Ciarkowski, J. Acta Biochim Pol 2004 51, 129-136) as a template. The 1-ns unconstrained molecular dynamics (MD) of receptor-dDAVP complexes immersed in the fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) membrane model was conducted in an Amber 7.0 force field. Highly conserved transmembrane residues have been proposed as being responsible for V2R activation and G protein coupling. Molecular mechanism of the dDAVP binding has been suggested. The internal water molecules involved in an intricate network of the hydrogen bonds inside the receptor cavity have been identified and their role in the stabilization of the agonist-bound state proposed.

Partial Agonism, Neutral Antagonism, and Inverse Agonism at the Human Wild-Type and Constitutively Active Cholecystokinin-2 Receptors

Cholecystokinin receptor-2 (CCK2R) is a G protein receptor that regulates a number of physiological functions. Activation of CCK2R and/or expression of a constitutively active CCK2R variant may contribute to human diseases, including digestive cancers. Search for antagonists of the CCK2R has been an important challenge during the last few years, leading to discovery of a set of chemically distinct compounds. However, several early-discovered antagonists turned out to be partial agonists. In this context, we carried out pharmacological characterization of six CCK2R antagonists using COS-7 cells expressing the human CCK2R or a CCK2R mutant having a robust constitutive activity on inositol phosphates production, and we investigated the molecular mechanisms which, at a CCK2R binding site, account for these features. Results indicated that three compounds, 3R(+)-N-(2,3-dihydro-1-methyl-2-oxo-5-phenyl-1H-1,4-benzodiazepin-3-yl)-N'-(3-methylphenyl)urea (L365,260), 4-{[2-[[3-(lH-indol-3-yl)-2-methyl-1-oxo-2-[[[1.7.7-trimethyl-bicyclo[2.2.1]hept-2-yl)-oxy]carbonyl]amino]propyl]amino]-1-phenylethyl]amino-4-oxo-[lS-la.2[S*(S*)]4a]}-butanoate N-methyl-D-glucamine (PD135,158), and (R)-1-naphthalenepropanoic acid, b-[2-[[2-(8-azaspiro-[4.5]dec-8-ylcarbonyl)-4,6-dimethylphenyl]amino]-2-oxoethyl] (CR2945), were partial agonists; one molecule, 1-[(R)-2,3-dihydro-1-(2,3-dihydro-1-(2-methylphenacyl)-2-oxo-5-phenyl-1H-1,4-benzodiazepin-3-yl]-3-(3-methylphenyl)urea (YM022), was a neutral antagonist; and two compounds, N-(+)-[1-(adamant-1-ylmethyl)-2,4-dioxo-5-phenyl2,3,4,5-tetrahydro-1H-1,5-benzodiazepin-3-yl]-N'-phenylurea (GV150,013X) and ([(N-[methoxy-3 phenyl] N-[N-methyl N-phenyl carbamoylmethyl], carbomoyl-methyl)-3 ureido]-3-phenyl)2-propionic acid (RPR101,048), were inverse agonists. Furthermore, target- and pharmacophore-based docking of ligands followed by molecular dynamic simulation experiments resulted in consistent motion of aromatic residues belonging to a network presumably important for activation, thus providing the first structural explanations for the different pharmacological profiles of tested compounds. This study confirms that several referenced so-called antagonists are in fact partial agonists, and because of this undesired activity, we suggest that newly generated molecules should be preferred to efficiently block CCK2R-related physiological effects. Furthermore, data on the structural basis for the different pharmacological features of CCK2R ligands will serve to further clarify CCK2R mechanism of activation.

A Model of Inverse Agonist Action at Thyrotropin-Releasing Hormone Receptor Type 1: Role of a Conserved Tryptophan in Helix 6

A binding pocket for thyrotropin-releasing hormone (TRH) within the transmembrane helices of the TRH receptor type 1 (TRH-R1) has been identified based on experimental evidence and computer simulations. To determine the binding site for a competitive inverse agonist, midazolam, three of the four residues that directly contact TRH and other residues that restrain TRH-R1 in an inactive conformation were screened by mutagenesis and binding assays. We found that two residues that directly contact TRH, Asn-110 in transmembrane helix 3 (3.37) and Arg-306 in transmembrane helix 7 (7.39), were important for midazolam binding but another, Tyr-282 in transmembrane helix 6 (6.51), was not. A highly conserved residue, Trp-279 in transmembrane helix 6 (6.48), which was reported to be critical in stabilizing TRH-R1 in an inactive state but not for TRH binding, was critical for midazolam binding. We used our previous model of the unoccupied TRH-R1 to generate a model of the TRH-R1/midazolam complex. The experimental results and the molecular model of the complex suggest that midazolam binds to TRH-R1 within a transmembrane helical pocket that partially overlaps the TRH binding pocket. This result is consistent with the competitive antagonism of midazolam binding. We suggest that the mechanism of inverse agonism effected by midazolam involves its direct interaction with Trp-279, which contributes to the stabilization of the inactive conformation of TRH-R1.

Homology model of the CB1 cannabinoid receptor: sites critical for nonclassical cannabinoid agonist interaction

Association of cannabimimetic compounds such as cannabinoids, aminoalkylindoles (AAIs), and arachidonylethanolamide (anandamide) with the brain cannabinoid (CB(1)) receptor activates G-proteins and relays signals to regulate neuronal functions. A CB(1) receptor homology model was constructed using the published x-ray crystal structure of bovine rhodopsin (Palczewski et al., Science, 2000, Vol. 289, pp. 739-745) in the conformation most likely to represent the "high-affinity" state for agonist binding to G-protein coupled receptors (GPCRs). A molecular docking approach that combined Monte Carlo and molecular dynamics simulations was used to identify the putative binding conformations of nonclassical cannabinoid agonists, including AC-bicyclic CP47497 and CP55940, and ACD-tricyclic CP55244. Placement of these ligands was based upon the assumption of a critical hydrogen bond between the A-ring OH and the side chain N of Lys192 in transmembrane helix 3. We evaluated two alternative binding conformations, C3-in and C3-out, denoting the directionality of the ligand C3 side chain within the receptor with respect to the inside or the outside of the cell. Assuming both the C3-in or C3-out conformation, the calculated ligand-receptor binding energy (DeltaE(bind)) was correlated with the experimentally observed binding affinity (K(i)) for a series of nonclassical cannabinoid agonists. The C3-in conformation was marginally better than the alternative C3-out conformation in predicting the rank order of the tested nonclassical cannabinoid analogs. Adopting the C3-in conformation due to the greater number of receptor interactions with known pharmacophoric elements of the ligand, key residues were identified comprising the presumed hydrophobic pocket that interacts with the C3 side chain of cannabinoid agonists. Key hydrogen bonds would form between both K3.28(192) and E(258) and the A-ring OH, and between Q(261) and the C-ring C-12 hydroxypropyl. In summary, the present study represents one of the first attempts to construct a homology model of the CB(1) cannabinoid receptor based upon the published bovine rhodopsin x-ray crystal structure and to elucidate the putative ligand binding site for nonclassical cannabinoid agonists. We postulated sites of the CB(1) receptor critical for the ligand interaction, including the hydrophobic pocket interacting with the key pharmacophoric moiety, the C3 side chain. More work is needed to delineate between two alternative (and possibly other) binding conformations of the nonclassical cannabinoid ligands within the CB(1) receptor. The present study provides a consistent framework for further investigation of the CB(1) receptor-ligand interaction and for the study of CB(1) receptor activation.

Random mutagenesis of the human adenosine A2B receptor followed by growth selection in yeast. Identification of constitutively active and gain of function mutations

To gain insight in spontaneous as well as agonist-induced activation of the human adenosine A2B receptor, we applied a random mutagenesis approach in yeast to create a large number of receptor mutants and selected mutants of interest with a robust screening assay based on growth. The amino acid sequence of 14 mutated receptors was determined. All these mutated receptors displayed constitutive activity. In particular, single-point mutations at T42A, V54L, and F84S and a triple-point mutation at N36S, T42A, and T66A resulted in high constitutive activity. In addition, a C-terminally truncated (after Lys269) mutant, Q214L I230N V240M V250M N254Y T257S K269stop, was highly constitutively active. The T42A, V54L, and F84S mutants showed a considerable decrease, 4.9- to 6.9-fold, in the EC50 value of 5'-N-ethylcarboxamidoadenosine (NECA), an adenosine analog. Combined mutation of I242T, K269R, V284A, and H302Q, as well as F84L together with S95G, resulted in an even greater potency of NECA of 10- and 18-fold, respectively. In fact, all constitutively active mutants had an increased potency for NECA. This suggests that the wild-type (wt) human A2B receptor itself is rather silent, which may explain the low affinity of agonists for this receptor. To verify the ability of the mutant receptors to activate mammalian second messenger systems, cAMP experiments were performed in CHO cells stably expressing the wt and T42A receptors. These experiments confirmed the increased sensitivity of T42A for NECA, because the EC50 values of T42A and the wt receptor were 0.15 +/- 0.04 and 1.3 +/- 0.4 microM, respectively.

Activation of CCR5 by chemokines involves an aromatic cluster between transmembrane helices 2 and 3

CCR5 is a G protein-coupled receptor responding to four natural agonists, the chemokines RANTES (regulated on activation normal T cell expressed and secreted), macrophage inflammatory protein (MIP)-1 alpha, MIP-1 beta, and monocyte chemotactic protein (MCP)-2, and is the main co-receptor for the macrophage-tropic human immunodeficiency virus strains. We have previously identified a structural motif in the second transmembrane helix of CCR5, which plays a crucial role in the mechanism of receptor activation. We now report the specific role of aromatic residues in helices 2 and 3 of CCR5 in this mechanism. Using site-directed mutagenesis and molecular modeling in a combined approach, we demonstrate that a cluster of aromatic residues at the extracellular border of these two helices are involved in chemokine-induced activation. These aromatic residues are involved in interhelical interactions that are key for the conformation of the helices and govern the functional response to chemokines in a ligand-specific manner. We therefore suggest that transmembrane helices 2 and 3 contain important structural elements for the activation mechanism of chemokine receptors, and possibly other related receptors as well.

Correlation between basal signaling and internalization of thyrotropin-releasing hormone receptors: evidence for involvement of similar receptor conformations

Previous studies have shown that rat thyrotropin-releasing hormone (TRH) receptor type 2 exhibits higher basal signaling activity and internalizes more rapidly upon agonist binding than rat TRH receptor type 1. The mouse TRH receptor type 2 (mR2) was recently cloned and, similar to its rat homolog, shows a higher basal signaling activity than mR1. Taking advantage of the high degree of sequence homology between mR1 and mR2, we used chimeras/mutants of these receptors to gain insight into the properties of the receptors that influence internalization and basal signaling. Chimeric receptors that have the mR1 extracellular and transmembrane domains with the carboxyl terminus and intracellular loops of mR2 (R1/R2-tail; R1/R2-I3,tail; R1/R2-I2,3,tail; R1/R2-I1,2,3,tail) exhibited internalization rates and basal activities that were similar to that of mR1. In contrast, a chimeric receptor with the extracellular and transmembrane domains of mR2 and the carboxyl terminus of mR1 exhibited the more rapid internalization rate and higher basal signaling activity characteristic of mR2. We showed previously that mutation of a highly conserved tryptophan to alanine caused mR1 to exhibit a high basal signaling activity and rapid internalization rate. In contrast, mutation of this tryptophan to alanine in mR2 decreased the rate of internalization and inhibited basal signaling activity. The rates of receptor internalization did not correlate with the binding affinities, coupling efficiencies, or potencies of the receptors. Thus, we observed that receptors with more rapid internalization rates showed relatively higher basal signaling activities, whereas receptors with lower basal signaling activities showed slower internalization rates. These data suggest that similar receptor conformations are required for productive coupling to signaling G proteins and to proteins involved in internalization.

A critical role for a tyrosine residue in the cannabinoid receptors for ligand recognition

Previous mutation and modeling studies have identified an aromatic cluster in the transmembrane helix (TMH) 3-4-5 region as important for ligand binding at the CB(1) and CB(2) cannabinoid receptors. Through novel mixed mode Monte Carlo/Stochastic Dynamics (MC/SD) calculations, we tested the importance of aromaticity at position 5.39(275) in CB(1). MC/SD calculations were performed on wild-type (WT) CB(1) and two mutants, Y5.39(275)F and Y5.39(275)I. Results indicated that while the CB(1) Y5.39(275)F mutant is very similar to WT, the Y5.39(275)I mutant shows pronounced topology changes in the TMH 3-4-5 region. Site-directed mutagenesis studies of tyrosine 5.39 to phenylalanine (Y-->F) or isoleucine (Y-->I) in both CB(1) and CB(2) were performed to determine the functional role of this amino acid in each receptor subtype. HEK 293 cells transfected with mutant receptor cDNAs were evaluated in radioligand binding and cyclic AMP assays. The CB(1) mutant and WT receptors were also co-expressed with G-protein-coupled inwardly rectifying channels (GIRK1 and GIRK4) in Xenopus oocytes to assess functional coupling. The Y-->F mutation resulted in cannnabinoid receptors with subtle differences in WT binding and signal transduction. In contrast, the Y-->I mutations produced receptors that could not produce signal transduction or bind to multiple cannabinoid compounds. However, immunofluorescence data indicate that the Y-->I mutation was compartmentalized and expressed at a level similar to that of the WT cannabinoid receptor. These results underscore the importance of aromaticity at position CB(1) 5.39(275) and CB(2) 5.39(191) for ligand recognition in the cannabinoid receptors.

Control of conformational equilibria in the human B2 bradykinin receptor. Modeling of nonpeptidic ligand action and comparison to the rhodopsin structure

A prototypic study of the molecular mechanisms of activation or inactivation of peptide hormone G protein-coupled receptors was carried out on the human B2 bradykinin receptor. A detailed pharmacological analysis of receptor mutants possessing either increased constitutive activity or impaired activation or ligand recognition allowed us to propose key residues participating in intramolecular interaction networks stabilizing receptor inactive or active conformations: Asn(113) and Tyr(115) (TM III), Trp(256) and Phe(259) (TM VI), Tyr(295) (TM VII) which are homologous of the rhodopsin residues Gly(120), Glu(122), Trp(265), Tyr(268), and Lys(296), respectively. An essential experimental finding was the spatial proximity between Asn(113), which is the cornerstone of inactive conformations, and Trp(256) which plays a subtle role in controlling the balance between active and inactive conformations. Molecular modeling and mutagenesis data showed that Trp(256) and Tyr(295) constitute, together with Gln(288), receptor contact points with original nonpeptidic ligands. It provided an explanation for the ligand inverse agonist behavior on the WT receptor, with underlying restricted motions of TMs III, VI, and VII, and its agonist behavior on the Ala(113) and Phe(256) constitutively activated mutants. These data on the B2 receptor emphasize that conformational equilibria are controlled in a coordinated fashion by key residues which are located at strategic positions for several G protein-coupled receptors. They are discussed in comparison with the recently determined rhodopsin crystallographic structure.

Characterization of the molecular motions of constitutively active G protein-coupled receptors for parathyroid hormone

The molecular mechanism of constitutive activity of the G protein-coupled receptor for human parathyroid hormone (PTH1) has been examined by molecular dynamics (MD) simulations. The single point mutations H223R, T410P, and I458R, of the PTH1 receptor result in ligand-independent receptor activation. Extensive MD simulations indicate that each of the mutations, through different mechanisms, lead to very similar conformational changes of the third intracellular loop. The structural changes, centered on K405 in the C-terminus of the third intracellular loop, can be traced back to the single-point mutations by calculation of the forces and torques responsible for the collective motions of the receptor. This analysis indicates a direct correlation between the conformational preferences of the cytoplasmic loop and the mutations in different locations of the receptor: TM2 (H223R), TM6 (T410P), and TM7 (1458R). Given the pivotal role of the third intracellular loop of PTH1 in coupling to the G proteins, the structural changes induced by these single-point mutations may be responsible for the ligand-free activation of the receptor. These results coupled with the high-resolution structure of the third cytoplasmic loop of PTH1, previously determined in our laboratory, provide unique insight into the mechanism of ligand free activation of the PTH1 receptor.

Charged residues at the intracellular boundary of transmembrane helices 2 and 3 independently affect constitutive activity of Kaposi's sarcoma-associated herpesvirus G protein-coupled receptor

Because charged residues at the intracellular ends of transmembrane helix (TMH) 2 and TMH3 of G protein-coupled receptors (GPCRs) affect signaling, we performed mutational analysis of these residues in the constitutively signaling Kaposi's sarcoma-associated herpesvirus GPCR (KSHV-GPCR). KSHV-GPCR contains the amino acid sequence Val-Arg-Tyr rather than the Asp/Glu-Arg-Tyr ((D/E)RY) motif at the intracellular end of TMH3. Mutation of Arg-143 to Ala (R143A) or Gln (R143Q) abolished constitutive signaling whereas R143K exhibited 50% of the basal activity of KSHV-GPCR. R143A was not stimulated by agonist, whereas R143Q was stimulated by growth-related oncogene-alpha, and R143K, similar to KSHV-GPCR, was stimulated further. These findings show that Arg-143 is critical for signal generation in KSHV-GPCR. In other GPCRs, Arg in this position may act as a signaling switch by movement of its sidechain from a hydrophilic pocket in the TMH bundle to a position outside the bundle. In rhodopsin, the Arg of Glu-Arg-Tyr interacts with the adjacent Asp to constrain Arg outside the TMH bundle. V142D was 70% more active than KSHV-GPCR, suggesting that an Arg residue, which is constrained outside the bundle by interacting with Asp-142, leads to a receptor that signals more actively. Because the usually conserved Asp in the middle of TMH2 is not present in KSHV-GPCR, we tested whether Asp-83 at the intracellular end of TMH2 was involved in signaling. D83N and D83A were 110 and 190% more active than KSHV-GPCR, respectively. The double mutant D83A/V142D was 510% more active than KSHV-GPCR. That is, cosubstitutions of Asp-83 by Ala and Val-142 by Asp act synergistically to increase basal signaling. A model of KSHV-GPCR predicts that Arg-143 interacts with residues in the TMH bundle and that the sidechain of Asp-83 does not interact with Arg-143. These data are consistent with the hypothesis that Arg-143 and Asp-83 independently affect the signaling activity of KSHV-GPCR.

Minireview: Insights into G protein-coupled receptor function using molecular models

G protein-coupled receptors (GPCRs) represent the largest family of signal-transducing molecules known. They convey signals for light and many extracellular regulatory molecules. GPCRs have been found to be dysfunctional/dysregulated in a growing number of human diseases and have been estimated to be the targets of more than 30% of the drugs used in clinical medicine today. Thus, understanding how GPCRs function at the molecular level is an important goal of biological research. In order to understand function at this level, it is necessary to delineate the 3D structure of these receptors. Recently, the 3D structure of rhodopsin has been resolved, but in the absence of experimentally determined 3D structures of other GPCRs, a powerful approach is to construct a theoretical model for the receptor and refine it based on experimental results. Computer-generated models for many GPCRs have been constructed. In this article, we will review these studies. We will place the greatest emphasis on an iterative, bi-directional approach in which models are used to generate hypotheses that are tested by experimentation and the experimental findings are, in turn, used to refine the model. The success of this approach is due to the synergistic interaction between theory and experiment.

Interaction between transmembrane domains five and six of the alpha -factor receptor

The alpha-factor pheromone receptor (STE2) activates a G protein signal pathway that induces conjugation of the yeast Saccharomyces cerevisiae. Previous studies implicated the third intracellular loop of this receptor in G protein activation. Therefore, the roles of transmembrane domains five and six (TMD5 and -6) that bracket the third intracellular loop were analyzed by scanning mutagenesis in which each residue was substituted with cysteine. Out of 42 mutants examined, four constitutive mutants and two strong loss-of-function mutants were identified. Double mutants combining Cys substitutions in TMD5 and TMD6 gave a broader range of phenotypes. Interestingly, a V223C mutation in TMD5 caused constitutive activity when combined with the L247C, L248C, or S251C mutations in TMD6. Also, the L226C mutation in TMD5 caused constitutive activity when combined with either the M250C or S251C mutations in TMD6. The residues affected by these mutations are predicted to fall on one side of their respective helices, suggesting that they may interact. In support of this, cysteines substituted at position 223 in TMD5 and position 247 in TMD6 formed a disulfide bond, providing the first direct evidence of an interaction between these transmembrane domains in the alpha-factor receptor. Altogether, these results identify an important region of interaction between conserved hydrophobic regions at the base of TMD5 and TMD6 that is required for the proper regulation of receptor signaling.

Conformational behavior of ionic self-complementary peptides

Several de novo designed ionic peptides that are able to undergo conformational change under the influence of temperature and pH were studied. These peptides have two distinct surfaces with regular repeats of alternating hydrophilic and hydrophobic side chains. This permits extensive ionic and hydrophobic interactions resulting in the formation of stable beta-sheet assemblies. The other defining characteristic of this type of peptide is a cluster of negatively charged aspartic or glutamic acid residues located toward the N-terminus and positively charged arginine or lysine residues located toward the C-terminus. This arrangement of charge balances the alpha-helical dipole moment (C --> N), resulting in a strong tendency to form stable alpha-helices as well. Therefore, these peptides can form both stable alpha-helices and beta-sheets. They are also able to undergo abrupt structural transformations between these structures induced by temperature and pH changes. The amino acid sequence of these peptides permits both stable beta-sheet and alpha-helix formation, resulting in a balance between these two forms as governed by the environment. Some segments in proteins may also undergo conformational changes in response to environmental changes. Analyzing the plasticity and dynamics of this type of peptide may provide insight into amyloid formation. Since these peptides have dynamic secondary structure, they will serve to refine our general understanding of protein structure.

The conformational switch in 7-transmembrane receptors: the muscarinic receptor paradigm

The rhodopsin-like superfamily of 7-transmembrane receptors is the largest class of signalling molecules in the mammalian genome. Recently, a combination of mutagenesis, biophysical and modelling studies have suggested a credible model for the alpha-carbon backbone in the transmembrane region of the 7-transmembrane receptors, and have started to reveal the structural basis of the conformational switch from the inactive to the active state. A key feature may be the replacement of a network of radial constraints, centred on transmembrane helix three, which stabilise the inactive ground state of the receptor by a new set of axial interactions which help to stabilise the activated state. Transmembrane helix three may act as a rotary switch in the activation mechanism.