Proposed minimum reporting standards for chemical analysis

There is a general consensus that supports the need for standardized reporting of metadata or information describing large-scale metabolomics and other functional genomics data sets. Reporting of standard metadata provides a biological and empirical context for the data, facilitates experimental replication, and enables the re-interrogation and comparison of data by others. Accordingly, the Metabolomics Standards Initiative is building a general consensus concerning the minimum reporting standards for metabolomics experiments of which the Chemical Analysis Working Group (CAWG) is a member of this community effort. This article proposes the minimum reporting standards related to the chemical analysis aspects of metabolomics experiments including: sample preparation, experimental analysis, quality control, metabolite identification, and data pre-processing. These minimum standards currently focus mostly upon mass spectrometry and nuclear magnetic resonance spectroscopy due to the popularity of these techniques in metabolomics. However, additional input concerning other techniques is welcomed and can be provided via the CAWG on-line discussion forum at or . Further, community input related to this document can also be provided via this electronic forum. The contents of this paper do not necessarily reflect any position of the Government or the opinion of the Food and Drug Administration Sponsor: Metabolomics Society http://www.metabolomicssociety.org/ Reference: http://msi-workgroups.sourceforge.net/bio-metadata/reporting/pbc/ http://msi-workgroups.sourceforge.net/chemical-analysis/ Version: Revision: 5.1 Date: 09 January, 2007     

Metabolic isotopomer labeling systems. Part II: structural flux identifiability analysis

Metabolic flux analysis using carbon labeling experiments (CLEs) is an important tool in metabolic engineering where the intracellular fluxes have to be computed from the measured extracellular fluxes and the partially measured distribution of 13C labeling within the intracellular metabolite pools. The relation between unknown fluxes and measurements is described by an isotopomer labeling system (ILS) (see Part I [Math. Biosci. 169 (2001) 173]). Part II deals with the structural flux identifiability of measured ILSs in the steady state. The central question is whether the measured data contains sufficient information to determine the unknown intracellular fluxes. This question has to be decided a priori, i.e. before the CLE is carried out. In structural identifiability analysis the measurements are assumed to be noise-free. A general theory of structural flux identifiability for measured ILSs is presented and several algorithms are developed to solve the identifiability problem. In the particular case of maximal measurement information, a symbolical algorithm is presented that decides the identifiability question by means of linear methods. Several upper bounds of the number of identifiable fluxes are derived, and the influence of the chosen inputs is evaluated. By introducing integer arithmetic this algorithm can even be applied to large networks. For the general case of arbitrary measurement information, identifiability is decided by a local criterion. A new algorithm based on integer arithmetic enables an a priori local identifiability analysis to be performed for networks of arbitrary size. All algorithms have been implemented and flux identifiability is investigated for the network of the central metabolic pathways of a microorganism. Moreover, several small examples are worked out to illustrate the influence of input metabolite labeling and the paradox of information loss due to network simplification.     

Bio-based succinate from sucrose: High-resolution 13C metabolic flux analysis and metabolic engineering of the rumen bacterium Basfia succiniciproducens

Succinic acid is a platform chemical of recognized industrial value and accordingly faces a continuous challenge to enable manufacturing from most attractive raw materials. It is mainly produced from glucose, using microbial fermentation. Here, we explore and optimize succinate production from sucrose, a globally applied substrate in biotechnology, using the rumen bacterium Basfia succiniciproducens DD1. As basis of the strain optimization, the yet unknown sucrose metabolism of the microbe was studied, using ¹³C metabolic flux analyses. When grown in batch culture on sucrose, the bacterium exhibited a high succinate yield of 1 mol mol⁻¹ and a by-product spectrum, which did not match the expected PTS-mediated sucrose catabolism. This led to the discovery of a fructokinase, involved in sucrose catabolism. The flux approach unraveled that the fructokinase and the fructose PTS both contribute to phosphorylation of the fructose part of sucrose. The contribution of the fructokinase reduces the undesired loss of the succinate precursor PEP into pyruvate and into pyruvate-derived by-products and enables increased succinate production, exclusively via the reductive TCA cycle branch. These findings were used to design superior producers. Mutants, which (i) overexpress the beneficial fructokinase, (II) lack the competing fructose PTS, and (iii) combine both traits, produce significantly more succinate. In a fed-batch process, B. succiniciproducens ΔfruA achieved a titer of 71 g L⁻¹ succinate and a yield of 2.5 mol mol⁻¹ from sucrose.     

Estradiol stimulates the biosynthetic pathways of breast cancer cells: detection by metabolic flux analysis

Selective estrogen receptor (ER) modulators are highly successful breast cancer therapies, but they are not effective in patients with ER negative and selective estrogen receptor modulator (SERM)-resistant tumors. Understanding the mechanisms of estrogen-stimulated proliferation may provide a route to design estrogen-independent therapies that would be effective in these patients. In this study, metabolic flux analysis was used to determine the intracellular fluxes that are significantly affected by estradiol stimulation in MCF-7 breast cancer cells. Intracellular fluxes were calculated from nuclear magnetic resonance (NMR)-generated isotope enrichment data and extracellular metabolite fluxes, using a specific flux analysis algorithm. The metabolic pathway model used by the algorithm includes glycolysis, the tricarboxylic acid cycle (TCA cycle), the pentose phosphate pathway, glutamine catabolism, pyruvate carboxylase, and malic enzyme. The pathway model also incorporates mitochondrial compartmentalization and reversible trans-mitochondrial membrane reactions to more accurately describe the role of mitochondria in cancer cell proliferation. Flux results indicate that estradiol significantly increases carbon flow through the pentose phosphate pathway and increases glutamine consumption. In addition, intra-mitochondrial malic enzyme was found to be inactive and the malate-aspartate shuttle (MAS) was only minimally active. The inactivity of these enzymes indicates that glutamine is not oxidized within mitochondria, but is consumed primarily to provide biosynthetic precursors. The excretion of glutamine carbons from the mitochondria has the secondary effect of limiting nicotinamide adenine dinucleotide (NADH) recycle, resulting in NADH buildup in the cytosol and the excretion of lactate. The observed dependence of breast cancer cells on pentose phosphate pathway activity and glutamine consumption for estradiol-stimulated biosynthesis suggests that these pathways may be targets for estrogen-independent breast cancer therapies.     

Storage and oxidation of long-chain fatty acids in the C57/BL6 mouse heart as measured by NMR spectroscopy

Triglyceride turnover in the isolated C57/BL6 mouse heart was measured by dynamic 13C edit-(1)H observe NMR and the rate of fatty acid oxidation was determined by 13C NMR isotopomer analysis. In the presence of a physiological mixture of substrates, energy was produced in the citric acid cycle by oxidation of long-chain fatty acids (18%), ketones (34%), lactate (24%), pyruvate (7%), and other sources (17%). Exogenous fatty acids appeared in the triglyceride pool at 0.24 micromol/g dry wt/min, similar to the rate of oxidation of long-chain fatty acids, 0.16 micromol/g dry wt/min. Isoproterenol decreased the rate of de novo triglyceride synthesis and increased the rate of fatty acid oxidation.     

Metabolic Analysis of Wild-type Escherichia coli and a Pyruvate Dehydrogenase Complex (PDHC)-deficient Derivative Reveals the Role of PDHC in the Fermentative Metabolism of Glucose

Pyruvate is located at a metabolic junction of assimilatory and dissimilatory pathways and represents a switch point between respiratory and fermentative metabolism. In Escherichia coli, the pyruvate dehydrogenase complex (PDHC) and pyruvate formate-lyase are considered the primary routes of pyruvate conversion to acetyl-CoA for aerobic respiration and anaerobic fermentation, respectively. During glucose fermentation, the in vivo activity of PDHC has been reported as either very low or undetectable, and the role of this enzyme remains unknown. In this study, a comprehensive characterization of wild-type E. coli MG1655 and a PDHC-deficient derivative (Pdh) led to the identification of the role of PDHC in the anaerobic fermentation of glucose. The metabolism of these strains was investigated by using a mixture of ¹³C-labeled and -unlabeled glucose followed by the analysis of the labeling pattern in protein-bound amino acids via two-dimensional ¹³C,¹H NMR spectroscopy. Metabolite balancing, biosynthetic ¹³C labeling of proteinogenic amino acids, and isotopomer balancing all indicated a large increase in the flux of the oxidative branch of the pentose phosphate pathway (ox-PPP) in response to the PDHC deficiency. Because both ox-PPP and PDHC generate CO₂ and the calculated CO₂ evolution rate was significantly reduced in Pdh, it was hypothesized that the role of PDHC is to provide CO₂ for cell growth. The similarly negative impact of either PDHC or ox-PPP deficiencies, and an even more pronounced impairment of cell growth in a strain lacking both ox-PPP and PDHC, provided further support for this hypothesis. The three strains exhibited similar phenotypes in the presence of an external source of CO₂, thus confirming the role of PDHC. Activation of formate hydrogen-lyase (which converts formate to CO₂ and H₂) rendered the PDHC deficiency silent, but its negative impact reappeared in a strain lacking both PDHC and formate hydrogen-lyase. A stoichiometric analysis of CO₂ generation via PDHC and ox-PPP revealed that the PDHC route is more carbon- and energy-efficient, in agreement with its beneficial role in cell growth.     

Kinetic isotope effects significantly influence intracellular metabolite 13C labeling patterns and flux determination

Rigorous mathematical modeling of carbon-labeling experiments allows estimation of fluxes through the pathways of central carbon metabolism, yielding powerful information for basic scientific studies as well as for a wide range of applications. However, the mathematical models that have been developed for flux determination from ¹³C labeling data have commonly neglected the influence of kinetic isotope effects on the distribution of ¹³C label in intracellular metabolites, as these effects have often been assumed to be inconsequential. We have used measurements of the ¹³C isotope effects on the pyruvate dehydrogenase enzyme from the literature to model isotopic fractionation at the pyruvate node and quantify the modeling errors expected to result from the assumption that isotope effects are negligible. We show that under some conditions kinetic isotope effects have a significant impact on the ¹³C labeling patterns of intracellular metabolites, and the errors associated with neglecting isotope effects in ¹³C-metabolic flux analysis models can be comparable in size to measurement errors associated with GC–MS. Thus, kinetic isotope effects must be considered in any rigorous assessment of errors in ¹³C labeling data, goodness-of-fit between model and data, confidence intervals of estimated metabolic fluxes, and statistical significance of differences between estimated metabolic flux distributions.     

Quantitative metabolome analysis using liquid chromatography-high-resolution mass spectrometry

In this report, we introduce a liquid chromatography single-mass spectrometry method for metabolome quantification, using the LTQ Orbitrap high-resolution mass spectrometer. Analytes were separated with hydrophilic interaction liquid chromatography. At a working resolution of 30,000 (at m/z 400), the limit of detection varied from 50 fmol to 5 pmol for 25 metabolites tested. In terms of metabolite concentration, the linearity was about 2 to 3 orders of magnitude for most compounds (R² > 0.99). To determine the accuracy of the system in complex sample matrices, the isotope dilution method was evaluated from mixtures of pure compounds and uniformly ¹³C-labeled cell extracts. With the application of this method, quantification was possible within single runs even when the pool sizes of individual metabolites varied from 0.13 to 55.6 μM. As a case study, intracellular concentrations of central metabolites were determined for Methylobacterium extorquens AM1 during growth on two different carbon sources, methanol and succinate. Reproducible results from technical and biological repetitions were obtained that revealed significant variations of intracellular metabolite pool sizes, depending on the carbon source. The LTQ Obitrap offers new perspectives and strategies for metabolome quantification.     

Metabolism of levulinate in perfused rat livers and live rats: conversion to the drug of abuse 4-hydroxypentanoate

Calcium levulinate (4-ketopentanoate) is used as an oral and parenteral source of calcium. We hypothesized that levulinate is converted in the liver to 4-hydroxypentanoate, a new drug of abuse, and that this conversion is accelerated by ethanol oxidation. We confirmed these hypotheses in live rats, perfused rat livers, and liver subcellular preparations. Levulinate is reduced to (R)-4-hydroxypentanoate by a cytosolic and a mitochondrial dehydrogenase, which are NADPH- and NADH-dependent, respectively. A mitochondrial dehydrogenase or racemase system also forms (S)-4-hydroxypentanoate. In livers perfused with [(13)C(5)]levulinate, there was substantial CoA trapping in levulinyl-CoA, 4-hydroxypentanoyl-CoA, and 4-phosphopentanoyl-CoA. This CoA trapping was increased by ethanol, with a 6-fold increase in the concentration of 4-phosphopentanoyl-CoA. Levulinate is catabolized by 3 parallel pathways to propionyl-CoA, acetyl-CoA, and lactate. Most intermediates of the 3 pathways were identified by mass isotopomer analysis and metabolomics. The production of 4-hydroxypentanoate from levulinate and its stimulation by ethanol is a potential public health concern.