Biosynthesis and NMR-studies of a double transmembrane domain from the Y4 receptor,a human GPCR

The human Y4 receptor, a class A G-protein coupled receptor (GPCR) primarily targeted by the pancreatic polypeptide (PP), is involved in a large number of physiologically important functions. This paper investigates a Y4 receptor fragment (N-TM1-TM2) comprising the N-terminal domain, the first two transmembrane (TM) helices and the first extracellular loop followed by a (His)₆ tag, and addresses synthetic problems encountered when recombinantly producing such fragments from GPCRs in Escherichia coli. Rigorous purification and usage of the optimized detergent mixture 28 mM dodecylphosphocholine (DPC)/118 mM% 1-palmitoyl-2-hydroxy-sn-glycero-3-[phospho-rac-(1-glycerol)] (LPPG) resulted in high quality TROSY spectra indicating protein conformational homogeneity. Almost complete assignment of the backbone, including all TM residue resonances was obtained. Data on internal backbone dynamics revealed a high secondary structure content for N-TM1-TM2. Secondary chemical shifts and sequential amide proton nuclear Overhauser effects defined the TM helices. Interestingly, the properties of the N-terminal domain of this large fragment are highly similar to those determined on the isolated N-terminal domain in the presence of DPC micelles.     

GPCR structures in drug design,emerging opportunities with new structures

In recent years, GPCR targets from diverse regions of phylogenetic space have been determined. This effort has culminated this year in the determination of representatives of all major classes of GPCRs (A, B, C, and F). Although much of the now well established knowledge on GPCR structures has been known for some years, the new high-resolution structures allow structural insight into the causes of ligand efficacy, biased signaling, and allosteric modulation. In this digest the structural basis for GPCR signaling in the light of the new structures is reviewed and the use of the new non-class A GPCRs for drug design is discussed.     

Advances in GPCR Modeling Evaluated by the GPCR Dock 2013 Assessment: Meeting New Challenges

1. Download : Download high-res image (559KB)
2. Download : Download full-size image

     

A model for the evaluation of domain based classification of GPCR

G-Protein Coupled Receptors (GPCR) are the largest family of membrane bound receptor and plays a vital role in various biological processes with their amenability to drug intervention. They are the spotlight for the pharmaceutical industry. Experimental methods are both time consuming and expensive so there is need to develop a computational approach for classification to expedite the drug discovery process. In the present study domain based classification model has been developed by employing and evaluating various machine learning approaches like Bagging, J48, Bayes net, and Naive Bayes. Various softwares are available for predicting domains. The result and accuracy of output for the same input varies for these software''s. Thus, there is dilemma in choosing any one of it. To address this problem, a simulation model has been developed using well known five softwares for domain prediction to explore the best predicted result with maximum accuracy. The classifier is developed for classification up to 3 levels for class A. An accuracy of 98.59% by Naïve Bayes for level I, 92.07% by J48 for level II and 82.14% by Bagging for level III has been achieved.     

Insights from the analysis of conserved motifs and permitted amino acid exchanges in the human, the fly and the worm GPCR clusters

G-protein coupled receptors (GPCRs) belong to biologically important and functionally diverse and largest super family of membrane proteins. GPCRs retain a characteristic membrane topology of seven alpha helices with three intracellular, three extracellular loops and flanking N' and C' terminal residues. Subtle differences do exist in the helix boundaries (TM-domain), loop lengths, sequence features such as conserved motifs, and substituting amino acid patterns and their physiochemical properties amongst these sequences (clusters) at intra-genomic and inter-genomic level (please re-phrase into 2 statements for clarity). In the current study, we employ prediction of helix boundaries and scores derived from amino acid substitution exchange matrices to identify the conserved amino acid residues (motifs) as consensus in aligned set of homologous GPCR sequences. Co-clustered GPCRs from human and other genomes, organized as 32 clusters, were employed to study the amino acid conservation patterns and species-specific or cluster-specific motifs. Critical analysis on sequence composition and properties provide clues to connect functional relevance within and across genome for vast practical applications such as design of mutations and understanding of disease-causing genetic abnormalities.     

Subunit-selective role of the M3 transmembrane domain of the nicotinic acetylcholine receptor in channel gating

The nicotinic acetylcholine receptor (AChR) can be either hetero-pentameric, composed of α and non-α subunits, or homo-pentameric, composed of α7 subunits. To explore the subunit-selective contributions of transmembrane domains to channel gating we analyzed single-channel activity of chimeric muscle AChRs. We exchanged M3 between α1 and ? or α7 subunits. The replacement of M3 in α1 by ?M3 significantly alters activation properties. Channel activity appears as bursts of openings whose durations are 20-fold longer than those of wild-type AChRs. In contrast, 7-fold briefer openings are observed in AChRs containing the reverse ? chimeric subunit. The duration of the open state decreases with the increase in the number of α1M3 segments, indicating additive contributions of M3 of all subunits to channel closing. Each α1M3 segment decreases the energy barrier of the closing process by ∼ 0.8 kcal/mol. Partial chimeric subunits show that small stretches of the M3 segment contribute additively to the open duration. The replacement of α1 sequence by α7 in M3 leads to 3-fold briefer openings whereas in M1 it leads to 10-fold prolonged openings, revealing that the subunit-selective role is unique to each transmembrane segment.     

Strong dimerization of wild-type ErbB2/Neu transmembrane domain and the oncogenic Val664Glu mutant in mammalian plasma membranes

Here, we study the homodimerization of the transmembrane domain of Neu, as well as an oncogenic mutant (V664E), in vesicles derived from the plasma membrane of mammalian cells. For the characterization, we use a Förster resonance energy transfer (FRET)-based method termed Quantitative Imaging-FRET (QI-FRET), which yields the donor and acceptor concentrations in addition to the FRET efficiencies in individual plasma membrane-derived vesicles. Our results demonstrate that both the wild-type and the mutant are 100% dimeric, suggesting that the Neu TM helix dimerizes more efficiently than other RTK TM domains in mammalian membranes. Furthermore, the data suggest that the V664E mutation causes a very small, but statistically significant change in dimer structure. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

Differences in Binding Sites of Two Melatonin Receptors Help to Explain Their Selectivity to Some Melatonin Analogs: A Molecular Modeling Study

Abstract

Numerous diseases have been linked to the malfunction of G-protein coupled receptors (GP-CRs). Their adequate treatment requires rational design of new high-affinity and high-selectivity drugs targeting these receptors. In this work, we report three-dimensional models of the human MT₁ and MT₂ melatonin receptors, members of the GPCR family. The models are based on the X-ray structure of bovine rhodopsin. The computational approach employs an original procedure for optimization of receptor-ligand structures. It includes rotation of one of the transmembrane α-helices around its axis with simultaneous assessment of quality of the resulting complexes according to a number of criteria we have developed for this purpose. The optimal geometry of the receptor-ligand binding is selected based on the analysis of complementarity of hydrophobic/hydrophilic properties between the ligand and its protein environment in the binding site. The elaborated “optimized” models are employed to explore the details of protein-ligand interactions for melatonin and a number of its analogs with known affinity to MT₁ and MT₂ receptors. The models permit rationalization of experimental data, including those that were not used in model building. The perspectives opened by the constructed models and by the optimization procedure in the design of new drugs are discussed.     

Molecular dynamics studies of the transmembrane domain of gp41 from HIV-1

Helix-helix interactions in the putative three-helix bundle formation of the gp41 transmembrane (TM) domain may contribute to the process of virus-cell membrane fusion in HIV-1 infection. In this study, molecular dynamics is used to analyze and compare the conformations of monomeric and trimeric forms of the TM domain in various solvent systems over the course of 4 to 23-ns simulations. The trimeric bundles of the TM domain were stable as helices and remained associated in a hydrated POPE lipid bilayer for the duration of the 23-ns simulation. Several stable inter-chain hydrogen bonds, mostly among the three deprotonated arginine residues located at the center of each of the three TM domains, formed in a right-handed bundle embedded in the lipid bilayer. No such bonds were observed when the bundle was left-handed or when the central arginine residue in each of the three TM helices was replaced with isoleucine (R_I mutant), suggesting that the central arginine residues may play an essential role in maintaining the integrity of the three-helix bundle. These observations suggest that formation of the three-helix bundle of the TM domain may play a role in the trimerization of gp41, thought to occur during the virus-cell membrane fusion process.     

Focused library design in GPCR projects on the example of 5-HT(2c) agonists: comparison of structure-based virtual screening with ligand-based search methods

The aim of this study was to investigate the usefulness of structure-based virtual screening (VS) for focused library design in G protein-coupled receptors (GPCR) projects on the example of 5-HT(2c) agonists. We compared the performance of structure-based VS against two different homology models using FRED for docking and ScreenScore, FlexX, and PMF for rescoring with the results of 12 ligand-based similarity searches using four different query compounds and three different similarity metrics (Daylight, FTree, Phacir). The result of the similarity search showed much variation, from an enrichment factor up to 3.2 to worse than random, whereas the structure-based VS gave a more stable result with a constant enrichment factor around 2. Additionally, actives retrieved by the structure-based approach were more diverse than the actives among the top scorers of the similarity searches. Based on these results, we suggest basing a focused library design for a GPCR project on a combination of a ligand-based similarity search and structure-based docking.     

TM-MOTIF: an alignment viewer to annotate predicted transmembrane helices and conserved motifs in aligned set of sequences

Multiple sequence alignments become biologically meaningful only if conserved and functionally important residues and secondary structural elements preserved can be identified at equivalent positions. This is particularly important for transmembrane proteins like G-protein coupled receptors (GPCRs) with seven transmembrane helices. TM-MOTIF is a software package and an effective alignment viewer to identify and display conserved motifs and amino acid substitutions (AAS) at each position of the aligned set of homologous sequences of GPCRs. The key feature of the package is to display the predicted membrane topology for seven transmembrane helices in seven colours (VIBGYOR colouring scheme) and to map the identified motifs on its respective helices /loop regions. It is an interactive package which provides options to the user to submit query or pre-aligned set of GPCR sequences to align with a reference sequence, like rhodopsin, whose structure has been solved experimentally. It also provides the possibility to identify the nearest homologue from the available inbuilt GPCR or Olfactory Receptor cluster dataset whose association is already known for its receptor type. AVAILABILITY: The database is available for free at mini@ncbs.res.in.     

The 3D structure of the defense-related rice protein Pir7b predicted by homology modeling and ligand binding studies

To better understand the ligand-binding mechanism of protein Pir7b, important part in detoxification of a pathogen-derived compound against Pyricularia oryzae, a 3D structure model of protein Pir7b was constructed based on the structure of the template SABP2. Three substrates were docking to this protein, two of them were proved to be active, and some critical residues are identified, which had not been confirmed by the experiments. His87 and Leu17 considered as ‘oxyanion hole’ contribute to initiating the Ser86 nucleophilic attack. Gln187 and Asp139 can form hydrogen bonds with the anilid group to maintain the active binding orientation with the substrates. The docking model can well interpret the specificity of protein Pir7b towards the anilid moiety of the substrates and provide valuable structure information about the ligand binding to protein Pir7b. Figure Ligand binding analysis based on the refined Pir7b model. Magenta dash line, hydrogen bond; Red dash line, distance label. (a) Docking of 2-naphthol AS-acetate to Pir7b model. A 3D figure of 2-naphthol AS-acetate-Pir7b complex is also attached (b) Docking of 2-naphthol AS-2-chlor-propionate to Pir7b model. (c) Docking of 2-naphthol-acetate to Pir7b model.     

Oligomerization of the fifth transmembrane domain from the adenosine A2A receptor

The human adenosine A_2A receptor (A_2AR) belongs to one of the largest family of membrane proteins, the G-protein coupled receptors (GPCRs), characterized by seven transmembrane (TM) helices. Little is known about the determinants of their structures, folding, assembly, activation mechanisms, and oligomeric states. Previous studies in our group showed that peptides corresponding to all seven TM domains form stable helical structures in detergent micelles and lipid vesicles. However, the peptides behave differently; TM5 is the only peptide to have a ratio [θ]₂₂₂/[θ]₂₀₈ obtained by circular dichroism (CD) spectroscopy>1. This finding suggested to us that TM5 might self-associate. In the present study, we investigate the unique properties of the TM5 domain. We performed detailed analyses of TM5 peptide behavior in membrane-mimetic environments using CD spectroscopy, fluorescence spectroscopy and Förster resonance energy transfer, and gel electrophoresis. We find that TM5 peptide has the ability to self-associate to form oligomeric structures in various hydrophobic milieus and that these oligomers are highly resistant to temperature and chemical denaturation. We also find that mutation of the full-length A_2AR at position M193, which is located in the fifth TM domain, noticeably alters A_2AR monomer: dimer ratio as observed on SDS-PAGE. Our results suggest that parallel association of TM5 dimers may play a role in the known adenosine A_2A receptor dimerization. This study represents the first evidence of an individual GPCR transmembrane domain self-association.     

Monte Carlo loop refinement and virtual screening of the thyroid-stimulating hormone receptor transmembrane domain

Metropolis Monte Carlo (MMC) loop refinement has been performed on the three extracellular loops (ECLs) of rhodopsin and opsin-based homology models of the thyroid-stimulating hormone receptor transmembrane domain, a class A type G protein-coupled receptor. The Monte Carlo sampling technique, employing torsion angles of amino acid side chains and local moves for the six consecutive backbone torsion angles, has previously reproduced the conformation of several loops with known crystal structures with accuracy consistently less than 2?Å. A grid-based potential map, which includes van der Waals, electrostatics, hydrophobic as well as hydrogen-bond potentials for bulk protein environment and the solvation effect, has been used to significantly reduce the computational cost of energy evaluation. A modified sigmoidal distance-dependent dielectric function has been implemented in conjunction with the desolvation and hydrogen-bonding terms. A long high-temperature simulation with 2?kcal/mol repulsion potential resulted in extensive sampling of the conformational space. The slow annealing leading to the low-energy structures predicted secondary structure by the MMC technique. Molecular docking with the reported agonist reproduced the binding site within 1.5?Å. Virtual screening performed on the three lowest structures showed that the ligand-binding mode in the inter-helical region is dependent on the ECL conformations.     

Structure prediction of the dimeric neu/ErbB-2 transmembrane domain from multi-nanosecond molecular dynamics simulations

Dimerization of the neu/ErbB-2 receptor tyrosine kinase is a necessary but not a sufficient step for signaling. Despite the efforts expended to identify the molecular interactions responsible for receptor-receptor contacts and particularly those involving the transmembrane domain, structural details are still unknown. In this work, molecular dynamics simulations of the helical transmembrane domain (TM) of neu and ErbB-2 receptors are used to predict their dimer structure both in the wild and oncogenic forms. A global conformational search method, applied to define the best orientations of parallel helices, showed an energetically favorable configuration with the specific mutation site within the interface, common for both the nontransforming and the transforming neu/ErbB-2 TM dimers. Starting from this configuration, a total of 10 simulations, about 1.4 ns each, performed in vacuum, without any constraints, show that the two helices preferentially wrap in left-handed interactions with a packing angle at about 20°. The resulting structures are nonsymmetric and the hydrogen bond network analysis shows that helices experience π local distortions that facilitate inter-helix hydrogen bond interactions and may result in a change in the helix packing, leading to a symmetric interface. For the mutated sequences, we show that the Glu side chain interacts directly with its cognate or with carbonyl groups of the facing backbone. We show that the connectivity between interfacial residues conforms to the knobs-into-holes packing mode of transmembrane helices. The dimeric interface described in our models is discussed with respect to mutagenesis studies. Received: 12 March 1999 / Revised version: 23 August 1999 / Accepted: 23 August 1999     

Proline residues link the active site to transmembrane domain movements in human nucleoside triphosphate diphosphohydrolase 3 (NTPDase3)

The active sites of the membrane-bound nucleoside triphosphate diphosphohydrolases (NTPDases) regulate and are regulated by coordinated and spatially distant movements of their transmembrane helices, modulating enzyme activity, and substrate specificity. Using site-directed mutagenesis, the roles of the conserved proline residues (N-terminal: P52 and P53; C-terminal: P472, P476, P481, P484, and P485) of human NTPDase3, located in the “linker regions” that connect the N- and C-terminal transmembrane helices with the extracellular active site, were examined. Single cysteine substitutions were strategically placed in the transmembrane domain (N-terminal helix: V42C; C-terminal helix: G489C) to serve as cross-linking “sensors” of helical interactions. These “sensor” background mutant proteins (V42C and G489C NTPDase3) are enzymatically active and are cross-linked by copper phenanthroline less efficiently in the presence of adenosine triphosphate (ATP). Proline to alanine substitutions at P53, P481, P484, and P485 in the V42C background, as well as P53, P481, and P484 in the G489C background, exhibited decreased nucleotidase activities. More importantly, alanine substitutions at P53 and P481 in the V42C background and P481 in the G489C background no longer exhibited the ATP-induced decrease in transmembrane cross-linking efficiency. Interestingly, the P485A mutation abolished oxidative cross-linking at G489C both in the presence and absence of ATP. Taken together, these results suggest a role for proline residues 53 and 481 in the linker regions of human NTPDase3 for coupling nucleotide binding at the enzyme active site to movements and/or rearrangements of the transmembrane helices necessary for optimal nucleotide hydrolysis.     

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

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