C5a receptor activation. Genetic identification of critical residues in four transmembrane helices.

Hormones and sensory stimuli activate serpentine receptors, transmembrane switches that relay signals to heterotrimeric guanine nucleotide-binding proteins (G proteins). To understand the switch mechanism, we subjected 93 amino acids in transmembrane helices III, V, VI, and VII of the human chemoattractant C5a receptor to random saturation mutagenesis. A yeast selection identified 121 functioning mutant receptors, containing a total of 523 amino acid substitutions. Conserved hydrophobic residues are located on helix surfaces that face other helices in a modeled seven-helix bundle (Baldwin, J. M., Schertler, G. F., and Unger, V. M. (1997) J. Mol. Biol. 272, 144-164), whereas surfaces predicted to contact the surrounding lipid tolerate many substitutions. Our analysis identified 25 amino acid positions resistant to nonconservative substitutions. These appear to comprise two distinct components of the receptor switch, a surface at or near the extracellular membrane interface and a core cluster in the cytoplasmic half of the bundle. Twenty-one of the 121 mutant receptors exhibit constitutive activity. Amino acids substitutions in these activated receptors predominate in helices III and VI; other activating mutations truncate the receptor near the extracellular end of helix VI. These results identify key elements of a general mechanism for the serpentine receptor switch.     

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

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

FoldGPCR: Structure prediction protocol for the transmembrane domain of G protein‐coupled receptors from class A

Building reliable structural models of G protein‐coupled receptors (GPCRs) is a difficult task because of the paucity of suitable templates, low sequence identity, and the wide variety of ligand specificities within the superfamily. Template‐based modeling is known to be the most successful method for protein structure prediction. However, refinement of homology models within 1–3 Å Cα RMSD of the native structure remains a major challenge. Here, we address this problem by developing a novel protocol (foldGPCR) for modeling the transmembrane (TM) region of GPCRs in complex with a ligand, aimed to accurately model the structural divergence between the template and target in the TM helices. The protocol is based on predicted conserved inter‐residue contacts between the template and target, and exploits an all‐atom implicit membrane force field. The placement of the ligand in the binding pocket is guided by biochemical data. The foldGPCR protocol is implemented by a stepwise hierarchical approach, in which the TM helical bundle and the ligand are assembled by simulated annealing trials in the first step, and the receptor‐ligand complex is refined with replica exchange sampling in the second step. The protocol is applied to model the human β₂‐adrenergic receptor (β₂AR) bound to carazolol, using contacts derived from the template structure of bovine rhodopsin. Comparison with the X‐ray crystal structure of the β₂AR shows that our protocol is particularly successful in accurately capturing helix backbone irregularities and helix‐helix packing interactions that distinguish rhodopsin from β₂AR. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Structure and pH sensitivity of the transmembrane segment 3 of rhodopsin

Activation of G protein-coupled receptors (GPCRs) originates in ligand-induced protein conformational changes that are transmitted to the cytosolic receptor surface. In the photoreceptor rhodopsin, and possibly other rhodopsin-like GPCRs, protonation of a carboxylic acid in the conserved E(D)RY motif at the cytosolic end of transmembrane helix 3 (TM3) is coupled to receptor activation. Here, we have investigated the structure of synthetic peptides derived from rhodopsin TM3. Polarized FTIR spectroscopy reveals a helical structure of a 31-mer TM3 peptide reconstituted into PC vesicles with a large tilt of 40-50 degrees of the helical axis relative to the membrane normal. Helical structure is also observed for the TM3 peptide in detergent micelles and depends on pH, especially in the C-terminal sequence. In addition, the fluorescence emission of the single tyrosine of the D(E)RY motif in the TM3 peptide exhibits a pronounced pH sensitivity that is abolished when Glu is replaced by Gln, demonstrating that protonation of the conserved Glu side chain affects the structure in the environment of the D(E)RY motif of TM3. The pH regulation of the C-terminal TM3 structure may be an intrinsic feature of the E(D)RY motif in other class I receptors, allowing the coupling of protonation and conformation of membrane-exposed residues in full-length GPCRs.     

The first 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 (TMs). Using the substituted-cysteine-accessibility method (SCAM), we are mapping the residues that contribute to the surface of this binding-site crevice. We have now mutated to cysteine, one at a time, 21 consecutive residues in TM1. Six of these mutants reacted with charged sulfhydryl reagents, whereas bound antagonist only protected N52(1.50)C from reaction. Except for A38(1.36)C, none of the substituted cysteine mutants in the extracellular half of TM1 appeared to be accessible. Pro(1.48) is highly conserved in opsins, but absent in catecholamine receptors, and the high-resolution rhodopsin structure showed that Pro(1.48) bends the extracellular portion of TM1 inward toward TM2 and TM7. Analysis of the conversation of residues in the extracellular portion of TM1 of opsins showed a pattern consistent with alpha-helical structure with a conserved face. In contrast, this region in catecholamine receptors is poorly conserved, suggesting a lack of critical contacts. Thus, in catecholamine receptors in the absence of Pro(1.48), TM1 may be straighter and therefore further from the helix bundle, consistent with the apparent lack of conserved contact residues. When examined in the context of a model of the D2 receptor, the accessible residues in the cytoplasmic half of TM1 are at the interface with TM7 and with helix 8 (H8). We propose the existence of critical contacts of TM1, TM7, and H8 that may stabilize the inactive state of the receptor.     

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.     

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.     

Exploring the structure of opioid receptors with homology modeling based on single and multiple templates and subsequent docking: A comparative study

Opioid receptors are the principal targets for opioids, which have been used as analgesics for centuries. Opioid receptors belong to the rhodopsin family of G-protein coupled receptors (GPCRs). In the absence of crystal structures of opioid receptors, 3D homology models have been reported with bovine rhodopsin as a template, though the sequence homology is low. Recently, it has been reported that use of multiple templates results in a better model for a target having low sequence identity with a single template. With the objective of carrying out a comparative study on the structural quality of the 3D models based on single and multiple templates, the homology models for opioid receptors (mu, delta and kappa) were generated using bovine rhodopsin as single template and the recently deposited crystal structures of squid rhodopsin, turkey β-1 and human β-2 adrenoreceptors along with bovine rhodopsin as multiple templates. In this paper we report the results of comparison between the refined 3D models based on multiple sequence alignment (MSA) and models built with bovine rhodopsin as template, using validation programs PROCHECK, PROSA, Verify 3D, Molprobity and docking studies. The results indicate that homology models of mu and kappa with multiple templates are better than those built with only bovine rhodopsin as template, whereas, in many aspects, the homology model of delta opioid receptor with single template is better with respect to the model based on multiple templates. Three nonselective ligands were docked to both the models of mu, delta and kappa opioid receptors using GOLD 3.1. The results of docking complied well with the pharamacophore, reported for nonspecific opioid ligands. The comparison of docking results for models with multiple templates and those with single template have been discussed in detail. Three selective ligands for each receptor were also docked. As the crystallographic structures are not yet known, this comparison will help in choosing better homology models of opioid receptors for studying ligand receptor interactions to design new potent opioid antagonists.     

Assembly of transmembrane helices of simple polytopic membrane proteins from sequence conservation patterns

The transmembrane (TM) domains of most membrane proteins consist of helix bundles. The seemingly simple task of TM helix bundle assembly has turned out to be extremely difficult. This is true even for simple TM helix bundle proteins, i.e., those that have the simple form of compact TM helix bundles. Herein, we present a computational method that is capable of generating native-like structural models for simple TM helix bundle proteins having modest numbers of TM helices based on sequence conservation patterns. Thus, the only requirement for our method is the presence of more than 30 homologous sequences for an accurate extraction of sequence conservation patterns. The prediction method first computes a number of representative well-packed conformations for each pair of contacting TM helices, and then a library of tertiary folds is generated by overlaying overlapping TM helices of the representative conformations. This library is scored using sequence conservation patterns, and a subsequent clustering analysis yields five final models. Assuming that neighboring TM helices in the sequence contact each other (but not that TM helices A and G contact each other), the method produced structural models of Calpha atom root-mean-square deviation (CA RMSD) of 3-5 A from corresponding crystal structures for bacteriorhodopsin, halorhodopsin, sensory rhodopsin II, and rhodopsin. In blind predictions, this type of contact knowledge is not available. Mimicking this, predictions were made for the rotor of the V-type Na(+)-adenosine triphosphatase without such knowledge. The CA RMSD between the best model and its crystal structure is only 3.4 A, and its contact accuracy reaches 55%. Furthermore, the model correctly identifies the binding pocket for sodium ion. These results demonstrate that the method can be readily applied to ab initio structure prediction of simple TM helix bundle proteins having modest numbers of TM helices.     

Use of an in situ disulfide cross-linking strategy to study the dynamic properties of the cytoplasmic end of transmembrane domain VI of the M3 muscarinic acetylcholine receptor

The ligand-induced activation of G protein-coupled receptors (GPCRs) is predicted to involve pronounced conformational changes on the intracellular surface or the receptor proteins. A reorientation of the cytoplasmic end of transmembrane domain VI (TM VI) is thought to play a key role in GPCR activation and productive receptor/G protein coupling. Disulfide cross-linking studies with solubilized, Cys-substituted mutant versions of bovine rhodopsin and the M3 muscarinic acetylcholine receptor suggested that the cytoplasmic end of TM VI is conformationally highly flexible, even in the absence of activating ligands (Farrens, D. L., et al. (1996) Science 274, 768-770; Zeng, F. Y., et al. (1999) J. Biol. Chem. 274, 16629-16640). To test the hypothesis that the promiscuous disulfide cross-linking pattern observed in these studies was caused by the use of solubilized receptor proteins endowed with increased conformational flexibility, we employed a recently developed in situ disulfide cross-linking strategy that allows the detection of disulfide bonds in Cys-substituted mutant M3 muscarinic receptors present in their native membrane environment. Specifically, we used membranes prepared from transfected COS-7 cells to analyze a series of double Cys mutant M3 receptors containing one Cys residue within the sequence K484(6.29) to S493(6.38) at the cytoplasmic end of TM VI and a second Cys residue at the cytoplasmic end of TM III (I169C(3.54)). This analysis revealed a disulfide cross-linking pattern that was strikingly more restricted than that observed previously with solubilized receptor proteins, both in the absence and in the presence of the muscarinic agonist, carbachol. Carbachol stimulated the formation of disulfide bonds in only two of the 10 analyzed mutant muscarinic receptors, I169C(3.54)/K484C(6.29) and I169C(3.54)/A488C(6.33), consistent with an agonist-induced rotation of the cytoplasmic end of TM VI. These findings underline the usefulness of analyzing the structural and dynamic properties of GPCRs in their native lipid environment.     

Tertiary interactions between the fifth and sixth transmembrane segments of rhodopsin.

We have used cysteine scanning mutagenesis and disulfide cross-linking in a split rhodopsin construct to investigate the secondary structure and tertiary contacts of the fifth (TM5) and sixth (TM6) transmembrane segments of rhodopsin. Using a simple increase in pH to promote disulfide bond formation, three cross-links between residues on the extracellular side of TM5 (at positions 198, 200, and 204) and TM6 (at position 276) have been identified and characterized. The helical pattern of cross-linking observed indicates that the fifth transmembrane helix extends through residue 200 but does not include residue 198. Rhodopsin mutants containing these disulfides demonstrate nativelike absorption spectra and light-dependent activation of transducin, suggesting that large movements on the extracellular side of TM5 with respect to TM6 are not required for receptor activation.     

Unique features of the structural model of 'hard' cuticle proteins: implications for chitin-protein interactions and cross-linking in cuticle

Cuticular proteins are one of the determinants of the physical properties of cuticle. A common consensus region (extended R&R Consensus) in these proteins binds to chitin, the other major component of cuticle. We previously predicted the preponderance of beta-pleated sheet in the consensus region and proposed its responsibility for the formation of helicoidal cuticle (Iconomidou et al., Insect Biochem. Mol. Biol. 29 (1999) 285). Subsequently, we verified experimentally the abundance of antiparallel beta-pleated sheet in the structure of cuticle proteins (Iconomidou et al., Insect Biochem. Mol. Biol. 31 (2001) 877). Homology modelling of soft (RR-1) cuticular proteins using bovine plasma retinol binding protein (RBP) as a template revealed an antiparallel beta-sheet half-barrel structure as the basic folding motif (Hamodrakas et al., Insect Biochem. Molec. Biol. 32 (2002) 1577). The RR-2 proteins characteristic of hard cuticle, have a far more conserved consensus and frequently more histidine residues. Extension of modelling to this class of consensus, in this work, reveals in detail several unique features of the proposed structural model to serve as a chitin binding structural motif, thus providing the basis for elucidating cuticle's overall architecture and chitin-protein interactions in cuticle.     

Improved model building and assessment of the Calcium-sensing receptor transmembrane domain

A three-dimensional model of the human Calcium-sensing receptor (CaSR) seven transmembrane domain was built via a novel sequence alignment method based on the conserved contacts in proteins using the crystal structure of bovine rhodopsin as the template. This model was tested by docking NPS 2143, the first identified allosteric antagonist of CaSR. In our model, Glu837 plays a critical role in anchoring the protonated nitrogen atom and hydroxy oxygen atom of NPS 2143. The phenyl moiety of the ligand contacts residues Phe668, Pro672, and Ile841. The naphthalene moiety is surrounded by several hydrophobic residues, including Phe684, Phe688, and Phe821. Our model appears to be consistent with all six residues that have been demonstrated to be critical for NPS 2143 binding, in contrast with existing homology models based on traditional sequence alignment of CaSR to rhodopsin. This provides validation of our sequence alignment method and the use of the rhodopsin backbone as the initial structure in homology modeling of other G protein-coupled receptors that are not members of the rhodopsin family.     

Structures of the transmembrane helices of the G-protein coupled receptor, rhodopsin.

An hypothesis is tested that individual peptides corresponding to the transmembrane helices of the membrane protein, rhodopsin, would form helices in solution similar to those in the native protein. Peptides containing the sequences of helices 1, 4 and 5 of rhodopsin were synthesized. Two peptides, with overlapping sequences at their termini, were synthesized to cover each of the helices. The peptides from helix 1 and helix 4 were helical throughout most of their length. The N- and C-termini of all the peptides were disordered and proline caused opening of the helical structure in both helix 1 and helix 4. The peptides from helix 5 were helical in the middle segment of each peptide, with larger disordered regions in the N- and C-termini than for helices 1 and 4. These observations show that there is a strong helical propensity in the amino acid sequences corresponding to the transmembrane domain of this G-protein coupled receptor. In the case of the peptides from helix 4, it was possible to superimpose the structures of the overlapping sequences to produce a construct covering the whole of the sequence of helix 4 of rhodopsin. As similar superposition for the peptides from helix 1 also produced a construct, but somewhat less successfully because of the disordering in the region of sequence overlap. This latter problem was more severe for helix 5 and therefore a single peptide was synthesized for the entire sequence of this helix, and its structure determined. It proved to be helical throughout. Comparison of all these structures with the recent crystal structure of rhodopsin revealed that the peptide structures mimicked the structures seen in the whole protein. Thus similar studies of peptides may provide useful information on the secondary structure of other transmembrane proteins built around helical bundles.     

Comparative model building of interleukin-7 using interleukin-4 as a template: a structural hypothesis that displays atypical surface chemistry in helix D important for receptor activation

Using a combination of theoretical sequence structure recognition predictions and experimental disulfide bond assignments, a three-dimensional (3D) model of human interleukin-7 (hIL-7) was constructed that predicts atypical surface chemistry in helix D that is important for receptor activation. A 3D model of hIL-7 was built using the X-ray crystal structure of interleukin-4 (IL-4) as a template (Walter MR et al., 1992, J Mol Biol. 224:1075-1085; Walter MR et al., 1992, J Biol Chem 267:20371-20376). Core secondary structures were constructed from sequences of hIL-7 predicted to form helices. The model was constructed by superimposing IL-7 helices onto the IL-4 template and connecting them together in an up-up down-down topology. The model was finished by incorporating the disulfide bond assignments (Cys3, Cys142), (Cys35, Cys130), and (Cys48, Cys93), which were determined by MALDI mass spectroscopy and site-directed mutagenesis (Cosenza L, Sweeney E, Murphy JR, 1997, J Biol Chem 272:32995-33000). Quality analysis of the hIL-7 model identified poor structural features in the carboxyl terminus that, when further studied using hydrophobic moment analysis, detected an atypical structural property in helix D, which contains Cys 130 and Cys142. This analysis demonstrated that helix D had a hydrophobic surface exposed to bulk solvent that accounted for the poor quality of the model, but was suggestive of a region in IL-7 that maybe important for protein interactions. Alanine (Ala) substitution scanning mutagenesis was performed to test if the predicted atypical surface chemistry of helix D in the hIL-7 model is important for receptor activation. This analysis resulted in the construction, purification, and characterization of four hIL-7 variants, hIL-7(K121A), hIL-7(L136A), hIL-7(K140A), and hIL-7(W143A), that displayed reduced or abrogated ability to stimulate a murine IL-7 dependent pre-B cell proliferation. The mutant hIL-7(W143A), which is biologically inactive and displaces [125I]-hIL-7, is the first reported IL-7R system antagonist.     

A conserved Asn in transmembrane helix 7 is an on/off switch in the activation of the thyrotropin receptor

The thyrotropin (TSH) receptor is an interesting model to study G protein-coupled receptor activation as many point mutations can significantly increase its basal activity. Here, we identified a molecular interaction between Asp(633) in transmembrane helix 6 (TM6) and Asn(674) in TM7 of the TSHr that is crucial to maintain the inactive state through conformational constraint of the Asn. We show that these residues are perfectly conserved in the glycohormone receptor family, except in one case, where they are exchanged, suggesting a direct interaction. Molecular modeling of the TSHr, based on the high resolution structure of rhodopsin, strongly favors this hypothesis. Our approach combining site-directed mutagenesis with molecular modeling shows that mutations disrupting this interaction, like the D633A mutation in TM6, lead to high constitutive activation. The strongly activating N674D (TM7) mutation, which in our modeling breaks the TM6-TM7 link, is reverted to wild type-like behavior by an additional D633N mutation (TM6), which would restore this link. Moreover, we show that the Asn of TM7 (conserved in most G protein-coupled receptors) is mandatory for ligand-induced cAMP accumulation, suggesting an active role of this residue in activation. In the TSHr, the conformation of this Asn residue of TM7 would be constrained, in the inactive state, by its Asp partner in TM6.     

Structure and function in rhodopsin: mapping light-dependent changes in distance between residue 65 in helix TM1 and residues in the sequence 306-319 at the cytoplasmic end of helix TM7 and in helix H8.

Spin-labeled double mutants of rhodopsin were produced containing a reference nitroxide at position 65, at the cytoplasmic termination of helix TM1, and a second nitroxide in the sequence of residues 306-319, which includes the cytoplasmic termination of helix TM7 and nearly the entire surface helix H8. Magnetic dipole-dipole interactions between the spins are analyzed to provide interspin distance distributions in both the dark and photoactivated states of rhodopsin. The distributions, apparently resulting from the conformational flexibility of the side chains, are found to be consistent with the structural model of rhodopsin in the dark state derived from crystallography. Photoactivation of the receptor triggers an increase in distance between residues in TM7, but not those in H8, relative to the reference at position 65 in TM1. The simplest interpretation of the result is a movement of the cytoplasmic portion of TM7 away from TM1 by 2-4 A.     

Comparison of the amino acid residues in the sixth transmembrane domains accessible in the binding-site crevices of mu, delta, and kappa opioid receptors.

We have mapped the residues in the sixth transmembrane domains (TMs 6) of the mu, delta, and kappa opioid receptors that are accessible in the binding-site crevices by the substituted cysteine accessibility method (SCAM). We previously showed that ligand binding to the C7.38S mutants of the mu and kappa receptors and the wild-type delta receptor was relatively insensitive to methanethiosulfonate ethylammonium (MTSEA), a positively charged sulfhydryl-specific reagent. These MTSEA-insensitive constructs were used as the templates, and 22 consecutive residues in TM6 (excluding C6.47) of each receptor were mutated to cysteine, 1 at a time. Most mutants retained binding affinities for [3H]diprenorphine, a nonselective opioid antagonist, similar to that of the template receptors. Treatment with MTSEA significantly inhibited [3H]diprenorphine binding to 11 of 22 mutants of the rat mu receptor and 9 of 22 mutants of the human delta receptor and 10 of 22 mutants of the human kappa receptor. Naloxone or diprenorphine protected all sensitive mutants, except the A6.42(287)C mu mutant. Thus, V6.40, F6.44, W6.48, I6.51, Y6.54, V6.55, I6.56, I6.57, K6.58, and A6.59 of the mu receptor; F6.44, I6.51, F6.54, V6.55, I6.56, V6.57, W6.58, T6.59, and L6.60 of the delta receptor; and F6.44, W6.48, I6.51, F6.54, I6.55, L6.56, V6.57, E6.58, A6.59, and L6.60 of the kappa receptor are on the water-accessible surface of the binding-site crevices. The accessibility patterns of residues in the TMs 6 of the mu, delta, and kappa opioid receptors are consistent with the notion that the TMs 6 are in alpha-helical conformations with a narrow strip of accessibility on the intracellular side of 6.54 and a wider area of accessibility on the extracellular side of 6.54, likely due to a proline kink at 6.50 that bends the helix in toward the binding pocket and enables considerable motion in this region. The wider exposure of residues 6.55-6.60 to the binding-site crevice, combined with the divergent amino acid sequences, is consistent with the inferred role of residues in this region in determining ligand binding selectivity. The conservation of the accessibility pattern on the cytoplasmic side of 6.54 suggests that this region may be important for receptor activation. This accessibility pattern is similar to that of the D2 dopamine receptor, the only other GPCR in which TM6 has been mapped by SCAM. That opioid receptors and the remotely related D2 dopamine receptor have similar accessibility patterns in TM6 suggest that these segments of GPCRs in the rhodopsin-like subfamily not only share secondary structure but also are packed similarly into the transmembrane bundle and thus have similar tertiary structure.     

Critical role of transmembrane segment zinc binding in the structure and function of rhodopsin

Zinc deficiency and retinitis pigmentosa are both important factors resulting in retinal dysfunction and night blindness. In this study, we address the critical biochemical and structural relevance of zinc ions in rhodopsin and examine whether zinc deficiency can lead to rhodopsin dysfunction. We report the identification of a high-affinity zinc coordination site within the transmembrane domain of rhodopsin, coordinated by the side chains of two highly conserved residues, Glu(122) in transmembrane helix III and His(211) in transmembrane helix V. We also demonstrate that this zinc coordination is critical for rhodopsin folding, 11-cis-retinal binding, and the stability of the chromophore-receptor interaction, defects of which are observed in retinitis pigmentosa. Furthermore, a cluster of retinitis pigmentosa mutations is localized within and around this zinc binding site. Based on these studies, we believe that improvement in zinc binding to rhodopsin at this site within the transmembrane domain may be a pharmacological approach for the treatment of select retinitis pigmentosa mutations. Transmembrane coordination of zinc may also be an important common principle across G-protein-coupled receptors.     

NMR structure and dynamics of a designed water-soluble transmembrane domain of nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is an important therapeutic target for a wide range of pathophysiological conditions, for which rational drug designs often require receptor structures at atomic resolution. Recent proof-of-concept studies demonstrated a water-solubilization approach to structure determination of membrane proteins by NMR (Slovic et al., PNAS, 101: 1828-1833, 2004; Ma et al., PNAS, 105: 16537-42, 2008). We report here the computational design and experimental characterization of WSA, a water-soluble protein with ~83% sequence identity to the transmembrane (TM) domain of the nAChR α1 subunit. Although the design was based on a low-resolution structural template, the resulting high-resolution NMR structure agrees remarkably well with the recent crystal structure of the TM domains of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), demonstrating the robustness and general applicability of the approach. NMR T(2) dispersion measurements showed that the TM2 domain of the designed protein was dynamic, undergoing conformational exchange on the NMR timescale. Photoaffinity labeling with isoflurane and propofol photolabels identified a common binding site in the immediate proximity of the anesthetic binding site found in the crystal structure of the anesthetic-GLIC complex. Our results illustrate the usefulness of high-resolution NMR analyses of water-solubilized channel proteins for the discovery of potential drug binding sites.