The retinal conformation and its environment in rhodopsin in light of a new 2.2 A crystal structure

Computer aided comparative analysis of the binding modes of the adenosine receptor agonists for all known subtypes of adenosine receptors

Molecular models of all known subtypes (A1, A2A, A2B, and A3) of the human adenosine receptors were built in homology with bovine rhodopsin. These models include the transmembrane domain as well as all extracellular and intracellular hydrophilic loops and terminal domains. The molecular docking of adenosine and 46 selected derivatives was performed for each receptor subtype. A binding mode common for all studied agonists was proposed, and possible explanations for differences in the ligand activities were suggested.

The Molecular Basis for Ligand Specificity in a Mouse Olfactory Receptor

Sequence differences between members of the mouse olfac-tory receptor MOR42 subfamily (MOR42-3 and MOR42-1) are likely to be the basis for variation in ligand binding preference among these receptors. We investigated the specificity of MOR42-3 for a variety of dicarboxylic acids. We used site-directed mutagenesis, guided by homology modeling and ligand docking studies, to locate functionally important residues. Receptors were expressed in Xenopus oocytes and assayed using high throughput electrophysiology. The importance of the Val-113 residue, located deep within the receptor, was analyzed in the context of interhelical interactions. We also screened additional residues predicted to be involved in ligand binding site, based on comparison of ortholog/paralog pairs from the mouse and human olfactory receptor genomes (Man, O., Gilad, Y., and Lancet, D. (2004) Protein Sci. 13, 240-254). A network of 8 residues in transmembrane domains III, V, and VI was identified. These residues form part of the ligand binding pocket of MOR42-3. C12 dicarboxylic acid did not activate the receptor in our functional assay, yet our docking simulations predicted its binding site in MOR42-3. Binding without activation implied that C12 dicarboxylic acid might act as an antagonist. In our functional assay, C12 dicarboxylic acid did indeed act as an antagonist of MOR42-3, in agreement with molecular docking studies. Our results demonstrate a powerful approach based on the synergy between computational predictions and physiological assays.

Effect of different treatments of long-range interactions and sampling conditions in molecular dynamic simulations of rhodopsin embedded in a dipalmitoyl phosphatidylcholine bilayer

The present study analyzes the effect of the simulation conditions on the results of molecular dynamics simulations of G-protein coupled receptors (GPCRs) performed with an explicit lipid bilayer. Accordingly, the present work reports the analysis of different simulations of bovine rhodopsin embedded in a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer using two different sampling conditions and two different approaches for the treatment of long-range electrostatic interactions. Specifically, sampling was carried out either by using the statistical ensembles NVT or NPT (constant number of atoms, a pressure of 1 atm in all directions and fixed temperature), and the electrostatic interactions were treated either by using a twin-cutoff, or the particle mesh Ewald summation method (PME). The results of the present study suggest that the use of the NPT ensemble in combination with the PME method provide more realistic simulations. The use of NPT during the equilibration avoids the need of an a priori estimation of the box dimensions, giving the correct area per lipid. However, once the system is equilibrated, the simulations are irrespective of the sampling conditions used. The use of an electrostatic cutoff induces artifacts on both lipid thickness and the ion distribution, but has no direct effect on the protein and water molecules.

Metals in membranes

In this critical review we discuss recent advances in understanding the modes of interaction of metal ions with membrane proteins, including channels, pumps, transporters, ATP-binding cassette proteins, G-protein coupled receptors, kinases and respiratory enzymes. Such knowledge provides a basis for elucidating the mechanism of action of some classes of metallodrugs, and a stimulus for the further exploration of the coordination chemistry of metal ions in membranes. Such research offers promise for the discovery of new drugs with unusual modes of action. The article will be of interest to bioinorganic chemists, chemical biologists, biochemists, pharmacologists and medicinal chemists. (247 references).

Human prostacyclin receptor structure and function from naturally-occurring and synthetic mutations

Prostacyclin (PGI2) is released by vascular endothelial cells and serves as a potent vasodilator, inhibitor of platelet aggregation (anti-thrombotic), and moderator of vascular smooth muscle cell proliferation-migration-differentiation (anti-atherosclerotic). These actions are mediated via a seven transmembrane-spanning G-protein coupled receptor (GPCR), known as the human prostacyclin receptor or hIP. Animal studies using prostacyclin receptor knock-out (IP-/-) mice have revealed increased propensities towards thrombosis, intimal hyperplasia, atherosclerosis, restenosis, as well as reperfusion injury. Of further importance has been the world-wide withdrawal of selective COX-2 inhibitors, due to their discriminating suppression of COX-2-derived PGI2 and its cardioprotective effects, leading to increased cardiovascular events, including myocardial infarction and thrombotic stroke. Over the last decade, mutagenesis studies of the IP receptor, in conjunction with in vitro functional assays and molecular modeling, have provided critical insights into the molecular mechanisms of both agonist binding and receptor activation. Most recently, the discovery of naturally-occurring and dysfunctional mutations within the hIP has provided additional insights into the proposed cardioprotective role of prostacyclin. The aim of this review is to summarize the most recent findings regarding hIP receptor structure-function that have developed through the study of both synthetic and naturally-occurring mutations.

The receptor-bound "empty pocket" state of the heterotrimeric G-protein alpha-subunit is conformationally dynamic

Heterotrimeric G-protein activation by a G-protein-coupled receptor (GPCR) requires the propagation of structural signals from the receptor-interacting surfaces to the guanine nucleotide-binding pocket. To probe conformational changes in the G-protein alpha-subunit (G(alpha)) associated with activated GPCR (R*) interactions and guanine nucleotide exchange, high-resolution solution NMR methods are being applied in studying signaling of the G-protein, transducin, by light-activated rhodopsin. Using these methods, we recently demonstrated that an isotope-labeled G(alpha) reconstituted heterotrimer forms functional complexes under NMR experimental conditions with light-activated, detergent-solubilized rhodopsin and a soluble mimic of R*, both of which trigger guanine nucleotide exchange [Ridge, K. D., et al. (2006) J. Biol. Chem. 281, 7635-7648]. Here, it is shown that both light-activated rhodopsin and the soluble mimic of R form trapped intermediate complexes with a GDP-released "empty pocket" state of the heterotrimer in the absence of GTP (or GTPgammaS). In contrast to guanine nucleotide-bound forms of G(alpha), the NMR spectra of the GDP-released, R-bound empty pocket state of G(alpha) display severe line broadening suggestive of a dynamic intermediate state. Interestingly, the conformation of a GDP-depleted, Mg(2+)-bound state of G(alpha) generated in a manner independent of R* does not exhibit a similar degree of line broadening but rather appears structurally similar to the GDP/Mg(2+)-bound form of the protein. Taken together, these results suggest that R*-mediated changes in the receptor-interacting regions of G(alpha), and not the absence of bound guanine nucleotide, are the predominant factors which dictate G(alpha) conformation and dynamics in this R*-bound state of the heterotrimer.

Activation of the CXCR3 Chemokine Receptor through Anchoring of a Small Molecule Chelator Ligand between TM-III, -IV, and -VI

Seven transmembrane segment (7TM) receptors are activated through a common, still rather unclear molecular mechanism by a variety of chemical messengers ranging from monoamines to large proteins. By introducing a His residue at position III:05 in the CXCR3 receptor a metal ion site was built between the extracellular ends of transmembrane (TM) III and TM-IV to anchor aromatic chelators at a location corresponding to the presumed binding pocket for adrenergic receptor agonists. In this construct, free metal ions had no agonistic effect in accordance with the optimal geometry of the metal ion site in molecular models built over the inactive form of rhodopsin. In contrast, the aromatic chelators bipyridine or phenanthrolene in complex with Zn(II) or Cu(II) acted as potent agonists displaying signaling efficacies similar to or even better than the endogenous chemokine agonists. Molecular modeling and molecular simulations combined with mutational analysis indicated that the metal ion site-anchored chelators act as agonists by establishing an aromatic-aromatic, second-site interaction with TyrVI:16 on the inner face of TM-VI. It is noteworthy that this interaction required that the extracellular segment of TM-VI moves inward in the direction of TM-III, whereby TyrVI:16 together with the chelators complete an "aromatic zipper" also comprising PheIII:08 (corresponding to the monoamine receptor anchoring point) and TyrVII:10 (corresponding to the retinal attachment site in rhodopsin). Chemokine agonism was independent of this aromatic zipper. It is proposed that in rhodopsin-like 7TM receptors, small-molecule compounds in general act as agonists in a similar manner as here demonstrated with the artificial, metal ion site anchored chelators, by holding TM-VI bent inward.

Identification of Leu276 of the S1P1 receptor and Phe263 of the S1P3 receptor in interaction with receptor specific agonists by molecular modeling, site-directed mutagenesis, and affinity studies

Sphingosine-1-phosphate (S1P) receptor agonists are novel immunosuppressive agents. The selectivity of S1P1 against S1P3 is strongly correlated with lymphocyte sequestration and minimum acute toxicity and bradycardia. This study describes molecular modeling, site-directed mutagenesis, and affinity studies exploring the molecular basis for selectivity between S1P1 and S1P3 receptors. Computational models of human S1P1 and S1P3 receptors bound with two nonselective agonists or two S1P1-selective agonists were developed based on the X-ray crystal structure of bovine rhodopsin. The models predict that S1P1 Leu276 and S1P3 Phe263 contribute to the S1P1/S1P3 selectivity of the two S1P1-selective agonists. These residues were subjected to site-directed mutagenesis. The wild-type and mutant S1P receptors were expressed in Chinese hamster ovary cells and examined for their abilities to bind to and be activated by agonists in vitro. The results indicate that the mutations have minimal effects on the activities of the two nonselective agonists, although they have dramatic effects on the S1P1-selective agonists. These studies provide a fundamental understanding of how these two receptor-selective agonists bind to the S1P1 and S1P3 receptors, which should aid development of more selective S1P1 receptor agonists with immunosuppressive properties and improved safety profiles.

The Role of the β-Ionone Ring in the Photochemical Reaction of Rhodopsin

Time-dependent density functional theory (TDDFT) calculations on the photoabsorption process of the 11-cis retinal protonated Schiff base (PSB) chromophore show that the Franck-Condon relaxation of the first excited state of the chromophore involves a torsional twist motion of the beta-ionone ring relative to the conjugated retinyl chain. For the ground state, the beta-ionone ring and the retinyl chain of the free retinal PSB chromophore form a -40 degrees dihedral angle as compared to -94 degrees for the first excited state. The double bonds of the retinal are shorter for the fully optimized structure of the excited state than for the ground state suggesting a higher cis-trans isomerization barrier for the excited state than for the ground state. According to the present TDDFT calculations, the excitation of the retinal PSB chromophore does not primarily lead to a reaction along the cis-trans torsional coordinate at the C11-C12 bond. The activation of the isomerization center seems to occur at a later stage of the photo reaction. The results obtained at the TDDFT level are supported by second-order Møller-Plesset (MP2) and approximate singles and doubles-coupled cluster (CC2) calculations on retinal chromophore models; the MP2 and CC2 calculations yield for them qualitatively the same ground state and excited-state structures as obtained in the density functional theory and TDDFT calculations.

Retinal Counterion Switch Mechanism in Vision Evaluated by Molecular Simulations

Photoisomerization of the retinylidene chromophore of rhodopsin is the starting point in the vision cascade. A counterion switch mechanism that stabilizes the retinal protonated Schiff base (PSB) has been proposed to be an essential step in rhodopsin activation. On the basis of vibrational and UV-visible spectroscopy, two counterion switch models have emerged. In the first model, the PSB is stabilized by Glu181 in the meta I state, while in the most recent proposal, it is stabilized by Glu113 as well as Glu181. We assess these models by conducting a pair of microsecond scale, all-atom molecular dynamics simulations of rhodopsin embedded in a 99-lipid bilayer of SDPC, SDPE, and cholesterol (2:2:1 ratio) varying the starting protonation state of Glu181. Theoretical simulations gave different orientations of retinal for the two counterion switch mechanisms, which were used to simulate experimental 2H NMR spectra for the C5, C9, and C13 methyl groups. Comparison of the simulated 2H NMR spectra with experimental data supports the complex-counterion mechanism. Hence, our results indicate that Glu113 and Glu181 stabilize the retinal PSB in the meta I state prior to activation of rhodopsin.

A methyl group at C7 of 11-cis-retinal allows chromophore formation but affects rhodopsin activation

The newly synthesized 11-cis-7-methylretinal can form an artificial visual pigment with kinetic and spectroscopic properties similar to the native pigment in the dark-state. However, its photobleaching behavior is altered, showing a Meta I-like photoproduct. This behavior reflects a steric constraint imposed by the 7-methyl group that affects the conformational change in the binding pocket as a result of retinal photoisomerization. Transducin activation is reduced, when compared to the native pigment with 11-cis-retinal. Molecular dynamics simulations suggest coupling of the C7 methyl group and the beta-ionone ring with Met207 in transmembrane helix 5 in agreement with recent experimental results.

NMR analysis of rhodopsin–transducin interactions

Heterotrimeric G-protein activation by an agonist-stimulated G-protein coupled receptor (R*) requires the propagation of structural signals from the receptor interacting surfaces to the guanine nucleotide-binding pocket. Employing high-resolution NMR methods, we are probing heterotrimer-associated and rhodopsin-stimulated changes in an isotope-labeled G-protein alpha-subunit (G(alpha)). A key aspect of the work involves the trapping and interrogation of discrete R*-bound conformations of G(alpha). Our results demonstrate that functionally important changes in G(alpha) structure and dynamics can be detected and characterized by NMR, enabling the generation of robust models for the global and local structural changes accompanying signal transfer from R* to the G-protein.

Improvements in G protein-coupled receptor purification yield light stable rhodopsin crystals

G protein-coupled receptors (GPCRs) represent the largest family of transmembrane signaling proteins and are the target of approximately half of all therapeutic agents. Agonist ligands bind their cognate GPCRs stabilizing the active conformation that is competent to bind G proteins, thus initiating a cascade of intracellular signaling events leading to modification of the cell activity. Despite their biomedical importance, the only known GPCR crystal structures are those of inactive rhodopsin forms. In order to understand how GPCRs are able to transduce extracellular signals across the plasma membrane, it is critical to determine the structure of these receptors in their ligand-bound, active state. Here, we report a novel combination of purification procedures that allowed the crystallization of rhodopsin in two new crystal forms and can be applicable to the purification and crystallization of other membrane proteins. Importantly, these new crystals are stable upon photoactivation and the preliminary X-ray diffraction analysis of both photoactivated and ground state rhodopsin crystals are also reported.

A comprehensive structure-function map of the intracellular surface of the human C5a receptor. I. Identification of critical residues

G protein-coupled receptors are one of the largest protein families in nature; however, the mechanisms by which they activate G proteins are still poorly understood. To identify residues on the intracellular face of the human C5a receptor that are involved in G protein activation, we performed a genetic analysis of each of the three intracellular loops and the carboxyl-terminal tail of the receptor. Amino acid substitutions were randomly incorporated into each loop, and functional receptors were identified in yeast. The third intracellular loop contains the largest number of preserved residues (positions resistant to amino acid substitutions), followed by the second loop, the first loop, and lastly the carboxyl terminus. Surprisingly, complete removal of the carboxyl-terminal tail did not impair C5a receptor signaling. When mapped onto a three-dimensional structural model of the inactive state of the C5a receptor, the preserved residues reside on one half of the intracellular surface of the receptor, creating a potential activation face. Together these data provide one of the most comprehensive functional maps of the intracellular surface of any G protein-coupled receptor to date.

Ligand specificity of odorant receptors

Odorant receptors belong to class A of the G protein-coupled receptors (GPCRs) and detect a large number of structurally diverse odorant molecules. A recent structural bioinformatic analysis suggests that structural features are conserved across class A of GPCRs in spite of their low sequence identity. Based on this work, we have aligned the sequences of 29 ORs for which ligand binding data are available. Recent site-directed mutagenesis experiments on one such receptor (MOR174-9) provide information that helped to identify nine amino-acid residues involved in ligand binding. Our modeling provides a rationale for amino acids in equivalent positions in most of the odorant receptors considered and helps to identify other amino acids that could be important for ligand binding. Our findings are consistent with most of the previous models and allow predictions for site-directed mutagenesis experiments, which could also validate our model.

Designer short peptide surfactants stabilize G protein-coupled receptor bovine rhodopsin

Membrane proteins play vital roles in every aspect of cellular activities. To study diverse membrane proteins, it is crucial to select the right surfactants to stabilize them for analysis. Despite much effort, little progress has been made in elucidating their structure and function, largely because of a lack of suitable surfactants. Here we report the stabilization of a G protein-coupled receptor bovine rhodopsin in solution, using a new class of designer short and simple peptide surfactants. These surfactants consist of seven amino acids with a hydrophilic head, aspartic acid or lysine, and a hydrophobic tail with six consecutive alanines. These peptide surfactants not only enhance the stability of bovine rhodopsin in the presence of lipids and the common surfactants n-dodecyl-beta-D-maltoside and octyl-D-glucoside, but they also significantly stabilize rhodopsin under thermal denaturation conditions, even after lipids are removed. These peptide surfactants are simple, versatile, effective, and affordable. They represent a designer molecular nanomaterial for use in studies of diverse elusive membrane proteins.

G protein coupled receptor structure and activation

G protein coupled receptors (GPCRs) are remarkably versatile signaling molecules. The members of this large family of membrane proteins are activated by a spectrum of structurally diverse ligands, and have been shown to modulate the activity of different signaling pathways in a ligand specific manner. In this manuscript I will review what is known about the structure and mechanism of activation of GPCRs focusing primarily on two model systems, rhodopsin and the beta(2) adrenoceptor.

The color of rhodopsins at the ab initio multiconfigurational perturbation theory resolution

We demonstrate that "brute force" quantum-mechanics/molecular-mechanics computations based on ab initio (i.e., first principles) multiconfigurational perturbation theory can reproduce the absorption maxima of a set of modified bovine rhodopsins with an accuracy allowing for the analysis of the factors determining their colors. In particular, we show that the theory accounts for the changes in excitation energy even when the proteins display the same charge distribution. Three color-tuning mechanisms, leading to changes of close magnitude, are demonstrated to operate in these conditions. The first is based on the change of the conformation of the conjugated backbone of the retinal chromophore. The second operates through the control of the distance between the positive charge residing on the chromophore and the carboxylate counterion. Finally, the third mechanism operates through the changes in orientation of the chromophore relative to the protein. These results offer perspectives for the unbiased computational design of mutants or chemically modified proteins with wanted optical properties.

Helix 8 of the Viral Chemokine Receptor ORF74 Directs Chemokine Binding

The constitutively active G-protein-coupled receptor and viral oncogene ORF74, encoded by Kaposi sarcoma-associated herpesvirus (human herpesvirus 8), binds a broad range of chemokines, including CXCL1 (agonist), CXCL8 (neutral ligand), and CXCL10 (inverse agonist). Although chemokines interact with the extracellular N terminus and loops of the receptor, we demonstrate that helix 8 (Hx8) in the intracellular carboxyl tail (C-tail) of ORF74 directs chemokine binding. Partial deletion of the C-tail resulted in a phenotype with reduced constitutive activity but intact regulation by ligands. Complete deletion of the C-tail, including Hx8, resulted in an inactive phenotype that lacks CXCL8 binding sites and has an increased number of binding sites for CXCL10. Similar effects were obtained with the single R7.61(322)W or Q7.62(323)P mutations in Hx8. We propose that the conserved charged or polar side chain at position 7.61 has a specific role in stabilizing the end of transmembrane domain 7 (TM7). Disruption of Hx8 by deletion or mutation distorts an H-bonding network, involving highly conserved amino acids within TM2, TM7, and Hx8, that is crucial for positioning of the TM domains, coupling to Galphaq, and CXCL8 binding. Thus, Hx8 appears to exert a key role in receptor stabilization through the conserved residue R7.61, directing the ligand binding profile of ORF74 and likely also that of other class A G-protein-coupled receptors.

Crystal structure of a photoactivated deprotonated intermediate of rhodopsin

The changes that lead to activation of G protein-coupled receptors have not been elucidated at the structural level. In this work we report the crystal structures of both ground state and a photoactivated deprotonated intermediate of bovine rhodopsin at a resolution of 4.15 A. In the photoactivated state, the Schiff base linking the chromophore and Lys-296 becomes deprotonated, reminiscent of the G protein-activating state, metarhodopsin II. The structures reveal that the changes that accompany photoactivation are smaller than previously predicted for the metarhodopsin II state and include changes on the cytoplasmic surface of rhodopsin that possibly enable the coupling to its cognate G protein, transducin. Furthermore, rhodopsin forms a potentially physiologically relevant dimer interface that involves helices I, II, and 8, and when taken with the prior work that implicates helices IV and V as the physiological dimer interface may account for one of the interfaces of the oligomeric structure of rhodopsin seen in the membrane by atomic force microscopy. The activation and oligomerization models likely extend to the majority of other G protein-coupled receptors.

Bacteriorhodopsin chimeras containing the third cytoplasmic loop of bovine rhodopsin activate transducin for GTP/GDP exchange

The mechanisms by which G-protein-coupled receptors (GPCRs) activate G-proteins are not well understood due to the lack of atomic structures of GPCRs in an active form or in GPCR/G-protein complexes. For study of GPCR/G-protein interactions, we have generated a series of chimeras by replacing the third cytoplasmic loop of a scaffold protein bacteriorhodopsin (bR) with various lengths of cytoplasmic loop 3 of bovine rhodopsin (Rh), and one such chimera containing loop 3 of the human beta2-adrenergic receptor. The chimeras expressed in the archaeon Halobacterium salinarum formed purple membrane lattices thus facilitating robust protein purification. Retinal was correctly incorporated into the chimeras, as determined by spectrophotometry. A 2D crystal (lattice) was evidenced by circular dichroism analysis, and proper organization of homotrimers formed by the bR/Rh loop 3 chimera Rh3C was clearly illustrated by atomic force microscopy. Most interestingly, Rh3C (and Rh3G to a lesser extent) was functional in activation of GTPgamma35S/GDP exchange of the transducin alpha subunit (Galphat) at a level 3.5-fold higher than the basal exchange. This activation was inhibited by GDP and by a high-affinity peptide analog of the Galphat C terminus, indicating specificity in the exchange reaction. Furthermore, a specific physical interaction between the chimera Rh3C loop 3 and the Galphat C terminus was demonstrated by cocentrifugation of transducin with Rh3C. This Galphat-activating bR/Rh chimera is highly likely to be a useful tool for studying GPCR/G-protein interactions.

Computational photochemistry of retinal proteins

High spectral tunability and quantum yield are the striking features of rhodopsin photochemistry. They rely on a strong and complex interaction of their chromophore, the protonated Schiff base of retinal, with its protein environment. In this article, we review the progress in the computational modeling of these systems, focusing on the optical properties and the excited state dynamics. While the earlier success of atomistic theoretical models was based on the breakthrough in X-ray crystallography and combined quantum mechanical molecular mechanical (QM/MM) methodology, recent advances point out the importance of high-level QM methods and the incorporation of effects that are neglected in conventional QM/MM or ONIOM schemes, like polarization and charge transfer.

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.

Ligand-supported homology modeling of the human angiotensin II type 1 (AT1) receptor: Insights into the molecular determinants of telmisartan binding

Angiotensin II type 1 (AT(1)) receptor belongs to the super-family of G-protein-coupled receptors, and antagonists of the AT(1) receptor are effectively used in the treatment of hypertension. To understand the molecular interactions of these antagonists, such as losartan and telmisartan, with the AT(1) receptor, a homology model of the human AT(1) (hAT(1)) receptor with all connecting loops was constructed from the 2.6 A resolution crystal structure (PDB i.d., 1L9H) of bovine rhodopsin. The initial model generated by MODELLER was subjected to a stepwise ligand-supported model refinement. This protocol involved initial docking of non-peptide AT(1) antagonists in the putative binding site, followed by several rounds of iterative energy minimizations and molecular dynamics simulations. The final model was validated based on its correlation with several structure-activity relationships and site-directed mutagenesis data. The final model was also found to be in agreement with a previously reported AT(1) antagonist pharmacophore model. Docking studies were performed for a series of non-peptide AT(1) receptor antagonists in the active site of the final hAT(1) receptor model. The docking was able to identify key molecular interactions for all the AT(1) antagonists studied. Reasonable correlation was observed between the interaction energy values and the corresponding binding affinities of these ligands, providing further validation for the model. In addition, an extensive unrestrained molecular dynamics simulation showed that the docking-derived bound pose of telmisartan is energetically stable. Knowledge gained from the present studies can be used in structure-based drug design for developing novel ligands for the AT(1) receptor.

Ab initio computational modeling of loops in G-protein-coupled receptors: lessons from the crystal structure of rhodopsin

With the help of the crystal structure of rhodopsin an ab initio method has been developed to calculate the three-dimensional structure of the loops that connect the transmembrane helices (TMHs). The goal of this procedure is to calculate the loop structures in other G-protein coupled receptors (GPCRs) for which only model coordinates of the TMHs are available. To mimic this situation a construct of rhodopsin was used that only includes the experimental coordinates of the TMHs while the rest of the structure, including the terminal domains, has been removed. To calculate the structure of the loops a method was designed based on Monte Carlo (MC) simulations which use a temperature annealing protocol, and a scaled collective variables (SCV) technique with proper structural constraints. Because only part of the protein is used in the calculations the usual approach of modeling loops, which consists of finding a single, lowest energy conformation of the system, is abandoned because such a single structure may not be a representative member of the native ensemble. Instead, the method was designed to generate structural ensembles from which the single lowest free energy ensemble is identified as representative of the native folding of the loop. To find the native ensemble a successive series of SCV-MC simulations are carried out to allow the loops to undergo structural changes in a controlled manner. To increase the chances of finding the native funnel for the loop, some of the SCV-MC simulations are carried out at elevated temperatures. The native ensemble can be identified by an MC search starting from any conformation already in the native funnel. The hypothesis is that native structures are trapped in the conformational space because of the high-energy barriers that surround the native funnel. The existence of such ensembles is demonstrated by generating multiple copies of the loops from their crystal structures in rhodopsin and carrying out an extended SCV-MC search. For the extracellular loops e1 and e3, and the intracellular loop i1 that were used in this work, the procedure resulted in dense clusters of structures with Calpha-RMSD approximately 0.5 angstroms. To test the predictive power of the method the crystal structure of each loop was replaced by its extended conformations. For e1 and i1 the procedure identifies native clusters with Calpha-RMSD approximately 0.5 angstroms and good structural overlap of the side chains; for e3, two clusters were found with Calpha-RMSD approximately 1.1 angstroms each, but with poor overlap of the side chains. Further searching led to a single cluster with lower Calpha-RMSD but higher energy than the two previous clusters. This discrepancy was found to be due to the missing elements in the constructs available from experiment for use in the calculations. Because this problem will likely appear whenever parts of the structural information are missing, possible solutions are discussed.

Retinitis pigmentosa associated with rhodopsin mutations: Correlation between phenotypic variability and molecular effects

Similar retinitis pigmentosa (RP) phenotypes can result from mutations affecting different rhodopsin regions, and distinct amino acid substitutions can cause different RP severity and progression rates. Specifically, both the R135L and R135W mutations (cytoplasmic end of H3) result in diffuse, severe disease (class A), but R135W causes more severe and more rapidly progressive RP than R135L. The P180A and G188R mutations (second intradiscal loop) exhibit a mild phenotype with regional variability (class B1) and diffuse disease of moderate severity (class B2), respectively. Computational and in vitro studies of these mutants provide molecular insights into this phenotypic variability.

Metarhodopsin-II stabilization by crosslinked Gtalpha C-terminal peptides and implications for the mechanism of GPCR-G protein coupling

Metarhodopsin-II is the light-excited form of rhodopsin that triggers G protein function. Metarhodopsin-II is stabilized when the N-terminus of the carboxyl (340-350) tail peptide of the alpha-subunit of transducin (Gtalpha) is crosslinked to rhodopsin cysteine 140 or the 340-350 peptide C-terminus of Gtalpha is crosslinked to rhodopsin cysteine 316. When the N-terminus of the peptide is coupled to C316 the MI<-->MII equilibrium is not affected. The evidence suggests that the N-terminus of the 340-350 region of Gtalpha is located near C140 when transducin stabilizes metarhodopsin-II and alternative explanations are suggested for the effectiveness of the 340-350 Gtalpha tail peptide when bound to C316.

Methyl Substituents at the 11 or 12 Position of Retinal Profoundly and Differentially Affect Photochemistry and Signalling Activity of Rhodopsin

The C-11=C-12 double bond of the retinylidene chromophore of rhodopsin holds a central position in its light-induced photoisomerization and hence the photosensory function of this visual pigment. To probe the local environment of the HC-11=C-12H element we have prepared the 11-methyl and 12-methyl derivatives of 11-Z retinal and incorporated these into opsin to generate the rhodopsin analogs 11-methyl and 12-methyl rhodopsin. These analog pigments form with much slower kinetics and lower efficiency than the native pigment. The initial photochemistry and the signaling activity of the analog pigments were investigated by UV-vis and FTIR spectroscopy, and by a G protein activation assay. Our data indicate that the ultrafast formation of the first photointermediate is strongly perturbed by the presence of an 11-methyl substituent, but much less by a 12-methyl substituent. These results support the current concept of the mechanism of the primary photoisomerization event in rhodopsin. An important stronghold of this concept is an out-of-plane movement of the C-12H element, which is facilitated by torsion as well as extended positive charge delocalization into the C-10-C-13 segment of the chromophore. We argue that this mechanism is maintained principally with a methyl substituent at C-12. In addition, we show that both an 11-methyl and a 12-methyl substitutent perturb the photointermediate cascade and finally yield a low-activity state of the receptor. The 11-methyl pigment retains about 30% of the G protein activation rate of native rhodopsin, while the 12-methyl chromophore behaves like an inverse agonist up to at least 20 degrees C, trapping the protein in a perturbed Meta-I-like conformation. We conclude that the isomerization region of the chromophore and the spatial structure of the binding site are finely tuned, in order to achieve a high photosensory potential with an efficient pathway to a high-activity state.

Signal Transfer from GPCRs to G Proteins

Catalysis of nucleotide exchange in heterotrimeric G proteins (Galphabetagamma) is a key step in cellular signal transduction mediated by G protein-coupled receptors. The Galpha N terminus with its helical stretch is thought to be crucial for G protein/activated receptor (R(*)) interaction. The N-terminal fatty acylation of Galpha is important for membrane targeting of G proteins. By applying biophysical techniques to the rhodopsin/transducin model system, we studied the effect of N-terminal truncations in Galpha. In Galphabetagamma, lack of the fatty acid and Galpha truncations up to 33 amino acids had little effect on R(*) binding and R(*)-catalyzed nucleotide exchange, implying that this region is not mandatory for R(*)/Galphabetagamma interaction. However, when the other hydrophobic modification of Galphabetagamma, the Ggamma C-terminal farnesyl moiety, is lacking, R(*) interaction requires the fatty acylated Galpha N terminus. This suggests that the two hydrophobic extensions can replace each other in the interaction of Galphabetagamma with R(*). We propose that in native Galphabetagamma, these two terminal regions are functionally redundant and form a microdomain that serves both to anchor the G protein to the membrane and to establish an initial docking complex with R(*). Accordingly, we find that the native fatty acylated Galpha is competent to interact with R(*) even in the absence of Gbetagamma, whereas nonacylated Galpha requires Gbetagamma for interaction. Experiments with N-terminally truncated Galpha subunits suggest that in the second step of the catalytic process, the receptor binds to the alphaN/beta1-loop region of Galpha to reduce nucleotide affinity and to make the Galpha C terminus available for subsequent interaction with R(*).

Inactive and active states and supramolecular organization of GPCRs: insights from computational modeling

Herein we make an overview of the results of our computational experiments aimed at gaining insight into the molecular mechanisms of GPCR functioning either in their normal conditions or when hit by gain-of-function or loss-of-function mutations. Molecular simulations of a number of GPCRs in their wild type and mutated as well as free and ligand-bound forms were instrumental in inferring the structural features, which differentiate the mutation- and ligand-induced active from the inactive states. These features essentially reside in the interaction pattern of the E/DRY arginine and in the degree of solvent exposure of selected cytosolic domains. Indeed, the active states differ from the inactive ones in the weakening of the interactions made by the highly conserved arginine and in the increase in solvent accessibility of the cytosolic interface between helices 3 and 6. Where possible, the structural hallmarks of the active and inactive receptor states are translated into molecular descriptors useful for in silico functional screening of novel receptor mutants or ligands. Computational modeling of the supramolecular organization of GPCRs and their intracellular partners is the current challenge toward a deep understanding of their functioning mechanisms.

A combined ligand-based and target-based drug design approach for G-protein coupled receptors: application to salvinorin A, a selective kappa opioid receptor agonist

Combined ligand-based and target-based drug design approaches provide a synergistic advantage over either method individually. Therefore, we set out to develop a powerful virtual screening model to identify novel molecular scaffolds as potential leads for the human KOP (hKOP) receptor employing a combined approach. Utilizing a set of recently reported derivatives of salvinorin A, a structurally unique KOP receptor agonist, a pharmacophore model was developed that consisted of two hydrogen bond acceptor and three hydrophobic features. The model was cross-validated by randomizing the data using the CatScramble technique. Further validation was carried out using a test set that performed well in classifying active and inactive molecules correctly. Simultaneously, a bovine rhodopsin based "agonist-bound" hKOP receptor model was also generated. The model provided more accurate information about the putative binding site of salvinorin A based ligands. Several protein structure-checking programs were used to validate the model. In addition, this model was in agreement with the mutation experiments carried out on KOP receptor. The predictive ability of the model was evaluated by docking a set of known KOP receptor agonists into the active site of this model. The docked scores correlated reasonably well with experimental pK (i) values. It is hypothesized that the integration of these two independently generated models would enable a swift and reliable identification of new lead compounds that could reduce time and cost of hit finding within the drug discovery and development process, particularly in the case of GPCRs.

Ab initio computational modeling of long loops in G-protein coupled receptors

A newly developed approach for predicting the structure of segments that connect known elements of secondary structure in proteins has been applied to some of the longer loops in the G-protein coupled receptors (GPCRs) rhodopsin and the dopamine receptor D2R. The algorithm uses Monte Carlo (MC) simulation in a temperature annealing protocol combined with a scaled collective variables (SCV) technique to search conformation space for loop structures that could belong to the native ensemble. Except for rhodopsin, structural information is only available for the transmembrane helices (TMHs), and therefore the usual approach of finding a single conformation of lowest energy has to be abandoned. Instead the MC search aims to find the ensemble located at the absolute minimum free energy, i.e., the native ensemble. It is assumed that structures in the native ensemble can be found by an MC search starting from any conformation in the native funnel. The hypothesis is that native structures are trapped in this part of conformational space because of the high-energy barriers that surround the native funnel. In this work it is shown that the crystal structure of the second extracellular loop (e2) of rhodopsin is a member of this loop's native ensemble. In contrast, the crystal structure of the third intracellular loop is quite different in the different crystal structures that have been reported. Our calculations indicate, that of three crystal structures examined, two show features characteristic of native ensembles while the other one does not. Finally the protocol is used to calculate the structure of the e2 loop in D2R. Here, the crystal structure is not known, but it is shown that several side chains that are involved in interaction with a class of substituted benzamides assume conformations that point into the active site. Thus, they are poised to interact with the incoming ligand.

Comparison of the retinitis pigmentosa mutations in rhodopsin with a functional map of the C5a receptor

We compare the known retinitis pigmentosa (RP) mutations in rhodopsin with mutational data obtained for the complement factor 5a receptor (C5aR), a member of the rhodopsin-like family of G protein-coupled receptors (GPCRs). We have performed genetic analyses that define residues that are required for C5aR folding and function. The cognate residues in rhodopsin are not preferentially mutated in RP, suggesting that the predominant molecular defect in RP involves more than simple misfolding or inactivation. Energy calculations are performed to elucidate the structural effects of the RP mutations. Many of these mutations specifically disrupt the environment of the retinal prosthetic group of rhodopsin, and these do not correspond to essential residues in C5aR. This may be because a retinal group is present in rhodopsin but not in C5aR. Another subset of RP mutations is more generally important for receptor structure, and these mutations correlate with essential residues of C5aR.

Cell-free expression as an emerging technique for the large scale production of integral membrane protein

Membrane proteins are highly underrepresented in structural data banks due to tremendous difficulties that occur upon approaching their structural analysis. Inefficient sample preparation from conventional cellular expression systems is in many cases the first major bottleneck. Preparative scale cell-free expression has now become an emerging alternative tool for the high level production of integral membrane proteins. Many toxic effects attributed to the overproduction of recombinant proteins are eliminated by cell-free expression as viable host cells are no longer required. A unique characteristic is the open nature of cell-free systems that offers a variety of options to manipulate the reaction conditions in order to protect or to stabilize the synthesized recombinant proteins. Detergents or lipids can easily be supplemented and membrane proteins can therefore be synthesized directly into a defined hydrophobic environment of choice that permits solubility and allows the functional folding of the proteins. Alternatively, cell-free produced precipitates of membrane proteins can efficiently be solubilized in mild detergents after expression. Highly valuable for structural approaches is the fast and efficient cell-free production of uniformly or specifically labeled proteins. A considerable number of membrane proteins from diverse families like prokaryotic small multidrug transporters or eukaryotic G-protein coupled receptors have been produced in cell-free systems in high amounts and in functionally active forms. We will give an overview about the current state of the art of this new approach with special emphasis on technical aspects as well as on the functional and structural characterization of cell-free produced membrane proteins.

The structural evolution of dopamine D3 receptor ligands: structure-activity relationships and selected neuropharmacological aspects

"Evolution consists largely of molecular tinkering."-Following the famous concept of the molecular geneticist and medicine Nobel laureate François Jacob, in this review we describe the structural evolution of dopamine D3 receptor ligands from the natural agonist dopamine (DA) to highly potent and subtype selective new agents by bioisosteric tinkering with well-established and privileged or novel and fancy chemical functionalities and scaffolds. Some of the more than 200 ligands presented herein have already achieved therapeutic or scientific value up to now, some will most likely achieve it in the future. Hence, great importance is not only attached to the relationship between structure and activity of the ligands, but also to their utility as pharmacological tools in animal models or as therapeutics in patients with neurological diseases or other disorders.

Local peptide movement in the photoreaction intermediate of rhodopsin

Photoactivation of the visual rhodopsin, a prototypical G protein-coupled receptor (GPCR), involves efficient conversion of the intrinsic inverse-agonist 11-cis-retinal to the all-trans agonist. This event leads to the rearrangement of the heptahelical transmembrane bundle, which is thought to be shared by hundreds of GPCRs. To examine this activation mechanism, we determined the x-ray crystallographic model of the photoreaction intermediate of rhodopsin, lumirhodopsin, which represents the conformational state having the nearly complete all-trans agonist form of the retinal. A difference electron density map clearly indicated that the distorted all-trans-retinal in the precedent intermediate bathorhodopsin relaxes by dislocation of the beta-ionone ring in lumirhodopsin, along with significant peptide displacement in the middle of helix III, including approximately two helical turns. This local movement results in the breaking of the electrostatic interhelical restraints mediated by many of the conserved residues among rhodopsin-like GPCRs, with consequent acquisition of full activity.

A cannabinoid receptor 1 mutation proximal to the DRY motif results in constitutive activity and reveals intramolecular interactions involved in receptor activation

Activation of a G-protein-coupled receptor involves changes in specific microdomain interactions within the transmembrane region of the receptor. Here, we have focused on the role of L207, proximal to the DRY motif of the human cannabinoid receptor 1 in the interconversion of the receptor resting and active states. Ligand binding analysis of the mutant receptor L207A revealed an enhanced affinity for agonists (three- to six-fold) and a diminished affinity for inverse agonists (19- to 35-fold) compared to the wild-type receptor, properties characteristic of constitutive activity. To further examine whether this mutant adopts a ligand-independent, active form, treatment with GTPgammaS was used to inhibit G protein coupling. Under these conditions, the L207A receptor exhibited a 10-fold increase in affinity for the inverse agonist SR141716A, consistent with a shift away from an enhanced precoupled state. Analysis of the cellular activity of the L207A receptor showed elevated basal cyclic AMP accumulation relative to the wild type that is inhibited by SR141716A, consistent with receptor-mediated Gs precoupling. Using toxins to selectively abrogate Gs or Gi coupling, we found that CP55940 nonetheless induced only a Gi response suggesting a strong preference of this ligand-bound form for Gi in this system. Molecular dynamics simulations reveal that the single residue change of L207A impacts the association of TM3 and TM6 in the receptor by altering hydrophobic interactions involving L207, the salt bridge involving the Arg of the DRY motif, and the helical structure of TM6, consistent with events leading to activation. The structural alterations parallel those observed in models of a mutant CB(1) receptor T210I, with established constitutive activity (D'Antona, A.M., Ahn, K.H. and Kendall, D.A., 2006. Mutations of CB1 T210 produce active and inactive receptor forms: correlations with ligand affinity, receptor stability, and cellular localization. Biochemistry, 45, 5606-5617).

Proposal of a New Binding Orientation for Non-Peptide AT1 Antagonists: Homology Modeling, Docking and Three-Dimensional Quantitative Structure−Activity Relationship Analysis

A three-dimensional model of the AT1 receptor was constructed by means of a homology modeling procedure, using the X-ray structure of bovine rhodopsin as the initial template and taking into account the available site-directed mutagenesis data. The docking of losartan and its active metabolite EXP3174, followed by 1 ns of molecular dynamics (MD) simulation inserted into the phospholipid bilayer, suggested a different binding orientation for these antagonists from those previously proposed. Furthermore, the docking of several non-peptide antagonists was used as an alignment tool for the development of a three-dimensional quantitative structure-activity relationship (3D-QSAR) model, and the good results confirmed our binding hypothesis and the reliability of the model.

Dynamic models of G-protein coupled receptor dimers: indications of asymmetry in the rhodopsin dimer from molecular dynamics simulations in a POPC bilayer

Based on the growing evidence that G-protein coupled receptors (GPCRs) form homo- and hetero-oligomers, models of GPCR signaling are now considering macromolecular assemblies rather than monomers, with the homo-dimer regarded as the minimal oligomeric arrangement required for functional coupling to the G-protein. The dynamic mechanisms of such signaling assemblies are unknown. To gain some insight into properties of GPCR dimers that may be relevant to functional mechanisms, we study their current structural prototype, rhodopsin. We have carried out nanosecond time-scale molecular dynamics (MD) simulations of a rhodopsin dimer and compared the results to the monomer simulated in the same type of bilayer membrane model composed of an equilibrated unit cell of hydrated palmitoyl-oleoyl-phosphatidyl choline (POPC). The dynamic representation of the homo-dimer reveals the location of structural changes in several regions of the monomeric subunits. These changes appear to be more pronounced at the dimerization interface that had been shown to be involved in the activation process [Proc Natl Acad Sci USA 102:17495, 2005]. The results are consistent with a model of GPCR activation that involves allosteric modulation through a single GPCR subunit per dimer.

Proton movement and photointermediate kinetics in rhodopsin mutants

The role of ionizable amino acid side chains in the bovine rhodopsin activation mechanism was studied in mutants E134Q, E134R/R135E, H211F, and E122Q. All mutants exhibited bathorhodopsin stability on the 30 ns to 1 micros time scale similar to that of the wild type. Lumirhodopsin decay was also similar to that of the wild type except for the H211F mutant where early decay (20 micros) to a second form of lumirhodopsin was seen, followed by formation of an extremely long-lived Meta I(480) product (34 ms), an intermediate which forms to a much reduced extent, if at all, in dodecyl maltoside suspensions of wild-type rhodopsin. A smaller amount of a similar long-lived Meta I(480) product was seen after photolysis of E122Q, but E134Q and E134R/R135Q displayed kinetics much more similar to those of the wild type under these conditions (i.e., no Meta I(480) product). These results support the idea that specific interaction of His211 and Glu122 plays a significant role in deprotonation of the retinylidene Schiff base and receptor activation. Proton uptake measurements using bromcresol purple showed that E122Q was qualitatively similar to wild-type rhodopsin, with at least one proton being released during lumirhodopsin decay per Meta I(380) intermediate formed, followed by uptake of at least two protons per rhodopsin bleached on a time scale of tens of milliseconds. Different results were obtained for H211F, E134Q, and E134R/R135E, which all released approximately two protons per rhodopsin bleached. These results show that several ionizable groups besides the Schiff base imine are affected by the structural changes involved in rhodopsin activation. At least two proton uptake groups and probably at least one proton release group in addition to the Schiff base are present in rhodopsin.

Metal Ion Site Engineering Indicates a Global Toggle Switch Model for Seven-transmembrane Receptor Activation

Much evidence indicates that, during activation of seven-transmembrane (7TM) receptors, the intracellular segments of the transmembrane helices (TMs) move apart with large amplitude, rigid body movements of especially TM-VI and TM-VII. In this study, AspIII:08 (Asp113), the anchor point for monoamine binding in TM-III, was used as the starting point to engineer activating metal ion sites between the extracellular segments of the beta2-adrenergic receptor. Cu(II) and Zn(II) alone and in complex with aromatic chelators acted as potent (EC50 decreased to 0.5 microm) and efficacious agonists in sites constructed between positions III:08 (Asp or His), VI:16 (preferentially Cys), and/or VII:06 (preferentially Cys). In molecular models built over the backbone conformation of the inactive rhodopsin structure, the heavy atoms that coordinate the metal ion were located too far away from each other to form high affinity metal ion sites in both the bidentate and potential tridentate settings. This indicates that the residues involved in the main ligand-binding pocket will have to move closer to each other during receptor activation. On the basis of the distance constraints from these activating metal ion sites, we propose a global toggle switch mechanism for 7TM receptor activation in which inward movement of the extracellular segments of especially TM-VI and, to some extent, TM-VII is coupled to the well established outward movement of the intracellular segments of these helices. We suggest that the pivots for these vertical seesaw movements are the highly conserved proline bends of the involved helices.