A low resolution model for the interaction of G proteins with G protein-coupled receptors

Proteome-wide analysis reveals G protein-coupled receptor-like proteins in rice (Oryza sativa)

G protein-coupled receptors (GPCRs) constitute the largest family of transmembrane proteins in metazoans that mediate the regulation of various physiological responses to discrete ligands through heterotrimeric G protein subunits. The existence of GPCRs in plant is contentious, but their comparable crucial role in various signaling pathways necessitates the identification of novel remote GPCR-like proteins that essentially interact with the plant G protein α subunit and facilitate the transduction of various stimuli. In this study, we identified three putative GPCR-like proteins (OsGPCRLPs) (LOC_Os06g09930.1, LOC_Os04g36630.1, and LOC_Os01g54784.1) in the rice proteome using a stringent bioinformatics workflow. The identified OsGPCRLPs exhibited a canonical GPCR 'type I' 7TM topology, patterns, and biologically significant sites for membrane anchorage and desensitization. Cluster-based interactome mapping revealed that the identified proteins interact with the G protein α subunit which is a characteristic feature of GPCRs. Computational results showing the interaction of identified GPCR-like proteins with G protein α subunit and its further validation by the membrane yeast-two-hybrid assay strongly suggest the presence of GPCR-like 7TM proteins in the rice proteome. The absence of a regulator of G protein signaling (RGS) box in the C- terminal domain, and the presence of signature motifs of canonical GPCR in the identified OsGPCRLPs strongly suggest that the rice proteome contains GPCR-like proteins that might be involved in signal transduction.

Identification and CRISPR-Cas9 validation of a novel β-adrenergic-like octopamine receptor mutation associated with amitraz resistance in Varroa destructor

Varroa destructor is widely recognized as a significant contributor to colony collapse disorder. Chemical acaricides, such as amitraz, have been extensively used for Varroa control due to their selectivity within beehives. However, the increasing number of cases of amitraz resistance across global V. destructor populations poses a significant challenge. In this study, we conducted a comprehensive molecular screening of the β-adrenergic-like octopamine receptor (Octβ2R), the target-site of amitraz, across 66 Turkish and 63 Belgian V. destructor populations. Although previously reported amitraz resistance mutations were not detected, the screening revealed a novel Y337F mutation located within transmembrane 7 (TM7) of Octβ2R in Turkish Varroa populations. Notably, this mutation was identified in the last residue of the highly conserved NPxxY motif associated with the activation of G-protein coupled receptors (GPCR). Among the 66 Varroa samples from Türkiye, twenty harbored the Y337F mutation, with eight samples exhibiting fixation of the mutation. Subsequent bioassays revealed over 8-fold resistance to amitraz in populations that contain the Y337F mutation. Genotyping of mites after exposure to 10 mg a.i./l amitraz demonstrated that all surviving mites were homozygous for the Y337F mutation, whereas dead mites carried susceptible alleles, providing genetic linkage between mutation and phenotype. Further, we used CRISPR-Cas9 editing to introduce the Y337F mutation in the orthologous Octβ2R of the model organism Tetranychus urticae. Crispants exhibited over threefold resistance to amitraz. In conclusion, this study identified and validated a novel amitraz resistance mutation. Additional research is required to further evaluate the phenotypic strength of Y337F in the context of operational resistance with current treatment strategies.

GPCR-BERT: Interpreting Sequential Design of G Protein-Coupled Receptors Using Protein Language Models

With the rise of transformers and large language models (LLMs) in chemistry and biology, new avenues for the design and understanding of therapeutics have been opened up to the scientific community. Protein sequences can be modeled as language and can take advantage of recent advances in LLMs, specifically with the abundance of our access to the protein sequence data sets. In this letter, we developed the GPCR-BERT model for understanding the sequential design of G protein-coupled receptors (GPCRs). GPCRs are the target of over one-third of Food and Drug Administration-approved pharmaceuticals. However, there is a lack of comprehensive understanding regarding the relationship among amino acid sequence, ligand selectivity, and conformational motifs (such as NPxxY, CWxP, and E/DRY). By utilizing the pretrained protein model (Prot-Bert) and fine-tuning with prediction tasks of variations in the motifs, we were able to shed light on several relationships between residues in the binding pocket and some of the conserved motifs. To achieve this, we took advantage of attention weights and hidden states of the model that are interpreted to extract the extent of contributions of amino acids in dictating the type of masked ones. The fine-tuned models demonstrated high accuracy in predicting hidden residues within the motifs. In addition, the analysis of embedding was performed over 3D structures to elucidate the higher-order interactions within the conformations of the receptors.

Melatonin receptor, GhCAND2-D5 motivated responding to NaCl signaling in cotton

As a receptor for plant melatonin, CAND2/PMTR plays an important role in melatonin signaling. Most of the CANDs are membrane proteins and play indispensable roles in signal transduction. In this study, the CANDs from four cotton species were characterized, and the phylogenetic relationships, expression patterns, stress responses of cotton CANDs were analyzed by bioinformatics. Through the analysis of phylogenetic and protein structure, it was found that the CANDs in clade Ⅱ might function as cotton melatonin receptors, and most of the GhCANDs in clade Ⅱ were induced by melatonin. A putative cotton melatonin receptor, GhCAND2-D5, was functionally probed by gene silencing. The plants with silenced expression of this gene exhibited decreased salt tolerance. Protein interaction prediction identified that GhCAND2-D5 interacted with several membrane proteins and played an important role in melatonin signaling. This study provided a theoretical reference for further investigation of melatonin signaling in cotton.

Life, death and resurrection of plant GPCRs

The activation of G-protein coupled receptors (GPCRs) by extracellular ligands constitutes the first step of heterotrimeric G-protein signalling in animals. In plants, canonical GPCRs have been known for over 25 years, often in association with agronomically important functions. But their role in plant G-protein signalling and even their annotation as GPCR was contested in the last decade, only to be revisited in the light of more recent evidences. In this first ever review on plant GPCRs, we catalogue all the plant GPCRs described to date and discuss the evidences for and against their role in plants in general and G-protein signalling in particular. We argue against writing off GPCRs and point to the missing links to be investigated to establish firm conclusions either way.

Research Advances in Heterotrimeric G-Protein α Subunits and Uncanonical G-Protein Coupled Receptors in Plants

As crucial signal transducers, G-proteins and G-protein-coupled receptors (GPCRs) have attracted increasing attention in the field of signal transduction. Research on G-proteins and GPCRs has mainly focused on animals, while research on plants is relatively rare. The mode of action of G-proteins is quite different from that in animals. The G-protein α (Gα) subunit is the most essential member of the G-protein signal cycle in animals and plants. The G-protein is activated when Gα releases GDP and binds to GTP, and the relationships with the GPCR and the downstream signal are also achieved by Gα coupling. It is important to study the role of Gα in the signaling pathway to explore the regulatory mechanism of G-proteins. The existence of a self-activated Gα in plants makes it unnecessary for the canonical GPCR to activate the G-protein by exchanging GDP with GTP. However, putative GPCRs have been found and proven to play important roles in G-protein signal transduction. The unique mode of action of G-proteins and the function of putative GPCRs in plants suggest that the same definition used in animal research cannot be used to study uncanonical GPCRs in plants. This review focuses on the different functions of the Gα and the mode of action between plants and animals as well as the functions of the uncanonical GPCR. This review employs a new perspective to define uncanonical GPCRs in plants and emphasizes the role of uncanonical GPCRs and Gα subunits in plant stress resistance and agricultural production.

Motions around conserved helical weak spots facilitate GPCR activation

G protein‐coupled receptors (GPCRs) participate in most physiological processes and are important drug targets in many therapeutic areas. Recently, many GPCR X‐ray structures became available, facilitating detailed studies of their sequence‐structure‐mobility‐function relations. We show that the functional role of many conserved GPCR sequence motifs is to create weak spots in the transmembrane helices that provide the structural plasticity necessary for ligand binding and signaling. Different receptor families use different conserved sequence motifs to obtain similar helix irregularities that allow for the same motions upon GPCR activation. These conserved motions come together to facilitate the timely release of the conserved sodium ion to the cytosol. Most GPCR crystal structures could be determined only after stabilization of the transmembrane helices by mutations that remove weak spots. These mutations often lead to diminished binding of agonists, but not antagonists, which logically agrees with the fact that large helix rearrangements occur only upon agonist binding. Upon activation, six of the seven TM helices in GPCRs undergo helix motions and/or deformations facilitated by weak spots in these helices. The location of these weak spots is much more conserved than the sequence motifs that cause them. Knowledge about these weak spots helps understand the activation process of GPCRs and thus helps design medicines.

Adaptations to High Salt in a Halophilic Protist: Differential Expression and Gene Acquisitions through Duplications and Gene Transfers

The capacity of halophiles to thrive in extreme hypersaline habitats derives partly from the tight regulation of ion homeostasis, the salt-dependent adjustment of plasma membrane fluidity, and the increased capability to manage oxidative stress. Halophilic bacteria, and archaea have been intensively studied, and substantial research has been conducted on halophilic fungi, and the green alga Dunaliella. By contrast, there have been very few investigations of halophiles that are phagotrophic protists, i.e., protozoa. To gather fundamental knowledge about salt adaptation in these organisms, we studied the transcriptome-level response of Halocafeteria seosinensis (Stramenopiles) grown under contrasting salinities. We provided further evolutionary context to our analysis by identifying genes that underwent recent duplications. Genes that were highly responsive to salinity variations were involved in stress response (e.g., chaperones), ion homeostasis (e.g., Na+/H+ transporter), metabolism and transport of lipids (e.g., sterol biosynthetic genes), carbohydrate metabolism (e.g., glycosidases), and signal transduction pathways (e.g., transcription factors). A significantly high proportion (43%) of duplicated genes were also differentially expressed, accentuating the importance of gene expansion in adaptation by H. seosinensis to high salt environments. Furthermore, we found two genes that were lateral acquisitions from bacteria, and were also highly up-regulated and highly expressed at high salt, suggesting that this evolutionary mechanism could also have facilitated adaptation to high salt. We propose that a transition toward high-salt adaptation in the ancestors of H. seosinensis required the acquisition of new genes via duplication, and some lateral gene transfers (LGTs), as well as the alteration of transcriptional programs, leading to increased stress resistance, proper establishment of ion gradients, and modification of cell structure properties like membrane fluidity.

Intrinsic Relative Stabilities of the Neutral Tautomers of Arginine Side-Chain Models

The specific protonation state of amino acids is crucial for the physicochemical properties of proteins and their biological functions. These protonation states influence, for instance, properties related to hydrogen bonding, solubility, and folding. pKa calculations for proteins are, therefore, important and require, in principle, a specification of the most stable protonated and deprotonated forms of each titratable group. This is complicated by the existence of multiple tautomers, like the five neutral tautomers of the guanidine moiety in arginine. In this study, the compounds N-methyl-guanidine and N-ethyl-guanidine were used to model the charged and all neutral protonation states of the arginine side chain. The relative stabilities of all five neutral tautomers were investigated systematically for the first time, using quantum-mechanical calculations. These relative stabilities were obtained in vacuo, water and chloroform, by combining the quantum-mechanical calculations with a continuum solvation model. The water model was used to represent arginines exposed to an aqueous solution, whereas the chloroform model has a polarity representative of a protein core or a membrane. This allowed determining the relative pKa's associated with each neutral tautomer in these environments. A key result is that significant differences in stability are found between the neutral tautomers, in both water and chloroform. Some tautomers are consistently found to be the most stable. These findings will be helpful to refine pKa calculations in proteins.

Generic GPCR residue numbers - aligning topology maps while minding the gaps

Generic residue numbers facilitate comparisons of, for example, mutational effects, ligand interactions, and structural motifs. The numbering scheme by Ballesteros and Weinstein for residues within the class A GPCRs (G protein-coupled receptors) has more than 1100 citations, and the recent crystal structures for classes B, C, and F now call for a community consensus in residue numbering within and across these classes. Furthermore, the structural era has uncovered helix bulges and constrictions that offset the generic residue numbers. The use of generic residue numbers depends on convenient access by pharmacologists, chemists, and structural biologists. We review the generic residue numbering schemes for each GPCR class, as well as a complementary structure-based scheme, and provide illustrative examples and GPCR database (GPCRDB) web tools to number any receptor sequence or structure.

Alpha-bulges in G protein-coupled receptors

Agonist binding is related to a series of motions in G protein-coupled receptors (GPCRs) that result in the separation of transmembrane helices III and VI at their cytosolic ends and subsequent G protein binding. A large number of smaller motions also seem to be associated with activation. Most helices in GPCRs are highly irregular and often contain kinks, with extensive literature already available about the role of prolines in kink formation and the precise function of these kinks. GPCR transmembrane helices also contain many α-bulges. In this article we aim to draw attention to the role of these α-bulges in ligand and G-protein binding, as well as their role in several aspects of the mobility associated with GPCR activation. This mobility includes regularization and translation of helix III in the extracellular direction, a rotation of the entire helix VI, an inward movement of the helices near the extracellular side, and a concerted motion of the cytosolic ends of the helices that makes their orientation appear more circular and that opens up space for the G protein to bind. In several cases, α-bulges either appear or disappear as part of the activation process.

Design and functional evaluation of an optically active μ-opioid receptor

The use of opioids, which achieve therapeutic analgesia through activation of μ-opioid receptors, are limited in the management of chronic pain by adverse effects including tolerance and addiction. Optogenetics is an emerging approach of designing molecular targets that can produce cell-specific receptor-mediated analgesia with minimal side effects. Here we report the design and functional characterization of a chimeric μ-opioid receptor that could be photoactivated to trigger intracellular signaling. A prototype optoactive μ-opioid receptor (optoMOR) was designed by replacing the intracellular domains from rhodopsin with those of the native μ-opioid receptor and was transiently expressed in human embryonic kidney (HEK293) cells. Expression and distribution of the protein were confirmed by immunocytochemistry. The signal-transduction mechanisms induced by photoactivation of the optoMOR were evaluated and compared with the native μ-opioid receptor stimulation by an agonist, D-Ala(2), N-MePhe(4), Gly-ol-enkephalin (DAMGO). Cells were depolarized by extracellular potassium and the depolarization-induced calcium (Ca(2+)) influx was quantified by using Fura-2 imaging. The forskolin-stimulated adenylate cyclase/cAMP cascade was evaluated by ELISA or western blotting of brain-derived neurotrophic factor (BDNF) and the phosphorylation of cAMP response element binding protein (CREB). The optoMOR protein distribution was observed intracellularly and on the plasma membrane similar to the native μ-opioid receptor in HEK293 cells. Photoactivation of optoMOR decreased the Ca(2+) influx and inhibited the forskolin-induced cAMP generation, activation of CREB, and BDNF levels in optoMOR-expressing cells similar to the activation of native μ-opioid receptor by DAMGO. Thus the current study has accomplished the design of a prototype optoMOR and characterized the cellular signaling mechanisms activated by light stimulation of this receptor.

Congenital hypogonadotropic hypogonadism due to GnRH receptor mutations in three brothers reveal sites affecting conformation and coupling

Congenital hypogonadotropic hypogonadism (CHH) is characterized by low gonadotropins and failure to progress normally through puberty. Mutations in the gene encoding the GnRH receptor (GNRHR1) result in CHH when present as compound heterozygous or homozygous inactivating mutations. This study identifies and characterizes the properties of two novel GNRHR1 mutations in a family in which three brothers display normosmic CHH while their sister was unaffected. Molecular analysis in the proband and the affected brothers revealed two novel non-synonymous missense GNRHR1 mutations, present in a compound heterozygous state, whereas their unaffected parents possessed only one inactivating mutation, demonstrating the autosomal recessive transmission in this kindred and excluding X-linked inheritance equivocally suggested by the initial pedigree analysis. The first mutation at c.845 C>G introduces an Arg substitution for the conserved Pro 282 in transmembrane domain (TMD) 6. The Pro282Arg mutant is unable to bind radiolabeled GnRH analogue. As this conserved residue is important in receptor conformation, it is likely that the mutation perturbs the binding pocket and affects trafficking to the cell surface. The second mutation at c.968 A>G introduces a Cys substitution for Tyr 323 in the functionally crucial N/DPxxY motif in TMD 7. The Tyr323Cys mutant has an increased GnRH binding affinity but reduced receptor expression at the plasma membrane and impaired G protein-coupling. Inositol phosphate accumulation assays demonstrated absent and impaired Gα_q/11 signal transduction by Pro282Arg and Tyr323Cys mutants, respectively. Pretreatment with the membrane permeant GnRHR antagonist NBI-42902, which rescues cell surface expression of many GNRHR1 mutants, significantly increased the levels of radioligand binding and intracellular signaling of the Tyr323Cys mutant but not Pro282Arg. Immunocytochemistry confirmed that both mutants are present on the cell membrane albeit at low levels. Together these molecular deficiencies of the two novel GNRHR1 mutations lead to the CHH phenotype when present as a compound heterozygote.

Drug design for ever, from hype to hope

In its first 25 years JCAMD has been disseminating a large number of techniques aimed at finding better medicines faster. These include genetic algorithms, COMFA, QSAR, structure based techniques, homology modelling, high throughput screening, combichem, and dozens more that were a hype in their time and that now are just a useful addition to the drug-designers toolbox. Despite massive efforts throughout academic and industrial drug design research departments, the number of FDA-approved new molecular entities per year stagnates, and the pharmaceutical industry is reorganising accordingly. The recent spate of industrial consolidations and the concomitant move towards outsourcing of research activities requires better integration of all activities along the chain from bench to bedside. The next 25 years will undoubtedly show a series of translational science activities that are aimed at a better communication between all parties involved, from quantum chemistry to bedside and from academia to industry. This will above all include understanding the underlying biological problem and optimal use of all available data.

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.

Rice G-protein coupled receptor (GPCR): in silico analysis and transcription regulation under abiotic stress

Majority of transmembrane signal transduction in response to diverse external stimuli is mediated by G-protein coupled receptors (GPCRs) and are the principal signal transducers. GPCRs are characterized by seven membrane-spanning domains with an extracellular N-terminus and a cytoplasmic C-terminus which functions along with GTP-binding protein in a highly coordinated fashion. Role of heterotrimeric G-proteins in abiotic stresses has been reported, but the response of GPCR is not yet well characterized. In the present study we report the isolation of one putative GPCR (966 bp) from Indica rice (Oryza sativa cv. Indica group Swarna) and described its transcriptional regulation under abiotic stresses. Amino acid sequence analyses shows the presence of typical heptahelical transmembrane spanning domains with extracellular N-terminus involved in ligand binding and cytoplasm facing C-terminus that binds with heterotrimeric G-protein. Sequence analysis also confirmed the presence of all signature motifs required for functional GPCR. Domain and site prediction shows the presence of myristoylation sites for membrane association, and protein kinase C sites for its desensitization. The transcript levels of rice GPCR was induced following NaCl and ABA treatments. However, in drought condition the expression profile of GPCR up regulated during early exposure which subsequently decreased. On the other hand it seems no significant effect due to cold and heat stress. These findings provide a direct evidence for transcriptional regulation of rice GPCR under abiotic stress conditions. These findings also suggest that GPCR can be exploited for promoting stress tolerance in plants.

Lipoic acid stimulates cAMP production via G protein-coupled receptor-dependent and -independent mechanisms

Lipoic acid (LA) is a naturally occurring fatty acid that exhibits anti-oxidant and anti-inflammatory properties and is being pursued as a therapeutic for many diseases including multiple sclerosis, diabetic polyneuropathy and Alzheimer's disease. We previously reported on the novel finding that racemic LA (50:50 mixture of R-LA and S-LA) stimulates cAMP production, activates prostanoid EP2 and EP4 receptors and adenylyl cyclases (AC), and suppresses activation and cytotoxicity in NK cells. In this study, we present evidence that furthers our understanding of the mechanisms of action of LA. Using various LA derivatives, such as dihydrolipoic acid (DHLA), S,S-dimethyl lipoic acid (DMLA) and lipoamide (LPM), we discovered that only LA is capable of stimulating cAMP production in NK cells. Furthermore, there is no difference in cAMP production after stimulation with either R-LA, S-LA or racemic LA. Competition and synergistic studies indicate that LA may also activate AC independent of the EP2 and EP4 receptors. Pretreatment of PBMCs with KH7 (a specific peptide inhibitor of soluble AC) and the calcium inhibitor (Bapta) prior to LA treatment resulted in reduced cAMP levels, suggesting that soluble AC and calcium signaling mediate LA stimulation of cAMP production. In addition, pharmacological inhibitor studies demonstrate that LA also activates other G protein-coupled receptors, including histamine and adenosine but not the β-adrenergic receptors. These novel findings provide information to better understand the mechanisms of action of LA, which can help facilitate the use of LA as a therapeutic for various diseases.

Cloning of an octopamine/tyramine receptor and plasticity of its expression as a function of adult sexual maturation in the male moth Agrotis ipsilon

In the male moth Agrotis ipsilon behavioural response and antennal lobe (AL) neuron sensitivity to the female-produced sex pheromone increase with age and juvenile hormone (JH) level. We recently showed that the neuromodulator, octopamine (OA), interacts with JH in this age-dependent olfactory plasticity. To further elucidate its role, we cloned a full cDNA encoding a protein that presents biochemical features essential to OA/tyramine receptor (AipsOAR/TAR) function. The AipsOAR/TAR transcript was detected predominantly in the antennae, the brain and, more specifically, in ALs where its expression level varied concomitantly with age. This expression plasticity indicates that AipsOAR/TAR might be involved in central processing of the pheromone signal during maturation of sexual behaviour in A. ipsilon.

Characterization of a novel octopamine receptor expressed in the surf clam Spisula solidissima

We have cloned and sequenced a cDNA from the surf clam (Spisula solidissima, a pelecypod mollusc) that encodes an octopamine receptor which we have named Spi-OAR. The sequence of Spi-OAR shares many similarities with two Aplysia and three Drosophila octopamine receptors belonging to a sub-group of beta-adrenergic-like octopamine receptors. Using an expression vector and transient transfections of Spi-OAR into HEK 293 cells, we observed an increase of cAMP upon addition of octopamine and, to a lesser extent, of tyramine, but not after addition of dopamine, serotonin, or histamine. Using a battery of known agonists and antagonists for octopamine receptors, we observed a rather unique pharmacological profile for Spi-OAR through measurements of cAMP. Spi-OAR exhibited some constitutive activity in HEK 293 cells and no Ca(2+) responses could be detected following addition of octopamine to Spi-OAR-transfected cells. RT-PCR analysis revealed ubiquitous expression of Spi-OAR mRNA in all adult tissues, oocytes and early embryos examined. While addition of serotonin to isolated clam oocytes resulted in meiotic activation, similar additions of octopamine had no effect, suggesting that its potential role in clam reproductive physiology differs significantly from that of serotonin. This work identifies Spi-OAR as a novel mollusc octopamine receptor closely related to other invertebrate beta-adrenergic-like octopamine receptors, with possible reproductive and other physiological functions. This initial characterization of Spi-OAR makes possible further investigations and comparisons with more studied and familiar insect or gastropod mollusc octopamine receptors.

Homology modelling and spectroscopy, a never-ending love story

Homology modelling is normally the technique of choice when experimental structure data are not available but three-dimensional coordinates are needed, for example, to aid with detailed interpretation of results of spectroscopic studies. Herein, the state of the art of homology modelling will be described in the light of a series of recent developments, and an overview will be given of the problems and opportunities encountered in this field. The major topic, the accuracy and precision of homology models, will be discussed extensively due to its influence on the reliability of conclusions drawn from the combination of homology models and spectroscopic data. Three real-world examples will illustrate how both homology modelling and spectroscopy can be beneficial for (bio)medical research.

Whole proteome identification of plant candidate G-protein coupled receptors in Arabidopsis, rice, and poplar: computational prediction and in-vivo protein coupling

Computational prediction and in vivo protein coupling experiments identify candidate plant G-protein coupled receptors in Arabidopsis, rice and poplar.

Background

The classic paradigm of heterotrimeric G-protein signaling describes a heptahelical, membrane-spanning G-protein coupled receptor that physically interacts with an intracellular Gα subunit of the G-protein heterotrimer to transduce signals. G-protein coupled receptors comprise the largest protein superfamily in metazoa and are physiologically important as they sense highly diverse stimuli and play key roles in human disease. The heterotrimeric G-protein signaling mechanism is conserved across metazoa, and also readily identifiable in plants, but the low sequence conservation of G-protein coupled receptors hampers the identification of novel ones. Using diverse computational methods, we performed whole-proteome analyses of the three dominant model plant species, the herbaceous dicot Arabidopsis thaliana (mouse-eared cress), the monocot Oryza sativa (rice), and the woody dicot Populus trichocarpa (poplar), to identify plant protein sequences most likely to be GPCRs.

Results

Our stringent bioinformatic pipeline allowed the high confidence identification of candidate G-protein coupled receptors within the Arabidopsis, Oryza, and Populus proteomes. We extended these computational results through actual wet-bench experiments where we tested over half of our highest ranking Arabidopsis candidate G-protein coupled receptors for the ability to physically couple with GPA1, the sole Gα in Arabidopsis. We found that seven out of eight tested candidate G-protein coupled receptors do in fact interact with GPA1. We show through G-protein coupled receptor classification and molecular evolutionary analyses that both individual G-protein coupled receptor candidates and candidate G-protein coupled receptor families are conserved across plant species and that, in some cases, this conservation extends to metazoans.

Conclusion

Our computational and wet-bench results provide the first step toward understanding the diversity, conservation, and functional roles of plant candidate G-protein coupled receptors.

How do receptors activate G proteins?

Heterotrimeric G proteins couple the activation of heptahelical receptors at the cell surface to the intracellular signaling cascades that mediate the physiological responses to extracellular stimuli. G proteins are molecular switches that are activated by receptor-catalyzed GTP for GDP exchange on the G protein alpha subunit, which is the rate-limiting step in the activation of all downstream signaling. Despite the important biological role of the receptor-G protein interaction, relatively little is known about the structure of the complex and how it leads to nucleotide exchange. This chapter will describe what is known about receptor and G protein structure and outline a strategy for assembling the current data into improved models for the receptor-G protein complex that will hopefully answer the question as to how receptors flip the G protein switch.

Endosomal colocalization of melanocortin-3 receptor and beta-arrestins in CAD cells with altered modification of AKT/PKB

The melanocortin 3-receptor is involved in regulating energy metabolism, body fluid composition and inflammatory responses. Melanocortin receptors function by activating membrane bound adenylate cyclase. However, the literature reports indicate that some G protein coupled receptors (GPCRs) can also activate mitogen activated protein kinase (MAPK) or phosphoinositide 3 kinase (PI3K) signaling pathways consequent to their endocytosis. These studies were undertaken to evaluate the role of these pathways in MC3R signaling in brain-stem neuronal cells. Recruitment of arrestins is implicated in the activation of secondary pathways by GPCRs and our data shows the colocalization of either arrestin B1 or B2 with MC3R in endosomes. An alteration in PKB phosphorylation pattern was observed in MC3R expressing cells independent of agonist stimulation. MC3R transfectants exhibited increased proliferation rates and inhibition of PKB pathway with triciribine abrogated cell proliferation in both vector control and MC3R transfectants. PKB is constitutively active in proliferating CAD cells but could be further activated by culturing the cells in differentiation medium. These studies suggest that the AKT/PKB pathway plays an important role in the proliferation of CAD cells and suggest a link between MC3R and cell growth pathways that may involve the alteration of AKT/PKB signaling pathway.

Computational studies of CXCR1, the receptor of IL-8/CXCL8, using molecular dynamics and electrostatics

The three-dimensional structure of IL-8/CXCL8 has been previously determined using NMR spectroscopy and X-ray crystallography, but the structure of the receptors for this chemokine has not been determined experimentally. We present here the development of a model for the structure of the IL-8/CXCL8 receptor CXCR1, using a combination of homology modeling and a molecular dynamics simulation. Based on this model, we discuss the analysis of structural, dynamic, and physicochemical properties of CXCR1. We focused on the role of pairwise ionic interactions in local structural stability of CXCR1 and the role of electrostatic potentials in recognition of CXCR1 with IL-8/CXCL8. We have performed theoretical mutations of six charged amino acids in CXCR1, which abolish binding as suggested by earlier experimental data, to shed light on the effect of charge on association ability. We propose that the observed loss of binding in the six CXCR1 mutants is owed to loss of local structural stability, rather than hindrance of the recognition process because of changes in the overall electrostatic properties of the receptor. Based on further structural analysis, we propose some mutations of charged residues involving ion pairs in different elements of transmembrane helices and extracellular loops, which are expected to alter the local structure and possibly affect binding.

Emerging concepts of guanine nucleotide-binding protein-coupled receptor (GPCR) function and implications for high throughput screening

Guanine nucleotide binding protein (G protein) coupled receptors (GPCRs) comprise one of the largest families of proteins in the human genome and are a target for 40% of all approved drugs. GPCRs have unique structural motifs that allow them to interact with a wide and diverse series of extracellular ligands, as well as intracellular proteins, G proteins, receptor activity-modifying proteins, arrestins, and indeed other receptors. This distinctive structure has led to numerous efforts to discover drugs against GPCRs with targeted therapeutic uses. Such "designer" drugs currently include allosteric regulators, inverse agonists, and drugs targeting hetero-oligomeric complexes. Moreover, the large family of orphan GPCRs provides a rich and novel field of targets to discover drugs with unique therapeutic properties. The numerous technologies to discover GPCR drugs have also greatly advanced over the years, facilitating compound screening against known and orphan GPCRs, as well as in the identification of unique designer GPCR drugs. Indeed, high throughput screening (HTS) technologies employing functional cell-based approaches are now widely used. These include measurement of second messenger accumulation such as cyclic AMP, calcium ions, and inositol phosphates, as well as mitogen-activated protein kinase activation, protein-protein interactions, and GPCR oligomerization. This review focuses on how the improved understanding of the molecular pharmacology of GPCRs, coupled with a plethora of novel HTS technologies, is leading to the discovery and development of an entirely new generation of GPCR-based therapeutics.

The Angiotensin II AT1 Receptor Structure-Activity Correlations in the Light of Rhodopsin Structure

The most prevalent physiological effects of ANG II, the main product of the renin-angiotensin system, are mediated by the AT1 receptor, a rhodopsin-like AGPCR. Numerous studies of the cardiovascular effects of synthetic peptide analogs allowed a detailed mapping of ANG II's structural requirements for receptor binding and activation, which were complemented by site-directed mutagenesis studies on the AT1 receptor to investigate the role of its structure in ligand binding, signal transduction, phosphorylation, binding to arrestins, internalization, desensitization, tachyphylaxis, and other properties. The knowledge of the high-resolution structure of rhodopsin allowed homology modeling of the AT1 receptor. The models thus built and mutagenesis data indicate that physiological (agonist binding) or constitutive (mutated receptor) activation may involve different degrees of expansion of the receptor's central cavity. Residues in ANG II structure seem to control these conformational changes and to dictate the type of cytosolic event elicited during the activation. 1) Agonist aromatic residues (Phe8 and Tyr4) favor the coupling to G protein, and 2) absence of these residues can favor a mechanism leading directly to receptor internalization via phosphorylation by specific kinases of the receptor's COOH-terminal Ser and Thr residues, arrestin binding, and clathrin-dependent coated-pit vesicles. On the other hand, the NH2-terminal residues of the agonists ANG II and [Sar1]-ANG II were found to bind by two distinct modes to the AT1 receptor extracellular site flanked by the COOH-terminal segments of the EC-3 loop and the NH2-terminal domain. Since the [Sar1]-ligand is the most potent molecule to trigger tachyphylaxis in AT1 receptors, it was suggested that its corresponding binding mode might be associated with this special condition of receptors.

Modeling GPCRs

Many GPCR models have been built over the years for many different purposes, of which drug-design undoubtedly has been the most frequent one. The release of the structure of bovine rhodopsin in August 2000 enabled us to analyze models built before that period to learn things for the models we build today. We conclude that the GPCR modeling field is riddled with "common knowledge". Several characteristics of the bovine rhodopsin structure came as a big surprise, and had obviously not been predicted, which led to large errors in the models. Some of these surprises, however, could have been predicted if the modelers had more rigidly stuck to the rule that holds for all models, namely that a model should explain all experimental facts, and not just those facts that agree with the modeler's preconceptions.

Mechanism of the receptor-catalyzed activation of heterotrimeric G proteins

Heptahelical receptors activate intracellular signaling pathways by catalyzing GTP for GDP exchange on the heterotrimeric G protein alpha subunit (G alpha). Despite the crucial role of this process in cell signaling, little is known about the mechanism of G protein activation. Here we explore the structural basis for receptor-mediated GDP release using electron paramagnetic resonance spectroscopy. Binding to the activated receptor (R*) causes an apparent rigid-body movement of the alpha5 helix of G alpha that would perturb GDP binding at the beta6-alpha5 loop. This movement was not observed when a flexible loop was inserted between the alpha5 helix and the R*-binding C terminus, which uncouples R* binding from nucleotide exchange, suggesting that this movement is necessary for GDP release. These data provide the first direct observation of R*-mediated conformational changes in G proteins and define the structural basis for GDP release from G alpha.

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.

Developmental expression of a tyramine receptor gene in the brain of the honey bee, Apis mellifera

This study reveals that the tyramine receptor gene, Amtyr1, is expressed in the developing brain, as well as in the brain of the adult worker honey bee. Changes in levels of Amtyr1 expression were examined using Northern analysis. Age-related increases in Amtyr1 transcript levels were observed not only during metamorphic adult development, but also in the brain of the adult worker bee. RNA in situ hybridization revealed the pattern of Amtyr1 expression. Cell bodies staining intensely for tyramine receptor-gene transcript were observed throughout the somata rind, with well-defined clusters of cells associated with developing mushroom bodies, optic lobes, and antennal lobes of the brain. Staining for Amtyr1 transcript was particularly intense within the three major divisions of mushroom body intrinsic neurons (outer compact, noncompact, and inner compact cells), suggesting that Amtyr1 is highly expressed in these structures. Activation of AmTYR1 receptors heterologously expressed in insect (Spodoptera frugiperda) cells led to a reduction in intracellular levels of cAMP similar to that reported for AmTYR1 receptors expressed in mammalian (HEK 293) cells (Blenau et al. [2000] J Neurochem 74:900-908). Taken together, these results suggest that AmTYR1 receptors may play a role in the developing brain as well as in the brain of the adult worker bee. The actions of tyramine are likely to be mediated, at least in part, via the cAMP-signaling pathway.

Isolation, cloning, and tissue expression of a putative octopamine/tyramine receptor from locust visceral muscle tissues

Octopamine has been shown to play major roles in invertebrate nervous systems as a neurotransmitter, neuromodulator, and neurohormone. Tyramine is the biochemical precursor of octopamine and its neuromodulatory role is now being investigated and clarified in invertebrates, particularly in insects. Both octopamine and tyramine mediate their actions via G protein-coupled receptors (GPCRs) and are believed to play important functions in the regulation of physiological processes in locust oviduct. Here we report the isolation, cloning, and tissue expression of a putative octopamine/tyramine receptor from the locust, Locusta migratoria. Degenerate oligonucleotides in PCR reactions were first used to obtain partial cDNA sequences and then these partial sequences were used in screens to obtain a full-length cDNA. The cloned cDNA is about 3.1 kb long and encodes a protein of 484 amino acid residues with typical characteristics of GPCRs including seven transmembrane domains and many signature residues. The amino acid sequence of the cloned cDNA displays sequence similarities with known GPCRs, particularly octopamine/tyramine receptors. Screening of the locust genomic DNA library resulted in isolation of a genomic DNA with the same size as the cDNA, indicating that the gene is intron-less. RT-PCR and Northern blot analyses revealed the expression of the receptor mRNA in brain, ventral nerve cord, oviduct, and midgut tissues. Southern blot analyses using EcoRI and HindIII restriction endonucleases recognized at least two distinct gene bands.

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.

Mutation of the DRY motif reveals different structural requirements for the CC chemokine receptor 5-mediated signaling and receptor endocytosis

CC chemokine receptor 5 (CCR5) is a G protein-coupled receptor that governs migration of leukocytes and serves as a coreceptor for the R5 tropic strains of human immunodeficiency virus (HIV). CCR5-mediated signaling in response to CC chemokines relies on G protein activation. Desensitization, which rapidly turns off G protein-dependent signaling, involves phosphorylation of CCR5 that promotes interaction of the receptor with beta-arrestins for endocytosis. Whether coupling to G proteins, desensitization, and endocytosis of CCR5 require the same structural determinants remains a matter of investigation. Here, we show that CCR5 displayed agonist-independent coupling to G proteins. This constitutive activity of the receptor was abrogated by TAK779 (N,N-dimethyl-N-[4-[[[2-(4-methylphenyl)-6,7-dihydro-5H-benzocyclohepten-8-yl]carbonyl]amino]benzyl]tetrahydro-2H-pyran-4-aminium chloride), a nonpeptidic CCR5 ligand that inhibits HIV infection and was found to depend on the integrity of the Asp-Arg-Tyr (DRY) motif. Changing Arg-126 by the neutral residue Asn (R126N-CCR5 mutant) abolished CCR5-mediated activation of G proteins, either constitutively or in response to agonists. In contrast, R126N-CCR5 not only retained agonist-promoted phosphorylation and beta-arrestin-dependent endocytosis but also displayed a higher basal phosphorylation than wild-type CCR5. Expression of beta-arrestin in R126N-CCR5-expressing cells resulted in receptor down-regulation, thereby suggesting that R126N-CCR5 spontaneously interacts with beta-arrestins. However, although expression of beta-arrestin favored wild-type CCR5-mediated chemotaxis, it failed to promote migration of cells expressing R126N-CCR5. Overall, these data indicate that structural requirements for CCR5-mediated activation of G proteins, albeit not involved in receptor desensitization and internalization, are needed for beta-arrestin-mediated chemotaxis. These results have implications for how distinct biological responses of CCR5 might rely on a different set of receptor conformations.

Identification of a novel aminergic-like G protein-coupled receptor in the cnidarian Renilla koellikeri

Biogenic amines exert various physiological effects in cnidarians, but the receptors involved in these responses are not known. We have cloned a novel G protein-coupled receptor cDNA from an anthozoan, the sea pansy Renilla koellikeri, that shows homology to mammalian catecholamine receptors and, to a lesser extent, to peptidergic receptors. This putative receptor, named Ren2, has a DRC pattern that replaces the well-conserved DRY motif on the cytoplasmic side of the transmembrane III and lacks the cysteine residues usually found in the second extracellular loop and C-terminus tail. Both the second extracellular loop and the N-terminal tail were seen to be short (six and three amino acids, respectively). Northern blot analysis suggests that the receptor gene codes for two transcripts. Localization of these transcripts by in situ hybridization demonstrated abundant expression in the epithelium of the pharyngeal wall, the oral disk and tentacles as well as in the endodermal epithelium lining the gastrovascular cavities.

Gαs proteinC-terminal α-helix at the interface: does the plasma membrane play a critical role in the Gαs protein functionality?

The heterotrimeric guanine nucleotide-binding regulatory proteins (G proteins, Galphabetagamma) mediate the signalling process of a large number of receptors, known as G protein-coupled receptors. The C-terminal domain of the heterotrimeric G protein alpha-subunit plays a key role in the selective activation of G proteins by their cognate receptors. The interaction of this domain can take place at the end of a cascade including several successive conformational modifications. Galpha(s)(350-394) is the 45-mer peptide corresponding to the C-terminal region of the Galpha(s) subunit. In the crystal structure of the Galpha(s) subunit it encompasses the alpha4/beta6 loop, the beta6 beta-sheet segment and the alpha5 helix region. Following a previous study based on the synthesis, biological activity and conformational analysis of shorter peptides belonging to the same Galpha(s) region, Galpha(s)(350-394) was synthesized and investigated. The present study outlines the central role played by the residues involved in the alpha4/beta6 loop and beta6/alpha5 loops in the stabilization of the C-terminal Galpha(s)alpha-helix. H(2)O/(2)H(2)O exchange experiments, and NMR diffusion experiments show interesting evidence concerning the interaction between the SDS micelles and the polypeptide. These data prompt intriguing speculations on the role of the intracellular environment/cellular membrane interface in the stabilization and functionality of the C-terminal Galpha(s) region.

Modulation of ligand selectivity associated with activation of the transmembrane region of the human follitropin receptor

Recently, three naturally occurring mutations in the serpentine region of the FSH receptor (FSHr) (D567N and T449I/A) have been identified in three families with spontaneous ovarian hyperstimulation syndrome (OHSS). All mutant receptors displayed abnormally high sensitivity to human chorionic gonadotropin and, in addition, D567N and T449A displayed concomitant increase in sensitivity to TSH and detectable constitutive activity. In the present study, we have used a combination of site-directed mutagenesis experiments and molecular modeling to explore the mechanisms responsible for the phenotype of the three OHSS FSHr mutants. Our results suggest that all mutations lead to weakening of interhelical locks between transmembrane helix (TM)-VI and TM-III, or TM-VI and TM-VII, which contributes to maintaining the receptor in the inactive state. They also indicate that broadening of the functional specificity of the mutant FSHr constructs is correlated to their increase in constitutive activity. This relation between basal activity and functional specificity is a characteristic of the FSHr, which is not shared by the other glycoprotein hormone receptors. It leads to the interesting suggestion that different pathways have been followed during primate evolution to avoid promiscuous stimulation of the TSHr and FSHr by human chorionic gonadotropin. In the hFSHr, specificity would be exerted both by the ectodomain and the serpentine portion.

Gonadotropin-releasing hormone receptors

GnRH and its analogs are used extensively for the treatment of hormone-dependent diseases and assisted reproductive techniques. They also have potential as novel contraceptives in men and women. A thorough delineation of the molecular mechanisms involved in ligand binding, receptor activation, and intracellular signal transduction is kernel to understanding disease processes and the development of specific interventions. Twenty-three structural variants of GnRH have been identified in protochordates and vertebrates. In many vertebrates, three GnRHs and three cognate receptors have been identified with distinct distributions and functions. In man, the hypothalamic GnRH regulates gonadotropin secretion through the pituitary GnRH type I receptor via activation of G(q). In-depth studies have identified amino acid residues in both the ligand and receptor involved in binding, receptor activation, and translation into intracellular signal transduction. Although the predominant coupling of the type I GnRH receptor in the gonadotrope is through productive G(q) stimulation, signal transduction can occur via other G proteins and potentially by G protein-independent means. The eventual selection of intracellular signaling may be specifically directed by variations in ligand structure. A second form of GnRH, GnRH II, conserved in all higher vertebrates, including man, is present in extrahypothalamic brain and many reproductive tissues. Its cognate receptor has been cloned from various vertebrate species, including New and Old World primates. The human gene homolog of this receptor, however, has a frame-shift and stop codon, and it appears that GnRH II signaling occurs through the type I GnRH receptor. There has been considerable plasticity in the use of different GnRHs, receptors, and signaling pathways for diverse functions. Delineation of the structural elements in GnRH and the receptor, which facilitate differential signaling, will contribute to the development of novel interventive GnRH analogs.

G Protein–Coupled Receptor Oligomerization

The cardiovascular system is richly endowed with G protein–coupled receptors (GPCRs), members of the largest family of plasma membrane-localized receptors. During the last 10 years, it has become increasingly clear that many, if not all, GPCRs function in oligomeric complexes, as either homo- or hetero-oligomers. This review explores the mechanistic implications of GPCR dimerization and/or oligomerization on receptor activation and interactions with G proteins. The effects of GPCR oligomerization on receptor pharmacology, GPCR-mediated signaling, and potential contributions to GPCR crosstalk will be considered in the context of receptors important in the cardiovascular system. Our evolving understanding of the structural and functional consequences of GPCR oligomerization may provide novel and more selective sites for pharmacological tuning of cardiovascular function.

A new G protein-coupled receptor from a primitive metazoan shows homology with vertebrate aminergic receptors and displays constitutive activity in mammalian cells

Biogenic amine receptors mediate wide-ranging hormonal and modulatory functions in vertebrates, but are largely unknown in primitive invertebrates. In a representative of the most basal multicellular animals possessing a nervous system, the cnidarian Renilla koellikeri, aminergic-like receptors were previously characterized pharmacologically and found to engender control of the animal's bioluminescent and peristaltic reactions. Using degenerate oligonucleotides in a RT-PCR strategy, we obtained a full-length cDNA encoding a polypeptide with typical G protein-coupled receptor (GPCR) characteristics and which displayed a significant degree of sequence similarity (up to 45%) to biogenic amine receptors, particularly dopamine and adrenergic receptors. The new receptor, named Ren1, did not resemble any one specific type of amine GPCR and thus could not be identified on the basis of sequence. Ren1 was expressed transiently and stably in cultured mammalian cells, as demonstrated by immunocytochemistry and western blotting. Functional analysis of transfected HEK293, LTK- and COS-7 cells, based on both cAMP and Ca2+ signalling assays, revealed that Ren1 was not activated by any of the known biogenic amines tested and several related metabolites. The results indicated, however, that cells stably expressing Ren1 contained, on average, an 11-fold higher level of cAMP than the controls, in the absence of agonist stimulation. The high basal cAMP levels were shown to be specific for Ren1 and to vary proportionally with the level of Ren1 expressed in the transfected cells. Taken together, the data suggested that Ren1 was expressed as a constitutively active receptor. Its identification provides a basis for examination of the early evolutionary emergence of GPCRs and their functional properties.

Sequence analysis reveals how G protein-coupled receptors transduce the signal to the G protein

Sequence entropy-variability plots based on alignments of very large numbers of sequences-can indicate the location in proteins of the main active site and modulator sites. In the previous article in this issue, we applied this observation to a series of well-studied proteins and concluded that it was possible to detect most of the residues with a known functional role. Here, we apply the method to rhodopsin-like G protein-coupled receptors. Our conclusion is that G protein binding is the main evolutionary constraint on these receptors, and that other ligands, such as agonists, act as modulators. The activation of the receptors can be described as a simple, two-step process, and the residues involved in signal transduction can be identified.

Identification of the binding sites and selectivity of sarpogrelate, a novel 5-HT2 antagonist, to human 5-HT2A, 5-HT2B and 5-HT2C receptor subtypes by molecular modeling

The aim of the present study was to investigate the binding sites interactions and the selectivity of sarpogrelate to human 5-HT(2) receptor family (5-HT(2A), 5-HT(2B) and 5-HT(2C) receptor subtypes) using molecular modeling. Rhodopsin (RH) crystal structures were used as template to build structural models of the human serotonin-2A and -2C receptors (5-HT(2A)R, 5-HT(2C)R), whereas for 5-HT(2B)R, we used our previously published three-dimensional (3D) models based on bacteriorhodopsin (BR). Sarpogrelate, a novel 5-HT(2)R antagonist, was docked to the receptors. Molecular dynamics (MD) simulations produced the strongest interaction for 5-HT(2A)R/sarpogrelate complex. Upon binding, sarpogrelate constraints aromatic residues network (Trp(3.28), Phe(5.47), Trp(6.48), Phe(6.51), Phe(6.52) in 5-HT(2A)R; Phe(3.35), Phe(6.51), Trp(7.40) in 5-HT(2B)R; Trp(3.28), Phe(3.35), Phe(5.47), Trp(6.48), Phe(6.51), Phe(6.52) in 5-HT(2C)R) in a stacked configuration, preventing activation of the receptor. The models suggest that the structural origin of the selectivity of sarpogrelate to 5-HT(2A)R vs both 5-HT(2B)R and 5-HT(2C)R comes from the following results: (1) The tight interaction between the antagonist and the transmembrane domain (TMD) 3. Asp(3.32) neutralizes the cationic head and interacts simultaneously with carboxylic group hydrogen of the antagonist molecule. (2) Due to steric hindrance, Ser(5.46) (vs Ala(5.46) in 5HT(2B) and 5HT(2C)) prevents sarpogrelate to enter deeply inside the hydrophobic core of the helix bundle and to interact with Pro(5.50). (3) The side chain of Ile(4.56) (vs Ile(4.56) in 5HT(2B)R and Val(4.56) in 5HT(2C)R) constraints sarpogrelate to adjust its position by translating toward the strongly attractive Asp(3.32). These results are in good agreement with binding affinities (pKi) of sarpogrelate for 5-HT(2) receptor family expressed in transfected cell.

Trifluoroethanol and binding to model membranes stabilize a predicted turn in a peptide corresponding to the first extracellular loop of the angiotensin II AT(1A) receptor

Homology modeling of the angiotensin II AT(1A) receptor based on rhodopsin's crystal structure has assigned the 92-100 (YRWPFGNHL) sequence of the receptor to its first extracellular loop. Solution and membrane-bound conformational properties of a peptide containing this sequence (EL1) were examined by CD, fluorescence, and (1)H-NMR. CD spectra in aqueous solution revealed an equilibrium between less organized and folded conformers. NMR spectra indicated the coexistence of trans and cis isomers of the Trp(3)-Pro(4) bond. A positive band at 226 nm in the CD spectra suggested aromatic ring stacking, modulated by EL1's ionization degree. CD spectra showed that trifluoroethanol (TFE), or binding to detergent micelles and phospholipid bilayers, shifted the equilibrium toward conformers with higher secondary structure content. Different media gave rise to spectra suggestive of different beta-turns. Chemical shift changes in the NMR spectra corroborated the stabilization of different conformations. Thus, environments of lower polarity or binding to interfaces probably favored the formation of hydrogen bonds, stabilizing beta-turns, predicted for this sequence in the whole receptor. Increases in Trp(3) fluorescence intensity and anisotropy, blue shifts of the maximum emission wavelength, and pK changes also evinced the interaction between EL1 and model membranes. Binding was seen to depend on both hydrophobic and electrostatic interactions, as well as lipid phase packing. Studies with water-soluble and membrane-bound fluorescence quenchers demonstrated that Trp(3) is located close to the water-membrane interface. The results are discussed with regard to possible implications in receptor folding and function.

Conformational analysis of the Galpha(s) protein C-terminal region

The C-terminal domain of the heterotrimeric G protein a-subunits plays a key role in selective activation of G proteins by their cognate receptors. Several C-terminal fragments of Galpha(s) (from 11 to 21 residues) were recently synthesized. The ability of these peptides to stimulate agonist binding was found to be related to their size. Galpha(s)(380-394) is a 15-mer peptide of intermediate length among those synthesized and tested that displays a biological activity surprisingly weak compared with that of the corresponding 21-mer peptide, shown to be the most active. In the present investigation, Galpha(s)(380-394) was subjected to a conformational NMR analysis in a fluorinated isotropic environment. An NMR structure, calculated on the basis of the data derived from conventional 1D and 2D homonuclear experiments, shows that the C-terminal residues of Galpha(s)(380-394) are involved in a helical arrangement whose length is comparable to that of the most active 21 -mer peptide. A comparative structural refinement of the NMR structures of Galpha(s)(380-394) and Galpha(s)(374-394)C379A was performed using molecular dynamics calculations. The results give structural elements to interpret the role played by both the backbone conformation and the side chain arrangement in determining the activity of the G protein C-terminal fragments. The orientation of the side chains allows the peptides to assume contacts crucial for the G protein/receptor interaction. In the 15-mer peptide the lack as well as the disorder of some N-terminal residues could explain the low biological activity observed.