mtbinfo.github.io

msinfo manuals

View the Project on GitHub systemsomicslab/mtbinfo.github.io

MS-FINDER tutorial

Last edited in Dec. 6, 2018

Abstract

The purpose of metabolomics is to perform the ‘comprehensive’ analysis for small biomolecules of living organisms. Gas chromatography coupled with electron ionization mass spectrometer (GC/MS) and liquid chromatography coupled with electrospray ionization- (ESI-) tandem mass spectrometer (LC/MS/MS) are the preferred tools for untargeted metabolomics. Currently, the main bottleneck of GC/MS- and LC/MS/MS based untargeted analysis is compound identification due to the limitation of EI-MS and MS/MS records of authentic standard.
    MS-FINDER was launched as a universal program for compound ‘annotation’ that supports EI-MS (GC/MS) and MS/MS spectral mining. First, MS-FINDER aims to provide solutions for 1) formula predictions, 2) fragment annotations, and 3) structure elucidations by means of unknown spectra. In addition, the program can annotate your unknowns by the public spectral databases such as MassBank, LipidBlast, and GNPS. MS-FINDER has been developed as the collaborative work between Prof. Masanori Arita team (RIKEN, Reifycs Inc.) and Prof. Oliver Fiehn team (UC Davis) supported by the JST/NSF SICORP ‘Metabolomics for the low carbon society’ project.

Hiroshi Tsugawa
RIKEN Center for Sustainable Resource Science
hiroshi.tsugawa@riken.jp


alt MS-FINDER screenshot

 

Section 1

Software environments

 

Section 2

Required programs

 

Section 3

Acceptable ASCII formats

This program accepts two file extensions, i.e. MSP or MAT formatted by the following explanations. Unknown queries should be separately stored in the ASCII file: the MSP or MAP file ‘CANNOT’ store multi compound records in the single file.

The format of MSP basically follows the NIST MS search manual.
Link: http://www.nist.gov/srd/upload/NIST1a11Ver2-0Man.pdf

 

Section 3-1

MSP format for MS/MS

Required fields

MSP example
alt  

Section 3-2

MSP format for EI-MS

Required fields

The fields are the minimum requirement for searching spectral databases. In the case that you want to perform formula predictions and structure elucidations in EI-MS data, two files ‘PRECURSORMZ:’ and ‘PRECURSORTYPE:’ must be required.

 

Section 3-3

MAT format

The MAT format was defined as the improved version of MSP in the MS-FIDNER program to store both MS1 and MS/MS spectra in the same file. The survey scan MS data should be required to calculate ‘isotopic ion score’ for formula predictions. Importantly, for EI-MS spectra, put your spectra into both MS1- and MS2 fields for the calculation of isotopic ratio and fragment ion similarities, respectively.

Required fields

✓ Three fields including MSTYPE, Num Peaks, and m/z intensity pair should be SERIALLY stored.

✓ If you type ‘MSTYPE: MS1’, the spectrum written from next field should be recognized as the survey scan MS (MS1). If you type ‘MSTYPE: MS2’, next spectrum should be recognized as the MS/MS spectrum.

✓ EI-MS spectra must be stored in both MS1 and MS2, which is the requirement of MS-FINDER.

✓ Both field (MSTYPE: MS1 and MSTYPE: MS2) is not necessary for this program, i.e. the users can import the ASCII file as only MS1 spectrum or as only MS/MS spectrum record.

✓ Users may prepare the MAT or MSP files without any spectrum record. In such case, the formula prediction will be performed by means of mass accuracy and database criteria.

✓ If you want to perform the MS/MS peak annotation with the known structure, prepare two fields including FORMULA and SMILES. The formula and SMILES of the neutralized structure should be made.

 
MAT example
alt

Section 3-4 

Adduct ion format: [M+Na]+, [M+2H]2+, [M-2H2O+H]+, [2M+FA-H]-, etc.

  1. The parentheses ‘[’ and ‘]’ must be used to bracket the ion information.
  2. The char + and - must be required after ‘]’ and the number must be written before + or -.
  3. When you want to define the organic formula like C6H12O5, you have to write it without any replicate elements or parentheses. For example, the descriptions like [M+C2H5COOH-H]- or [M+H+(CH3)3SiOH]+ are not accepted.
  4. The beginning figure of organic formula like ‘2’H2O is recognized as the H2O × 2. Again, never use 2(H2O) for that.
  5. Sequential equations are acceptable: [2M+H-C6H12O5+Na]2+ (very apt.)
  6. Radical ion can be described by ‘.(dot)’ after + or – like [M]+. And [M-CH3]+. as adduct format.
  7. MS-FINDER accepts some abbreviations or common organic formulas for adduct types as follows.
    For Acetonitrile: ACN, CH3CN
    For Methanol: CH3OH
    For Isopropanol: IsoProp, C3H7OH
    For Dimethyl sulfoxide: DMSO
    For Formic acid: FA, HCOOH
    For Acetic acid: Hac, CH3COOH
    For Trifluoroacetic acid: TFA, CFCOOH

 

Section 3-5 

User defined structure database format

MS-FINDER supports the structure elucidations from the candidates that users provide. The following format file should be prepared as tab-delimited text file. The identifiers of InChIKey, short InChIKey, and database ID are not required, but the values must be filled by some mimic values. If you set the correct InChIKeys, the database scoring is also utilized for ranking structures. The files of exact mass, formula, and SMILES must be prepared.

alt
 

Section 3-6 

Specific field to fix the formula element count

MS-FINDER recently accepts the filtering field for molecular formula prediction. For example, if you can perform fully labeled stable isotope experiment, the element count for CHNOS can be determined by checking the mass shift between non-labeled and labeled sample data. The fields described in the below figure can be written in MSP and MAT files.

alt
 

Section 4

Import queries

There are three ways to import unknown queries.

Section 4-1

From a folder which includes MSP or MAT format files

  1. File -> import
  2. Select a folder containing MSP or MAT files.

alt

 

Section 4-2

From the graphical user interface of the MS-FINDER program

  1. File -> Create a query
  2. Fill in the form to make a query.

Required files

alt
 

Section 4-3

From the MS-DIAL program

The MS-DIAL program which has been reported as software for data processing of LC/MS/MS can call the MS-FINDER program directly. On the first time when you call the MS-FINDER program at MS-DIAL, please select the file path of MS-FINDER via GUI.

alt

Section 5

Parameter setting

Section 5-1

Method tab

alt

MS-FINDER provides two options for compound annotation: one is by spectral databases, and the other is by formula- and structure finder programs using in silico fragmenter. You can simultaneously check both ‘spectral database search’ and ‘formula prediction and structure elucidation by in silico fragmenter’ options, and the result of spectral database search has priority for ranking structures.
    The ‘internal experimental library’ is stored in ‘EIMS-DBs-vs*.egm’ and ‘MSMS-DBs-vs*.etm’ of Resources folder as NIST MSP format. If ‘Formula finder > TMS-MeOX derivative compound’ is checked, EIMS database will be used; otherwise, MSMS database is used. The ‘in silico library for lipids (LipidBlast)’ is stored in ‘MSDIAL-LipidDBs-vs*.lbm’ of Resources folder. Select the appropriate solvent condition for searching your unknowns. The user-defined spectral database must be formatted by NIST MSP.
    If ‘Precursor oriented spectral search’ is checked, the structure candidates will be filtered out by the precursor m/z of spectral records in combination with MS1 tolerance value; otherwise, all of spectral records will be used. Uncheck this option if you want to search EI spectral databases.
 

Section 5-2

Mass spectrum tab

alt

Mass tolerance (MS1): the mass tolerance to generate formula candidates.

Mass tolerance (MS/MS): the mass tolerance for matching experimental- and reference fragments.

Relative abundance cut off: The product ions more than this parameter on the basis of base peak ion are utilized for the product ion matching.

* For EI-MS spectra, set the same tolerance into MS1 and MS2.

 

Section 5-3

Formula finder tab

alt

Formula calculation setting: You can set the parameters for formula calculation.

LEWIS and SENIOR check: to generate formula candidates that match the valence rules of formula elements. The valences of hetero atoms, i.e. N, O, S, and P are currently set to 3, 2, 6, and 5, respectively.

Isotopic ratio tolerance: to calculate the isotopic score. The tolerance should be utilized as the sigma value for the Gaussian scoring as described in the MS-FINDER paper.

Element ratio check: to generate formula candidates that satisfy every element ratios (ex. H/C ratio should be between 0 and 3.33 for ‘Common range (99.7%)’ restriction. ) as described in the MS-FINDER paper.

Element probability check: to generate formula candidates that satisfy the heuristic rules as described in the Seven Golden Rules paper. For example, if a formula candidate contains the following element counts, i.e. NOPS all > 1, the element counts of N, O, P, and S should be less than 9, 19, 3, and 2, respectively.

Element selection: to generate formula candidates that just contain the elements selected by the users. Check ‘TMS-MEOX derivative compound’ if you want to annotate EI-MS spectra.

Result cut off: formula candidates ranked by the MS-FINDER program will be reported within up to this number.

Section 5-4

Structure finder tab

alt

Here is the parameter setting for in silico fragmenter.

Tree depth: the limitation of in silico cleavages, i.e. if the user sets ‘2’, the MS-FINDER program generates the fragments until product ions of a product ion.

The current MS-FINDER program can utilize the fragment ion library for EI-MS spectral mining which is stored in ‘*.eif’ (recommended to use).

Result cut off: structure candidates ranked by the MS-FINDER program will be reported within up to this number.

XLogP based RT prediction and cut off setting for structure elucidation
Recently, MS-FINDER provides a simple RT prediction function using XLogP calculated by CDK. If you prepare the tab-delimited text format file containing (first column) metabolite name, (second column) retention time (min), and (third column) SMILES code as described below, the predicted retention time is calculated for searching structure candidates.

Retention time setting for spectral searching
This function is for spectral searching. If your MSP or MAT files contain the retention time or retention index information, you can use the RT or RI filtering by means of this checkbox.

alt
 

Section 5-5

Data source tab

alt

Local databases: currently, total 14 metabolome databases are prepared in the MS-FINDER program which is stored in ‘*.esd’. The local databases selected by users will be used to retrieve the structure data. Please see ‘user defined database format’ section for searching your own structure candidates.

MINEs (Metabolic In silico Network Expansions) and PubChem online settings:


Section 6

Compound annotation by in silico fragmenter

The general workflow of MS-FINDER is described here. (Batch analysis is shown below)

1. The formula prediction is executed by double click at the title name of ‘File navigator’. The detail of formula calculations is described in our paper. *The point of formula prediction is to correctly select the precursor type (adduct ion) and to set parameters for picking up the correct formulas.

alt

2. The structure finder will be executed by right click at the formula result table followed by clicking ‘Search the structure’.

alt

Section 7

Compound annotation by searching spectral databases

It’s very simple. Double click the unknown record that you want to annotate. The below is the examples for searching EI-MS spectral database and LipidBlast database.
alt

 

Section 8

FSEA: fragment set enrichment analysis

This function recommends the ‘metabolite class’ for unknown MS/MS spectrum. The classification is performed in formula prediction process. Therefore, you can see the result once the formula candidates are generated. To check the result of FSEA (for example, for C20H18O10 as shown below), please click ‘FSEA’ button after the row of formula finder results is activated. Because the correct structure is ‘Kaempferol-3-O-arabinoside’, the second hit of recommended metabolite classes, i.e. Flavonol O-glycosides, is the correct metabolite class for this query: this is actually the current ‘accuracy’ of metabolite class recommendation, i.e. several candidates will be suggested with the information of estimated p-values. The p-value estimation can be executed by three different ‘decoy’ criteria.

alt

alt

alt

The challenging issue to use over-representation analysis (ORA) for metabolite class recommendation is how to determine the ‘non-significant’ peaks in the MS/MS spectrum set. Currently, there are three options to define the non-significant peaks. You can change the setting at ‘Setting’ -> ‘FSEA parameter setting’. If you choose ‘Low abundance ions’ with 5% relative abundance cut off, the ions having less than 5% relative abundances are recognized as ‘non-significant’ peaks. The result of FSEA is shown below: 29 metabolite classes are recommended with less than 1% estimated p-value.

alt

Moreover, if you choose the ‘reverse spectrum’, all of ‘significant’ peaks in the reverse spectrum is defined as ‘non-significant’ peaks. The reversed spectrum is created by mirroring the m/z nominal places with a 5-mass shift. For example, the m/z sequence [101.023, 305.035, 421.098, 634.201, 754.235] is changed to [106.023, 226.035, 439.098, 555.201, 759.235]. The FSEA result is shown below: 4 classes are recommended with 1% p-value cut off.

alt

As the third option (default), the non-significant peak is automatically calculated by undetected ontologies in the FSEA set (n = 459), i.e. the current knowledge-space of fragment ontologies. The FSEA result is shown below: 29 candidates are recommended.
    Generally, the option of ‘reverse spectrum’ is the most conservative method to reduce false positive recommendation for metabolite classes although the false negatives are increased in this option. In order to improve the accuracy of metabolite class recommendation, the program will be further updated in future.

alt

 

Section 9

Compound annotation (batch analysis)

1. Analysis -> Compound annotation (batch job)

A) If you want to perform the batch job for both formula predictions and structure elucidations, please check ‘Both processes. Also, add the number in ‘Top N hits’ textbox where the formula candidates generated by the formula finder program are supposed to be searched.

alt

B) If you want to perform the batch job for formula predictions, please just check ‘molecular formula finder’.

alt

C) If you want to perform the structure finder program, please just check ‘structure finder’. Here, the formula candidates checked by the users are supposed to be examined. Also, if the formula finder is not executed before this analysis, the queries will be passed.

alt
 

Section 10

Peak assignment (single)

The MS-FINDER program can be used as the peak assignment tool to assign substructures in the MS/MS spectrum from user-defined structure.

  1. Analysis -> Peak assignment (single)
  2. Select the query file that you want to analyze.
  3. Add both formula and SMILES into the textboxes as the neutralized form.
  4. The result will be generated as shown below.

alt
 

Section 11

Peak assignment (batch job)

Analysis-> Peak assignment (batch job)

alt

 

Section 12

Molecular spectrum networking

The integrated networking method using MS/MS similarity, ontology similarity, and formula based bioreaction program can be executed. The requirement to fully use this function is that

  1. Import the queries where you obtained with the same analytical condition.
  2. Execute ‘batch job’ from ‘Analysis’ -> ‘Compound annotation (batch job)’.
  3. You may modify the biotransformation rule at the ‘Resource’ folder -> ‘Biotransformation-VS1.fbt’ (this file can be opened by a text editor like notepad++).

To execute the network analysis, please follow as ‘Analysis’ -> ‘Molecular spectrum networking’. Each query will be linked with others by using MS/MS similarity (red), ontology similarity (blue), and formula oriented bioreaction (yellow). The top hit result is used for the network analysis. Due to the limitation of cytoscape.js (to be used in a local browser), not all of nodes and edges are shown in the network. Therefore, you can export all of nodes- and edges by ‘Export’ -> ‘Export nodes and edges for molecular spectrum network’ as a tab delimited file to be used in locally available cytoscape program (http://www.cytoscape.org/cy3.html).

alt

alt
 

Section 13

Mouse function

alt

 

Section 14

Export

The result of formula and structure finders can be exported from this option. Currently, the top 10 candidates will be automatically exported. However, you can also check the details of result as the ASCII file. It means that the MS-FINDER program is supposed to generate the FGT file containing formula results in the same directory as the project folder. Also, the program generates the respective folder containing the SFD file (per formula candidate) which stores the result of structure finder program.

alt
 

Section 15

Help

You can check the version of MS-FINDER.

alt