Quantification and identification of components in solution mixtures from 1D proton NMR spectra using singular value decomposition

One-dimensional proton NMR spectra of complex solutions provide rich molecular information, but limited chemical shift dispersion creates peak overlap that often leads to difficulty in peak identification and analyte quantification. Modern high-field NMR spectrometers provide high digital resolution with improved peak dispersion. We took advantage of these spectral qualities and developed a quantification method based on linear least-squares fitting using singular value decomposition (SVD). The linear least-squares fitting of a mixture spectrum was performed on the basis of reference spectra from individual small-molecule analytes. Each spectrum contained an internal quantitative reference (e.g., DSS-d6 or other suitable small molecules) by which the intensity of the spectrum was scaled. Normalization of the spectrum facilitated quantification based on peak intensity using linear least-squares fitting analysis. This methodology provided quantification of individual analytes as well as chemical identification. The analysis of small-molecule analytes over a wide concentration range indicated the accuracy and reproducibility of the SVD-based quantification. To account for the contribution from residual protein, lipid or polysaccharide in solution, a reference spectrum showing the macromolecules or aggregates was obtained using a diffusion-edited 1D proton NMR analysis. We demonstrated this approach with a mixture of small-molecule analytes in the presence of macromolecules (e.g., protein). The results suggested that this approach should be applicable to the quantification and identification of small-molecule analytes in complex biological samples.     

Web server based complex mixture analysis by NMR

Comprehensive metabolite identification and quantification of complex biological mixtures are central aspects of metabolomics. NMR shows excellent promise for these tasks. An automated fingerprinting strategy is presented, termed COLMAR query, which screens NMR chemical shift lists or raw 1D NMR cross sections taken from covariance total correlation spectroscopy (TOCSY) spectra or other multidimensional NMR spectra against an NMR spectral database. Cross peaks are selected using local clustering to avoid ambiguities between chemical shifts and scalar J-coupling effects. With the use of three different algorithms, the corresponding chemical shift list is then screened against chemical shift lists extracted from 1D spectra of a NMR database. The resulting query scores produced by forward assignment, reverse assignment, and bipartite weighted-matching algorithms are combined into a consensus score, which provides a robust means for identifying the correct compound. The approach is demonstrated for a metabolite model mixture that is screened against the metabolomics BioMagResDatabank (BMRB). This NMR-based compound identification approach has been implemented in a public Web server that allows the efficient analysis of a wide range of metabolite mixtures.     

Virtual chromatographic resolution enhancement in cryoflow LC-NMR experiments via statistical total correlation spectroscopy

A new approach to enhancing information recovery from cryogenic probe "on-flow" LC-NMR spectroscopic analyses of complex biological mixtures is demonstrated using a variation on the statistical total correlation spectroscopy (STOCSY) method. Cryoflow probe technology enables sensitive and efficient NMR detection of metabolites on-flow, and the rapid spectral scanning allows multiple spectra to be collected over chromatographic peaks containing several species with similar, but nonidentical, retention times. This enables 1H NMR signal connectivities between close-eluting metabolites to be identified resulting in a "virtual" chromatographic resolution enhancement visualized directly in the NMR spectral projection. We demonstrate the applicability of the approach for structure assignment of drug and endogenous metabolites in urine. This approach is of wide general applicability to any complex mixture analysis problem involving chromatographic peak overlap and with particular application in metabolomics and metabonomics.     

Between-person comparison of metabolite fitting for NMR-based quantitative metabolomics

Nuclear magnetic resonance (NMR) spectroscopy is widely used as an analytical platform for metabolomics. Many studies make use of 1D spectra, which have the advantages of relative simplicity and rapid acquisition times. The spectral data can then be analyzed either with a chemometric workflow or by an initial deconvolution or fitting step to generate a list of identified metabolites and associated sample concentrations. Various software tools exist to simplify the fitting process, but at least for 1D spectra, this still requires a degree of skilled operator input. It is of critical importance that we know how much person-to-person variability affects the results, in order to be able to judge between different studies. Here we tested a commercially available software package (Chenomx' NMR Suite) for fitting metabolites to a set of NMR spectra of yeast extracts and compared the output of five different people for both metabolite identification and quantitation. An initial comparison showed good agreement for a restricted set of common metabolites with characteristic well-resolved resonances but wide divergence in the overall identities and number of compounds fitted; refitting according to an agreed set of metabolites and spectral processing approach increased the total number of metabolites fitted but did not dramatically increase the quality of the metabolites that could be fitted without prior knowledge about peak identity. Hence, robust peak assignments are required in advance of manual deconvolution, when the widest range of metabolites is desired. However, very low concentration metabolites still had high coefficients of variation even with shared information on peak assignment. Overall, the effect of the person was less than the experimental group (in this case, sampling method) for almost all of the metabolites.     

"Add to subtract": a simple method to remove complex background signals from the 1H nuclear magnetic resonance spectra of mixtures

Because of its highly reproducible and quantitative nature and minimal requirements for sample preparation or separation, (1)H nuclear magnetic resonance (NMR) spectroscopy is widely used for profiling small-molecule metabolites in biofluids. However (1)H NMR spectra contain many overlapped peaks. In particular, blood serum/plasma and diabetic urine samples contain high concentrations of glucose, which produce strong peaks between 3.2 ppm and 4.0 ppm. Signals from most metabolites in this region are overwhelmed by the glucose background signals and become invisible. We propose a simple "Add to Subtract" background subtraction method and show that it can reduce the glucose signals by 98% to allow retrieval of the hidden information. This procedure includes adding a small drop of concentrated glucose solution to the sample in the NMR tube, mixing, waiting for an equilibration time, and acquisition of a second spectrum. The glucose-free spectra are then generated by spectral subtraction using Bruker Topspin software. Subsequent multivariate statistical analysis can then be used to identify biomarker candidate signals for distinguishing different types of biological samples. The principle of this approach is generally applicable for all quantitative spectral data and should find utility in a variety of NMR-based mixture analyses as well as in metabolite profiling.     

Statistical correlation and projection methods for improved information recovery from diffusion-edited NMR spectra of biological samples

Although NMR spectroscopic techniques coupled with multivariate statistics can yield much useful information for classifying biological samples based on metabolic profiles, biomarker identification remains a time-consuming and complex procedure involving separation methods, two-dimensional NMR, and other spectroscopic tools. We present a new approach to aid complex biomixture analysis that combines diffusion ordered (DO) NMR spectroscopy with statistical total correlation spectroscopy (STOCSY) and demonstrate its application in the characterization of urinary biomarkers and enhanced information recovery from plasma NMR spectra. This method relies on calculation and display of the covariance of signal intensities from the various nuclei on the same molecule across a series of spectra collected under different pulsed field gradient conditions that differentially attenuate the signal intensities according to translational molecular diffusion rates. We term this statistical diffusion-ordered spectroscopy (S-DOSY). We also have developed a new visualization tool in which the apparent diffusion coefficients from DO spectra are projected onto a 1D NMR spectrum (diffusion-ordered projection spectroscopy, DOPY). Both methods either alone or in combination have the potential for general applications to any complex mixture analysis where the sample contains compounds with a range of diffusion coefficients.     

Hierarchical Alignment and Full Resolution Pattern Recognition of 2D NMR Spectra: Application to Nematode Chemical Ecology

Nuclear magnetic resonance (NMR) is the most widely used nondestructive technique in analytical chemistry. In recent years, it has been applied to metabolic profiling due to its high reproducibility, capacity for relative and absolute quantification, atomic resolution, and ability to detect a broad range of compounds in an untargeted manner. While one-dimensional (1D) (1)H NMR experiments are popular in metabolic profiling due to their simplicity and fast acquisition times, two-dimensional (2D) NMR spectra offer increased spectral resolution as well as atomic correlations, which aid in the assignment of known small molecules and the structural elucidation of novel compounds. Given the small number of statistical analysis methods for 2D NMR spectra, we developed a new approach for the analysis, information recovery, and display of 2D NMR spectral data. We present a native 2D peak alignment algorithm we term HATS, for hierarchical alignment of two-dimensional spectra, enabling pattern recognition (PR) using full-resolution spectra. Principle component analysis (PCA) and partial least squares (PLS) regression of full resolution total correlation spectroscopy (TOCSY) spectra greatly aid the assignment and interpretation of statistical pattern recognition results by producing back-scaled loading plots that look like traditional TOCSY spectra but incorporate qualitative and quantitative biological information of the resonances. The HATS-PR methodology is demonstrated here using multiple 2D TOCSY spectra of the exudates from two nematode species: Pristionchus pacificus and Panagrellus redivivus. We show the utility of this integrated approach with the rapid, semiautomated assignment of small molecules differentiating the two species and the identification of spectral regions suggesting the presence of species-specific compounds. These results demonstrate that the combination of 2D NMR spectra with full-resolution statistical analysis provides a platform for chemical and biological studies in cellular biochemistry, metabolomics, and chemical ecology.     

Strategy for automated analysis of dynamic metabolic mixtures by NMR. Application to an insect venom

Elucidation of the composition of chemical-biological samples is a main focus of systems biology and metabolomics. Due to the inherent complexity of these mixtures, reliable, efficient, and potentially automatable methods are needed to identify the underlying metabolites and natural products. Because of its rich chemical information content, nuclear magnetic resonance (NMR) spectroscopy has a unique potential for this task. Here we present a generalization and application of a recently introduced NMR data collection, processing, and analysis strategy that circumvents the need for extensive purification and hyphenation prior to analysis. It uses covariance TOCSY NMR spectra measured on a 1-mm high-temperature cryogenic probe that are analyzed by a spectral trace clustering algorithm yielding 1D NMR spectra of the individual components for their unambiguous identification. The method is demonstrated on a metabolic model mixture and is then applied to the unpurified venom mixture of an individual walking stick insect that contains several slowly interconverting and closely related metabolites.     

Metabolite projection analysis for fast identification of metabolites in metabonomics. Application in an amiodarone study

The assignment of significantly changed NMR signals, which were identified with the help of multivariate models, to individual metabolites in biofluids is a manual and tedious task requiring knowledge in chemometrics and NMR spectroscopy. Metabolite projection analysis, introduced in this work, allows automatic linking of multivariate models with metabolites by skipping the level of manual NMR signal identification. The method depends on the projection of sets of metabolite NMR spectra from a database into PCA or PLS models of NMR spectra of biofluid samples. Metabolites that are significantly changed can be identified graphically in metabolite projection plots or numerically as projected virtual concentration. The method is demonstrated together with a newly introduced algorithm for refined nonequidistant binning using a metabonomics study with amiodarone as administered drug. Amiodarone can induce phospholipidosis in the lung and liver, which is accompanied by associated organ toxicity in these organs. It is shown how metabolite projection analysis allows easy and fast tentative assignment of all structures of metabolites whose concentrations in the urine samples significantly changed upon dosage. These metabolites had also been identified previously by manually interpreting the multivariate models and spectra. Among these metabolites, phenylacetylglycine was also identified as being significantly increased. This metabolite has recently been proposed as urinary biomarker for phospholipidosis.     

Method for determining molar concentrations of metabolites in complex solutions from two-dimensional 1H-13C NMR spectra

One-dimensional (1D) (1)H nuclear magnetic resonance (NMR) spectroscopy is used extensively for high-throughput analysis of metabolites in biological fluids and tissue extracts. Typically, such spectra are treated as multivariate statistical objects rather than as collections of quantifiable metabolites. We report here a two-dimensional (2D) (1)H-(13)C NMR strategy (fast metabolite quantification, FMQ, by NMR) for identifying and quantifying the approximately 40 most abundant metabolites in biological samples. To validate this technique, we prepared mixtures of synthetic compounds and extracts from Arabidopsis thaliana, Saccharomyces cerevisiae, and Medicago sativa. We show that accurate (technical error 2.7%) molar concentrations can be determined in 12 min using our quantitative 2D (1)H-(13)C NMR strategy. In contrast, traditional 1D (1)H NMR analysis resulted in 16.2% technical error under nearly ideal conditions. We propose FMQ by NMR as a practical alternative to 1D (1)H NMR for metabolomics studies in which 50-mg (extract dry weight) samples can be obtained.     

Investigations of the effects of gender, diurnal variation, and age in human urinary metabolomic profiles

Metabolomics may have the capacity to revolutionize disease diagnosis through the identification of scores of metabolites that vary during environmental, pathogenic, or toxicological insult. NMR spectroscopy has become one of the main tools for measuring these changes since an NMR spectrum can accurately identify metabolites and their concentrations. The predominant approach in analyzing NMR data has been through the technique of spectral binning. However, identification of spectral areas in an NMR spectrum is insufficient for diagnostic evaluation, since it is unknown whether areas of interest are strictly caused by metabolic changes or are simply artifacts. In this paper, we explore differences in gender, diurnal variation, and age in a human population. We use the example of gender differences to compare traditional spectral binning techniques (NMR spectral areas) to novel targeted profiling techniques (metabolites and their concentrations). We show that targeted profiling produces robust models, generates accurate metabolite concentration data, and provides data that can be used to help understand metabolic differences in a healthy population. Metabolites relating to mitochondrial energy metabolism were found to differentiate gender and age. Dietary components and some metabolites related to circadian rhythms were found to differentiate time of day urine collection. The mechanisms by which these differences arise will be key to the discovery of new diagnostic tests and new understandings of the mechanism of disease.     

Measurement of absolute concentrations of individual compounds in metabolite mixtures by gradient-selective time-zero 1H-13C HSQC with two concentration references and fast maximum likelihood reconstruction analysis

Time-zero 2D (13)C HSQC (HSQC(0)) spectroscopy offers advantages over traditional 2D NMR for quantitative analysis of solutions containing a mixture of compounds because the signal intensities are directly proportional to the concentrations of the constituents. The HSQC(0) spectrum is derived from a series of spectra collected with increasing repetition times within the basic HSQC block by extrapolating the repetition time to zero. Here we present an alternative approach to data collection, gradient-selective time-zero (1)H-(13)C HSQC(0) in combination with fast maximum likelihood reconstruction (FMLR) data analysis and the use of two concentration references for absolute concentration determination. Gradient-selective data acquisition results in cleaner spectra, and NMR data can be acquired in both constant-time and non-constant-time mode. Semiautomatic data analysis is supported by the FMLR approach, which is used to deconvolute the spectra and extract peak volumes. The peak volumes obtained from this analysis are converted to absolute concentrations by reference to the peak volumes of two internal reference compounds of known concentration: DSS (4,4-dimethyl-4-silapentane-1-sulfonic acid) at the low concentration limit (which also serves as chemical shift reference) and MES (2-(N-morpholino)ethanesulfonic acid) at the high concentration limit. The linear relationship between peak volumes and concentration is better defined with two references than with one, and the measured absolute concentrations of individual compounds in the mixture are more accurate. We compare results from semiautomated gsHSQC(0) with those obtained by the original manual phase-cycled HSQC(0) approach. The new approach is suitable for automatic metabolite profiling by simultaneous quantification of multiple metabolites in a complex mixture.     

Self-consistent metabolic mixture analysis by heteronuclear NMR. Application to a human cancer cell line

Elucidation of the chemical composition of biological samples is a main focus of systems biology and metabolomics. In order to comprehensively study these complex mixtures, reliable, efficient, and automatable methods are needed to identify and quantify the underlying metabolites and natural products. Because of its rich information content, nuclear magnetic resonance (NMR) spectroscopy has a unique potential for this task. Here we present a generalization of the recently introduced homonuclear TOCSY-based DemixC method to heteronuclear HSQC-TOCSY NMR spectroscopy. The approach takes advantage of the high resolution afforded along the (13)C dimension due to the narrow (13)C line widths for the identification of spin systems and compounds. The method combines information from both 1D (13)C and (1)H traces by querying them against an NMR spectral database using our COLMAR query web server. The complementarity of (13)C and (1)H spectral information improves the robustness of compound identification. The method is demonstrated for a metabolic model mixture and is then applied to an extract from DU145 human prostate cancer cells.     

Deconvolution of two-dimensional NMR spectra by fast maximum likelihood reconstruction: application to quantitative metabolomics

We have developed an algorithm called fast maximum likelihood reconstruction (FMLR) that performs spectral deconvolution of 1D-2D NMR spectra for the purpose of accurate signal quantification. FMLR constructs the simplest time-domain model (e.g., the model with the fewest number of signals and parameters) whose frequency spectrum matches the visible regions of the spectrum obtained from identical Fourier processing of the acquired data. We describe the application of FMLR to quantitative metabolomics and demonstrate the accuracy of the method by analysis of complex, synthetic mixtures of metabolites and liver extracts. The algorithm demonstrates greater accuracy (0.5-5.0% error) than peak height analysis and peak integral analysis with greatly reduced operator intervention. FMLR has been implemented in a Java-based framework that is available for download on multiple platforms and is interoperable with popular NMR display and processing software. Two-dimensional (1)H-(13)C spectra of mixtures can be acquired with acquisition times of 15 min and analyzed by FMLR in the range of 2-5 min per spectrum to identify and quantify constituents present at concentrations of 0.2 mM or greater.     

Heteronuclear 1H-31P statistical total correlation NMR spectroscopy of intact liver for metabolic biomarker assignment: application to galactosamine-induced hepatotoxicity

As part of our ongoing development of methods for enhanced biomarker information recovery from spectroscopic data we present the first example of a new hetero-nuclear statistical total correlation spectroscopy (HET-STOCSY) approach applied to intact tissue samples collected as part of a toxicological study. One-dimensional 1H and 31P-{1H} magic angle spinning (MAS) NMR spectra of intact liver samples after galactosamine (galN) treatment to rats and after cotreatment of galN plus uridine were collected at 275 K. Individual samples were also followed by 1H and 31P-{1H} MAS NMR through time generating time dependent modulations in metabolite signatures relating to toxicity. High-resolution 1H NMR spectra of urine and plasma and clinical chemical data were also collected to establish a biological framework in which to place these novel statistical heterospectroscopic data. In HET-STOCSY, calculation of the covariance between the 31P-{1H} and 1H NMR signals of phosphorus containing metabolites allows their molecular connectivities to be established and the construction of virtual two-dimensional heteronuclear correlation spectra that connect all protons on the molecule to the heteroatom. We show how HET-STOCSY applied to MAS NMR spectra of liver samples can be used to augment biomarker detection. This approach is generic and can be applied to correlate the covarying signals for any spin-active nuclei where there is parallel or serial collection of data.     

Statistical heterospectroscopy, an approach to the integrated analysis of NMR and UPLC-MS data sets: application in metabonomic toxicology studies

Statistical heterospectroscopy (SHY) is a new statistical paradigm for the coanalysis of multispectroscopic data sets acquired on multiple samples. This method operates through the analysis of the intrinsic covariance between signal intensities in the same and related molecules measured by different techniques across cohorts of samples. The potential of SHY is illustrated using both 600-MHz 1H NMR and UPLC-TOFMS data obtained from control rat urine samples (n = 54) and from a corresponding hydrazine-treated group (n = 58). We show that direct cross-correlation of spectral parameters, viz. chemical shifts from NMR and m/z data from MS, is readily achievable for a variety of metabolites, which leads to improved efficiency of molecular biomarker identification. In addition to structure, higher level biological information can be obtained on metabolic pathway activity and connectivities by examination of different levels of the NMR to MS correlation and anticorrelation matrixes. The SHY approach is of general applicability to complex mixture analysis, if two or more independent spectroscopic data sets are available for any sample cohort. Biological applications of SHY as demonstrated here show promise as a new systems biology tool for biomarker recovery.     

Chemicalome and metabolome matching approach to elucidating biological metabolic networks of complex mixtures

Global metabolite identification of complex compound mixtures in biological systems is a very challenging task. Herein, we developed and validated a chemicalome to metabolome matching approach by taking herbal medicine as an example to delineate the metabolic networks of complex systems. This approach consists of five steps of data processing including raw data output, endogenous background subtraction, parent compound and metabolite differentiation, chemicalome to metabolome correlation, and the final validation via manual fragment comparison. Chemicalome to metabolome correlation, the core step of this approach, was performed based on matching the accurate mass differences of pseudomolecular ions between them with the accurate mass changes of known metabolic pathways and validating the matches by validation ions. A step-forward approach that confers a gradual identification of metabolites generated from different steps (1-4) and types (degradation, phase I/II, or mixed) of metabolic reactions was further proposed for chemicalome to metabolome matching. This approach was validated to be very useful and powerful for the metabolite identification of a single compound, a homologous compound mixture, and a complex herbal system. Using this approach, all metabolites (162) detected from urine samples of rats treated with Mai-Luo-Ning injection could be linked to their respective parent compounds, and 143 of them were supported by the final validation via manual fragment analysis. In most cases, more than 80% of the automatic matching results could be supported by the manual fragment validations. A complex metabolic network showing all the possible links between precursors and metabolites was successfully constructed. This study provides a generally applicable approach to global metabolite identification of complex compound mixtures in complex matrixes.     

Targeted profiling: quantitative analysis of 1H NMR metabolomics data

Extracting meaningful information from complex spectroscopic data of metabolite mixtures is an area of active research in the emerging field of "metabolomics", which combines metabolism, spectroscopy, and multivariate statistical analysis (pattern recognition) methods. Chemometric analysis and comparison of 1H NMR1 spectra is commonly hampered by intersample peak position and line width variation due to matrix effects (pH, ionic strength, etc.). Here a novel method for mixture analysis is presented, defined as "targeted profiling". Individual NMR resonances of interest are mathematically modeled from pure compound spectra. This database is then interrogated to identify and quantify metabolites in complex spectra of mixtures, such as biofluids. The technique is validated against a traditional "spectral binning" analysis on the basis of sensitivity to water suppression (presaturation, NOESY-presaturation, WET, and CPMG), relaxation effects, and NMR spectral acquisition times (3, 4, 5, and 6 s/scan) using PCA pattern recognition analysis. In addition, a quantitative validation is performed against various metabolites at physiological concentrations (9 microM-8 mM). "Targeted profiling" is highly stable in PCA-based pattern recognition, insensitive to water suppression, relaxation times (within the ranges examined), and scaling factors; hence, direct comparison of data acquired under varying conditions is made possible. In particular, analysis of metabolites at low concentration and overlapping regions are well suited to this analysis. We discuss how targeted profiling can be applied for mixture analysis and examine the effect of various acquisition parameters on the accuracy of quantification.     

MASPIC: intensity-based tandem mass spectrometry scoring scheme that improves peptide identification at high confidence

Algorithmic search engines bridge the gap between large tandem mass spectrometry data sets and the identification of proteins associated with biological samples. Improvements in these tools can greatly enhance biological discovery. We present a new scoring scheme for comparing tandem mass spectra with a protein sequence database. The MASPIC (Multinomial Algorithm for Spectral Profile-based Intensity Comparison) scorer converts an experimental tandem mass spectrum into a m/z profile of probability and then scores peak lists from potential candidate peptides using a multinomial distribution model. The MASPIC scoring scheme incorporates intensity, spectral peak density variations, and m/z error distribution associated with peak matches into a multinomial distribution. The scoring scheme was validated on two standard protein mixtures and an additional set of spectra collected on a complex ribosomal protein mixture from Rhodopseudomonas palustris. The results indicate a 5-15% improvement over Sequest for high-confidence identifications. The performance gap grows as sequence database size increases. Additional tests on spectra from proteinase-K digest data showed similar performance improvements demonstrating the advantages in using MASPIC for studying proteins digested with less specific proteases. All these investigations show MASPIC to be a versatile and reliable system for peptide tandem mass spectral identification.     

Structural elucidation and quantification of phenolic conjugates present in human urine after tea intake

In dietary polyphenol exposure studies, annotation and identification of urinary metabolites present at low (micromolar) concentrations are major obstacles. To determine the biological activity of specific components, it is necessary to have the correct structures and the quantification of the polyphenol-derived conjugates present in the human body. We present a procedure for identification and quantification of metabolites and conjugates excreted in human urine after single bolus intake of black or green tea. A combination of a solid-phase extraction (SPE) preparation step and two high pressure liquid chromatography (HPLC)-based analytical platforms was used, namely, accurate mass fragmentation (HPLC-FTMS(n)) and mass-guided SPE-trapping of selected compounds for nuclear magnetic resonance spectroscopy (NMR) measurements (HPLC-TOFMS-SPE-NMR). HPLC-FTMS(n) analysis led to the annotation of 138 urinary metabolites, including 48 valerolactone and valeric acid conjugates. By combining the results from MS(n) fragmentation with the one-dimensional (1D)-(1)H NMR spectra of HPLC-TOFMS-SPE-trapped compounds, we elucidated the structures of 36 phenolic conjugates, including the glucuronides of 3',4'-di- and 3',4',5'-trihydroxyphenyl-γ-valerolactone, three urolithin glucuronides, and indole-3-acetic acid glucuronide. We also obtained 26 h-quantitative excretion profiles for specific valerolactone conjugates. The combination of the HPLC-FTMS(n) and HPLC-TOFMS-SPE-NMR platforms results in the efficient identification and quantification of less abundant phenolic conjugates down to nanomoles of trapped amounts of metabolite corresponding to micromolar metabolite concentrations in urine.