Assessment of metabolome annotation quality: a method for evaluating the false discovery rate of elemental composition searches.

BACKGROUND: In metabolomics researches using mass spectrometry (MS), systematic searching of high-resolution mass data against compound databases is often the first step of metabolite annotation to determine elemental compositions possessing similar theoretical mass numbers. However, incorrect hits derived from errors in mass analyses will be included in the results of elemental composition searches. To assess the quality of peak annotation information, a novel methodology for false discovery rates (FDR) evaluation is presented in this study. Based on the FDR analyses, several aspects of an elemental composition search, including setting a threshold, estimating FDR, and the types of elemental composition databases most reliable for searching are discussed. METHODOLOGY/PRINCIPAL FINDINGS: The FDR can be determined from one measured value (i.e., the hit rate for search queries) and four parameters determined by Monte Carlo simulation. The results indicate that relatively high FDR values (30-50%) were obtained when searching time-of-flight (TOF)/MS data using the KNApSAcK and KEGG databases. In addition, searches against large all-in-one databases (e.g., PubChem) always produced unacceptable results (FDR >70%). The estimated FDRs suggest that the quality of search results can be improved not only by performing more accurate mass analysis but also by modifying the properties of the compound database. A theoretical analysis indicates that FDR could be improved by using compound database with smaller but higher completeness entries. CONCLUSIONS/SIGNIFICANCE: High accuracy mass analysis, such as Fourier transform (FT)-MS, is needed for reliable annotation (FDR <10%). In addition, a small, customized compound database is preferable for high-quality annotation of metabolome data.

Metabolomics approach for determining growth-specific metabolites based on Fourier transform ion cyclotron resonance mass spectrometry.

Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR/MS) is the best MS technology for obtaining exact mass measurements owing to its great resolution and accuracy, and several outstanding FT-ICR/MS-based metabolomics approaches have been reported. A reliable annotation scheme is needed to deal with direct-infusion FT-ICR/MS metabolic profiling. Correlation analyses can help us not only uncover relations between the ions but also annotate the ions originated from identical metabolites (metabolite derivative ions). In the present study, we propose a procedure for metabolite annotation on direct-infusion FT-ICR/MS by taking into consideration the classification of metabolite-derived ions using correlation analyses. Integrated analysis based on information of isotope relations, fragmentation patterns by MS/MS analysis, co-occurring metabolites, and database searches (KNApSAcK and KEGG) can make it possible to annotate ions as metabolites and estimate cellular conditions based on metabolite composition. A total of 220 detected ions were classified into 174 metabolite derivative groups and 72 ions were assigned to candidate metabolites in the present work. Finally, metabolic profiling has been able to distinguish between the growth stages with the aid of PCA. The constructed model using PLS regression for OD(600) values as a function of metabolic profiles is very useful for identifying to what degree the ions contribute to the growth stages. Ten phospholipids which largely influence the constructed model are highly abundant in the cells. Our analyses reveal that global modification of those phospholipids occurs as E. coli enters the stationary phase. Thus, the integrated approach involving correlation analyses, metabolic profiling, and database searching is efficient for high-throughput metabolomics.

Differential metabolomics unraveling light/dark regulation of metabolic activities in Arabidopsis cell culture.

Differential metabolomics based on a non-targeted FT-ICR/MS analysis demonstrated metabolite accumulation patterns reflecting light/dark conditions in Arabidopsis T87 cell culture. First, FT-ICR/MS data sets were converted into metabolome information using the Dr.DMASS software (http://kanaya.naist.jp/DrDMASS/). A quick search of a metabolite-species database, KNApSAcK (http://kanaya.naist.jp/KNApSAcK/), was implemented to assign metabolite candidates to each accurate MS data (<1 ppm) through the prediction of molecular formulas, and the candidate structures were further studied using MS/MS analyses. Specific metabolites representing the culture conditions included sugars, phenylpropanoid derivatives, flavonol aglycons, and a plastid nonmevalonate pathway intermediate. Transcriptomics data were obtained in parallel and analyzed using a transcriptome analysis tool, KaPPA-View (http://kpv.kazusa.or.jp/kappa-view/). The specific accumulation patterns of flavonol aglycons were in good agreement with the light/dark regulation of a cytochrome P450 gene, CYP75B, and the build-up of 2-C-methyl-D-erythritol 4-phosphate, a nonmevalonate pathway intermediate, in the light grown cells was also consistent with a gene expression profile. The differential metablomics scheme based on the FT-ICR/MS metabolomics can serve as an evaluation system of metabolic activities contributing to successful identification and proper manipulation of key enzymatic steps in metabolic engineering studies.

Clarification of pathway-specific inhibition by Fourier transform ion cyclotron resonance/mass spectrometry-based metabolic phenotyping studies.

We have developed a metabolic profiling scheme based on direct-infusion Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR/MS). The scheme consists of: (1) reproducible data collection under optimized FT-ICR/MS analytical conditions; (2) automatic mass-error correction and multivariate analyses for metabolome characterization using a newly developed metabolomics tool (DMASS software); (3) identification of marker metabolite candidates by searching a species-metabolite relationship database, KNApSAcK; and (4) structural analyses by an MS/MS method. The scheme was applied to metabolic phenotyping of Arabidopsis (Arabidopsis thaliana) seedlings treated with different herbicidal chemical classes for pathway-specific inhibitions. Arabidopsis extracts were directly infused into an electrospray ionization source on an FT-ICR/MS system. Acquired metabolomics data were comprised of mass-to-charge ratio values with ion intensity information subjected to principal component analysis, and metabolic phenotypes from the herbicide treatments were clearly differentiated from those of the herbicide-free treatment. From each herbicide treatment, candidate metabolites representing such metabolic phenotypes were found through the KNApSAcK database search. The database search and MS/MS analyses suggested dose-dependent accumulation patterns of specific metabolites including several flavonoid glycosides. The metabolic phenotyping scheme on the basis of FT-ICR/MS coupled with the DMASS program is discussed as a general tool for high throughput metabolic phenotyping studies.

Development and implementation of an algorithm for detection of protein complexes in large interaction networks.

BACKGROUND: After complete sequencing of a number of genomes the focus has now turned to proteomics. Advanced proteomics technologies such as two-hybrid assay, mass spectrometry etc. are producing huge data sets of protein-protein interactions which can be portrayed as networks, and one of the burning issues is to find protein complexes in such networks. The enormous size of protein-protein interaction (PPI) networks warrants development of efficient computational methods for extraction of significant complexes. RESULTS: This paper presents an algorithm for detection of protein complexes in large interaction networks. In a PPI network, a node represents a protein and an edge represents an interaction. The input to the algorithm is the associated matrix of an interaction network and the outputs are protein complexes. The complexes are determined by way of finding clusters, i. e. the densely connected regions in the network. We also show and analyze some protein complexes generated by the proposed algorithm from typical PPI networks of Escherichia coli and Saccharomyces cerevisiae. A comparison between a PPI and a random network is also performed in the context of the proposed algorithm. CONCLUSION: The proposed algorithm makes it possible to detect clusters of proteins in PPI networks which mostly represent molecular biological functional units. Therefore, protein complexes determined solely based on interaction data can help us to predict the functions of proteins, and they are also useful to understand and explain certain biological processes.