PipeAlign: A new toolkit for protein family analysis

ESG: extended similarity group method for automated protein function prediction

MOTIVATION

Importance of accurate automatic protein function prediction is ever increasing in the face of a large number of newly sequenced genomes and proteomics data that are awaiting biological interpretation. Conventional methods have focused on high sequence similarity-based annotation transfer which relies on the concept of homology. However, many cases have been reported that simple transfer of function from top hits of a homology search causes erroneous annotation. New methods are required to handle the sequence similarity in a more robust way to combine together signals from strongly and weakly similar proteins for effectively predicting function for unknown proteins with high reliability.

RESULTS

We present the extended similarity group (ESG) method, which performs iterative sequence database searches and annotates a query sequence with Gene Ontology terms. Each annotation is assigned with probability based on its relative similarity score with the multiple-level neighbors in the protein similarity graph. We will depict how the statistical framework of ESG improves the prediction accuracy by iteratively taking into account the neighborhood of query protein in the sequence similarity space. ESG outperforms conventional PSI-BLAST and the protein function prediction (PFP) algorithm. It is found that the iterative search is effective in capturing multiple-domains in a query protein, enabling accurately predicting several functions which originate from different domains.

AVAILABILITY

ESG web server is available for automated protein function prediction at http://dragon.bio.purdue.edu/ESG/.

Structural bases for 16 S rRNA methylation catalyzed by ArmA and RmtB methyltransferases

Aminoglycosides are used extensively for the treatment of severe infections due to Gram-negative bacteria. However, certain species have become highly resistant after acquisition of genes for methyltransferases which catalyze post-transcriptional methylation of N7-G1405 in 16 S rRNA of 30 S ribosomal subunits. Inactivation of this enzymatic activity is therefore an important challenge for development of an effective therapy. The present work describes the crystallographic structures of methyltransferases RmtB and ArmA from clinical isolates. Together with biochemical experiments, the 3D structures indicate that the N-terminal domain specific for this family of methyltransferases is required for enzymatic activity. Site-directed mutagenesis has enabled important residues for catalysis and RNA binding to be identified. These high-resolution structures should underpin the design of potential inhibitors of these enzymes, which could be used to restore the activity of aminoglycosides against resistant pathogens.

Structure of the archaeal pab87 peptidase reveals a novel self-compartmentalizing protease family

Self-compartmentalizing proteases orchestrate protein turnover through an original architecture characterized by a central catalytic chamber. Here we report the first structure of an archaeal member of a new self-compartmentalizing protease family forming a cubic-shaped octamer with D₄ symmetry and referred to as CubicO. We solved the structure of the Pyrococcus abyssi Pab87 protein at 2.2 Å resolution using the anomalous signal of the high-phasing-power lanthanide derivative Lu-HPDO3A. A 20 Å wide channel runs through this supramolecular assembly of 0.4 MDa, giving access to a 60 Å wide central chamber holding the eight active sites. Surprisingly, activity assays revealed that Pab87 degrades specifically d-amino acid containing peptides, which have never been observed in archaea. Genomic context of the Pab87 gene showed that it is surrounded by genes involved in the amino acid/peptide transport or metabolism. We propose that CubicO proteases are involved in the processing of d-peptides from environmental origins.

Knowledge-based expert systems and a proof-of-concept case study for multiple sequence alignment construction and analysis

The traditional approach to bioinformatics analyses relies on independent task-specific services and applications, using different input and output formats, often idiosyncratic, and frequently not designed to inter-operate. In general, such analyses were performed by experts who manually verified the results obtained at each step in the process. Today, the amount of bioinformatics information continuously being produced means that handling the various applications used to study this information presents a major data management and analysis challenge to researchers. It is now impossible to manually analyse all this information and new approaches are needed that are capable of processing the large-scale heterogeneous data in order to extract the pertinent information. We review the recent use of integrated expert systems aimed at providing more efficient knowledge extraction for bioinformatics research. A general methodology for building knowledge-based expert systems is described, focusing on the unstructured information management architecture, UIMA, which provides facilities for both data and process management. A case study involving a multiple alignment expert system prototype called AlexSys is also presented.

Function and ribosomal localization of aIF6, a translational regulator shared by archaea and eukarya

The translation factor IF6 is shared by the Archaea and the Eukarya, but is not found in Bacteria. The properties of eukaryal IF6 (eIF6) have been extensively studied, but remain somewhat elusive. eIF6 behaves as a ribosome-anti-association factor and is involved in miRNA-mediated gene silencing; however, it also seems to participate in ribosome synthesis and export. Here we have determined the function and ribosomal localization of the archaeal (Sulfolobus solfataricus) IF6 homologue (aIF6). We find that aIF6 binds specifically to the 50S ribosomal subunits, hindering the formation of 70S ribosomes and strongly inhibiting translation. aIF6 is uniformly expressed along the cell cycle, but it is upregulated following both cold- and heat shock. The aIF6 ribosomal binding site lies in the middle of the 30-S interacting surface of the 50S subunit, including a number of critical RNA and protein determinants involved in subunit association. The data suggest that the IF6 protein evolved in the archaeal–eukaryal lineage to modulate translational efficiency under unfavourable environmental conditions, perhaps acquiring additional functions during eukaryotic evolution.

Ortho-proteogenomics: multiple proteomes investigation through orthology and a new MS-based protocol

The progress in sequencing technologies irrigates biology with an ever-increasing number of genome sequences. In most cases, the gene repertoire is predicted in silico and conceptually translated into proteins. As recently highlighted, the predicted genes exhibit frequent errors, particularly in start codons, with a serious impact on subsequent biological studies. A new "ortho-proteogenomic" approach is presented here for the annotation refinement of multiple genomes at once. It combines comparative genomics with an original proteomic protocol that allows the characterization of both N-terminal and internal peptides in a single experiment. This strategy was applied to the Mycobacterium genus with Mycobacterium smegmatis as the reference, and identified 946 distinct proteins, including 443 characterized N termini. These experimental data allowed the correction of 19% of the characterized start codons, the identification of 29 proteins missed during the annotation process, and the curation, thanks to comparative genomics, of 4328 sequences of 16 other Mycobacterium proteomes.

PtdIns5P regulation through evolution: roles in membrane trafficking?

Phosphoinositides are lipid second messengers that are essential for many cellular processes, including signal transduction and cell compartmentalization. Among them, phosphatidylinositol 5-phosphate (PtdIns5P) is the least characterized, although several proteins involved in its regulation are implicated in human diseases. We studied the distribution of 32 PtdIns5P-metabolizing proteins in 39 eukaryotic genomes. Phylogenetic profiles identify four groups of co-evolving proteins, confirming known protein complexes and revealing new ones. The complexes comprise a phosphatase, a kinase and a regulator; this indicates that physical interactions between the three partners are necessary for the acute spatial regulation of PtdIns5P turnover. By examining PtdIns5P metabolism in this new perspective, we propose a role for PtdIns5P in membrane trafficking from late endosomal compartments to the plasma membrane.

RPS19 mutations in patients with Diamond-Blackfan anemia

Diamond-Blackfan anemia (DBA) is an inherited disease characterized by pure erythroid aplasia. Thirty percent (30%) of patients display malformations, especially of the hands, face, heart, and urogenital tract. DBA has an autosomal dominant pattern of inheritance. De novo mutations are common and familial cases display wide clinical heterogeneity. Twenty-five percent (25%) of patients carry a mutation in the ribosomal protein (RP) S19 gene, whereas mutations in RPS24, RPS17, RPL35A, RPL11, and RPL5 are rare. These genes encode for structural proteins of the ribosome. A link between ribosomal functions and erythroid aplasia is apparent in DBA, but its etiology is not clear. Most authors agree that a defect in protein synthesis in a rapidly proliferating tissue, such as the erythroid bone marrow, may explain the defective erythropoiesis. A total of 77 RPS19 mutations have been described. Most are whole gene deletions, translocations, or truncating mutations (nonsense or frameshift), suggesting that haploinsufficiency is the basis of DBA pathology. A total of 22 missense mutations have also been described and several works have provided in vitro functional data for the mutant proteins. This review looks at the data on all these mutations, proposes a functional classification, and describes six new mutations. It is shown that patients with RPS19 mutations display a poorer response to steroids and a worse long-term prognosis compared to other DBA patients.

The repertoire of G protein-coupled receptors in the sea squirt Ciona intestinalis

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.

RESULTS

We have identified 169 gene products in the Ciona genome that code for putative GPCRs. Phylogenetic analyses reveal that Ciona GPCRs have homologous representatives from the five major GRAFS (Glutamate, Rhodopsin, Adhesion, Frizzled and Secretin) families concomitant with other vertebrate GPCR repertoires. Nearly 39% of Ciona GPCRs have unambiguous orthologs of vertebrate GPCR families, as defined for the human, mouse, puffer fish and chicken genomes. The Rhodopsin family accounts for ~68% of the Ciona GPCR repertoire wherein the LGR-like subfamily exhibits a lineage specific gene expansion of a group of receptors that possess a novel domain organisation hitherto unobserved in metazoan genomes.

CONCLUSION

Comparison of GPCRs in Ciona to that in human reveals a high level of orthology of a protochordate repertoire with that of vertebrate GPCRs. Our studies suggest that the ascidians contain the basic ancestral complement of vertebrate GPCR genes. This is evident at the subfamily level comparisons since Ciona GPCR sequences are significantly analogous to vertebrate GPCR subfamilies even while exhibiting Ciona specific genes. Our analysis provides a framework to perform future experimental and comparative studies to understand the roles of the ancestral chordate versions of GPCRs that predated the divergence of the urochordates and the vertebrates.

A new protein linear motif benchmark for multiple sequence alignment software

BACKGROUND

Linear motifs (LMs) are abundant short regulatory sites used for modulating the functions of many eukaryotic proteins. They play important roles in post-translational modification, cell compartment targeting, docking sites for regulatory complex assembly and protein processing and cleavage. Methods for LM detection are now being developed that are strongly dependent on scores for motif conservation in homologous proteins. However, most LMs are found in natively disordered polypeptide segments that evolve rapidly, unhindered by structural constraints on the sequence. These regions of modular proteins are difficult to align using classical multiple sequence alignment programs that are specifically optimised to align the globular domains. As a consequence, poor motif alignment quality is hindering efforts to detect new LMs.

RESULTS

We have developed a new benchmark, as part of the BAliBASE suite, designed to assess the ability of standard multiple alignment methods to detect and align LMs. The reference alignments are organised into different test sets representing real alignment problems and contain examples of experimentally verified functional motifs, extracted from the Eukaryotic Linear Motif (ELM) database. The benchmark has been used to evaluate and compare a number of multiple alignment programs. With distantly related proteins, the worst alignment program correctly aligns 48% of LMs compared to 73% for the best program. However, the performance of all the programs is adversely affected by the introduction of other sequences containing false positive motifs. The ranking of the alignment programs based on LM alignment quality is similar to that observed when considering full-length protein alignments, however little correlation was observed between LM and overall alignment quality for individual alignment test cases.

CONCLUSION

We have shown that none of the programs currently available is capable of reliably aligning LMs in distantly related sequences and we have highlighted a number of specific problems. The results of the tests suggest possible ways to improve program accuracy for difficult, divergent sequences.

Mutagenesis in the alpha3alpha4 GyrA helix and in the Toprim domain of GyrB refines the contribution of Mycobacterium tuberculosis DNA gyrase to intrinsic resistance to quinolones

The replacement of M74 in GyrA, A83 in GyrA, and R447 in GyrB of Mycobacterium tuberculosis gyrase by their Escherichia coli homologs resulted in active enzymes as quinolone susceptible as the E. coli gyrase. This demonstrates that the primary structure of gyrase determines intrinsic quinolone resistance and was supported by a three-dimensional model of N-terminal GyrA.

ADCK3, an ancestral kinase, is mutated in a form of recessive ataxia associated with coenzyme Q10 deficiency

Muscle coenzyme Q(10) (CoQ(10) or ubiquinone) deficiency has been identified in more than 20 patients with presumed autosomal-recessive ataxia. However, mutations in genes required for CoQ(10) biosynthetic pathway have been identified only in patients with infantile-onset multisystemic diseases or isolated nephropathy. Our SNP-based genome-wide scan in a large consanguineous family revealed a locus for autosomal-recessive ataxia at chromosome 1q41. The causative mutation is a homozygous splice-site mutation in the aarF-domain-containing kinase 3 gene (ADCK3). Five additional mutations in ADCK3 were found in three patients with sporadic ataxia, including one known to have CoQ(10) deficiency in muscle. All of the patients have childhood-onset cerebellar ataxia with slow progression, and three of six have mildly elevated lactate levels. ADCK3 is a mitochondrial protein homologous to the yeast COQ8 and the bacterial UbiB proteins, which are required for CoQ biosynthesis. Three out of four patients tested showed a low endogenous pool of CoQ(10) in their fibroblasts or lymphoblasts, and two out of three patients showed impaired ubiquinone synthesis, strongly suggesting that ADCK3 is also involved in CoQ(10) biosynthesis. The deleterious nature of the three identified missense changes was confirmed by the introduction of them at the corresponding positions of the yeast COQ8 gene. Finally, a phylogenetic analysis shows that ADCK3 belongs to the family of atypical kinases, which includes phosphoinositide and choline kinases, suggesting that ADCK3 plays an indirect regulatory role in ubiquinone biosynthesis possibly as part of a feedback loop that regulates ATP production.

Homology and phylogeny and their automated inference

The analysis of the ever-increasing amount of biological and biomedical data can be pushed forward by comparing the data within and among species. For example, an integrative analysis of data from the genome sequencing projects for various species traces the evolution of the genomes and identifies conserved and innovative parts. Here, I review the foundations and advantages of this "historical" approach and evaluate recent attempts at automating such analyses. Biological data is comparable if a common origin exists (homology), as is the case for members of a gene family originating via duplication of an ancestral gene. If the family has relatives in other species, we can assume that the ancestral gene was present in the ancestral species from which all the other species evolved. In particular, describing the relationships among the duplicated biological sequences found in the various species is often possible by a phylogeny, which is more informative than homology statements. Detecting and elaborating on common origins may answer how certain biological sequences developed, and predict what sequences are in a particular species and what their function is. Such knowledge transfer from sequences in one species to the homologous sequences of the other is based on the principle of 'my closest relative looks and behaves like I do', often referred to as 'guilt by association'. To enable knowledge transfer on a large scale, several automated 'phylogenomics pipelines' have been developed in recent years, and seven of these will be described and compared. Overall, the examples in this review demonstrate that homology and phylogeny analyses, done on a large (and automated) scale, can give insights into function in biology and biomedicine.

Tex19, a mammalian-specific protein with a restricted expression in pluripotent stem cells and germ line

Although the properties of embryonic stem (ES) cells make these cells very attractive in the field of replacement therapy, the molecular mechanisms involved in the maintenance of their pluripotency are not fully characterized. Starting from the observation that most pluripotent markers are also expressed by spermatogonia stem cells, we identified Tex19 as a new potential pluripotency marker. We show that Tex19 is a mammalian-specific protein duplicated in mouse and rat, renamed Tex19.1 and Tex19.2, whereas only one form is found in human. In mouse, both forms are localized on chromosome 11 and transcribed in opposite directions. Tex19 proteins are well conserved, showing two highly conserved domains that do not present any similarity with any other known domains. We show that Tex19.2 is specifically detected in the male somatic gonad lineage, whereas Tex19.1 expression is very similar to that of Oct4. Transcripts are maternally inherited, and expression starts as soon as the early embryo and later is limited to the germ line. Tex19.1 transcripts were also detected in mouse pluripotent stem cells, and expression of Tex19.1, like that of Oct4, decreases after murine embryonic stem and germ cell differentiation. Human TEX19 was more closely related to murine Tex19.1 and was also detected in adult testis and in undifferentiated ES cells. By immunofluorescence, we found that Tex19.1 protein localizes to the nucleus of mouse ES and inner cell mass cells. All these results suggest that Tex19.1, as well as human TEX19, could be a new factor involved in the maintenance of self-renewal or pluripotency of stem cells.

Structured mRNAs regulate translation initiation by binding to the platform of the ribosome

Gene expression can be regulated at the level of initiation of protein biosynthesis via structural elements present at the 5' untranslated region of mRNAs. These folded mRNA segments may bind to the ribosome, thus blocking translation until the mRNA unfolds. Here, we report a series of cryo-electron microscopy snapshots of ribosomal complexes directly visualizing either the mRNA structure blocked by repressor protein S15 or the unfolded, active mRNA. In the stalled state, the folded mRNA prevents the start codon from reaching the peptidyl-tRNA (P) site inside the ribosome. Upon repressor release, the mRNA unfolds and moves into the mRNA channel allowing translation initiation. A comparative structure and sequence analysis suggests the existence of a universal stand-by site on the ribosome (the 30S platform) dedicated for binding regulatory 5' mRNA elements. Different types of mRNA structures may be accommodated during translation preinitiation and regulate gene expression by transiently stalling the ribosome.

Rod-derived Cone Viability Factor-2 is a novel bifunctional-thioredoxin-like protein with therapeutic potential

Background

Cone degeneration is the hallmark of the inherited retinal disease retinitis pigmentosa. We have previously identified a trophic factor "Rod-derived Cone Viability Factor (RdCVF) that is secreted by rods and promote cone viability in a mouse model of the disease.

Results

Here we report the bioinformatic identification and the experimental analysis of RdCVF2, a second trophic factor belonging to the Rod-derived Cone Viability Factor family. The mouse RdCVF gene is known to be bifunctional, encoding both a long thioredoxin-like isoform (RdCVF-L) and a short isoform with trophic cone photoreceptor viability activity (RdCVF-S). RdCVF2 shares many similarities with RdCVF in terms of gene structure, expression in a rod-dependent manner and protein 3D structure. Furthermore, like RdCVF, the RdCVF2 short isoform exhibits cone rescue activity that is independent of its putative thiol-oxydoreductase activity.

Conclusion

Taken together, these findings define a new family of bifunctional genes which are: expressed in vertebrate retina, encode trophic cone viability factors, and have major therapeutic potential for human retinal neurodegenerative diseases such as retinitis pigmentosa.

Molecular basis of Diamond-Blackfan anemia: structure and function analysis of RPS19

Diamond–Blackfan anemia (DBA) is a rare congenital disease linked to mutations in the ribosomal protein genes rps19, rps24 and rps17. It belongs to the emerging class of ribosomal disorders. To understand the impact of DBA mutations on RPS19 function, we have solved the crystal structure of RPS19 from Pyrococcus abyssi. The protein forms a five α-helix bundle organized around a central amphipathic α-helix, which corresponds to the DBA mutation hot spot. From the structure, we classify DBA mutations relative to their respective impact on protein folding (class I) or on surface properties (class II). Class II mutations cluster into two conserved basic patches. In vivo analysis in yeast demonstrates an essential role for class II residues in the incorporation into pre-40S ribosomal particles. This data indicate that missense mutations in DBA primarily affect the capacity of the protein to be incorporated into pre-ribosomes, thus blocking maturation of the pre-40S particles.

Automated protein subfamily identification and classification

Function prediction by homology is widely used to provide preliminary functional annotations for genes for which experimental evidence of function is unavailable or limited. This approach has been shown to be prone to systematic error, including percolation of annotation errors through sequence databases. Phylogenomic analysis avoids these errors in function prediction but has been difficult to automate for high-throughput application. To address this limitation, we present a computationally efficient pipeline for phylogenomic classification of proteins. This pipeline uses the SCI-PHY (Subfamily Classification in Phylogenomics) algorithm for automatic subfamily identification, followed by subfamily hidden Markov model (HMM) construction. A simple and computationally efficient scoring scheme using family and subfamily HMMs enables classification of novel sequences to protein families and subfamilies. Sequences representing entirely novel subfamilies are differentiated from those that can be classified to subfamilies in the input training set using logistic regression. Subfamily HMM parameters are estimated using an information-sharing protocol, enabling subfamilies containing even a single sequence to benefit from conservation patterns defining the family as a whole or in related subfamilies. SCI-PHY subfamilies correspond closely to functional subtypes defined by experts and to conserved clades found by phylogenetic analysis. Extensive comparisons of subfamily and family HMM performances show that subfamily HMMs dramatically improve the separation between homologous and non-homologous proteins in sequence database searches. Subfamily HMMs also provide extremely high specificity of classification and can be used to predict entirely novel subtypes. The SCI-PHY Web server at http://phylogenomics.berkeley.edu/SCI-PHY/ allows users to upload a multiple sequence alignment for subfamily identification and subfamily HMM construction. Biologists wishing to provide their own subfamily definitions can do so. Source code is available on the Web page. The Berkeley Phylogenomics Group PhyloFacts resource contains pre-calculated subfamily predictions and subfamily HMMs for more than 40,000 protein families and domains at http://phylogenomics.berkeley.edu/phylofacts/.

The structure of the CstF-77 homodimer provides insights into CstF assembly

The cleavage stimulation factor (CstF) is essential for the first step of poly(A) tail formation at the 3' ends of mRNAs. This heterotrimeric complex is built around the 77-kDa protein bridging both CstF-64 and CstF-50 subunits. We have solved the crystal structure of the 77-kDa protein from Encephalitozoon cuniculi at a resolution of 2 Å. The structure folds around 11 Half-a-TPR repeats defining two domains. The crystal structure reveals a tight homodimer exposing phylogenetically conserved areas for interaction with protein partners. Mapping experiments identify the C-terminal region of Rna14p, the yeast counterpart of CstF-77, as the docking domain for Rna15p, the yeast CstF-64 homologue.

Cloning and early expression pattern of two melatonin biosynthesis enzymes in the turbot (Scophthalmus maximus)

Melatonin biosynthesis from serotonin involves the sequential activation of the arylalkylamine N-acetyltransferase (AANAT) and hydroxyindole-O-methyltransferase (HIOMT). Photoperiod synchronizes a daily rhythm in pineal and retinal melatonin secretion through controlling AANAT activity. Teleost fish possess two Aanat, one expressed in the retina (AANAT1) and the other expressed in the pineal gland (AANAT2). We report here the full-length cloning of Aanat1, Aanat2, SmHiomt and Otx5 (orthodenticle homeobox homolog 5) in the turbot (Scophthalmus maximus, Sm), a flatfish belonging to an evolutionary recent group of Teleost. The temporal expression pattern of the genes investigated is consistent with the idea that OTX5 is needed for photoreceptor specification, and that the pineal gland differentiates before the retina. SmAanat2 expression remained pineal specific during the period of time investigated, whereas SmOtx5 and SmHiomt expressions were seen in both the retina and pineal gland. Our results do not support the existence of a second SmHiomt, as is the case for SmAanat. Neither SmAanat2 nor SmHiomt mRNAs displayed cyclic accumulation in the pineal organ of embryos and larvae maintained under a light-dark cycle from fertilization onward. This is in marked contrast with the situation observed with zebrafish Aanat2, indicating that the molecular mechanisms controlling the development of the pineal melatonin system have been modified during the evolution of Teleost.

Effects of single nucleotide polymorphisms on Toll-like receptor 3 activity and expression in cultured cells

Recognition of double-stranded RNA by Toll-like receptor 3 (TLR3) will increase the production of cytokines and chemokines through transcriptional activation by the NF-kappaB protein. Over 136 single-nucleotide polymorphisms (SNPs) in TLR3 have been identified in the human population. Of these, four alter the sequence of the TLR3 protein. Molecular modeling suggests that two of the SNPs, N284I and L412F, could affect the packing of the leucine-rich repeating units in TLR3. Notably, L412F is reported to be present in 20% of the population and is higher in the asthmatic population. To examine whether the four SNPs affect TLR3 function, each were cloned and tested for their ability to activate the expression of TLR3-dependent reporter constructs. SNP N284I was nearly completely defective for activating reporter activity, and L412F was reduced in activity. These two SNPs did not obviously affect the level of TLR3 expression or their intracellular location in vesicles. However, N284I and L412F were underrepresented on the cell surface, as determined by flow cytometry analysis, and were not efficiently secreted into the culture medium when expressed as the soluble ectodomain. They were also reduced in their ability to act in a dominant negative fashion on the wild type TLR3 allele. These observations suggest that N284I and L412F affect the activities of TLR3 needed for proper signaling.

Blast sampling for structural and functional analyses

BACKGROUND

The post-genomic era is characterised by a torrent of biological information flooding the public databases. As a direct consequence, similarity searches starting with a single query sequence frequently lead to the identification of hundreds, or even thousands of potential homologues. The huge volume of data renders the subsequent structural, functional and evolutionary analyses very difficult. It is therefore essential to develop new strategies for efficient sampling of this large sequence space, in order to reduce the number of sequences to be processed. At the same time, it is important to retain the most pertinent sequences for structural and functional studies.

RESULTS

An exhaustive analysis on a large scale test set (284 protein families) was performed to compare the efficiency of four different sampling methods aimed at selecting the most pertinent sequences. These four methods sample the proteins detected by BlastP searches and can be divided into two categories: two customisable methods where the user defines either the maximal number or the percentage of sequences to be selected; two automatic methods in which the number of sequences selected is determined by the program. We focused our analysis on the potential information content of the sampled sets of sequences using multiple alignment of complete sequences as the main validation tool. The study considered two criteria: the total number of sequences in BlastP and their associated E-values. The subsequent analyses investigated the influence of the sampling methods on the E-value distributions, the sequence coverage, the final multiple alignment quality and the active site characterisation at various residue conservation thresholds as a function of these criteria.

CONCLUSION

The comparative analysis of the four sampling methods allows us to propose a suitable sampling strategy that significantly reduces the number of homologous sequences required for alignment, while at the same time maintaining the relevant information concerning the active site residues.

MAGOS: multiple alignment and modelling server

MAGOS is a web server allowing automated protein modelling coupled to the creation of a hierarchical and annotated multiple alignment of complete sequences. MAGOS is designed for an interactive approach of structural information within the framework of the evolutionary relevance of mined and predicted sequence information.

AVAILABILITY

The web server is freely available at http://pig-pbil.ibcp.fr/magos.

MACSIMS: multiple alignment of complete sequences information management system

BACKGROUND

In the post-genomic era, systems-level studies are being performed that seek to explain complex biological systems by integrating diverse resources from fields such as genomics, proteomics or transcriptomics. New information management systems are now needed for the collection, validation and analysis of the vast amount of heterogeneous data available. Multiple alignments of complete sequences provide an ideal environment for the integration of this information in the context of the protein family.

RESULTS

MACSIMS is a multiple alignment-based information management program that combines the advantages of both knowledge-based and ab initio sequence analysis methods. Structural and functional information is retrieved automatically from the public databases. In the multiple alignment, homologous regions are identified and the retrieved data is evaluated and propagated from known to unknown sequences with these reliable regions. In a large-scale evaluation, the specificity of the propagated sequence features is estimated to be >99%, i.e. very few false positive predictions are made. MACSIMS is then used to characterise mutations in a test set of 100 proteins that are known to be involved in human genetic diseases. The number of sequence features associated with these proteins was increased by 60%, compared to the features available in the public databases. An XML format output file allows automatic parsing of the MACSIM results, while a graphical display using the JalView program allows manual analysis.

CONCLUSION

MACSIMS is a new information management system that incorporates detailed analyses of protein families at the structural, functional and evolutionary levels. MACSIMS thus provides a unique environment that facilitates knowledge extraction and the presentation of the most pertinent information to the biologist. A web server and the source code are available at http://bips.u-strasbg.fr/MACSIMS/.

Structural genomics of eukaryotic targets at a laboratory scale

Structural genomics programs are distributed worldwide and funded by large institutions such as the NIH in United-States, the RIKEN in Japan or the European Commission through the SPINE network in Europe. Such initiatives, essentially managed by large consortia, led to technology and method developments at the different steps required to produce biological samples compatible with structural studies. Besides specific applications, method developments resulted mainly upon miniaturization and parallelization. The challenge that academic laboratories faces to pursue structural genomics programs is to produce, at a higher rate, protein samples. The Structural Biology and Genomics Department (IGBMC - Illkirch - France) is implicated in a structural genomics program of high eukaryotes whose goal is solving crystal structures of proteins and their complexes (including large complexes) related to human health and biotechnology. To achieve such a challenging goal, the Department has established a medium-throughput pipeline for producing protein samples suitable for structural biology studies. Here, we describe the setting up of our initiative from cloning to crystallization and we demonstrate that structural genomics may be manageable by academic laboratories by strategic investments in robotic and by adapting classical bench protocols and new developments, in particular in the field of protein expression, to parallelization.

Double hexameric ring assembly of the type III protein translocase ATPase HrcN

The specialized type III secretion (T3S) apparatus of pathogenic and symbiotic Gram-negative bacteria comprises a complex transmembrane organelle and an ATPase homologous to the F1-ATPase beta subunit. The T3S ATPase HrcN of Pseudomonas syringae associates with the inner membrane, and its ATP hydrolytic activity is stimulated by dodecamerization. The structure of dodecameric HrcN (HrcN12) determined to 1.6 nm by cryo-electron microscopy is presented. HrcN12 comprises two hexameric rings that are probably stacked face-to-face by the association of their C-terminal domains. It is 11.5 +/- 1.0 nm in diameter, 12.0 +/- 2.0 nm high and has a 2.0-3.8 nm wide inner channel. This structure is compared to a homology model based on the structure of the F1-beta-ATPase. A model for its incorporation within the T3S apparatus is presented.

Distribution of cannabinoid receptor 1 in the CNS of zebrafish

The cannabinoid receptor 1 (Cb1) mediates the psychoactive effect of marijuana. In mammals, there is abundant evidence advocating the importance of cannabinoid signaling; activation of Cb1 exerts diverse functions, chiefly by its ability to modulate neurotransmission. Thus, much attention has been devoted to understand its role in health and disease and to evaluate its therapeutic potential. Here, we have cloned zebrafish cb1 and investigated its expression in developing and adult zebrafish brain. Sequence analysis showed that there is a high degree of conservation, especially in residues demonstrated to be critical for function in mammals. In situ hybridization revealed that zebrafish cb1 appears first in the preoptic area at 24 hours post-fertilization. Subsequently, transcripts are detected in the dorsal telencephalon, hypothalamus, pretectum and torus longitudinalis. A similar pattern of expression is recapitulated in the adult brain. While cb1 is intensively stained in the medial zone of the dorsal telencephalon, expression elsewhere is weak by comparison. In particular, localization of cb1 in the telencephalic periventricular matrix is suggestive of the involvement of Cb1 in neurogenesis, bearing strong resemblance in terms of expression and function to the proliferative mammalian hippocampal formation. In addition, a gradient-like expression of cb1 is detected in the torus longitudinalis, a teleost specific neural tissue. In relation to dopaminergic neurons in the diencephalic posterior tuberculum (considered to be the teleostean homologue of the mammalian midbrain dopaminergic system), both cb1 and tyrosine hydroxylase-expressing cells occupy non-overlapping domains. However there is evidence that they are co-localized in the caudal zone of the hypothalamus, implying a direct modulation of dopamine release in this particular region. Collectively, our data indicate the propensity of zebrafish cb1 to participate in multiple neurological processes.

vALId: validation of protein sequence quality based on multiple alignment data

The validation of sequences is essential to perform accurate phylogeny and structure/function analysis. However among the thousands of protein sequences available in the public databases, most have been predicted in silico and have not systematically undergone a quality verification. It has recently become evident that they often contain sequence errors. To address the problem of automatic protein quality control, we have developed vALId, an interactive web interfaced software. Taking advantage of high quality multiple alignments of complete protein sequences (MACS), vALId first warns about the presence of suspicious insertions, deletions (indels) and divergent segments, and second, proposes corrections based on transcripts and genome contigs. In a first evaluation test, hundreds of indels and divergent segments were randomly generated in a manually refined MACS. The sensitivity (Sn) and specificity (Sp) of indel detection were excellent (0.96) while the mean Sn(0.49) and Sp(0.56) of divergent segment delineation depended on the percent identity between sequence neighbors. In a second test, 6195 sequences in 100 MACS corresponding to different functional and structural protein families were analyzed. 65% of the sequences were in silico predictions and 44% of eukaryote predicted proteins were partially incorrect with at least one suspicious indel or divergent segment.

Conformational transition of initiation factor 2 from the GTP- to GDP-bound state visualized on the ribosome

Initiation of protein synthesis is a universally conserved event that requires initiation factors IF1, IF2 and IF3 in prokaryotes. IF2 is a GTPase essential for binding initiator transfer RNA to the 30S ribosomal subunit and recruiting the 50S subunit into the 70S initiation complex. We present two cryo-EM structures of the assembled 70S initiation complex comprising mRNA, fMet-tRNA(fMet) and IF2 with either a non-hydrolyzable GTP analog or GDP. Transition from the GTP-bound to the GDP-bound state involves substantial conformational changes of IF2 and of the entire ribosome. In the GTP analog-bound state, IF2 interacts mostly with the 30S subunit and extends to the initiator tRNA in the peptidyl (P) site, whereas in the GDP-bound state IF2 steps back and adopts a 'ready-to-leave' conformation. Our data also provide insights into the molecular mechanism guiding release of IF1 and IF3.

Cloning, purification and crystallization of a Walker-type Pyrococcus abyssi ATPase family member

Several ATPase proteins play essential roles in the initiation of chromosomal DNA replication in archaea. Walker-type ATPases are defined by their conserved Walker A and B motifs, which are associated with nucleotide binding and ATP hydrolysis. A family of 28 ATPase proteins with non-canonical Walker A sequences has been identified by a bioinformatics study of comparative genomics in Pyrococcus genomes. A high-throughput structural study on P. abyssi has been started in order to establish the structure of these proteins. 16 genes have been cloned and characterized. Six out of the seven soluble constructs were purified in Escherichia coli and one of them, PABY2304, has been crystallized. X-ray diffraction data were collected from selenomethionine-derivative crystals using synchrotron radiation. The crystals belong to the orthorhombic space group C2, with unit-cell parameters a = 79.41, b = 48.63, c = 108.77 A, and diffract to beyond 2.6 A resolution.

Cloning, expression, and characterization of pollen allergens from Humulus scandens (Lour) Merr and Ambrosia artemisiifolia L

AIM

To clone the pollen allergen genes in Humulus scandens (Lour) Merr (LvCao in Chinese) and short ragweed (Ambrosia artemisiifolia L) for recombinant allergen production and immunotherapy.

METHODS

The allergen genes were selectively amplified in the weed pollen cDNA pool by using a special PCR profile, with the primers designed by a modeling procedure. Following truncated gene cloning and confirmation of the pollen source, unknown 3'cDNA ends were identified by using the 3'-RACE method. The gene function conferred by the full-length coding region was evaluated by a homologue search in the GenBank database. Recombinant proteins expressed in Escherichia coli pET-44 RosettaBlue cells were subsequently characterized by N-terminal end sequencing, IgE binding, and cross-reactivity.

RESULTS

Three full-length cDNAs were obtained in each weed. Multiple alignment analysis revealed that the deduced amino acid sequences were 83% identical to each other and 56%-90% identical to panallergen profilins from other species. Five recombinant proteins were abundantly expressed in non-fusion forms and were confirmed by using the N-terminal end sequence identity. Sera from patients who were allergic to A artemisiifolia reacted not only with rAmb a 8(D03) derived from A artemisiifolia, but also with recombinant protein rHum s 1(LCM9) derived from H scandens, which confirmed the allergenicity and cross-reactivity of the recombinant proteins from the 2 sources. Comparison of the degenerate primers used for truncated gene cloning with the full-length cDNA demonstrated that alternative nucleotide degeneracy occurred.

CONCLUSION

This study demonstrates a useful method for cloning homologous allergen genes across different species, particularly for little-studied species. The recombinant allergens obtained might be useful for the immunotherapeutic treatment of H scandens and/or A artemisiifolia pollen allergies.

Structural Basis for tRNA-Dependent Amidotransferase Function

Besides direct charging of tRNAs by aminoacyl-tRNA synthetases, indirect routes also ensure attachment of some amino acids onto tRNA. Such routes may explain how new amino acids entered into protein synthesis. In archaea and in most bacteria, tRNA(Gln) is first misaminoacylated by glutamyl-tRNA synthetase. Glu-tRNA(Gln) is then matured into Gln-tRNA(Gln) by a tRNA-dependent amidotransferase. We report the structure of a tRNA-dependent amidotransferase-that of GatDE from Pyrococcus abyssi. The 3.0 A resolution crystal structure shows a tetramer with two GatD molecules as the core and two GatE molecules at the periphery. The fold of GatE cannot be related to that of any tRNA binding enzyme. The ammonium donor site on GatD and the tRNA site on GatE are markedly distant. Comparison of GatD and L-asparaginase structures shows how the motion of a beta hairpin region containing a crucial catalytic threonine may control the overall reaction cycle of GatDE.

Sequence and comparative genomic analysis of actin-related proteins

Actin-related proteins (ARPs) are key players in cytoskeleton activities and nuclear functions. Two complexes, ARP2/3 and ARP1/11, also known as dynactin, are implicated in actin dynamics and in microtubule-based trafficking, respectively. ARP4 to ARP9 are components of many chromatin-modulating complexes. Conventional actins and ARPs codefine a large family of homologous proteins, the actin superfamily, with a tertiary structure known as the actin fold. Because ARPs and actin share high sequence conservation, clear family definition requires distinct features to easily and systematically identify each subfamily. In this study we performed an in depth sequence and comparative genomic analysis of ARP subfamilies. A high-quality multiple alignment of approximately 700 complete protein sequences homologous to actin, including 148 ARP sequences, allowed us to extend the ARP classification to new organisms. Sequence alignments revealed conserved residues, motifs, and inserted sequence signatures to define each ARP subfamily. These discriminative characteristics allowed us to develop ARPAnno (http://bips.u-strasbg.fr/ARPAnno), a new web server dedicated to the annotation of ARP sequences. Analyses of sequence conservation among actins and ARPs highlight part of the actin fold and suggest interactions between ARPs and actin-binding proteins. Finally, analysis of ARP distribution across eukaryotic phyla emphasizes the central importance of nuclear ARPs, particularly the multifunctional ARP4.

An unusual primary sigma factor in the Bacteroidetes phylum

The presence of housekeeping gene promoters with a unique consensus sequence in Bacteroides fragilis, previously described by Bayley et al. (2000, FEMS Microbiol Lett 193: 149-154), suggested the existence of a particular primary sigma factor. The single rpoD-like gene observed in the B. fragilis genome, and similarly in those of other members of the Bacteroidetes phylum, was found to be essential. It encodes a protein, sigma(ABfr), of only 32.7 kDa that is produced with equal abundance during all phases of growth and was concluded to be the primary sigma factor. sigma(ABfr) and its orthologues in the Bacteroidetes are unusual primary sigma factors in that they lack region 1.1, have a unique signature made up of 29 strictly identical amino acids and are the only RpoD factors that cluster with the RpoS factors. Although binding to the Escherichia coli core RNA polymerase, sigma(ABfr) does not support transcription initiation from any promoter when it is part of the heterologous holoenzyme, while in the reconstituted homologous holoenzyme it does so only from typical B. fragilis, including rrs, promoters but not from the lacUV5 or RNA I promoters.

Identifying sigma factors in Mycobacterium smegmatis by comparative genomic analysis

Mycobacterium smegmatis is a saprophytic species that has been used for 15 years as a model to perform heterologous regulation and virulence studies of Mycobacterium tuberculosis. Members of the extracytoplasmic sigma factors family, which are required for adaptive responses to various environmental stresses, are responsible for some of the virulence traits of M. tuberculosis. A bioinformatic search on the genome of M. smegmatis has predicted the existence of 26 sigma factors, which is twice the number that are present in M. tuberculosis. A phylogenetic analysis has shown that despite this high number of sigma factors the orthologs of the genes sigC, sigI and sigK of M. tuberculosis are absent in the M. smegmatis genome. Several sigma factors are specific for M. smegmatis, with a special enrichment in the sigH and, to a lesser extent, in the sigJ and sigL subfamily, pinpointing the potential variability of the repertoire of adaptive response in this saprophytic species.

MAO: a Multiple Alignment Ontology for nucleic acid and protein sequences

The application of high-throughput techniques such as genomics, proteomics or transcriptomics means that vast amounts of heterogeneous data are now available in the public databases. Bioinformatics is responding to the challenge with new integrated management systems for data collection, validation and analysis. Multiple alignments of genomic and protein sequences provide an ideal environment for the integration of this mass of information. In the context of the sequence family, structural and functional data can be evaluated and propagated from known to unknown sequences. However, effective integration is being hindered by syntactic and semantic differences between the different data resources and the alignment techniques employed. One solution to this problem is the development of an ontology that systematically defines the terms used in a specific domain. Ontologies are used to share data from different resources, to automatically analyse information and to represent domain knowledge for non-experts. Here, we present MAO, a new ontology for multiple alignments of nucleic and protein sequences. MAO is designed to improve interoperation and data sharing between different alignment protocols for the construction of a high quality, reliable multiple alignment in order to facilitate knowledge extraction and the presentation of the most pertinent information to the biologist.

GOAnno: GO annotation based on multiple alignment

GOAnno is a web tool that automatically annotates proteins according to the Gene Ontology (GO) using evolutionary information available in hierarchized multiple alignments. GO terms present in the aligned functional subfamily can be cross-validated and propagated to obtain highly reliable predicted GO annotation based on the GOAnno algorithm.

AVAILABILITY

The web tool and a reduced version for local installation are freely available at http://igbmc.u-strasbg.fr/GOAnno/GOAnno.html

SUPPLEMENTARY INFORMATION

The website supplies a detailed explanation and illustration of the algorithm at http://igbmc.u-strasbg.fr/GOAnno/GOAnnoHelp.html.

Crystal structures of Weissella viridescens FemX and its complex with UDP-MurNAc-pentapeptide: insights into FemABX family substrates recognition

Members of the FemABX protein family are novel therapeutic targets, as they are involved in the synthesis of the bacterial cell wall. They catalyze the addition of amino acid(s) on the peptidoglycan precursor using aminoacylated tRNA as a substrate. We report here the high-resolution structure of Weissella viridescens L-alanine transferase FemX and its complex with the UDP-MurNAc-pentapeptide. This is the first structure example of a FemABX family member that does not possess a coiled-coil domain. FemX consists of two structurally equivalent domains, separated by a cleft containing the binding site of the UDP-MurNAc-pentapeptide and a long channel that traverses one of the two domains. Our structural studies bring new insights into the evolution of the FemABX and the related GNAT superfamilies, shed light on the recognition site of the aminoacylated tRNA in Fem proteins, and allowed manual docking of the acceptor end of the alanyl-tRNAAla.

Association of the breast cancer protein MLN51 with the exon junction complex via its speckle localizer and RNA binding module

MLN51 is a nucleocytoplasmic shuttling protein that is overexpressed in breast cancer. The function of MLN51 in mammals remains elusive. Its fly homolog, named barentsz, as well as the proteins mago nashi and tsunagi have been shown to be required for proper oskar mRNA localization to the posterior pole of the oocyte. Magoh and Y14, the human homologs of mago nashi and tsunagi, are core components of the exon junction complex (EJC). The EJC is assembled on spliced mRNAs and plays important roles in post-splicing events including mRNA export, nonsense-mediated mRNA decay, and translation. In the present study, we show that human MLN51 is an RNA-binding protein present in ribonucleo-protein complexes. By co-immunoprecipitation assays, endogenous MLN51 protein is found to be associated with EJC components, including Magoh, Y14, and NFX1/TAP, and subcellular localization studies indicate that MLN51 transiently co-localizes with Magoh in nuclear speckles. Moreover, we demonstrate that MLN51 specifically associates with spliced mRNAs in co-precipitation experiments, both in the nucleus and in the cytoplasm, at the position where the EJC is deposited. Most interesting, we have identified a region within MLN51 sufficient to bind RNA, to interact with Magoh and spliced mRNA, and to address the protein to nuclear speckles. This conserved region of MLN51 was therefore named SELOR for speckle localizer and RNA binding module. Altogether our data demonstrate that MLN51 associates with EJC in the nucleus and remains stably associated with mRNA in the cytoplasm, suggesting that its overexpression might alter mRNA metabolism in cancer.

Visualization of release factor 3 on the ribosome during termination of protein synthesis

Termination of protein synthesis by the ribosome requires two release factor (RF) classes. The class II RF3 is a GTPase that removes class I RFs (RF1 or RF2) from the ribosome after release of the nascent polypeptide. RF3 in the GDP state binds to the ribosomal class I RF complex, followed by an exchange of GDP for GTP and release of the class I RF. As GTP hydrolysis triggers release of RF3 (ref. 4), we trapped RF3 on Escherichia coli ribosomes using a nonhydrolysable GTP analogue. Here we show by cryo-electron microscopy that the complex can adopt two different conformational states. In 'state 1', RF3 is pre-bound to the ribosome, whereas in 'state 2' RF3 contacts the ribosome GTPase centre. The transfer RNA molecule translocates from the peptidyl site in state 1 to the exit site in state 2. This translocation is associated with a large conformational rearrangement of the ribosome. Because state 1 seems able to accommodate simultaneously both RF3 and RF2, whose position is known from previous studies, we can infer the release mechanism of class I RFs.

LEON: multiple aLignment Evaluation Of Neighbours

Sequence alignments are fundamental to a wide range of applications, including database searching, functional residue identification and structure prediction techniques. These applications predict or propagate structural/functional/evolutionary information based on a presumed homology between the aligned sequences. If the initial hypothesis of homology is wrong, no subsequent application, however sophisticated, can be expected to yield accurate results. Here we present a novel method, LEON, to predict homology between proteins based on a multiple alignment of complete sequences (MACS). In MACS, weak signals from distantly related proteins can be considered in the overall context of the family. Intermediate sequences and the combination of individual weak matches are used to increase the significance of low-scoring regions. Residue composition is also taken into account by incorporation of several existing methods for the detection of compositionally biased sequence segments. The accuracy and reliability of the predictions is demonstrated in large-scale comparisons with structural and sequence family databases, where the specificity was shown to be >99% and the sensitivity was estimated to be approximately 76%. LEON can thus be used to reliably identify the complex relationships between large multidomain proteins and should be useful for automatic high-throughput genome annotations, 2D/3D structure predictions, protein-protein interaction predictions etc.

The 'PINIT' motif, of a newly identified conserved domain of the PIAS protein family, is essential for nuclear retention of PIAS3L

PIAS proteins, cytokine-dependent STAT-associated repressors, exhibit intrinsic E3-type SUMO ligase activities and form a family of transcriptional modulators. Three conserved domains have been identified so far in this protein family, the SAP box, the MIZ-Zn finger/RING module and the acidic C-terminal domain, which are essential for protein interactions, DNA binding or SUMO ligase activity. We have identified a novel conserved domain of 180 residues in PIAS proteins and shown that its 'PINIT' motif as well as other conserved motifs (in the SAP box and in the RING domain) are independently involved in nuclear retention of PIAS3L, the long form of PIAS3, that we have characterized in mouse embryonic stem cells.