Combined biophysical and biochemical information confirms arrangement of transmembrane helices visible from the three-dimensional map of frog rhodopsin

Unconventional Roles of Opsins

Rhodopsin is the classical light sensor. Although rhodopsin has long been known to be important for image formation in the eye, the requirements for opsins in non-image formation and in extraocular light sensation were revealed much later. Most recent is the demonstration that an opsin in the fruit fly, Drosophila melanogaster, is expressed in pacemaker neurons in the brain and functions in light entrainment of circadian rhythms. However, the biggest surprise is that opsins have light-independent roles, countering more than a century of dogma that they function exclusively as light sensors. Through studies in Drosophila, light-independent roles of opsins have emerged in temperature sensation and hearing. Although these findings have been uncovered in the fruit fly, there are hints that opsins have light-independent roles in a wide array of animals, including mammals. Thus, despite the decades of focus on opsins as light detectors, they represent an important new class of polymodal sensory receptor.

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 A C alpha 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 beta(2)-adrenergic receptor (beta(2)AR) bound to carazolol, using contacts derived from the template structure of bovine rhodopsin. Comparison with the X-ray crystal structure of the beta(2)AR shows that our protocol is particularly successful in accurately capturing helix backbone irregularities and helix-helix packing interactions that distinguish rhodopsin from beta(2)AR.

Critical review of general guidelines for membrane proteins model building and analysis

Membrane proteins play major roles in many biological processes such as signalling, transport, etc. They have been shown to be involved in the development of many diseases and have become important drug targets per se. The understanding of their functional properties may be facilitated if a 3D structure is available. However, in the case of membrane proteins, only a few 3D structures have been solved to date. Bioinformatics and molecular modelling approaches are thus powerful alternatives to fill the gap between the sequence and the structure. Here, a review of the most recent approaches is proposed together with guidelines on how to use them. In addition, insofar as important biological processes require conformational changes, we discuss some interesting methods aimed at exploring the dynamic behaviour of proteins in their membrane environment. The paper ends with a brief description of useful approaches for determining oligomerisation or ligand binding sites.

Structural refinement of membrane proteins by restrained molecular dynamics and solvent accessibility data

We present an approach for incorporating solvent accessibility data from electron paramagnetic resonance experiments in the structural refinement of membrane proteins through restrained molecular dynamics simulations. The restraints have been parameterized from oxygen (PiO(2)) and nickel-ethylenediaminediacetic acid (PiNiEdda) collision frequencies, as indicators of lipid or aqueous exposed spin-label sites. These are enforced through interactions between a pseudoatom representation of the covalently attached Nitroxide spin-label and virtual "solvent" particles corresponding to O(2) and NiEdda in the surrounding environment. Interactions were computed using an empirical potential function, where the parameters have been optimized to account for the different accessibilities of the spin-label pseudoatoms to the surrounding environment. This approach, "pseudoatom-driven solvent accessibility refinement", was validated by refolding distorted conformations of the Streptomyces lividans potassium channel (KcsA), corresponding to a range of 2-30 A root mean-square deviations away from the native structure. Molecular dynamics simulations based on up to 58 electron paramagnetic resonance restraints derived from spin-label mutants were able to converge toward the native structure within 1-3 A root mean-square deviations with minimal computational cost. The use of energy-based ranking and structure similarity clustering as selection criteria helped in the convergence and identification of correctly folded structures from a large number of simulations. This approach can be applied to a variety of integral membrane protein systems, regardless of oligomeric state, and should be particularly useful in calculating conformational changes from a known reference crystal structure.

Packing of transmembrane helices in bacteriorhodopsin folding: structure and thermodynamics

We propose a coarse-grained (CG) model to study the native structure and physical properties of helical membrane proteins (HMPs) using off-lattice computer simulations. Instead of considering sequence heterogeneity explicitly, we model its effect on the packing of helices by employing a mean packing parameter r(0), which is calculated from an all-atom (AA) model. Specifically, this CG model is applied to investigate the packing of helices in bacteriorhodopsin (BR), and predicts the seven helix bundle structure of BR with a root mean square deviation (RMSD) in coordinates of helix backbone atoms (N, C, C(alpha)) of 3.99 A from its crystal structure. This predicted structure is further refined in an AA model by Amber and the refined structure has a RMSD (in coordinates of helix backbone atoms) of 2.64 A. The predicted packing position, tilting angle, and orientation angle of each helix in the refined structure are consistent with experimental data and their physical origins can be well understood in our model. Our results show that a reasonably good structure of BR can be predicted by using such a dual-scale approach, provided that its secondary structure is known. Starting from a random initial configuration, the folded structure can be obtained in days using a regular desktop computer. Various thermodynamic properties of helix packing of BR are also investigated in this CG model.

Method To Assess Packing Quality of Transmembrane α-Helices in Proteins. 1. Parametrization Using Structural Data

Integral membrane proteins (MPs) are pharmaceutical targets of exceptional importance. Modern methods of three-dimensional protein structure determination often fail to supply the fast growing field of structure-based drug design with the requested MPs' structures. That is why computational modeling techniques gain a special importance for these objects. Among the principal difficulties limiting application of these methods is the low quality of the MPs' models built in silico. In this series of two papers we present a computational approach to the assessment of the packing "quality" of transmembrane (TM) alpha-helical domains in proteins. The method is based on the concept of protein environment classes, whereby each amino acid residue is described in terms of its environment polarity and accessibility to the membrane. In the first paper we analyze a nonredundant set of 26 TM alpha-helical domains and compute the residues' propensities to five predefined classes of membrane-protein environments. Here we evaluate the proposed approach only by various test sets, cross-validation protocols and ability of the method to delimit the crystal structure of visual rhodopsin, and a number of its erroneous theoretical models. More advanced validation of the method is given in the second article of this series. We assume that the developed "membrane score" method will be helpful in optimizing computer models of TM domains of MPs, especially G-protein coupled receptors.

Modeling GPCRs

Many GPCR models have been built over the years for many different purposes, of which drug-design undoubtedly has been the most frequent one. The release of the structure of bovine rhodopsin in August 2000 enabled us to analyze models built before that period to learn things for the models we build today. We conclude that the GPCR modeling field is riddled with "common knowledge". Several characteristics of the bovine rhodopsin structure came as a big surprise, and had obviously not been predicted, which led to large errors in the models. Some of these surprises, however, could have been predicted if the modelers had more rigidly stuck to the rule that holds for all models, namely that a model should explain all experimental facts, and not just those facts that agree with the modeler's preconceptions.

G-protein coupled receptor structure

Because of their central role in regulation of cellular function, structure/function relationships for G-protein coupled receptors (GPCR) are of vital importance, yet only recently have sufficient data been obtained to begin mapping those relationships. GPCRs regulate a wide range of cellular processes, including the senses of taste, smell, and vision, and control a myriad of intracellular signaling systems in response to external stimuli. Many diseases are linked to GPCRs. A critical need exists for structural information to inform studies on mechanism of receptor action and regulation. X-ray crystal structures of only one GPCR, in an inactive state, have been obtained to date. However considerable structural information for a variety of GPCRs has been obtained using non-crystallographic approaches. This review begins with a review of the very earliest GPCR structural information, mostly derived from rhodopsin. Because of the difficulty in crystallizing GPCRs for X-ray crystallography, the extensive published work utilizing alternative approaches to GPCR structure is reviewed, including determination of three-dimensional structure from sparse constraints. The available X-ray crystallographic analyses on bovine rhodopsin are reviewed as the only available high-resolution structures for any GPCR. Structural information available on ligand binding to several receptors is included. The limited information on excited states of receptors is also reviewed. It is concluded that while considerable basic structural information has been obtained, more data are needed to describe the molecular mechanism of activation of a GPCR.

A strategy using NMR peptide structures of thromboxane A2 receptor as templates to construct ligand-recognition pocket of prostacyclin receptor

Background:

Prostacyclin receptor (IP) and thromboxane A2 receptor (TP) belong to rhodopsin-type G protein-coupling receptors and respectively bind to prostacyclin and thromboxane A2 derived from arachidonic acid. Recently, we have determined the extracellular loop (eLP) structures of the human TP receptor by 2-D 1H NMR spectroscopy using constrained peptides mimicking the individual eLP segments. The studies have identified the segment along with several residues in the eLP domains important to ligand recognition, as well as proposed a ligand recognition pocket for the TP receptor.

Results:

The IP receptor shares a similar primary structure in the eLPs with those of the TP receptor. Forty percent residues in the second eLPs of the receptors are identical, which is the major region involved in forming the ligand recognition pocket in the TP receptor. Based on the high homology score, the eLP domains of the IP receptor were constructed by the homology modeling approach using the NMR structures of the TP eLPs as templates, and then configured to the seven transmembrane (TM) domains model constructed using the crystal structure of the bovine rhodopsin as a template. A NMR structure of iloprost was docked into the modeled IP ligand recognition pocket. After dynamic studies, the segments and residues involved in the IP ligand recognition were proposed. A key residue, Arg173 involved in the ligand recognition for the IP receptor, as predicted from the modeling, was confirmed by site-directed mutagenesis.

Conclusion:

A 3-D model of the human IP receptor was constructed by homology modeling using the crystal structure of bovine rhodopsin TM domains and the NMR structures of the synthetic constrained peptides of the eLP domains of the TP receptor as templates. This strategy can be applied to molecular modeling and the prediction of ligand recognition pockets for other prostanoid receptors.

Biosynthesis and NMR Analysis of a 73-Residue Domain of a Saccharomyces cerevisiae G Protein-Coupled Receptor

The yeast Saccharomyces cerevisiae alpha-factor pheromone receptor (Ste2p) was used as a model G protein-coupled receptor (GPCR). A 73-mer multidomain fragment of Ste2p (residues 267-339) containing the third extracellular loop, the seventh transmembrane domain, and 40 residues of the cytosolic tail (E3-M7-24-T40) was biosynthesized fused to a carrier protein. The multidomain fusion protein (designated M7FP) was purified to near homogeneity as judged by HPLC and characterized by mass spectrometry. In minimal medium, 30-40 mg of M7FP were obtained per liter of culture. The 73-residue peptide was released from its carrier by CNBr and obtained in wild-type, (15)N, and (13)C/(15)N forms. The E3-M7-24-T40 peptide integrated into 1-palmitoyl-2-hydroxy-sn-glycero-3-[phospho-rac-(1-glycerol)] and dodecylphosphocholine micelles at concentrations (200-500 microM) suitable for NMR investigations. HSQC experiments performed in organic solvents and detergent micelles on (15)N-labeled E3-M7-24-T40 showed a clear dispersion of the nitrogen-amide proton correlation cross-peaks indicative of a pure, uniformly labeled molecule that assumed a partially ordered structure. NOE connectivities, chemical shift indices, J-coupling analysis, and structural modeling suggested that in trifluoroethanol/water (1:1) helical subdomains existed in both the transmembrane and cytoslic tail of the multidomain peptide. Similar conclusions were reached in chloroform/methanol/water (4:4:1). As the cytosolic tail participates in down-regulation of Ste2p, the helical regions in the Ste2p tail may play a role in protein-protein interactions involved in endocytosis.

Heavier-than-air flying machines are impossible

Many G protein-coupled receptor (GPCR) models have been built over the years. The release of the structure of bovine rhodopsin in August 2000 enabled us to analyze models built before that period to learn more about the models we build today. We conclude that the GPCR modelling field is riddled with 'common knowledge' similar to Lord Kelvin's remark in 1895 that "heavier-than-air flying machines are impossible", and we summarize what we think are the (im)possibilities of modelling GPCRs using the coordinates of bovine rhodopsin as a template. Associated WWW pages: www.gpcr.org/articles/2003_mod

Biosynthesis and biophysical analysis of domains of a yeast G protein-coupled receptor

The alpha-factor receptor(Ste2p) is required for the sexual conjugation of the yeast Saccharomyces cerevisiae. Ste2p belongs to the G protein-coupled receptor (GPCR) family sharing a common heptahelical transmembrane structure. Biological synthesis of regions of Ste2p fused to a leader protein in Escherichia coli yielded milligram quantities of polypeptides that corresponded to one or two transmembrane domains. Fusion proteins were characterized by polyacrylamide gel electrophoresis, high performance liquid chromatography, and mass spectrometry. CD studies on the fusion proteins in trifluoroethanol:water mixtures indicated the existence of alpha-helical structures in the single- and the double-transmembrane domains. NMR experiments were performed in CDCl(3):CD(3)OH:H(2)O (4:4:1) on the (15)N-labeled single-transmembrane peptide showing a clear dispersion of the nitrogen-amide proton correlation cross peaks indicative of a high-purity, uniformly labeled molecule. The results indicate that single- and double-transmembrane domains of a GPCR can be produced by biosynthetic methods in quantities and purity sufficient for biophysical studies.

The three-dimensional structure of bovine rhodopsin determined by electron cryomicroscopy

G-protein-coupled receptors are integral membrane proteins that respond to environmental signals and initiate signal transduction pathways, which activate cellular processes. Rhodopsin, a well known member of the G-protein-coupled receptor family, is located in the disk membranes of the rod outer segment, where it is responsible for the visualization of dim light. Rhodopsin is the most extensively studied G-protein-coupled receptor, and knowledge about its structure serves as a template for other related receptors. We have gained detailed structural knowledge from the crystal structure (1), which was solved by x-ray crystallography in 2000 using three-dimensional crystals. Here we report a three-dimensional density map of bovine rhodopsin determined by electron cryomicroscopy of two-dimensional crystals with p22(1)2(1) symmetry. The usage of relatively small and disordered crystals made the process of structure determination challenging. Special attention was paid to the extraction of amplitudes and phases, since usable raw data were limited to a maximum tilt of 45 degrees. In the refinement process, an improved unbending procedure was applied. This led to a final resolution of 5.5 A in the membrane plane and approximately 13 A perpendicular to it, making our electron density map the most accurate map of a G-protein-coupled receptor currently available by electron microscopy. Most important is the information we gain about the center of the membrane plane and the orientation of the molecule relative to the bilayer. This information cannot be retrieved from the three-dimensional crystals. In our electron density map, all seven transmembrane helices were identified, and their arrangement is in agreement with the arrangement known from the crystal structure (1). In the retinal binding pocket, a density peak adjacent to helix 3 suggests the position of the beta-ionine ring of the chromophore, and in its vicinity several of the bigger amino acids can be identified.

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.

Recent progress in discovery of small-molecule CCR5 chemokine receptor ligands as HIV-1 inhibitors

This review addresses key pharmacology and virology issues relevant in discovery and development of CCR5 antagonists as anti-HIV drugs, such as target validation, receptor internalization, allosterism, viral resistance and tropism. Recent progress in the discovery and development of CCR5 antagonists, SAR and clinical status are reviewed. Finally, modeling-based structure of CCR5 is discussed in the context of a small-molecule antagonism of the CCR5 receptor.

The crystallographic model of rhodopsin and its use in studies of other G protein-coupled receptors

G protein-coupled receptors (GPCRs) are integral membrane proteins that respond to environmental signals and initiate signal transduction pathways activating cellular processes. Rhodopsin is a GPCR found in rod cells in retina where it functions as a photopigment. Its molecular structure is known from cryo-electron microscopic and X-ray crystallographic studies, and this has reshaped many structure/function questions important in vision science. In addition, this first GPCR structure has provided a structural template for studies of other GPCRs, including many known drug targets. After presenting an overview of the major structural elements of rhodopsin, recent literature covering the use of the rhodopsin structure in analyzing other GPCRs will be summarized. Use of the rhodopsin structural model to understand the structure and function of other GPCRs provides strong evidence validating the structural model.

Molecular modelling and site-directed mutagenesis of human GALR1 galanin receptor defines determinants of receptor subtype specificity

Human galanin is a 30 amino acid neuropeptide that elicits a range of biological activities by interaction with G protein-coupled receptors. We have generated a model of the human GALR1 galanin receptor subtype (hGALR1) based on the alpha carbon maps of frog rhodopsin and investigated the significance of potential contact residues suggested by the model using site-directed mutagenesis. Mutation of Phe186 within the second extracellular loop to Ala resulted in a 6-fold decrease in affinity for galanin, representing a change in free energy consistent with hydrophobic interaction. Our model suggests interaction between Phe186 of hGALR1 and Ala7 or Leu11 of galanin. Receptor subtype specificity was investigated by replacement of residues in hGALR1 with the corresponding residues in hGALR2 and use of the hGALR2-specific ligands hGalanin(2-30) and [D-Trp2]hGalanin(1-30). The His267Ile mutant receptor exhibited a pharmacological profile corresponding to that of hGALR1, suggesting that His267 is not involved in a receptor-ligand interaction. The mutation Phe115Ala resulted in a decreased binding affinity for hGalanin and for hGALR2-specific analogues, indicating Phe115 to be of structural importance to the ligand binding pocket of hGALR1 but not involved in direct ligand interaction. Analysis of Glu271Trp suggested that Glu271 of hGALR1 interacts with the N-terminus of galanin and that the Trp residue in the corresponding position in hGALR2 is involved in receptor subtype specificity of binding. Our model supports previous reports of Phe282 of hGALR1 interacting with Trp2 of galanin and His264 of hGALR1 interacting with Tyr9 of galanin.

Docking ligands to vasopressin and oxytocin receptors via genetic algorithm

The aim of the study was to computer-dock selected ligands to neurophyseal receptors in order to identify amino acid residues responsible for ligand-receptor interactions. To this aim, reliable oxytocin receptor (OTR) and arginine-vasopressin receptor (V1aR/V2R) models were built. The OTR-selective agonist [Thr4,Gly7]OT, the OTR-selective cyclohexapeptide antagonist L-366,948 and OT itself were docked via genetic algorithm to OTR, V1aR, and V2R and relaxed using a constrained simulated annealing protocol. For the analysis of receptor/ligand interactions a subset of initial conformations was chosen using energetic and steric criteria. All three ligands seem to prefer similar modes of binding to the receptors, manifested by repetitive residues of the receptors which directly interact with the ligands. Taking into account that many aspects of mechanisms of G protein-coupled receptor (GPCR) action are still unsolved, the results obtained with the docking simulations may propose future experimental research, especially in site-directed mutagenesis analysis and searching for key amino acid residues responsible for drug activities.

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.

Receptor activation: what does the rhodopsin structure tell us?

G-protein-coupled receptors (GPCRs) are a large family of seven-transmembrane-helix proteins that mediate responses to hormones, neurotransmitters and, in the case of rhodopsin, photons. The recent determination of the structure of rhodopsin at atomic resolution opens avenues to a deeper understanding of GPCR activation and transmembrane signaling. Data from previous crosslinking, spin labeling and scanning accessibility experiments on rhodopsin have been mapped onto the high-resolution structure. These data correlate well and are consistent with the structure, and suggest that activation by light opens a cleft at the cytoplasmic end of the seven-helix bundle of rhodopsin. Furthermore, lessons learned from rhodopsin might also apply to other members of this essential family of receptors. (For an animation of the crystal structure of rhodopsin see http://archive.bmn.com/supp/tips/tips2211a.html)

Confronting complexity: the interlink of phototransduction and retinoid metabolism in the vertebrate retina

Absorption of light by rhodopsin or cone pigments in photoreceptors triggers photoisomerization of their universal chromophore, 11-cis-retinal, to all-trans-retinal. This photoreaction is the initial step in phototransduction that ultimately leads to the sensation of vision. Currently, a great deal of effort is directed toward elucidating mechanisms that return photoreceptors to the dark-adapted state, and processes that restore rhodopsin and counterbalance the bleaching of rhodopsin. Most notably, enzymatic isomerization of all-trans-retinal to 11-cis-retinal, called the visual cycle (or more properly the retinoid cycle), is required for regeneration of these visual pigments. Regeneration begins in rods and cones when all-trans-retinal is reduced to all-trans-retinol. The process continues in adjacent retinal pigment epithelial cells (RPE), where a complex set of reactions converts all-trans-retinol to 11-cis-retinal. Although remarkable progress has been made over the past decade in understanding the phototransduction cascade, our understanding of the retinoid cycle remains rudimentary. The aim of this review is to summarize recent developments in our current understanding of the retinoid cycle at the molecular level, and to examine the relevance of these reactions to phototransduction.

An activation switch in the ligand binding pocket of the C5a receptor

Although agonists are thought to occupy binding pockets within the seven-helix core of serpentine receptors, the topography of these binding pockets and the conformational changes responsible for receptor activation are poorly understood. To identify the ligand binding pocket in the receptor for complement factor 5a (C5aR), we assessed binding affinities of hexapeptide ligands, each mutated at a single position, for seven mutant C5aRs, each mutated at a single position in the putative ligand binding site. In ChaW (an antagonist) and W5Cha (an agonist), the side chains at position 5 are tryptophan and cyclohexylalanine, respectively. Comparisons of binding affinities indicated that the hexapeptide residue at this position interacts with two C5aR residues, Ile-116 (helix III) and Val-286 (helix VII); in a C5aR model these two side chains point toward one another. Both the I116A and the V286A mutations markedly increased binding affinity of W5Cha but not that of ChaW. Moreover, ChaW, the antagonist hexapeptide, acted as a full agonist on the I116A mutant. These results argue that C5aR residues Ile-116 and Val-286 interact with the side chain at position 5 of the hexapeptide ligand to form an activation switch. Based on this and previous work, we present a docking model for the hexapeptide within the C5aR binding pocket. We propose that agonists induce a small change in the relative orientations of helices III and VII and that these helices work together to allow movement of helix VI away from the receptor core, thereby triggering G protein activation.

Functionally discrete mimics of light-activated rhodopsin identified through expression of soluble cytoplasmic domains

Numerous studies on the seven-helix receptor rhodopsin have implicated the cytoplasmic loops and carboxyl-terminal region in the binding and activation of proteins involved in visual transduction and desensitization. In our continuing studies on rhodopsin folding, assembly, and structure, we have attempted to reconstruct the interacting surface(s) for these proteins by inserting fragments corresponding to the cytoplasmic loops and/or the carboxyl-terminal tail of bovine opsin either singly, or in combination, onto a surface loop in thioredoxin. The purpose of the thioredoxin fusion is to provide a soluble scaffold for the cytoplasmic fragments thereby allowing them sufficient conformational freedom to fold to a structure that mimics the protein-binding sites on light-activated rhodopsin. All of the fusion proteins are expressed to relatively high levels in Escherichia coli and can be purified using a two- or three-step chromatography procedure. Biochemical studies show that some of the fusion proteins effectively mimic the activated conformation(s) of rhodopsin in stimulating G-protein or competing with the light-activated rhodopsin/G-protein interaction, in supporting phosphorylation of the carboxyl-terminal opsin fragment by rhodopsin kinase, and/or phosphopeptide-stimulated arrestin binding. These results suggest that specific segments of the cytoplasmic surface of rhodopsin can adopt functionally discrete conformations in the absence of the connecting transmembrane helices and retinal chromophore.

Modeling and docking the endothelin G-protein-coupled receptor

A model of the endothelin G-protein-coupled receptor (ET(A)) has been constructed using a segmented approach. The model was produced using a bovine rhodopsin model as a template for the seven transmembrane alpha-helices. The three cytoplasmic loop regions and the C-terminal region were modeled on NMR structures of corresponding segments from bovine rhodopsin. The three extracellular loops were modeled on homologous loop regions in other proteins of known structure. The N-terminal region was modeled as a three-helix domain based on its homology with a hydrolase protein. To test the model, the FTDOCK algorithm was used to predict the ligand-binding site for the crystal structure of human endothelin. The site of docking is consistent with mutational and biochemical data. The principal sites of interaction in the endothelin ligand all lie on one face of a helix that has been implicated by structure-activity relationship studies as being essential for binding. As further support for the model, attempts to dock bigET, an inactive precursor to endothelin that does not bind to the receptor, found no sites for tight binding. The model of the receptor-ligand complex produced forms a basis for rational drug design of agonists and antagonists for this G-protein-coupled receptor.

A triangle lattice model that predicts transmembrane helix configuration using a polar jigsaw puzzle

We developed a method of predicting the tertiary structures of seven transmembrane helical proteins in triangle lattice models, assuming that the configuration of helices is stabilized by polar interactions. Triangle lattice models having 12 or 11 nearest neighbor pairs were used as general templates of a seven-helix system, then the orientation angles of all helices were varied at intervals of 15 degrees. The polar interaction energy for all possible positions of each helix was estimated using the calculated polar indices of transmembrane helices. An automated system was constructed and applied to bacteriorhodopsin, a typical membrane protein with seven transmembrane helices. The predicted optimal and actual structures were similar. The top 100 predicted helical configurations indicated that the helix-triangle, CFG, occurred at the highest frequency. In fact, this helix-triangle of bacteriorhodopsin forms an active proton-pumping site, suggesting that the present method can identify functionally important helices in membrane proteins. The possibility of studying the structure change of bacteriorhodopsin during the functional process by this method is discussed, and may serve to explain the experimental structures of photointermediate states.

Structure. Rhodopsin sees the light

Members of the seven transmembrane receptor superfamily bind a remarkable variety of ligands, from neurotransmitters to odorants, and activate a spectacular array of G protein signaling molecules. These G-protein coupled receptors (GPCRs) are important in many cellular functions and so there has been great interest in elucidating how they transmit their signals to the interior of the cell after activation by ligand. As Bourne and Meng explain in their Perspective, the molecular movements of activated GPCRs are becoming clear now that the first crystal structure of a GPCR (rhodopsin, the light-trapping receptor found in the retina of the eye) has been reported (Palczweski et al.).

Crystal structure of rhodopsin: A G protein-coupled receptor

Heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) respond to a variety of different external stimuli and activate G proteins. GPCRs share many structural features, including a bundle of seven transmembrane alpha helices connected by six loops of varying lengths. We determined the structure of rhodopsin from diffraction data extending to 2.8 angstroms resolution. The highly organized structure in the extracellular region, including a conserved disulfide bridge, forms a basis for the arrangement of the seven-helix transmembrane motif. The ground-state chromophore, 11-cis-retinal, holds the transmembrane region of the protein in the inactive conformation. Interactions of the chromophore with a cluster of key residues determine the wavelength of the maximum absorption. Changes in these interactions among rhodopsins facilitate color discrimination. Identification of a set of residues that mediate interactions between the transmembrane helices and the cytoplasmic surface, where G-protein activation occurs, also suggests a possible structural change upon photoactivation.

The molecular origin of the inhibition of transducin activation in rhodopsin lacking the 9-methyl group of the retinal chromophore: a UV-Vis and FTIR spectroscopic study

The formation of the active rhodopsin state metarhodopsin II (MII) is believed to be partially governed by specific steric constraints imposed onto the protein by the 9-methyl group of the retinal chromophore. We studied the properties of the synthetic pigment 9-demethyl rhodopsin (9dm-Rho), consisting of the rhodopsin apoprotein regenerated with synthetic retinal lacking the 9-methyl group, by UV-vis and Fourier transform infrared difference spectroscopy. Low activation rates of the visual G-protein transducin by the modified pigment reported in previous studies are shown to not be caused by the reduced activity of its MII state, but to be due to a dramatic equilibrium shift from MII to its immediate precursor, MI. The MII state of 9dm-Rho displays only a partial deprotonation of the retinal Schiff base, leading to the formation of two MII subspecies absorbing at 380 and 470 nm, both of which seem to be involved in transducin activation. The rate of MII formation is slowed by 2 orders of magnitude compared to rhodopsin. The dark state and the MI state of 9dm-Rho are distinctly different from their respective states in the native pigment, pointing to a more relaxed fit of the retinal chromophore in its binding pocket. The shifted equilibrium between MI and MII is therefore discussed in terms of an increased entropy of the 9dm-Rho MI state due to changed steric interactions.

The nociceptin (ORL1) receptor: molecular cloning and functional architecture

Nociceptin and the ORL1 receptor share high sequence similarity with opioid peptides, particularly dynorphin A, and their receptors. However, nociceptin and dynorphin A may use distinct molecular pathways to bind and activate their cognate receptors. Activation of the kappa-opioid receptor by dynorphin A is thought to require interactions of its N-terminal hydrophobic domain (Y(1)GGF) with the receptor opioid binding pocket, located within the transmembrane helix bundle, while activation of the ORL1 receptor appears to require interactions of the positively charged core (R(8)KSARK) of nociceptin with the negatively charged second extracellular receptor loop.

On the spatial disposition of the fifth transmembrane helix and the structural integrity of the transmembrane binding site in the opioid and ORL1 G protein-coupled receptor family

Evidence from statistical cluster analyses of a multiple sequence alignment of G protein-coupled receptor seven-helix folds supports the existence of structurally conserved transmembrane (TM) ligand binding sites in the opioid/opioid receptor-like (ORL1) and amine receptor families. Based on the expectation that functionally conserved regions in homologous proteins will display locally higher levels of sequence identity compared with global sequence similarities that pertain to the overall fold, this approach may have wider applications in functional genomics to annotate sequence data. Binding sites in models of the kappa-opioid receptor seven-helix bundle built from the rhodopsin templates of Baldwin et al. (1997) [J. Mol. Biol., 272, 144-164] and Herzyk and Hubbard (1998) [J. Mol. Biol., 281, 742-751] are compared. The Herzyk and Hubbard template is found to be in better accord with experimental studies of amine, opioid and rhodopsin receptors owing to the reduced physical separation of the extracellular parts of TM helices V and VI and differences in the rotational orientation of the N-terminal of helix V that reveal side chain accessibilities in the Baldwin et al. structure to be out of phase with relative alkylation rates of engineered cysteine residues in the TM binding site of the alpha(2A)-adrenergic receptor. TM helix V in the Baldwin et al. template has been remodelled with a different proline kink to satisfy experimental constraints. A recent proposal that rotation of helix V is associated with receptor activation is critically discussed.

Movement of retinal along the visual transduction path

Movement of the ligand/receptor complex in rhodopsin (Rh) has been traced. Bleaching of diazoketo rhodopsin (DK-Rh) containing 11-cis-3-diazo-4-oxo-retinal yields batho-, lumi-, meta-I-, and meta-II-Rh intermediates corresponding to those of native Rh but at lower temperatures. Photoaffinity labeling of DK-Rh and these bleaching intermediates shows that the ionone ring cross-links to tryptophan-265 on helix F in DK-Rh and batho-Rh, and to alanine-169 on helix D in lumi-, meta-I-, and meta-II-Rh intermediates. It is likely that these movements involving a flip-over of the chromophoric ring trigger changes in cytoplasmic membrane loops resulting in heterotrimeric guanine nucleotide-binding protein (G protein) activation.

Observations of light-induced structural changes of retinal within rhodopsin

Photo-isomerization of the 11-cis retinal chromophore activates the mammalian light-receptor rhodopsin, a representative member of a major superfamily of transmembrane G-protein-coupled receptor proteins (GPCRs) responsible for many cell signal communication pathways. Although low-resolution (5 A) electron microscopy studies confirm a seven transmembrane helix bundle as a principal structural component of rhodopsin, the structure of the retinal within this helical bundle is not known in detail. Such information is essential for any theoretical or functional understanding of one of the fastest occurring photoactivation processes in nature, as well as the general mechanism behind GPCR activation. Here we determine the three-dimensional structure of 11-cis retinal bound to bovine rhodopsin in the ground state at atomic level using a new high-resolution solid-state NMR method. Significant structural changes are observed in the retinal following activation by light to the photo-activated M(I) state of rhodopsin giving the all-trans isomer of the chromophore. These changes are linked directly to the activation of the receptor, providing an insight into the activation mechanism of this class of receptors at a molecular level.

X-Ray diffraction analysis of three-dimensional crystals of bovine rhodopsin obtained from mixed micelles

Rhodopsin, a prototypic G protein-coupled receptor responsible for absorption of photons in retinal rod photoreceptor cells, was selectively extracted from bovine rod outer segment membranes, employing mixed micelles of nonyl beta-d-glucoside and heptanetriol. Highly purified rhodopsin was crystallized from solutions containing varying amounts of detergent and amphiphile. The crystals contained ground state rhodopsin molecules as judged by their red color and the linear dichroism originating from the 11-cis-retinal chromophore. However, when exposed to visible light, even at 4 degrees C, rhodopsin was bleached and the crystals decomposed. Reflections in the diffraction pattern were observed out to 3.5-A resolution at 100 K for the most ordered crystals. Diffraction data have been processed to 3.85-A resolution. The symmetry of the diffraction pattern and the systematic absences indicate that the crystals have tetragonal symmetry, space group P4(1)22 or P4(3)22, a = b = 96.51 A, c = 148.55 A. A value of 4.12 A(3)/Da for V(M) was obtained for one monomer in the asymmetric unit (eight molecules per unit cell). Our study is the first characterization of a three-dimensional crystal of a G protein-coupled receptor and may be valuable for future structural studies on related receptors of this important superfamily.

Modelling G-protein-coupled receptors for drug design

The G-protein coupled receptors form a large and diverse multi-gene superfamily with many important physiological functions. As such, they have become important targets in pharmaceutical research. Molecular modelling and site-directed mutagenesis have played an important role in our increasing understanding of the structural basis of drug action at these receptors. Aspects of this understanding, how these techniques can be used within a drug-design programme, and remaining challenges for the future are reviewed.

Conversion of agonist site to metal-ion chelator site in the beta(2)-adrenergic receptor

Previously metal-ion sites have been used as structural and functional probes in seven transmembrane receptors (7TM), but as yet all the engineered sites have been inactivating. Based on presumed agonist interaction points in transmembrane III (TM-III) and -VII of the beta(2)-adrenergic receptor, in this paper we construct an activating metal-ion site between the amine-binding Asp-113 in TM-III-or a His residue introduced at this position-and a Cys residue substituted for Asn-312 in TM-VII. No increase in constitutive activity was observed in the mutant receptors. Signal transduction was activated in the mutant receptors not by normal catecholamine ligands but instead either by free zinc ions or by zinc or copper ions in complex with small hydrophobic metal-ion chelators. Chelation of the metal ions by small hydrophobic chelators such as phenanthroline or bipyridine protected the cells from the toxic effect of, for example Cu(2+), and in several cases increased the affinity of the ions for the agonistic site. Wash-out experiments and structure-activity analysis indicated, that the high-affinity chelators and the metal ions bind and activate the mutant receptor as metal ion guided ligand complexes. Because of the well-understood binding geometry of the small metal ions, an important distance constraint has here been imposed between TM-III and -VII in the active, signaling conformation of 7TM receptors. It is suggested that atoxic metal-ion chelator complexes could possibly in the future be used as generic, pharmacologic tools to switch 7TM receptors with engineered metal-ion sites on or off at will.

Is the protein/lipid hydrophobic matching principle relevant to membrane organization and functions?

Biological membranes are complex and well-organized multimolecular assemblies composed of a wide variety of protein and lipid molecular species. If such a diversity in protein and lipid polar headgroup structures may easily be related to a large panel of functions, the wide dispersion in acyl chain length and structure which the lipids display is more difficult to understand. It is not required for maintaining bilayer assembly and fluidity. Direct information on the lateral distribution of these various molecular species, on their potential specificity for interaction between themselves and with proteins and on the functional implications of these interactions is also still lacking. Because hydrophobic interactions play a major role in stabilizing membrane structures, we suggest considering the problem from the point of view of the matching of the hydrophobic surface of proteins by the acyl chains of the lipids. After a brief introduction to the hydrophobic matching principle, we will present experimental results which demonstrate the predictive power of the current theories and then, we will introduce the new and important concept of protein/lipid sorting in membranes. Finally, we will show how the hydrophobic matching condition may play a key role in the membrane organization and function.

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.

Kinetics of the light-induced proton translocation associated with the pH-dependent formation of the metarhodopsin I/II equilibrium of bovine rhodopsin

The kinetics of the formation of the metaII (MII) state of bovine rhodopsin was investigated by time-resolved electrical and absorption measurements with rod outer segment (ROS) fragments. Photoexcitation leads to proton transfer in the direction from the cytosolic to the intradiscal side of the membrane, probably from the Schiff base to the acceptor glutamate 113. Two components of comparable amplitude are required to describe the charge movement with exponential times of 1.1 (45%) and 3.0 ms (55%) (pH 7.8, 22 degreesC, 150 mM KCl). The corresponding activation energies are 86 and 123 kJ/mol, respectively (150 mM KCl). The time constants and amplitudes depend strongly on pH. Between pH 7.1 and 3.8 the kinetics becomes much faster, with the faster and slower components accelerating by factors of about 8 and 2, respectively. Complementary single-flash absorption experiments at 380 nm and 10 degreesC show that the formation of MII also occurs with two components with similar time constants and pH dependence. This suggests that both signals monitor the same molecular events. The pH dependence of the two apparent time constants and amplitudes of the optical data can be described well over the pH range 4-7.5 by two coupled equilibria between MI and two isochromic MII species MIIa and MIIb: MI MIIa(380) MIIb(380), with k0 proportional to the proton concentration. This model implies that deprotonation of the Schiff base and proton uptake are tightly coupled in ROS membranes. Models with k2 proportional to the proton concentration cannot describe the data. Photoreversal of MII by blue flashes (420 nm) leads to proton transfer in a direction opposite to that of the signal associated with MII formation. In this transition the Schiff base is reprotonated, most likely from glutamate 113. At pH 7.3, 150 mM KCl, 22 degreesC, this electrical charge reversal has an exponential time constant of about 30 ms and is about 10 times slower than the forward charge motion.