Good practices and recommendations for using and benchmarking computational metabolomics metabolite annotation tools

Highly reliable LC-MS lipidomics database for efficient human plasma profiling based on NIST SRM 1950

Liquid chromatography coupled to high resolution mass spectrometry (LC-HRMS)-based methods have become the gold standard methodology for the comprehensive profiling of the human plasma lipidome. However, both the complexity of lipid chemistry and LC-HRMS-associated data pose challenges to the characterization of this biological matrix. In accordance with the current consensus of quality requirements for LC-HRMS lipidomics data, we aimed to characterize the NIST® Standard Reference Material for Human Plasma (SRM 1950) using an LC-ESI(+/-)-MS method compatible with high-throughput lipidome profiling. We generated a highly curated lipid database with increased coverage, quality, and consistency, including additional quality assurance procedures involving adduct formation, within-method m/z evaluation, retention behavior of species within lipid chain isomers, and expert-driven resolution of isomeric and isobaric interferences. As a proof-of-concept, we showed the utility of our in-house LC-MS lipidomic database -consisting of 592 lipid entries- for the fast, comprehensive, and reliable lipidomic profiling of the human plasma from healthy human volunteers. We are confident that the implementation of this robust resource and methodology will have a significant impact by reducing data redundancy and the current delays and bottlenecks in untargeted plasma lipidomic studies.

Combined LC-MS/MS feature grouping, statistical prioritization, and interactive networking in msFeaST

SUMMARY

Computational metabolomics workflows have revolutionized the untargeted metabolomics field. However, the organization and prioritization of metabolite features remains a laborious process. Organizing metabolomics data is often done through mass fragmentation-based spectral similarity grouping, resulting in feature sets that also represent an intuitive and scientifically meaningful first stage of analysis in untargeted metabolomics. Exploiting such feature sets, feature-set testing has emerged as an approach that is widely used in genomics and targeted metabolomics pathway enrichment analyses. It allows for formally combining groupings with statistical testing into more meaningful pathway enrichment conclusions. Here, we present msFeaST (mass spectral Feature Set Testing), a feature-set testing and visualization workflow for LC-MS/MS untargeted metabolomics data. Feature-set testing involves statistically assessing differential abundance patterns for groups of features across experimental conditions. We developed msFeaST to make use of spectral similarity-based feature groupings generated using k-medoids clustering, where the resulting clusters serve as a proxy for grouping structurally similar features with potential biosynthesis pathway relationships. Spectral clustering done in this way allows for feature group-wise statistical testing using the globaltest package, which provides high power to detect small concordant effects via joint modeling and reduced multiplicity adjustment penalties. Hence, msFeaST provides interactive integration of the semi-quantitative experimental information with mass-spectral structural similarity information, enhancing the prioritization of features and feature sets during exploratory data analysis.

AVAILABILITY AND IMPLEMENTATION

The msFeaST workflow is freely available through https://github.com/kevinmildau/msFeaST and built to work on MacOS and Linux systems.

Statistical analysis of feature-based molecular networking results from non-targeted metabolomics data

Feature-based molecular networking (FBMN) is a popular analysis approach for liquid chromatography-tandem mass spectrometry-based non-targeted metabolomics data. While processing liquid chromatography-tandem mass spectrometry data through FBMN is fairly streamlined, downstream data handling and statistical interrogation are often a key bottleneck. Especially users new to statistical analysis struggle to effectively handle and analyze complex data matrices. Here we provide a comprehensive guide for the statistical analysis of FBMN results, focusing on the downstream analysis of the FBMN output table. We explain the data structure and principles of data cleanup and normalization, as well as uni- and multivariate statistical analysis of FBMN results. We provide explanations and code in two scripting languages (R and Python) as well as the QIIME2 framework for all protocol steps, from data clean-up to statistical analysis. All code is shared in the form of Jupyter Notebooks ( https://github.com/Functional-Metabolomics-Lab/FBMN-STATS ). Additionally, the protocol is accompanied by a web application with a graphical user interface ( https://fbmn-statsguide.gnps2.org/ ) to lower the barrier of entry for new users and for educational purposes. Finally, we also show users how to integrate their statistical results into the molecular network using the Cytoscape visualization tool. Throughout the protocol, we use a previously published environmental metabolomics dataset for demonstration purposes. Together, the protocol, code and web application provide a complete guide and toolbox for FBMN data integration, cleanup and advanced statistical analysis, enabling new users to uncover molecular insights from their non-targeted metabolomics data. Our protocol is tailored for the seamless analysis of FBMN results from Global Natural Products Social Molecular Networking and can be easily adapted to other mass spectrometry feature detection, annotation and networking tools.

Critical review on in silico methods for structural annotation of chemicals detected with LC/HRMS non-targeted screening

Non-targeted screening with liquid chromatography coupled to high-resolution mass spectrometry (LC/HRMS) is increasingly leveraging in silico methods, including machine learning, to obtain candidate structures for structural annotation of LC/HRMS features and their further prioritization. Candidate structures are commonly retrieved based on the tandem mass spectral information either from spectral or structural databases; however, the vast majority of the detected LC/HRMS features remain unannotated, constituting what we refer to as a part of the unknown chemical space. Recently, the exploration of this chemical space has become accessible through generative models. Furthermore, the evaluation of the candidate structures benefits from the complementary empirical analytical information such as retention time, collision cross section values, and ionization type. In this critical review, we provide an overview of the current approaches for retrieving and prioritizing candidate structures. These approaches come with their own set of advantages and limitations, as we showcase in the example of structural annotation of ten known and ten unknown LC/HRMS features. We emphasize that these limitations stem from both experimental and computational considerations. Finally, we highlight three key considerations for the future development of in silico methods.

An Ensemble Spectral Prediction (ESP) model for metabolite annotation

MOTIVATION

A key challenge in metabolomics is annotating measured spectra from a biological sample with chemical identities. Currently, only a small fraction of measurements can be assigned identities. Two complementary computational approaches have emerged to address the annotation problem: mapping candidate molecules to spectra, and mapping query spectra to molecular candidates. In essence, the candidate molecule with the spectrum that best explains the query spectrum is recommended as the target molecule. Despite candidate ranking being fundamental in both approaches, limited prior works incorporated rank learning tasks in determining the target molecule.

RESULTS

We propose a novel machine learning model, Ensemble Spectral Prediction (ESP), for metabolite annotation. ESP takes advantage of prior neural network-based annotation models that utilize multilayer perceptron (MLP) networks and Graph Neural Networks (GNNs). Based on the ranking results of the MLP- and GNN-based models, ESP learns a weighting for the outputs of MLP and GNN spectral predictors to generate a spectral prediction for a query molecule. Importantly, training data is stratified by molecular formula to provide candidate sets during model training. Further, baseline MLP and GNN models are enhanced by considering peak dependencies through label mixing and multi-tasking on spectral topic distributions. When trained on the NIST 2020 dataset and evaluated on the relevant candidate sets from PubChem, ESP improves average rank by 23.7% and 37.2% over the MLP and GNN baselines, respectively, demonstrating performance gain over state-of-the-art neural network approaches. However, MLP approaches remain strong contenders when considering top five ranks. Importantly, we show that annotation performance is dependent on the training dataset, the number of molecules in the candidate set and candidate similarity to the target molecule.

AVAILABILITY AND IMPLEMENTATION

The ESP code, a trained model, and a Jupyter notebook that guide users on using the ESP tool is available at https://github.com/HassounLab/ESP.

NMR Fingerprinting of Conventional and Genetically Modified Soybean Plants with AtAREB1 Transcription Factors

Drought stress impacts soybean yields and physiological processes. However, the insertion of the activated form of the AtAREB1 gene in the soybean cultivar BR16, which is sensitive to water deficit, improved the drought response of the genetically modified plants. Thus, in this study, we used 1H NMR in solution and solid-state NMR to investigate the response of genetically modified soybean overexpressing AtAREB1 under water deficiency conditions. We achieved that drought-tolerant soybean yields high content of amino acids isoleucine, leucine, threonine, valine, proline, glutamate, aspartate, asparagine, tyrosine, and phenylalanine after 12 days of drought stress conditions, as compared to drought-sensitive soybean under the same conditions. Specific target compounds, including sugars, organic acids, and phenolic compounds, were identified as involved in controlling sensitive soybean during the vegetative stage. Solid-state NMR was used to study the impact of drought stress on starch and cellulose contents in different soybean genotypes. The findings provide insights into the metabolic adjustments of soybean overexpressing AREB transcription factors in adapting to dry climates. This study presents NMR techniques for investigating the metabolome of transgenic soybean plants in response to the water deficit. The approach allowed for the identification of physiological and morphological changes in drought-resistant and drought-tolerant soybean tissues. The findings indicate that drought stress significantly alters micro- and macromolecular metabolism in soybean plants. Differential responses were observed among roots and leaves as well as drought-tolerant and drought-sensitive cultivars, highlighting the complex interplay between overexpressed transcription factors and drought stress in soybean plants.

Reproducible MS/MS library cleaning pipeline in matchms

Mass spectral libraries have proven to be essential for mass spectrum annotation, both for library matching and training new machine learning algorithms. A key step in training machine learning models is the availability of high-quality training data. Public libraries of mass spectrometry data that are open to user submission often suffer from limited metadata curation and harmonization. The resulting variability in data quality makes training of machine learning models challenging. Here we present a library cleaning pipeline designed for cleaning tandem mass spectrometry library data. The pipeline is designed with ease of use, flexibility, and reproducibility as leading principles.Scientific contributionThis pipeline will result in cleaner public mass spectral libraries that will improve library searching and the quality of machine-learning training datasets in mass spectrometry. This pipeline builds on previous work by adding new functionality for curating and correcting annotated libraries, by validating structure annotations. Due to the high quality of our software, the reproducibility, and improved logging, we think our new pipeline has the potential to become the standard in the field for cleaning tandem mass spectrometry libraries.

Common data models to streamline metabolomics processing and annotation, and implementation in a Python pipeline

To standardize metabolomics data analysis and facilitate future computational developments, it is essential to have a set of well-defined templates for common data structures. Here we describe a collection of data structures involved in metabolomics data processing and illustrate how they are utilized in a full-featured Python-centric pipeline. We demonstrate the performance of the pipeline, and the details in annotation and quality control using large-scale LC-MS metabolomics and lipidomics data and LC-MS/MS data. Multiple previously published datasets are also reanalyzed to showcase its utility in biological data analysis. This pipeline allows users to streamline data processing, quality control, annotation, and standardization in an efficient and transparent manner. This work fills a major gap in the Python ecosystem for computational metabolomics.

Comprehensive mass spectrometric metabolomic profiling of a chemically diverse collection of plants of the Celastraceae family

Natural products exhibit interesting structural features and significant biological activities. The discovery of new bioactive molecules is a complex process that requires high-quality metabolite profiling data to properly target the isolation of compounds of interest and enable their complete structural characterization. The same metabolite profiling data can also be used to better understand chemotaxonomic links between species. This Data Descriptor details a dataset resulting from the untargeted liquid chromatography-mass spectrometry metabolite profiling of 76 natural extracts of the Celastraceae family. The spectral annotation results and related chemical and taxonomic metadata are shared, along with proposed examples of data reuse. This data can be further studied by researchers exploring the chemical diversity of natural products. This can serve as a reference sample set for deep metabolome investigation of this chemically rich plant family.

Tailored Mass Spectral Data Exploration Using the SpecXplore Interactive Dashboard

Untargeted metabolomics promises comprehensive characterization of small molecules in biological samples. However, the field is hampered by low annotation rates and abstract spectral data. Despite recent advances in computational metabolomics, manual annotations and manual confirmation of in-silico annotations remain important in the field. Here, exploratory data analysis methods for mass spectral data provide overviews, prioritization, and structural hypothesis starting points to researchers facing large quantities of spectral data. In this research, we propose a fluid means of dealing with mass spectral data using specXplore, an interactive Python dashboard providing interactive and complementary visualizations facilitating mass spectral similarity matrix exploration. Specifically, specXplore provides a two-dimensional t-distributed stochastic neighbor embedding embedding as a jumping board for local connectivity exploration using complementary interactive visualizations in the form of partial network drawings, similarity heatmaps, and fragmentation overview maps. SpecXplore makes use of state-of-the-art ms2deepscore pairwise spectral similarities as a quantitative backbone while allowing fast changes of threshold and connectivity limitation settings, providing flexibility in adjusting settings to suit the localized node environment being explored. We believe that specXplore can become an integral part of mass spectral data exploration efforts and assist users in the generation of structural hypotheses for compounds of interest.

Hydrophilic Interaction Liquid Chromatography-Hydrogen/Deuterium Exchange-Mass Spectrometry (HILIC-HDX-MS) for Untargeted Metabolomics

Liquid chromatography with mass spectrometry (LC-MS)-based metabolomics detects thousands of molecular features (retention time-m/z pairs) in biological samples per analysis, yet the metabolite annotation rate remains low, with 90% of signals classified as unknowns. To enhance the metabolite annotation rates, researchers employ tandem mass spectral libraries and challenging in silico fragmentation software. Hydrogen/deuterium exchange mass spectrometry (HDX-MS) may offer an additional layer of structural information in untargeted metabolomics, especially for identifying specific unidentified metabolites that are revealed to be statistically significant. Here, we investigate the potential of hydrophilic interaction liquid chromatography (HILIC)-HDX-MS in untargeted metabolomics. Specifically, we evaluate the effectiveness of two approaches using hypothetical targets: the post-column addition of deuterium oxide (D2O) and the on-column HILIC-HDX-MS method. To illustrate the practical application of HILIC-HDX-MS, we apply this methodology using the in silico fragmentation software MS-FINDER to an unknown compound detected in various biological samples, including plasma, serum, tissues, and feces during HILIC-MS profiling, subsequently identified as N1-acetylspermidine.

An assessment of AcquireX and Compound Discoverer software 3.3 for non-targeted metabolomics

We used the Exploris 240 mass spectrometer for non-targeted metabolomics on Saccharomyces cerevisiae strain BY4741 and tested AcquireX software for increasing the number of detectable compounds and Compound Discoverer 3.3 software for identifying compounds by MS2 spectral library matching. AcquireX increased the number of potentially identifiable compounds by 50% through six iterations of MS2 acquisition. On the basis of high-scoring MS2 matches made by Compound Discoverer, there were 483 compounds putatively identified from nearly 8000 candidate spectra. Comparisons to 20 amino acid standards, however, revealed instances whereby compound matches could be incorrect despite strong scores. Situations included the candidate with the top score not being the correct compound, matching the same compound at two different chromatographic peaks, assigning the highest score to a library compound much heavier than the mass for the parent ion, and grouping MS2 isomers to a single parent ion. Because the software does not calculate false positive and false discovery rates at these multiple levels where such errors can propagate, we conclude that manual examination of findings will be required post software analysis. These results will interest scientists who may use this platform for metabolomics research in diverse disciplines including medical science, environmental science, and agriculture.

Discovery of Glycation Products: Unraveling the Unknown Glycation Space Using a Mass Spectral Library from In Vitro Model Systems

The nonenzymatic reaction between amino acids (AAs) and reducing sugars, also known as the Maillard reaction, is the primary source of free glycation products (GPs) in vivo and in vitro. The limited number of MS/MS records for GPs in public libraries hinders the annotation and investigation of nonenzymatic glycation. To address this issue, we present a mass spectral library containing the experimental MS/MS spectra of diverse GPs from model systems. Based on the conceptional reaction processes and structural characteristics of products, we classified GPs into common GPs (CGPs) and modified AAs (MAAs). A workflow for annotating GPs was established based on the structural and fragmentation patterns of each GP type. The final spectral library contains 157 CGPs, 499 MAAs, and 2426 GP spectra with synthetic model system information, retention time, precursor m/z, MS/MS, and annotations. As a proof-of-concept, we demonstrated the use of the library for screening GPs in unidentified spectra of human plasma and urine. The AAs with the C6H10O5 modification, fructosylation from Amadori rearrangement, were the most found GPs. With the help of the model system, we confirmed the existence of C6H10O5-modified Valine in human plasma by matching both retention time, MS1, and MS/MS without reference standards. In summary, our GP library can serve as an online resource to quickly screen possible GPs in an untargeted metabolomics workflow, furthermore with the model system as a practical synthesis method to confirm their identity.

Common data models to streamline metabolomics processing and annotation, and implementation in a Python pipeline

To standardize metabolomics data analysis and facilitate future computational developments, it is essential is have a set of well-defined templates for common data structures. Here we describe a collection of data structures involved in metabolomics data processing and illustrate how they are utilized in a full-featured Python-centric pipeline. We demonstrate the performance of the pipeline, and the details in annotation and quality control using large-scale LC-MS metabolomics and lipidomics data and LC-MS/MS data. Multiple previously published datasets are also reanalyzed to showcase its utility in biological data analysis. This pipeline allows users to streamline data processing, quality control, annotation, and standardization in an efficient and transparent manner. This work fills a major gap in the Python ecosystem for computational metabolomics.

Mapping the metabolic responses to oxaliplatin-based chemotherapy with in vivo spatiotemporal metabolomics

Adjuvant chemotherapy improves the survival outlook for patients undergoing operations for lung metastases caused by colorectal cancer (CRC). However, a multidisciplinary approach that evaluates several factors related to patient and tumor characteristics is necessary for managing chemotherapy treatment in metastatic CRC patients with lung disease, as such factors dictate the timing and drug regimen, which may affect treatment response and prognosis. In this study, we explore the potential of spatial metabolomics for evaluating metabolic phenotypes and therapy outcomes during the local delivery of the anticancer drug, oxaliplatin, to the lung. 12 male Yorkshire pigs underwent a 3 h left lung in vivo lung perfusion (IVLP) with various doses of oxaliplatin (7.5, 10, 20, 40, and 80 mg/L), which were administered to the perfusion circuit reservoir as a bolus. Biocompatible solid-phase microextraction (SPME) microprobes were combined with global metabolite profiling to obtain spatiotemporal information about the activity of the drug, determine toxic doses that exceed therapeutic efficacy, and conduct a mechanistic exploration of associated lung injury. Mild and subclinical lung injury was observed at 40 mg/L of oxaliplatin, and significant compromise of the hemodynamic lung function was found at 80 mg/L. This result was associated with massive alterations in metabolic patterns of lung tissue and perfusate, resulting in a total of 139 discriminant compounds. Uncontrolled inflammatory response, abnormalities in energy metabolism, and mitochondrial dysfunction next to accelerated kynurenine and aldosterone production were recognized as distinct features of dysregulated metabolipidome. Spatial pharmacometabolomics may be a promising tool for identifying pathological responses to chemotherapy.

Machine learning-assisted structure annotation of natural products based on MS and NMR data

Covering: up to March 2023Machine learning (ML) has emerged as a popular tool for analyzing the structures of natural products (NPs). This review presents a summary of the recent advancements in ML-assisted mass spectrometry (MS) and nuclear magnetic resonance (NMR) data analysis to establish the chemical structures of NPs. First, ML-based MS/MS analyses that rely on library matching are discussed, which involves the utilization of ML algorithms to calculate similarity, predict the MS/MS fragments, and form molecular fingerprint. Then, ML assisted MS/MS structural annotation without library matching is reviewed. Furthermore, the cases of ML algorithms in assisting structural studies of NPs based on NMR are discussed from four perspectives: NMR prediction, functional group identification, structural categorization and quantum chemical calculation. Finally, the review concludes with a discussion of the challenges and the trends associated with the structural establishment of NPs based on ML algorithms.

Persistent Biofluid Small-Molecule Alterations Induced by Trypanosoma cruzi Infection Are Not Restored by Parasite Elimination

Chagas disease (CD), caused by Trypanosoma cruzi (T. cruzi) protozoa, is a complicated parasitic illness with inadequate medical measures for diagnosing infection and monitoring treatment success. To address this gap, we analyzed changes in the metabolome of T. cruzi-infected mice via liquid chromatography tandem mass spectrometry of clinically accessible biofluids: saliva, urine, and plasma. Urine was the most indicative of infection status across mouse and parasite genotypes. Metabolites perturbed by infection in urine include kynurenate, acylcarnitines, and threonylcarbamoyladenosine. Based on these results, we sought to implement urine as a tool for the assessment of CD treatment success. Strikingly, it was found that mice with parasite clearance following benznidazole antiparasitic treatment had an overall urine metabolome comparable to that of mice that failed to clear parasites. These results provide a complementary hypothesis to explain clinical trial data in which benznidazole treatment did not improve patient outcomes in late-stage disease, even in patients with successful parasite clearance. Overall, this study provides insights into new small-molecule-based CD diagnostic methods and a new approach to assess functional responses to treatment.

Localized cardiac small molecule trajectories and persistent chemical sequelae in experimental Chagas disease

Post-infectious conditions present major health burdens but remain poorly understood. In Chagas disease (CD), caused by Trypanosoma cruzi parasites, antiparasitic agents that successfully clear T. cruzi do not always improve clinical outcomes. In this study, we reveal differential small molecule trajectories between cardiac regions during chronic T. cruzi infection, matching with characteristic CD apical aneurysm sites. Incomplete, region-specific, cardiac small molecule restoration is observed in animals treated with the antiparasitic benznidazole. In contrast, superior restoration of the cardiac small molecule profile is observed for a combination treatment of reduced-dose benznidazole plus an immunotherapy, even with less parasite burden reduction. Overall, these results reveal molecular mechanisms of CD treatment based on simultaneous effects on the pathogen and on host small molecule responses, and expand our understanding of clinical treatment failure in CD. This link between infection and subsequent persistent small molecule perturbation broadens our understanding of infectious disease sequelae.

Integration of Untargeted Metabolomics and Microbial Community Analyses to Characterize Distinct Deep-Sea Methane Seeps

Deep-sea methane seeps host highly diverse microbial communities whose biological diversity is distinct from other marine habitats. Coupled with microbial community analysis, untargeted metabolomics of environmental samples using high resolution tandem mass spectrometry provides unprecedented access to the unique specialized metabolisms of these chemosynthetic microorganisms. In addition, the diverse microbial natural products are of broad interest due to their potential applications for human and environmental health and well-being. In this exploratory study, sediment cores were collected from two methane seeps (-1000 m water depth) with very different gross geomorphologies, as well as a non-seep control site. Cores were subjected to parallel metabolomic and microbial community analyses to assess the feasibility of representative metabolite detection and identify congruent patterns between metabolites and microbes. Metabolomes generated using high resolution liquid chromatography tandem mass spectrometry were annotated with predicted structure classifications of the majority of mass features using SIRIUS and CANOPUS. The microbiome was characterized by analysis of 16S rRNA genes and analyzed both at the whole community level, as well as the small subgroup of Actinobacteria, which are known to produce societally useful compounds. Overall, the younger Dagorlad seep possessed a greater abundance of metabolites while there was more variation in abundance, number, and distribution of metabolites between samples at the older Emyn Muil seep. Lipid and lipid-like molecules displayed the greatest variation between sites and accounted for a larger proportion of metabolites found at the older seep. Overall, significant differences in composition of the microbial community mirrored the patterns of metabolite diversity within the samples; both varied greatly as a function of distance from methane seep, indicating a deterministic role of seepage. Interdisciplinary research to understand microbial and metabolic diversity is essential for understanding the processes and role of ubiquitous methane seeps in global systems and here we increase understanding of these systems by visualizing some of the chemical diversity that seeps add to marine systems.

Illuminating the oral microbiome and its host interactions: recent advancements in omics and bioinformatics technologies in the context of oral microbiome research

The oral microbiota has an enormous impact on human health, with oral dysbiosis now linked to many oral and systemic diseases. Recent advancements in sequencing, mass spectrometry, bioinformatics, computational biology, and machine learning are revolutionizing oral microbiome research, enabling analysis at an unprecedented scale and level of resolution using omics approaches. This review contains a comprehensive perspective of the current state-of-the-art tools available to perform genomics, metagenomics, phylogenomics, pangenomics, transcriptomics, proteomics, metabolomics, lipidomics, and multi-omics analysis on (all) microbiomes, and then provides examples of how the techniques have been applied to research of the oral microbiome, specifically. Key findings of these studies and remaining challenges for the field are highlighted. Although the methods discussed here are placed in the context of their contributions to oral microbiome research specifically, they are pertinent to the study of any microbiome, and the intended audience of this includes researchers would simply like to get an introduction to microbial omics and/or an update on the latest omics methods. Continued research of the oral microbiota using omics approaches is crucial and will lead to dramatic improvements in human health, longevity, and quality of life.

Computation-aided studies related to the induction of specialized metabolite biosynthesis in microbial co-cultures: An introductory overview

Co-cultivation is an effective method of inducing the production of specialized metabolites (SMs) in microbial strains. By mimicking the ecological interactions that take place in natural environment, this approach enables to trigger the biosynthesis of molecules which are not formed under monoculture conditions. Importantly, microbial co-cultivation may lead to the discovery of novel chemical entities of pharmaceutical interest. The experimental efforts aimed at the induction of SMs are greatly facilitated by computational techniques. The aim of this overview is to highlight the relevance of computational methods for the investigation of SM induction via microbial co-cultivation. The concepts related to the induction of SMs in microbial co-cultures are briefly introduced by addressing four areas associated with the SM induction workflows, namely the detection of SMs formed exclusively under co-culture conditions, the annotation of induced SMs, the identification of SM producer strains, and the optimization of fermentation conditions. The computational infrastructure associated with these areas, including the tools of multivariate data analysis, molecular networking, genome mining and mathematical optimization, is discussed in relation to the experimental results described in recent literature. The perspective on the future developments in the field, mainly in relation to the microbiome-related research, is also provided.

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.

Metabolome-guided genome mining of RiPP natural products

Ribosomally synthesized and post-translationally modified peptides (RiPPs) are a chemically diverse class of metabolites. Many RiPPs show potent biological activities that make them attractive starting points for drug development. A promising approach for the discovery of new classes of RiPPs is genome mining. However, the accuracy of genome mining is hampered by the lack of signature genes shared across different RiPP classes. One way to reduce false-positive predictions is by complementing genomic information with metabolomics data. In recent years, several new approaches addressing such integrative genomics and metabolomics analyses have been developed. In this review, we provide a detailed discussion of RiPP-compatible software tools that integrate paired genomics and metabolomics data. We highlight current challenges in data integration and identify opportunities for further developments targeting new classes of bioactive RiPPs.

Recent advances in mass spectrometry-based computational metabolomics

The computational metabolomics field brings together computer scientists, bioinformaticians, chemists, clinicians, and biologists to maximize the impact of metabolomics across a wide array of scientific and medical disciplines. The field continues to expand as modern instrumentation produces datasets with increasing complexity, resolution, and sensitivity. These datasets must be processed, annotated, modeled, and interpreted to enable biological insight. Techniques for visualization, integration (within or between omics), and interpretation of metabolomics data have evolved along with innovation in the databases and knowledge resources required to aid understanding. In this review, we highlight recent advances in the field and reflect on opportunities and innovations in response to the most pressing challenges. This review was compiled from discussions from the 2022 Dagstuhl seminar entitled "Computational Metabolomics: From Spectra to Knowledge".

MS2Query: reliable and scalable MS2 mass spectra-based analogue search

Metabolomics-driven discoveries of biological samples remain hampered by the grand challenge of metabolite annotation and identification. Only few metabolites have an annotated spectrum in spectral libraries; hence, searching only for exact library matches generally returns a few hits. An attractive alternative is searching for so-called analogues as a starting point for structural annotations; analogues are library molecules which are not exact matches but display a high chemical similarity. However, current analogue search implementations are not yet very reliable and relatively slow. Here, we present MS2Query, a machine learning-based tool that integrates mass spectral embedding-based chemical similarity predictors (Spec2Vec and MS2Deepscore) as well as detected precursor masses to rank potential analogues and exact matches. Benchmarking MS2Query on reference mass spectra and experimental case studies demonstrate improved reliability and scalability. Thereby, MS2Query offers exciting opportunities to further increase the annotation rate of metabolomics profiles of complex metabolite mixtures and to discover new biology.