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.     

Ratio analysis nuclear magnetic resonance spectroscopy for selective metabolite identification in complex samples

Metabolite identification in the complex NMR spectra of biological samples is a challenging task due to significant spectral overlap and limited signal-to-noise. In this study we present a new approach, RANSY (ratio analysis NMR spectroscopy), which identifies all the peaks of a specific metabolite on the basis of the ratios of peak heights or integrals. We show that the spectrum for an individual metabolite can be generated by exploiting the fact that the peak ratios for any metabolite in the NMR spectrum are fixed and proportional to the relative numbers of magnetically distinct protons. When the peak ratios are divided by their coefficients of variation derived from a set of NMR spectra, the generation of an individual metabolite spectrum is enabled. We first tested the performance of this approach using one-dimensional (1D) and two-dimensional (2D) NMR data of mixtures of synthetic analogues of common body fluid metabolites. Subsequently, the method was applied to (1)H NMR spectra of blood serum samples to demonstrate the selective identification of a number of metabolites. The RANSY approach, which does not need any additional NMR experiments for spectral simplification, is easy to perform and has the potential to aid in the identification of unknown metabolites using 1D or 2D NMR spectra in virtually any complex biological mixture.     

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.     

Timing and Spectral Diagrams of Magnetic Video-Pulse Excitation Influence on NMR Spin–Echo in Magnets

NMR methods are successfully used in obtaining the microscopic information on the internal magnetic structure and dynamical processes in magnets. The important application of NMR in magnets is also the investigation of structure and dynamics of domain walls. In this work, it is accomplished by observation of nuclear spin echo signals of nuclei arranged in domain walls of magnets at excitation by an additional magnetic video-pulse field. We present the first systematic study of timing and spectral diagrams of magnetic video-pulse influence on the NMR two-pulse echo in a number of magnets (ferromagnets, ferrites, half metals, intermetals). It is shown that the timing diagrams showing the dependence of two-pulse echo intensity on the temporal location of a magnetic video-pulse in respect to radio-frequency pulses and spectral diagram of this influence are defined mainly by local hyperfine field anisotropy and domain walls mobility. These diagrams could be used for the identification of the nature of NMR spectra in multidomain magnetic materials and to improve the resolution capacity of the NMR method in magnets.     

Use of Raman microscopy and band-target entropy minimization analysis to identify dyes in a commercial stamp. Implications for authentication and counterfeit detection

Raman microscopy was used in mapping mode to collect more than 1000 spectra in a 100 microm x 100 microm area from a commercial stamp. Band-target entropy minimization (BTEM) was then employed to unmix the mixture spectra in order to extract the pure component spectra of the samples. Three pure component spectral patterns with good signal-to-noise ratios were recovered, and their spatial distributions were determined. The three pure component spectral patterns were then identified as copper phthalocyanine blue, calcite-like material, and yellow organic dye material by comparison to known spectral libraries. The present investigation, consisting of (1) advanced curve resolution (blind-source separation) followed by (2) spectral data base matching, readily suggests extensions to authenticity and counterfeit studies of other types of commercial objects. The presence or absence of specific observable components form the basis for assessment. The present spectral analysis (BTEM) is applicable to highly overlapping spectral information. Since a priori information such as the number of components present and spectral libraries are not needed in BTEM, and since minor signals arising from trace components can be reconstructed, this analysis offers a robust approach to a wide variety of material problems involving authenticity and counterfeit issues.     

Kinetic and J-resolved statistical total correlation NMR spectroscopy approaches to structural information recovery in complex reacting mixtures: application to acyl glucuronide intramolecular transacylation reactions

We demonstrate here a new variant on a statistical spectroscopic method for recovering structural information on unstable intermediates formed in reaction mixtures. We exemplify this approach with respect to the internal acyl migration reactions of 1-beta-O-acyl glucuronides (AGs), which rearrange at neutral or slightly alkaline pH on a minute to hour time scale to yield a series of positional glucuronide ring isomers and alpha/beta anomers from the 1-beta (starting material), i.e. 2-beta, 2-alpha, 1-alpha, 3-beta, 3-alpha, and 4-beta, 4-alpha isomers together with the aglycon and alpha- and beta-glucuronic acid hydrolysis products. Multiple sequential 800 MHz cryoprobe (1)H NMR spectra (1D and 2D J-resolved, JRES) were collected on a 5.1 mM solution of a synthetic model drug glucuronide, 1-beta-O-acyl (S)-alpha-methyl phenylacetyl glucuronide (MPG) in 0.1 M sodium phosphate buffer in D2O at pD 7.4 over 18 h to monitor the reaction which leads to the formation of the eight positional isomers and hydrolysis products. As the reaction proceeds and new isomers form, the NMR signal intensities vary accordingly allowing the application of a novel kinetic variant on statistical total correlation spectroscopy (K-STOCSY) method to recover the connectivities between proton signals on the same reacting molecule based on their intensity covariance through time. We performed K-STOCSY analysis on both the standard 1D NMR spectra and the skyline projected singlets of the (1)H-(1)H JRES NMR spectra through time, i.e. the K-JRES-STOCSY experimental variant, which increases the effective spectral dispersion and is ideally suited for the analysis of heavily overlapped spin systems. High statistical correlations were observed between mutarotated alpha- and beta-anomers of individual positional isomers, as well as directly acyl migrated products and anticorrelation observed between signals from compounds that were being depleted as others increased, e.g. between the 1-beta and 2-alpha/2-beta isomers. This statistical kinetic approach enabled the recovery of structural connectivity information on all isomers allowing unequivocal resonance assignment, and this approach to spectroscopic information recovery has wider potential uses in the study of reactions that occur on the second-to-minute time scale in conditions where multiple sequential NMR spectra can be collected. JRES-STOCSY is also of potential use as a method for recovering spectroscopic information in highly overlapped NMR signals and spin systems in other types of complex mixture analysis.     

Development of 2D band-target entropy minimization and application to the deconvolution of multicomponent 2D nuclear magnetic resonance spectra

Spectral reconstruction from multicomponent spectroscopic data is the frequent primary goal in chemical system identification and exploratory chemometric studies. Various methods and techniques have been reported in the literature. However, few algorithms/methods have been devised for spectral recovery without the use of any a priori information. In the present studies, a higher dimensional entropy minimization method based on the BTEM algorithm (Widjaja, E.; Li, C.; Garland, M. Organometallics 2002, 21, 1991-1997.) and related techniques were extended to large-scale arrays, namely, 2D NMR spectroscopy. The performance of this novel method had been successfully verified on various real experimental mixture spectra from a series of randomized 2D NMR mixtures (COSY NMR and HSQC NMR). With the new algorithm and raw multicomponent NMR alone, it was possible to reconstruct the pure spectroscopic patterns and calculate the relative concentration of each species without recourse to any libraries or any other a priori information. The potential advantages of this novel algorithm and its implications for general chemical system identification of unknown mixtures are discussed.     

"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.     

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.     

Variable reference alignment: an improved peak alignment protocol for NMR spectral data with large intersample variation

In an effort to address the variable correspondence problem across large sample cohorts common in metabolomic/metabonomic studies, we have developed a prealignment protocol that aims to generate spectral segments sharing a common target spectrum. Under the assumption that a single reference spectrum will not correctly represent all spectra of a data set, the goal of this approach is to perform local alignment corrections on spectral regions which share a common "most similar" spectrum. A natural beneficial outcome of this procedure is the automatic definition of spectral segments, a feature that is not common to all alignment methods. This protocol is shown to specifically improve the quality of alignment in (1)H NMR data sets exhibiting large intersample compositional variation (e.g., pH, ionic strength). As a proof-of-principle demonstration, we have utilized two recently developed alignment algorithms specific to NMR data, recursive segment-wise peak alignment and interval correlated shifting, and applied them to two data sets composed of 15 aqueous cell line extract and 20 human urine (1)H NMR profiles. Application of this protocol represents a fundamental shift from current alignment methodologies that seek to correct misalignments utilizing a single representative spectrum, with the added benefit that it can be appended to any alignment algorithm.     

Quantification of high-resolution (1)H NMR spectra from rat brain extracts

Extracting quantitative information about absolute concentrations from high-resolution (1)H NMR spectra of complex mixtures such as brain extracts remains challenging. Partial overlap of resonances complicates integration, whereas simple line fitting algorithms cannot accommodate the spectral complexity of coupled spin systems. Here, it is shown that high-resolution (1)H NMR spectra of rat brain extracts from 11 distinct brain regions can be reproducibly quantified using a basis set of 29 compounds. The basis set is simulated with the density matrix formalism using complete prior knowledge of chemical shifts and scalar couplings. A crucial aspect to obtain reproducible results was the inclusion of a line shape distortion common among all 73 resonances of the 29 compounds. All metabolites could be quantified with <10% and <3% inter- and intrasubject variation, respectively.     

Two-dimensional correlation analysis of nuclear magnetic resonance metabonomics data

Two-dimensional (2D) correlation analysis has been used in this study to identify changes in complex nuclear magnetic resonance (NMR) metabonomics spectra of rat urine samples obtained during a study in which vasculitis (vascular injury), an important safety element in preclinical trials, was induced. Two types of correlation analysis were performed, along the variables and along the samples, and both 2D covariance and correlation coefficient maps were calculated. The binned and 'raw' NMR spectra were analyzed (0.04 and 0.001 ppm resolution, respectively). Good correlation was found among the major peaks of the binned spectra, and two groups of samples were identified using sample-sample 2D correlation maps. Much more complex correlation features were obtained from the 'raw' spectra, in which the specific, butterfly-like patterns were obtained in the covariance map but with only a few significant correlation coefficients in the corresponding 2D correlation maps. In terms of classification, the same group of the last nine spectra that indicated the end of the process and clustered in the 2D sample-sample covariance map of the binned data was also found in the 2D sample-sample covariance map of the raw NMR spectra but, again, not in the 2D correlation coefficient map. A discussion is given on the details of the application of the correlation analysis with regard to spectral data resolution, alignment, the effect of actual intensities of the NMR signal, and reference to various results from 2D correlation analysis of vibrational spectra.     

Effects of binary encoding on pattern recognition and library matching of spectral data

Scott, D.R., 1988. Effects of binary encoding on pattern recognition and library matching of spectral data. Chemometrics and Intelligent Laboratory Systems, 4: 47–63.Binary encoding is frequently employed in pattern recognition studies and matching of unknown against library spectral data. In this study the effect of binary encoding on pattern recognition and library searches is determined using the Hamming and Euclidean distance metrics. The effect on a full intensity spectrum is to compress the total information into a qualitative spectral data vector, the most basic information in the full spectrum. Geometrically, binary encoding of unit normalized spectral data can be visualized as shifting spectral points on the faces of the measurement space hypercube to the corners of a Hamming hypercube. A new classification scheme for comparison of analytical spectra based on their binary encoded spectra is introduced. Quantitative expressions for the effect of binary encoding on general Euclidean distances between spectral points are derived and shown to depend upon their spectral classification. Generally binary encoding increases the interclass distances in pattern recognition and may decrease the intraclass distances. This effect is illustrated with a mass spectral pattern recognition example. The effect of binary encoding on library searches is to produce possible false compound identification in identity searches and to flag spectrally similar compounds in structure searches. A scheme which uses both the Hamming and the Euclidean metrics is proposed for improved library searches. This scheme is illustrated with searches of a small mass spectral library for benzene and p-dioxane spectra.     

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.     

Effect of magnetic field strength on NMR-based metabonomic human urine data. Comparative study of 250, 400, 500, and 800 MHz

Metabonomic analysis of urine utilizing high-resolution NMR spectroscopy and chemometric techniques has proven valuable in characterizing the biochemical response to an intervention. To assess the effect of magnetic field strength on information contained in NMR-based metabonomic data sets, 1H NMR spectra were acquired on 250-, 400-, 500-, and 800-MHz instruments, respectively, on the same set of human urine samples collected before and after dietary interventions with milk and with meat proteins. Partial least-squares regression discriminant analyses (PLS-DA) were performed in order to elucidate the ability of the 1H spectra acquired at various field strengths to identify possible spectral differences and discriminate between pre- and postintervention samples. The loadings from PLS-DA contained the same spectral regions, implying that the same metabolites were involved in the discrimination independent of magnetic field strength. The investigation revealed a strong increase in prediction performance and thereby spectral information content when increasing the magnetic field strength from 250 to 500 MHz, while from 500 to 800 MHz the increase was less pronounced.     

Optimization of human plasma 1H NMR spectroscopic data processing for high-throughput metabolic phenotyping studies and detection of insulin resistance related to type 2 diabetes

Optimizing NMR experimental parameters for high-throughput metabolic phenotyping requires careful examination of the total biochemical information obtainable from (1)H NMR data, which includes concentration and molecular dynamics information. Here we have applied two different types of mathematical transformation (calculation of the first derivative of the NMR spectrum and Gaussian shaping of the free-induction decay) to attenuate broad spectral features from macromolecules and enhance the signals of small molecules. By application of chemometric methods such as principal component analysis (PCA), orthogonal projections to latent structures discriminant analysis (O-PLS-DA) and statistical spectroscopic tools such as statistical total correlation spectroscopy (STOCSY), we show that these methods successfully identify the same potential biomarkers as spin-echo (1)H NMR spectra in which broad lines are suppressed via T2 relaxation editing. Finally, we applied these methods for identification of the metabolic phenotype of patients with type 2 diabetes. This "virtual" relaxation-edited spectroscopy (RESY) approach can be particularly useful for high-throughput screening of complex mixtures such as human plasma and may be useful for extraction of latent biochemical information from legacy or archived NMR data sets for which only standard 1D data sets exist.     

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.     

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.     

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.