Automated Annotation of Untargeted All-Ion Fragmentation LC-MS Metabolomics Data with MetaboAnnotatoR

B. Trivella D. NP3 MS Workflow: An Open-Source Software System to Empower Natural Product-Based Drug Discovery Using Untargeted Metabolomics

Natural products (or specialized metabolites) are historically the main source of new drugs. However, the current drug discovery pipelines require miniaturization and speeds that are incompatible with traditional natural product research methods, especially in the early stages of the research. This article introduces the NP3 MS Workflow, a robust open-source software system for liquid chromatography-tandem mass spectrometry (LC-MS/MS) untargeted metabolomic data processing and analysis, designed to rank bioactive natural products directly from complex mixtures of compounds, such as bioactive biota samples. NP3 MS Workflow allows minimal user intervention as well as customization of each step of LC-MS/MS data processing, with diagnostic statistics to allow interpretation and optimization of LC-MS/MS data processing by the user. NP3 MS Workflow adds improved computing of the MS2 spectra in an LC-MS/MS data set and provides tools for automatic [M + H]+ ion deconvolution using fragmentation rules; chemical structural annotation against MS2 databases; and relative quantification of the precursor ions for bioactivity correlation scoring. The software will be presented with case studies and comparisons with equivalent tools currently available. NP3 MS Workflow shows a robust and useful approach to select bioactive natural products from complex mixtures, improving the set of tools available for untargeted metabolomics. It can be easily integrated into natural product-based drug-discovery pipelines and to other fields of research at the interface of chemistry and biology.

MetaboAnalystR 4.0: a unified LC-MS workflow for global metabolomics

The wide applications of liquid chromatography - mass spectrometry (LC-MS) in untargeted metabolomics demand an easy-to-use, comprehensive computational workflow to support efficient and reproducible data analysis. However, current tools were primarily developed to perform specific tasks in LC-MS based metabolomics data analysis. Here we introduce MetaboAnalystR 4.0 as a streamlined pipeline covering raw spectra processing, compound identification, statistical analysis, and functional interpretation. The key features of MetaboAnalystR 4.0 includes an auto-optimized feature detection and quantification algorithm for LC-MS1 spectra processing, efficient MS2 spectra deconvolution and compound identification for data-dependent or data-independent acquisition, and more accurate functional interpretation through integrated spectral annotation. Comprehensive validation studies using LC-MS1 and MS2 spectra obtained from standards mixtures, dilution series and clinical metabolomics samples have shown its excellent performance across a wide range of common tasks such as peak picking, spectral deconvolution, and compound identification with good computing efficiency. Together with its existing statistical analysis utilities, MetaboAnalystR 4.0 represents a significant step toward a unified, end-to-end workflow for LC-MS based global metabolomics in the open-source R environment.

Addressing the untargeted lipidomics challenge in urine samples: Comparative study of extraction methods by UHPLC-ESI-QTOF-MS

Urine analysis has remained a fundamental and widely used method in clinical diagnostics for over a century. With its minimal invasive nature and comprehensive range of analytes, urine has established itself as a clinical diagnostic tool for various disorders, including renal, urological, metabolic, and endocrine diseases. Furthermore, urine's unique attributes make it an attractive matrix for biomarker discovery, as well as in assessing the metabolic and physiological states of patients and healthy individuals alike. However, limitations in our knowledge of average values and sources of urinary lipids decrease the wider clinical application of urinary lipidomics. In this context, untargeted lipidomics analysis relies heavily on the extraction and analysis of lipids in biological samples. Nevertheless, this type of analysis presents challenges in lipid identification due to the diverse nature of lipids. Therefore, proper sample treatment before analysis is crucial to obtain robust and reproducible lipidomic profiles. To address this gap, we conducted a comparative study of a urine pool sample collected from twenty healthy volunteers using four different lipid extraction methods: one biphasic and three monophasic protocols. The extracted lipids were then analyzed using UHPLC-MS and MS/MS, and the semi-quantification of all the accurately annotated lipid species was performed for each extraction method.

LC-MS/DIA-based strategy for comprehensive flavonoid profiling: an Ocotea spp. applicability case

We introduce a liquid chromatography - mass spectrometry with data-independent acquisition (LC-MS/DIA)-based strategy, specifically tailored to achieve comprehensive and reliable glycosylated flavonoid profiling. This approach facilitates in-depth and simultaneous exploration of all detected precursors and fragments during data processing, employing the widely-used open-source MZmine 3 software. It was applied to a dataset of six Ocotea plant species. This framework suggested 49 flavonoids potentially newly described for these plant species, alongside 45 known features within the genus. Flavonols kaempferol and quercetin, both exhibiting O-glycosylation patterns, were particularly prevalent. Gas-phase fragmentation reactions further supported these findings. For the first time, the apigenin flavone backbone was also annotated in most of the examined Ocotea species. Apigenin derivatives were found mainly in the C-glycoside form, with O. porosa displaying the highest flavone : flavonol ratio. The approach also allowed an unprecedented detection of kaempferol and quercetin in O. porosa species, and it has underscored the untapped potential of LC-MS/DIA data for broad and reliable flavonoid profiling. Our study annotated more than 50 flavonoid backbones in each species, surpassing the current literature.

Enhanced efficiency of MS/MS all-ion fragmentation for non-targeted analysis of trace contaminants in surface water using multivariate curve resolution and data fusion

Data-independent acquisition-all-ion fragmentation (DIA-AIF) mode of mass spectrometry can facilitate wide-scope non-target analysis of contaminants in surface water due to comprehensive spectral identification. However, because of the complexity of the resulting MS2 AIF spectra, identifying unknown pollutants remains a significant challenge, with a significant bottleneck in translating non-targeted chemical signatures into environmental impacts. The present study proposes to process fused MS1 and MS2 data sets obtained from LC-HRMS/MS measurements in non-targeted AIF workflows on surface water samples using multivariate curve resolution-alternating least squares (MCR-ALS). This enables straightforward assignment between precursor ions obtained from resolved MS1 spectra and their corresponding MS2 spectra. The method was evaluated for two sets of tap water and surface water contaminated with 14 target chemicals as a proof of concept. The data set of surface water samples consisting of 3506 MS1 and 2170 MS2 AIF mass spectral features was reduced to 81 components via a fused MS1-MS2 MCR model that describes at least 98.8% of the data. Each component summarizes the distinct chromatographic elution of components together with their corresponding MS1 and MS2 spectra. MS2 spectral similarity of more than 82% was obtained for most target chemicals. This highlights the potential of this method for unraveling the composition of MS/MS complex data in a water environment. Ultimately, the developed approach was applied to the retrospective non-target analysis of an independent set of surface water samples.

Reactivity-based identification of oxygen containing functional groups of chemicals applied as potential classifier in non-target analysis

In this work, we developed a reactivity-based strategy to identify functional groups of unknown analytes, which can be applied as classifier in non-target analysis with gas chromatography. The aim of this strategy is to reduce the number of potential candidate structures generated for a molecular formula determined by high resolution mass spectrometry. We selected an example of 18 isomers with the molecular formula C12H10O2 to test the performance of different derivatization reagents, whereas our aim was to select mild and fast reaction conditions. Based on the results for the isomers, we developed a four-step workflow for the identification of functional groups containing oxygen.

Cumulative Neutral Loss Model for Fragment Deconvolution in Electrospray Ionization High-Resolution Mass Spectrometry Data

Clean high-resolution mass spectra (HRMS) are essential to a successful structural elucidation of an unknown feature during nontarget analysis (NTA) workflows. This is a crucial step, particularly for the spectra generated during data-independent acquisition or during direct infusion experiments. The most commonly available tools only take advantage of the time domain for spectral cleanup. Here, we present an algorithm that combines the time domain and mass domain information to perform spectral deconvolution. The algorithm employs a probability-based cumulative neutral loss (CNL) model for fragment deconvolution. The optimized model, with a mass tolerance of 0.005 Da and a scoreCNL threshold of 0.00, was able to achieve a true positive rate (TPr) of 95.0%, a false discovery rate (FDr) of 20.6%, and a reduction rate of 35.4%. Additionally, the CNL model was extensively tested on real samples containing predominantly pesticides at different concentration levels and with matrix effects. Overall, the model was able to obtain a TPr above 88.8% with FD rates between 33 and 79% and reduction rates between 9 and 45%. Finally, the CNL model was compared with the retention time difference method and peak shape correlation analysis, showing that a combination of correlation analysis and the CNL model was the most effective for fragment deconvolution, obtaining a TPr of 84.7%, an FDr of 54.4%, and a reduction rate of 51.0%.

IDSL.CSA: Composite Spectra Analysis for Chemical Annotation of Untargeted Metabolomics Datasets

Poor chemical annotation of high-resolution mass spectrometry data limits applications of untargeted metabolomics datasets. Our new software, the Integrated Data Science Laboratory for Metabolomics and Exposomics─Composite Spectra Analysis (IDSL.CSA) R package, generates composite mass spectra libraries from MS1-only data, enabling the chemical annotation of high-resolution mass spectrometry coupled with liquid chromatography peaks regardless of the availability of MS2 fragmentation spectra. We demonstrate comparable annotation rates for commonly detected endogenous metabolites in human blood samples using IDSL.CSA libraries versus MS/MS libraries in validation tests. IDSL.CSA can create and search composite spectra libraries from any untargeted metabolomics dataset generated using high-resolution mass spectrometry coupled to liquid or gas chromatography instruments. The cross-applicability of these libraries across independent studies may provide access to new biological insights that may be missed due to the lack of MS2 fragmentation data. The IDSL.CSA package is available in the R-CRAN repository at https://cran.r-project.org/package=IDSL.CSA. Detailed documentation and tutorials are provided at https://github.com/idslme/IDSL.CSA.

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".

Oligonucleotide Catabolism-Derived Gluconucleosides in Caenorhabditis elegans

Nucleosides are essential cornerstones of life, and nucleoside derivatives and synthetic analogues have important biomedical applications. Correspondingly, production of non-canonical nucleoside derivatives in animal model systems is of particular interest. Here, we report the discovery of diverse glucose-based nucleosides in Caenorhabditis elegans and related nematodes. Using a mass spectrometric screen based on all-ion fragmentation in combination with total synthesis, we show that C. elegans selectively glucosylates a series of modified purines but not the canonical purine and pyrimidine bases. Analogous to ribonucleosides, the resulting gluconucleosides exist as phosphorylated and non-phosphorylated forms. The phosphorylated gluconucleosides can be additionally decorated with diverse acyl moieties from amino acid catabolism. Syntheses of representative variants, facilitated by a novel 2'-O- to 3'-O-dibenzyl phosphoryl transesterification reaction, demonstrated selective incorporation of different nucleobases and acyl moieties. Using stable-isotope labeling, we further show that gluconucleosides incorporate modified nucleobases derived from RNA and possibly DNA breakdown, revealing extensive recycling of oligonucleotide catabolites. Gluconucleosides are conserved in other nematodes, and biosynthesis of specific subsets is increased in germline mutants and during aging. Bioassays indicate that gluconucleosides may function in stress response pathways.

IDSL.CSA: Composite Spectra Analysis for Chemical Annotation of Untargeted Metabolomics Datasets

Poor chemical annotation of high-resolution mass spectrometry data limit applications of untargeted metabolomics datasets. Our new software, the Integrated Data Science Laboratory for Metabolomics and Exposomics - Composite Spectra Analysis (IDSL.CSA) R package, generates composite mass spectra libraries from MS1-only data, enabling the chemical annotation of LC/HRMS peaks regardless of the availability of MS2 fragmentation spectra. We demonstrate comparable annotation rates for commonly detected endogenous metabolites in human blood samples using IDSL.CSA libraries versus MS/MS libraries in validation tests. IDSL.CSA can create and search composite spectra libraries from any untargeted metabolomics dataset generated using high-resolution mass spectrometry coupled to liquid or gas chromatography instruments. The cross-applicability of these libraries across independent studies may provide access to new biological insights that may be missed due to the lack of MS2 fragmentation data. The IDSL.CSA package is available in the R CRAN repository at https://cran.r-project.org/package=IDSL.CSA . Detailed documentation and tutorials are provided at https://github.com/idslme/IDSL.CSA .

For Table of Contents Only

Factors affecting untargeted detection of doping agents in biological samples

Doping control is essential for sports, and untargeted detection of doping agents (UDDA) is the holy grail for anti-doping strategies. The present study examined major factors impacting UDDA with metabolomic data processing, including the use of blank samples, signal-to-noise ratio thresholds, and the minimum chromatographic peak intensity. Contrary to data processing in metabolomics studies, both blank sample use (either blank solvent or plasma) and marking of background compounds were found to be unnecessary for UDDA in biological samples, the first such report to the authors' knowledge. The minimum peak intensity required to detect chromatographic peaks affected the limit of detection (LOD) and data processing time for untargeted detection of 57 drugs spiked into equine plasma. The ratio of the mean (ROM) of the extracted ion chromatographic peak area of a compound in the sample group (SG) to that in the control group (CG) impacted its LOD, and a small ROM value such as 2 is recommended for UDDA. Mathematical modeling of the required signal-to-noise ratio (S/N) for UDDA provided insights into the effect of the number of samples in the SG, the number of positive samples, and the ROM on the required S/N, highlighting the power of mathematics in addressing issues in analytical chemistry. The UDDA method was validated by its successful identification of untargeted doping agents in real-world post-competition equine plasma samples. This advancement in UDDA methodology will be a useful addition to the arsenal of approaches used to combat doping in sports.

MAW: the reproducible Metabolome Annotation Workflow for untargeted tandem mass spectrometry

Mapping the chemical space of compounds to chemical structures remains a challenge in metabolomics. Despite the advancements in untargeted liquid chromatography-mass spectrometry (LC-MS) to achieve a high-throughput profile of metabolites from complex biological resources, only a small fraction of these metabolites can be annotated with confidence. Many novel computational methods and tools have been developed to enable chemical structure annotation to known and unknown compounds such as in silico generated spectra and molecular networking. Here, we present an automated and reproducible Metabolome Annotation Workflow (MAW) for untargeted metabolomics data to further facilitate and automate the complex annotation by combining tandem mass spectrometry (MS²) input data pre-processing, spectral and compound database matching with computational classification, and in silico annotation. MAW takes the LC-MS² spectra as input and generates a list of putative candidates from spectral and compound databases. The databases are integrated via the R package Spectra and the metabolite annotation tool SIRIUS as part of the R segment of the workflow (MAW-R). The final candidate selection is performed using the cheminformatics tool RDKit in the Python segment (MAW-Py). Furthermore, each feature is assigned a chemical structure and can be imported to a chemical structure similarity network. MAW is following the FAIR (Findable, Accessible, Interoperable, Reusable) principles and has been made available as the docker images, maw-r and maw-py. The source code and documentation are available on GitHub ( https://github.com/zmahnoor14/MAW ). The performance of MAW is evaluated on two case studies. MAW can improve candidate ranking by integrating spectral databases with annotation tools like SIRIUS which contributes to an efficient candidate selection procedure. The results from MAW are also reproducible and traceable, compliant with the FAIR guidelines. Taken together, MAW could greatly facilitate automated metabolite characterization in diverse fields such as clinical metabolomics and natural product discovery.

Simulated-to-real benchmarking of acquisition methods in untargeted metabolomics

Data-Dependent and Data-Independent Acquisition modes (DDA and DIA, respectively) are both widely used to acquire MS2 spectra in untargeted liquid chromatography tandem mass spectrometry (LC-MS/MS) metabolomics analyses. Despite their wide use, little work has been attempted to systematically compare their MS/MS spectral annotation performance in untargeted settings due to the lack of ground truth and the costs involved in running a large number of acquisitions. Here, we present a systematic in silico comparison of these two acquisition methods in untargeted metabolomics by extending our Virtual Metabolomics Mass Spectrometer (ViMMS) framework with a DIA module. Our results show that the performance of these methods varies with the average number of co-eluting ions as the most important factor. At low numbers, DIA outperforms DDA, but at higher numbers, DDA has an advantage as DIA can no longer deal with the large amount of overlapping ion chromatograms. Results from simulation were further validated on an actual mass spectrometer, demonstrating that using ViMMS we can draw conclusions from simulation that translate well into the real world. The versatility of the Virtual Metabolomics Mass Spectrometer (ViMMS) framework in simulating different parameters of both Data-Dependent and Data-Independent Acquisition (DDA and DIA) modes is a key advantage of this work. Researchers can easily explore and compare the performance of different acquisition methods within the ViMMS framework, without the need for expensive and time-consuming experiments with real experimental data. By identifying the strengths and limitations of each acquisition method, researchers can optimize their choice and obtain more accurate and robust results. Furthermore, the ability to simulate and validate results using the ViMMS framework can save significant time and resources, as it eliminates the need for numerous experiments. This work not only provides valuable insights into the performance of DDA and DIA, but it also opens the door for further advancements in LC-MS/MS data acquisition methods.

NMR-Based Chromatography Readouts: Indispensable Tools to “Translate” Analytical Features into Molecular Structures

Gaining structural information is a must to allow the unequivocal structural characterization of analytes from natural sources. In liquid state, NMR spectroscopy is almost the only possible alternative to HPLC-MS and hyphenating the effluent of an analyte separation device to the probe head of an NMR spectrometer has therefore been pursued for more than three decades. The purpose of this review article was to demonstrate that, while it is possible to use mass spectrometry and similar methods to differentiate, group, and often assign the differentiating variables to entities that can be recognized as single molecules, the structural characterization of these putative biomarkers usually requires the use of NMR spectroscopy.

Antibacterial Activity from Momordica charantia L. Leaves and Flavones Enriched Phase

Momordica charantia L. (Cucurbitaceae) is a plant known in Brazil as “melão de São Caetano”, which has been related to many therapeutic applications in folk medicine. Herein, we describe antibacterial activities and related metabolites for an extract and fractions obtained from the leaves of that species. An ethanolic extract and its three fractions were used to perform in vitro antibacterial assays. In addition, liquid chromatography coupled to mass spectrometry and the molecular networking approach were used for the metabolite annotation process. Overall, 25 compounds were annotated in the ethanolic extract from M. charantia leaves, including flavones, terpenes, organic acids, and inositol pyrophosphate derivatives. The ethanolic extract exhibited low activity against Proteus mirabilis (MIC 312.5 µg·mL⁻¹) and Klebsiella pneumoniae (MIC 625 µg·mL⁻¹). The ethyl acetate phase showed interesting antibacterial activity (MIC 156.2 µg·mL⁻¹) against Klebsiella pneumoniae, and it was well justified by the high content of glycosylated flavones. Therefore, based on the ethyl acetate phase antibacterial result, we suggest that M. charantia leaves could be considered as an alternative antibacterial source against K. pneumoniae and can serve as a pillar for future studies as well as pharmacological application against the bacteria.