Identifying the Unknowns by Aligning Fragmentation Trees

Self Cite

840

829

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.mass spectrometry | small compound identification | metabolomics | bioinformatics | machine learning M etabolites, small molecules that are involved in cellular reactions, can provide detailed information about cellular state. Untargeted metabolomic studies may use NMR or MS technologies, but liquid chromatography followed by MS (LC/MS) can detect the highest number of metabolites from minimal amounts of sample (1, 2). Untargeted metabolomics comprehensively compares the mass spectral intensities of metabolite signals (peaks) between two or more samples (3, 4). Advances in MS instrumentation allow us to simultaneously detect thousands of metabolites in a biological sample. Identification of these compounds relies on tandem MS (MS/MS) data, produced by fragmenting the compound and recording the masses of the fragments. Structural elucidation remains a challenging problem, in particular for compounds that cannot be found in any spectral library (1): In total, all available spectral MS/MS libraries of pure chemical standards cover fewer than 20,000 compounds (5). Growth of spectral libraries is limited by the unavailability of pure reference standards for many compounds.In contrast, molecular structure databases such as PubChem (6) and ChemSpider (7) contain millions of compounds, with PubChem alone having surpassed 50 million entries. Searching in molecular structure databases using MS/MS data is therefore considered a powerful tool for assisting an expert in the elucidation of a compound. This problem is considerably harder than the fundamental analysis step in the shotgun proteomics workflow, namely, searching peptide MS/MS data in a peptide sequence database (8): Unlike proteins and peptides, metabolites show a large structural variability and, consequently, also large variations in MS/MS fragmentation. Computational approaches for interpreting and predicting MS/MS data of small molecules date back to the 1960s (9): Due to the unavailability of molecular structure databases at that time, structure libraries were combinatorially generated and then "searched" using the experimental MS/MS data. "Modern" methods for this question have been developed since mid-2000. Particular progress has been made for restricted metabolite cl...

Section: Resultsmentioning

confidence: 99%

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

Dührkop

Shen

Meusel

et al. 2015

Self Cite

840

829

“…for genome mining for new RiPPs because of the direct link between final structure and the gene-encoded precursor peptide, the method in principle can also be used for structural analysis of other compounds such as nonribosomal peptides, and is complementary to other tandem MS-based methods that aim to ultimately determine structures in high throughput, including blind search (38), network analysis (39)(40)(41), and fragmentation tree construction (42,43). Moreover, HSEE greatly facilitates investigations of RiPP biosynthetic mechanisms, as illustrated by several examples in this study.…”

Section: Resultsmentioning

confidence: 99%

Structural investigation of ribosomally synthesized natural products by hypothetical structure enumeration and evaluation using tandem MS

Zhang

Ortega

Shi

et al. 2014

Ribosomally synthesized and posttranslationally modified peptides (RiPPs) are a growing class of natural products that are found in all domains of life. These compounds possess vast structural diversity and have a wide range of biological activities, promising a fertile ground for exploring novel natural products. One challenging aspect of RiPP research is the difficulty of structure determination due to their architectural complexity. We here describe a method for automated structural characterization of RiPPs by tandem mass spectrometry. This method is based on the combined analysis of multiple mass spectra and evaluation of a collection of hypothetical structures predicted based on the biosynthetic gene cluster and molecular weight. We show that this method is effective in structural characterization of complex RiPPs, including lanthipeptides, glycopeptides, and azole-containing peptides. Using this method, we have determined the structure of a previously structurally uncharacterized lanthipeptide, prochlorosin 1.2, and investigated the order of the posttranslational modifications in three biosynthetic systems.dehydration | genome mining | lantibiotics | directionality R ibosomally synthesized and posttranslationally modified peptides (RiPPs) are a major class of natural products as revealed by the genome-sequencing efforts of the past decade (1). RiPPs are biosynthesized from genetically encoded and ribosomally produced precursor peptides, which typically consist of a core peptide that is transformed to the final product and an N-terminal extension called the leader peptide that is usually important for recognition by the posttranslational modification (PTM) enzymes (1). Because of the highly diverse PTMs, these compounds possess vast structural diversity and have a wide range of biological activities, thus representing a fertile ground for exploration. Furthermore, the ribosomal origin of RiPPs makes them particularly well suited for genome mining efforts. By using genome mining to explicitly avoid species harboring biosynthetic gene clusters identical to those that produce known compounds, a combination of strain prioritization and mass spectrometry (MS)-based analysis offers a new route to discovering natural products that can overcome the burden of rediscovery that has increasingly hampered discovery efforts (2, 3). One challenging aspect of high-throughput genome mining for new natural products is the difficulty to determine their molecular structures in high throughput. We present here a method that allows automated RiPP structure elucidation.In contrast to nonribosomal peptides that have an average molecular weight of less than 1,000 Da, as documented in the NORINE database (4), RiPPs in many cases have molecular weights larger than 2,500 Da. Molecules of this size are difficult to rapidly analyze by NMR spectroscopy, rendering MS the most convenient tool for RiPP structural characterization. Even when the precursor peptide sequences are known and the types of PTMs can be predicted based on the sequen...

“…Several recent studies have highlighted the importance of considering neutral losses (NLs) for aligning spectra and constructing similarity clusters overlapping with compound familial groupings. NL analysis has been especially well-implemented in the context of MS/MS fragmentation tree studies (53,54). We therefore combined these two types of information into a bidimensional clustering method.…”

Section: A Biclustering Classification Of the Idms/ms Landscape Highlmentioning

confidence: 99%

Navigating natural variation in herbivory-induced secondary metabolism in coyote tobacco populations using MS/MS structural analysis

Baldwin

Gaquerel

2015

Natural variation can be extremely useful in unraveling the determinants of phenotypic trait evolution but has rarely been analyzed with unbiased metabolic profiling to understand how its effects are organized at the level of biochemical pathways. Native populations of Nicotiana attenuata, a wild tobacco species, have been shown to be highly genetically diverse for traits important for their interactions with insects. To resolve the chemodiversity existing in these populations, we developed a metabolomics and computational pipeline to annotate leaf metabolic responses to Manduca sexta herbivory. We selected seeds from 43 accessions of different populations from the southwestern United Statesincluding the well-characterized Utah 30th generation inbred accession-and grew 183 plants in the glasshouse for standardized herbivory elicitation. Metabolic profiles were generated from elicited leaves of each plant using a high-throughput ultra HPLC (UHPLC)-quadrupole TOFMS (qTOFMS) method, processed to systematically infer covariation patterns among biochemically related metabolites, as well as unknown ones, and finally assembled to map natural variation. Navigating this map revealed metabolic branch-specific variations that surprisingly only partly overlapped with jasmonate accumulation polymorphisms and deviated from canonical jasmonate signaling. Fragmentation analysis via indiscriminant tandem mass spectrometry (idMS/MS) was conducted with 10 accessions that spanned a large proportion of the variance found in the complete accession dataset, and compound spectra were computationally assembled into spectral similarity networks. The biological information captured by this networking approach facilitates the mining of the mass spectral data of unknowns with high natural variation, as demonstrated by the annotation of a strongly herbivory-inducible phenolic derivative, and can guide pathway analysis.plant-insect interactions | metabolomics | mass spectrometry | natural variation E lucidating the structure of metabolites underlying complex traits and the factors that maintain their variation in natural populations are important challenges in plant ecological studies (1). Many studies have notably shown that stress-responsive pathways that produce secondary metabolites are sporadically found across different plant taxa with extensive diversification (2). This important diversification suggests that particular metabolic systems have been recruited through natural selection when the set of compounds that they produce address specific ecological needs. Interactions with insects are important selection pressures that have sculpted plant metabolism, and many plant metabolites protect against herbivore attack and physical damage (3-5). The timely production of particular secondary metabolites in response to insect attack benefits plants by decreasing the costs of constitutive metabolite production. Trade-offs between defense metabolite productions and the intrinsic growth-related functions of central metabolic pathways likely prov...