Deconvolution of Two-Dimensional NMR Spectra by Fast Maximum Likelihood Reconstruction: Application to Quantitative Metabolomics

We recently developed two databases and a laboratory information system as resources for the metabolomics community. These tools are freely available and are intended to ease data analysis in both MS and NMR based metabolomics studies. The first database is a metabolomics extension to the BioMagResBank (BMRB, http://www.bmrb.wisc.edu), which currently contains experimental spectral data on over 270 pure compounds. Each small molecule entry consists of five or six one- and two-dimensional NMR data sets, along with information about the source of the compound, solution conditions, data collection protocol and the NMR pulse sequences. Users have free access to peak lists, spectra, and original time-domain data. The BMRB database can be queried by name, monoisotopic mass and chemical shift. We are currently developing a deposition tool that will enable people in the community to add their own data to this resource. Our second database, the Madison Metabolomics Consortium Database (MMCD, available from http://mmcd.nmrfam.wisc.edu/), is a hub for information on over 10,000 metabolites. These data were collected from a variety of sites with an emphasis on metabolites found in Arabidopsis. The MMC database supports extensive search functions and allows users to make bulk queries using experimental MS and/or NMR data. In addition to these databases, we have developed a new module for the Sesame laboratory information management system (http://www.sesame.wisc.edu) that captures all of the experimental protocols, background information, and experimental data associated with metabolomics samples. Sesame was designed to help coordinate research efforts in laboratories with high sample throughput and multiple investigators and to track all of the actions that have taken place in a particular study.

A general theory has been developed for the application of the maximum likelihood (ML) principle to the estimation of NMR parameters (frequency and amplitudes) from multidimensional time-domain NMR data. A computer program (ChiFit) has been written that carries out ML parameter estimation in the D-1 indirectly detected dimensions of a D-dimensional NMR data set. The performance of this algorithm has been tested with experimental three-dimensional (HNCO) and four-dimensional (HN(CO)-CAHA) data from a small protein labeled with 13C and 15N. These data sets, with different levels of digital resolution, were processed using ChiFit for ML analysis and employing conventional Fourier transform methods with prior extrapolation of the time-domain dimensions by linear prediction. Comparison of the results indicates that the ML approach provides superior frequency resolution compared to conventional methods, particularly under conditions of limited digital resolution in the time-domain input data, as is characteristic of D-dimensional NMR data of biomolecules. Close correspondence is demonstrated between the results of analyzing multidimensional time-domain NMR data by Fourier transformation, Bayesian probability theory [Chylla, R.A. and Markley, J.L. (1993) J. Biomol. NMR, 3, 515-533], and the ML principle.