GRIP: A server for predicting interfaces for GPCR oligomerization

G-Protein Coupled Receptors (GPCRs) are one of the most important pharmaceutical targets. Recent studies have revealed that many GPCRs form homo- and/or hetero-oligomers. The molecular mechanisms of oligomerization are not fully understood yet, due to the lack of structural data for GPCR complexes. Therefore, accurate interface prediction would accelerate investigations of the molecular mechanisms of oligomerization and signaling via GPCRs. However, interface prediction for GPCR oligomerization is difficult, because the various GPCR subtypes often use different structural regions as their interfaces, even when the subtypes belong to the same subfamily. Previously, we developed a method to predict the interfaces for GPCR oligomerization, which overcomes the difficulty described above. We have now launched a web service, named G-protein coupled Receptors Interaction Partners (GRIP) (http://grip.cbrc.jp/GRIP/index.html), to predict the interfaces for GPCR oligomerization. As far as we know, it is the only service to predict the interfaces for GPCR oligomerization.     

Prediction of interfaces for oligomerizations of G-protein coupled receptors

Several lines of biochemical and pharmacological evidence have suggested that some G-protein-coupled receptors (GPCRs) form homo oligomers, hetero oligomers or both. The GPCRs oligomerizations are considered to be related to signal transduction and some diseases. Therefore, an accurate prediction of the residues that interact upon oligomerization interface would further our understanding of signal transduction and the diseases in which GPCRs are involved. One of the complications for such a prediction is that the interfaces differ with the subtypes, even within the same GPCR family. Focusing on the distribution of residues conserved on the molecular surface in a particular subtype, we developed a new method to predict the interface for the GPCR oligomers, and applied it to several subtypes of known GPCRs to check the sensitivity. Subsequently, we found that predicted interfaces of rhodopsin, D(2) dopamine receptor and beta(2) adrenergic receptor agreed with the experimentally suggested interfaces, despite difference in the interface region among the three subtypes. Moreover, a highly conserved residue detected from the D(2) dopamine receptor corresponded to a residue involved in a missense change found in the large family of myoclonus dystonia. Our observation suggests the possibility that the disease is caused by the disorder of the oligomerization, although the molecular mechanism of the disease has not been revealed yet. The benefits and the pitfalls of the new method will be discussed, based on the results of the applications.     

Influence of oligomerization on the dynamics of G-protein coupled receptors as assessed by normal mode analysis

The recently discovered impact of oligomerization on G-protein coupled receptor (GPCR) function further complicates the already challenging goal of unraveling the molecular and dynamic mechanisms of these receptors. To help understand the effect of oligomerization on the dynamics of GPCRs, we have compared the motion of monomeric, dimeric, and tetrameric arrangements of the prototypic GPCR rhodopsin, using an approximate-yet powerful-normal mode analysis (NMA) technique termed elastic network model (ENM). Moreover, we have used ENM to discriminate between putative dynamic mechanisms likely to account for the recently observed conformational rearrangement of the TM4,5-TM4,5 dimerization interface of GPCRs that occurs upon activation. Our results indicate: (1) significant perturbation of the normal modes (NMs) of the rhodopsin monomer upon oligomerization, which is mainly manifested at interfacial regions; (2) increased positive correlation among the transmembrane domains (TMs) and between the extracellular loop (EL) and TM regions of the rhodopsin protomer; (3) highest interresidue positive correlation at the interfaces between protomers; and (4) experimentally testable hypotheses of differential motional changes within different putative oligomeric arrangements.     

The study of G-protein coupled receptor oligomerization with computational modeling and bioinformatics

To achieve a structural context for the analysis of G-protein coupled receptor (GPCR) oligomers, molecular modeling must be used to predict the corresponding interaction interfaces. The task is complicated by the paucity of detailed structural data at atomic resolution, and the large number of possible modes in which the bundles of seven transmembrane (TM) segments of the interacting GPCR monomers can be packed together into dimers and/or higher-order oligomers. Approaches and tools offered by bioinformatics can be used to reduce the complexity of this task and, combined with computational modeling, can serve to yield testable predictions for the structural properties of oligomers. Most of the bioinformatics methods take advantage of the evolutionary relation that exists among GPCRs, as expressed in their sequences and measurable in the common elements of their structural and functional features. These common elements are responsible for the presence of detectable patterns of motifs and correlated mutations evident from the alignment of the sequences of these complex biological systems. The decoding of these patterns in terms of structural and functional determinants can provide indications about the most likely interfaces of dimerization/oligomerization of GPCRs. We review here the main approaches from bioinformatics, enhanced by computational molecular modeling, that have been used to predict likely interfaces of dimerization/oligomerization of GPCRs, and compare results from their application to rhodopsin-like GPCRs. A compilation of the most frequently predicted GPCR oligomerization interfaces points to specific regions of TMs 4-6.     

Structural models for dimerization of G-protein coupled receptors: the opioid receptor homodimers

Among the most exciting functional features of G-protein coupled receptors (GPCRs) that are coming into focus lately are those relating to the role and structural characteristics of their oligomerization (mostly homo- and heterodimers). The structural underpinnings of these novel functional insights are still not clear, as current experimental techniques have not yet succeeded in identifying the dimerization interfaces between GPCR monomers. Two computational approaches have recently been designed in our lab to provide reasonable three-dimensional (3D) molecular models of the transmembrane (TM) regions of GPCR dimers based on a combination of the structural information of receptor monomers and analyses of correlated mutations in receptor families. The modeling of GPCR heterodimers has been described recently. We present here a related approach for modeling of GPCR homodimers that identifies the interfaces in the most likely configurations of the complexes. The approach is illustrated for the three cloned opioid receptor subtypes (OPRD, OPRM, and OPRK).     

Factor analysis for gene regulatory networks and transcription factor activity profiles

Background

G Protein-Coupled Receptors (GPCRs) are a large and diverse family of membrane proteins whose members participate in the regulation of most cellular and physiological processes and therefore represent key pharmacological targets. Although several bioinformatics resources support research on GPCRs, most of these have been designed based on the traditional assumption that monomeric GPCRs constitute the functional receptor unit. The increase in the frequency and number of reports about GPCR dimerization/oligomerization and the implication of oligomerization in receptor function makes necessary the ability to store and access information about GPCR dimers/oligomers electronically.

Results

We present here the requirements and ontology (the information scheme to describe oligomers and associated concepts and their relationships) for an information system that can manage the elements of information needed to describe comprehensively the phenomena of both homo- and hetero-oligomerization of GPCRs. The comprehensive information management scheme that we plan to use for the development of an intuitive and user-friendly GPCR-Oligomerization Knowledge Base (GPCR-OKB) is the result of a community dialog involving experimental and computational colleagues working on GPCRs.

Conclusion

Our long term goal is to disseminate to the scientific community organized, curated, and detailed information about GPCR dimerization/oligomerization and its related structural context. This information will be reported as close to the data as possible so the user can make his own judgment on the conclusions drawn for a particular study. The requirements and ontology described here will facilitate the development of future information systems for GPCR oligomers that contain both computational and experimental information about GPCR oligomerization. This information is freely accessible at http://www.gpcr-okb.org.     

GGIP: Structure and sequence‐based GPCR–GPCR interaction pair predictor

G Protein‐Coupled Receptors (GPCRs) are important pharmaceutical targets. More than 30% of currently marketed pharmaceutical medicines target GPCRs. Numerous studies have reported that GPCRs function not only as monomers but also as homo‐ or hetero‐dimers or higher‐order molecular complexes. Many GPCRs exert a wide variety of molecular functions by forming specific combinations of GPCR subtypes. In addition, some GPCRs are reportedly associated with diseases. GPCR oligomerization is now recognized as an important event in various biological phenomena, and many researchers are investigating this subject. We have developed a support vector machine (SVM)‐based method to predict interacting pairs for GPCR oligomerization, by integrating the structure and sequence information of GPCRs. The performance of our method was evaluated by the Receiver Operating Characteristic (ROC) curve. The corresponding area under the curve was 0.938. As far as we know, this is the only prediction method for interacting pairs among GPCRs. Our method could accelerate the analyses of these interactions, and contribute to the elucidation of the global structures of the GPCR networks in membranes. Proteins 2016; 84:1224–1233. © 2016 Wiley Periodicals, Inc.     

Use of fluorescence resonance energy transfer to analyze oligomerization of G-protein-coupled receptors expressed in yeast

Oligomerization or dimerization of G-protein-coupled receptors (GPCRs) has emerged as an important theme in signal transduction. This concept has recently gained widespread interest due to the application of direct and noninvasive biophysical techniques such as fluorescence resonance energy transfer (FRET), which have shown unequivocally that several types of GPCR can form dimers or oligomers in living cells. Current challenges are to determine which GPCRs can self-associate and/or interact with other GPCRs, to define the molecular principles that govern these specific interactions, and to establish which aspects of GPCR function require oligomerization. Although these questions ultimately must be addressed by using GPCRs expressed endogenously in their native cell types, analysis of GPCR oligomerization in heterologous expression systems will be useful to survey which GPCRs can interact, to conduct structure-function studies, and to identify peptides or small molecules that disrupt GPCR oligomerization and function. Here, we describe methods employing scanning fluorometry to detect FRET between GPCRs tagged with enhanced cyan and yellow fluorescent proteins (CFP and YFP) in living yeast cells. This approach provides a powerful means to analyze oligomerization of a variety of GPCRs that can be expressed in yeast, such as adrenergic, adenosine, C5a, muscarinic acetylcholine, vasopressin, opioid, and somatostatin receptors.     

Prediction of heterodimerization interfaces of G-protein coupled receptors with a new subtractive correlated mutation method

Recent studies employing differential epitope tagging, selective immunoprecipitation of receptor complexes and fluorescence or bioluminescence resonance energy transfer techniques provide direct evidence for heterodimerization between both closely and distantly related members of the G-protein coupled receptor (GPCR) family. Since heterodimerization appears to play a role in modulating agonist affinity, efficacy and/or trafficking properties, the molecular models of GPCRs required to understand receptor function must consider these oligomerization hypotheses. To advance knowledge in this field, we present here a computational approach based on correlated mutation analysis and the structural information contained in three-dimensional molecular models of the transmembrane regions of GPCRs built using the rhodopsin crystal structure as a template. The new subtractive correlated mutation method reveals likely heterodimerization interfaces amongst the different alternatives for the positioning of two tightly packed bundles of seven transmembrane domains next to each other in contact heterodimers of GPCRs. Predictions are applied to GPCRs in the class of opioid receptors. However, in the absence of a known structure of any GPCR dimer, the features of the method and predictions are also illustrated and analyzed for a dimeric complex of known structure.     

Increasingly accurate dynamic molecular models of G-protein coupled receptor oligomers: Panacea or Pandora's box for novel drug discovery?

For years, conventional drug design at G-protein coupled receptors (GPCRs) has mainly focused on the inhibition of a single receptor at a usually well-defined ligand-binding site. The recent discovery of more and more physiologically relevant GPCR dimers/oligomers suggests that selectively targeting these complexes or designing small molecules that inhibit receptor–receptor interactions might provide new opportunities for novel drug discovery. To uncover the fundamental mechanisms and dynamics governing GPCR dimerization/oligomerization, it is crucial to understand the dynamic process of receptor–receptor association, and to identify regions that are suitable for selective drug binding. This minireview highlights current progress in the development of increasingly accurate dynamic molecular models of GPCR oligomers based on structural, biochemical, and biophysical information that has recently appeared in the literature. In view of this new information, there has never been a more exciting time for computational research into GPCRs than at present. Information-driven modern molecular models of GPCR complexes are expected to efficiently guide the rational design of GPCR oligomer-specific drugs, possibly allowing researchers to reach for the high-hanging fruits in GPCR drug discovery, i.e. more potent and selective drugs for efficient therapeutic interventions.     

The extracellular N-terminal domain and transmembrane domains 1 and 2 mediate oligomerization of a yeast G protein-coupled receptor

G protein-coupled receptors (GPCRs) can form homodimers/oligomers and/or heterodimers/oligomers. The mechanisms used to form specific GPCR oligomers are poorly understood because the domains that mediate such interactions and the step(s) in the secretory pathway where oligomerization occurs have not been well characterized. Here we have used subcellular fractionation and fluorescence resonance energy transfer (FRET) experiments to show that oligomerization of a GPCR (alpha-factor receptor; STE2 gene product) of the yeast Saccharomyces cerevisiae occurs in the endoplasmic reticulum. To identify domains of this receptor that mediate oligomerization, we used FRET and endocytosis assays of oligomerization in vivo to analyze receptor deletion mutants. A mutant lacking the N-terminal extracellular domain and transmembrane (TM) domain 1 was expressed at the cell surface but did not self-associate. In contrast, a receptor fragment containing only the N-terminal extracellular domain and TM1 could self-associate and heterodimerize with wild type receptors. Analysis of other mutants suggested that oligomerization is facilitated by the N-terminal extracellular domain and TM2. Therefore, the N-terminal extracellular domain, TM1, and TM2 appear to stabilize alpha-factor receptor oligomers. These domains may form an interface in contact or domain-swapped oligomers. Similar domains may mediate dimerization of certain mammalian GPCRs.     

The expanding repertoire of receptor activity modifying protein (RAMP) function

Receptor activity modifying proteins (RAMPs) associate with G-protein-coupled receptors (GPCRs) at the plasma membrane and together bind a variety of peptide ligands, serving as a communication interface between the extracellular and intracellular environments. The collection of RAMP-interacting GPCRs continues to expand and now consists of GPCRs from families A, B and C, suggesting that RAMP activity is extremely prevalent. RAMP association with GPCRs can regulate GPCR function by altering ligand binding, receptor trafficking and desensitization, and downstream signaling pathways. Here, we elaborate on these RAMP-dependent mechanisms of GPCR regulation, which provide opportunities for pharmacological intervention.     

Homo- and hetero-oligomerization of G protein-coupled receptors

G protein-coupled receptors (GPCRs) form homo-oligomeric and hetero-oligomeric complexes. This understanding has prompted a re-evaluation of many aspects of GPCR biology, however the concept of receptor complexes has not been fully integrated into the current thinking about GPCR structure and function. Nevertheless, receptor oligomerization is a pivotal aspect of the structure and function of GPCRs that has been shown to have implications for receptor trafficking, signaling, and pharmacology and more intricate models for understanding the physiological roles of these receptors are emerging. Here, we summarize some of the advances made in understanding the structural basis and the functional roles of homo- and hetero- oligomerization in this important group of receptors. Although this discussion focuses primarily on the dopamine receptors, particularly the D2 dopamine receptor, and the opioid and serotonin receptors, we discuss the principles governing the oligomerization of all rhodopsin-like GPCRs and potentially of the entire superfamily of these receptors.     

Structural features of the G-protein/GPCR interactions

Background

The details of the functional interaction between G proteins and the G protein coupled receptors (GPCRs) have long been subjected to extensive investigations with structural and functional assays and a large number of computational studies.

Scope of review

The nature and sites of interaction in the G-protein/GPCR complexes, and the specificities of these interactions selecting coupling partners among the large number of families of GPCRs and G protein forms, are still poorly defined.

Major conclusions

Many of the contact sites between the two proteins in specific complexes have been identified, but the three dimensional molecular architecture of a receptor-Gα interface is only known for one pair. Consequently, many fundamental questions regarding this macromolecular assembly and its mechanism remain unanswered.

General significance

In the context of current structural data we review the structural details of the interfaces and recognition sites in complexes of sub-family A GPCRs with cognate G-proteins, with special emphasis on the consequences of activation on GPCR structure, the prevalence of preassembled GPCR/G-protein complexes, the key structural determinants for selective coupling and the possible involvement of GPCR oligomerization in this process.     

A crystal clear solution for determining G-protein-coupled receptor structures

G-protein-coupled receptors (GPCRs) are medically important membrane proteins that are targeted by over 30% of small molecule drugs. At the time of writing, 15 unique GPCR structures have been determined, with 77 structures deposited in the PDB database, which offers new opportunities for drug development and for understanding the molecular mechanisms of GPCR activation. Many different factors have contributed to this success, but if there is one single factor that can be singled out as the foundation for producing well-diffracting GPCR crystals, it is the stabilisation of the detergent-solubilised receptor-ligand complex. This review will focus predominantly on one of the successful strategies for the stabilisation of GPCRs, namely the thermostabilisation of GPCRs using systematic mutagenesis coupled with thermostability assays. Structures of thermostabilised GPCRs bound to a wide variety of ligands have been determined, which has led to an understanding of ligand specificity; why some ligands act as agonists as opposed to partial or inverse agonists; and the structural basis for receptor activation.     

Structure modeling of all identified G protein-coupled receptors in the human genome

G protein-coupled receptors (GPCRs), encoded by about 5% of human genes, comprise the largest family of integral membrane proteins and act as cell surface receptors responsible for the transduction of endogenous signal into a cellular response. Although tertiary structural information is crucial for function annotation and drug design, there are few experimentally determined GPCR structures. To address this issue, we employ the recently developed threading assembly refinement (TASSER) method to generate structure predictions for all 907 putative GPCRs in the human genome. Unlike traditional homology modeling approaches, TASSER modeling does not require solved homologous template structures; moreover, it often refines the structures closer to native. These features are essential for the comprehensive modeling of all human GPCRs when close homologous templates are absent. Based on a benchmarked confidence score, approximately 820 predicted models should have the correct folds. The majority of GPCR models share the characteristic seven-transmembrane helix topology, but 45 ORFs are predicted to have different structures. This is due to GPCR fragments that are predominantly from extracellular or intracellular domains as well as database annotation errors. Our preliminary validation includes the automated modeling of bovine rhodopsin, the only solved GPCR in the Protein Data Bank. With homologous templates excluded, the final model built by TASSER has a global C(alpha) root-mean-squared deviation from native of 4.6 angstroms, with a root-mean-squared deviation in the transmembrane helix region of 2.1 angstroms. Models of several representative GPCRs are compared with mutagenesis and affinity labeling data, and consistent agreement is demonstrated. Structure clustering of the predicted models shows that GPCRs with similar structures tend to belong to a similar functional class even when their sequences are diverse. These results demonstrate the usefulness and robustness of the in silico models for GPCR functional analysis. All predicted GPCR models are freely available for noncommercial users on our Web site (http://www.bioinformatics.buffalo.edu/GPCR).     

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.     

G protein activation by serotonin type 4 receptor dimers: evidence that turning on two protomers is more efficient

The discovery that class C G protein-coupled receptors (GPCRs) function as obligatory dimeric entities has generated major interest in GPCR oligomerization. Oligomerization now appears to be a common feature among all GPCR classes. However, the functional significance of this process remains unclear because, in vitro, some monomeric GPCRs, such as rhodopsin and β(2)-adrenergic receptors, activate G proteins. By using wild type and mutant serotonin type 4 receptors (5-HT(4)Rs) (including a 5-HT(4)-RASSL) expressed in COS-7 cells as models of class A GPCRs, we show that activation of one protomer in a dimer was sufficient to stimulate G proteins. However, coupling efficiency was 2 times higher when both protomers were activated. Expression of combinations of 5-HT(4), in which both protomers were able to bind to agonists but only one could couple to G proteins, suggested that upon agonist occupancy, protomers did not independently couple to G proteins but rather that only one G protein was activated. Coupling of a single heterotrimeric G(s) protein to a receptor dimer was further confirmed in vitro, using the purified recombinant WT RASSL 5-HT(4)R obligatory heterodimer. These results, together with previous findings, demonstrate that, differently from class C GPCR dimers, class A GPCR dimers have pleiotropic activation mechanisms.     

Formation of new genes explains lower intron density in mammalian Rhodopsin G protein-coupled receptors

Mammalian G protein-coupled receptor (GPCR) genes are characterised by a large proportion of intronless genes or a lower density of introns when compared with GPCRs of invertebrates. It is unclear which mechanisms have influenced intron density in this protein family, which is one of the largest in the mammalian genomes. We used a combination of Hidden Markov Models (HMM) and BLAST searches to establish the comprehensive repertoire of Rhodopsin GPCRs from seven species and performed overall alignments and phylogenetic analysis using the maximum parsimony method for over 1400 receptors in 12 subgroups. We identified 14 different Ancestral Receptor Groups (ARGs) that have members in both vertebrate and invertebrate species. We found that there exists a remarkable difference in the intron density among ancestral and new Rhodopsin GPCRs. The intron density among ARGs members was more than 3.5-fold higher than that within non-ARG members and more than 2-fold higher when considering only the 7TM region. This suggests that the new GPCR genes have been predominantly formed intronless while the ancestral receptors likely accumulated introns during their evolution. Many of the intron positions found in mammalian ARG receptor sequences were found to be present in orthologue invertebrate receptors suggesting that these intron positions are ancient. This analysis also revealed that one intron position is much more frequent than any other position and it is common for a number of phylogenetically different Rhodopsin GPCR groups. This intron position lies within a functionally important, conserved, DRY motif which may form a proto-splice site that could contribute to positional intron insertion. Moreover, we have found that other receptor motifs, similar to DRY, also contain introns between the second and third nucleotide of the arginine codon which also forms a proto-splice site. Our analysis presents compelling evidence that there was not a major loss of introns in mammalian GPCRs and formation of new GPCRs among mammals explains why these have fewer introns compared to invertebrate GPCRs. We also discuss and speculate about the possible role of different RNA- and DNA-based mechanisms of intron insertion and loss.     

On the expanding terminology in the GPCR field: The meaning of receptor mosaics and receptor heteromers

The oligomerization of G protein-coupled receptors (GPCRs) is a fact that deserves further attention as increases both the complexity and diversity of the receptor-mediated signal transduction, thus enriching the cell signaling. Consequently, in the present review we tackle among others the problems concerning the terminology used to describe aspects surrounding the GPCRs oligomerization phenomenon. Therefore, the theoretical implications of the GPCR oligomerization will be briefly discussed together with possible implications of this phenomenon especially for new strategies in drug development.