Curve-fitting method for direct quantitation of compounds in complex biological mixtures using 1H NMR: application in metabonomic toxicology studies

NMR in Metabolomics: From Conventional Statistics to Machine Learning and Neural Network Approaches

NMR measurements combined with chemometrics allow achieving a great amount of information for the identification of potential biomarkers responsible for a precise metabolic pathway. These kinds of data are useful in different fields, ranging from food to biomedical fields, including health science. The investigation of the whole set of metabolites in a sample, representing its fingerprint in the considered condition, is known as metabolomics and may take advantage of different statistical tools. The new frontier is to adopt self-learning techniques to enhance clustering or classification actions that can improve the predictive power over large amounts of data. Although machine learning is already employed in metabolomics, deep learning and artificial neural networks approaches were only recently successfully applied. In this work, we give an overview of the statistical approaches underlying the wide range of opportunities that machine learning and neural networks allow to perform with accurate metabolites assignment and quantification.Various actual challenges are discussed, such as proper metabolomics, deep learning architectures and model accuracy.

Complete Reduction to Amplitude Frequency Table (CRAFT)-A perspective

The CRAFT (Complete Reduction to Amplitude Frequency Table) technique, based on Bayesian analysis approach, converts FID and/or interferogram (time domain) to a frequency-amplitude table (tabular domain) in a robust, automated, and time-efficient fashion. This mini review/perspective presents an introduction to CRAFT as a processing workflow followed by a discussion of several practical 1D and 2D examples of its applicability and associated benefit. CRAFT provides high quality quantitative results for complex systems without any need for conventional preprocessing steps, such as phase and baseline corrections. Two-dimensional time domain data are typically truncated, particularly in the evolution dimension, and conventional processing after zero-filling and t_1max -matched apodization masks potentially available peak resolution. The line broadening introduced by extensive zero-filling and severe apodization functions leads to the lack of clear resolution of cross peaks. CRAFT decimation of interferograms, on the other hand, requires minimal or no apodization prior to extraction of the NMR parameters and significantly improves the spectral linewidth of the cross peaks along F₁ dimension compared to conventional (FT) processing. The tabular representation of the CRAFT2d cross peaks information can be visualized in a variety of frequency domain formats for conventional spectral interpretation as well as quantitative applications. A simple workflow to generate in silico oversampled interferogram (iSOS) is presented, and its potential benefit in CRAFT decimation of highly crowded 2D NMR is demonstrated. This report is meant as a collective thesis to present a potentially new paradigm in data processing that questions the need for hitherto unchallenged preprocessing steps, such as phase and baseline correction in 1D and zero-fill/severe apodization in 2D.

Practical guide for selection of 1 H qNMR acquisition and processing parameters confirmed by automated spectra evaluation

In our recent paper, a new technique for automated spectra integration and quality control of the acquired results in qNMR was developed and validated (Monakhova & Diehl, Magn. Res. Chem. 2017, doi: 10.1002/mrc.4591). The present study is focused on the influence of acquisition and postacquisition parameters on the developed automated routine in particular, and on the quantitative NMR (qNMR) results in general, which has not been undertaken previously in a systematic and automated manner. Results are presented for a number of model mixtures and authentic pharmaceutical products measured on 500- and 600-MHz NMR spectrometers. The influence of the most important acquisition (spectral width, transmitter [frequency] offset, number of scans, and time domain) and processing (size of real spectrum, deconvolution, Gaussian window multiplication, and line broadening) parameters for qNMR was automatically investigated. Moderate modification of the majority of the investigated parameters from default instrument settings within evaluated ranges does not significantly affect the trueness and precision of the qNMR. Lite Gaussian window multiplication resulted in accuracy improvement of the qNMR output and is recommended for routine measurements. In general, given that the acquisition and processing parameters were selected based on the presented guidelines, automated qNMR analysis can be employed for reproducible high-precision concentration measurements in practice.

Quantitative NMR analysis of intra- and extracellular metabolism of mammalian cells: A tutorial

Metabolomics analysis of body fluids as well as cells is depended on many factors. While several well-accepted standard operating procedures for the analysis of body fluids are available, the NMR based quantitative analysis of cellular metabolites is less well standardized. Experimental designs depend on the cell type, the quenching protocol and the applied post-acquisition workflow. Here, we provide a tutorial for the quantitative description of the metabolic phenotype of mammalian cells using NMR spectroscopy. We discuss all key steps of the process, starting from the selection of the appropriate culture medium, quenching techniques to arrest metabolism in a reproducible manner, the extraction of the intracellular components and the profiling of the culture medium. NMR data acquisition and methods for both qualitative and quantitative analysis are also provided. The suggested methods cover experiments for adherent cells and cells in suspension. We ultimately describe the application of the discussed workflow to a thyroid cancer cell line. Although this tutorial focuses on mammalian cells, the given guidelines and procedures may be adjusted for the analysis of other cell types.

Purity assessment of ginsenoside Rg1 using quantitative 1H nuclear magnetic resonance

Ginseng herbs comprise a group of the most popular herbs, including Panax ginseng, P. notoginseng and P. quinquefolius (Family Araliaceae), which are used as traditional Chinese medicine (TCM) and are some of the best-selling natural products in the world. The accurate quantification of ginsenoside Rg1 is one of the major aspects of its quality control. However, the purity of the commercial Rg1 chemical reference substance (CRS) is often measured with high-performance chromatography coupled with an ultraviolet detector (HPLC-UV), which is a selective detector with unequal responses to different compounds; thus, this detector introduces probable error to purity assessments. In the present study, quantitative nuclear magnetic resonance (qNMR), due to its absolute quantification ability, was applied to accurately assess the purity of Rg1 CRS. Phenylmethyl phthalate was used as the internal standard (IS) to calibrate the purity of Rg1 CRS. The proton signal of Rg1 CRS in methanol-d4 at 4.37ppm was selected to avoid interfering signals, enabling accurate quantitative analysis. The relaxation delay, number of scans, and NMR windowing were optimized for data acquisition. For post-processing, the Lorentz/Gauss deconvolution method was employed to increase the signal accuracy by separating the impurities and noise in the integrated region of the quantitative proton. The method validation showed that the developed method has acceptable sensitivity, linearity, precision, and accuracy. The purity of the commercial Rg1 CRS examined with the method developed in this research was 90.34±0.21%, which was obviously lower than that reported by the manufacturer (>98.0%, HPLC-UV). The cross-method validation shows that the commonly used HPLC-UV, HPLC-ELSD (evaporative light scattering detector) and even LC-MS (mass spectrometry) methods provide significantly higher purity values of Rg1 CRS compared with the qNMR method, and the accuracy of these LC-based methods largely depend on the amount of the sample that was loaded and the properties of the impurities.

Natural Products Diversity of Marine Ascidians (Tunicates; Ascidiacea) and Successful Drugs in Clinical Development

This present study reviewed the chemical diversity of marine ascidians and their pharmacological applications, challenges and recent developments in marine drug discovery reported during 1994-2014, highlighting the structural activity of compounds produced by these specimens. Till date only 5% of living ascidian species were studied from <3000 species, this study represented from family didemnidae (32%), polyclinidae (22%), styelidae and polycitoridae (11-12%) exhibiting the highest number of promising MNPs. Close to 580 compound structures are here discussed in terms of their occurrence, structural type and reported biological activity. Anti-cancer drugs are the main area of interest in the screening of MNPs from ascidians (64%), followed by anti-malarial (6%) and remaining others. FDA approved ascidian compounds mechanism of action along with other compounds status of clinical trials (phase 1 to phase 3) are discussed here in. This review highlights recent developments in the area of natural products chemistry and biotechnological approaches are emphasized.

Precision high-throughput proton NMR spectroscopy of human urine, serum, and plasma for large-scale metabolic phenotyping

Proton nuclear magnetic resonance (NMR)-based metabolic phenotyping of urine and blood plasma/serum samples provides important prognostic and diagnostic information and permits monitoring of disease progression in an objective manner. Much effort has been made in recent years to develop NMR instrumentation and technology to allow the acquisition of data in an effective, reproducible, and high-throughput approach that allows the study of general population samples from epidemiological collections for biomarkers of disease risk. The challenge remains to develop highly reproducible methods and standardized protocols that minimize technical or experimental bias, allowing realistic interlaboratory comparisons of subtle biomarker information. Here we present a detailed set of updated protocols that carefully consider major experimental conditions, including sample preparation, spectrometer parameters, NMR pulse sequences, throughput, reproducibility, quality control, and resolution. These results provide an experimental platform that facilitates NMR spectroscopy usage across different large cohorts of biofluid samples, enabling integration of global metabolic profiling that is a prerequisite for personalized healthcare.

Potential of metabolomics in preclinical and clinical drug development

Metabolomics is an upcoming technology system which involves detailed experimental analysis of metabolic profiles. Due to its diverse applications in preclinical and clinical research, it became an useful tool for the drug discovery and drug development process. This review covers the brief outline about the instrumentation and interpretation of metabolic profiles. The applications of metabolomics have a considerable scope in the pharmaceutical industry, almost at each step from drug discovery to clinical development. These include finding drug target, potential safety and efficacy biomarkers and mechanisms of drug action, the validation of preclinical experimental models against human disease profiles, and the discovery of clinical safety and efficacy biomarkers. As we all know, nowadays the drug discovery and development process is a very expensive, and risky business. Failures at any stage of drug discovery and development process cost millions of dollars to the companies. Some of these failures or the associated risks could be prevented or minimized if there were better ways of drug screening, drug toxicity profiling and monitoring adverse drug reactions. Metabolomics potentially offers an effective route to address all the issues associated with the drug discovery and development.

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.

NMR methods for metabolomics of mammalian cell culture bioreactors

Metabolomics has become an important tool for measuring pools of small molecules in mammalian cell cultures expressing therapeutic proteins. NMR spectroscopy has played an important role, largely because it requires minimal sample preparation, does not require chromatographic separation, and is quantitative. The concentrations of large numbers of small molecules in the extracellular media or within the cells themselves can be measured directly on the culture supernatant and on the supernatant of the lysed cells, respectively, and correlated with endpoints such as titer, cell viability, or glycosylation patterns. The observed changes can be used to generate hypotheses by which these parameters can be optimized. This chapter focuses on the sample preparation, data acquisition, and analysis to get the most out of NMR metabolomics data from CHO cell cultures but could easily be extended to other in vitro culture systems.

Reproducibility of NMR analysis of urine samples: impact of sample preparation, storage conditions, and animal health status

Introduction. Spectroscopic analysis of urine samples from laboratory animals can be used to predict the efficacy and side effects of drugs. This employs methods combining ¹H NMR spectroscopy with quantification of biomarkers or with multivariate data analysis. The most critical steps in data evaluation are analytical reproducibility of NMR data (collection, storage, and processing) and the health status of the animals, which may influence urine pH and osmolarity. Methods. We treated rats with a solvent, a diuretic, or a nephrotoxicant and collected urine samples. Samples were titrated to pH 3 to 9, or salt concentrations increased up to 20-fold. The effects of storage conditions and freeze-thaw cycles were monitored. Selected metabolites and multivariate data analysis were evaluated after ¹H NMR spectroscopy. Results. We showed that variation of pH from 3 to 9 and increases in osmolarity up to 6-fold had no effect on the quantification of the metabolites or on multivariate data analysis. Storage led to changes after 14 days at 4°C or after 12 months at −20°C, independent of sample composition. Multiple freeze-thaw cycles did not affect data analysis. Conclusion. Reproducibility of NMR measurements is not dependent on sample composition under physiological or pathological conditions.

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

Metabolomics techniques in nanotoxicology studies

The rapid growth in the development of nanoparticles for uses in a variety of applications including targeted drug delivery, cancer therapy, imaging, and as biological sensors has led to questions about potential toxicity of such particles to humans. High-throughput methods are necessary to evaluate the potential toxicity of nanoparticles. The omics technologies are particularly well suited to evaluate toxicity in both in vitro and in vivo systems. Metabolomics, specifically, can rapidly screen for biomarkers related to predefined pathways or processes in biofluids and tissues. Specifically, oxidative stress has been implicated as a potential mechanism of toxicity in nanoparticles and is generally difficult to measure by conventional methods. Furthermore, metabolomics can provide mechanistic insight into nanotoxicity. This chapter focuses on the application of both LC/MS and NMR-based metabolomics approaches to study the potential toxicity of nanoparticles.

MetaboHunter: an automatic approach for identification of metabolites from 1H-NMR spectra of complex mixtures

BACKGROUND

One-dimensional 1H-NMR spectroscopy is widely used for high-throughput characterization of metabolites in complex biological mixtures. However, the accurate identification of individual compounds is still a challenging task, particularly in spectral regions with higher peak densities. The need for automatic tools to facilitate and further improve the accuracy of such tasks, while using increasingly larger reference spectral libraries becomes a priority of current metabolomics research.

RESULTS

We introduce a web server application, called MetaboHunter, which can be used for automatic assignment of 1H-NMR spectra of metabolites. MetaboHunter provides methods for automatic metabolite identification based on spectra or peak lists with three different search methods and with possibility for peak drift in a user defined spectral range. The assignment is performed using as reference libraries manually curated data from two major publicly available databases of NMR metabolite standard measurements (HMDB and MMCD). Tests using a variety of synthetic and experimental spectra of single and multi metabolite mixtures show that MetaboHunter is able to identify, in average, more than 80% of detectable metabolites from spectra of synthetic mixtures and more than 50% from spectra corresponding to experimental mixtures. This work also suggests that better scoring functions improve by more than 30% the performance of MetaboHunter's metabolite identification methods.

CONCLUSIONS

MetaboHunter is a freely accessible, easy to use and user friendly 1H-NMR-based web server application that provides efficient data input and pre-processing, flexible parameter settings, fast and automatic metabolite fingerprinting and results visualization via intuitive plotting and compound peak hit maps. Compared to other published and freely accessible metabolomics tools, MetaboHunter implements three efficient methods to search for metabolites in manually curated data from two reference libraries.

Sparse non-negative generalized PCA with applications to metabolomics

MOTIVATION

Nuclear magnetic resonance (NMR) spectroscopy has been used to study mixtures of metabolites in biological samples. This technology produces a spectrum for each sample depicting the chemical shifts at which an unknown number of latent metabolites resonate. The interpretation of this data with common multivariate exploratory methods such as principal components analysis (PCA) is limited due to high-dimensionality, non-negativity of the underlying spectra and dependencies at adjacent chemical shifts.

RESULTS

We develop a novel modification of PCA that is appropriate for analysis of NMR data, entitled Sparse Non-Negative Generalized PCA. This method yields interpretable principal components and loading vectors that select important features and directly account for both the non-negativity of the underlying spectra and dependencies at adjacent chemical shifts. Through the reanalysis of experimental NMR data on five purified neural cell types, we demonstrate the utility of our methods for dimension reduction, pattern recognition, sample exploration and feature selection. Our methods lead to the identification of novel metabolites that reflect the differences between these cell types.

AVAILABILITY

www.stat.rice.edu/~gallen/software.html.

CONTACT

gallen@rice.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification and quantification of metabolites in (1)H NMR spectra by Bayesian model selection

MOTIVATION

Nuclear magnetic resonance (NMR) spectroscopy is widely used for high-throughput characterization of metabolites in complex biological mixtures. However, accurate interpretation of the spectra in terms of identities and abundances of metabolites can be challenging, in particular in crowded regions with heavy peak overlap. Although a number of computational approaches for this task have recently been proposed, they are not entirely satisfactory in either accuracy or extent of automation.

RESULTS

We introduce a probabilistic approach Bayesian Quantification (BQuant), for fully automated database-based identification and quantification of metabolites in local regions of (1)H NMR spectra. The approach represents the spectra as mixtures of reference profiles from a database, and infers the identities and the abundances of metabolites by Bayesian model selection. We show using a simulated dataset, a spike-in experiment and a metabolomic investigation of plasma samples that BQuant outperforms the available automated alternatives in accuracy for both identification and quantification.

AVAILABILITY

The R package BQuant is available at: http://www.stat.purdue.edu/~ovitek/BQuant-Web/.

Application of a Bayesian deconvolution approach for high-resolution (1)H NMR spectra to assessing the metabolic effects of acute phenobarbital exposure in liver tissue

High-resolution (1)H NMR spectroscopy is frequently used in the field of metabolomics to assess the metabolites found in biofluids or tissue extracts to define a metabolic profile that describes a given biological process. In this study, we aimed to increase the utility of NMR-based metabolomics by using advanced Bayesian modeling of the time-domain high-resolution 1D NMR free induction decay (FID). The improvement over traditional nonparametric binning is twofold and associated with enhanced resolution of the analysis and automation of the signal processing stage. The automation is achieved by using a Bayesian formalism for all parameters of the model including the number of components. The approach is illustrated with a study of early markers of acute exposure to different doses of a well-characterized nongenotoxic hepatocarcinogen, phenobarbital, in rats. The results demonstrate that Bayesian deconvolution produces a better model for the NMR spectra that allows the identification of subtle changes in metabolic concentrations and a decrease in the expected false discovery rate compared with approaches based on "binning". These properties suggest that Bayesian deconvolution could facilitate the biomarker discovery process and improve information extraction from high-resolution NMR spectra.

NMR-based metabolic profiling identifies biomarkers of liver regeneration following partial hepatectomy in the rat

Tissue injury and repair are often overlapping consequences of disease or toxic exposure, but are not often considered as distinct processes in molecular studies. To establish the systemic metabolic response to liver regeneration, the partial hepatectomy (PH) model has been studied in the rat by an integrated metabonomics strategy, utilizing (1)H NMR spectroscopy of urine, liver and serum. Male Sprague-Dawley rats were subjected to either surgical removal of approximately two-thirds of the liver, sham operated (SO) surgery, or no treatment (n = 10/group) and samples collected over a 7 day period. A number of urinary metabolic perturbations were observed in PH rats compared with SO and control animals, including elevated levels of taurine, hypotaurine, creatine, guanidinoacetic acid, betaine, dimethylglycine and bile acids. Serum betaine and creatine were also elevated after PH, while levels of triglyceride were reduced. In the liver, triglycerides, cholesterol, alanine and betaine were elevated after PH, while choline and its derivatives were reduced. Upon examining the dynamic pattern of urinary response (the 'metabolic trajectory'), several metabolites could be categorized into groups likely to reflect perturbations to different processes such as dietary intake or hepatic 1-carbon metabolism. Several of the urinary perturbations observed during the regenerative phase of the PH model have also been observed after exposure to liver toxins, indicating that hepatic regeneration may make a contribution to the systemic alterations in metabolism associated with hepatotoxicity. The observed changes in 1-carbon and lipid metabolism are consistent with the proposed role of these pathways in the activation of a regenerative response and provide further evidence regarding the utility of urinary NMR profiles in the detection of liver-specific pathology. Biofluid (1)H NMR-based metabolic profiling provides new insight into the role of metabolism of liver regeneration, and suggests putative biomarkers for the noninvasive monitoring of the regeneration process.

Opening up the "Black Box": metabolic phenotyping and metabolome-wide association studies in epidemiology

BACKGROUND

Metabolic phenotyping of humans allows information to be captured on the interactions between dietary, xenobiotic, other lifestyle and environmental exposures, and genetic variation, which together influence the balance between health and disease risks at both individual and population levels.

OBJECTIVES

We describe here the main procedures in large-scale metabolic phenotyping and their application to metabolome-wide association (MWA) studies.

METHODS

By use of high-throughput technologies and advanced spectroscopic methods, application of metabolic profiling to large-scale epidemiologic sample collections, including metabolome-wide association (MWA) studies for biomarker discovery and identification.

DISCUSSION

Metabolic profiling at epidemiologic scale requires optimization of experimental protocol to maximize reproducibility, sensitivity, and quantitative reliability, and to reduce analytical drift. Customized multivariate statistical modeling approaches are needed for effective data visualization and biomarker discovery with control for false-positive associations since 100s or 1,000s of complex metabolic spectra are being processed.

CONCLUSION

Metabolic profiling is an exciting addition to the armamentarium of the epidemiologist for the discovery of new disease-risk biomarkers and diagnostics, and to provide novel insights into etiology, biological mechanisms, and pathways.

Effects of the Yangtze River source of drinking water on metabolites of Mus musculus

The effects of the Yangtze River source of drinking water on metabolites of mouse (Mus musculus) were implemented to observe the environmental health issue of the water by use of the 1H nuclear magnetic resonance (NMR) spectroscopy-based metabonomics. All the sampled mice were treated for 90 days. There were 20 organic pollutants discovered in the water with total concentration of 9.41 microg/l. The NMR spectra for the sampled mice were different at delta2.06, delta2.24 and delta5.22. The concentrations of alanine, glycoprotein, acetone and trimethylamine-N-oxide in the source water group mice were decreased but that of creatinine and glucose were increased, which indicated that hepatotoxicity and kidney dysfunction occurred. There were six parameters for the source water group mice were different or extremely different from that of the control group. And these metabolites are responsible for separation of the data along PC1 and PC2 which may be used as biomarkers to indicate the source water pollution. The results indicate that 1H NMR-based metabonomic approach is a useful technique to test toxicity of xenobiotics on metabolites for observation of the environmental health issue of source water.

Metabolic profiling and population screening of analgesic usage in nuclear magnetic resonance spectroscopy-based large-scale epidemiologic studies

The application of a (1)H nuclear magnetic resonance (NMR) spectroscopy-based screening method for determining the use of two widely available analgesics (acetaminophen and ibuprofen) in epidemiologic studies has been investigated. We used samples and data from the cross-sectional INTERMAP Study involving participants from Japan (n = 1145), China (n = 839), U.K. (n = 501), and the U.S. (n = 2195). An orthogonal projection to latent structures discriminant analysis (OPLS-DA) algorithm with an incorporated Monte Carlo resampling function was applied to the NMR data set to determine which spectra contained analgesic metabolites. OPLS-DA preprocessing parameters (normalization, bin width, scaling, and input parameters) were assessed systematically to identify an optimal acetaminophen prediction model. Subsets of INTERMAP spectra were examined to verify and validate the presence/absence of acetaminophen/ibuprofen based on known chemical shift and coupling patterns. The optimized and validated acetaminophen model correctly predicted 98.2%, and the ibuprofen model correctly predicted 99.0% of the urine specimens containing these drug metabolites. The acetaminophen and ibuprofen models were subsequently used to predict the presence/absence of these drug metabolites for the remaining INTERMAP specimens. The acetaminophen model identified 415 out of 8436 spectra as containing acetaminophen metabolite signals while the ibuprofen model identified 245 out of 8604 spectra as containing ibuprofen metabolite signals from the global data set after excluding samples used to construct the prediction models. The NMR-based metabolic screening strategy provides a new objective approach for evaluation of self-reported medication data and is extendable to other aspects of population xenometabolome profiling.

MetaboMiner--semi-automated identification of metabolites from 2D NMR spectra of complex biofluids

Background

One-dimensional (1D) ¹H nuclear magnetic resonance (NMR) spectroscopy is widely used in metabolomic studies involving biofluids and tissue extracts. There are several software packages that support compound identification and quantification via 1D ¹H NMR by spectral fitting techniques. Because 1D ¹H NMR spectra are characterized by extensive peak overlap or spectral congestion, two-dimensional (2D) NMR, with its increased spectral resolution, could potentially improve and even automate compound identification or quantification. However, the lack of dedicated software for this purpose significantly restricts the application of 2D NMR methods to most metabolomic studies.

Results

We describe a standalone graphics software tool, called MetaboMiner, which can be used to automatically or semi-automatically identify metabolites in complex biofluids from 2D NMR spectra. MetaboMiner is able to handle both ¹H-¹H total correlation spectroscopy (TOCSY) and ¹H-¹³C heteronuclear single quantum correlation (HSQC) data. It identifies compounds by comparing 2D spectral patterns in the NMR spectrum of the biofluid mixture with specially constructed libraries containing reference spectra of ~500 pure compounds. Tests using a variety of synthetic and real spectra of compound mixtures showed that MetaboMiner is able to identify >80% of detectable metabolites from good quality NMR spectra.

Conclusion

MetaboMiner is a freely available, easy-to-use, NMR-based metabolomics tool that facilitates automatic peak processing, rapid compound identification, and facile spectrum annotation from either 2D TOCSY or HSQC spectra. Using comprehensive reference libraries coupled with robust algorithms for peak matching and compound identification, the program greatly simplifies the process of metabolite identification in complex 2D NMR spectra.

Metabolic profiling of an Echinostoma caproni infection in the mouse for biomarker discovery

BACKGROUND

Metabolic profiling holds promise with regard to deepening our understanding of infection biology and disease states. The objectives of our study were to assess the global metabolic responses to an Echinostoma caproni infection in the mouse, and to compare the biomarkers extracted from different biofluids (plasma, stool, and urine) in terms of characterizing acute and chronic stages of this intestinal fluke infection.

METHODOLOGY/PRINCIPAL FINDINGS

Twelve female NMRI mice were infected with 30 E. caproni metacercariae each. Plasma, stool, and urine samples were collected at 7 time points up to day 33 post-infection. Samples were also obtained from non-infected control mice at the same time points and measured using (1)H nuclear magnetic resonance (NMR) spectroscopy. Spectral data were subjected to multivariate statistical analyses. In plasma and urine, an altered metabolic profile was already evident 1 day post-infection, characterized by reduced levels of plasma choline, acetate, formate, and lactate, coupled with increased levels of plasma glucose, and relatively lower concentrations of urinary creatine. The main changes in the urine metabolic profile started at day 8 post-infection, characterized by increased relative concentrations of trimethylamine and phenylacetylglycine and lower levels of 2-ketoisocaproate and showed differentiation over the course of the infection.

CONCLUSION/SIGNIFICANCE

The current investigation is part of a broader NMR-based metabonomics profiling strategy and confirms the utility of this approach for biomarker discovery. In the case of E. caproni, a diagnosis based on all three biofluids would deliver the most comprehensive fingerprint of an infection. For practical purposes, however, future diagnosis might aim at a single biofluid, in which case urine would be chosen for further investigation, based on quantity of biomarkers, ease of sampling, and the degree of differentiation from the non-infected control group.

Human metabolic phenotype diversity and its association with diet and blood pressure

Metabolic phenotypes are the products of interactions among a variety of factors-dietary, other lifestyle/environmental, gut microbial and genetic. We use a large-scale exploratory analytical approach to investigate metabolic phenotype variation across and within four human populations, based on 1H NMR spectroscopy. Metabolites discriminating across populations are then linked to data for individuals on blood pressure, a major risk factor for coronary heart disease and stroke (leading causes of mortality worldwide). We analyse spectra from two 24-hour urine specimens for each of 4,630 participants from the INTERMAP epidemiological study, involving 17 population samples aged 40-59 in China, Japan, UK and USA. We show that urinary metabolite excretion patterns for East Asian and western population samples, with contrasting diets, diet-related major risk factors, and coronary heart disease/stroke rates, are significantly differentiated (P < 10(-16)), as are Chinese/Japanese metabolic phenotypes, and subgroups with differences in dietary vegetable/animal protein and blood pressure. Among discriminatory metabolites, we quantify four and show association (P < 0.05 to P < 0.0001) of mean 24-hour urinary formate excretion with blood pressure in multiple regression analyses for individuals. Mean 24-hour urinary excretion of alanine (direct) and hippurate (inverse), reflecting diet and gut microbial activities, are also associated with blood pressure of individuals. Metabolic phenotyping applied to high-quality epidemiological data offers the potential to develop an area of aetiopathogenetic knowledge involving discovery of novel biomarkers related to cardiovascular disease risk.

NMR-based characterization of metabolic alterations in hypertension using an adaptive, intelligent binning algorithm

As with every -omics technology, metabolomics requires new methodologies for data processing. Due to the large spectral size, a standard approach in NMR-based metabolomics implies the division of spectra into equally sized bins, thereby simplifying subsequent data analysis. Yet, disadvantages are the loss of information and the occurrence of artifacts caused by peak shifts. Here, a new binning algorithm, Adaptive Intelligent Binning (AI-Binning), which largely circumvents these problems, is presented. AI-Binning recursively identifies bin edges in existing bins, requires only minimal user input, and avoids the use of arbitrary parameters or reference spectra. The performance of AI-Binning is demonstrated using serum spectra from 40 hypertensive and 40 matched normotensive subjects from the Asklepios study. Hypertension is a major cardiovascular risk factor characterized by a complex biochemistry and, in most cases, an unknown origin. The binning algorithm resulted in an improved classification of hypertensive status compared with that of standard binning and facilitated the identification of relevant metabolites. Moreover, since the occurrence of noise variables is largely avoided, AI-Binned spectra can be unit-variance scaled. This enables the detection of relevant, low-intensity metabolites. These results demonstrate the power of AI-Binning and suggest the involvement of alpha-1 acid glycoproteins and choline biochemistry in hypertension.

Improved identification of metabolites in complex mixtures using HSQC NMR spectroscopy

The automated and robust identification of metabolites in a complex biological sample remains one of the greatest challenges in metabolomics. In our experiments, HSQC carbon-proton correlation NMR data with a model that takes intensity information into account improves upon the identification of metabolites that was achieved using COSY proton-proton correlation NMR data with the binary model of [Y. Xi, J.S. de Ropp, M.R. Viant, D.L. Woodruff, P. Yu, Metabolomics, 2 (2006) 221-233]. In addition, using intensity information results in easier-to-interpret "grey areas" for cases where it is not clear if the compound might be present. We report on highly successful experiments that identify compounds in chemically defined mixtures as well as in biological samples, and compare our two-dimensional HSQC analyses against quantification of metabolites in the corresponding one-dimensional proton NMR spectra. We show that our approach successfully employs a fully automated algorithm for identifying the presence or absence of predefined compounds (held within a library) in biological HSQC spectra, and in addition calculates upper bounds on the compound intensities.

Integrative top-down system metabolic modeling in experimental disease states via data-driven Bayesian methods

Multivariate metabolic profiles from biofluids such as urine and plasma are highly indicative of the biological fitness of complex organisms and can be captured analytically in order to derive top-down systems biology models. The application of currently available modeling approaches to human and animal metabolic pathway modeling is problematic because of multicompartmental cellular and tissue exchange of metabolites operating on many time scales. Hence, novel approaches are needed to analyze metabolic data obtained using minimally invasive sampling methods in order to reconstruct the patho-physiological modulations of metabolic interactions that are representative of whole system dynamics. Here, we show that spectroscopically derived metabolic data in experimental liver injury studies (induced by hydrazine and alpha-napthylisothiocyanate treatment) can be used to derive insightful probabilistic graphical models of metabolite dependencies, which we refer to as metabolic interactome maps. Using these, system level mechanistic information on homeostasis can be inferred, and the degree of reversibility of induced lesions can be related to variations in the metabolic network patterns. This approach has wider application in assessment of system level dysfunction in animal or human studies from noninvasive measurements.

Species variation in the fecal metabolome gives insight into differential gastrointestinal function

The metabolic composition of fecal extracts provides a window for elucidating the complex metabolic interplay between mammals and their intestinal ecosystems, and these metabolite profiles can yield information on a range of gut diseases. Here, the metabolites present in aqueous fecal extracts of humans, mice and rats were characterized using high-resolution (1)H NMR spectroscopy coupled with multivariate pattern recognition techniques. Additionally, the effects of sample storage and preparation methods were evaluated in order to assess the stability of fecal metabolite profiles, and to optimize information recovery from fecal samples. Finally, variations in metabolite profiles were investigated in healthy mice as a function of time. Interspecies variation was found to be greater than the variation due to either time or sample preparation. Although many fecal metabolites were common to the three species, such as short chain fatty acids and branched chain amino acids, each species generated a unique profile. Relatively higher levels of uracil, hypoxanthine, phenylacetic acid, glucose, glycine, and tyrosine amino acids were present in the rat, with beta-alanine being unique to the rat, and glycerol and malonate being unique to the human. Human fecal extracts showed a greater interindividual variation than the two rodent species, reflecting the natural genetic and environmental diversity in human populations. Fecal composition in healthy mice was found to change over time, which might be explained by altered gut microbial presence or activity. The systematic characterization of fecal composition across humans, mice, and rats, together with the evaluation of inherent variation, provides a benchmark for future studies seeking to determine fecal biomarkers of disease and/or response to dietary or therapeutic interventions.

High-precision frequency measurements: indispensable tools at the core of the molecular-level analysis of complex systems

This perspective article provides an assessment of the state-of-the-art in the molecular-resolution analysis of complex organic materials. These materials can be divided into biomolecules in complex mixtures (which are amenable to successful separation into unambiguously defined molecular fractions) and complex nonrepetitive materials (which cannot be purified in the conventional sense because they are even more intricate). Molecular-level analyses of these complex systems critically depend on the integrated use of high-performance separation, high-resolution organic structural spectroscopy and mathematical data treatment. At present, only high-precision frequency-derived data exhibit sufficient resolution to overcome the otherwise common and detrimental effects of intrinsic averaging, which deteriorate spectral resolution to the degree of bulk-level rather than molecular-resolution analysis. High-precision frequency measurements are integral to the two most influential organic structural spectroscopic methods for the investigation of complex materials-NMR spectroscopy (which provides unsurpassed detail on close-range molecular order) and FTICR mass spectrometry (which provides unrivalled resolution)-and they can be translated into isotope-specific molecular-resolution data of unprecedented significance and richness. The quality of this standalone de novo molecular-level resolution data is of unparalleled mechanistic relevance and is sufficient to fundamentally advance our understanding of the structures and functions of complex biomolecular mixtures and nonrepetitive complex materials, such as natural organic matter (NOM), aerosols, and soil, plant and microbial extracts, all of which are currently poorly amenable to meaningful target analysis. The discrete analytical volumetric pixel space that is presently available to describe complex systems (defined by NMR, FT mass spectrometry and separation technologies) is in the range of 10(8-14) voxels, and is therefore capable of providing the necessary detail for a meaningful molecular-level analysis of very complex mixtures. Nonrepetitive complex materials exhibit mass spectral signatures in which the signal intensity often follows the number of chemically feasible isomers. This suggests that even the most strongly resolved FTICR mass spectra of complex materials represent simplified (e.g. isomer-filtered) projections of structural space.

1H NMR metabonomics of plasma lipoprotein subclasses: elucidation of metabolic clustering by self-organising maps

(1)H NMR spectra of plasma are known to provide specific information on lipoprotein subclasses in the form of complex overlapping resonances. A combination of (1)H NMR and self-organising map (SOM) analysis was applied to investigate if automated characterisation of subclass-related metabolic interactions can be achieved. To reliably assess the intrinsic capability of (1)H NMR for resolving lipoprotein subclass profiles, sum spectra representing the pure lipoprotein subclass part of actual plasma were simulated with the aid of experimentally derived model signals for 11 distinct lipoprotein subclasses. Two biochemically characteristic categories of spectra, representing normolipidaemic and metabolic syndrome status, were generated with corresponding lipoprotein subclass profiles. A set of spectra representing a metabolic pathway between the two categories was also generated. The SOM analysis, based solely on the aliphatic resonances of these simulated spectra, clearly revealed the lipoprotein subclass profiles and their changes. Comparable SOM analysis in a group of 69 experimental (1)H NMR spectra of serum samples, which according to biochemical analyses represented a wide range of lipoprotein lipid concentrations, corroborated the findings based on the simulated data. Interestingly, the choline-N(CH(3))(3) region seems to provide more resolved clustering of lipoprotein subclasses in the SOM analyses than the methyl-CH(3) region commonly used for subclass quantification. The results illustrate the inherent suitability of (1)H NMR metabonomics for automated studies of lipoprotein subclass-related metabolism and demonstrate the power of SOM analysis in an extensive and representative case of (1)H NMR metabonomics.

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.

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.

Analytical aspects of plant metabolite profiling platforms: current standings and future aims

Over the past years, metabolic profiling has been established as a comprehensive systems biology tool. Mass spectrometry or NMR spectroscopy-based technology platforms combined with unsupervised or supervised multivariate statistical methodologies allow a deep insight into the complex metabolite patterns of plant-derived samples. Within this review, we provide a thorough introduction to the analytical hard- and software requirements of metabolic profiling platforms. Methodological limitations are addressed, and the metabolic profiling workflow is exemplified by summarizing recent applications ranging from model systems to more applied topics.

Weighted least-squares deconvolution method for discovery of group differences between complex biofluid 1H NMR spectra

Biomarker discovery through analysis of high-throughput NMR data is a challenging, time-consuming process due to the requirement of sophisticated, dataset specific preprocessing techniques and the inherent complexity of the data. Here, we demonstrate the use of weighted, constrained least-squares for fitting a linear mixture of reference standard data to complex urine NMR spectra as an automated way of utilizing current assignment knowledge and the ability to deconvolve confounded spectral regions. Following the least-squares fit, univariate statistics were used to identify metabolites associated with group differences. This method was evaluated through applications on simulated datasets and a murine diabetes dataset. Furthermore, we examined the differential ability of various weighting metrics to correctly identify discriminative markers. Our findings suggest that the weighted least-squares approach is effective for identifying biochemical discriminators of varying physiological states. Additionally, the superiority of specific weighting metrics is demonstrated in particular datasets. An additional strength of this methodology is the ability for individual investigators to couple this analysis with laboratory specific preprocessing techniques.

Metabolic profiling technologies for biomarker discovery in biomedicine and drug development

The state-of-the-art of nuclear magnetic resonance spectroscopy, mass spectrometry and statistical tools for the acquisition and evaluation of complex multidimensional spectroscopic data in metabolic profiling is reviewed in this article. The continuous evolution of the sensitivity, precision and throughput has made these technologies powerful and extremely robust tools for application in systems biology, pharmaceutical and diagnostics research. Particular emphasis is also given to the collection and storage of biological samples that are subjected to metabolite profiling. Selected examples from preclinical and clinical applications are paradigmatically shown. These illustrate the power of the profiling technologies for characterizing the metabolic phenotype of healthy, diseased and treated subjects. The complexity of disease and drug treatment is asking for an adequate response by integrated and comprehensive metabolite profiling approaches that allow the discovery of new combinations of metabolic biomarkers.