New modes of data partitioning based on PARS peak alignment for improved multivariate biomarker/biopattern detection in 1H-NMR spectroscopic metabolic profiling of urine

Recent Advancement in NMR Based Plant Metabolomics: Techniques, Tools, and Analytical Approaches

Plant metabolomics, a rapidly advancing field within plant biology, is dedicated to comprehensively exploring the intricate array of small molecules in plant systems. This entails precisely gathering comprehensive chemical data, detecting numerous metabolites, and ensuring accurate molecular identification. Nuclear magnetic resonance (NMR) spectroscopy, with its detailed chemical insights, is crucial in obtaining metabolite profiles. Its widespread application spans various research disciplines, aiding in comprehending chemical reactions, kinetics, and molecule characterization. Biotechnological advancements have further expanded NMR's utility in metabolomics, particularly in identifying disease biomarkers across diverse fields such as agriculture, medicine, and pharmacology. This review covers the stages of NMR-based metabolomics, including historical aspects and limitations, with sample preparation, data acquisition, spectral processing, analysis, and their application parts.

Multi-scale Gaussian/Haar wavelet strategies coupled with sub-window factor analysis for an accurate alignment in nontargeted metabolic profiling to enhance herbal origin discrimination capability

Metabolic dataset can provide an overview of different herbal origin, which is conducted by some statistical procedures. Such results often deviate to a certain degree, due to peaks shifts in chromatographic signals. In order to solve this problem, an improved algorithm of combining sub-window factor analysis with the mass spectrum information is proposed. The algorithm uses a peak detection approach derived either from multi-scale Gaussian function or Haar wavelet to locate the peaks with different application scope; the candidate drift points at each peak are estimated by Fast Fourier transform cross correlation; Specifically, the best drift points at each candidate peaks are confirmed by sub-window factor analysis and mass spectrum information in nontargeted metabolic profiling. Finally, the peak regions were aligned against a reference chromatogram, and the non-peak regions were used linear interpolation. The chromatographic signals of 30 Bupleurum samples were aligned as an illustration of this algorithm, and they could be well distinguished using some statistical procedures. The result demonstrates that the presented method is stronger than other mass-spectra based algorithms, when facing the alignment of some co-eluted peaks.

Quantitative Quantum Mechanical Spectral Analysis (qQMSA) of (1)H NMR spectra of complex mixtures and biofluids

The quantitative interpretation of (1)H NMR spectra of mixtures like the biofluids is a demanding task due to spectral complexity and overlap. Complications may arise also from water suppression, T2-editing, protein interactions, relaxation differences of the species, experimental artifacts and, furthermore, the spectra may contain unknown components and macromolecular background which cannot be easily separated from baseline. In this work, tools and strategies for quantitative Quantum Mechanical Spectral Analysis (qQMSA) of (1)H NMR spectra from complex mixtures were developed and systematically assessed. In the present approach, the signals of well-defined, stoichiometric components are described by a QM model, while the background is described by a multiterm baseline function and the unknown signals using optimizable and adjustable lines, regular multiplets or any spectral structures which can be composed from spectral lines. Any prior knowledge available from the spectrum can also be added to the model. Fitting strategies for weak and strongly overlapping spectral systems were developed and assessed using two basic model systems, the metabolite mixtures without and with macromolecular (serum) background. The analyses show that if the spectra are measured in high-throughput manner, the consistent absolute quantification demands some calibration to compensate the different response factors of the protons and compounds. On the other hand, the results show that also the T2-edited spectra can be measured so that they obey well the QM rules. In general, qQMSA exploits and interprets the spectral information in maximal way taking full advantage from the QM properties of the spectra and, at the same time, offers chemical confidence which means that individual components can be identified with high confidence on the basis of their accurate spectral parameters.

Automated annotation and quantification of metabolites in 1H NMR data of biological origin

In (1)H NMR metabolomic datasets, there are often over a thousand peaks per spectrum, many of which change position drastically between samples. Automatic alignment, annotation, and quantification of all the metabolites of interest in such datasets have not been feasible. In this work we propose a fully automated annotation and quantification procedure which requires annotation of metabolites only in a single spectrum. The reference database built from that single spectrum can be used for any number of (1)H NMR datasets with a similar matrix. The procedure is based on the generalized fuzzy Hough transform (GFHT) for alignment and on Principal-components analysis (PCA) for peak selection and quantification. We show that we can establish quantities of 21 metabolites in several (1)H NMR datasets and that the procedure is extendable to include any number of metabolites that can be identified in a single spectrum. The procedure speeds up the quantification of previously known metabolites and also returns a table containing the intensities and locations of all the peaks that were found and aligned but not assigned to a known metabolite. This enables both biopattern analysis of known metabolites and data mining for new potential biomarkers among the unknowns.

Evaluation of standard and advanced preprocessing methods for the univariate analysis of blood serum 1H-NMR spectra

Proton nuclear magnetic resonance ((1)H-NMR)-based metabolomics enables the high-resolution and high-throughput assessment of a broad spectrum of metabolites in biofluids. Despite the straightforward character of the experimental methodology, the analysis of spectral profiles is rather complex, particularly due to the requirement of numerous data preprocessing steps. Here, we evaluate how several of the most common preprocessing procedures affect the subsequent univariate analyses of blood serum spectra, with a particular focus on how the standard methods perform compared to more advanced examples. Carr-Purcell-Meiboom-Gill 1D (1)H spectra were obtained for 240 serum samples from healthy subjects of the Asklepios study. We studied the impact of different preprocessing steps--integral (standard method) and probabilistic quotient normalization; no, equidistant (standard), and adaptive-intelligent binning; mean (standard) and maximum bin intensity data summation--on the resonance intensities of three different types of metabolites: triglycerides, glucose, and creatinine. The effects were evaluated by correlating the differently preprocessed NMR data with the independently measured metabolite concentrations. The analyses revealed that the standard methods performed inferiorly and that a combination of probabilistic quotient normalization after adaptive-intelligent binning and maximum intensity variable definition yielded the best overall results (triglycerides, R = 0.98; glucose, R = 0.76; creatinine, R = 0.70). Therefore, at least in the case of serum metabolomics, these or equivalent methods should be preferred above the standard preprocessing methods, particularly for univariate analyses. Additional optimization of the normalization procedure might further improve the analyses.

icoshift: A versatile tool for the rapid alignment of 1D NMR spectra

The increasing scientific and industrial interest towards metabonomics takes advantage from the high qualitative and quantitative information level of nuclear magnetic resonance (NMR) spectroscopy. However, several chemical and physical factors can affect the absolute and the relative position of an NMR signal and it is not always possible or desirable to eliminate these effects a priori. To remove misalignment of NMR signals a posteriori, several algorithms have been proposed in the literature. The icoshift program presented here is an open source and highly efficient program designed for solving signal alignment problems in metabonomic NMR data analysis. The icoshift algorithm is based on correlation shifting of spectral intervals and employs an FFT engine that aligns all spectra simultaneously. The algorithm is demonstrated to be faster than similar methods found in the literature making full-resolution alignment of large datasets feasible and thus avoiding down-sampling steps such as binning. The algorithm uses missing values as a filling alternative in order to avoid spectral artifacts at the segment boundaries. The algorithm is made open source and the Matlab code including documentation can be downloaded from www.models.life.ku.dk.

Application of NMR in plant metabolomics: techniques, problems and prospects

The present state-of-the-art of NMR in plant metabolomics is reviewed. Attention is paid to the different practical aspects of the application of NMR. The sample preparation, the measurement of the spectrum, quantitative aspects and data analysis are discussed. Each stage has its specific problems, which are pointed out and recommendations are made.

A solution to the 1D NMR alignment problem using an extended generalized fuzzy Hough transform and mode support

This paper approaches the problem of intersample peak correspondence in the context of later applying statistical data analysis techniques to 1D 1H-nuclear magnetic resonance (NMR) data. Any data analysis methodology will fail to produce meaningful results if the analyzed data table is not synchronized, i.e., each analyzed variable frequency (Hz) does not originate from the same chemical source throughout the entire dataset. This is typically the case when dealing with NMR data from biological samples. In this paper, we present a new state of the art for solving this problem using the generalized fuzzy Hough transform (GFHT). This paper describes significant improvements since the method was introduced for NMR datasets of plasma in Csenki et al. (Anal Bioanal Chem 389:875-885, 15) and is now capable of synchronizing peaks from more complex datasets such as urine as well as plasma data. We present a novel way of globally modeling peak shifts using principal component analysis, a new algorithm for calculating the transform and an effective peak detection algorithm. The algorithm is applied to two real metabonomic 1H-NMR datasets and the properties of the method are compared to bucketing. We implicitly prove that GFHT establishes the objectively true correspondence. Desirable features of the GFHT are: (1) intersample peak correspondence even if peaks change order on the frequency axis and (2) the method is symmetric with respect to the samples.

Proof of principle of a generalized fuzzy Hough transform approach to peak alignment of one-dimensional 1H NMR data

In metabolic profiling, multivariate data analysis techniques are used to interpret one-dimensional (1D) 1H NMR data. Multivariate data analysis techniques require that peaks are characterised by the same variables in every spectrum. This location constraint is essential for correct comparison of the intensities of several NMR spectra. However, variations in physicochemical factors can cause the locations of the peaks to shift. The location prerequisite may thus not be met, and so, to solve this problem, alignment methods have been developed. However, current state-of-the-art algorithms for data alignment cannot resolve the inherent problems encountered when analysing NMR data of biological origin, because they are unable to align peaks when the spatial order of the peaks changes-a commonly occurring phenomenon. In this paper a new algorithm is proposed, based on the Hough transform operating on an image representation of the NMR dataset that is capable of correctly aligning peaks when existing methods fail. The proposed algorithm was compared with current state-of-the-art algorithms operating on a selected plasma dataset to demonstrate its potential. A urine dataset was also processed using the algorithm as a further demonstration. The method is capable of successfully aligning the plasma data but further development is needed to address more challenging applications, for example urine data.

Application of metabonomics in drug development

Metabolic profiling (metabonomics/metabolomics) is the untargeted analysis of metabolic composition in a biological sample, and is principally aimed at biomarker discovery. The frequent use of noninvasive biofluid analysis in metabonomics is suited to the clinic and facilitates dynamic monitoring. Analytical protocols for metabolic biomarkers are potentially robust because a metabolite is the same chemical entity irrespective of its origin, facilitating ‘bench-to-bedside’ translational research. Metabonomics can make an impact at several points in the drug-development process: target identification; lead discovery and optimization; preclinical efficacy and safety assessment; mode-of-action and mechanistic toxicology; patient stratification; and clinical pharmacological monitoring. This review describes and exemplifies the latest developments in each of these areas, including the impact of new data and chemical analytical techniques. The future goals for metabonomics are the validation of existing biomarkers, in terms of mechanism and translation to man, together with a focus on characterizing the individual (‘personalized healthcare’).