Peak fitting in 2D 1H-13C HSQC NMR spectra for metabolomic studies

A modified Lorentzian distribution function is used to model peaks in two-dimensional (2D) H-1-C-13 heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra. The model fit is used to determine accurate chemical shifts from genuine signals in complex metabolite mixtures such as blood. The algorithm can be used to extract features from a set of spectra from different samples for exploratory metabolomics. First a reference spectrum is created in which the peak intensities are given by the median value over all samples at each point in the 2D spectra so that H-1-C-13 correlations in any spectra are accounted for. The mathematical model provides a footprint for each peak in the reference spectrum, which can be used to bin the H-1-C-13 correlations in each HSQC spectrum. The binned intensities are then used as variables in multivariate analyses and those found to be discriminatory are rapidly identified by cross referencing the chemical shifts of the bins with a database of C-13 and H-1 chemical shift correlations from known metabolites.