Membrane Environment Imposes Unique Selection Pressures on Transmembrane Domains of G Protein-Coupled Receptors

We have investigated the influence of the plasma membrane environment on the molecular evolution of G protein-coupled receptors (GPCRs), the largest receptor family in Metazoa. In particular, we have analyzed the site-specific rate variation across the two primary structural partitions, transmembrane (TM) and extramembrane (EM), of these membrane proteins. We find that TM domains evolve more slowly than do EM domains, though TM domains display increased rate heterogeneity relative to their EM counterparts. Although the majority of residues across GPCRs experience strong to weak purifying selection, many GPCRs experience positive selection at both TM and EM residues, albeit with a slight bias towards the EM. Further, a subset of GPCRs, chemosensory receptors (including olfactory and taste receptors), exhibit increased rates of evolution relative to other GPCRs, an effect which is more pronounced in their TM spans. Although it has been previously suggested that the TM’s low evolutionary rate is caused by their high percentage of buried residues, we show that their attenuated rate seems to stem from the strong biophysical constraints of the membrane itself, or by functional requirements. In spite of the strong evolutionary constraints acting on the TM spans of GPCRs, positive selection and high levels of evolutionary rate variability are common. Thus, biophysical constraints should not be presumed to preclude a protein’s ability to evolve.     

Formation of new genes explains lower intron density in mammalian Rhodopsin G protein-coupled receptors

Mammalian G protein-coupled receptor (GPCR) genes are characterised by a large proportion of intronless genes or a lower density of introns when compared with GPCRs of invertebrates. It is unclear which mechanisms have influenced intron density in this protein family, which is one of the largest in the mammalian genomes. We used a combination of Hidden Markov Models (HMM) and BLAST searches to establish the comprehensive repertoire of Rhodopsin GPCRs from seven species and performed overall alignments and phylogenetic analysis using the maximum parsimony method for over 1400 receptors in 12 subgroups. We identified 14 different Ancestral Receptor Groups (ARGs) that have members in both vertebrate and invertebrate species. We found that there exists a remarkable difference in the intron density among ancestral and new Rhodopsin GPCRs. The intron density among ARGs members was more than 3.5-fold higher than that within non-ARG members and more than 2-fold higher when considering only the 7TM region. This suggests that the new GPCR genes have been predominantly formed intronless while the ancestral receptors likely accumulated introns during their evolution. Many of the intron positions found in mammalian ARG receptor sequences were found to be present in orthologue invertebrate receptors suggesting that these intron positions are ancient. This analysis also revealed that one intron position is much more frequent than any other position and it is common for a number of phylogenetically different Rhodopsin GPCR groups. This intron position lies within a functionally important, conserved, DRY motif which may form a proto-splice site that could contribute to positional intron insertion. Moreover, we have found that other receptor motifs, similar to DRY, also contain introns between the second and third nucleotide of the arginine codon which also forms a proto-splice site. Our analysis presents compelling evidence that there was not a major loss of introns in mammalian GPCRs and formation of new GPCRs among mammals explains why these have fewer introns compared to invertebrate GPCRs. We also discuss and speculate about the possible role of different RNA- and DNA-based mechanisms of intron insertion and loss.     

Comparison of class A and D G protein-coupled receptors: common features in structure and activation

All G protein-coupled receptors (GPCRs) share a common seven TM helix architecture and the ability to activate heterotrimeric G proteins. Nevertheless, these receptors have widely divergent sequences with no significant homology. We present a detailed structure-function comparison of the very divergent Class A and D receptors to address whether there is a common activation mechanism across the GPCR superfamily. The Class A and D receptors are represented by the vertebrate visual pigment rhodopsin and the yeast alpha-factor pheromone receptor Ste2, respectively. Conserved amino acids within each specific receptor class and amino acids where mutation alters receptor function were located in the structures of rhodopsin and Ste2 to assess whether there are functionally equivalent positions or regions within these receptors. We find several general similarities that are quite striking. First, strongly polar amino acids mediate helix interactions. Their mutation generally leads to loss of function or constitutive activity. Second, small and weakly polar amino acids facilitate tight helix packing. Third, proline is essential at similar positions in transmembrane helices 6 and 7 of both receptors. Mapping the specific location of the conserved amino acids and sites of constitutively active mutations identified conserved microdomains on transmembrane helices H3, H6, and H7, suggesting that there are underlying similarities in the mechanism of the widely divergent Class A and Class D receptors.     

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.     

The role of receptor oligomerization in modulating the expression and function of leukocyte adhesion-G protein-coupled receptors

The human leukocyte adhesion-G protein-coupled receptors (GPCRs), the epidermal growth factor (EGF)-TM7 proteins, are shown here to function as homo- and hetero-oligomers. Using cell surface cross-linking, co-immunoprecipitation, and fluorescence resonance energy transfer analysis of EMR2, an EGF-TM7 receptor predominantly expressed in myeloid cells, we demonstrate that it forms dimers in a reaction mediated exclusively by the TM7 moiety. We have also identified a naturally occurring but structurally unstable EMR2 splice variant that acts as a dominant negative modulator by dimerizing with the wild type receptor and down-regulating its expression. Additionally, heterodimerization between closely related EGF-TM7 members is shown to result in the modulation of expression and ligand binding properties of the receptors. These findings suggest that receptor homo- and hetero-oligomerization play a regulatory role in modulating the expression and function of leukocyte adhesion-GPCRs.     

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

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

The G protein-coupled receptors in the pufferfish <Emphasis Type="Italic">Takifugu rubripes</Emphasis>

Background

Guanine protein-coupled receptors (GPCRs) constitute a eukaryotic transmembrane protein family and function as “molecular switches” in the second messenger cascades and are found in all organisms between yeast and humans. They form the single, biggest drug-target family due to their versatility of action and their role in several physiological functions, being active players in detecting the presence of light, a variety of smells and tastes, amino acids, nucleotides, lipids, chemicals etc. in the environment of the cell. Comparative genomic studies on model organisms provide information on target receptors in humans and their function. The Japanese teleost Fugu has been identified as one of the smallest vertebrate genomes and a compact model to study the human genome, owing to the great similarity in its gene repertoire with that of human and other vertebrates. Thus the characterization of the GPCRs of Fugu would provide insights to the evolution of the vertebrate genome.

Results

We classified the GPCRs in the Fugu genome and our analysis of its 316 membrane-bound receptors, available on the public databases as well as from literature, detected 298 GPCRs that were grouped into five main families according to the GRAFS classification system (namely, Glutamate, Rhodopsin, Adhesion, Frizzled and Secretin). We also identified 18 other GPCRs that could not be grouped under the GRAFS family and hence were classified as ‘Other 7TM’ receptors. On comparison of the GPCR information from the Fugu genome with those in the human and chicken genomes, we detected 96.83% (306/316) and 96.51% (305/316) orthology in GPCRs among the Fugu-human genomes and Fugu-chicken genomes, respectively.

Conclusions

This study reveals the position of pisces in vertebrate evolution from the GPCR perspective. Fugu can act as a reference model for the human genome for other protein families as well, going by the high orthology observed for GPCRs between Fugu and human. The evolutionary comparison of GPCR sequences between key vertebrate classes of mammals, birds and fish will help in identifying key functional residues and motifs so as to fill in the blanks in the evolution of GPCRs in vertebrates.

     

GPCRs heterodimerization: a new way towards the discovery of function for the orphan receptors?

G protein-coupled receptors (GPCRs), also called seven transmembrane domain (7TM) proteins, represent the largest family of cell surface receptors. GPCRs control a variety of physiological processes, are involved in multiple diseases and are major drug targets. Despite a vast effort of academic and industrial research, more than one hundred receptors remain orphans. These orphan GPCRs offer a great potential for drug discovery, as almost 60% of currently prescribed drugs target GPCRs. Deorphenization strategies have concentrated mainly on the identification of the natural ligands of these proteins. Recent advances have shown that orphan GPCRs, similar to orphan nuclear receptors, can regulate the function of non-orphan receptors by heterodimerization. These findings not only help to better understand the extraordinary diversity of GPCRs, but also open new perspectives for the identification of the function of these orphan receptors that hold great therapeutic potential.     

There exist at least 30 human G-protein-coupled receptors with long Ser/Thr-rich N-termini

We report six novel members of the superfamily of human G-protein coupled receptors (GPCRs) found by searches in the human genome databases, termed GPR123, GPR124, GPR125, GPR126, GPR127, and GPR128. Phylogenetic analysis demonstrates that these are additional members of the family of GPCRs with long N-termini, previously termed EGF-7TM, LNB-7TM, B2 or LN-7TM, showing that there exist at least 30 such GPCRs in the human genome. Three of these receptors form their own phylogenetic cluster, while two other places in a cluster with the previously reported HE6 and GPR56 (TM7XN1) and one with EMR1-3. All the novel receptors have a GPS domain in their N-terminus, except GPR123, as well as long Ser/Thr rich regions forming mucin-like stalks. GPR124 and GPR125 have a leucine rich repeat (LRR), an immunoglobulin (Ig) domain, and a hormone-binding domain (HBD). The Ig domain shows similarities to motilin and titin, while the LRR domain shows similarities to LRIG1 and SLIT1-2. GPR127 has one EGF domain while GPR126 and GPR128 do not contain domains that are readily recognized in other proteins beyond the GPS domain. We found several human EST sequences for most of the receptors showing differential expression patterns, which may indicate that some of these receptors participate in central functions while others are more likely to have a role in the immune or reproductive systems.     

The fourth transmembrane segment of the dopamine D2 receptor: accessibility in the binding-site crevice and position in the transmembrane bundle

The binding site of the dopamine D2 receptor, like that of homologous G-protein-coupled receptors (GPCRs), is contained within a water-accessible crevice formed among its seven transmembrane segments (TMSs). Using the substituted-cysteine-accessibility method (SCAM), we are mapping the residues that contribute to the surface of this binding-site crevice. We have mutated to cysteine, one at a time, 21 consecutive residues in the fourth TMS (TM4). Eleven of these mutants reacted with charged sulfhydryl-specific reagents, and bound antagonist protected nine of these from reaction. For the mutants in which cysteine was substituted for residues in the cytoplasmic half of TM4, treatment with the reagents had no effect on binding, consistent with these residues being inaccessible and with the low-resolution structure of the homologous rhodopsin, in which TM3 and TM5 occlude the cytoplasmic half of TM4. Although hydrophobicity analysis positions the C-terminus of TM4 at 4.64, Pro-Pro and Pro-X-Pro motifs, which are known to disrupt alpha-helices, occur at position 4.59 in a number of homologous GPCRs. The SCAM data were consistent with a C-terminus at 4.58, but it is also possible that the alpha-helix extends one additional turn to 4.62 in the D2 receptor, which has a single Pro at 4.59. In homologous GPCRs, the high degree of sequence variation between 4.59 and 4.68 is more characteristic of a loop domain than a helical segment. This region is shown here to be very conserved within functionally related receptors, suggesting an important functional role for this putative nonhelical domain. This inference is supported by observed ligand-specific effects of mutations in this region and by the predicted spatial proximity of this segment to known ligand binding sites in other TMs.     

3D modeling, ligand binding and activation studies of the cloned mouse delta, mu; and kappa opioid receptors

Refined 3D models of the transmembrane domains of the cloned delta, mu and kappa opioid receptors belonging to the superfamily of G-protein coupled receptors (GPCRs) were constructed from a multiple sequence alignment using the alpha carbon template of rhodopsin recently reported. Other key steps in the procedure were relaxation of the 3D helix bundle by unconstrained energy optimization and assessment of the stability of the structure by performing unconstrained molecular dynamics simulations of the energy optimized structure. The results were stable ligand-free models of the TM domains of the three opioid receptors. The ligand-free delta receptor was then used to develop a systematic and reliable procedure to identify and assess putative binding sites that would be suitable for similar investigation of the other two receptors and GPCRs in general. To this end, a non-selective, 'universal' antagonist, naltrexone, and agonist, etorphine, were used as probes. These ligands were first docked in all sites of the model delta opioid receptor which were sterically accessible and to which the protonated amine of the ligands could be anchored to a complementary proton-accepting residue. Using these criteria, nine ligand-receptor complexes with different binding pockets were identified and refined by energy minimization. The properties of all these possible ligand-substrate complexes were then examined for consistency with known experimental results of mutations in both opioid and other GPCRs. Using this procedure, the lowest energy agonist-receptor and antagonist-receptor complexes consistent with these experimental results were identified. These complexes were then used to probe the mechanism of receptor activation by identifying differences in receptor conformation between the agonist and the antagonist complex during unconstrained dynamics simulation. The results lent support to a possible activation mechanism of the mouse delta opioid receptor similar to that recently proposed for several other GPCRs. They also allowed the selection of candidate sites for future mutagenesis experiments.     

Novel human G protein-coupled receptors with long N-terminals containing GPS domains and Ser/Thr-rich regions

We report eight novel members of the superfamily of human G protein-coupled receptors (GPCRs) found by searches in the human genome databases, termed GPR97, GPR110, GPR111, GPR112, GPR113, GPR114, GPR115 and GPR116. Phylogenetic analysis shows that these are additional members of a family of GPCRs with long N-termini, previously termed EGF-7TM, LNB-7TM, B2 or LN-7TM. Five of the receptors form their own phylogenetic cluster, while three others form a cluster with the previously reported HE6 and GPR56 (TM7XN1). All the receptors have a GPS domain in their N-terminus and long Ser/Thr-rich regions forming mucin-like stalks. GPR113 has a hormone binding domain and one EGF domain. GPR112 has over 20 Ser/Thr repeats and a pentraxin domain. GPR116 has two immunoglobulin-like repeats and a SEA box. We found several human EST sequences for most of the receptors showing differential expression patterns, which may indicate that some of these receptors participate in reproductive functions while others are more likely to have a role in the immune system.     

Conserved activation pathways in G-protein-coupled receptors

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

Constitutive activation of angiotensin II type 1 receptor alters the orientation of transmembrane Helix-2

A key step in transmembrane (TM) signal transduction by G-protein-coupled receptors (GPCRs) is the ligand-induced conformational change of the receptor, which triggers the activation of a guanine nucleotide-binding protein. GPCRs contain a seven-TM helical structure essential for signal transduction in response to a large variety of sensory and hormonal signals. Primary structure comparison of GPCRs has shown that the second TM helix contains a highly conserved Asp residue, which is critical for agonist activation in these receptors. How conformational changes in TM2 relate to signal transduction by a GPCR is not known, because activation-induced conformational changes in TM2 helix have not been measured. Here we use modification of reporter cysteines to measure water accessibility at specific residues in TM2 of the type 1 receptor for the octapeptide hormone angiotensin II. Activation-dependent changes in the accessibility of Cys76 on TM2 were measured in constitutively activated mutants. These changes were directly correlated with measurement of function, establishing the link between physical changes in TM2 and function. Accessibility changes were measured at several consecutive residues on TM2, which suggest that TM2 undergoes a transmembrane movement in response to activation. This is the first report of in situ measurement of TM2 movement in a GPCR.     

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.     

Bioinformatics and type II G-protein-coupled receptors

The best known family B, or Type II, G-protein-coupled receptors (GPCRs) recognize peptides as ligands. The receptors for corticotrophin-releasing factor, parathyroid hormone and secretin typify this group. However, there are only 15 such GPCRs. Many other receptors share sequence homology and have been assigned to this family. The ten 'Frizzled' and one 'Smoothened' receptors show the lowest sequence homology and are not necessarily G-protein coupled. Drosophila genetics have enabled our understanding of their biology. In contrast, relatively little is known about the largest group with family B, the 33 'large amino termini' or large N-terminal family B seven-transmembrane (LNB 7TM) receptors. This review highlights the similarities found between family B receptors and provides a classification of LNB 7TM receptors.     

Structure-based identification of binding sites,native ligands and potential inhibitors for G-protein coupled receptors

G-protein coupled receptors (GPCRs) are the largest family of cell-surface receptors involved in signal transmission. Drugs associated with GPCRs represent more than one fourth of the 100 top-selling drugs and are the targets of more than half of the current therapeutic agents on the market. Our methodology based on the internal coordinate mechanics (ICM) program can accurately identify the ligand-binding pocket in the currently available crystal structures of seven transmembrane (7TM) proteins [bacteriorhodopsin (BR) and bovine rhodopsin (bRho)]. The binding geometry of the ligand can be accurately predicted by ICM flexible docking with and without the loop regions, a useful finding for GPCR docking because the transmembrane regions are easier to model. We also demonstrate that the native ligand can be identified by flexible docking and scoring in 1.5% and 0.2% (for bRho and BR, respectively) of the best scoring compounds from two different types of compound database. The same procedure can be applied to the database of available chemicals to identify specific GPCR binders. Finally, we demonstrate that even if the sidechain positions in the bRho binding pocket are entirely wrong, their correct conformation can be fully restored with high accuracy (0.28 A) through the ICM global optimization with and without the ligand present. These binding site adjustments are critical for flexible docking of new ligands to known structures or for docking to GPCR homology models. The ICM docking method has the potential to be used to "de-orphanize" orphan GPCRs (oGPCRs) and to identify antagonists-agonists for GPCRs if an accurate model (experimentally and computationally validated) of the structure has been constructed or when future crystal structures are determined.     

Prediction of the odorant binding site of olfactory receptor proteins by human-mouse comparisons

Olfactory receptors (ORs) are a large family of proteins involved in the recognition and discrimination of numerous odorants. These receptors belong to the G-protein coupled receptor (GPCR) hyperfamily, for which little structural data are available. In this study we predict the binding site residues of OR proteins by analyzing a set of 1441 OR protein sequences from mouse and human. The central insight utilized is that functional contact residues would be conserved among pairs of orthologous receptors, but considerably less conserved among paralogous pairs. Using judiciously selected subsets of 218 ortholog pairs and 518 paralog pairs, we have identified 22 sequence positions that are both highly conserved among the putative orthologs and variable among paralogs. These residues are disposed on transmembrane helices 2 to 7, and on the second extracellular loop of the receptor. Strikingly, although the prediction makes no assumption about the location of the binding site, these amino acid positions are clustered around a pocket in a structural homology model of ORs, mostly facing the inner lumen. We propose that the identified positions constitute the odorant binding site. This conclusion is supported by the observation that all but one of the predicted binding site residues correspond to ligand-contact positions in other rhodopsin-like GPCRs.     

Role of the extracellular amino terminus and first membrane-spanning helix of dopamine D1 and D5 receptors in shaping ligand selectivity and efficacy

Transmembrane (TM) helices of human D1-like dopaminergic receptors (hD1R and hD5R) harbor the same residues implicated in ligand binding and activation of catecholamine G protein-coupled receptors (GPCRs). Yet, hD1R and hD5R naturally display the distinct functional properties shared by wild type and constitutively active mutant GPCRs, respectively. Interestingly, we show in the present study that a class of synthetic phenylbenzazepine agonists containing a methyl on the azepine ring exhibited lower affinity for the more constitutively activated hD5R. These results cannot be explained by the “allosteric ternary complex model” postulating a higher agonist affinity for constitutively active GPCRs. We have also explored the functional role of distinct extracellular amino terminus (NT) and TM1 regions of hD1R and hD5R using a chimerical approach. Of these two regions, our studies suggest that TM1 predominantly shapes D1-like ligand affinity and selectivity. Additionally, NT and TM1 of hD1R and hD5R play no role in receptor constitutive activity but differentially modulate dopamine-mediated responsiveness. The TM1 exchange mediated drastic changes in intrinsic efficacy and activity of phenylbenzazepine drugs displaying partial agonism at hD1R and hD5R. Phenylbenzazepines were converted into strong partial agonists or full agonists in cells expressing hD1R-TM1^D5 chimera while being switched from full agonists to partial agonists and partial agonists to antagonists in cells harboring hD5R-TM1^D1 chimera. TM1 exchange had no effect on antipsychotic-mediated inverse agonism. In summary, our study shows that NT and TM1 of D1-like receptors control ligand binding and agonist-induced activation, poising these regions as important structural determinants for catecholamine GPCR function.     

Random mutagenesis of the M3 muscarinic acetylcholine receptor expressed in yeast: identification of second-site mutations that restore function to a coupling-deficient mutant M3 receptor

The M(3) muscarinic receptor is a prototypical member of the class A family of G protein-coupled receptors (GPCRs). To gain insight into the structural mechanisms governing agonist-mediated M(3) receptor activation, we recently developed a genetically modified yeast strain (Saccharomyces cerevisiae) which allows the efficient screening of large libraries of mutant M(3) receptors to identify mutant receptors with altered/novel functional properties. Class A GPCRs contain a highly conserved Asp residue located in transmembrane domain II (TM II; corresponding to Asp-113 in the rat M(3) muscarinic receptor) which is of fundamental importance for receptor activation. As observed previously with other GPCRs analyzed in mammalian expression systems, the D113N point mutation abolished agonist-induced receptor/G protein coupling in yeast. We then subjected the D113N mutant M(3) receptor to PCR-based random mutagenesis followed by a yeast genetic screen to recover point mutations that can restore G protein coupling to the D113N mutant receptor. A large scale screening effort led to the identification of three such second-site suppressor mutations, R165W, R165M, and Y250D. When expressed in the wild-type receptor background, these three point mutations did not lead to an increase in basal activity and reduced the efficiency of receptor/G protein coupling. Similar results were obtained when the various mutant receptors were expressed and analyzed in transfected mammalian cells (COS-7 cells). Interestingly, like Asp-113, Arg-165 and Tyr-250, which are located at the cytoplasmic ends of TM III and TM V, respectively, are also highly conserved among class A GPCRs. Our data suggest a conformational link between the highly conserved Asp-113, Arg-165, and Tyr-250 residues which is critical for receptor activation.