Extraction of multiple pure component 1H and 13C NMR spectra from two mixtures: Novel solution obtained by sparse component analysis-based blind decomposition

Non-negative Least Squares Approach to Quantification of 1H Nuclear Magnetic Resonance Spectra of Human Urine

Because of its quantitative character and capability for high-throughput screening, ¹H nuclear magnetic resonance (NMR) spectroscopy is used extensively in the profiling of biofluids such as urine and blood plasma. However, the narrow frequency bandwidth of ¹H NMR spectroscopy leads to a severe overlap of the spectra of components present in the complex mixtures such as biofluids. Therefore, ¹H NMR-based metabolomics analysis is focused on targeted studies related to concentrations of the small number of metabolites. Here, we propose a library-based approach to quantify proportions of overlapping metabolites from ¹H NMR mixture spectra. The method boils down to the linear non-negative least squares (NNLS) problem, whereas proportions of the pure components contained in the library stand for the unknowns. The method is validated on an estimation of the proportions of (i) the 78 pure spectra, presumably related to type 2 diabetes mellitus (T2DM), from their synthetic linear mixture; (ii) metabolites present in 62 ¹H NMR spectra of urine of subjects with T2DM and 62 ¹H NMR spectra of urine of control subjects. In both cases, the in-house library of 210 pure component ¹H NMR spectra represented the design matrix in the related NNLS problem. The proposed method pinpoints 63 metabolites that in a statistically significant way discriminate the T2DM group from the control group and 46 metabolites discriminating control from the T2DM group. For several T2DM-discriminative metabolites, we prove their presence by independent analytical determination or by pointing out the corresponding findings in the published literature.

Library-assisted nonlinear blind separation and annotation of pure components from a single 1H nuclear magnetic resonance mixture spectra

Due to its capability for high-throughput screening ¹H nuclear magnetic resonance (NMR) spectroscopy is commonly used for metabolite research. The key problem in ¹H NMR spectroscopy of multicomponent mixtures is overlapping of component signals and that is increasing with the number of components, their complexity and structural similarity. It makes metabolic profiling, that is carried out through matching acquired spectra with metabolites from the library, a hard problem. Here, we propose a method for nonlinear blind separation of highly correlated components spectra from a single ¹H NMR mixture spectra. The method transforms a single nonlinear mixture into multiple high-dimensional reproducible kernel Hilbert Spaces (mRKHSs). Therein, highly correlated components are separated by sparseness constrained nonnegative matrix factorization in each induced RKHS. Afterwards, metabolites are identified through comparison of separated components with the library comprised of 160 pure components. Thereby, a significant number of them are expected to be related with diabetes type 2. Conceptually similar methodology for nonlinear blind separation of correlated components from two or more mixtures is presented in the Supplementary material. Single-mixture blind source separation is exemplified on: (i) annotation of five components spectra separated from one ¹H NMR model mixture spectra; (ii) annotation of fifty five metabolites separated from one ¹H NMR mixture spectra of urine of subjects with and without diabetes type 2. Arguably, it is for the first time a method for blind separation of a large number of components from a single nonlinear mixture has been proposed. Moreover, the proposed method pinpoints urinary creatine, glutamic acid and 5-hydroxyindoleacetic acid as the most prominent metabolites in samples from subjects with diabetes type 2, when compared to healthy controls.

Challenges in the decomposition of 2D NMR spectra of mixtures of small molecules

Analytical methods for mixtures of small molecules require specificity (is a certain molecule present in the mix?) and speciation capabilities. NMR spectroscopy has been a tool of choice for both of these issues since its early days, due to its quantitative (linear) response, sufficiently high resolving power and capabilities of inferring molecular structures from spectral features (even in the absence of a reference database). However, the analytical performances of NMR spectroscopy are being stretched by the increased complexity of the samples, the dynamic range of the components, and the need for a reasonable turnover time. One approach that has been actively pursued for disentangling the composition complexity is the use of 2D NMR spectroscopy. While any of the many experiments from this family will increase the spectral resolution, some are more apt for mixtures, as they are capable of unveiling signals belonging to whole molecules or fragments of it. Among the most popular ones, one can enumerate HSQC-TOCSY, DOSY and Maximum-Quantum (MaxQ) NMR spectroscopy. For multicomponent samples, the development of robust mathematical methods of signal decomposition would provide a clear edge towards identification. We have been pursuing, along these lines, Blind Source Separation (BSS). Here, the un-mixing of the spectra is achieved relying on correlations detected on a series of datasets. The series could be associated with samples of different relative composition or in a classically acquired 2D experiment by the mathematical laws underlying the construction of the indirect dimension, the one not recorded by the spectrometer. Many algorithms have been proposed for BSS in NMR spectroscopy since the seminal work of Nuzillard. In this paper, we use rather standard algorithms in BSS in order to disentangle NMR spectra. We show on simulated data (both 1D and 2D HSQC) that these approaches enable us to accurately disentangle multiple components, and provide good estimates for the concentrations of compounds. Furthermore, we show that after proper realignment of the signals, the same algorithms are able to disentangle real 1D NMR spectra. We obtain similar results on 2D HSQC spectra, where the BSS algorithms are able to successfully disentangle components, and provide even better estimates for concentrations.

Underdetermined blind separation of three-way fluorescence spectra of PAHs in water

In this work, underdetermined blind decomposition method is developed to recognize individual components from the three-way fluorescent spectra of their mixtures by using sparse component analysis (SCA). The mixing matrix is estimated from the mixtures using fuzzy data clustering algorithm together with the scatters corresponding to local energy maximum value in the time-frequency domain, and the spectra of object components are recovered by pseudo inverse technique. As an example, using this method three and four pure components spectra can be blindly extracted from two samples of their mixture, with similarities between resolved and reference spectra all above 0.80. This work opens a new and effective path to realize monitoring PAHs in water by three-way fluorescence spectroscopy technique.

A review of blind source separation in NMR spectroscopy

Fourier transform is the data processing naturally associated to most NMR experiments. Notable exceptions are Pulse Field Gradient and relaxation analysis, the structure of which is only partially suitable for FT. With the revamp of NMR of complex mixtures, fueled by analytical challenges such as metabolomics, alternative and more apt mathematical methods for data processing have been sought, with the aim of decomposing the NMR signal into simpler bits. Blind source separation is a very broad definition regrouping several classes of mathematical methods for complex signal decomposition that use no hypothesis on the form of the data. Developed outside NMR, these algorithms have been increasingly tested on spectra of mixtures. In this review, we shall provide an historical overview of the application of blind source separation methodologies to NMR, including methods specifically designed for the specificity of this spectroscopy.