MS2Tox is a machine learning tool for predicting ecotoxicity of unidentified chemicals in water by nontarget LC-HRMS. Calculated values are in unit log-mM (milli-molar in logaritmic scale)
For installing R-package:
remotes::install_github("kruvelab/MS2Tox", INSTALL_opts="--no-multiarch")
library(MS2Tox)
Before installing MS2Tox package Rdisop is needed to be installed. To install this package, start R (version "4.1") and enter:
if (!require("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("Rdisop")
PS! Calculations are currently limited to [M+H]+, [M+Na]+, [M+K]+, [M-H]-, [M+Cl]- adducts
Before using function FishLC50Prediction(folderwithSIRIUSfiles, LC50mode) fingerprints must be calculated with SIRIUS+CSI:FingerID application. Additionally LC50 can be calculated from structure using function LC50fromSMILES(compoundslistwithSMILES, LC50mode). In this case input table needs two columns with names SMILES and exactMass. PS! For code to work, column names should be written exactly like that.
PS! If SIRIUS version 5 and higher is used, function UnZip_SIRIUS5(folderwithSIRIUSfiles) must be used before toxicity calculations.
Folder "ModelsData" contains train, test and validation sets with calculated results.
Folder "Example data" contains SIRIUS+CSI:FingerID calculated fingerprints for testing if the code works.
Step by step instructions from MS data to LC50 values
Before calculating toxicity values, SIRIUS+CSI:FingerID must be used for calculations of fingerprints. Fingerprint calculations with SIRIUS+CSI:FingerID can be done all at once and in this case all the calculated fingerprints will be in this one folder. This folder is used as path to directory in the input for function FishLC50Prediction(). For too many peaks, the SIRIUS+CSI:FingerID calculation may crash. In this case we suggest loading the peaks for calculations in the batch of 200 at a time. Fingerprints can be calculated using SIRIUS GUI, command line or patRoon. More detailed instructions for SIRIUS file can be found from (https://boecker-lab.github.io/docs.sirius.github.io/io/). Here short desription of workflow used by the Kruvelab group while conducting the experiments is given. For this workflow SIRIUS (4.9.5) GUI was used.
For fingerprint calculation it is important that correct input files are dragged into SIRIUS+CSI:FingerID GUI. Example of the .ms file is shown in Figure 1. In order for the toxicity prediction function to give out correct table, SIRIUS+CSI:FingerID input file should start with unique ID. Function uses this ID and later subfolder number to differentiate one compound from another. The first row in the file must start with >compound and then after space can contain either name of the compound or unique id for the analysed peak/spectrum. >parentmass and >ionization are not mandatory but helpful for better identification. Row >formula can be added if known. Information from MS2 spectra are added after row >collision collision value. As SIRIUS+CSI:FingerID adds all the different MS2 spectra together before calculations it is not extremely important to add exactly correct collision energy values or even separate the spectra (RAMP can be used as well). In this case instead of >collision you can use >ms2 and add the MS2 information same way below. For last part after >ms1peaks (or >ms1) isotope pattern from MS1 spectrum can be added; however, this is not mandatory for calculations but can improve the formula prediction.
Figure 1. Example input file for SIRIUS+CSI:FingerID
For calculations .ms files can be dragged into the GUI. Before that new project directory must be established by clicking to icon New (see Fig 2A). This path to the directory will be used in final R function for calculations. After dragging .ms files into the application all the information about each file is listed on the left side in the application (Fig 2B). By clicking on Compute All button (Fig 2C) compute window (Fig 2D) opens. In this window parameters for calculations can be chosen. In this work for Orbitrap measurements mass deviation 5 ppm was chosen. Calculations may take time from few minutes to hours depending on the number of .ms files and molecular mass of the compounds. After calculations a subfolder per .ms file has been created and it contains calculated fingerprints, scores, fragmentation trees, etc. All output folders start with order number that is shown in the final table in column named “foldernumber”. See Fig 3.
Figure 2. SIRIUS+CSI:FingerID application with chosen parameters for calculations
Table 1. Selected elements and their maximum amount for SIRIUS+CSI:FingerID calculations
| Element | Amount | Element | Amount | Element | Amount | Element | Amount | |||
|---|---|---|---|---|---|---|---|---|---|---|
| H | Inf | C | Inf | N | Inf | O | Inf | |||
| P | 8 | B | 11 | Si | 9 | S | 12 | |||
| Cl | 18 | Se | 2 | Br | 10 | F | Inf | |||
| I | 6 | K | 1 | Na | 1 | As | 2 |
Figure 3. Folder after fingerprint calculations that directory must be put into the function
R package MS2Tox can be then used for toxicity value predictions. For installing R-package:
remotes::install_github("kruvelab/MS2Tox")
library(MS2Tox)
Before installing MS2Tox package Rdisop is needed to be installed. To install this package, start R (version "4.1") and enter:
if (!require("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("Rdisop")
FishLC50Prediction() function takes in the path to the directory containing all the calculated fingerprints and information about the LC50 mode. Right now there are two possible LC50 modes for which predictions can be made: static and flow through. If LC50mode is not chosen, calculations will be done automatically in static mode. Directory that will be input for prediction function is shown in Fig 3. It returns a summary table with calculated toxicities for all of the peaks for which fingerprints were calculated by SIRIUS+CSI:FingerID. The calculation is based on the fingerprints associated with the first ranked molecular formula. SIRIUS+CSI:FingerID predicted formula, prediced ion and folder number are also displayed in the summary table. In order for the code to work properly input files for SIRIUS+CSI:FingerID must contain the compounds id/name in the first part of the file name (e.g AU231458_restofthename.ms). This first part of the name is used as id column in the summary table and can be used to link a specific peak with a specific toxicity value.
If data sets are very big, in SIRIUS+CSI:FingerID calculations results could be saved so that one directory contains up to 1000 subfolders, otherwise calculation with FishLC50Prediction() might exceed the memory limits of R.
The main function links to other functions, which gather the fingerprint data (FingerPrintTable()), process it, filter out only first rank chemical formulas (FpTableForPredictions()) based on SiriusScore (SiriusScoreRank1()) and calculate the toxicity value. PS! If for spectrum, two formulas are calculated with same SiriusScore, both of them are considered as first rank and in this case two different toxicity values may be seen for one peak/id. In final table, possible chemical formulas are given as well.
In this work SIRIUS+CSI:FingerID version 4.9.5 was used so the prediction method is limited to these fingerprints. If newer or older versions do not calculate some of the fingerprints used here, the FishLC50Prediction() might not work . In this case please notify the author for retraining the model.
PS! Versions up to 5.5.7 have been tested and work! If SIRIUS 5 is used, fingerprints and score must be unzipped beforehand. For that run this before LC50 predictions:
folderwithSIRIUSfiles <- "C:/Desktop/Folder" #Add your directory with SIRIUS calculated fingerprints
UnZip_SIRIUS5(folderwithSIRIUSfiles)
Code example for calculating LC50 values:
setwd(C:/Desktop/MS2Tox) #Directory must be set to downloaded MS2Tox folder as trained model is located there
folderwithSIRIUSfiles <- "C:/Desktop/Folder" #Add your directory with SIRIUS calculated fingerprints
chosen_mode <- “static” #or "flow"
results <- FishLC50Prediction(folderwithSIRIUSfiles, chosen_mode)
Addition to SIRIUS fingerprints it is also possible to calculate LC50 values using SMILES as an input. For that use function SMILESFingerprints(compoundslistwithSMILES). PS! Input needs to be a table containing columns named "SMILES" and "exactMass". Use exactly this format for naming.