In silico enzymatic synthesis of a 400,000 compound biochemical database for nontargeted metabolomics

Knowledge-based in silico fragmentation and annotation of mass spectra for natural products with MassKG

Liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) is a potent analytical technique utilized for identifying natural products from complex sources. However, due to the structural diversity, annotating LC-MS/MS data of natural products efficiently remains challenging, hindering the discovery process of novel active structures. Here, we introduce MassKG, an algorithm that combines a knowledge-based fragmentation strategy and a deep learning-based molecule generation model to aid in rapid dereplication and the discovery of novel NP structures. Specifically, MassKG has compiled 407,720 known NP structures and, based on this, generated 266,353 new structures using chemical language models for the discovery of potential novel compounds. Furthermore, MassKG demonstrates exceptional performance in spectra annotation compared to state-of-the-art algorithms. To enhance usability, MassKG has been implemented as a web server for annotating tandem mass spectral data (MS/MS, MS2) with a user-friendly interface, automatic reporting, and fragment tree visualization. Lastly, the interpretive capability of MassKG is comprehensively validated through composition analysis and MS annotation of Panax notoginseng, Ginkgo biloba, Codonopsis pilosula, and Astragalus membranaceus. MassKG is now accessible at https://xomics.com.cn/masskg.

Unlocking The Mysteries of DNA Adducts with Artificial Intelligence

Cellular genome is considered a dynamic blueprint of a cell since it encodes genetic information that gets temporally altered due to various endogenous and exogenous insults. Largely, the extent of genomic dynamicity is controlled by the trade-off between DNA repair processes and the genotoxic potential of the causative agent (genotoxins or potential carcinogens). A subset of genotoxins form DNA adducts by covalently binding to the cellular DNA, triggering structural or functional changes that lead to significant alterations in cellular processes via genetic (e. g., mutations) or non-genetic (e. g., epigenome) routes. Identification, quantification, and characterization of DNA adducts are indispensable for their comprehensive understanding and could expedite the ongoing efforts in predicting carcinogenicity and their mode of action. In this review, we elaborate on using Artificial Intelligence (AI)-based modeling in adducts biology and present multiple computational strategies to gain advancements in decoding DNA adducts. The proposed AI-based strategies encompass predictive modeling for adduct formation via metabolic activation, novel adducts' identification, prediction of biochemical routes for adduct formation, adducts' half-life predictions within biological ecosystems, and, establishing methods to predict the link between adducts chemistry and its location within the genomic DNA. In summary, we discuss some futuristic AI-based approaches in DNA adduct biology.

Combining Experimental with Computational Infrared and Mass Spectra for High-Throughput Nontargeted Chemical Structure Identification

The inability to identify the structures of most metabolites detected in environmental or biological samples limits the utility of nontargeted metabolomics. The most widely used analytical approaches combine mass spectrometry and machine learning methods to rank candidate structures contained in large chemical databases. Given the large chemical space typically searched, the use of additional orthogonal data may improve the identification rates and reliability. Here, we present results of combining experimental and computational mass and IR spectral data for high-throughput nontargeted chemical structure identification. Experimental MS/MS and gas-phase IR data for 148 test compounds were obtained from NIST. Candidate structures for each of the test compounds were obtained from PubChem (mean = 4444 candidate structures per test compound). Our workflow used CSI:FingerID to initially score and rank the candidate structures. The top 1000 ranked candidates were subsequently used for IR spectra prediction, scoring, and ranking using density functional theory (DFT-IR). Final ranking of the candidates was based on a composite score calculated as the average of the CSI:FingerID and DFT-IR rankings. This approach resulted in the correct identification of 88 of the 148 test compounds (59%). 129 of the 148 test compounds (87%) were ranked within the top 20 candidates. These identification rates are the highest yet reported when candidate structures are used from PubChem. Combining experimental and computational MS/MS and IR spectral data is a potentially powerful option for prioritizing candidates for final structure verification.

FT-ICR-MS-based metabolomics: A deep dive into plant metabolism

Metabolomics involves the identification and quantification of metabolites to unravel the chemical footprints behind cellular regulatory processes and to decipher metabolic networks, opening new insights to understand the correlation between genes and metabolites. In plants, it is estimated the existence of hundreds of thousands of metabolites and the majority is still unknown. Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS) is a powerful analytical technique to tackle such challenges. The resolving power and sensitivity of this ultrahigh mass accuracy mass analyzer is such that a complex mixture, such as plant extracts, can be analyzed and thousands of metabolite signals can be detected simultaneously and distinguished based on the naturally abundant elemental isotopes. In this review, FT-ICR-MS-based plant metabolomics studies are described, emphasizing FT-ICR-MS increasing applications in plant science through targeted and untargeted approaches, allowing for a better understanding of plant development, responses to biotic and abiotic stresses, and the discovery of new natural nutraceutical compounds. Improved metabolite extraction protocols compatible with FT-ICR-MS, metabolite analysis methods and metabolite identification platforms are also explored as well as new in silico approaches. Most recent advances in MS imaging are also discussed.

Metabolomics: small molecules that matter more

Metabolomics, an analytical study with high-throughput profiling, helps to understand interactions within a biological system. Small molecules, called metabolites or metabolomes with the size of <1500 Da, depict the status of a biological system in a different manner. Currently, we are in need to globally analyze the metabolome and the pathways involved in healthy, as well as diseased conditions, for possible therapeutic applications. Metabolome analysis has revealed high-abundance molecules during different conditions such as diet, environmental stress, microbiota, and disease and treatment states. As a result, it is hard to understand the complete and stable network of metabolites of a biological system. This review helps readers know the available techniques to study metabolomics in addition to other major omics such as genomics, transcriptomics, and proteomics. This review also discusses the metabolomics in various pathological conditions and the importance of metabolomics in therapeutic applications.

Seeing the forest for the trees: Retrieving plant secondary biochemical pathways from metabolome networks

Over the last decade, a giant leap forward has been made in resolving the main bottleneck in metabolomics, i.e., the structural characterization of the many unknowns. This has led to the next challenge in this research field: retrieving biochemical pathway information from the various types of networks that can be constructed from metabolome data. Searching putative biochemical pathways, referred to as biotransformation paths, is complicated because several flaws occur during the construction of metabolome networks. Multiple network analysis tools have been developed to deal with these flaws, while in silico retrosynthesis is appearing as an alternative approach. In this review, the different types of metabolome networks, their flaws, and the various tools to trace these biotransformation paths are discussed.

MAGI: A Method for Metabolite Annotation and Gene Integration

Metabolomics is a widely used technology for obtaining direct measures of metabolic activities from diverse biological systems. However, ambiguous metabolite identifications are a common challenge and biochemical interpretation is often limited by incomplete and inaccurate genome-based predictions of enzyme activities (that is, gene annotations). Metabolite Annotation and Gene Integration (MAGI) generates a metabolite-gene association score using a biochemical reaction network. This is calculated by a method that emphasizes consensus between metabolites and genes via biochemical reactions. To demonstrate the potential of this method, we applied MAGI to integrate sequence data and metabolomics data collected from Streptomyces coelicolor A3(2), an extensively characterized bacterium that produces diverse secondary metabolites. Our findings suggest that coupling metabolomics and genomics data by scoring consensus between the two increases the quality of both metabolite identifications and gene annotations in this organism. MAGI also made biochemical predictions for poorly annotated genes that were consistent with the extensive literature on this important organism. This limited analysis suggests that using metabolomics data has the potential to improve annotations in sequenced organisms and also provides testable hypotheses for specific biochemical functions. MAGI is freely available for academic use both as an online tool at https://magi.nersc.gov and with source code available at https://github.com/biorack/magi .

Advancing internal exposure and physiologically-based toxicokinetic modeling for 21st-century risk assessments

Scientifically sound, risk-informed evaluation of chemicals is essential to protecting public health. Systematically leveraging information from exposure, toxicology, and epidemiology studies can provide a holistic understanding of how real-world exposure to chemicals may impact the health of populations, including sensitive and vulnerable individuals and life-stages. Increasingly, public health policy makers are employing toxicokinetic (TK) modeling tools to integrate these data streams and predict potential human health impact. Development of a suite of tools for predicting internal exposure, including physiologically-based toxicokinetic (PBTK) models, is being driven by needs to address large numbers of data-poor chemicals efficiently, translate bioactivity, and mechanistic information from new in vitro test systems, and integrate multiple lines of evidence to enable scientifically sound, risk-informed decisions. New modeling approaches are being designed "fit for purpose" to inform specific decision contexts, with applications ranging from rapid screening of hundreds of chemicals, to improved prediction of risks during sensitive stages of development. New data are being generated experimentally and computationally to support these models. Progress to meet the demand for internal exposure and PBTK modeling tools will require transparent publication of models and data to build credibility in results, as well as opportunities to partner with decision makers to evaluate and build confidence in use of these for improved decisions that promote safe use of chemicals.

Integrating tools for non-targeted analysis research and chemical safety evaluations at the US EPA

Tens-of-thousands of chemicals are registered in the U.S. for use in countless processes and products. Recent evidence suggests that many of these chemicals are measureable in environmental and/or biological systems, indicating the potential for widespread exposures. Traditional public health research tools, including in vivo studies and targeted analytical chemistry methods, have been unable to meet the needs of screening programs designed to evaluate chemical safety. As such, new tools have been developed to enable rapid assessment of potentially harmful chemical exposures and their attendant biological responses. One group of tools, known as "non-targeted analysis" (NTA) methods, allows the rapid characterization of thousands of never-before-studied compounds in a wide variety of environmental, residential, and biological media. This article discusses current applications of NTA methods, challenges to their effective use in chemical screening studies, and ways in which shared resources (e.g., chemical standards, databases, model predictions, and media measurements) can advance their use in risk-based chemical prioritization. A brief review is provided of resources and projects within EPA's Office of Research and Development (ORD) that provide benefit to, and receive benefits from, NTA research endeavors. A summary of EPA's Non-Targeted Analysis Collaborative Trial (ENTACT) is also given, which makes direct use of ORD resources to benefit the global NTA research community. Finally, a research framework is described that shows how NTA methods will bridge chemical prioritization efforts within ORD. This framework exists as a guide for institutions seeking to understand the complexity of chemical exposures, and the impact of these exposures on living systems.

Identification of small molecules using accurate mass MS/MS search

Tandem mass spectral library search (MS/MS) is the fastest way to correctly annotate MS/MS spectra from screening small molecules in fields such as environmental analysis, drug screening, lipid analysis, and metabolomics. The confidence in MS/MS-based annotation of chemical structures is impacted by instrumental settings and requirements, data acquisition modes including data-dependent and data-independent methods, library scoring algorithms, as well as post-curation steps. We critically discuss parameters that influence search results, such as mass accuracy, precursor ion isolation width, intensity thresholds, centroiding algorithms, and acquisition speed. A range of publicly and commercially available MS/MS databases such as NIST, MassBank, MoNA, LipidBlast, Wiley MSforID, and METLIN are surveyed. In addition, software tools including NIST MS Search, MS-DIAL, Mass Frontier, SmileMS, Mass++, and XCMS2 to perform fast MS/MS search are discussed. MS/MS scoring algorithms and challenges during compound annotation are reviewed. Advanced methods such as the in silico generation of tandem mass spectra using quantum chemistry and machine learning methods are covered. Community efforts for curation and sharing of tandem mass spectra that will allow for faster distribution of scientific discoveries are discussed.

Cardiovascular Metabolomics

Disturbances in cardiac metabolism underlie most cardiovascular diseases. Metabolomics, one of the newer omics technologies, has emerged as a powerful tool for defining changes in both global and cardiac-specific metabolism that occur across a spectrum of cardiovascular disease states. Findings from metabolomics studies have contributed to better understanding of the metabolic changes that occur in heart failure and ischemic heart disease and have identified new cardiovascular disease biomarkers. As technologies advance, the metabolomics field continues to evolve rapidly. In this review, we will discuss the current state of metabolomics technologies, including consideration of various metabolomics platforms and elements of study design; the emerging utility of stable isotopes for metabolic flux studies; and the use of metabolomics to better understand specific cardiovascular diseases, with an emphasis on recent advances in the field.

Designed inhibitors with hetero linkers for gastric proton pump H+,K+-ATPase: Steered molecular dynamics and metadynamics studies

Acid suppressant SCH28080 and its derivatives reversibly reduce acid secretion activity of the H⁺,K⁺-ATPase in a K⁺ competitive manner. The results on homologation of the SCH28080 by varying the linker chain length suggested the improvement in efficacy. However, the pharmacokinetic studies reveal that the hydrophobic nature of the CH₂ linker units may not help it to function as a better acid suppressant. We have exploited the role of linker unit to enhance the efficacy of such reversible acid suppressant drug molecules using hetero linker, i.e., disulfide and peroxy linkers. The logarithm of partition coefficient defined for a drug molecule relates to the partition coefficient, which allows the optimum solubility characteristics to reach the active site. The logarithm of partition coefficient calculated for the designed inhibitors suggests that inhibitors would possibly reach the active site in sufficient concentration like in the case of SCH28080. The steered molecular dynamics studies have revealed that the Inhibitor-1 with disulfide linker unit is more stable at the active site due to greater noncovalent interactions compared to the SCH28080. Centre of mass distance analysis suggests that the Cysteine-813 amino acid residue selectively plays an important role in the inhibition of H⁺,K⁺-ATPase for Inhibitor-1. Furthermore, the quantum chemical calculations with M11L/6-31+G(d,p) level of theory have been performed to account the noncovalent interactions responsible for the stabilization of inhibitor molecules in the active site gorge of the gastric proton pump at different time scale. The hydrogen bonding and hydrophobic interaction studies corroborate the center of mass distance analysis as well. Well-tempered metadynamics free energy surface and center of mass separation analysis for the Inhibitor-1 is in good agreement with the steered molecular dynamics results. The torsional angle of the linker units seems to be crucial for better efficacy of drug molecules. The torsional angle of linker units of SCH28080 (COCH₂C) and of Inhibitor 1 (CSSC) prefers to lie within ∼60°-90° for a longer time during the simulations, whereas, the peroxy linker (COOC) of Inhibitor 2 prefers to adopt ∼120-160°. Therefore, it appears that the smaller torsion angle of linker units can achieve better interactions with the active site residues of H⁺,K⁺-ATPase to inhibit the acid secretion activity. The reversible drug molecules with disulfide linker unit would be a promising candidate as proton pump antagonist to H⁺,K⁺-ATPase.

Mining molecular structure databases: Identification of small molecules based on fragmentation mass spectrometry data

Mass spectrometry (MS) is a key technology for the analysis of small molecules. For the identification and structural elucidation of novel molecules, new approaches beyond straightforward spectral comparison are required. In this review, we will cover computational methods that help with the identification of small molecules by analyzing fragmentation MS data. We focus on the four main approaches to mine a database of metabolite structures, that is rule-based fragmentation spectrum prediction, combinatorial fragmentation, competitive fragmentation modeling, and molecular fingerprint prediction. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 36:624-633, 2017.

Quantifying confidence in the reporting of metabolic biotransformations

How confident can we be in the assignment of metabolite structures? Are the analytical techniques used sufficient to support hypotheses about what is being formed? In this Feature, we discuss the results of an extensive survey into the analytical techniques used, and their value in the characterisation of metabolites. The survey covers the structures of over 16000 metabolites formed from 1732 query compounds, covering over 35 years of the literature and a variety of journals. The value of different characterisation techniques is considered, alongside or in the absence of synthetic standards. The changes in analytical techniques used over time are briefly considered, and a metric for the confidence that a claimed metabolite has been confirmed is proposed.

Integrating ion mobility spectrometry into mass spectrometry-based exposome measurements: what can it add and how far can it go?

Measuring the exposome remains a challenge due to the range and number of anthropogenic molecules that are encountered in our daily lives, as well as the complex systemic responses to these exposures. One option for improving the coverage, dynamic range and throughput of measurements is to incorporate ion mobility spectrometry (IMS) into current MS-based analytical methods. The implementation of IMS in exposomics studies will lead to more frequent observations of previously undetected chemicals and metabolites. LC-IMS-MS will provide increased overall measurement dynamic range, resulting in detections of lower abundance molecules. Alternatively, the throughput of IMS-MS alone will provide the opportunity to analyze many thousands of longitudinal samples over lifetimes of exposure, capturing evidence of transitory accumulations of chemicals or metabolites. The volume of data corresponding to these new chemical observations will almost certainly outpace the generation of reference data to enable their confident identification. In this perspective, we briefly review the state-of-the-art in measuring the exposome, and discuss the potential use for IMS-MS and the physico-chemical property of collisional cross section in both exposure assessment and molecular identification.

Searching molecular structure databases with tandem mass spectra using CSI:FingerID

Metabolites provide a direct functional signature of cellular state. Untargeted metabolomics experiments usually rely on tandem MS to identify the thousands of compounds in a biological sample. Today, the vast majority of metabolites remain unknown. We present a method for searching molecular structure databases using tandem MS data of small molecules. Our method computes a fragmentation tree that best explains the fragmentation spectrum of an unknown molecule. We use the fragmentation tree to predict the molecular structure fingerprint of the unknown compound using machine learning. This fingerprint is then used to search a molecular structure database such as PubChem. Our method is shown to improve on the competing methods for computational metabolite identification by a considerable margin.

MINEs: open access databases of computationally predicted enzyme promiscuity products for untargeted metabolomics

BACKGROUND

In spite of its great promise, metabolomics has proven difficult to execute in an untargeted and generalizable manner. Liquid chromatography-mass spectrometry (LC-MS) has made it possible to gather data on thousands of cellular metabolites. However, matching metabolites to their spectral features continues to be a bottleneck, meaning that much of the collected information remains uninterpreted and that new metabolites are seldom discovered in untargeted studies. These challenges require new approaches that consider compounds beyond those available in curated biochemistry databases.

DESCRIPTION

Here we present Metabolic In silico Network Expansions (MINEs), an extension of known metabolite databases to include molecules that have not been observed, but are likely to occur based on known metabolites and common biochemical reactions. We utilize an algorithm called the Biochemical Network Integrated Computational Explorer (BNICE) and expert-curated reaction rules based on the Enzyme Commission classification system to propose the novel chemical structures and reactions that comprise MINE databases. Starting from the Kyoto Encyclopedia of Genes and Genomes (KEGG) COMPOUND database, the MINE contains over 571,000 compounds, of which 93% are not present in the PubChem database. However, these MINE compounds have on average higher structural similarity to natural products than compounds from KEGG or PubChem. MINE databases were able to propose annotations for 98.6% of a set of 667 MassBank spectra, 14% more than KEGG alone and equivalent to PubChem while returning far fewer candidates per spectra than PubChem (46 vs. 1715 median candidates). Application of MINEs to LC-MS accurate mass data enabled the identity of an unknown peak to be confidently predicted.

CONCLUSIONS

MINE databases are freely accessible for non-commercial use via user-friendly web-tools at http://minedatabase.mcs.anl.gov and developer-friendly APIs. MINEs improve metabolomics peak identification as compared to general chemical databases whose results include irrelevant synthetic compounds. Furthermore, MINEs complement and expand on previous in silico generated compound databases that focus on human metabolism. We are actively developing the database; future versions of this resource will incorporate transformation rules for spontaneous chemical reactions and more advanced filtering and prioritization of candidate structures. Graphical abstractMINE database construction and access methods. The process of constructing a MINE database from the curated source databases is depicted on the left. The methods for accessing the database are shown on the right.

Metabolic pathway predictions for metabolomics: a molecular structure matching approach

Metabolic pathways are composed of a series of chemical reactions occurring within a cell. In each pathway, enzymes catalyze the conversion of substrates into structurally similar products. Thus, structural similarity provides a potential means for mapping newly identified biochemical compounds to known metabolic pathways. In this paper, we present TrackSM, a cheminformatics tool designed to associate a chemical compound to a known metabolic pathway based on molecular structure matching techniques. Validation experiments show that TrackSM is capable of associating 93% of tested structures to their correct KEGG pathway class and 88% to their correct individual KEGG pathway. This suggests that TrackSM may be a valuable tool to aid in associating previously unknown small molecules to known biochemical pathways and improve our ability to link metabolomics, proteomic, and genomic data sets. TrackSM is freely available at http://metabolomics.pharm.uconn.edu/?q=Software.html .

Analysis of phenolic and cyclic compounds in plants using derivatization techniques in combination with GC-MS-based metabolite profiling

Metabolite profiling has been established as a modern technology platform for the description of complex chemical matrices and compound identification in biological samples. Gas chromatography coupled with mass spectrometry (GC-MS) in particular is a fast and accurate method widely applied in diagnostics, functional genomics and for screening purposes. Following solvent extraction and derivatization, hundreds of metabolites from different chemical groups can be characterized in one analytical run. Besides sugars, acids, and polyols, diverse phenolic and other cyclic metabolites can be efficiently detected by metabolite profiling. The review describes own results from plant research to exemplify the applicability of GC-MS profiling and concurrent detection and identification of phenolics and other cyclic structures.

Efficient searching and annotation of metabolic networks using chemical similarity

MOTIVATION

The urgent need for efficient and sustainable biological production of fuels and high-value chemicals has elicited a wave of in silico techniques for identifying promising novel pathways to these compounds in large putative metabolic networks. To date, these approaches have primarily used general graph search algorithms, which are prohibitively slow as putative metabolic networks may exceed 1 million compounds. To alleviate this limitation, we report two methods--SimIndex (SI) and SimZyme--which use chemical similarity of 2D chemical fingerprints to efficiently navigate large metabolic networks and propose enzymatic connections between the constituent nodes. We also report a Byers-Waterman type pathway search algorithm for further paring down pertinent networks.

RESULTS

Benchmarking tests run with SI show it can reduce the number of nodes visited in searching a putative network by 100-fold with a computational time improvement of up to 10(5)-fold. Subsequent Byers-Waterman search application further reduces the number of nodes searched by up to 100-fold, while SimZyme demonstrates ∼ 90% accuracy in matching query substrates with enzymes. Using these modules, we have designed and annotated an alternative to the methylerythritol phosphate pathway to produce isopentenyl pyrophosphate with more favorable thermodynamics than the native pathway. These algorithms will have a significant impact on our ability to use large metabolic networks that lack annotation of promiscuous reactions.

AVAILABILITY AND IMPLEMENTATION

Python files will be available for download at http://tyolab.northwestern.edu/tools/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The food metabolome: a window over dietary exposure

The food metabolome is defined as the part of the human metabolome directly derived from the digestion and biotransformation of foods and their constituents. With >25,000 compounds known in various foods, the food metabolome is extremely complex, with a composition varying widely according to the diet. By its very nature it represents a considerable and still largely unexploited source of novel dietary biomarkers that could be used to measure dietary exposures with a high level of detail and precision. Most dietary biomarkers currently have been identified on the basis of our knowledge of food compositions by using hypothesis-driven approaches. However, the rapid development of metabolomics resulting from the development of highly sensitive modern analytic instruments, the availability of metabolite databases, and progress in (bio)informatics has made agnostic approaches more attractive as shown by the recent identification of novel biomarkers of intakes for fruit, vegetables, beverages, meats, or complex diets. Moreover, examples also show how the scrutiny of the food metabolome can lead to the discovery of bioactive molecules and dietary factors associated with diseases. However, researchers still face hurdles, which slow progress and need to be resolved to bring this emerging field of research to maturity. These limits were discussed during the First International Workshop on the Food Metabolome held in Glasgow. Key recommendations made during the workshop included more coordination of efforts; development of new databases, software tools, and chemical libraries for the food metabolome; and shared repositories of metabolomic data. Once achieved, major progress can be expected toward a better understanding of the complex interactions between diet and human health.

In silico prediction and automatic LC-MS(n) annotation of green tea metabolites in urine

The colonic breakdown and human biotransformation of small molecules present in food can give rise to a large variety of potentially bioactive metabolites in the human body. However, the absence of reference data for many of these components limits their identification in complex biological samples, such as plasma and urine. We present an in silico workflow for automatic chemical annotation of metabolite profiling data from liquid chromatography coupled with multistage accurate mass spectrometry (LC-MS(n)), which we used to systematically screen for the presence of tea-derived metabolites in human urine samples after green tea consumption. Reaction rules for intestinal degradation and human biotransformation were systematically applied to chemical structures of 75 green tea components, resulting in a virtual library of 27,245 potential metabolites. All matching precursor ions in the urine LC-MS(n) data sets, as well as the corresponding fragment ions, were automatically annotated by in silico generated (sub)structures. The results were evaluated based on 74 previously identified urinary metabolites and lead to the putative identification of 26 additional green tea-derived metabolites. A total of 77% of all annotated metabolites were not present in the Pubchem database, demonstrating the benefit of in silico metabolite prediction for the automatic annotation of yet unknown metabolites in LC-MS(n) data from nutritional metabolite profiling experiments.