Molecular interactions between the photoreceptor G protein and rhodopsin

1. The visual transduction system of the vertebrate retina is a well-studied model for biochemical and molecular studies of signal transduction. The structure and function of rhodopsin, a prototypical G protein-coupled receptor, and transducin or Gt, the photoreceptor G protein, have been particularly well studied. Mechanisms of rhodopsin-Gt interaction are discussed in this review. 2. The visual pigment rhodopsin contains a chromophore, and thus conformational changes leading to activation can be monitored spectroscopically. A model of the conformational changes in the activated receptor is presented based on biophysical and biochemical data. 3. The current information on sites of interaction on receptors and cognate G proteins is summarized. Studies using synthetic peptides from amino acid sequences corresponding to Gt and rhodopsin have provided information on the sites of rhodopsin-Gt interaction. Synthetic peptides from the carboxyl terminal region of alpha t mimic Gt by stabilizing the active conformation of rhodopsin, Metarhodopsin II. 4. The conformation of one such peptide when it is bound to Metarhodopsin II was determined by 2D NMR. The model based on the NMR data was tested using peptide analogs predicted to stabilize or break the structure. These studies yield molecular insight into why toxin-treated and mutant G proteins are uncoupled from receptors.     

Identification of CC chemokine receptor 7 residues important for receptor activation

The binding pocket of family A GPCRs that bind small biogenic amines is well characterized. In this study we identify residues on CC chemokine receptor 7 (CCR-7) that are involved in agonist-mediated receptor activation but not in high affinity ligand binding. The mutations also affect the ability of the ligands to induce chemotaxis. Two of the residues, Lys3.33(137) and Gln5.42(227), are consistent with the binding pocket described for biogenic amines, while Lys3.26(130) and Asn7.32(305), are found at, or close to, the cell surface. Our observations are in agreement with findings from other peptide and chemokine receptors, which indicate that receptors that bind larger ligands contain contact sites closer to the cell surface in addition to the conventional transmembrane binding pocket. These findings also support the theory that chemokine receptors require different sets of interactions for high affinity ligand binding and receptor activation.     

Lipid-facing correlated mutations and dimerization in G-protein coupled receptors.

A correlated mutation analysis has been performed on the aligned protein sequences of a number of class A G-protein coupled receptor families, including the chemokine, neurokinin, opioid, somatostatin, thyrotrophin and the whole biogenic amine family. Many of the correlated mutations are observed flanking or neighbouring conserved residues. The correlated residues have been plotted onto the transmembrane portion of the rhodopsin crystal structure. The structure shows that a significant proportion of the correlated mutations are located on the external (lipid-facing) region of the helices. The occurrence of these highly correlated patterns of change amongst the external residues suggest that they are sites for protein-protein interactions. In particular, it is suggested that the correlated residues may be involved in either large conformational changes, the formation of heterodimers or homodimers (which may be domain swapped) or oligomers required for activation or internalization. The results are discussed in the light of the subtype-specific heterodimerization observed for the chemokine, opioid and somatostatin receptors.     

Structural determinants for selective recognition of peptide ligands for endothelin receptor subtypes ETA and ETB

The molecular basis for recognition of peptide ligands endothelin‐1, ‐2 and ‐3 in endothelin receptors is poorly understood. Especially the origin of ligand selectivity for ET_A or ET_B is not clearly resolved. We derived sequence‐structure‐function relationships of peptides and receptors from mutational data and homology modeling. Our major findings are the dissection of peptide ligands into four epitopes and the delineation of four complementary structural portions on receptor side explaining ligand recognition in both endothelin receptor subtypes. In addition, structural determinants for ligand selectivity could be described. As a result, we could improve the selectivity of BQ3020 about 10‐fold by a single amino acid substitution, validating our hypothesis for ligand selectivity caused by different entrances to the receptors' transmembrane binding sites. A narrow tunnel shape in ET_A is restrictive for a selected group of peptide ligands' N‐termini, whereas a broad funnel‐shaped entrance in ET_B accepts a variety of different shapes and properties of ligands. Copyright © 2009 European Peptide Society and John Wiley & Sons, Ltd.     

Photoaffinity scanning in the mapping of the peptide receptor interface of class II G protein-coupled receptors.

The family of G protein-coupled receptors constitutes about 50% of the therapeutic drug targets used in clinical medicine today, although the mechanisms of ligand binding, activation and signal transduction for G protein-coupled receptors are not yet well defined. This review discusses ongoing research using the photoaffinity scanning method to map the bimolecular interface between class II G protein-coupled receptors and their ligands. Furthermore the available computer model of class II peptide ligand docking into the receptor, based on the positional constraints imposed by the photoaffinity scanning analyses, will be discussed briefly. The ultimate goal of these efforts is to understand the molecular basis of receptor binding and therefore to generate a template for rational drug design.     

Sequence,Structure and Ligand Binding Evolution of Rhodopsin-Like G Protein-Coupled Receptors: A Crystal Structure-Based Phylogenetic Analysis

G protein-coupled receptors (GPCRs) form the largest family of membrane receptors in the human genome. Advances in membrane protein crystallization so far resulted in the determination of 24 receptors available as high-resolution atomic structures. We performed the first phylogenetic analysis of GPCRs based on the available set of GPCR structures. We present a new phylogenetic tree of known human rhodopsin-like GPCR sequences based on this structure set. We can distinguish the three separate classes of small-ligand binding GPCRs, peptide binding GPCRs, and olfactory receptors. Analyzing different structural subdomains, we found that small molecule binding receptors most likely have evolved from peptide receptor precursors, with a rhodopsin/S1PR1 ancestor, most likely an ancestral opsin, forming the link between both classes. A light-activated receptor therefore seems to be the origin of the small molecule hormone receptors of the central nervous system. We find hints for a common evolutionary path of both ligand binding site and central sodium/water binding site. Surprisingly, opioid receptors exhibit both a binding cavity and a central sodium/water binding site similar to the one of biogenic amine receptors instead of peptide receptors, making them seemingly prone to bind small molecule ligands, e.g. opiates. Our results give new insights into the relationship and the pharmacological properties of rhodopsin-like GPCRs.     

Ligand binding properties of human galanin receptors

Abstract

The galanin receptor family comprises of three members, GalR1, GalR2 and GalR3, all belonging to the G-protein-couple receptor superfamily. All three receptors bind the peptide hormone galanin, but show distinctly different binding properties to other molecules and effects on intracellular signaling. To gain insight on the molecular basis of receptor subtype specificity, we have generated a three-dimensional model for each of the galanin receptors based on its homologs in the same family. We found significant differences in the organization of the binding pockets among the three types of receptors, which might be the key for specific molecular recognition of ligands. Through docking of fragments of the galanin peptide and a number of ligands, we investigated the involvement of transmembrane and loop residues in ligand interaction.     

The unique ligand-binding pocket for the human prostacyclin receptor. Site-directed mutagenesis and molecular modeling

The human prostacyclin receptor is a seven-transmembrane alpha-helical G-protein coupled receptor, which plays important roles in both vascular smooth muscle relaxation as well as prevention of blood coagulation. The position of the native ligand-binding pocket for prostacyclin as well as other derivatives of the 20-carbon eicosanoid, arachidonic acid, has yet to be determined. Through the use of prostanoid receptor sequence alignments, site-directed mutagenesis, and the 2.8-A x-ray crystallographic structure of bovine rhodopsin, we have developed a three-dimensional model of the agonist-binding pocket within the seven-transmembrane (TM) domains of the human prostacyclin receptor. Upon mutation to alanine, 11 of 29 candidate residues within TM domains II, III, IV, V, and VII exhibited a marked decrease in agonist binding. Of this group, four amino acids, Arg-279 (TMVII), Phe-278 (TMVII), Tyr-75 (TMII), and Phe-95 (TMIII), were identified (via receptor amino acid sequence alignment, ligand structural comparison, and computer-assisted homology modeling) as having direct molecular interactions with ligand side-chain constituents. This binding pocket is distinct from that of the biogenic amine receptors and rhodopsin where the native ligands (also composed of a carbon ring and a carbon chain) are accommodated in an opposing direction. These findings should assist in the development of novel and highly specific ligands including selective antagonists for further molecular pharmacogenetic studies of the human prostacyclin receptor.     

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

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

A computer modeling postulated mechanism for angiotensin II receptor activation

The angiotensin II receptor of the AT₁-type has been modeled starting from the experimentally determined three-dimensional structure of bacteriorhodopsin as the template. Intermediate 3D structures of rhodopsin and ₂-adrenergic receptors were built because no direct sequence alignment is possible between the AT₁ receptor and bacteriorhodopsin. Docking calculations were carried out on the complex of the modeled receptor with AII, and the results were used to analyze the binding possibilities of DuP753-type antagonistic non-peptide ligands. We confirm that the positively charged Lys¹⁹⁹ on helix 5 is crucial for ligand binding, as in our model; the charged side chain of this amino acid interacts strongly with the C-terminal carboxyl group of peptide agonists or with the acidic group at the 2-position of the biphenyl moiety of DuP753-type antagonists. Several other receptor residues which are implicated in the binding of ligands and the activation of receptor by agonists are identified, and their functional role is discussed. Therefore, a plausible mechanism of receptor activation is proposed. The three-dimensional docking model integrates most of the available experimental observations and helps to plan pertinent site-directed mutagenesis experiments which in turn may validate or modify the present model and the proposed mechanism of receptor activation.     

Structure and regulation of opioid receptors

Significant advances have been made in understanding the structure, function, and regulation of opioid receptors and endogenous opioid peptides since their discovery approximately 25 years ago. This review summarizes recent studies aimed at identifying key amino acids that confer ligand selectivity to the opioid receptors and that are critical constituents of the ligand binding sites. A molecular model of the delta receptor based on the crystal structure of rhodopsin is presented. Agonist-induced down regulation of opioid receptors is discussed, highlighting recent evidence for the involvement of the ubiquitin/proteasome system in this process.     

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

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

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

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

Characterization of 5-HT transporter and receptor system in HeLaS3 cells by [(3)H]8-OH-DPAT and other serotonergic ligands

[(3)H]8-OH-DPAT is a selective ligand for labeling 5-HT(1A) receptor sites. In competition binding experiments, we found that classic biogenic amine transporter inhibitors displaced [(3)H]8-OH-DPAT binding at its high-affinity binding sites in HeLaS3 cells. [(125)I]RTI-55 and [(3)H]paroxetine are known to specifically label amine transporter sites, and this was observed in our cells. Displacement studies showed that 8-OH-DPAT displayed affinity in a dose-dependent manner for the labeled amine transporter sites. These data suggest that [(3)H]8-OH-DPAT binds to amine uptake sites in HeLaS3 cells. A variety of drugs targeting different classes of receptors did not significantly affect [(3)H]8-OH-DPAT binding. Moreover, we determined the specific binding effects of various serotonergic ligands (i.e. [(125)I]cyanopindolol, [(3)H]ketanserin/[(3)H]mesulergine, [(3)H]GR-65630, [(3)H]GR-113808 and [(3)H]LSD) that specifically labeled 5-HT(1), 5-HT(2), 5-HT(3), 5-HT(4) and 5-HT(5-7) receptors, respectively. It is suggested that HeLaS3 cells contain distinct types of the related to 5-HT receptor recognition binding sites. These observations could help elucidate the relevant characteristics of different types of 5-HT receptors and 5-HT membrane transporters in tumor cells and their role in tumorigenesis.     

Multiple molecular recognition properties of the lipocalin protein family

The lipocalins, a diverse family of small extracellular ligand proteins, display a remarkable range of different molecular properties. While their binding of small hydrophobic molecules, and to a lesser extent their binding to cell surface receptors, is well known, it is shown here that formation of macromolecular complexes is also a common feature of this family. Analysis of known crystallographic structures reveals that the lipocalins process a conserved common structure: an antiparallel β-barrel with a repeated +1 topology. Comparisons show that within this overall similarity the structure of individual proteins is specifically adapted to bind their particular ligands, forming a binding site from an internal cavity (within the barrel) and/or an external loop scaffold, which gives rise to different binding modes that reflects the need to accommodate ligands of different shape, size, and chemical structure. The architecture of the lipocalin fold suggests that the both the ends and sides of this barrel are topologically distinct, differences also apparent in analyses of structural and sequence variation within the family. These different can be linked to experimental evidence suggesting a possible functional dichotomy between the two ends of the lipocalin fold. The structurally invariant end of the molecule may be implicated in general binding small ligands and forming macromolecular complexes via an exposed binding surface.     

Construction of a sequence motif characteristic of aminergic G protein-coupled receptors

An approach to discover sequence patterns characteristic of ligand classes is described and applied to aminergic G protein-coupled receptors (GPCRs). Putative ligand-binding residue positions were inferred from considering three lines of evidence: conservation in the subfamily absent or underrepresented in the superfamily, any available mutation data, and the physicochemical properties of the ligand. For aminergic GPCRs, the motif is composed of a conserved aspartic acid in the third transmembrane (TM) domain (rhodopsin position 117) and a conserved tryptophan in the seventh TM domain (rhodopsin position 293); the roles of each are readily justified by molecular modeling of ligand-receptor interactions. This minimally defined motif is an appropriate computational tool for identifying additional, potentially novel aminergic GPCRs from a set of experimentally uncharacterized "orphan" GPCRs, complementing existing sequence matching, clustering, and machine-learning techniques. Motif sensitivity stems from the stepwise addition of residues characteristic of an entire class of ligand (and not tailored for any particular biogenic amine). This sensitivity is balanced by careful consideration of residues (evidence drawn from mutation data, correlation of ligand properties to residue properties, and location with respect to the extracellular face), thereby maintaining specificity for the aminergic class. A number of orphan GPCRs assigned to the aminergic class by this motif were later discovered to be a novel subfamily of trace amine GPCRs, as well as the successful classification of the histamine H4 receptor.     

Characterization of 5-HT transporter and receptor system in HeLaS3 cells by [H]8-OH-DPAT and other serotonergic ligands

[³H]8-OH-DPAT is a selective ligand for labeling 5-HT_1A receptor sites. In competition binding experiments, we found that classic biogenic amine transporter inhibitors displaced [³H]8-OH-DPAT binding at its high-affinity binding sites in HeLaS3 cells. [¹²⁵I]RTI-55 and [³H]paroxetine are known to specifically label amine transporter sites, and this was observed in our cells. Displacement studies showed that 8-OH-DPAT displayed affinity in a dose-dependent manner for the labeled amine transporter sites. These data suggest that [³H]8-OH-DPAT binds to amine uptake sites in HeLaS3 cells. A variety of drugs targeting different classes of receptors did not significantly affect [³H]8-OH-DPAT binding. Moreover, we determined the specific binding effects of various serotonergic ligands (i.e. [¹²⁵I]cyanopindolol, [³H]ketanserin/[³H]mesulergine, [³H]GR-65630, [³H]GR-113808 and [³H]LSD) that specifically labeled 5-HT₁, 5-HT₂, 5-HT₃, 5-HT₄ and 5-HT_5–7 receptors, respectively. It is suggested that HeLaS3 cells contain distinct types of the related to 5-HT receptor recognition binding sites. These observations could help elucidate the relevant characteristics of different types of 5-HT receptors and 5-HT membrane transporters in tumor cells and their role in tumorigenesis.     

mRNA display selection and solid‐phase synthesis of Fc‐binding cyclic peptide affinity ligands

Cyclic peptides are attractive candidates for synthetic affinity ligands due to their favorable properties, such as resistance to proteolysis, and higher affinity and specificity relative to linear peptides. Here we describe the discovery, synthesis and characterization of novel cyclic peptide affinity ligands that bind the Fc portion of human Immunoglobulin G (IgG; hFc). We generated an mRNA display library of cyclic pentapeptides wherein peptide cyclization was achieved with high yield and selectivity, using a solid‐phase crosslinking reaction between two primary amine groups, mediated by a homobifunctional linker. Subsequently, a pool of cyclic peptide binders to hFc was isolated from this library and chromatographic resins incorporating the selected cyclic peptides were prepared by on‐resin solid‐phase peptide synthesis and cyclization. Significantly, this approach results in resins that are resistant to harsh basic conditions of column cleaning and regeneration. Further studies identified a specific cyclic peptide—cyclo[Link‐M‐WFRHY‐K]—as a robust affinity ligand for purification of IgG from complex mixtures. The cyclo[Link‐M‐WFRHY‐K] resin bound selectively to the Fc fragment of IgG, with no binding to the Fab fragment, and also bound immunoglobulins from a variety of mammalian species. Notably, while the recovery of IgG using the cyclo[Link‐M‐WFRHY‐K] resin was comparable to a Protein A resin, elution of IgG could be achieved under milder conditions (pH 4 vs. pH 2.5). Thus, cyclo[Link‐M‐WFRHY‐K] is an attractive candidate for developing a cost‐effective and robust chromatographic resin to purify monoclonal antibodies (mAbs). Finally, our approach can be extended to efficiently generate and evaluate cyclic peptide affinity ligands for other targets of interest. Biotechnol. Bioeng. 2013; 110: 857–870. © 2012 Wiley Periodicals, Inc.     

Multiple residues in the second extracellular loop are critical for M3 muscarinic acetylcholine receptor activation

Recent studies suggest that the second extracellular loop (o2 loop) of bovine rhodopsin and other class I G protein-coupled receptors (GPCRs) targeted by biogenic amine ligands folds deeply into the transmembrane receptor core where the binding of cis-retinal and biogenic amine ligands is known to occur. In the past, the potential role of the o2 loop in agonist-dependent activation of biogenic amine GPCRs has not been studied systematically. To address this issue, we used the M(3) muscarinic acetylcholine receptor (M3R), a prototypic class I GPCR, as a model system. Specifically, we subjected the o2 loop of the M3R to random mutagenesis and subsequently applied a novel yeast genetic screen to identity single amino acid substitutions that interfered with M3R function. This screen led to the recovery of about 20 mutant M3Rs containing single amino acid changes in the o2 loop that were inactive in yeast. In contrast, application of the same strategy to the extracellular N-terminal domain of the M3R did not yield any single point mutations that disrupted M3R function. Pharmacological characterization of many of the recovered mutant M3Rs in mammalian cells, complemented by site-directed mutagenesis studies, indicated that the presence of several o2 loop residues is important for efficient agonist-induced M3R activation. Besides the highly conserved Cys(220) residue, Gln(207), Gly(211), Arg(213), Gly(218), Ile(222), Phe(224), Leu(225), and Pro(228) were found to be of particular functional importance. In general, mutational modification of these residues had little effect on agonist binding affinities. Our findings are therefore consistent with a model in which multiple o2 loop residues are involved in stabilizing the active state of the M3R. Given the high degree of structural homology found among all biogenic amine GPCRs, our findings should be of considerable general relevance.     

Protection of all cleavable sites of DNA with the multiple CGCG or continuous CGG sites from the restriction enzyme,indicative of simultaneous binding of small ligands

The anthracenone ligands (1–12) with a keto-phenol and a hydroxamic acid unit were synthesized and evaluated by a restriction enzyme inhibition assay. DNA substrates composed of multiple CGCG or CGG sites are fully hydrolyzed by a restriction enzyme that is selective for each sequence. Under such conditions, the full-length DNA substrate remains only when the ligand binds to all binding sites and protects it from hydrolysis by the restriction enzymes. In the assay using AccII and the 50-mer DNA substrates containing a different number of CGCG sites at different non-binding AT base pair intervals, the more the CGCG sites, the more the full-length DNA increased. Namely, simultaneous binding of the ligand (5) to the CGCG sites increased in the order of (CGCG)₅>(CGCG)₂>(CGCG)₁. Furthermore, the length of the spacer of the hydroxamic acid to the anthracenone skeleton played an important role in the preference for the number of the d(A/T) base pairs between the CGCG sites. The long spacer-ligand (5) showed a preference to the CGCG sites with five AT pairs, and the short spacer-ligand (10) to that with two AT pairs. The ligand (12) with the shortest spacer showed a preference in simultaneous binding to the 54-mer DNA composed of 16 continuous CGG sites in the assay using the restriction enzyme Fnu4HI that hydrolyzes the d(GCGGC)/d(CGCCG) site. Application of these ligands to biological systems including the repeat DNA sequence should be of significant interest.