Structure of rhodopsin

A new era of GPCR structural and chemical biology

G protein-coupled receptors (GPCRs) are versatile molecular machines that regulate the majority of physiological responses to chemically diverse hormones and neurotransmitters. Recent breakthroughs in structural studies have advanced our understanding of GPCR signaling, particularly the selectivity of ligand recognition and receptor activation of G proteins.

Identification of anti-malarial compounds as novel antagonists to chemokine receptor CXCR4 in pancreatic cancer cells

Despite recent advances in targeted therapies, patients with pancreatic adenocarcinoma continue to have poor survival highlighting the urgency to identify novel therapeutic targets. Our previous investigations have implicated chemokine receptor CXCR4 and its selective ligand CXCL12 in the pathogenesis and progression of pancreatic intraepithelial neoplasia and invasive pancreatic cancer; hence, CXCR4 is a promising target for suppression of pancreatic cancer growth. Here, we combined in silico structural modeling of CXCR4 to screen for candidate anti-CXCR4 compounds with in vitro cell line assays and identified NSC56612 from the National Cancer Institute's (NCI) Open Chemical Repository Collection as an inhibitor of activated CXCR4. Next, we identified that NSC56612 is structurally similar to the established anti-malarial drugs chloroquine and hydroxychloroquine. We evaluated these compounds in pancreatic cancer cells in vitro and observed specific antagonism of CXCR4-mediated signaling and cell proliferation. Recent in vivo therapeutic applications of chloroquine in pancreatic cancer mouse models have demonstrated decreased tumor growth and improved survival. Our results thus provide a molecular target and basis for further evaluation of chloroquine and hydroxychloroquine in pancreatic cancer. Historically safe in humans, chloroquine and hydroxychloroquine appear to be promising agents to safely and effectively target CXCR4 in patients with pancreatic cancer.

Structure prediction of G protein-coupled receptors and their ensemble of functionally important conformations

G protein-coupled receptors (GPCRs) are integral membrane proteins whose "pleiotropic" nature enables transmembrane (TM) signal transduction, amplification, and diversification via G protein-coupled and β arrestin-coupled pathways. GPCRs appear to enable this by being structurally flexible and by existing in different conformational states with potentially different signaling and functional consequences. We describe a method for the prediction of the three-dimensional structures of these different conformations of GPCRs starting from their amino acid sequence. It combines a unique protocol of computational methods that first predict the TM regions of these receptors and TM helix shapes based on those regions, which is followed by a locally complete sampling of TM helix packings and their scoring that results in a few (~10-20) lowest energy conformations likely to play a role in binding to different ligands and signaling events. Prediction of the structures for multiple conformations of a GPCR is starting to enable the testing of multiple hypotheses related to GPCR activation and binding to ligands with different signaling profiles.

Predicted structures and dynamics for agonists and antagonists bound to serotonin 5-HT2B and 5-HT2C receptors

Subtype 2 serotonin (5-hydroxytryptamine, 5-HT) receptors are major drug targets for schizophrenia, feeding disorders, perception, depression, migraines, hypertension, anxiety, hallucinogens, and gastrointestinal dysfunctions. (1) We report here the predicted structure of 5-HT2B and 5-HT2C receptors bound to highly potent and selective 5-HT2B antagonist PRX-08066 3, (pKi: 30 nM), including the key binding residues [V103 (2.53), L132 (3.29), V190 (4.60), and L347 (6.58)] determining the selectivity of binding to 5-HT2B over 5-HT2A. We also report structures of the endogenous agonist (5-HT) and a HT2B selective antagonist 2 (1-methyl-1-1,6,7,8-tetrahydro-pyrrolo[2,3-g]quinoline-5-carboxylic acid pyridine-3-ylamide). We examine the dynamics for the agonist- and antagonist-bound HT2B receptors in explicit membrane and water finding dramatically different patterns of water migration into the NPxxY motif and the binding site that correlates with the stability of ionic locks in the D(E)RY region.

Importance of receptor flexibility in binding of cyclam compounds to the chemokine receptor CXCR4

We have elucidated the binding sites of four moncyclam and one bicyclam antagonist AMD3100, in the human chemokine receptor CXCR4. Using the predicted structural models of CXCR4, we have further predicted the binding sites of these cyclam compounds. We used the computational method LITiCon to map the differences in receptor structure stabilized by the mono and bicyclam compounds. Accounting for the receptor flexibility lead to a single binding mode for the cyclam compounds, that has not been possible previously using a single receptor structural model and fixed receptor docking algorithms. There are several notable differences in the receptor conformations stabilized by monocyclam antagonist compared to a bicylam antagonist. The loading of the Cu(2+) ions in the cyclam compounds, shrinks the size of the cyclam rings and the residue D262(6.58) plays an important role in bonding to the copper ion in the monocylam compounds while residue E288(7.39) is important for the bicyclam compound.

Prediction of the three-dimensional structure for the rat urotensin II receptor, and comparison of the antagonist binding sites and binding selectivity between human and rat receptors from atomistic simulations

Urotensin-II (U-II) has been shown to be the most potent mammalian vasoconstrictor known. Thus, a U-II antagonist might be of therapeutic value in a number of cardiovascular disorders. However, interspecies variability of several nonpeptidic ligands complicates the interpretation of in vivo studies of such antagonists in preclinical animal disease models. ACT058362 is a selective antagonist for the human U-II receptor (hUT2R) with a reported K(d) value of approximately 4 nM in a molecular binding assay, but it is reported to bind weakly to rat UT2R (rUT2R), with a K(d) value of approximately 1 500 nM. In contrast, the arylsulphonamide SB706375 is a selective antagonist against both hUT2R (K(d)= approximately 9 nM) and rUT2R (K(d)= approximately 21 nM). To understand the species selectivity of the UT2R, we investigated the binding site of ACT058362 and SB706375 in both hUT2R and rUT2R to explain the dramatically lower (approximately 400-fold) affinity of ACT058362 for rUT2R and the similar affinity (approximately 10 nM) of SB706375 for both UT2Rs. These studies used MembStruk and MSCDock to predict the UT2R structure and the binding site of ACT058362 and SB706375. Based on binding energies, we found two binding modes each with D130(3.32) as the crucial anchoring point (Ballesteros-Weinstein numbering given in superscript). We predict that ACT058362 (an aryl-amine-aryl or ANA ligand) binds in the transmembrane (TM) 3456 region, while SB706375 (an aryl-aryl-amine or AAN ligand) binds in the TM 1237 region. These predicted sites explain the known differences in binding of the ANA ligand to rat and human receptors, while explaining the similar binding of the AAN compound to rat and human receptors. Moreover the predictions explain currently available structure-activity relationship (SAR) data. To further validate the predicted binding sites of these ligands in hUT2R and rUT2R, we propose several mutations that would help define the structural origins of differential responses between UT2R of different species, potentially indicating novel UT2R antagonists with cross-species high affinity.

Discrimination of saturated aldehydes by the rat I7 olfactory receptor

The discrimination of n-alkyl-saturated aldehydes during the early stage of odorant recognition by the rat I7 olfactory receptor (OR-I7) is investigated. The concentrations of odorants necessary for 50% activation (or inhibition) of the OR-I7 are measured by calcium imaging recordings of dissociated rat olfactory sensory neurons, expressing the recombinant OR-I7 from an adenoviral vector. These are correlated with the corresponding binding free energies computed for a homology structural model of the OR-I7 built from the crystal structure of bovine visual rhodopsin at 2.2 A resolution.

Allosteric antagonist binding sites in class B GPCRs: corticotropin receptor 1

The 41 amino acid neuropeptide, corticotropin-releasing factor (CRF) and its associated receptors CRF(1)-R and CRF(2)-R have been targeted for treating stress related disorders. Both CRF(1)-R and CRF(2)-R belong to the class B G-protein coupled receptors for which little information is known regarding the small molecule antagonist binding characteristics. However, it has been shown recently that different non-peptide allosteric ligands stabilize different receptor conformations for CRF(1)-R and hence an understanding of the ligand induced receptor conformational changes is important in the pharmacology of ligand binding. In this study, we modeled the receptor and identified the binding sites of representative small molecule allosteric antagonists for CRF(1)-R. The predicted binding sites of the investigated compounds are located within the transmembrane (TM) domain encompassing TM helices 3, 5 and 6. The docked compounds show strong interactions with H228 on TM3 and M305 on TM5 that have also been implicated in the binding by site directed mutation studies. H228 forms a hydrogen bond of varied strengths with all the antagonists in this study and this is in agreement with the decreased binding affinity of several compounds with H228F mutation. Also mutating M305 to Ile showed a sharp decrease in the calculated binding energy whereas the binding energy loss on M305 to Leu was less significant. These results are in qualitative agreement with the decrease in binding affinities observed experimentally. We further predicted the conformational changes in CRF(1)-R induced by the allosteric antagonist NBI-27914. Movement of TM helices 3 and 5 are dominant and generates three degenerate conformational states two of which are separated by an energy barrier from the third, when bound to NBI-27914. Binding of NBI-27914 was predicted to improve the interaction of the ligand with M305 and also enhanced the aromatic stacking between the ligand and F232 on TM3. A virtual ligand screening of ~13,000 compounds seeded with ~350 CRF(1)-R specific active antagonists performed on the NBI-27914 stabilized conformation of CRF(1)-R yielded a 44% increase in enrichment compared to the initially modeled receptor conformation at a 10% cutoff. The NBI-27914 stabilized conformation also shows a high enrichment for high affinity antagonists compared to the weaker ones. Thus, the conformational changes induced by NBI-27914 improved the ligand screening efficiency of the CRF(1)-R model and demonstrate a generalized application of the method in drug discovery.

The year in G protein-coupled receptor research

I (R.P.M.) presented "The Year In G Protein-Coupled Receptor Research" at ENDO 2009. I first described the diversity of ligands and the five families into which the approximately 800 G protein-coupled receptors (GPCRs) are grouped, their basic structural architectures, their preeminent role in signaling, and the enormous scope for developing drugs targeted at GPCRs. I then spoke about some of the exciting breakthroughs in solving the atomic level structures of the active state of rhodopsin, beta(2)-adrenergic, beta(1)-adrenergic, and A(2A)-adenosine receptors. I also described studies on the structural changes accompanying the activation of the rhodopsin family of GPCRs. From these recent technical advances, we can anticipate that many more GPCR structures will emerge, which will afford us greater insight into their common and unique structural features and, particularly, the mechanisms underlying their activation. These insights will guide us in our understanding of how GPCRs operate, both in the normal and pathological situation. Although these crystal structures are highly informative, it is important to recognize that they represent static frozen conformations of a single GPCR state. New biophysical techniques are therefore being utilized to facilitate the dynamic monitoring of GPCR structural changes in relation to ligand activation. Solving of the crystal structures of GPCRs has also presented the real possibility of using the information of the ligand-binding pocket to allow in silico screening for novel small-molecule ligands. I then reviewed the concept of ligand-induced selective signaling of GPCRs, which is opening up new insights into more selective drug development. The assembly of GPCRs as homo- and heterooligomers and their phosphorylation and association with a vast array of trafficking and signal-modulating proteins are emerging as major mechanisms underlying the functioning of GPCRs. Differential expression and recruitment of these proteins provide a mechanism for subtle physiological regulation of cellular activity. Finally, I mentioned some of the GPCRs that have lately come to the fore as novel regulators in endocrinology. These included fatty acid-specific GPCRs expressed in pancreatic beta-cells and novel neuroendocrine GPCRs regulating reproduction.

Elucidation of binding sites of dual antagonists in the human chemokine receptors CCR2 and CCR5

Design of dual antagonists for the chemokine receptors CCR2 and CCR5 will be greatly facilitated by knowledge of the structural differences of their binding sites. Thus, we computationally predicted the binding site of the dual CCR2/CCR5 antagonist N-dimethyl-N-[4-[[[2-(4-methylphenyl)-6,7-dihydro-5H-benzohepten-8-yl] carbonyl]amino]benzyl]tetrahydro-2H-pyran-4-aminium (TAK-779), and a CCR2-specific antagonist N-(carbamoylmethyl)-3-trifluoromethyl benzamido-parachlorobenzyl 3-aminopyrrolidine (Teijin compound 1) in an ensemble of predicted structures of human CCR2 and CCR5. Based on our predictions of the protein-ligand interactions, we examined the activity of the antagonists for cells expressing thirteen mutants of CCR2 and five mutants of CCR5. The results show that residues Trp98(2.60) and Thr292(7.40) contribute significantly to the efficacy of both TAK-779 and Teijin compound 1, whereas His121(3.33) and Ile263(6.55) contribute significantly only to the antagonistic effect of Teijin compound 1 at CCR2. Mutation of residues Trp86(2.60) and Tyr108(3.32) adversely affected the efficacy of TAK-779 in antagonizing CCR5-mediated chemotaxis. Y49A(1.39) and E291A(7.39) mutants of CCR2 showed a complete loss of CCL2 binding and chemotaxis, despite robust cell surface expression, suggesting that these residues are critical in maintaining the correct receptor architecture. Modeling studies support the hypothesis that the residues Tyr49(1.39), Trp98(2.60), Tyr120(3.32), and Glu291(7.39) of CCR2 form a tight network of aromatic cluster and polar contacts between transmembrane helices 1, 2, 3, and 7.

The structure of human serotonin 2c G-protein-coupled receptor bound to agonists and antagonists

We used the MembStruk computational procedure to predict the three-dimensional (3D) structure for the serotonin 5-HT(2C) G-protein-coupled receptor (GPCR). Using this structure, we used the MSCDock computational procedure to predict the 3D structures for bound ligand-protein complexes for agonists such as serotonin and antagonists such as ritanserin, metergoline, and methiothepin. In addition, we predicted the SAR data for a series of psilocybin analogs, both agonists and antagonists. We performed molecular dynamics (MD) on serotonin bound to 5-HT(2C) and we find the protein and binding site to be stable after 5ns. We find good agreement with the currently known experimental data and we predict a number of new mutations which could be used to validate further our predicted structures. This agreement between theory and experiment suggests that our 3D structure is sufficiently accurate for use in drug design. We also compare a preliminary prediction for 5-HT(2B) with our prediction for 5-HT(2C) and find a difference in TM5 that contributes to different serotonin binding modes in 5-HT(2B) and 5-HT(2C).

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.

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 predicted 3D structures of the human M1 muscarinic acetylcholine receptor with agonist or antagonist bound

The muscarinic acetylcholine G-protein-coupled receptors are implicated in diseases ranging from cognitive dysfunctions to smooth-muscle disorders. To provide a structural basis for drug design, we used the MembStruk computational method to predict the 3D structure of the human M1 muscarinic receptor. We validated this structure by using the HierDock method to predict the binding sites for three agonists and four antagonists. The intermolecular ligand-receptor contacts at the predicted binding sites agree well with deductions from available mutagenesis experiments, and the calculated relative binding energies correlate with measured binding affinities. The predicted binding site of all four antagonists is located between transmembrane (TM) helices 3, 4, 5, 6, and 7, whereas the three agonists prefer a site involving residues from TM3, TM6, and TM7. We find that Trp 157(4) contributes directly to antagonist binding, whereas Pro 159(4) provides an indirect conformational switch to position Trp 157(4) in the binding site (the number in parentheses indicates the TM helix). This explains the large decrease in ligand binding affinity and signaling efficacy by mutations of Trp 157(4) and Pro 159(4) not previously explained by homology models. We also found that Asp 105(3) and aromatic residues Tyr 381(6), Tyr 404(7), and Tyr 408(7) are critical for binding the quaternary ammonium head group of the ligand through cation-pi interactions. For ligands with a charged tertiary amine head group, we suggest that proton transfer from the ligand to Asp 105(3) occurs upon binding. Furthermore, we found that an extensive aromatic network involving Tyr 106(3), Trp 157(4), Phe 197(5), Trp 378(6), and Tyr 381(6) is important in stabilizing antagonist binding. For antagonists with two terminal phenyl rings, this aromatic network extends to Trp 164(4), Tyr 179(extracellular loop 2), and Phe 390(6) located at the extracellular end of the TMs. We find that Asn 382(6) forms hydrogen bonds with selected antagonists. Tyr381(6) and Ser 109(3) form hydrogen bonds with the ester moiety of acetylcholine, which binds in the gauche conformation.

Structural activation pathways from dynamic olfactory receptor-odorant interactions

We have simulated an odor ligand's dynamic behavior in the binding region of an olfactory receptor (OR). Our short timescale computational studies (up to 200 ps) have helped identify unprecedented postdocking ligand behavior of ligands. From in vacuo molecular dynamics simulations of interactions between models of rat OR I7 and 10 aldehyde ligands, we have identified a dissociative pathway along which the ligand exits and enters the OR-binding pocket--a transit event. The ligand's transit through the receptor's binding region may mark the beginning of a signal transduction cascade leading to odor recognition. We have graphically traced the rotameric changes in key OR amino acid side chains during the transit. Our results have helped substantiate or refute previously held notions of amino acid contribution to ligand stability in the binding pocket. Our observations of ligand activity when compared to those of experimental (electroolfactogram response) OR-activation studies provide a view to predicting the stability of ligands in the binding pocket as a precursor to OR activation by the ligand.

Structural basis for a broad but selective ligand spectrum of a mouse olfactory receptor: mapping the odorant-binding site

The olfactory receptor (OR) superfamily provides a basis for the remarkable ability to recognize and discriminate a large number of odorants. In mice, the superfamily includes approximately 1000 members, and they recognize overlapping sets of odorants with distinct affinities and specificities. To address the molecular basis of odor discrimination by the mammalian OR superfamily, we performed functional analysis on a series of site-directed mutants and performed ligand docking simulation studies to define the odorant-binding site of a mouse OR. Our results indicate that several amino acids in the transmembrane domains formed a ligand-binding pocket. Although other G-protein-coupled receptors (GPCRs) recognize biogenic ligands mainly with ionic or hydrogen bonding interactions, ORs recognize odorants mostly via hydrophobic and van der Waals interactions. This accounts for the broad but selective binding by ORs as well as their relatively low ligand-binding affinities. Furthermore, we succeeded in rational receptor design, inserting point mutations in the odorant-binding site that resulted in predicted changes in ligand specificity and antagonist activity. This ability to rationally design the receptor validated the binding site structure that was deduced with our mutational and ligand docking studies. Such broad and specific sensitivity suggests an evolutionary process during which mutations in the active site led to an enormous number of ORs with a wide range of ligand specificity. The current study reveals the molecular environment of the odorant-binding site, and it further advances the understanding of GPCR pharmacology.

Test of the Binding Threshold Hypothesis for olfactory receptors: explanation of the differential binding of ketones to the mouse and human orthologs of olfactory receptor 912-93

We tested the Binding Threshold Hypothesis (BTH) for activation of olfactory receptors (ORs): To activate an OR, the odorant must bind to the OR with binding energy above some threshold value. The olfactory receptor (OR) 912-93 is known experimentally to be activated by ketones in mouse, but is inactive to ketones in human, despite an amino acid sequence identity of approximately 66%. To investigate the origins of this difference, we used the MembStruk first-principles method to predict the tertiary structure of the mouse OR 912-93 (mOR912-93), and the HierDock first-principles method to predict the binding site for ketones to this receptor. We found that the strong binding of ketones to mOR912-93 is dominated by a hydrogen bond of the ketone carbonyl group to Ser105. All ketones predicted to have a binding energy stronger than EBindThresh = 26 kcal/mol were observed experimentally to activate this OR, while the two ketones predicted to bind more weakly do not. In addition, we predict that 2-undecanone and 2-dodecanone both bind sufficiently strongly to activate mOR912-93. A similar binding site for ketones was predicted in hOR912-93, but the binding is much weaker because the human ortholog has a Gly at the position of Ser105. We predict that mutating this Gly to Ser in human should lead to activation of hOR912-93 by these ketones. Experimental substantiations of the above predictions would provide further tests of the validity of the BTH, our predicted 3D structures, and our predicted binding sites for these ORs.

Essential role for the second extracellular loop in C5a receptor activation

More than 90% of G protein-coupled receptors (GPCRs) contain a disulfide bridge that tethers the second extracellular loop (EC2) to the third transmembrane helix. To determine the importance of EC2 and its disulfide bridge in receptor activation, we subjected this region of the complement factor 5a receptor (C5aR) to random saturation mutagenesis and screened for functional receptors in yeast. The cysteine forming the disulfide bridge was the only conserved residue in the EC2-mutated receptors. Notably, approximately 80% of the functional receptors exhibited potent constitutive activity. These results demonstrate an unexpected role for EC2 as a negative regulator of C5a receptor activation. We propose that in other GPCRs, EC2 might serve a similar role by stabilizing the inactive state of the receptor.

Molecular mechanisms of ligand binding, signaling, and regulation within the superfamily of G-protein-coupled receptors: molecular modeling and mutagenesis approaches to receptor structure and function

The superfamily of G-protein-coupled receptors (GPCRs) could be subclassified into 7 families (A, B, large N-terminal family B-7 transmembrane helix, C, Frizzled/Smoothened, taste 2, and vomeronasal 1 receptors) among mammalian species. Cloning and functional studies of GPCRs have revealed that the superfamily of GPCRs comprises receptors for chemically diverse native ligands including (1) endogenous compounds like amines, peptides, and Wnt proteins (i.e., secreted proteins activating Frizzled receptors); (2) endogenous cell surface adhesion molecules; and (3) photons and exogenous compounds like odorants. The combined use of site-directed mutagenesis and molecular modeling approaches have provided detailed insight into molecular mechanisms of ligand binding, receptor folding, receptor activation, G-protein coupling, and regulation of GPCRs. The vast majority of family A, B, C, vomeronasal 1, and taste 2 receptors are able to transduce signals into cells through G-protein coupling. However, G-protein-independent signaling mechanisms have also been reported for many GPCRs. Specific interaction motifs in the intracellular parts of these receptors allow them to interact with scaffold proteins. Protein engineering techniques have provided information on molecular mechanisms of GPCR-accessory protein, GPCR-GPCR, and GPCR-scaffold protein interactions. Site-directed mutagenesis and molecular dynamics simulations have revealed that the inactive state conformations are stabilized by specific interhelical and intrahelical salt bridge interactions and hydrophobic-type interactions. Constitutively activating mutations or agonist binding disrupts such constraining interactions leading to receptor conformations that associates with and activate G-proteins.

First principles predictions of the structure and function of g-protein-coupled receptors: validation for bovine rhodopsin

G-protein-coupled receptors (GPCRs) are involved in cell communication processes and with mediating such senses as vision, smell, taste, and pain. They constitute a prominent superfamily of drug targets, but an atomic-level structure is available for only one GPCR, bovine rhodopsin, making it difficult to use structure-based methods to design receptor-specific drugs. We have developed the MembStruk first principles computational method for predicting the three-dimensional structure of GPCRs. In this article we validate the MembStruk procedure by comparing its predictions with the high-resolution crystal structure of bovine rhodopsin. The crystal structure of bovine rhodopsin has the second extracellular (EC-II) loop closed over the transmembrane regions by making a disulfide linkage between Cys-110 and Cys-187, but we speculate that opening this loop may play a role in the activation process of the receptor through the cysteine linkage with helix 3. Consequently we predicted two structures for bovine rhodopsin from the primary sequence (with no input from the crystal structure)-one with the EC-II loop closed as in the crystal structure, and the other with the EC-II loop open. The MembStruk-predicted structure of bovine rhodopsin with the closed EC-II loop deviates from the crystal by 2.84 A coordinate root mean-square (CRMS) in the transmembrane region main-chain atoms. The predicted three-dimensional structures for other GPCRs can be validated only by predicting binding sites and energies for various ligands. For such predictions we developed the HierDock first principles computational method. We validate HierDock by predicting the binding site of 11-cis-retinal in the crystal structure of bovine rhodopsin. Scanning the whole protein without using any prior knowledge of the binding site, we find that the best scoring conformation in rhodopsin is 1.1 A CRMS from the crystal structure for the ligand atoms. This predicted conformation has the carbonyl O only 2.82 A from the N of Lys-296. Making this Schiff base bond and minimizing leads to a final conformation only 0.62 A CRMS from the crystal structure. We also used HierDock to predict the binding site of 11-cis-retinal in the MembStruk-predicted structure of bovine rhodopsin (closed loop). Scanning the whole protein structure leads to a structure in which the carbonyl O is only 2.85 A from the N of Lys-296. Making this Schiff base bond and minimizing leads to a final conformation only 2.92 A CRMS from the crystal structure. The good agreement of the ab initio-predicted protein structures and ligand binding site with experiment validates the use of the MembStruk and HierDock first principles' methods. Since these methods are generic and applicable to any GPCR, they should be useful in predicting the structures of other GPCRs and the binding site of ligands to these proteins.

The predicted 3D structure of the human D2 dopamine receptor and the binding site and binding affinities for agonists and antagonists

Dopamine neurotransmitter and its receptors play a critical role in the cell signaling process responsible for information transfer in neurons functioning in the nervous system. Development of improved therapeutics for such disorders as Parkinson's disease and schizophrenia would be significantly enhanced with the availability of the 3D structure for the dopamine receptors and of the binding site for dopamine and other agonists and antagonists. We report here the 3D structure of the long isoform of the human D2 dopamine receptor, predicted from primary sequence using first-principles theoretical and computational techniques (i.e., we did not use bioinformatic or experimental 3D structural information in predicting structures). The predicted 3D structure is validated by comparison of the predicted binding site and the relative binding affinities of dopamine, three known dopamine agonists (antiparkinsonian), and seven known antagonists (antipsychotic) in the D2 receptor to experimentally determined values. These structures correctly predict the critical residues for binding dopamine and several antagonists, identified by mutation studies, and give relative binding affinities that correlate well with experiments. The predicted binding site for dopamine and agonists is located between transmembrane (TM) helices 3, 4, 5, and 6, whereas the best antagonists bind to a site involving TM helices 2, 3, 4, 6, and 7 with minimal contacts to TM helix 5. We identify characteristic differences between the binding sites of agonists and antagonists.

Predicted 3D structure for the human beta 2 adrenergic receptor and its binding site for agonists and antagonists

We report the 3D structure of human beta2 adrenergic receptor (AR) predicted by using the MembStruk first principles method. To validate this structure, we use the HierDock first principles method to predict the ligand-binding sites for epinephrine and norepinephrine and for eight other ligands, including agonists and antagonists to beta 2 AR and ligands not observed to bind to beta 2 AR. The binding sites agree well with available mutagenesis data, and the calculated relative binding energies correlate reasonably with measured binding affinities. In addition, we find characteristic differences in the predicted binding sites of known agonists and antagonists that allow us to infer the likely activity of other ligands. The predicted ligand-binding properties validate the methods used to predict the 3D structure and function. This validation is a successful step toward applying these procedures to predict the 3D structures and function of the other eight subtypes of ARs, which should enable the development of subtype-specific antagonists and agonists with reduced side effects.

Prediction of structure and function of G protein-coupled receptors

G protein-coupled receptors (GPCRs) mediate our sense of vision, smell, taste, and pain. They are also involved in cell recognition and communication processes, and hence have emerged as a prominent superfamily for drug targets. Unfortunately, the atomic-level structure is available for only one GPCR (bovine rhodopsin), making it difficult to use structure-based methods to design drugs and mutation experiments. We have recently developed first principles methods (MembStruk and HierDock) for predicting structure of GPCRs, and for predicting the ligand binding sites and relative binding affinities. Comparing to the one case with structural data, bovine rhodopsin, we find good accuracy in both the structure of the protein and of the bound ligand. We report here the application of MembStruk and HierDock to beta1-adrenergic receptor, endothelial differential gene 6, mouse and rat I7 olfactory receptors, and human sweet receptor. We find that the predicted structure of beta1-adrenergic receptor leads to a binding site for epinephrine that agrees well with the mutation experiments. Similarly the predicted binding sites and affinities for endothelial differential gene 6, mouse and rat I7 olfactory receptors, and human sweet receptor are consistent with the available experimental data. These predicted structures and binding sites allow the design of mutation experiments to validate and improve the structure and function prediction methods. As these structures are validated they can be used as targets for the design of new receptor-selective antagonists or agonists for GPCRs.

Mutations at position 125 in transmembrane helix III of rhodopsin affect the structure and signalling of the receptor

Mutation of L125R in trasmembrane helix III of rhodopsin, associated with the retinal degenerative disease retinitis pigmentosa, was previously shown to cause structural misfolding of the mutant protein. Also, conservative mutations at this position were found to cause partial misfolding of the mutant receptors. We report here on a series of mutations at position 125 to further investigate the role of Leu125 in the correct folding and function of rhodopsin. In particular, the effect of the size of the substituted amino-acid side chain in the functionality of the receptor, measured as the ability of the mutant rhodopsins to activate the G protein transducin, has been analysed. The following mutations have been studied: L125G, L125N, L125I, L125H, L125P, L125T, L125D, L125E, L125Y and L125W. Most of the mutant proteins, expressed in COS-1 cells, showed reduced 11-cis-retinal binding, red-shifts in the wavelength of the visible absorbance maximum, and increased reactivity towards hydroxylamine in the dark. Thermal stability in the dark was reduced, particularly for L125P, L125Y and L125W mutants. The ability of the mutant rhodopsins to activate the G protein transducin was significantly reduced in a size dependent manner, especially in the case of the bulkier L125Y and L125W substitutions, suggesting a steric effect of the substituted amino acid. On the basis of the present and previous results, Leu125 in transmembrane helix III of rhodopsin, in the vicinity of the beta-ionone ring of 11-cis-retinal, is proposed to be an important residue in maintaining the correct structure of the chromophore binding pocket. Thus, bulky substitutions at this position may affect the structure and signallling of the receptor by altering the optimal conformation of the retinal binding pocket, rather than by direct interaction with the chromophore, as seen from the recent crystallographic structure of rhodopsin.

Modeling the 3D structure of GPCRs from sequence

G-protein-coupled receptors (GPCRs) are a large and functionally diverse protein superfamily, which form a seven transmembrane (TM) helices bundle with alternating extra-cellular and intracellular loops. GPCRs are considered to be one of the most important groups of drug targets because they are involved in a broad range of body functions and processes and are related to major diseases. In this paper we present a new technology, named PREDICT, for modeling the 3D structure of any GPCR from its amino acid sequence. This approach takes into account both internal protein properties (i.e., the amino acid sequence) and the properties of the membrane environment. Unlike competing approaches, the new technology does not rely on the single known structure of rhodopsin, and is thus capable of predicting novel GPCR conformations. We demonstrate the capabilities of PREDICT in reproducing the known experimental structure of rhodopsin. In principle, PREDICT-generated models offer new opportunities for structure-based drug discovery towards GPCR targets.

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.

Receptors coupling to G proteins: is there a signal behind the sequence?

Upon the binding of their ligands, G protein-coupled receptors couple to the heterotrimeric G proteins to transduce a signal. One receptor family may couple to a single G protein subtype and another family to several ones. Is there a signal in the receptor sequence that can give an indication of the G protein subtype selectivity? We used a sequence analysis method on biogenic amine and adenosine receptors and concluded that a weak signal can be detected in receptor families where specialization for coupling to a given G protein occurred during a recent divergent evolutionary process. Proteins 2000;41:448-459.

Alpha 1-adrenergic receptor regulation: basic science and clinical implications

Adrenergic receptors (ARs) are members of the G-protein-coupled receptor family, which includes alpha 1ARs, alpha 2ARs, beta 1ARs, beta 2ARs, beta 3ARs, adenosine, muscarinic, angiotensin, endothelin receptors, and many others that are responsible for a large variety of physiologic effects through G-protein coupling. This review focuses on alpha 1ARs and their regulation at both the mRNA and protein levels. Currently, three alpha 1AR subtypes have been characterized both pharmacologically and at the gene level: alpha 1aAR, alpha 1bAR, and alpha 1dAR. These are expressed in a species- and tissue-dependent manner. Mutagenesis approaches have been extremely valuable in the identification of key residues that govern alpha 1AR ligand binding and signaling. These studies reveal that alpha 1ARs have evolved an exquisitely sensitive regulation of their activity in which any disruption of the native structure has profound effects on subsequent function and effector coupling. Significant advances have also been made in the elucidation of signaling pathway components, resulting in the identification of novel pathways that can lead to pathologic conditions. Specific topics include mitogen-activated protein kinase, phosphatidylinositol 3-kinase, and G-protein-coupled receptor cross-talk pathways. Within this context, recent studies identifying underlying transcriptional mechanisms involved in the regulation of the alpha 1AR subtypes are also discussed. Finally, given the potentially important role of alpha 1ARs in the vasculature, as well as in the pathology of many diseases, such as myocardial hypertrophy and benign prostatic hyperplasia, the clinical relevance of alpha 1AR distribution, pharmacology, and therapeutic intervention is reviewed.

Dopamine receptors: from structure to behavior

The responses obtained with drugs that act at dopamine receptors depend on the spectrum of receptors stimulated, the pattern of stimulation and the neuronal signal-transduction pathways that are activated. In the absence of drugs that reliably discriminate between the various cloned receptors, elucidating the role of these receptors has largely relied on molecular genetic approaches that include expression of genes for receptors in cell lines and manipulation of this expression in animal models. Connecting molecular events that occur consequent to receptor stimulation with the resulting physiological effects entails bridging a complex network of interactions. This article reviews the current understanding of the molecular, cellular and systemic consequences of activation of the different dopamine receptors.

Molecular mechanisms underlying differential odor responses of a mouse olfactory receptor

The prevailing paradigm for G protein-coupled receptors is that each receptor is narrowly tuned to its ligand and closely related agonists. An outstanding problem is whether this paradigm applies to olfactory receptor (ORs), which is the largest gene family in the genome, in which each of 1,000 different G protein-coupled receptors is believed to interact with a range of different odor molecules from the many thousands that comprise "odor space." Insights into how these interactions occur are essential for understanding the sense of smell. Key questions are: (i) Is there a binding pocket? (ii) Which amino acid residues in the binding pocket contribute to peak affinities? (iii) How do affinities change with changes in agonist structure? To approach these questions, we have combined single-cell PCR results [Malnic, B., Hirono, J., Sato, T. & Buck, L. B. (1999) Cell 96, 713-723] and well-established molecular dynamics methods to model the structure of a specific OR (OR S25) and its interactions with 24 odor compounds. This receptor structure not only points to a likely odor-binding site but also independently predicts the two compounds that experimentally best activate OR S25. The results provide a mechanistic model for olfactory transduction at the molecular level and show how the basic G protein-coupled receptor template is adapted for encoding the enormous odor space. This combined approach can significantly enhance the identification of ligands for the many members of the OR family and also may shed light on other protein families that exhibit broad specificities, such as chemokine receptors and P450 oxidases.

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.

Encephalopsin: a novel mammalian extraretinal opsin discretely localized in the brain

We have identified a mammalian opsin, encephalopsin, that shows strong and specific expression in the brain. Encephalopsin defines a new family of opsins and shows highest homology to vertebrate retinal and pineal opsins. Encephalopsin is highly expressed in the preoptic area and paraventricular nucleus of the hypothalamus, both regions implicated in encephalic photoreception in nonmammalian vertebrates. In addition, encephalopsin shows highly patterned expression in other regions of the brain, being enriched in selected regions of the cerebral cortex, cerebellar Purkinje cells, a subset of striatal neurons, selected thalamic nuclei, and a subset of interneurons in the ventral horn of the spinal cord. Rostrocaudal gradients of encephalopsin expression are present in the cortex, cerebellum, and striatum. Radial stripes of encephalopsin expression are seen in the cerebellum. In the cortex and cerebellum, encephalopsin expression is considerably higher and more highly patterned in the adult than in the neonate. Encephalopsin is the first putative extraocular opsin identified in mammals and may play a role in encephalic photoreception.