Neutral Loss Mass Spectral Data Enhances Molecular Similarity Analysis in METLIN

Continuous collective analysis of chemical reactions

The automated synthesis of small organic molecules from modular building blocks has the potential to transform our capacity to create medicines and materials1-3. Disruptive acceleration of this molecule-building strategy broadly unlocks its functional potential and requires the integration of many new assembly chemistries. Although recent advances in high-throughput chemistry4-6 can speed up the development of appropriate synthetic methods, for example, in selecting appropriate chemical reaction conditions from the vast range of potential options, equivalent high-throughput analytical methods are needed. Here we report a streamlined approach for the rapid, quantitative analysis of chemical reactions by mass spectrometry. The intrinsic fragmentation features of chemical building blocks generalize the analyses of chemical reactions, allowing sub-second readouts of reaction outcomes. Central to this advance was identifying that starting material fragmentation patterns function as universal barcodes for downstream product analysis by mass spectrometry. Combining these features with acoustic droplet ejection mass spectrometry7,8 we could eliminate slow chromatographic steps and continuously evaluate chemical reactions in multiplexed formats. This enabled the assignment of reaction conditions to molecules derived from ultrahigh-throughput chemical synthesis experiments. More generally, these results indicate that fragmentation features inherent to chemical synthesis can empower rapid data-rich experimentation.

LC-HRMS/MS-Guided Isolation of Unusual Diarylheptanoids from the Rhizomes of Alpinia officinarum

LC-HRMS/MS analysis facilitated the precise targeting, isolation, and identification of unusual dimeric diarylheptanoids from Alpinia officinarum (A. officinarum). The tandem MS data for (4E)-1,7-diphenyl-4-hepten-3-one (7) revealed fragment ions at m/z 91, 105, and 117, which are fragmentation patterns specific to diarylheptanoids. In the tandem MS data, peaks with m/z values ranging from 450 to 600 that exhibited these specific fragment ions were selected and isolated. Consequently, two previously undescribed dimeric diarylheptanoids (1 and 2) and four unusual diarylheptanoids (3-6) along with 10 monomeric diarylheptanoids (7-16) were isolated from the rhizomes of A. officinarum using various chromatographic techniques. The structures of the isolates were elucidated by an analysis of 1D/2D NMR and HRESIMS data, and a combination of DP4+ probability analysis and ECD calculations. To evaluate the anti-inflammatory effects of the isolated compounds, their inhibitory activity against nitric oxide production in LPS-induced RAW 264.7 cells was assessed. Compounds 1, 7, and 9 exhibited remarkable inhibitory effects with IC50 values of 14.7, 6.6, and 5.0 μM, respectively.

Substrate Specificity of a Methyltransferase Involved in the Biosynthesis of the Lantibiotic Cacaoidin

Modification of the N- and C-termini of peptides enhances their stability against degradation by exopeptidases. The biosynthetic pathways of many peptidic natural products feature enzymatic modification of their termini, and these enzymes may represent a valuable pool of biocatalysts. The lantibiotic cacaoidin carries an N,N-dimethylated N-terminal amine group. Its biosynthetic gene cluster encodes the putative methyltransferase Cao4. In this work, we present reconstitution of the activity of the enzyme, which we termed CaoSC following standardized lanthipeptide nomenclature, using a heterologously produced peptide as the model substrate. In vitro methylation of diverse lanthipeptides revealed the substrate requirements of CaoSC. The enzyme accepts peptides of varying lengths and C-terminal sequences but requires dehydroalanine or dehydrobutyrine at the second position. CaoSC-mediated dimethylation of natural lantibiotics resulted in modestly enhanced antimicrobial activity of the lantibiotic haloduracin compared to that of the native compound. Improved activity and/or metabolic stability as a result of methylation illustrates the potential future application of CaoSC in the bioengineering of therapeutic peptides.

MetaboAnalyst 6.0: towards a unified platform for metabolomics data processing, analysis and interpretation

We introduce MetaboAnalyst version 6.0 as a unified platform for processing, analyzing, and interpreting data from targeted as well as untargeted metabolomics studies using liquid chromatography - mass spectrometry (LC-MS). The two main objectives in developing version 6.0 are to support tandem MS (MS2) data processing and annotation, as well as to support the analysis of data from exposomics studies and related experiments. Key features of MetaboAnalyst 6.0 include: (i) a significantly enhanced Spectra Processing module with support for MS2 data and the asari algorithm; (ii) a MS2 Peak Annotation module based on comprehensive MS2 reference databases with fragment-level annotation; (iii) a new Statistical Analysis module dedicated for handling complex study design with multiple factors or phenotypic descriptors; (iv) a Causal Analysis module for estimating metabolite - phenotype causal relations based on two-sample Mendelian randomization, and (v) a Dose-Response Analysis module for benchmark dose calculations. In addition, we have also improved MetaboAnalyst's visualization functions, updated its compound database and metabolite sets, and significantly expanded its pathway analysis support to around 130 species. MetaboAnalyst 6.0 is freely available at https://www.metaboanalyst.ca.

TransExION: a transformer based explainable similarity metric for comparing IONS in tandem mass spectrometry

Small molecule identification is a crucial task in analytical chemistry and life sciences. One of the most commonly used technologies to elucidate small molecule structures is mass spectrometry. Spectral library search of product ion spectra (MS/MS) is a popular strategy to identify or find structural analogues. This approach relies on the assumption that spectral similarity and structural similarity are correlated. However, popular spectral similarity measures, usually calculated based on identical fragment matches between the MS/MS spectra, do not always accurately reflect the structural similarity. In this study, we propose TransExION, a Transformer based Explainable similarity metric for IONS. TransExION detects related fragments between MS/MS spectra through their mass difference and uses these to estimate spectral similarity. These related fragments can be nearly identical, but can also share a substructure. TransExION also provides a post-hoc explanation of its estimation, which can be used to support scientists in evaluating the spectral library search results and thus in structure elucidation of unknown molecules. Our model has a Transformer based architecture and it is trained on the data derived from GNPS MS/MS libraries. The experimental results show that it improves existing spectral similarity measures in searching and interpreting structural analogues as well as in molecular networking. SCIENTIFIC CONTRIBUTION: We propose a transformer-based spectral similarity metrics that improves the comparison of small molecule tandem mass spectra. We provide a post hoc explanation that can serve as a good starting point for unknown spectra annotation based on database spectra.

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.

Weighting Low-Intensity MS/MS Ions and m/z Frequency for Spectral Library Annotation

Calculating spectral similarity is a fundamental step in MS/MS data analysis in untargeted metabolomics experiments, as it facilitates the identification of related spectra and the annotation of compounds. To improve matching accuracy when querying an experimental mass spectrum against a spectral library, previous approaches have proposed increasing peak intensities for high m/z ranges. These high m/z values tend to be smaller in magnitude, yet they offer more crucial information for identifying the chemical structure. Here, we evaluate the impact of using these weights for identifying structurally related compounds and mass spectral library searches. Additionally, we propose a weighting approach that (i) takes into account the frequency of the m/z values within a spectral library in order to assign higher importance to the most common peaks and (ii) increases the intensity of lower peaks, similar to previous approaches. To demonstrate our approach, we applied weighting preprocessing to modified cosine, entropy, and fidelity distance metrics and benchmarked it against previously reported weights. Our results demonstrate how weighting-based preprocessing can assist in annotating the structure of unknown spectra as well as identifying structurally similar compounds. Finally, we examined scenarios in which the utilization of weights resulted in diminished performance, pinpointing spectral features where the application of weights might be detrimental.

Open access repository-scale propagated nearest neighbor suspect spectral library for untargeted metabolomics

Despite the increasing availability of tandem mass spectrometry (MS/MS) community spectral libraries for untargeted metabolomics over the past decade, the majority of acquired MS/MS spectra remain uninterpreted. To further aid in interpreting unannotated spectra, we created a nearest neighbor suspect spectral library, consisting of 87,916 annotated MS/MS spectra derived from hundreds of millions of MS/MS spectra originating from published untargeted metabolomics experiments. Entries in this library, or "suspects," were derived from unannotated spectra that could be linked in a molecular network to an annotated spectrum. Annotations were propagated to unknowns based on structural relationships to reference molecules using MS/MS-based spectrum alignment. We demonstrate the broad relevance of the nearest neighbor suspect spectral library through representative examples of propagation-based annotation of acylcarnitines, bacterial and plant natural products, and drug metabolism. Our results also highlight how the library can help to better understand an Alzheimer's brain phenotype. The nearest neighbor suspect spectral library is openly available for download or for data analysis through the GNPS platform to help investigators hypothesize candidate structures for unknown MS/MS spectra in untargeted metabolomics data.

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.

The Integration of the Metabolome and Transcriptome for Dendrobium nobile Lindl. in Response to Methyl Jasmonate

Dendrobium nobile Lindl., as an endangered medicinal plant within the genus Dendrobium, is widely distributed in southwestern China and has important ecological and economic value. There are a variety of metabolites with pharmacological activity in D. nobile. The alkaloids and polysaccharides contained within D. nobile are very important active components, which mainly have antiviral, anti-tumor, and immunity improvement effects. However, the changes in the compounds and functional genes of D. nobile induced by methyl jasmonate (MeJA) are not clearly understood. In this study, the metabolome and transcriptome of D. nobile were analyzed after exposure to MeJA. A total of 377 differential metabolites were obtained through data analysis, of which 15 were related to polysaccharide pathways and 35 were related to terpenoids and alkaloids pathways. Additionally, the transcriptome sequencing results identified 3256 differentially expressed genes that were discovered in 11 groups. Compared with the control group, 1346 unigenes were differentially expressed in the samples treated with MeJA for 14 days (TF14). Moreover, the expression levels of differentially expressed genes were also significant at different growth and development stages. According to GO and KEGG annotations, 189 and 99 candidate genes were identified as being involved in terpenoid biosynthesis and polysaccharide biosynthesis, respectively. In addition, the co-expression analysis indicated that 238 and 313 transcription factors (TFs) may contribute to the regulation of terpenoid and polysaccharide biosynthesis, respectively. Through a heat map analysis, fourteen terpenoid synthetase genes, twenty-three cytochrome P450 oxidase genes, eight methyltransferase genes, and six aminotransferase genes were identified that may be related to dendrobine biosynthesis. Among them, one sesquiterpene synthase gene was found to be highly expressed after the treatment with MeJA and was positively correlated with the content of dendrobine. This study provides important and valuable metabolomics and transcriptomic information for the further understanding of D. nobile at the metabolic and molecular levels and provides candidate genes and possible intermediate compounds for the dendrobine biosynthesis pathway, which lays a certain foundation for further research on and application of Dendrobium.

Flash entropy search to query all mass spectral libraries in real time

Public repositories of metabolomics mass spectra encompass more than 1 billion entries. With open search, dot product or entropy similarity, comparisons of a single tandem mass spectrometry spectrum take more than 8 h. Flash entropy search speeds up calculations more than 10,000 times to query 1 billion spectra in less than 2 s, without loss in accuracy. It benefits from using multiple threads and GPU calculations. This algorithm can fully exploit large spectral libraries with little memory overhead for any mass spectrometry laboratory.

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

Identification of Pseudomonas aeruginosa exopolysaccharide Psl in biofilms using 3D OrbiSIMS

Secondary ion mass spectrometry (SIMS) offers advantages over both liquid extraction mass spectrometry and matrix assisted laser desorption mass spectrometry in that it provides the direct in situ analysis of molecules and has the potential to preserve the 3D location of an analyte in a sample. Polysaccharides are recognized as challenging analytes in the mass spectrometry of liquids and are also difficult to identify and assign using SIMS. Psl is an exopolysaccharide produced by Pseudomonas aeruginosa, which plays a key role in biofilm formation and maturation. In this Letter, we describe the use of the OrbiTrap analyzer with SIMS (3D OrbiSIMS) for the label-free mass spectrometry of Psl, taking advantage of its high mass resolving power for accurate secondary ion assignment. We study a P. aeruginosa biofilm and compare it with purified Psl to enable the assignment of secondary ions specific to the Psl structure. This resulted in the identification of 17 peaks that could confidently be ascribed to Psl fragments within the biofilm matrix. The complementary approach of the following neutral loss sequences is also shown to identify multiple oligosaccharide fragments without the requirement of a biological reference sample.

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

BACKGROUND

Untargeted metabolomics approaches based on mass spectrometry obtain comprehensive profiles of complex biological samples. However, on average only 10% of the molecules can be annotated. This low annotation rate hampers biochemical interpretation and effective comparison of metabolomics studies. Furthermore, de novo structural characterization of mass spectral data remains a complicated and time-intensive process. Recently, the field of computational metabolomics has gained traction and novel methods have started to enable large-scale and reliable metabolite annotation. Molecular networking and machine learning-based in-silico annotation tools have been shown to greatly assist metabolite characterization in diverse fields such as clinical metabolomics and natural product discovery.

AIM OF REVIEW

We highlight recent advances in computational metabolite annotation workflows with a special focus on their evaluation and comparison with other tools. Whilst the progress is substantial and promising, we also argue that inconsistencies in benchmarking different tools hamper users from selecting the most appropriate and promising method for their research. We summarize benchmarking strategies of the different tools and outline several recommendations for benchmarking and comparing novel tools.

KEY SCIENTIFIC CONCEPTS OF REVIEW

This review focuses on recent advances in mass spectral library-based and machine learning-supported metabolite annotation workflows. We discuss large-scale library matching and analogue search, the current bloom of mass spectral similarity scores, and how molecular networking has changed the field. In addition, the potentials and challenges of machine learning-supported metabolite annotation workflows are highlighted. Overall, recent developments in computational metabolomics have started to fundamentally change metabolomics workflows, and we expect that as a community we will be able to overcome current method performance ambiguities and annotation bottlenecks.

Comparison of Cosine, Modified Cosine, and Neutral Loss Based Spectrum Alignment For Discovery of Structurally Related Molecules

Spectrum alignment of tandem mass spectrometry (MS/MS) data using the modified cosine similarity and subsequent visualization as molecular networks have been demonstrated to be a useful strategy to discover analogs of molecules from untargeted MS/MS-based metabolomics experiments. Recently, a neutral loss matching approach has been introduced as an alternative to MS/MS-based molecular networking with an implied performance advantage in finding analogs that cannot be discovered using existing MS/MS spectrum alignment strategies. To comprehensively evaluate the scoring properties of neutral loss matching, the cosine similarity, and the modified cosine similarity, similarity measures of 955 228 peptide MS/MS spectrum pairs and 10 million small molecule MS/MS spectrum pairs were compared. This comparative analysis revealed that the modified cosine similarity outperformed neutral loss matching and the cosine similarity in all cases. The data further indicated that the performance of MS/MS spectrum alignment depends on the location and type of the modification, as well as the chemical compound class of fragmented molecules.

CFM-ID 4.0 - a web server for accurate MS-based metabolite identification

The CFM-ID 4.0 web server (https://cfmid.wishartlab.com) is an online tool for predicting, annotating and interpreting tandem mass (MS/MS) spectra of small molecules. It is specifically designed to assist researchers pursuing studies in metabolomics, exposomics and analytical chemistry. More specifically, CFM-ID 4.0 supports the: 1) prediction of electrospray ionization quadrupole time-of-flight tandem mass spectra (ESI-QTOF-MS/MS) for small molecules over multiple collision energies (10 eV, 20 eV, and 40 eV); 2) annotation of ESI-QTOF-MS/MS spectra given the structure of the compound; and 3) identification of a small molecule that generated a given ESI-QTOF-MS/MS spectrum at one or more collision energies. The CFM-ID 4.0 web server makes use of a substantially improved MS fragmentation algorithm, a much larger database of experimental and in silico predicted MS/MS spectra and improved scoring methods to offer more accurate MS/MS spectral prediction and MS/MS-based compound identification. Compared to earlier versions of CFM-ID, this new version has an MS/MS spectral prediction performance that is ∼22% better and a compound identification accuracy that is ∼35% better on a standard (CASMI 2016) testing dataset. CFM-ID 4.0 also features a neutral loss function that allows users to identify similar or substituent compounds where no match can be found using CFM-ID’s regular MS/MS-to-compound identification utility. Finally, the CFM-ID 4.0 web server now offers a much more refined user interface that is easier to use, supports molecular formula identification (from MS/MS data), provides more interactively viewable data (including proposed fragment ion structures) and displays MS mirror plots for comparing predicted with observed MS/MS spectra. These improvements should make CFM-ID 4.0 much more useful to the community and should make small molecule identification much easier, faster, and more accurate.