Towards high-throughput metabolomics using ultrahigh-field Fourier transform ion cyclotron resonance mass spectrometry

With unmatched mass resolution, mass accuracy, and exceptional detection sensitivity, Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FTICR-MS) has the potential to be a powerful new technique for high-throughput metabolomic analysis. In this study, we examine the properties of an ultrahigh-field 12-Tesla (12T) FTICR-MS for the identification and absolute quantitation of human plasma metabolites, and for the untargeted metabolic fingerprinting of inbred-strain mouse serum by direct infusion (DI). Using internal mass calibration (mass error ≤1 ppm), we determined the rational elemental compositions (incorporating unlimited C, H, N and O, and a maximum of two S, three P, two Na, and one K per formula) of approximately 250 out of 570 metabolite features detected in a 3-min infusion analysis of aqueous extract of human plasma, and were able to identify more than 100 metabolites. Using isotopically-labeled internal standards, we were able to obtain excellent calibration curves for the absolute quantitation of choline with sub-pmol sensitivity, using 500 times less sample than previous LC/MS analyses. Under optimized serum dilution conditions, chemical compounds spiked into mouse serum as metabolite mimics showed a linear response over a 600-fold concentration range. DI/FTICR-MS analysis of serum from 26 mice from 2 inbred strains, with and without acute trichloroethylene (TCE) treatment, gave a relative standard deviation (RSD) of 4.5%. Finally, we extended this method to the metabolomic fingerprinting of serum samples from 49 mice from 5 inbred strains involved in an acute alcohol toxicity study, using both positive and negative electrospray ionization (ESI). Using these samples, we demonstrated the utility of this method for high-throughput metabolomics, with more than 400 metabolites profiled in only 24 h. Our experiments demonstrate that DI/FTICR-MS is well-suited for high-throughput metabolomic analysis. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ab initio prediction of metabolic networks using Fourier transform mass spectrometry data

Fourier transform mass spectrometry has recently been introduced into the field of metabolomics as a technique that enables the mass separation of complex mixtures at very high resolution and with ultra high mass accuracy. Here we show that this enhanced mass accuracy can be exploited to predict large metabolic networks ab initio, based only on the observed metabolites without recourse to predictions based on the literature. The resulting networks are highly information-rich and clearly non-random. They can be used to infer the chemical identity of metabolites and to obtain a global picture of the structure of cellular metabolic networks. This represents the first reconstruction of metabolic networks based on unbiased metabolomic data and offers a breakthrough in the systems-wide analysis of cellular metabolism.     

Description and analysis of metabolic connectivity and dynamics in the human red blood cell

The human red blood cell (hRBC) metabolic network is relatively simple compared with other whole cell metabolic networks, yet too complicated to study without the aid of a computer model. Systems science techniques can be used to uncover the key dynamic features of hRBC metabolism. Herein, we have studied a full dynamic hRBC metabolic model and developed several approaches to identify metabolic pools of metabolites. In particular, we have used phase planes, temporal decomposition, and statistical analysis to show hRBC metabolism is characterized by the formation of pseudoequilibrium concentration states. Such equilibria identify metabolic "pools" or aggregates of concentration variables. We proceed to define physiologically meaningful pools, characterize them within the hRBC, and compare them with those derived from systems engineering techniques. In conclusion, systems science methods can decipher detailed information about individual enzymes and metabolites within metabolic networks and provide further understanding of complex biological networks.     

Use of reconstituted metabolic networks to assist in metabolomic data visualization and mining

Metabolomics experiments seldom achieve their aim of comprehensively covering the entire metabolome. However, important information can be gleaned even from sparse datasets, which can be facilitated by placing the results within the context of known metabolic networks. Here we present a method that allows the automatic assignment of identified metabolites to positions within known metabolic networks, and, furthermore, allows automated extraction of sub-networks of biological significance. This latter feature is possible by use of a gap-filling algorithm. The utility of the algorithm in reconstructing and mining of metabolomics data is shown on two independent datasets generated with LC–MS LTQ-Orbitrap mass spectrometry. Biologically relevant metabolic sub-networks were extracted from both datasets. Moreover, a number of metabolites, whose presence eluded automatic selection within mass spectra, could be identified retrospectively by virtue of their inferred presence through gap filling.     

NMR-based stable isotope resolved metabolomics in systems biochemistry

An important goal of metabolomics is to characterize the changes in metabolic networks in cells or various tissues of an organism in response to external perturbations or pathologies. The profiling of metabolites and their steady state concentrations does not directly provide information regarding the architecture and fluxes through metabolic networks. This requires tracer approaches. NMR is especially powerful as it can be used not only to identify and quantify metabolites in an unfractionated mixture such as biofluids or crude cell/tissue extracts, but also determine the positional isotopomer distributions of metabolites derived from a precursor enriched in stable isotopes such as (13)C and (15)N via metabolic transformations. In this article we demonstrate the application of a variety of 2-D NMR editing experiments to define the positional isotopomers of compounds present in polar and non-polar extracts of human lung cancer cells grown in either [U-(13)C]-glucose or [U-(13)C,(15)N]-glutamine as source tracers. The information provided by such experiments enabled unambiguous reconstruction of metabolic pathways, which is the foundation for further metabolic flux modeling.     

Biosynthetic potentials of metabolites and their hierarchical organization

A major challenge in systems biology is to understand how complex and highly connected metabolic networks are organized. The structure of these networks is investigated here by identifying sets of metabolites that have a similar biosynthetic potential. We measure the biosynthetic potential of a particular compound by determining all metabolites than can be produced from it and, following a terminology introduced previously, call this set the scope of the compound. To identify groups of compounds with similar scopes, we apply a hierarchical clustering method. We find that compounds within the same cluster often display similar chemical structures and appear in the same metabolic pathway. For each cluster we define a consensus scope by determining a set of metabolites that is most similar to all scopes within the cluster. This allows for a generalization from scopes of single compounds to scopes of a chemical family. We observe that most of the resulting consensus scopes overlap or are fully contained in others, revealing a hierarchical ordering of metabolites according to their biosynthetic potential. Our investigations show that this hierarchy is not only determined by the chemical complexity of the metabolites, but also strongly by their biological function. As a general tendency, metabolites which are necessary for essential cellular processes exhibit a larger biosynthetic potential than those involved in secondary metabolism. A central result is that chemically very similar substances with different biological functions may differ significantly in their biosynthetic potentials. Our studies provide an important step towards understanding fundamental design principles of metabolic networks determined by the structural and functional complexity of metabolites.     

Stability of Metabolic Correlations under Changing Environmental Conditions in Escherichia coli – A Systems Approach

Background

Biological systems adapt to changing environments by reorganizing their cellular and physiological program with metabolites representing one important response level. Different stresses lead to both conserved and specific responses on the metabolite level which should be reflected in the underlying metabolic network.

Methodology/Principal Findings

Starting from experimental data obtained by a GC-MS based high-throughput metabolic profiling technology we here develop an approach that: (1) extracts network representations from metabolic condition-dependent data by using pairwise correlations, (2) determines the sets of stable and condition-dependent correlations based on a combination of statistical significance and homogeneity tests, and (3) can identify metabolites related to the stress response, which goes beyond simple observations about the changes of metabolic concentrations. The approach was tested with Escherichia coli as a model organism observed under four different environmental stress conditions (cold stress, heat stress, oxidative stress, lactose diauxie) and control unperturbed conditions. By constructing the stable network component, which displays a scale free topology and small-world characteristics, we demonstrated that: (1) metabolite hubs in this reconstructed correlation networks are significantly enriched for those contained in biochemical networks such as EcoCyc, (2) particular components of the stable network are enriched for functionally related biochemical pathways, and (3) independently of the response scale, based on their importance in the reorganization of the correlation network a set of metabolites can be identified which represent hypothetical candidates for adjusting to a stress-specific response.

Conclusions/Significance

Network-based tools allowed the identification of stress-dependent and general metabolic correlation networks. This correlation-network-based approach does not rely on major changes in concentration to identify metabolites important for stress adaptation, but rather on the changes in network properties with respect to metabolites. This should represent a useful complementary technique in addition to more classical approaches.     

The connectivity structure,giant strong component and centrality of metabolic networks

MOTIVATION: Structural and functional analysis of genome-based large-scale metabolic networks is important for understanding the design principles and regulation of the metabolism at a system level. The metabolic network is conventionally considered to be highly integrated and very complex. A rational reduction of the metabolic network to its core structure and a deeper understanding of its functional modules are important. RESULTS: In this work, we show that the metabolites in a metabolic network are far from fully connected. A connectivity structure consisting of four major subsets of metabolites and reactions, i.e. a fully connected sub-network, a substrate subset, a product subset and an isolated subset is found to exist in metabolic networks of 65 fully sequenced organisms. The largest fully connected part of a metabolic network, called 'the giant strong component (GSC)', represents the most complicated part and the core of the network and has the feature of scale-free networks. The average path length of the whole network is primarily determined by that of the GSC. For most of the organisms, GSC normally contains less than one-third of the nodes of the network. This connectivity structure is very similar to the 'bow-tie' structure of World Wide Web. Our results indicate that the bow-tie structure may be common for large-scale directed networks. More importantly, the uncovered structure feature makes a structural and functional analysis of large-scale metabolic network more amenable. As shown in this work, comparing the closeness centrality of the nodes in the GSC can identify the most central metabolites of a metabolic network. To quantitatively characterize the overall connection structure of the GSC we introduced the term 'overall closeness centralization index (OCCI)'. OCCI correlates well with the average path length of the GSC and is a useful parameter for a system-level comparison of metabolic networks of different organisms. SUPPLEMENTARY INFORMATION: http://genome.gbf.de/bioinformatics/     

Reconstruction of metabolic networks from genome data and analysis of their global structure for various organisms

MOTIVATION: Information from fully sequenced genomes makes it possible to reconstruct strain-specific global metabolic network for structural and functional studies. These networks are often very large and complex. To properly understand and analyze the global properties of metabolic networks, methods for rationally representing and quantitatively analyzing their structure are needed. RESULTS: In this work, the metabolic networks of 80 fully sequenced organisms are in silico reconstructed from genome data and an extensively revised bioreaction database. The networks are represented as directed graphs and analyzed by using the 'breadth first searching algorithm to identify the shortest pathway (path length) between any pair of the metabolites. The average path length of the networks are then calculated and compared for all the organisms. Different from previous studies the connections through current metabolites and cofactors are deleted to make the path length analysis physiologically more meaningful. The distribution of the connection degree of these networks is shown to follow the power law, indicating that the overall structure of all the metabolic networks has the characteristics of a small world network. However, clear differences exist in the network structure of the three domains of organisms. Eukaryotes and archaea have a longer average path length than bacteria. AVAILABILITY: The reaction database in excel format and the programs in VBA (Visual Basic for Applications) are available upon request. SUPPLEMENTARY MATERIAL: Bioinformatics Online.     

FTICR-mass spectrometry for high-resolution analysis in combinatorial chemistry

The diversity of compound collections required for finding lead structures in pharmaceutical research can be provided by means of combinatorial organic chemistry. The resultant enormous number of single compounds but also of compound mixtures represents a challenge for the analyst. With the introduction of Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS or FT-MS), a new and, as yet, not widespread mass spectrometric technique (a means of analysis of such compound libraries with a very high mass resolution) high mass accuracy and high sensitivity has become available. Moreover, in combination with electrospray ionization (ESI), not only high-throughput measurements via flow-injection analysis (FIA) but also coupling with separation techniques such as high-performance liquid chromatography (HPLC) or capillary electrophoresis (CE) is possible. Structural verification by way of decomposing ions (MS(n); n > or = 2) using a variety of different dissociation techniques can be performed by FTICR-MS. This is the first review specifically covering applications of FTICR-MS in the field of combinatorial chemistry.     

Mass spectrometry for metabolic flux analysis

Mass spectrometry in combination with tracer experiments based on 13C substrates can serve as a powerful tool for the modeling and analysis of intracellular fluxes and the investigation of biochemical networks. The theoretical background for the application of mass spectrometry to metabolic flux analysis is discussed. Mass spectrometry methods are especially useful to determine mass distribution of metabolites. Additional information gained from fragmentation of metabolites, e.g., by electron impact ionization, allows further localization of labeling positions, up to complete resolution of isotopomer pools. To effectively handle mass distributions in simulation experiments, a matrix based general methodology is formulated. The natural isotope distribution of carbon, oxygen, hydrogen and nitrogen in the target metabolites is considered by introduction of correction matrices. It is shown by simulation results for the central carbon metabolism that neglecting natural isotope distributions causes significant errors in intracellular flux distributions. By varying relative fluxes into pentosephosphate pathway and pyruvate carboxylation reaction, marked changes in the mass distributions of metabolites result, which are illustrated for pyruvate, oxaloacetate, and alpha-ketoglutarate. In addition mass distributions of metabolites are significantly influenced over a broad range by the degree of reversibility of transaldolase and transketolase reactions in the pentosephosphate pathway. The mass distribution of metabolites is very sensitive towards intracellular flux patterns and can be measured with high accuracy by routine mass spectrometry methods. Copyright 1999 John Wiley & Sons, Inc.     

Profiling of Arabidopsis secondary metabolites by capillary liquid chromatography coupled to electrospray ionization quadrupole time-of-flight mass spectrometry

Large-scale metabolic profiling is expected to develop into an integral part of functional genomics and systems biology. The metabolome of a cell or an organism is chemically highly complex. Therefore, comprehensive biochemical phenotyping requires a multitude of analytical techniques. Here, we describe a profiling approach that combines separation by capillary liquid chromatography with the high resolution, high sensitivity, and high mass accuracy of quadrupole time-of-flight mass spectrometry. About 2000 different mass signals can be detected in extracts of Arabidopsis roots and leaves. Many of these originate from Arabidopsis secondary metabolites. Detection based on retention times and exact masses is robust and reproducible. The dynamic range is sufficient for the quantification of metabolites. Assessment of the reproducibility of the analysis showed that biological variability exceeds technical variability. Tools were optimized or established for the automatic data deconvolution and data processing. Subtle differences between samples can be detected as tested with the chalcone synthase deficient tt4 mutant. The accuracy of time-of-flight mass analysis allows to calculate elemental compositions and to tentatively identify metabolites. In-source fragmentation and tandem mass spectrometry can be used to gain structural information. This approach has the potential to significantly contribute to establishing the metabolome of Arabidopsis and other model systems. The principles of separation and mass analysis of this technique, together with its sensitivity and resolving power, greatly expand the range of metabolic profiling.     

Does deamidation cause protein unfolding? A top‐down tandem mass spectrometry study

Deamidation is a nonenzymatic post-translational modification of asparagine to aspartic acid or glutamine to glutamic acid, converting an uncharged amino acid to a negatively charged residue. It is plausible that deamidation of asparagine and glutamine residues would result in disruption of a proteins'' hydrogen bonding network and thus lead to protein unfolding. To test this hypothesis Calmodulin and B2M were deamidated and analyzed using tandem mass spectrometry on a Fourier transform ion cyclotron resonance mass spectrometer (FTICR-MS). The gas phase hydrogen bonding networks of deamidated and nondeamidated protein isoforms were probed by varying the infra-red multi-photon dissociation laser power in a linear fashion and plotting the resulting electron capture dissociation fragment intensities as a melting curve at each amino acid residue. Analysis of the unfolding maps highlighted increased fragmentation at lower laser powers localized around heavily deamidated regions of the proteins. In addition fragment intensities were decreased across the rest of the proteins which we propose is because of the formation of salt-bridges strengthening the intramolecular interactions of the central regions. These results were supported by a computational flexibility analysis of the mutant and unmodified proteins, which would suggest that deamidation can affect the global structure of a protein via modification of the hydrogen bonding network near the deamidation site and that top down FTICR-MS is an appropriate technique for studying protein folding.     

Thermodynamics-based Metabolite Sensitivity Analysis in metabolic networks

The increasing availability of large metabolomics datasets enhances the need for computational methodologies that can organize the data in a way that can lead to the inference of meaningful relationships. Knowledge of the metabolic state of a cell and how it responds to various stimuli and extracellular conditions can offer significant insight in the regulatory functions and how to manipulate them. Constraint based methods, such as Flux Balance Analysis (FBA) and Thermodynamics-based flux analysis (TFA), are commonly used to estimate the flow of metabolites through genome-wide metabolic networks, making it possible to identify the ranges of flux values that are consistent with the studied physiological and thermodynamic conditions. However, unless key intracellular fluxes and metabolite concentrations are known, constraint-based models lead to underdetermined problem formulations. This lack of information propagates as uncertainty in the estimation of fluxes and basic reaction properties such as the determination of reaction directionalities. Therefore, knowledge of which metabolites, if measured, would contribute the most to reducing this uncertainty can significantly improve our ability to define the internal state of the cell. In the present work we combine constraint based modeling, Design of Experiments (DoE) and Global Sensitivity Analysis (GSA) into the Thermodynamics-based Metabolite Sensitivity Analysis (TMSA) method. TMSA ranks metabolites comprising a metabolic network based on their ability to constrain the gamut of possible solutions to a limited, thermodynamically consistent set of internal states. TMSA is modular and can be applied to a single reaction, a metabolic pathway or an entire metabolic network. This is, to our knowledge, the first attempt to use metabolic modeling in order to provide a significance ranking of metabolites to guide experimental measurements.     

Use of the Local False Discovery Rate for Identification of Metabolic Biomarkers in Rat Urine Following Genkwa Flos-Induced Hepatotoxicity

Metabolomics is concerned with characterizing the large number of metabolites present in a biological system using nuclear magnetic resonance (NMR) and HPLC/MS (high-performance liquid chromatography with mass spectrometry). Multivariate analysis is one of the most important tools for metabolic biomarker identification in metabolomic studies. However, analyzing the large-scale data sets acquired during metabolic fingerprinting is a major challenge. As a posterior probability that the features of interest are not affected, the local false discovery rate (LFDR) is a good interpretable measure. However, it is rarely used to when interrogating metabolic data to identify biomarkers. In this study, we employed the LFDR method to analyze HPLC/MS data acquired from a metabolomic study of metabolic changes in rat urine during hepatotoxicity induced by Genkwa flos (GF) treatment. The LFDR approach was successfully used to identify important rat urine metabolites altered by GF-stimulated hepatotoxicity. Compared with principle component analysis (PCA), LFDR is an interpretable measure and discovers more important metabolites in an HPLC/MS-based metabolomic study.     

Chemical basis of metabolic network organization

Although the metabolic networks of the three domains of life consist of different constituents and metabolic pathways, they exhibit the same scale-free organization. This phenomenon has been hypothetically explained by preferential attachment principle that the new-recruited metabolites attach preferentially to those that are already well connected. However, since metabolites are usually small molecules and metabolic processes are basically chemical reactions, we speculate that the metabolic network organization may have a chemical basis. In this paper, chemoinformatic analyses on metabolic networks of Kyoto Encyclopedia of Genes and Genomes (KEGG), Escherichia coli and Saccharomyces cerevisiae were performed. It was found that there exist qualitative and quantitative correlations between network topology and chemical properties of metabolites. The metabolites with larger degrees of connectivity (hubs) are of relatively stronger polarity. This suggests that metabolic networks are chemically organized to a certain extent, which was further elucidated in terms of high concentrations required by metabolic hubs to drive a variety of reactions. This finding not only provides a chemical explanation to the preferential attachment principle for metabolic network expansion, but also has important implications for metabolic network design and metabolite concentration prediction.     

Predicting Network Activity from High Throughput Metabolomics

The functional interpretation of high throughput metabolomics by mass spectrometry is hindered by the identification of metabolites, a tedious and challenging task. We present a set of computational algorithms which, by leveraging the collective power of metabolic pathways and networks, predict functional activity directly from spectral feature tables without a priori identification of metabolites. The algorithms were experimentally validated on the activation of innate immune cells.     

Metabolomics and microarrays for improved understanding of phenotypic characteristics controlled by both genomics and environmental constraints

Advances in our understanding of functional genomics are best addressed by integrative studies that include measurements of mRNA, proteins, and low molecular weight metabolites over time and varied conditions. Bioinformatics can then be used to relate this data to the genome. Current technology allows for comprehensive and rapid mRNA expression profiling and mass spectrophotometric measurement of low molecular weight intermediates and metabolic products. In prokaryotic organisms, this combination provides a potentially powerful tool for identifying gene function and regulatory networks even in the absence of a combined proteomic approach.     

Multi-Organ Contribution to the Metabolic Plasma Profile Using Hierarchical Modelling

Hierarchical modelling was applied in order to identify the organs that contribute to the levels of metabolites in plasma. Plasma and organ samples from gut, kidney, liver, muscle and pancreas were obtained from mice. The samples were analysed using gas chromatography time-of-flight mass spectrometry (GC TOF-MS) at the Swedish Metabolomics centre, Umeå University, Sweden. The multivariate analysis was performed by means of principal component analysis (PCA) and orthogonal projections to latent structures (OPLS). The main goal of this study was to investigate how each organ contributes to the metabolic plasma profile. This was performed using hierarchical modelling. Each organ was found to have a unique metabolic profile. The hierarchical modelling showed that the gut, kidney and liver demonstrated the greatest contribution to the metabolic pattern of plasma. For example, we found that metabolites were absorbed in the gut and transported to the plasma. The kidneys excrete branched chain amino acids (BCAAs) and fatty acids are transported in the plasma to the muscles and liver. Lactic acid was also found to be transported from the pancreas to plasma. The results indicated that hierarchical modelling can be utilized to identify the organ contribution of unknown metabolites to the metabolic profile of plasma.