RTExtract: time-series NMR spectra quantification based on 3D surface ridge tracking

Abstract Motivation Time-series nuclear magnetic resonance (NMR) has advanced our knowledge about metabolic dynamics. Before analyzing compounds through modeling or statistical methods, chemical features need to be tracked and quantified. However, because of peak overlap and peak shifting, the available protocols are time consuming at best or even impossible for some regions in NMR spectra. Results We introduce Ridge Tracking-based Extract (RTExtract), a computer vision-based algorithm, to quantify time-series NMR spectra. The NMR spectra of multiple time points were formulated as a 3D surface. Candidate points were first filtered using local curvature and optima, then connected into ridges by a greedy algorithm. Interactive steps were implemented to refine results. Among 173 simulated ridges, 115 can be tracked (RMSD < 0.001). For reproducing previous results, RTExtract took less than 2 h instead of ∼48 h, and two instead of seven parameters need tuning. Multiple regions with overlapping and changing chemical shifts are accurately tracked. Availability and implementation Source code is freely available within Metabolomics toolbox GitHub repository (https://github.com/artedison/Edison_Lab_Shared_Metabolomics_UGA/tree/master/metabolomics_toolbox/code/ridge_tracking) and is implemented in MATLAB and R. Supplementary information Supplementary data are available at Bioinformatics online.

[1]  M. D. Joshi,et al.  Complete measurement of the pKa values of the carboxyl and imidazole groups in Bacillus circulans xylanase , 1997, Protein science : a publication of the Protein Society.

[2]  Timothy M. D. Ebbels,et al.  Modelling the acid/base 1H NMR chemical shift limits of metabolites in human urine , 2016, Metabolomics.

[4]  A. Jenne,et al.  In-Vivo NMR Spectroscopy: A Powerful and Complimentary Tool for Understanding Environmental Toxicity , 2018, Metabolites.

[5]  Suchendra M. Bhandarkar,et al.  Three-Dimensional Object Recognition from Range Images , 1992, Computer Science Workbench.

[6]  Maurice Berk,et al.  Modelling short time series in metabolomics: a functional data analysis approach. , 2011, Advances in experimental medicine and biology.

[7]  Johannes Schindelin,et al.  TrackMate: An open and extensible platform for single-particle tracking. , 2017, Methods.

[8]  Ralf J. O. Torgrip,et al.  Proof of principle of a generalized fuzzy Hough transform approach to peak alignment of one-dimensional 1H NMR data , 2007, Analytical and bioanalytical chemistry.

[9]  Ramesh C. Jain,et al.  Invariant surface characteristics for 3D object recognition in range images , 1985, Comput. Vis. Graph. Image Process..

[10]  I. Lindberg,et al.  Structural studies of a neuropeptide precursor protein with an RGD proteolytic site. , 2001, Biochemistry.

[11]  U. Sauer,et al.  Advancing metabolic models with kinetic information. , 2014, Current opinion in biotechnology.

[12]  R. Haralick,et al.  The Topographic Primal Sketch , 1983 .

[13]  K. Laukens,et al.  Getting Your Peaks in Line: A Review of Alignment Methods for NMR Spectral Data , 2013, Metabolites.

[14]  M. Simpson,et al.  Reducing impacts of organism variability in metabolomics via time trajectory in vivo NMR , 2018, Magnetic resonance in chemistry : MRC.

[15]  S. Withers,et al.  The pKa values of the catalytic residues in the retaining glycoside hydrolase T26H mutant of T4 lysozyme , 2018, Protein science : a publication of the Protein Society.

[16]  M. Fey,et al.  NMR flow tube for online NMR reaction monitoring. , 2014, Analytical chemistry.

[17]  U. Günther,et al.  Metabolic plasticity in CLL: adaptation to the hypoxic niche , 2015, Leukemia.

[18]  G. Hägele,et al.  1H/31P NMR pH indicator series to eliminate the glass electrode in NMR spectroscopic pKa determinations , 2004 .

[19]  Christian Ludwig,et al.  MetaboLab - advanced NMR data processing and analysis for metabolomics , 2011, BMC Bioinformatics.

[20]  Michael E. Taylor,et al.  Differential Geometry I , 1994 .

[21]  Ramesh C. Jain,et al.  Segmentation through Variable-Order Surface Fitting , 1988, IEEE Trans. Pattern Anal. Mach. Intell..

[22]  Carsten Steger,et al.  An Unbiased Detector of Curvilinear Structures , 1998, IEEE Trans. Pattern Anal. Mach. Intell..

[23]  Piotr Klukowski,et al.  NMRNet: a deep learning approach to automated peak picking of protein NMR spectra , 2018, Bioinform..

[24]  Jie Hao,et al.  Combining spectral ordering with peak fitting for one-dimensional NMR quantitative metabolomics. , 2013, Analytical chemistry.

[25]  A. Edison,et al.  Conformational Ensembles: The Role of Neuropeptide Structures in Receptor Binding , 1999, The Journal of Neuroscience.

[26]  Piotr Klukowski,et al.  Computer vision-based automated peak picking applied to protein NMR spectra , 2015, Bioinform..

[27]  Manfred Spraul,et al.  Deconvoluting interrelationships between concentrations and chemical shifts in urine provides a powerful analysis tool , 2017, Nature Communications.

[28]  K. Krishnamurthy.  CRAFT (complete reduction to amplitude frequency table) – robust and time‐efficient Bayesian approach for quantitative mixture analysis by NMR , 2013, Magnetic resonance in chemistry : MRC.

[29]  U. Sauer,et al.  Real-time metabolome profiling of the metabolic switch between starvation and growth , 2015, Nature Methods.

[30]  Maria De Iorio,et al.  Bayesian deconvolution and quantification of metabolites in complex 1D NMR spectra using BATMAN , 2014, Nature Protocols.