Learning from Decoys to Improve the Sensitivity and Specificity of Proteomics Database Search Results

The statistical validation of database search results is a complex issue in bottom-up proteomics. The correct and incorrect peptide spectrum match (PSM) scores overlap significantly, making an accurate assessment of true peptide matches challenging. Since the complete separation between the true and false hits is practically never achieved, there is need for better methods and rescoring algorithms to improve upon the primary database search results. Here we describe the calibration and False Discovery Rate (FDR) estimation of database search scores through a dynamic FDR calculation method, FlexiFDR, which increases both the sensitivity and specificity of search results. Modelling a simple linear regression on the decoy hits for different charge states, the method maximized the number of true positives and reduced the number of false negatives in several standard datasets of varying complexity (18-mix, 49-mix, 200-mix) and few complex datasets (E. coli and Yeast) obtained from a wide variety of MS platforms. The net positive gain for correct spectral and peptide identifications was up to 14.81% and 6.2% respectively. The approach is applicable to different search methodologies- separate as well as concatenated database search, high mass accuracy, and semi-tryptic and modification searches. FlexiFDR was also applied to Mascot results and showed better performance than before. We have shown that appropriate threshold learnt from decoys, can be very effective in improving the database search results. FlexiFDR adapts itself to different instruments, data types and MS platforms. It learns from the decoy hits and sets a flexible threshold that automatically aligns itself to the underlying variables of data quality and size.     

Defining components of the chromosomal origin of replication of the hyperthermophilic archaeon Pyrococcus furiosus needed for construction of a stable replicating shuttle vector

We report the construction of a series of replicating shuttle vectors that consist of a low-copy-number cloning vector for Escherichia coli and functional components of the origin of replication (oriC) of the chromosome of the hyperthermophilic archaeon Pyrococcus furiosus. In the process of identifying the minimum replication origin sequence required for autonomous plasmid replication in P. furiosus, we discovered that several features of the origin predicted by bioinformatic analysis and in vitro binding studies were not essential for stable autonomous plasmid replication. A minimum region required to promote plasmid DNA replication was identified, and plasmids based on this sequence readily transformed P. furiosus. The plasmids replicated autonomously and existed in a single copy. In contrast to shuttle vectors based on a plasmid from the closely related hyperthermophile Pyrococcus abyssi for use in P. furiosus, plasmids based on the P. furiosus chromosomal origin were structurally unchanged after transformation and were stable without selection for more than 100 generations.     

Microbial Metalloproteomes Explored Using MIRAGE

Metalloproteomics is a rapidly developing field of science that involves the comprehensive analysis of all metal-containing or metal-binding proteins in a biological sample. The purpose of this review is to offer an overview of the research involving Metal Isotope native RadioAutography in Gel Electrophoresis (MIRAGE), a powerful new method to visualize and study the proteome of a particular metal ion. MIRAGE involves four steps: i) labelling of target proteins with a radioisotope; ii) separation of intact holo-proteins using native isoelectric focusing (1D) combined with Blue Native PAGE (2D); iii) spot visualization and quantification using autoradiography; and iv) protein identification by tandem mass spectrometry. MIRAGE Investigations of the soluble Cu, Zn, and Fe metalloproteomes of Escherichia coli, and of the soluble Mo and W proteomes of the hyperthermophilic archaeon Pyrococcus furiosus are reviewed.     

The GINS complex from Pyrococcus furiosus stimulates the MCM helicase activity

Pyrococcus furiosus, a hyperthermophilic Archaea, has homologs of the eukaryotic MCM (mini-chromosome maintenance) helicase and GINS complex. The MCM and GINS proteins are both essential factors to initiate DNA replication in eukaryotic cells. Many biochemical characterizations of the replication-related proteins have been reported, but it has not been proved that the homologs of each protein are also essential for replication in archaeal cells. Here, we demonstrated that the P. furiosus GINS complex interacts with P. furiosus MCM. A chromatin immunoprecipitation assay revealed that the GINS complex is detected preferentially at the oriC region on Pyrococcus chromosomal DNA during the exponential growth phase but not in the stationary phase. Furthermore, the GINS complex stimulates both the ATPase and DNA helicase activities of MCM in vitro. These results strongly suggest that the archaeal GINS is involved in both the initiation and elongation processes of DNA replication in P. furiosus, as observed in eukaryotic cells.     

The tungsten metallome of Pyrococcus furiosus

The tungsten metallome of the hyperthermophilic archaeon Pyrococcus furiosus has been investigated using electroanalytical metal analysis and native-native 2D-PAGE with the radioactive tungsten isotope (187)W (t(1/2) = 23.9 h). P. furiosus cells have an intracellular tungsten concentration of 29 μM, of which ca. 30% appears to be free tungsten, probably in the form of tungstate or polytungstates. The remaining 70% is bound by five different tungsten enzymes: formaldehyde ferredoxin oxidoreductase, aldehyde ferredoxin oxidoreductase, glyceraldehyde-3-phosphate ferredoxin oxidoreductase and the tungsten-containing oxidoreductases WOR4 and WOR5. The membrane proteome of P. furiosus is devoid of tungsten. The differential expression, as measured by the tungsten level, of the five soluble tungsten enzymes when the cells are subjected to a cold-shock shows a strong correlation with previously published DNA microarray analyses.     

A computational method for assessing peptide- identification reliability in tandem mass spectrometry analysis with SEQUEST

High-throughput protein identification in mass spectrometry is predominantly achieved by first identifying tryptic peptides by a database search and then by combining the peptide hits for protein identification. One of the popular tools used for the database search is SEQUEST. Peptide identification is carried out by selecting SEQUEST hits above a specified threshold, the value of which is typically chosen empirically in an attempt to separate true identifications from false ones. These SEQUEST scores are not normalized with respect to the composition, length and other parameters of the peptides. Furthermore, there is no rigorous reliability estimate assigned to the protein identifications derived from these scores. Hence, the interpretation of SEQUEST hits generally requires human involvement, making it difficult to scale up the identification process for genome-scale applications. To overcome these limitations, we have developed a method, which combines a neural network and a statistical model, for normalizing SEQUEST scores, and also for providing a reliability estimate for each SEQUEST hit. This method improves the sensitivity and specificity of peptide identification compared to the standard filtering procedure used in the SEQUEST package, and provides a basis for estimating the reliability of protein identifications.     

DNA protection mechanisms are not involved in the radioresistance of the hyperthermophilic archaea<Emphasis Type="Italic"> Pyrococcus abyssi</Emphasis> and<Emphasis Type="Italic"> P. furiosus</Emphasis>

Hyperthermophilic archaea of the genus Pyrococcus are resistant to gamma radiation, suggesting that efficient mechanisms for DNA repair exist in these organisms. To determine whether protective mechanisms might also be implicated in this radioresistance, we have estimated the linear density of DNA double-stranded breaks caused by gamma irradiation in the genomic DNA of two Pyrococcus species, using Escherichia coli and the radioresistant bacterium Deinococcus radiodurans as controls. The linear density of double-stranded breaks was essentially the same in all four microorganisms when irradiation was carried under similar anaerobic conditions, indicating that no specific DNA protection mechanisms exist in Pyrococcus species. Using one- and two-dimensional gel electrophoresis we compared the protein patterns from Pyrococcus abyssi and P. furiosus cells that had or had not been exposed to gamma rays. We did not detect any significant protein induction following DNA damage in either species.     

Integrated data management and validation platform for phosphorylated tandem mass spectrometry data

MS/MS is a widely used method for proteome‐wide analysis of protein expression and PTMs. The thousands of MS/MS spectra produced from a single experiment pose a major challenge for downstream analysis. Standard programs, such as MASCOT, provide peptide assignments for many of the spectra, including identification of PTM sites, but these results are plagued by false‐positive identifications. In phosphoproteomic experiments, only a single peptide assignment is typically available to support identification of each phosphorylation site, and hence minimizing false positives is critical. Thus, tedious manual validation is often required to increase confidence in the spectral assignments. We have developed phoMSVal, an open‐source platform for managing MS/MS data and automatically validating identified phosphopeptides. We tested five classification algorithms with 17 extracted features to separate correct peptide assignments from incorrect ones using over 2600 manually curated spectra. The naïve Bayes algorithm was among the best classifiers with an AUC value of 97% and PPV of 97% for phosphotyrosine data. This classifier required only three features to achieve a 76% decrease in false positives as compared with MASCOT while retaining 97% of true positives. This algorithm was able to classify an independent phosphoserine/threonine data set with AUC value of 93% and PPV of 91%, demonstrating the applicability of this method for all types of phospho‐MS/MS data. PhoMSVal is available at http://csbi.ltdk.helsinki.fi/phomsval .     

Stabilization of Taq DNA polymerase at high temperature by protein folding pathways from a hyperthermophilic archaeon, Pyrococcus furiosus

Pyrococcus furiosus, a hyperthermophilic archaeon growing optimally at 100 degrees C, encodes three protein chaperones, a small heat shock protein (sHsp), a prefoldin (Pfd), and a chaperonin (Cpn). In this study, we report that the passive chaperones sHsp and Pfd from P. furiosus can boost the protein refolding activity of the ATP-dependent Cpn from the same hyperthermophile. The thermo-stability of Taq polymerase was significantly improved by combinations of P. furiosus chaperones, showing ongoing protein folding activity at elevated temperatures and during thermal cycling. Based on these results, we propose that the protein folding apparatus in the hyperthermophilic archaeon, P. furiosus can be utilized to enhance the durability and cost effectiveness of high temperature biocatalysts.     

The crystal structure of a virus-like particle from the hyperthermophilic archaeon Pyrococcus furiosus provides insight into the evolution of viruses

Pyrococcus furiosus is a hyperthermophilic archaeal microorganism found near deep-sea thermal vents and its optimal growth temperature of 100 degrees C. Recently, a 38.8-kDa protein from P. furiosus DSM 3638 was isolated and characterized. Electron microscopy revealed that this protein aggregated as spheres of approximately 30 nm in diameter, which we designated P. furiosus virus-like particles (PfVs). X-ray crystallographic analysis at 3.6-A resolution revealed that each PfV consisted of 180 copies of the 38.8-kDa protein and retained T=3 icosahedral symmetry, as is often the case in spherical viruses. The total molecular mass of each particle was approximately 7 MDa. An examination of capsid structures suggested strong evolutionary links among PfV, tailed double-stranded DNA bacteriophages, and herpes viruses. The similar three-dimensional structures of the various coat proteins indicate that these viral capsids might have originated and evolved from a common ancestor. The structure of PfV provides a previously undescribed example of viral relationships across the three domains of life (Eukarya, Bacteria, and Archaea).     

A proposal to rename the hyperthermophile Pyrococcus woesei as Pyrococcus furiosus subsp. woesei

Pyrococcus species are hyperthermophilic members of the order Thermococcales, with optimal growth temperatures approaching 100 degrees C. All species grow heterotrophically and produce H2 or, in the presence of elemental sulfur (S(o)), H2S. Pyrococcus woesei and P. furiosus were isolated from marine sediments at the same Vulcano Island beach site and share many morphological and physiological characteristics. We report here that the rDNA operons of these strains have identical sequences, including their intergenic spacer regions and part of the 23S rRNA. Both species grow rapidly and produce H2 in the presence of 0.1% maltose and 10-100 microM sodium tungstate in S(o)-free medium. However, P. woesei shows more extensive autolysis than P. furiosus in the stationary phase. Pyrococcus furiosus and P. woesei share three closely related families of insertion sequences (ISs). A Southern blot performed with IS probes showed extensive colinearity between the genomes of P. woesei and P. furiosus. Cloning and sequencing of ISs that were in different contexts in P. woesei and P. furiosus revealed that the napA gene in P. woesei is disrupted by a type III IS element, whereas in P. furiosus, this gene is intact. A type I IS element, closely linked to the napA gene, was observed in the same context in both P. furiosus and P. woesei genomes. Our results suggest that the IS elements are implicated in genomic rearrangements and reshuffling in these closely related strains. We propose to rename P. woesei a subspecies of P. furiosus based on their identical rDNA operon sequences, many common IS elements that are shared genomic markers, and the observation that all P. woesei nucleotide sequences deposited in GenBank to date are > 99% identical to P. furiosus sequences.     

Increasing specificity in high-throughput yeast two-hybrid experiments

Since its inception, the yeast two-hybrid (Y2H) system has proven to be an efficient system to identify novel protein-protein interactions. However, Y2H screens are sometimes criticized for generating high rates of false-positives. Minimizing false-positive interactions is especially important in proteome wide high-throughput (HT) Y2H. Here, we summarize various approaches that reduce false-positives in HT-Y2H projects. We evaluated the potential of examining putative positives after removing the prey encoding plasmid by negative selection. We found that this method reliably identifies false-positives caused by spontaneous conversion of baits into auto-activators and provides significant time-savings in HT screens. In addition, we present a method to eliminate an important source of false-positives: contaminating prey plasmids. Y2H interactors can be wrongly identified due to the presence of two or more different plasmids in the cells of a single yeast colony. Of these independent plasmids, only one encodes a genuine interactor. Contaminating plasmids are eliminated by extended culture of yeast cells under positive selection for the interaction, allowing the identification of the true interaction partner.     

Properties of the prolyl oligopeptidase homologue from Pyrococcus furiosus

Prolyl oligopeptidase (POP), the paradigm of a serine peptidase family, hydrolyses peptides, but not proteins. The thermophilic POP from Pyrococcus furiosus (Pfu) appeared to be an exception, since it hydrolysed large proteins. Here we demonstrate that the Pfu POP does not display appreciable activity against azocasein. The autolysis observed earlier was an artefact. We have also found that the pH-rate profile is different from that of the mammalian enzyme and the low pK(a) extracted from the curve represents the ionization of the catalytic histidine. We conclude that some oligopeptidases may be true endopeptidases, cleaving at disordered segments of proteins, but with very low efficacy.     

Prefractionation of proteome by liquid isoelectric focusing prior to two-dimensional liquid chromatography mass spectrometric identification

Due to the complexity of proteomes, developing methods of sample fractionation, separation, concentration, and detection have become urgent to the identification of large numbers of proteins, as well as the acquisition of those proteins in low abundance. In this work, liquid isoelectric focusing (LIEF) combined with 2D-LC-MS/MS was applied to the proteome of Saccharomyces cerevisiae. This yielded a total of 1795 proteins that were detected and identified by 30 fractions of protein prefractionation. Categorization of these hits demonstrated the ability of this technology to detect and identify proteins rarely seen in proteome analysis without protein fractionation. LIEF-2D-LC-MS/MS also produced improved resolution of low-abundance proteins. Furthermore, we analyzed the characteristics of proteins obtained by LIEF-2D-LC-MS/MS. 1103 proteins with CAI under 0.2 were identified, allowing us to specifically obtain detailed biochemical information on these kind proteins. It was observed that LIEF-2D-LC-MS/MS is useful for large-scale proteome analysis and may be specifically applied to systems with wide dynamic ranges.     

Production of a thermostable archaeal superoxide reductase in plant cells

Pyrococcus furiosus superoxide reductase (SOR) is a thermostable archaeal enzyme that reduces superoxide without producing oxygen. When produced as a fusion protein with the green fluorescent protein in plant cells, P. furiosus SOR is located in the cytosol and nucleus. The recombinant SOR enzyme retains its function and heat stability when assayed in vitro. Importantly, expressing SOR in plant cells enhances their survival at high temperature indicating that it functions in vivo. The archaeal SOR provides a novel mechanism to reduce superoxide and demonstrates the potential for using archaeal genes to alter eukaryotic metabolism.     

Electroanalytical determination of tungsten and molybdenum in proteins

Recent crystal structure determinations accelerated the progress in the biochemistry of tungsten-containing enzymes. In order to characterize these enzymes, a sensitive determination of this metal in protein-containing samples is necessary. An electroanalytical tungsten determination has successfully been adapted to determine the tungsten and molybdenum content in enzymes. The tungsten and molybdenum content can be measured simultaneously from 1 to 10 microg of purified protein with little or no sample handling. More crude protein samples require precipitation of interfering surface active material with 10% perchloric acid. This method affords the isolation of novel molybdenum- and tungsten-containing proteins via molybdenum and tungsten monitoring of column fractions, without using radioactive isotopes. A screening of soluble proteins from Pyrococcus furiosus for tungsten, using anion-exchange column chromatography to separate the proteins, has been performed. The three known tungsten-containing enzymes from P. furiosus were recovered with this screening.