-
Notifications
You must be signed in to change notification settings - Fork 25
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Find all isotope peaks in a spectrum #185
Comments
@jorainer I have highly accurate isotope distributions and background proportions if you need them and so code to simulate the isotope distribution given the sequence |
@andreavicini is currently calculating distributions based on all chemical formulas of metabolites from HMDB (human metabolome database). On what did you calculate that? Our approach is currently simpler than isotope distribution simulation - we're essentially looking for peaks with a difference in m/z that matches the expected difference for an isotope (e.g. C12, C13) allowing a user-defined ppm and checking that the intensity is lower than a certain threshold. |
So currently, I'm using, for example 12C and 13C that their masses are It looks like your use cases is slightly different, but thought I'd share in case its useful to discuss |
OK, if I get you correctly, in your case the sequence (=chemical formula) and the charge is known beforehand. That's definitely also a good use case. Is that somewhat similar to what My use case at present is a completely unsupervised one, given that I have a spectrum, identify groups of peaks that could represent isotope peaks of a (yet unknown) compound. |
Yep, exactly, mine is more simialr to envipat, just it returns a spectra object so its easier to use. Though, would also be cool in your unsupervised approach to be able to identify a glyco or phospho group (because that is unknown for us). Thanks for clarification - look forward to the development! |
Similar issue: rformassspectrometry/MetaboCoreUtils#10
|
Thanks @sgibb ! I completely forgot about that one! |
Maybe you can get inspired here: https://github.com/RECETOX/recetox-xMSannotator/blob/main/xmsannotator/R/compute_isotopes.R The rdkit chem library gives you the pattern, so with some spectral matching you could maybe identify those peaks. |
Given a spectrum, find all sets of peaks that could represent isotope groups (e.g. C12, C13 peaks). This functionality could then be used e.g. in a
filterIsotopes
function or another function to extract just isotope peaks from aSpectra
(e.g. to pass it to functions to predict the formula based on the isotope pattern).The text was updated successfully, but these errors were encountered: