Skip to main content

Auto MS mass checker

Project description

MScheck for locating target compound masses in mass spectra

PyPI version

build test

MScheck is a python package that hunts for a target compound mass + given ion mass (eg. H+, Na+). MScheck was created to assist with the automated mass spectrum analysis of target compounds synthesised using a high throughput approach.

MSCheck uses the excellent Python library pyOpenMS for reading and handling .mzML mass spectra files. The spectra are stored and handled as a MassSpectrum class object - see spectrum.py.

The AnalyseSpectum class - see analyse.py - use Scipy's signal peak finding algorithms (find_peaks and peak_widths) to find peaks and calculate the full width at half maximum height (FWHM) of the peaks found. Mass spectrum data points are analysed in the area of the peak above the FWHM height by searching for the sum of the parent mass of the target molecule and ion. Different ions can be included in the search - see the example below.

MScheck installation

MScheck relies on rdkit for generating molecule SVG images and for calculating molecular weights
Installing rdkit using conda works best followed by a pip install of MScheck

  1. Create a conda environment with rdkit

    conda create -c conda-forge -n MScheck rdkit

  2. Activate the MScheck conda environment created

    conda activate MScheck

  3. Pip install MScheck

    pip install mscheck

Preparing vendor files

MScheck has been tested on using Agilent LCMS files (.D) as the starting file format. One challenge is to convert vendor file formats into a format with the binary decoded.

  1. Convert .D Agilent folder to a .d MassHunter format using Agilent's ChemStation to MassHunter Translator (B.04.00)
  2. Convert .d format into .mzML format using ProteoWizard's MSConvert tool
  3. Finally - we have an file format that we can use!
  4. In your favourite IDE or Jupyter notebook - as basic example of using MSCheck is provided below:
from mscheck.analyse import AnalyseSpectrum

# Create MS scptrum object and find peaks
test = AnalyseSpectrum("<path to .mzML file>", mode="Positive")

# Analyse test spectrum
test.analyse(compoundsmiles=target_compound,
             ionstoadd=["[H]", "[Na]", "[K]", "[NH4+]"],
             tolerance=1)test.analyse

# Create a .svg report - if you do not give a compound_name
# the ending leaf of the file name will be used
test.create_report(compound_name="Test")

The .svg report will be in a folder called Reports

Example of report output:

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for mscheck, version 0.1.9b0
Filename, size File type Python version Upload date Hashes
Filename, size mscheck-0.1.9b0.tar.gz (10.5 kB) File type Source Python version None Upload date Hashes View

Supported by

AWS AWS Cloud computing Datadog Datadog Monitoring DigiCert DigiCert EV certificate Facebook / Instagram Facebook / Instagram PSF Sponsor Fastly Fastly CDN Google Google Object Storage and Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Salesforce Salesforce PSF Sponsor Sentry Sentry Error logging StatusPage StatusPage Status page