A reinforced merging methodology for mapping unique peptide motifs in members of protein families

Background  

Members of a protein family often have highly conserved sequences; most of these sequences carry identical biological functions and possess similar three-dimensional (3-D) structures. However, enzymes with high sequence identity may acquire differential functions other than the common catalytic ability. It is probable that each of their variable regions consists of a unique peptide motif (UPM), which selectively interacts with other cellular proteins, rendering additional biological activities. The ability to identify and localize such UPMs is paramount in recognizing the characteristic role of each member of a protein family.     

Human RNase 7: a new cationic ribonuclease of the RNase A superfamily

Here we report on the expression and function of RNase 7, one of the final RNase A superfamily ribonucleases identified in the human genome sequence. The human RNase 7 gene is expressed in various somatic tissues including the liver, kidney, skeletal muscle and heart. Recombinant RNase 7 is ribonucleolytically active against yeast tRNA, as expected from the presence of eight conserved cysteines and the catalytic histidine–lysine– histidine triad which are signature motifs of this superfamily. The protein is atypically cationic with an isoelectric point (pI) of 10.5. Expression of recombinant RNase 7 in Escherichia coli completely inhibits the growth of the host bacteria, similar to what has been observed for the cationic RNase, eosinophil cationic protein (ECP/RNase 3, pI 11.4). An in vitro assay demonstrates dose-dependent cytotoxicity of RNase 7 against bacteria E.coli, Pseudomonas aeruginosa and Staphylococcus aureus. While RNase 7 and ECP/RNase 3 are both cationic and share this particular aspect of functional similarity, their protein sequence identity is only 40%. Of particular interest, ECP/RNase 3’s cationicity is based on an (over)abundance of arginine residues, whereas RNase 7 includes an excess of lysine. This difference, in conjunction with the independent origins and different expression patterns, suggests that RNase 7 and ECP/RNase 3 may have been recruited to target different pathogens in vivo, if their physiological functions are indeed host defenses.     

Getting more from less: algorithms for rapid protein identification with multiple short peptide sequences

We describe two novel sequence similarity search algorithms, FASTS and FASTF, that use multiple short peptide sequences to identify homologous sequences in protein or DNA databases. FASTS searches with peptide sequences of unknown order, as obtained by mass spectrometry-based sequencing, evaluating all possible arrangements of the peptides. FASTF searches with mixed peptide sequences, as generated by Edman sequencing of unseparated mixtures of peptides. FASTF deconvolutes the mixture, using a greedy heuristic that allows rapid identification of high scoring alignments while reducing the total number of explored alternatives. Both algorithms use the heuristic FASTA comparison strategy to accelerate the search but use alignment probability, rather than similarity score, as the criterion for alignment optimality. Statistical estimates are calculated using an empirical correction to a theoretical probability. These calculated estimates were accurate within a factor of 10 for FASTS and 1000 for FASTF on our test dataset. FASTS requires only 15-20 total residues in three or four peptides to robustly identify homologues sharing 50% or greater protein sequence identity. FASTF requires about 25% more sequence data than FASTS for equivalent sensitivity, but additional sequence data are usually available from mixed Edman experiments. Thus, both algorithms can identify homologues that diverged 100 to 500 million years ago, allowing proteomic identification from organisms whose genomes have not been sequenced.     

Analysis of the resolution limitations of peptide identification algorithms

Proteome identification using peptide-centric proteomics techniques is a routinely used analysis technique. One of the most powerful and popular methods for the identification of peptides from MS/MS spectra is protein database matching using search engines. Significance thresholding through false discovery rate (FDR) estimation by target/decoy searches is used to ensure the retention of predominantly confident assignments of MS/MS spectra to peptides. However, shortcomings have become apparent when such decoy searches are used to estimate the FDR. To study these shortcomings, we here introduce a novel kind of decoy database that contains isobaric mutated versions of the peptides that were identified in the original search. Because of the supervised way in which the entrapment sequences are generated, we call this a directed decoy database. Since the peptides found in our directed decoy database are thus specifically designed to look quite similar to the forward identifications, the limitations of the existing search algorithms in making correct calls in such strongly confusing situations can be analyzed. Interestingly, for the vast majority of confidently identified peptide identifications, a directed decoy peptide-to-spectrum match can be found that has a better or equal match score than the forward match score, highlighting an important issue in the interpretation of peptide identifications in present-day high-throughput proteomics.     

A data bank merging related protein structures and sequences.

A data collection which merges protein structural and sequence information is described. Structural superpositions amongst proteins with similar main-chain fold were performed or collected from the literature. Sequences taken from the protein primary structure databases were associated with the multiple structural alignments providing they were at least 50% homologous in residue identity to one of the structural sequences and at least 50% of the structural sequence residues were alignable. Such restrictions allow reasonable confidence that the primary sequences share the conformation of the tertiary structural templates, except in the less conserved loop regions. Multiple structural superpositions were collected for 38 familial groups containing a total of 209 tertiary structures; 45 structures had no superposable mates and were used individually. Other information is also provided as main-chain and side-chain conformational angles, secondary structural assignments and the like. Wedding the primary and tertiary structural data resulted in an 8-fold increase of data bank sequence entries over those associated with the known three-dimensional architectures alone.     

Angiogenin is a cytotoxic, tRNA-specific ribonuclease in the RNase A superfamily.

Angiogenin is a 14.4-kDa human plasma protein with 65% homology to RNase A that retains the key active site residues and three of the four RNase A disulfide bonds. We demonstrate that recombinant angiogenin functions as a cytotoxic tRNA-specific RNase in cell-free lysates and when injected into Xenopus oocytes. Inhibition of protein synthesis by angiogenin correlates with degradation of endogenous oocyte tRNA. Exogenous, radiolabeled tRNA is also hydrolyzed by angiogenin, whereas oocyte rRNA and mRNA are not detectably degraded by angiogenin. Protein synthesis was restored to angiogenin-injected oocytes by injecting the RNase inhibitor RNasin plus total Xenopus or calf liver tRNAs, thereby demonstrating that the tRNA degradation induced by angiogenin was the sole cause of cytotoxicity. A similar tRNA-reversible inhibition of protein synthesis was seen in rabbit reticulocyte lysates. Angiogenin therefore appears to be a specific cellular tRNase, whereas five homologues in the RNase A superfamily lack angiogenin's specificity for tRNA. One of these homologues purified from human eosinophils, eosinophil-derived neurotoxin, nonspecifically degrades oocyte RNA similar to RNase A and is also cytotoxic at very low concentrations.     

Seminal RNase: a unique member of the ribonuclease superfamily

The RNase found in bull semen, although a member of the mammalian superfamily of ribonucleases, possesses some unusual properties. Besides its unique structure and enzymic properties, it displays antispermatogenic, antitumor and immunosuppressive activities. Seminal RNase belongs to an interesting group of RNases, the RISBASES (RIbonucleases with Special, i.e. non catalytic, Biological Actions) other members of which include angiogenin, selectively neurotoxic RNases, a lectin and the self-incompatibility factors from a flowering plant.     

Zebrafish ribonucleases are bactericidal: implications for the origin of the vertebrate RNase A superfamily

Understanding the evolutionary origin of the ribonuclease (RNase) A superfamily is of great interest because the superfamily is the sole vertebrate-specific enzyme family known to date. Although mammalian RNases have a diverse array of biochemical and physiological functions, the original function of the superfamily at its birth is enigmatic. Such information may be obtained by studying basal lineages of the vertebrate phylogeny and is necessary for discerning how and why this superfamily originated. Here, we clone and characterize 3 RNase genes from the zebrafish, the most basal vertebrate examined for RNases. We report 1) that all the 3 zebrafish RNases are ribonucleolytically active, with one of them having an RNase activity comparable to that of bovine RNase A, the prototype of the superfamily; 2) that 2 zebrafish RNases have prominent expressions in adult liver and gut, whereas the 3rd is expressed in adult eye and heart; and 3) that all 3 RNases have antibacterial activities in vitro. These results, together with the presence of antibacterial and/or antiviral activities in multiple distantly related mammalian RNases, strongly suggest that the superfamily started as a host-defense mechanism in vertebrate evolution.     

LOC 390443 (RNase 9) on chromosome 14q11.2 is related to the RNase A superfamily and contains a unique amino-terminal preproteinlike sequence

A new member of the human RNase A superfamily is reported. Identified in the human genome assembly as LOC 390443, this locus is located 128 kb telomeric to the established RNase A gene family cluster on chromosome 14q11.2. The amino acid sequence of this locus is sufficiently similar to the eight previously identified gene family members to warrant a designation as RNase 9. RNase 9 is expressed in a wide range of human tissues. In addition, a 30-amino acid sequence lying between a 26-amino acid putative signal peptide and the last 148 amino acids that align with the other RNases A is not seen in other members of the RNase A superfamily in any species. Nucleotide and amino acid sequences of RNase 9 in 13 nonhuman primate species were determined and indicate several conserved sites but, also, an excess of nonsynonymous substitutions, about one-third of which are radical substitutions. This suggests that RNase 9, similar to several other human RNases A, has been under diversifying selection in the primates. Data from the mouse and rat genomes indicate that RNase 9 is also present in rodents, thus making it older than most of the established members of the human RNase A superfamily. Many of the human RNases A have been shown to have antimicrobial, antiviral, or antiparasitic functions involved in host-defense mechanisms. The features of RNase 9 described here suggest that it, too, may be involved in host defense and that it, along with the rest of the superfamily, may prove to have played an important role in anthropoid evolution.     

Unique repetitive sequences of 170 bp inChlorella

Summary Cot analysis ofChlorella DNA revealed that the genome of the unicellular green alga contained a small amount of repetitive sequences (at most 15% of the total DNA). Short repetitive sequences (SRS) of 170 bp produced by enzymatic digestion of algal DNA with eitherHaeIII,HinfI, orPstI, were found by polyacrylamide gel electrophoresis, and their copy number was estimated to be a few hundred (about 2% of the total repetitive sequences). All three members showed high sequence homology and could be be unified into one family, HaeIII family. The family was divided further into two subfamilies,HinfI- (HaeIII-andHinfI-SRS) andPstI-(PstI-SRS) subfamilies, based on small sequence differences among the members. TheHaeIII family had characteristic structural features, including a considerable number of small unique sequence units (purine-CC) and both direct and inverted repeats, and were organized in tandem arrays in the genome.     

Phylogenetic study and identification of Vibrio splendidus-related strains based on gyrB gene sequences

Different strains related to Vibrio splendidus have been associated with infection of aquatic animals. An epidemiological study of V. splendidus strains associated with Crassostrea gigas mortalities demonstrated genetic diversity within this group and suggested its polyphyletic nature. Recently 4 species, V. lentus, V. chagasii, V. pomeroyi and V. kanaloae, phenotypically related to V. splendidus, have been described, although biochemical methods do not clearly discriminate species within this group. Here, we propose a polyphasic approach to investigate their taxonomic relationships. Phylogenetic analysis of V. splendidus-related strains was carried out using the nucleotide sequences of 16S ribosomal DNA (16S rDNA) and gyrase B subunit (gyrB) genes. Species delineation based on 16S rDNA-sequencing is limited because of divergence between cistrons, roughly equivalent to divergence between strains. Despite a high level of sequence similarity, strains were separated into 2 clades. In the phylogenetic tree constructed on the basis of gyrB gene sequences, strains were separated into 5 independent clusters containing V. splendidus, V. lentus, V. chagasii-type strains and a putative new genomic species. This phylogenetic grouping was almost congruent with that based on DNA-DNA hybridisation analysis. V. pomeroyi, V. kanaloae and V. tasmaniensis-type strains clustered together in a fifth clade. The gyrB gene-sequencing approach is discussed as an alternative for investigating the taxonomy of Vibrio species.     

A new superfamily of protein-O-fucosyltransferases,alpha2-fucosyltransferases,and alpha6-fucosyltransferases: phylogeny and identification of conserved peptide motifs

The presence of three conserved peptide motifs shared by alpha2-fucosyltransferases, alpha6-fucosyltransferases, the protein-O-fucosyltransferase family 1 (POFUT1) and a newly identified protein-O-fucosyltransferase family 2 (POFUT2), together with evidence that the present genes encoding for these enzymes have originated from a common ancestor by duplication and divergent evolution, suggests that they constitute a new superfamily of fucosyltransferases.     

Protein and peptide identification algorithms using MS for use in high-throughput, automated pipelines

Current proteomics experiments can generate vast quantities of data very quickly, but this has not been matched by data analysis capabilities. Although there have been a number of recent reviews covering various aspects of peptide and protein identification methods using MS, comparisons of which methods are either the most appropriate for, or the most effective at, their proposed tasks are not readily available. As the need for high-throughput, automated peptide and protein identification systems increases, the creators of such pipelines need to be able to choose algorithms that are going to perform well both in terms of accuracy and computational efficiency. This article therefore provides a review of the currently available core algorithms for PMF, database searching using MS/MS, sequence tag searches and de novo sequencing. We also assess the relative performances of a number of these algorithms. As there is limited reporting of such information in the literature, we conclude that there is a need for the adoption of a system of standardised reporting on the performance of new peptide and protein identification algorithms, based upon freely available datasets. We go on to present our initial suggestions for the format and content of these datasets.     

New light on bacterial carbonic anhydrases phylogeny based on the analysis of signal peptide sequences

Among protein families, carbonic anhydrases (CAs, EC 4.2.1.1) are metalloenzymes characterized by a common reaction mechanism in all life domains: the carbon dioxide hydration to bicarbonate and protons (CO₂+H₂O ? HCO₃^?+H⁺). Six genetically distinct CA families are known to date, the α-, β-, γ-, δ-, ζ- and η-CAs. The last CA class was recently discovered analyzing the amino acid sequences of CAs from Plasmodia. Bacteria encode for enzymes belonging to the α-, β-, and γ-CA classes and recently, phylogenetic analysis revealed an interesting relationship regarding the evolution of bacterial CA classes. This result evidenced that the three bacterial CA classes, in spite of the high level of the structural similarity, are evolutionarily distinct, but we noted that the primary structure of some β-CAs identified in the genome of Gram-negative bacteria present a pre-sequence of 18 or more amino acid residues at the N-terminal part. These observations and subsequent phylogenetic data presented here prompted us to propose that the β-CAs found in Gram-negative bacteria with a periplasmic space and characterized by the presence of a signal peptide might have a periplasmic localization and a role similar to that described previously for the α-CAs.     

A proposal on metrics for identification using nucleic acid sequences

Abstract The need is stressed for attempts to be made to permit diagnostic nucleic acid sequences to be used in a quantitative manner. Sequence differences or binding values should be converted to a distance measure and from this an ultrametric tree should be constructed. A single quantitative determination can yield considerable information about the likely identity of an unknown microorganism when the distance obtained from the sequence is compared with the tree. The concept is illustrated by hypothetical species and genus subsequences, and it is suitable both for successive use of hierarchical subsequences and for automated identification. It is pointed out that entirely specific subsequences for higher taxa may be difficult to discover. These principles will be useful for the future design of diagnostic sequences, including possible application to DNA-DNA pairing.     

A DNA microarray for identification of selected Korean birds based on mitochondrial cytochrome c oxidase I gene sequences

DNA barcoding with the gene encoding cytochrome c oxidase I (COI) in the mitochondrial genome has been proposed as a standard marker to identify and discover animal species. Some migratory wild birds are suspected of transmitting avian influenza and pose a threat to aircraft safety because of bird strikes. We have previously reported the COI gene sequences of 92 Korean bird species. In the present study, we developed a DNA microarray to identify 17 selected bird species on the basis of nucleotide diversity. We designed and synthesized 19 specific oligonucleotide probes; these probes were arrayed on a silylated glass slide. The length of the probes was 19-24 bps. The COI sequences amplified from the tissues of the selected birds were labeled with a fluorescent probe for microarray hybridization, and unique hybridization patterns were detected for each selected species. These patterns may be considered diagnostic patterns for species identification. This microarray system will provide a sensitive and a high-throughput method for identification of Korean birds.