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

Structure and dynamics determine G protein coupling specificity at a class A GPCR

G protein coupled receptors (GPCRs) exhibit varying degrees of selectivity for different G protein isoforms. Despite the abundant structures of GPCR-G protein complexes, little is known about the mechanism of G protein coupling specificity. The β2-adrenergic receptor is an example of GPCR with high selectivity for Gαs, the stimulatory G protein for adenylyl cyclase, and much weaker for the Gαi family of G proteins inhibiting adenylyl cyclase. By developing a new Gαi-biased agonist (LM189), we provide structural and biophysical evidence supporting that distinct conformations at ICL2 and TM6 are required for coupling of the different G protein subtypes Gαs and Gαi. These results deepen our understanding of G protein specificity and bias and can accelerate the design of ligands that select for preferred signaling pathways.

Biophysical investigations of class A GPCRs

G protein-coupled receptors (GPCRs) play a central role in cellular communication, converting external stimuli into intracellular responses. GPCRs bind a very broad panel of ligands, such as hormones, neurotransmitters, peptides and lipids. Ligand binding triggers a series of receptor conformational rearrangements, enabling the coupling to intracellular partners and the activation of signaling cascades. The major breakthrough in GPCRs structural biology of the past decade has considerably advanced our understanding of GPCR activation. However, structural information cannot fully explain the molecular details of GPCRs pharmacology. Biophysical investigations reveal that GPCRs are very dynamic proteins, capable of exploring a wide range of conformational states. Binding to ligands of various pharmacological classes, as well as intracellular effectors and allosteric modulators, can shift the equilibrium between these states and the kinetic of interconversions among the different conformers. Investigation of GPCR dynamic interplay is therefore important to better understand the complex pharmacology and signaling profile of these receptors.

Understanding G Protein Selectivity of Muscarinic Acetylcholine Receptors Using Computational Methods

The neurotransmitter molecule acetylcholine is capable of activating five muscarinic acetylcholine receptors, M1 through M5, which belong to the superfamily of G-protein-coupled receptors (GPCRs). These five receptors share high sequence and structure homology; however, the M1, M3, and M5 receptor subtypes signal preferentially through the Gαq/11 subset of G proteins, whereas the M2 and M4 receptor subtypes signal through the Gαi/o subset of G proteins, resulting in very different intracellular signaling cascades and physiological effects. The structural basis for this innate ability of the M1/M3/M5 set of receptors and the highly homologous M2/M4 set of receptors to couple to different G proteins is poorly understood. In this study, we used molecular dynamics (MD) simulations coupled with thermodynamic analyses of M1 and M2 receptors coupled to both Gαi and Gαq to understand the structural basis of the M1 receptor’s preference for the Gαq protein and the M2 receptor’s preference for the Gαi protein. The MD studies showed that the M1 and M2 receptors can couple to both Gα proteins such that the M1 receptor engages with the two Gα proteins in slightly different orientations and the M2 receptor engages with the two Gα proteins in the same orientation. Thermodynamic studies of the free energy of binding of the receptors to the Gα proteins showed that the M1 and M2 receptors bind more strongly to their cognate Gα proteins compared to their non-cognate ones, which is in line with previous experimental studies on the M3 receptor. A detailed analysis of receptor–G protein interactions showed some cognate-complex-specific interactions for the M2:Gαi complex; however, G protein selectivity determinants are spread over a large overlapping subset of residues. Conserved interaction between transmembrane helices 5 and 6 far away from the G-protein-binding receptor interface was found only in the two cognate complexes and not in the non-cognate complexes. An analysis of residues implicated previously in G protein selectivity, in light of the cognate and non-cognate structures, shaded a more nuanced role of those residues in affecting G protein selectivity. The simulation of both cognate and non-cognate receptor–G protein complexes fills a structural gap due to difficulties in determining non-cognate complex structures and provides an enhanced framework to probe the mechanisms of G protein selectivity exhibited by most GPCRs.

Illuminating G-Protein-Coupling Selectivity of GPCRs

Heterotrimetic G proteins consist of four subfamilies (G_s, G_i/o, G_q/11, and G_12/13) that mediate signaling via G-protein-coupled receptors (GPCRs), principally by receptors binding Gα C termini. G-protein-coupling profiles govern GPCR-induced cellular responses, yet receptor sequence selectivity determinants remain elusive. Here, we systematically quantified ligand-induced interactions between 148 GPCRs and all 11 unique Gα subunit C termini. For each receptor, we probed chimeric Gα subunit activation via a transforming growth factor-α (TGF-α) shedding response in HEK293 cells lacking endogenous G_q/11 and G_12/13 proteins, and complemented G-protein-coupling profiles through a NanoBiT-G-protein dissociation assay. Interrogation of the dataset identified sequence-based coupling specificity features, inside and outside the transmembrane domain, which we used to develop a coupling predictor that outperforms previous methods. We used the predictor to engineer designer GPCRs selectively coupled to G₁₂. This dataset of fine-tuned signaling mechanisms for diverse GPCRs is a valuable resource for research in GPCR signaling.

Light control of G protein signaling pathways by a novel photopigment

Channelopsins and photo-regulated ion channels make it possible to use light to control electrical activity of cells. This powerful approach has lead to a veritable explosion of applications, though it is limited to changing membrane voltage of the target cells. An enormous potential could be tapped if similar opto-genetic techniques could be extended to the control of chemical signaling pathways. Photopigments from invertebrate photoreceptors are an obvious choice—as they do not bleach upon illumination -however, their functional expression has been problematic. We exploited an unusual opsin, pScop2, recently identified in ciliary photoreceptors of scallop. Phylogenetically, it is closer to vertebrate opsins, and offers the advantage of being a bi-stable photopigment. We inserted its coding sequence and a fluorescent protein reporter into plasmid vectors and demonstrated heterologous expression in various mammalian cell lines. HEK 293 cells were selected as a heterologous system for functional analysis, because wild type cells displayed the largest currents in response to the G-protein activator, GTP-γ-S. A line of HEK cells stably transfected with pScop2 was generated; after reconstitution of the photopigment with retinal, light responses were obtained in some cells, albeit of modest amplitude. In native photoreceptors pScop2 couples to G_o; HEK cells express poorly this G-protein, but have a prominent Gq/PLC pathway linked to internal Ca mobilization. To enhance pScop2 competence to tap into this pathway, we swapped its third intracellular loop—important to confer specificity of interaction between 7TMDRs and G-proteins—with that of a G_q-linked opsin which we cloned from microvillar photoreceptors present in the same retina. The chimeric construct was evaluated by a Ca fluorescence assay, and was shown to mediate a robust mobilization of internal calcium in response to illumination. The results project pScop2 as a potentially powerful optogenetic tool to control signaling pathways.

Evidence That Calcium Entry Into Calcium-Transporting Dental Enamel Cells Is Regulated by Cholecystokinin, Acetylcholine and ATP

Dental enamel is formed by specialized epithelial cells which handle large quantities of Ca²⁺ while producing the most highly mineralized tissue. However, the mechanisms used by enamel cells to handle bulk Ca²⁺ safely remain unclear. Our previous work contradicted the dogma that Ca²⁺ is ferried through the cytosol of Ca²⁺-transporting cells and instead suggested an organelle-based route across enamel cells. This new paradigm involves endoplasmic reticulum (ER)-associated Ca²⁺ stores and their concomitant refilling by store-operated Ca²⁺ entry (SOCE) mediated by Ca²⁺ release activated Ca²⁺ (CRAC) channels. Given that Ca²⁺ handling is maximal during the enamel-mineralization stage (maturation), we anticipated that SOCE would also be elevated then. Confirmation was obtained here using single-cell recordings of cytosolic Ca²⁺ concentration ([Ca²⁺]_cyt) in rat ameloblasts. A candidate SOCE agonist, cholecystokinin (CCK), was found to be upregulated during maturation, with Cck transcript abundance reaching 30% of that in brain. CCK-receptor transcripts were also detected and Ca²⁺ imaging showed that CCK stimulation increased [Ca²⁺]_cyt in a dose-responsive manner that was sensitive to CRAC-channel inhibitors. Similar effects were observed with two other SOCE activators, acetylcholine and ATP, whose receptors were also found in enamel cells. These results provide the first evidence of a potential regulatory system for SOCE in enamel cells and so strengthen the Ca²⁺ transcytosis paradigm for ER-based transport of bulk Ca²⁺. Our findings also implicate enamel cells as a new physiological target of CCK and raise the possibility of an auto/paracrine system for regulating Ca²⁺ transport.

Selectivity determinants of GPCR-G-protein binding

The selective coupling of G-protein-coupled receptors (GPCRs) to specific G proteins is critical to trigger the appropriate physiological response. However, the determinants of selective binding have remained elusive. Here we reveal the existence of a selectivity barcode (that is, patterns of amino acids) on each of the 16 human G proteins that is recognized by distinct regions on the approximately 800 human receptors. Although universally conserved positions in the barcode allow the receptors to bind and activate G proteins in a similar manner, different receptors recognize the unique positions of the G-protein barcode through distinct residues, like multiple keys (receptors) opening the same lock (G protein) using non-identical cuts. Considering the evolutionary history of GPCRs allows the identification of these selectivity-determining residues. These findings lay the foundation for understanding the molecular basis of coupling selectivity within individual receptors and G proteins.

Novel insights on thyroid-stimulating hormone receptor signal transduction

The TSH receptor (TSHR) is a member of the glycoprotein hormone receptors, a subfamily of family A G protein-coupled receptors. The TSHR is of great importance for the growth and function of the thyroid gland. The TSHR and its endogenous ligand TSH are pivotal proteins with respect to a variety of physiological functions and malfunctions. The molecular events of TSHR regulation can be summarized as a process of signal transduction, including signal reception, conversion, and amplification. The steps during signal transduction from the extra- to the intracellular sites of the cell are not yet comprehensively understood. However, essential new insights have been achieved in recent years on the interrelated mechanisms at the extracellular region, the transmembrane domain, and intracellular components. This review contains a critical summary of available knowledge of the molecular mechanisms of signal transduction at the TSHR, for example, the key amino acids involved in hormone binding or in the structural conformational changes that lead to G protein activation or signaling regulation. Aspects of TSHR oligomerization, signaling promiscuity, signaling selectivity, phenotypes of genetic variations, and potential extrathyroidal receptor activity are also considered, because these are relevant to an understanding of the overall function of the TSHR, including physiological, pathophysiological, and pharmacological perspectives. Directions for future research are discussed.

Non-neuronal functions of the m2 muscarinic acetylcholine receptor

Acetylcholine is an important neurotransmitter whose effects are mediated by two classes of receptors. The nicotinic acetylcholine receptors are ion channels, whereas the muscarinic receptors belong to the large family of G protein coupled seven transmembrane helix receptors. Beyond its function in neuronal systems, it has become evident that acetylcholine also plays an important role in non-neuronal cells such as epithelial and immune cells. Furthermore, many cell types in the periphery are capable of synthesizing acetylcholine and express at least some of the receptors. In this review, we summarize the non-neuronal functions of the muscarinic acetylcholine receptors, especially those of the M₂ muscarinic receptor in epithelial cells. We will review the mechanisms of signaling by the M₂ receptor but also the cellular trafficking and ARF6 mediated endocytosis of this receptor, which play an important role in the regulation of signaling events. In addition, we provide an overview of the M₂ receptor in human pathological conditions such as autoimmune diseases and cancer.

Pathophysiology of GPCR Homo- and Heterodimerization: Special Emphasis on Somatostatin Receptors

G-protein coupled receptors (GPCRs) are cell surface proteins responsible for translating >80% of extracellular reception to intracellular signals. The extracellular information in the form of neurotransmitters, peptides, ions, odorants etc is converted to intracellular signals via a wide variety of effector molecules activating distinct downstream signaling pathways. All GPCRs share common structural features including an extracellular N-terminal, seven-transmembrane domains (TMs) linked by extracellular/intracellular loops and the C-terminal tail. Recent studies have shown that most GPCRs function as dimers (homo- and/or heterodimers) or even higher order of oligomers. Protein-protein interaction among GPCRs and other receptor proteins play a critical role in the modulation of receptor pharmacology and functions. Although ~50% of the current drugs available in the market target GPCRs, still many GPCRs remain unexplored as potential therapeutic targets, opening immense possibility to discover the role of GPCRs in pathophysiological conditions. This review explores the existing information and future possibilities of GPCRs as tools in clinical pharmacology and is specifically focused for the role of somatostatin receptors (SSTRs) in pathophysiology of diseases and as the potential candidate for drug discovery.

More than two decades of research on insect neuropeptide GPCRs: an overview

This review focuses on the state of the art on neuropeptide receptors in insects. Most of these receptors are G protein-coupled receptors (GPCRs) and are involved in the regulation of virtually all physiological processes during an insect's life. More than 20 years ago a milestone in invertebrate endocrinology was achieved with the characterization of the first insect neuropeptide receptor, i.e., the Drosophila tachykinin-like receptor. However, it took until the release of the Drosophila genome in 2000 that research on neuropeptide receptors boosted. In the last decade a plethora of genomic information of other insect species also became available, leading to a better insight in the functions and evolution of the neuropeptide signaling systems and their intracellular pathways. It became clear that some of these systems are conserved among all insect species, indicating that they fulfill crucial roles in their physiological processes. Meanwhile, other signaling systems seem to be lost in several insect orders or species, suggesting that their actions were superfluous in those insects, or that other neuropeptides have taken over their functions. It is striking that the deorphanization of neuropeptide GPCRs gets much attention, but the subsequent unraveling of the intracellular pathways they elicit, or their physiological functions are often hardly examined. Especially in insects besides Drosophila this information is scarce if not absent. And although great progress made in characterizing neuropeptide signaling systems, even in Drosophila several predicted neuropeptide receptors remain orphan, awaiting for their endogenous ligand to be determined. The present review gives a précis of the insect neuropeptide receptor research of the last two decades. But it has to be emphasized that the work done so far is only the tip of the iceberg and our comprehensive understanding of these important signaling systems will still increase substantially in the coming years.

Classification of G proteins and prediction of GPCRs-G proteins coupling specificity using continuous wavelet transform and information theory

The coupling between G protein-coupled receptors (GPCRs) and guanine nucleotide-binding proteins (G proteins) regulates various signal transductions from extracellular space into the cell. However, the coupling mechanism between GPCRs and G proteins is still unknown, and experimental determination of their coupling specificity and function is both expensive and time consuming. Therefore, it is significant to develop a theoretical method to predict the coupling specificity between GPCRs and G proteins as well as their function using their primary sequences. In this study, a novel four-layer predictor (GPCRsG_CWTIT) based on support vector machine (SVM), continuous wavelet transform (CWT) and information theory (IT) is developed to classify G proteins and predict the coupling specificity between GPCRs and G proteins. SVM is used for construction of models. CWT and IT are used to characterize the primary structure of protein. Performance of GPCRsG_CWTIT is evaluated with cross-validation test on various working dataset. The overall accuracy of the G proteins at the levels of class and family is 98.23 and 85.42%, respectively. The accuracy of the coupling specificity prediction varies from 74.60 to 94.30%. These results indicate that the proposed predictor is an effective and feasible tool to predict the coupling specificity between GPCRs and G proteins as well as their functions using only the protein full sequence. The establishment of such an accurate prediction method will facilitate drug discovery by improving the ability to identify and predict protein-protein interactions. GPCRsG_CWTIT and dataset can be acquired freely on request from the authors.

Chemogenomic analysis of G-protein coupled receptors and their ligands deciphers locks and keys governing diverse aspects of signalling

Understanding the molecular mechanism of signalling in the important super-family of G-protein-coupled receptors (GPCRs) is causally related to questions of how and where these receptors can be activated or inhibited. In this context, it is of great interest to unravel the common molecular features of GPCRs as well as those related to an active or inactive state or to subtype specific G-protein coupling. In our underlying chemogenomics study, we analyse for the first time the statistical link between the properties of G-protein-coupled receptors and GPCR ligands. The technique of mutual information (MI) is able to reveal statistical inter-dependence between variations in amino acid residues on the one hand and variations in ligand molecular descriptors on the other. Although this MI analysis uses novel information that differs from the results of known site-directed mutagenesis studies or published GPCR crystal structures, the method is capable of identifying the well-known common ligand binding region of GPCRs between the upper part of the seven transmembrane helices and the second extracellular loop. The analysis shows amino acid positions that are sensitive to either stimulating (agonistic) or inhibitory (antagonistic) ligand effects or both. It appears that amino acid positions for antagonistic and agonistic effects are both concentrated around the extracellular region, but selective agonistic effects are cumulated between transmembrane helices (TMHs) 2, 3, and ECL2, while selective residues for antagonistic effects are located at the top of helices 5 and 6. Above all, the MI analysis provides detailed indications about amino acids located in the transmembrane region of these receptors that determine G-protein signalling pathway preferences.

Aminergic control and modulation of honeybee behaviour

Biogenic amines are important messenger substances in the central nervous system and in peripheral organs of vertebrates and of invertebrates. The honeybee, Apis mellifera, is excellently suited to uncover the functions of biogenic amines in behaviour, because it has an extensive behavioural repertoire, with a number of biogenic amine receptors characterised in this insect.In the honeybee, the biogenic amines dopamine, octopamine, serotonin and tyramine modulate neuronal functions in various ways. Dopamine and serotonin are present in high concentrations in the bee brain, whereas octopamine and tyramine are less abundant. Octopamine is a key molecule for the control of honeybee behaviour. It generally has an arousing effect and leads to higher sensitivity for sensory inputs, better learning performance and increased foraging behaviour. Tyramine has been suggested to act antagonistically to octopamine, but only few experimental data are available for this amine. Dopamine and serotonin often have antagonistic or inhibitory effects as compared to octopamine.Biogenic amines bind to membrane receptors that primarily belong to the large gene-family of GTP-binding (G) protein coupled receptors. Receptor activation leads to transient changes in concentrations of intracellular second messengers such as cAMP, IP(3) and/or Ca(2+). Although several biogenic amine receptors from the honeybee have been cloned and characterised more recently, many genes still remain to be identified. The availability of the completely sequenced genome of Apis mellifera will contribute substantially to closing this gap.In this review, we will discuss the present knowledge on how biogenic amines and their receptor-mediated cellular responses modulate different behaviours of honeybees including learning processes and division of labour.

Principles and determinants of G-protein coupling by the rhodopsin-like thyrotropin receptor

In this study we wanted to gain insights into selectivity mechanisms between G-protein-coupled receptors (GPCR) and different subtypes of G-proteins. The thyrotropin receptor (TSHR) binds G-proteins promiscuously and activates both Gs (cAMP) and Gq (IP). Our goal was to dissect selectivity patterns for both pathways in the intracellular region of this receptor. We were particularly interested in the participation of poorly investigated receptor parts.

We systematically investigated the amino acids of intracellular loop (ICL) 1 and helix 8 using site-directed mutagenesis alongside characterization of cAMP and IP accumulation. This approach was guided by a homology model of activated TSHR in complex with heterotrimeric Gq, using the X-ray structure of opsin with a bound G-protein peptide as a structural template.

We provide evidence that ICL1 is significantly involved in G-protein activation and our model suggests potential interactions with subunits Gα as well as Gβγ. Several amino acid substitutions impaired both IP and cAMP accumulation. Moreover, we found a few residues in ICL1 (L440, T441, H443) and helix 8 (R687) that are sensitive for Gq but not for Gs activation. Conversely, not even one residue was found that selectively affects cAMP accumulation only.

Together with our previous mutagenesis data on ICL2 and ICL3 we provide here the first systematically completed map of potential interfaces between TSHR and heterotrimeric G-protein. The TSHR/Gq-heterotrimer complex is characterized by more selective interactions than the TSHR/Gs complex. In fact the receptor interface for binding Gs is a subset of that for Gq and we postulate that this may be true for other GPCRs coupling these G-proteins. Our findings support that G-protein coupling and preference is dominated by specific structural features at the intracellular region of the activated GPCR but is completed by additional complementary recognition patterns between receptor and G-protein subtypes.

Observation of "ionic lock" formation in molecular dynamics simulations of wild-type beta 1 and beta 2 adrenergic receptors

G protein coupled receptors (GPCRs) are a large family of integral membrane proteins involved in signal transduction pathways, making them appealing drug targets for a wide spectrum of diseases. The recently crystallized structures of two engineered adrenergic receptors have opened new avenues for the understanding of the molecular mechanisms of action of GPCRs. Taking the two crystal structures as a starting point, we carried out submicrosecond molecular dynamics simulations of wild-type beta(1) and beta(2) adrenergic receptors in a lipid bilayer under physiological conditions. These simulations give access to structural and dynamic properties of the receptors in pseudo in vivo conditions. For both systems the overall fold properties of the transmembrane region as well as the binding pocket remain close to the crystal structure of the engineered systems, thus suggesting that the ligand binding mode is not affected by the introduced modifications. Both simulations indicate the presence of one or two internal water molecules absent in both crystal structures and essential for the stabilization of the binding pocket at the interface between transmembrane helices III, IV, and V. The different interactions arising from the substitution of Tyr308 in beta(2)AR into Phe325 in beta(1)AR induce different conformations of the homologous Asn(6.55) inside the binding pockets of the two receptors, suggesting a possible origin of receptor specificity in agonist binding. The equilibrated structures of both receptors recover all of the previously suggested features of inactive GPCRs including formation of a salt bridge between the cytoplasmatic moieties of helices III and VI ("ionic lock") that is absent in the crystal structures.

Immunolocalization of a Galphaq protein to the chemosensory organs of Dipolydora quadrilobata (polychaeta: spionidae)

Chemoreception in marine invertebrates mediates a variety of ecologically important behaviors including defense, reproduction, larval settlement and recruitment, and feeding. The sensory pathways that regulate deposit-feeding activity by polychaetes living in sedimentary habitats are of particular interest because such feeding has profound effects on the physical and chemical properties of the habitat. Nevertheless, little is known concerning the molecular mechanisms of chemical signal transduction associated with deposit feeding and other behaviors in polychaetes. Chemosensory-based feeding behaviors are typically regulated by G-protein-coupled signal transduction pathways. However, the presence and role of such pathways have not been demonstrated in marine polychaetes. Methodologies involving degenerate primer-based reverse transcription with the polymerase chain reaction and rapid amplification of cDNA ends were used to identify and characterize a Galphaq subunit expressed in the feeding palps of the spionid polychaete Dipolydora quadrilobata. The D. quadrilobata Galphaq protein had high sequence similarity with previously reported Galphaq subunits from both invertebrate and vertebrate taxa. Immunhistochemistry and immunocytochemistry were used with confocal laser scanning microscopy and transmission electron microscopy to visualize the distribution of a Galphaq antibody in whole worms and in cilia of the feeding palps. Galphaq immunoreactivity was concentrated in the nuchal organs, food-groove cilia, and lateral/abfrontal cilia of the feeding palps. Because these structures are known to be involved in chemoreception, we propose that Galphaq isolated from D. quadrilobata is a key component of chemosensory signal transduction pathways in this species.

Predicting G-protein coupled receptors-G-protein coupling specificity based on autocross-covariance transform

Determining G-protein coupled receptors (GPCRs) coupling specificity is very important for further understanding the functions of receptors. A successful method in this area will benefit both basic research and drug discovery practice. Previously published methods rely on the transmembrane topology prediction at training step, even at prediction step. However, the transmembrane topology predicted by even the best algorithm is not of high accuracy. In this study, we developed a new method, autocross-covariance (ACC) transform based support vector machine (SVM), to predict coupling specificity between GPCRs and G-proteins. The primary amino acid sequences are translated into vectors based on the principal physicochemical properties of the amino acids and the data are transformed into a uniform matrix by applying ACC transform. SVMs for nonpromiscuous coupled GPCRs and promiscuous coupled GPCRs were trained and validated by jackknife test and the results thus obtained are very promising. All classifiers were also evaluated by the test datasets with good performance. Besides the high prediction accuracy, the most important feature of this method is that it does not require any transmembrane topology prediction at either training or prediction step but only the primary sequences of proteins. The results indicate that this relatively simple method is applicable. Academic users can freely download the prediction program at http://www.scucic.net/group/database/Service.asp.

Molecular mechanisms of constitutive activity: mutations at position 111 of the angiotensin AT1 receptor

A possible molecular mechanism for the constitutive activity of mutants of the angiotensin type 1 receptor (AT1) at position 111 was suggested by molecular modeling. This involves a cascade of conformational changes in spatial positions of side chains along transmembrane helix (TM3) from L112 to Y113 to F117, which in turn, results in conformational changes in TM4 (residues I152 and M155) leading to the movement of TM4 as a whole. The mechanism is consistent with the available data of site-directed mutagenesis, as well as with correct predictions of constitutive activity of mutants L112F and L112C. It was also predicted that the double mutant N111G/L112A might possess basal constitutive activity comparable with that of the N111G mutant, whereas the double mutants N111G/Y113A, N111G/F117A, and N111G/I152A would have lower levels of basal activity. Experimental studies of the above double mutants showed significant constitutive activity of N111G/L112A and N111G/F117A. The basal activity of N111G/I152A was higher than expected, and that of N111G/Y113A was not determined due to poor expression of the mutant. The proposed mechanism of constitutive activity of the AT(1) receptor reveals a novel nonsimplistic view on the general problem of constitutive activity, and clearly demonstrates the inherent complexity of the process of G protein-coupled receptor (GPCR) activation.

AmphiD1/β, a dopamine D1/β-adrenergic receptor from the amphioxus Branchiostoma floridae: evolutionary aspects of the catecholaminergic system during development

Catecholamine receptors mediate wide-ranging functions in vertebrates and invertebrates but are largely unknown in invertebrate chordates such as amphioxus. Catecholaminergic cells have been described in amphioxus adults, but few data are known about the transmembrane signal transduction pathways and the expression pattern of related receptors during development. In Branchiostoma floridae, we cloned a full-length cDNA (AmphiD1/beta) that corresponds to the dopamine D1/beta receptor previously cloned from a related species of amphioxus, Branchiostoma lanceolatum, but no expression studies have been performed for such receptor in amphioxus. In B. floridae, AmphiD1/beta encodes a polypeptide with typical G-protein-coupled receptor features, characterized by highest sequence similarity with D1 dopamine and beta-adrenergic receptors. The expression of AmphiD1/beta mRNA in different regions of the cerebral vesicle corresponds to that of D1-like receptors in vertebrate homologous structures. Furthermore, in situ experiments show that during development, the expression in the nervous system is restricted to cells located anteriorly. A further expression was found in larvae at the level of the endostyle, but it has no counterpart in the predominant expression domains of vertebrate dopamine and/or adrenergic receptor genes. At the same time, we compared the dopaminergic system, consisting of AmphiTH-expressing cells, with the AmphiD1/beta expression. In conclusion, the identification of the AmphiD1/beta receptor provides further basis for understanding the evolutionary history of the dopaminergic system at the transition from invertebrates and vertebrates.

Location and nature of the residues important for ligand recognition in G-protein coupled receptors

The overall structure of the biogenic amine subclass of the G-protein-coupled receptors, and of their ligand binding sites, is discussed with the aim of highlighting the major structural features of these receptors that are responsible for ligand recognition. A comparison is made between biogenic amine receptors, peptide receptors of the rhodopsin class, and the secretin receptors which all have peptide ligands. The question of where the peptide ligands bind, whether at extracellular sites or within the transmembrane helix bundle, is discussed. The suitability of the rhodopsin crystal structure as a template for construction of homology models is discussed and it is concluded that there are many reasons why a caution should be issued against using it uncritically.

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.

The serotonin1A receptor: a representative member of the serotonin receptor family

1. Serotonin is an intrinsically fluorescent biogenic amine that acts as a neurotransmitter and is found in a wide variety of sites in the central and peripheral nervous system. Serotonergic signaling appears to play a key role in the generation and modulation of various cognitive and behavioral functions. 2. Serotonin exerts its diverse actions by binding to distinct cell surface receptors which have been classified into many groups. The serotonin1A (5-HT1A) receptor is the most extensively studied of the serotonin receptors and belongs to the large family of seven transmembrane domain G-protein coupled receptors. 3. The tissue and sub-cellular distribution, structural characteristics, signaling of the serotonin1A receptor and its interaction with G-proteins are discussed. 4. The pharmacology of serotonin1A receptors is reviewed in terms of binding of agonists and antagonists and sensitivity of their binding to guanine nucleotides. 5. Membrane biology of 5-HT1A receptors is presented using the bovine hippocampal serotonin1A receptor as a model system. The ligand binding activity and G-protein coupling of the receptor is modulated by membrane cholesterol thereby indicating the requirement of cholesterol in maintaining the receptor organization and function. This, along with the reported detergent resistance characteristics of the receptor, raises important questions on the role of membrane lipids and domains in the function of this receptor.

A method for the prediction of GPCRs coupling specificity to G-proteins using refined profile Hidden Markov Models

BACKGROUND

G- Protein coupled receptors (GPCRs) comprise the largest group of eukaryotic cell surface receptors with great pharmacological interest. A broad range of native ligands interact and activate GPCRs, leading to signal transduction within cells. Most of these responses are mediated through the interaction of GPCRs with heterotrimeric GTP-binding proteins (G-proteins). Due to the information explosion in biological sequence databases, the development of software algorithms that could predict properties of GPCRs is important. Experimental data reported in the literature suggest that heterotrimeric G-proteins interact with parts of the activated receptor at the transmembrane helix-intracellular loop interface. Utilizing this information and membrane topology information, we have developed an intensive exploratory approach to generate a refined library of statistical models (Hidden Markov Models) that predict the coupling preference of GPCRs to heterotrimeric G-proteins. The method predicts the coupling preferences of GPCRs to Gs, Gi/o and Gq/11, but not G12/13 subfamilies.

RESULTS

Using a dataset of 282 GPCR sequences of known coupling preference to G-proteins and adopting a five-fold cross-validation procedure, the method yielded an 89.7% correct classification rate. In a validation set comprised of all receptor sequences that are species homologues to GPCRs with known coupling preferences, excluding the sequences used to train the models, our method yields a correct classification rate of 91.0%. Furthermore, promiscuous coupling properties were correctly predicted for 6 of the 24 GPCRs that are known to interact with more than one subfamily of G-proteins.

CONCLUSION

Our method demonstrates high correct classification rate. Unlike previously published methods performing the same task, it does not require any transmembrane topology prediction in a preceding step. A web-server for the prediction of GPCRs coupling specificity to G-proteins available for non-commercial users is located at http://bioinformatics.biol.uoa.gr/PRED-COUPLE.

Modeling and simulation of the human delta opioid receptor

A model for the human delta opioid receptor has been generated via sequence alignment, structure building using the crystal structure of bovine rhodopsin as a template, and refinement by molecular dynamics simulation. The model building suggested that, in addition to the previously postulated interaction between D128 and Y308, an internal salt bridge also exists between residues D128 and R192, both of which are conserved in all the opioid receptors. The model and salt bridge were then shown to be stable during a 20-nsec simulation in a lipid bilayer. It is therefore proposed that both of these interactions play a role in stabilizing the inactive state of the receptor. The model is also used in an effort to rationalize many of the mutational studies performed on delta opioid receptors, and to suggest a plausible explanation for the differences between known delta opioid agonists and antagonists.

Novel approaches for modeling of the A1adenosine receptor and its agonist binding site

The present work describes the building of a human A(1) adenosine receptor (hA(1)AR) model, based on the X-ray crystal structure of bovine rhodopsin, and its use as a basis for the investigation of some important structural characteristics of the receptor. One of the issues investigated was the protonation position of two histidine residues known to influence ligand binding, with protonation of His251 (6.52) in epsilon position and His278 (7.43) in delta position showing the best agreement with experimental evidence. The model was also used to study the position and structural role of water molecules present in the helical bundle. Finally, the binding site location and the ligand docking were investigated using an objective strategy. A suitable site for the binding of the ribose moiety of adenosine was first postulated and further confirmed by means of a novel chemometric strategy based on GRIND descriptors. Using this position as an anchor point, the binding of adenosine was studied by docking and molecular dynamics simulations obtaining two putative binding positions in good agreement with experimental data.

Depicting a protein's two faces: GPCR classification by phylogenetic tree-based HMMs

Related proteins with similar biological functions generally share common features, allowing us to extract the common sequence features. These common features enable us to build statistical models that can be used to classify proteins, to predict new members, and to study the sequence-function relationship of this protein function group. Although evolution underlies the basis of multiple sequence analysis methods, most methods ignore phylogenetic relationships and the evolutionary process in building these statistical models. Previously we have shown that a phylogenetic tree-based profile hidden Markov model (T-HMM) is superior in generating a profile for a group of similar proteins. In this study we used the method to generate common features of G protein-coupled receptors (GPCRs). The profile generated by T-HMM gives high accuracy in GPCR function classification, both by ligand and by coupled G protein.

CASTOR: clustering algorithm for sequence taxonomical organization and relationships

Given a set of related proteins, two important problems in biology are the inference of protein subsets such that members of one subset share a common function and the identification of protein regions that possess functional significance. The former is typically approached by hierarchical bottom-up clustering based on pairwise sequence similarity and various linkage rules. The latter is typically approached in a supervised manner, based on global multiple sequence alignment. However, the two problems are inextricably linked, since functional subsets are usually characterized by distinctive functional regions. This paper introduces CASTOR, an automatic and unsupervised system that addresses both problems simultaneously and efficiently. It identifies protein regions that are likely to have functional significance by discovering and refining statistically significant motifs. It infers likely functional protein subsets and their relationships based on the presence of the discovered motifs in a top-down and recursive manner, allowing the identification of both hierarchical and nonhierarchical subset relationships. This is, to our knowledge, the first system that approaches both problems simultaneously in a top-down, systematic manner. CASTOR's performance is evaluated against the G-protein coupled receptor superfamily. The identified protein regions lead to a taxonomical organization of this superfamily that is in remarkable agreement with a biologically motivated one and which outperforms those produced by bottom-up clustering methods. We also find that conventional hierarchical representations may fail to accurately describe the complexity of evolutionary development responsible for the final organization of a complex protein family. In particular, many functional relationships governing distant subfamilies of such a protein family may not be represented hierarchically.

From subgenome analysis to protein structure

Groups of related genes abound in large eukaryotic genomes. In such 'subgenomes', homology modeling carried out for a few genes will probably have relevance to the entire group. Subgenomes also afford unique ways of determining protein structural information. In addition to analyses based on the quantification of residue variability in paralogs, two-way comparisons, both within and among species, help to disclose functional amino acids. Comparative studies of gene families throughout the mammalian genome will also help elucidate the functional significance of single nucleotide polymorphisms in coding regions.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

3D model for TM region of the AT-1 receptor in complex with angiotensin II independently validated by site-directed mutagenesis data

A three-dimensional model of the complex of angiotensin II (AII) with the transmembrane (TM) region of the angiotensin II receptor of type 1 (the AT-1 receptor) was obtained by molecular modeling procedures employing structural homology to the X-ray structure of rhodopsin. Since the modeling procedure considered only steric and energy considerations without prior knowledge of the experimental results of site-directed mutagenesis, the results with receptor mutants could be used for independent validation of the model. Indeed, the model brings in contact the residues of AII responsible for agonistic activity, Tyr(4), His(6), and Phe(8), with many residues of AT-1 involved in signal transduction according to site-directed mutagenesis. The model also predicts the existence of several possible conformational pathways for transferring the binding signal through the TM region of AT-1 to the intracellular loops interacting with the G-protein.