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.     

Diverse Cell Type-Specific Mechanisms Localize G Protein-Coupled Receptors to Caenorhabditis elegans Sensory Cilia

The localization of signaling molecules such as G protein-coupled receptors (GPCRs) to primary cilia is essential for correct signal transduction. Detailed studies over the past decade have begun to elucidate the diverse sequences and trafficking mechanisms that sort and transport GPCRs to the ciliary compartment. However, a systematic analysis of the pathways required for ciliary targeting of multiple GPCRs in different cell types in vivo has not been reported. Here we describe the sequences and proteins required to localize GPCRs to the cilia of the AWB and ASK sensory neuron types in Caenorhabditis elegans. We find that GPCRs expressed in AWB or ASK utilize conserved and novel sequences for ciliary localization, and that the requirement for a ciliary targeting sequence in a given GPCR is different in different neuron types. Consistent with the presence of multiple ciliary targeting sequences, we identify diverse proteins required for ciliary localization of individual GPCRs in AWB and ASK. In particular, we show that the TUB-1 Tubby protein is required for ciliary localization of a subset of GPCRs, implying that defects in GPCR localization may be causal to the metabolic phenotypes of tub-1 mutants. Together, our results describe a remarkable complexity of mechanisms that act in a protein- and cell-specific manner to localize GPCRs to cilia, and suggest that this diversity allows for precise regulation of GPCR-mediated signaling as a function of external and internal context.     

3D Modeling for TM Receptors: Algorithms and Validations

This paper outlines the basic strategy to build 3D models of transmembrane G-protein coupled receptors (GPCRs) starting from their amino acid sequences in a block-by-block manner: (i) locate possible TM helical fragments in the GPCR sequence; (ii) build 3D structures for these helices; (iii) arrange isolated helices across the membrane; (iv) calculate all pairwise helix-helix interactions; (v) assemble helical bundle(s); (vi) restore interhelical loops and N- and C-termini; and (vii) refine the entire 3D structure(s). Computer algorithms and preliminary results for most of the steps are discussed.     

GCR2 is a new member of the eukaryotic lanthionine synthetase component C-like protein family

GCR2 was recently proposed to represent a G-protein-coupled receptor (GPCR) for the plant hormone, abscisic acid (ABA). We and others provided evidence that GCR2 is unlikely to be a bona fide GPCR because it is not clearly predicted to contain seven transmembrane domains, a structural hallmark for classical GPCRs. Instead, GCR2 shows significant sequence similarity to homologs of bacterial lanthionine synthetase component C (LanC). Here, we provide additional analysis of GCR2 and LanC-like (LANCL) proteins in plants, and propose that GCR2 is a new member of the eukaryotic LANCL protein family.Key words: GCR2, G-protein-coupled receptor, abscisic acid (ABA), lanthionine synthetaseSeven transmembrane (7TM) G-protein-coupled receptors (GPCRs) comprise the largest protein family in mammals, and are the most pharmacologically important receptor family, being the target of approximately half of all modern medicinal drugs. All canonical GPCRs are integral membrane proteins and are predicted to contain 7TM-spanning domains as their structural hallmark, a pattern confirmed by the high-resolution crystal structure of human β2-adrenergic GPCR.^1,2 GPCRs sense extracellular molecules and activate intracellular cell signaling via coupling with heterotrimeric G-proteins. Heterotrimeric G-protein subunits are conserved in plants, but the repertoire of heterotrimeric G-protein complexes to which they contribute in plants is much simpler than in mammals.^3,4 Liu et al. (2007) proposed that GCR2 is a GPCR for the plant hormone abscisic acid (ABA) in Arabidopsis.⁵ However, GCR2 was predicted not to be a 7TM protein when its amino acid sequence was analyzed in robust transmembrane prediction systems.^6,7 On the other hand, GCR2 has significant sequence similarity to homologs of bacterial lanthionine synthetase component C (LanC) that are found in diverse eukaryotes and which have predicted structural similarity to prokaryotic LanC.^6,7 These findings raise the possibility that GCR2 belongs to the LanC protein superfamily, rather than the GPCR superfamily.     

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

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

SFESA: a web server for pairwise alignment refinement by secondary structure shifts

Background

Protein sequence alignment is essential for a variety of tasks such as homology modeling and active site prediction. Alignment errors remain the main cause of low-quality structure models. A bioinformatics tool to refine alignments is needed to make protein alignments more accurate.

Results

We developed the SFESA web server to refine pairwise protein sequence alignments. Compared to the previous version of SFESA, which required a set of 3D coordinates for a protein, the new server will search a sequence database for the closest homolog with an available 3D structure to be used as a template. For each alignment block defined by secondary structure elements in the template, SFESA evaluates alignment variants generated by local shifts and selects the best-scoring alignment variant. A scoring function that combines the sequence score of profile-profile comparison and the structure score of template-derived contact energy is used for evaluation of alignments. PROMALS pairwise alignments refined by SFESA are more accurate than those produced by current advanced alignment methods such as HHpred and CNFpred. In addition, SFESA also improves alignments generated by other software.

Conclusions

SFESA is a web-based tool for alignment refinement, designed for researchers to compute, refine, and evaluate pairwise alignments with a combined sequence and structure scoring of alignment blocks. To our knowledge, the SFESA web server is the only tool that refines alignments by evaluating local shifts of secondary structure elements. The SFESA web server is available at http://prodata.swmed.edu/sfesa.     

UMARS: Un-MAppable Reads Solution

Background

Few high-resolution structures of integral membranes proteins are available, as crystallization of such proteins needs yet to overcome too many technical limitations. Nevertheless, prediction of their transmembrane (TM) structure by bioinformatics tools provides interesting insights on the topology of these proteins.

Methods

We describe here how to extract new information from the analysis of hydrophobicity variations or hydrophobic pulses (HPulses) in the sequence of integral membrane proteins using the Hydrophobic Pulse Predictor, a new tool we developed for this purpose. To analyze the primary sequence of 70 integral membrane proteins we defined two levels of analysis: G1-HPulses for sliding windows of n = 2 to 6 and G2-HPulses for sliding windows of n = 12 to 16.

Results

The G2-HPulse analysis of 541 transmembrane helices allowed the definition of the new concept of transmembrane unit (TMU) that groups together transmembrane helices and segments with potential adjacent structures. In addition, the G1-HPulse analysis identified helix irregularities that corresponded to kinks, partial helices or unannotated structural events. These irregularities could represent key dynamic elements that are alternatively activated depending on the channel status as illustrated by the crystal structures of the lactose permease in different conformations.

Conclusions

Our results open a new way in the understanding of transmembrane secondary structures: hydrophobicity through hydrophobic pulses strongly impacts on such embedded structures and is not confined to define the transmembrane status of amino acids.     

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.     

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.

     

Novel approach to computer modeling of seven-helical transmembrane proteins: current progress in the test case of bacteriorhodopsin.

G-protein coupled receptors (GPCRs) are thought to be proteins with 7-membered transmembrane helical bundles (7TM proteins). Recently, the X-ray structures have been solved for two such proteins, namely for bacteriorhodopsin (BR) and rhodopsin (Rh), the latter being a GPCR. Despite similarities, the structures are different enough to suggest that 3D models for different GPCRs cannot be obtained directly employing 3D structures of BR or Rh as a unique template. The approach to computer modeling of 7TM proteins developed in this work was capable of reproducing the experimental X-ray structure of BR with great accuracy. A combination of helical packing and low-energy conformers for loops most close to the X-ray structure possesses the r.m.s.d. value of 3.13 A. Such a level of accuracy for the 3D-structure prediction for a 216-residue protein has not been achieved, so far, by any available ab initio procedure of protein folding. The approach may produce also other energetically consistent combinations of helical bundles and loop conformers, creating a variety of possible templates for 3D structures of 7TM proteins, including GPCRs. These templates may provide experimentalists with various plausible options for 3D structure of a given GPCR; in our view, only experiments will determine the final choice of the most reasonable 3D template.     

PROCAIN: protein profile comparison with assisting information

Detection of remote sequence homology is essential for the accurate inference of protein structure, function and evolution. The most sensitive detection methods involve the comparison of evolutionary patterns reflected in multiple sequence alignments (MSAs) of protein families. We present PROCAIN, a new method for MSA comparison based on the combination of ‘vertical’ MSA context (substitution constraints at individual sequence positions) and ‘horizontal’ context (patterns of residue content at multiple positions). Based on a simple and tractable profile methodology and primitive measures for the similarity of horizontal MSA patterns, the method achieves the quality of homology detection comparable to a more complex advanced method employing hidden Markov models (HMMs) and secondary structure (SS) prediction. Adding SS information further improves PROCAIN performance beyond the capabilities of current state-of-the-art tools. The potential value of the method for structure/function predictions is illustrated by the detection of subtle homology between evolutionary distant yet structurally similar protein domains. ProCAIn, relevant databases and tools can be downloaded from: http://prodata.swmed.edu/procain/download. The web server can be accessed at http://prodata.swmed.edu/procain/procain.php.     

3D Modeling for TM Receptors: Algorithms and Validations

Summary This paper outlines the basic strategy to build 3D models of transmembrane G-protein coupled receptors (GPCRs) starting from their amino acid sequences in a ‘block-by-block’ manner: (i) locate possible TM helical fragments in the GPCR sequence; (ii) build 3D structures for these helices; (iii) arrange isolated helices across the membrane; (iv) calculate all pairwise helix-helix interactions; (v) assemble helical bundle(s); (vi) restore interhelical loops and N- and C-termini; and (vii) refine the entire 3D structure(s). Computer algorithms and preliminary results for most of the steps are discussed.     

Translational Control of the Oogenic Program by Components of OMA Ribonucleoprotein Particles in Caenorhabditis elegans

The oocytes of most sexually reproducing animals arrest in meiotic prophase I. Oocyte growth, which occurs during this period of arrest, enables oocytes to acquire the cytoplasmic components needed to produce healthy progeny and to gain competence to complete meiosis. In the nematode Caenorhabditis elegans, the major sperm protein hormone promotes meiotic resumption (also called meiotic maturation) and the cytoplasmic flows that drive oocyte growth. Prior work established that two related TIS11 zinc-finger RNA-binding proteins, OMA-1 and OMA-2, are redundantly required for normal oocyte growth and meiotic maturation. We affinity purified OMA-1 and identified associated mRNAs and proteins using genome-wide expression data and mass spectrometry, respectively. As a class, mRNAs enriched in OMA-1 ribonucleoprotein particles (OMA RNPs) have reproductive functions. Several of these mRNAs were tested and found to be targets of OMA-1/2-mediated translational repression, dependent on sequences in their 3′-untranslated regions (3′-UTRs). Consistent with a major role for OMA-1 and OMA-2 in regulating translation, OMA-1-associated proteins include translational repressors and activators, and some of these proteins bind directly to OMA-1 in yeast two-hybrid assays, including OMA-2. We show that the highly conserved TRIM-NHL protein LIN-41 is an OMA-1-associated protein, which also represses the translation of several OMA-1/2 target mRNAs. In the accompanying article in this issue, we show that LIN-41 prevents meiotic maturation and promotes oocyte growth in opposition to OMA-1/2. Taken together, these data support a model in which the conserved regulators of mRNA translation LIN-41 and OMA-1/2 coordinately control oocyte growth and the proper spatial and temporal execution of the meiotic maturation decision.     

Comprehensive FISH Probe Design Tool Applied to Imaging Human Immunoglobulin Class Switch Recombination

We present a web engine boosted fluorescence in-situ hybridization (webFISH) algorithm using a genome-wide sequence similarity search to design target-specific single-copy and repetitive DNA FISH probes. The webFISH algorithm featuring a user-friendly interface (http://www.webfish2.org/) maximizes the coverage of the examined sequences with FISH probes by considering locally repetitive sequences absent from the remainder of the genome. The highly repetitive human immunoglobulin heavy chain sequence was analyzed using webFISH to design three sets of FISH probes. These allowed direct simultaneous detection of class switch recombination in both immunoglobulin-heavy chain alleles in single cells from a population of cultured primary B cells. It directly demonstrated asynchrony of the class switch recombination in the two alleles in structurally preserved nuclei while permitting parallel readout of protein expression by immunofluorescence staining. This novel technique offers the possibility of gaining unprecedented insight into the molecular mechanisms involved in class switch recombination.     

Functional Genomic Identification of Genes Required for Male Gonadal Differentiation in Caenorhabditis elegans

The Caenorhabditis elegans somatic gonad develops from a four-cell primordium into a mature organ that differs dramatically between the sexes in overall morphology (two arms in hermaphrodites and one in males) and in the cell types comprising it. Gonadal development in C. elegans is well studied, but regulation of sexual differentiation, especially later in gonadal development, remains poorly elucidated. To identify genes involved in this process, we performed a genome-wide RNAi screen using sex-specifically expressed gonadal GFP reporters. This screen identified several phenotypic classes, including ∼70 genes whose depletion feminized male gonadal cells. Among the genes required for male cell fate specification are Wnt/β-catenin pathway members, cell cycle regulators, and genes required for mitotic spindle function and cytokinesis. We find that a Wnt/β-catenin pathway independent of extracellular Wnt ligand is essential for asymmetric cell divisions and male differentiation during gonadal development in larvae. We also find that the cell cycle regulators cdk-1 and cyb-3 and the spindle/cytokinesis regulator zen-4 are required for Wnt/β-catenin pathway activity in the developing gonad. After sex is determined in the gonadal primordium the global sex determination pathway is dispensable for gonadal sexual fate, suggesting that male cell fates are promoted and maintained independently of the global pathway during this period.THE Caenorhabditis elegans gonad derives from a simple primordium of four cells that coalesces during embryogenesis and contains two somatic gonad precursors (SGPs), Z1 and Z4, flanking two germline precursors, Z2 and Z3 (Kimble and Hirsh 1979). The SGPs undergo very different developmental programs in each sex, involving sexually dimorphic cell lineages and migrations and sex-specific cellular differentiation. The result is a two-armed bilaterally symmetrical gonad in the adult hermaphrodite or a single-armed asymmetric gonad in the adult male. The high degree of sexual dimorphism of the mature organ and variety of cellular events that occur sex specifically during its development make the C. elegans gonad an outstanding model for sex-specific organogenesis.Development of the somatic gonad occurs in two phases. The early phase defines the gonadal axes and establishes the precursors of the major gonadal cell types. This takes place during the first larval stage (L1), beginning shortly after hatching with the first division of the SGPs. In both sexes SGP division is asymmetric in terms of both the sizes and the fates of the daughter cells, and establishes the proximal/distal axis of the gonad (Hirsh et al. 1976; Kimble and Hirsh 1979). The global sex determination pathway establishes the future sex of the gonad around the time of hatching (Klass et al. 1976; Nelson et al. 1978), and sexual dimorphism is already apparent when the SGPs divide: the size asymmetry of the SGP daughters is much more pronounced in males than hermaphrodites. In both sexes the asymmetry of the first SGP division requires a Wnt/β-catenin pathway. Mutations compromising this pathway cause a “symmetrical sisters” phenotype in which both daughters adopt the same fate (Miskowski et al. 2001; Siegfried and Kimble 2002; Phillips and Kimble 2009). Sex specificity is imposed on the SGPs by the global sex determining gene tra-1 (Hodgkin 1987) and the gonad-specific sex determining gene fkh-6 (Chang et al. 2004). These genes play opposing roles in SGP sex determination, with tra-1 feminizing and fkh-6 masculinizing the somatic gonad, and they also act redundantly to promote mitotic proliferation of the SGP lineage (Chang et al. 2004). SGP sex determination is linked to cell cycle progression by cyclin D, which is required to overcome repression of fkh-6 expression in the SGPs by E2F (Tilmann and Kimble 2005).The later phase of gonadal development involves the elongation of the gonad, together with cellular proliferation and differentiation, and lasts from L2 to adulthood. During L2 the somatic cells enlarge and leader cells (distal tip cells in the hermaphrodite, linker cell in the male) begin long-range migrations that extend the gonad. During L3, somatic gonad cell division resumes in both sexes, leading to the formation of differentiated somatic cell types by the end of L3 or beginning of L4. Gonadal morphogenesis is completed and gametogenesis begins during L4 (Kimble and Hirsh 1979).Although SGP division and much of hermaphrodite gonadal development have been well studied (Hubbard and Greenstein 2000), sexual cell fate specification in the somatic gonad is more poorly understood, particularly after the L1 stage. Despite the importance of fkh-6 in promoting male differentiation, it is expressed in males only during early L1 and null mutants have incomplete gonadal sex reversal. We have therefore performed a genome-wide RNAi screen to identify additional genes required after hatching for gonadal development in each sex. Among the advantages of this approach is the ability to identify gonadal regulators that also are essential for embryonic development. To our knowledge this is the first functional genomic study of gonadal sex differentiation.The screen identified many genes whose depletion disrupts gonadogenesis in each sex and nearly 70 genes whose depletion causes gonadal feminization in males. Prominent among this latter class were components of a Wnt/β-catenin pathway, cell cycle regulators, and genes involved in mitotic spindle function and cytokinesis. We find that Wnt/β-catenin activity continues in both sexes after SGP division and is required for male cell fate commitment in the gonad. We also find that the cyclin-dependent kinase cdk-1 and its cognate cyclin cyb-3 as well as the mitotic spindle regulator zen-4 are required for gonadal Wnt/β-catenin pathway activity, providing a potential new link between the cell cycle, asymmetric division, and sexual differentiation. The feminization caused by depletion of Wnt/β-catenin pathway components or cdk-1 is independent of the global sex determination pathway, suggesting that sexual fates in the male gonad remain plastic after the primary sex determination decision.     

High quality draft genome sequence of Staphylococcus cohnii subsp. cohnii strain hu-01

Staphylococcus cohnii subsp. cohnii belongs to the family Staphylococcaceae in the order Bacillales, class Bacilli and phylum Firmicutes. The increasing relevance of S. cohnii to human health prompted us to determine the genomic sequence of Staphylococcus cohnii subsp. cohnii strain hu-01, a multidrug-resistant isolate from a hospital in China. Here we describe the features of S. cohnii subsp. cohnii strain hu-01, together with the genome sequence and its annotation. This is the first genome sequence of the species Staphylococcus cohnii.