Delineating a Ca2+ binding pocket within the venus flytrap module of the human calcium-sensing receptor

The Ca(2+)-sensing receptor (CaSR) belongs to the class III G-protein-coupled receptors (GPCRs), which include receptors for pheromones, amino acids, sweeteners, and the neurotransmitters glutamate and gamma-aminobutyric acid (GABA). These receptors are characterized by a long extracellular amino-terminal domain called a Venus flytrap module (VFTM) containing the ligand binding pocket. To elucidate the molecular determinants implicated in Ca(2+) recognition by the CaSR VFTM, we developed a homology model of the human CaSR VFTM from the x-ray structure of the metabotropic glutamate receptor type 1 (mGluR1), and a phylogenetic analysis of 14 class III GPCR VFTMs. We identified critical amino acids delineating a Ca(2+) binding pocket predicted to be adjacent to, but distinct from, a cavity reminiscent of the binding site described for amino acids in mGluRs, GABA-B receptor, and GPRC6a. Most interestingly, these Ca(2+)-contacting residues are well conserved within class III GPCR VFTMs. Our model was validated by mutational and functional analysis, including the characterization of activating and inactivating mutations affecting a single amino acid, Glu-297, located within the proposed Ca(2+) binding pocket of the CaSR and associated with autosomal dominant hypocalcemia and familial hypocalciuric hypercalcemia, respectively, genetic diseases characterized by perturbations in Ca(2+) homeostasis. Altogether, these data define a Ca(2+) binding pocket within the CaSR VFTM that may be conserved in several other class III GPCRs, thereby providing a molecular basis for extracellular Ca(2+) sensing by these receptors.     

Factor analysis for gene regulatory networks and transcription factor activity profiles

Background

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

Results

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

Conclusion

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

Identification of a gonadotropin-releasing hormone receptor orthologue in Caenorhabditis elegans

Background  

The Caenorhabditis elegans genome is known to code for at least 1149 G protein-coupled receptors (GPCRs), but the GPCR(s) critical to the regulation of reproduction in this nematode are not yet known. This study examined whether GPCRs orthologous to human gonadotropin-releasing hormone receptor (GnRHR) exist in C. elegans.     

The serendipitous origin of chordate secretin peptide family members

Background  

The secretin family is a pleotropic group of brain-gut peptides with affinity for class 2 G-protein coupled receptors (secretin family GPCRs) proposed to have emerged early in the metazoan radiation via gene or genome duplications. In human, 10 members exist and sequence and functional homologues and ligand-receptor pairs have been characterised in representatives of most vertebrate classes. Secretin-like family GPCR homologues have also been isolated in non-vertebrate genomes however their corresponding ligands have not been convincingly identified and their evolution remains enigmatic.     

The G protein-coupled receptors in the pufferfish <Emphasis Type="Italic">Takifugu rubripes</Emphasis>

Background

Guanine protein-coupled receptors (GPCRs) constitute a eukaryotic transmembrane protein family and function as “molecular switches” in the second messenger cascades and are found in all organisms between yeast and humans. They form the single, biggest drug-target family due to their versatility of action and their role in several physiological functions, being active players in detecting the presence of light, a variety of smells and tastes, amino acids, nucleotides, lipids, chemicals etc. in the environment of the cell. Comparative genomic studies on model organisms provide information on target receptors in humans and their function. The Japanese teleost Fugu has been identified as one of the smallest vertebrate genomes and a compact model to study the human genome, owing to the great similarity in its gene repertoire with that of human and other vertebrates. Thus the characterization of the GPCRs of Fugu would provide insights to the evolution of the vertebrate genome.

Results

We classified the GPCRs in the Fugu genome and our analysis of its 316 membrane-bound receptors, available on the public databases as well as from literature, detected 298 GPCRs that were grouped into five main families according to the GRAFS classification system (namely, Glutamate, Rhodopsin, Adhesion, Frizzled and Secretin). We also identified 18 other GPCRs that could not be grouped under the GRAFS family and hence were classified as ‘Other 7TM’ receptors. On comparison of the GPCR information from the Fugu genome with those in the human and chicken genomes, we detected 96.83% (306/316) and 96.51% (305/316) orthology in GPCRs among the Fugu-human genomes and Fugu-chicken genomes, respectively.

Conclusions

This study reveals the position of pisces in vertebrate evolution from the GPCR perspective. Fugu can act as a reference model for the human genome for other protein families as well, going by the high orthology observed for GPCRs between Fugu and human. The evolutionary comparison of GPCR sequences between key vertebrate classes of mammals, birds and fish will help in identifying key functional residues and motifs so as to fill in the blanks in the evolution of GPCRs in vertebrates.

     

Comparative Sequence and Structural Analyses of G-Protein-Coupled Receptor Crystal Structures and Implications for Molecular Models

Background

Up until recently the only available experimental (high resolution) structure of a G-protein-coupled receptor (GPCR) was that of bovine rhodopsin. In the past few years the determination of GPCR structures has accelerated with three new receptors, as well as squid rhodopsin, being successfully crystallized. All share a common molecular architecture of seven transmembrane helices and can therefore serve as templates for building molecular models of homologous GPCRs. However, despite the common general architecture of these structures key differences do exist between them. The choice of which experimental GPCR structure(s) to use for building a comparative model of a particular GPCR is unclear and without detailed structural and sequence analyses, could be arbitrary. The aim of this study is therefore to perform a systematic and detailed analysis of sequence-structure relationships of known GPCR structures.

Methodology

We analyzed in detail conserved and unique sequence motifs and structural features in experimentally-determined GPCR structures. Deeper insight into specific and important structural features of GPCRs as well as valuable information for template selection has been gained. Using key features a workflow has been formulated for identifying the most appropriate template(s) for building homology models of GPCRs of unknown structure. This workflow was applied to a set of 14 human family A GPCRs suggesting for each the most appropriate template(s) for building a comparative molecular model.

Conclusions

The available crystal structures represent only a subset of all possible structural variation in family A GPCRs. Some GPCRs have structural features that are distributed over different crystal structures or which are not present in the templates suggesting that homology models should be built using multiple templates. This study provides a systematic analysis of GPCR crystal structures and a consistent method for identifying suitable templates for GPCR homology modelling that will help to produce more reliable three-dimensional models.     

Virtual screening of GPCRs: An <Emphasis Type="Italic">in silico</Emphasis> chemogenomics approach

Background  

The G-protein coupled receptor (GPCR) superfamily is currently the largest class of therapeutic targets. In silico prediction of interactions between GPCRs and small molecules in the transmembrane ligand-binding site is therefore a crucial step in the drug discovery process, which remains a daunting task due to the difficulty to characterize the 3D structure of most GPCRs, and to the limited amount of known ligands for some members of the superfamily. Chemogenomics, which attempts to characterize interactions between all members of a target class and all small molecules simultaneously, has recently been proposed as an interesting alternative to traditional docking or ligand-based virtual screening strategies.     

Towards Improved Quality of GPCR Models by Usage of Multiple Templates and Profile-Profile Comparison

G-protein coupled receptors (GPCRs) are targets of nearly one third of the drugs at the current pharmaceutical market. Despite their importance in many cellular processes the crystal structures are available for less than 20 unique GPCRs of the Rhodopsin-like class. Fortunately, even though involved in different signaling cascades, this large group of membrane proteins has preserved a uniform structure comprising seven transmembrane helices that allows quite reliable comparative modeling. Nevertheless, low sequence similarity between the GPCR family members is still a serious obstacle not only in template selection but also in providing theoretical models of acceptable quality. An additional level of difficulty is the prediction of kinks and bulges in transmembrane helices. Usage of multiple templates and generation of alignments based on sequence profiles may increase the rate of success in difficult cases of comparative modeling in which the sequence similarity between GPCRs is exceptionally low. Here, we present GPCRM, a novel method for fast and accurate generation of GPCR models using averaging of multiple template structures and profile-profile comparison. In particular, GPCRM is the first GPCR structure predictor incorporating two distinct loop modeling techniques: Modeller and Rosetta together with the filtering of models based on the Z-coordinate. We tested our approach on all unique GPCR structures determined to date and report its performance in comparison with other computational methods targeting the Rhodopsin-like class. We also provide a database of precomputed GPCR models of the human receptors from that class.

Availability

GPCRM server and database: http://gpcrm.biomodellab.eu     

An improved classification of G-protein-coupled receptors using sequence-derived features

Background  

G-protein-coupled receptors (GPCRs) play a key role in diverse physiological processes and are the targets of almost two-thirds of the marketed drugs. The 3 D structures of GPCRs are largely unavailable; however, a large number of GPCR primary sequences are known. To facilitate the identification and characterization of novel receptors, it is therefore very valuable to develop a computational method to accurately predict GPCRs from the protein primary sequences.     

Expression analysis of G Protein-Coupled Receptors in mouse macrophages

Background

Monocytes and macrophages express an extensive repertoire of G Protein-Coupled Receptors (GPCRs) that regulate inflammation and immunity. In this study we performed a systematic micro-array analysis of GPCR expression in primary mouse macrophages to identify family members that are either enriched in macrophages compared to a panel of other cell types, or are regulated by an inflammatory stimulus, the bacterial product lipopolysaccharide (LPS).

Results

Several members of the P2RY family had striking expression patterns in macrophages; P2ry6 mRNA was essentially expressed in a macrophage-specific fashion, whilst P2ry1 and P2ry5 mRNA levels were strongly down-regulated by LPS. Expression of several other GPCRs was either restricted to macrophages (e.g. Gpr84) or to both macrophages and neural tissues (e.g. P2ry12, Gpr85). The GPCR repertoire expressed by bone marrow-derived macrophages and thioglycollate-elicited peritoneal macrophages had some commonality, but there were also several GPCRs preferentially expressed by either cell population.

Conclusion

The constitutive or regulated expression in macrophages of several GPCRs identified in this study has not previously been described. Future studies on such GPCRs and their agonists are likely to provide important insights into macrophage biology, as well as novel inflammatory pathways that could be future targets for drug discovery.     

Structural features of the G-protein/GPCR interactions

Background

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

Scope of review

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

Major conclusions

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

General significance

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

On the origins of arrestin and rhodopsin

Background  

G protein coupled receptors (GPCRs) are the most numerous proteins in mammalian genomes, and the most common targets of clinical drugs. However, their evolution remains enigmatic. GPCRs are intimately associated with trimeric G proteins, G protein receptor kinases, and arrestins. We conducted phylogenetic studies to reconstruct the history of arrestins. Those findings, in turn, led us to investigate the origin of the photosensory GPCR rhodopsin.     

Identifying dramatic selection shifts in phylogenetic trees

Background

The rate of evolution varies spatially along genomes and temporally in time. The presence of evolutionary rate variation is an informative signal that often marks functional regions of genomes and historical selection events. There exist many tests for temporal rate variation, or heterotachy, that start by partitioning sampled sequences into two or more groups and testing rate homogeneity among the groups. I develop a Bayesian method to infer phylogenetic trees with a divergence point, or dramatic temporal shifts in selection pressure that affect many nucleotide sites simultaneously, located at an unknown position in the tree.

Results

Simulation demonstrates that the method is most able to detect divergence points when rate variation and the number of affected sites is high, but not beyond biologically relevant values. The method is applied to two viral data sets. A divergence point is identified separating the B and C subtypes, two genetically distinct variants of HIV that have spread into different human populations with the AIDS epidemic. In contrast, no strong signal of temporal rate variation is found in a sample of F and H genotypes, two genetic variants of HBV that have likely evolved with humans during their immigration and expansion into the Americas.

Conclusion

Temporal shifts in evolutionary rate of sufficient magnitude are detectable in the history of sampled sequences. The ability to detect such divergence points without the need to specify a prior hypothesis about the location or timing of the divergence point should help scientists identify historically important selection events and decipher mechanisms of evolution.

     

Allosteric activation of metabotropic glutamate receptor 5

Abstract

Metabotropic glutamate receptor 5 (mGluR5) is a class C G protein-coupled receptor (GPCR) with both an extracellular ligand binding site and an allosteric intrahelical chamber located similarly to the orthosteric ligand binding site of Class A GPCRs. Ligands binding to this ancestral site of mGluR5 can act as positive (PAM), negative (NAM) or silent (SAM) allosteric modulators, and their medicinal chemistry optimization is notoriously difficult, as subtle structural changes may cause significant variation in activity and switch in the functional response. Here we present all atom molecular dynamics simulations of NAM, SAM and PAM complexes formed by closely related ligands and analyse the structural differences of the complexes. Several residues involved in the activation are identified and the formation of a continuous water channel in the active complex but not in the inactive ones is recognized. Our results suggest that the mechanism of mGluR5 activation is similar to that of class A GPCRs.

Communicated by Ramaswamy H. Sarma     

Novel Structural and Functional Insights into M3 Muscarinic Receptor Dimer/Oligomer Formation

Class A G protein-coupled receptors (GPCRs) are able to form homodimers and/or oligomeric arrays. We recently proposed, based on bioluminescence resonance energy transfer studies with the M₃ muscarinic receptor (M3R), a prototypic class A GPCR, that the M3R is able to form multiple, structurally distinct dimers that are probably transient in nature (McMillin, S. M., Heusel, M., Liu, T., Costanzi, S., and Wess, J. (2011) J. Biol. Chem. 286, 28584–28598). To provide more direct experimental support for this concept, we employed a disulfide cross-linking strategy to trap various M3R dimeric species present in a native lipid environment (transfected COS-7 cells). Disulfide cross-linking studies were carried out with many mutant M3Rs containing single cysteine (Cys) substitutions within two distinct cytoplasmic M3R regions, the C-terminal portion of the second intracellular loop (i2) and helix H8 (H8). The pattern of cross-links that we obtained, in combination with molecular modeling studies, was consistent with the existence of two structurally distinct M3R dimer interfaces, one involving i2/i2 contacts (TM4-TM5-i2 interface) and the other one characterized by H8-H8 interactions (TM1-TM2-H8 interface). Specific H8-H8 disulfide cross-links led to significant impairments in M3R-mediated G protein activation, suggesting that changes in the structural orientation or mobility of H8 are critical for efficient receptor-G protein coupling. Our findings provide novel structural and functional insights into the mechanisms involved in M3R dimerization (oligomerization). Because the M3R shows a high degree of sequence similarity with many other class A GPCRs, our findings should be of considerable general interest.     

The repertoire of G protein-coupled receptors in the sea squirt <Emphasis Type="Italic">Ciona intestinalis</Emphasis>

Background  

G protein-coupled receptors (GPCRs) constitute a large family of integral transmembrane receptor proteins that play a central role in signal transduction in eukaryotes. The genome of the protochordate Ciona intestinalis has a compact size with an ancestral complement of many diversified gene families of vertebrates and is a good model system for studying protochordate to vertebrate diversification. An analysis of the Ciona repertoire of GPCRs from a comparative genomic perspective provides insight into the evolutionary origins of the GPCR signalling system in vertebrates.     

Molecular evolution and functional divergence of the cytochrome P450 3 (CYP3) Family in Actinopterygii (ray-finned fish)

Background

The cytochrome P450 (CYP) superfamily is a multifunctional hemethiolate enzyme that is widely distributed from Bacteria to Eukarya. The CYP3 family contains mainly the four subfamilies CYP3A, CYP3B, CYP3C and CYP3D in vertebrates; however, only the Actinopterygii (ray-finned fish) have all four subfamilies and detailed understanding of the evolutionary relationship of Actinopterygii CYP3 family members would be valuable.

Methods and Findings

Phylogenetic relationships were constructed to trace the evolutionary history of the Actinopterygii CYP3 family genes. Selection analysis, relative rate tests and functional divergence analysis were combined to interpret the relationship of the site-specific evolution and functional divergence in the Actinopterygii CYP3 family. The results showed that the four CYP3 subfamilies in Actinopterygii might be formed by gene duplication. The first gene duplication event was responsible for divergence of the CYP3B/C clusters from ancient CYP3 before the origin of the Actinopterygii, which corresponded to the fish-specific whole genome duplication (WGD). Tandem repeat duplication in each of the homologue clusters produced stable CYP3B, CYP3C, CYP3A and CYP3D subfamilies. Acceleration of asymmetric evolutionary rates and purifying selection together were the main force for the production of new subfamilies and functional divergence in the new subset after gene duplication, whereas positive selection was detected only in the retained CYP3A subfamily. Furthermore, nearly half of the functional divergence sites appear to be related to substrate recognition, which suggests that site-specific evolution is closely related with functional divergence in the Actinopterygii CYP3 family.

Conclusions

The split of fish-specific CYP3 subfamilies was related to the fish-specific WGD, and site-specific acceleration of asymmetric evolutionary rates and purifying selection was the main force for the origin of the new subfamilies and functional divergence in the new subset after gene duplication. Site-specific evolution in substrate recognition was related to functional divergence in the Actinopterygii CYP3 family.     

A database for G proteins and their interaction with GPCRs

Background  

G protein-coupled receptors (GPCRs) transduce signals from extracellular space into the cell, through their interaction with G proteins, which act as switches forming hetero-trimers composed of different subunits (α,β,γ). The α subunit of the G protein is responsible for the recognition of a given GPCR. Whereas specialised resources for GPCRs, and other groups of receptors, are already available, currently, there is no publicly available database focusing on G Proteins and containing information about their coupling specificity with their respective receptors.