Metabolomic and mass isotopomer analysis of liver gluconeogenesis and citric acid cycle. I. Interrelation between gluconeogenesis and cataplerosis; formation of methoxamates from aminooxyacetate and ketoacids

We conducted a study coupling metabolomics and mass isotopomer analysis of liver gluconeogenesis and citric acid cycle. Rat livers were perfused with lactate or pyruvate +/- aminooxyacetate or mercaptopicolinate in the presence of 40% enriched NaH(13)CO(3). Other livers were perfused with dimethyl [1,4-(13)C(2)]succinate +/- mercaptopicolinate. In this first of two companion articles, we show that a substantial fraction of gluconeogenic carbon leaves the liver as citric acid cycle intermediates, mostly alpha-ketoglutarate. The efflux of gluconeogenic carbon ranges from 10 to 200% of the rate of liver gluconeogenesis. This cataplerotic efflux of gluconeogenic carbon may contribute to renal gluconeogenesis in vivo. Multiple crossover analyses of concentrations of gluconeogenic intermediates and redox measurements expand previous reports on the regulation of gluconeogenesis and the effects of inhibitors. We also demonstrate the formation of adducts from the condensation, in the liver, of (i) aminooxyacetate with pyruvate, alpha-ketoglutarate, and oxaloacetate and (ii) mercaptopicolinate and pyruvate. These adducts may exert metabolic effects unrelated to their effect on gluconeogenesis.     

Metabolic flux analysis for a ppc mutant Escherichia coli based on 13C-labelling experiments together with enzyme activity assays and intracellular metabolite measurements

The physiology and central metabolism of a ppc mutant Escherichia coli were investigated based on the metabolic flux distribution obtained by (13)C-labelling experiments using gas chromatography-mass spectrometry (GC-MS) and 2-dimensional nuclear magnetic resonance (2D NMR) strategies together with enzyme activity assays and intracellular metabolite concentration measurements. Compared to the wild type, its ppc mutant excreted little acetate and produced less carbon dioxide at the expense of a slower growth rate and a lower glucose uptake rate. Consequently, an improvement of the biomass yield on glucose was observed in the ppc mutant. Enzyme activity measurements revealed that isocitrate lyase activity increased by more than 3-fold in the ppc mutant. Some TCA cycle enzymes such as citrate synthase, aconitase and malate dehydrogenase were also upregulated, but enzymes of glycolysis and the pentose phosphate pathway were downregulated. The intracellular intermediates in the glycolysis and the pentose phosphate pathway, therefore, accumulated, while acetyl coenzyme A and oxaloacetate concentrations decreased in the ppc mutant. The intracellular metabolic flux analysis uncovered that deletion of ppc resulted in the appearance of the glyoxylate shunt, with 18.9% of the carbon flux being channeled via the glyoxylate shunt. However, the flux of the pentose phosphate pathway significantly decreased in the ppc mutant.     

Hierarchy in Pentose Sugar Metabolism in Clostridium acetobutylicum

Bacterial metabolism of polysaccharides from plant detritus into acids and solvents is an essential component of the terrestrial carbon cycle. Understanding the underlying metabolic pathways can also contribute to improved production of biofuels. Using a metabolomics approach involving liquid chromatography-mass spectrometry, we investigated the metabolism of mixtures of the cellulosic hexose sugar (glucose) and hemicellulosic pentose sugars (xylose and arabinose) in the anaerobic soil bacterium Clostridium acetobutylicum. Simultaneous feeding of stable isotope-labeled glucose and unlabeled xylose or arabinose revealed that, as expected, glucose was preferentially used as the carbon source. Assimilated pentose sugars accumulated in pentose phosphate pathway (PPP) intermediates with minimal flux into glycolysis. Simultaneous feeding of xylose and arabinose revealed an unexpected hierarchy among the pentose sugars, with arabinose utilized preferentially over xylose. The phosphoketolase pathway (PKP) provides an alternative route of pentose catabolism in C. acetobutylicum that directly converts xylulose-5-phosphate into acetyl-phosphate and glyceraldehyde-3-phosphate, bypassing most of the PPP. When feeding the mixture of pentose sugars, the labeling patterns of lower glycolytic intermediates indicated more flux through the PKP than through the PPP and upper glycolysis, and this was confirmed by quantitative flux modeling. Consistent with direct acetyl-phosphate production from the PKP, growth on the pentose mixture resulted in enhanced acetate excretion. Taken collectively, these findings reveal two hierarchies in clostridial pentose metabolism: xylose is subordinate to arabinose, and the PPP is used less than the PKP.     

Sodium-Dependent Growth of Azotobacter chroococcum

Selected ion-monitoring gas chromatography-mass spectrometry was used for detection of beta-hydroxy fatty acids as an independent assay for the presence or absence of endotoxin in materials claimed to induce nonspecific activation of Limulus amoebocyte lysate. To this end, suspensions of gram-negative and -positive bacteria, one fungal species, cerebrospinal fluid, and hollow-fiber hemodialyzer rinses were assayed for endotoxin by gas chromatography-mass spectrometry and the Limulus amoebocyte lysate assay. Good qualitative agreement was shown for both methods when suspensions of test organisms were assayed. Two false-negative results were obtained by gas chromatography-mass spectrometry assays of cerebrospinal fluid and were shown to be a result of insufficient endotoxin in the cerebrospinal fluid specimens for detection by gas chromatography-mass spectrometry. Hemodialyzer rinses were Limulus assay positive; however, no beta-hydroxy fatty acids were detected by gas chromatography-mass spectrometry. These data were compared with data obtained from USP rabbit pyrogen tests of the rinse materials (nonpyrogenic) and chemical characterization of the Limulus assay-reactive rinses, which showed the rinses to be cellulosic in nature. It is suggested that beta-hydroxy fatty acids, as assayed by selected ion-monitoring gas chromatography-mass spectrometry, be used as chemical marker molecules for the presence or absence of endotoxin in materials reported to cause nonspecific activation of Limulus amoebocyte lysate.     

Sodium-Dependent Growth of Azotobacter chroococcum

Selected ion-monitoring gas chromatography-mass spectrometry was used for detection of beta-hydroxy fatty acids as an independent assay for the presence or absence of endotoxin in materials claimed to induce nonspecific activation of Limulus amoebocyte lysate. To this end, suspensions of gram-negative and -positive bacteria, one fungal species, cerebrospinal fluid, and hollow-fiber hemodialyzer rinses were assayed for endotoxin by gas chromatography-mass spectrometry and the Limulus amoebocyte lysate assay. Good qualitative agreement was shown for both methods when suspensions of test organisms were assayed. Two false-negative results were obtained by gas chromatography-mass spectrometry assays of cerebrospinal fluid and were shown to be a result of insufficient endotoxin in the cerebrospinal fluid specimens for detection by gas chromatography-mass spectrometry. Hemodialyzer rinses were Limulus assay positive; however, no beta-hydroxy fatty acids were detected by gas chromatography-mass spectrometry. These data were compared with data obtained from USP rabbit pyrogen tests of the rinse materials (nonpyrogenic) and chemical characterization of the Limulus assay-reactive rinses, which showed the rinses to be cellulosic in nature. It is suggested that beta-hydroxy fatty acids, as assayed by selected ion-monitoring gas chromatography-mass spectrometry, be used as chemical marker molecules for the presence or absence of endotoxin in materials reported to cause nonspecific activation of Limulus amoebocyte lysate.     

Comparison of the limulus amebocyte lysate test and gas chromatography-mass spectrometry for measuring lipopolysaccharides (endotoxins) in airborne dust from poultry-processing industries

The lipopolysaccharide (endotoxin) content in airborne dust samples from three different poultry slaughterhouses was determined with both the chromogenic Limulus amebocyte lysate assay and gas chromatography-mass spectrometry analysis of lipopolysaccharide-derived 3-hydroxy fatty acids. Gram-negative cell walls were also measured by using two-dimensional gas chromatography/electron-capture analysis of diaminopimelic acid originating from the peptidoglycan. The correlation between the results of the Limulus assay and those of gas chromatography-mass spectrometry for determination of the lipopolysaccharide content in the dust samples was poor, whereas a good correlation was obtained between lipopolysaccharide and diaminopimelic acid concentrations with the gas chromatographic methods. The results suggest that it is predominantly cell-wall-dissociated lipopolysaccharides that are measured with the Limulus assay, whereas the gas chromatographic methods allow determination of total concentrations of lipopolysaccharide, including Limulus-inactive lipopolysaccharide, gram-negative cells, and cellular debris.     

Comparison of the limulus amebocyte lysate test and gas chromatography-mass spectrometry for measuring lipopolysaccharides (endotoxins) in airborne dust from poultry-processing industries.

The lipopolysaccharide (endotoxin) content in airborne dust samples from three different poultry slaughterhouses was determined with both the chromogenic Limulus amebocyte lysate assay and gas chromatography-mass spectrometry analysis of lipopolysaccharide-derived 3-hydroxy fatty acids. Gram-negative cell walls were also measured by using two-dimensional gas chromatography/electron-capture analysis of diaminopimelic acid originating from the peptidoglycan. The correlation between the results of the Limulus assay and those of gas chromatography-mass spectrometry for determination of the lipopolysaccharide content in the dust samples was poor, whereas a good correlation was obtained between lipopolysaccharide and diaminopimelic acid concentrations with the gas chromatographic methods. The results suggest that it is predominantly cell-wall-dissociated lipopolysaccharides that are measured with the Limulus assay, whereas the gas chromatographic methods allow determination of total concentrations of lipopolysaccharide, including Limulus-inactive lipopolysaccharide, gram-negative cells, and cellular debris.     

Understanding in vivo benzenoid metabolism in petunia petal tissue

In vivo stable isotope labeling and computer-assisted metabolic flux analysis were used to investigate the metabolic pathways in petunia (Petunia hybrida) cv Mitchell leading from Phe to benzenoid compounds, a process that requires the shortening of the side chain by a C(2) unit. Deuterium-labeled Phe ((2)H(5)-Phe) was supplied to excised petunia petals. The intracellular pools of benzenoid/phenylpropanoid-related compounds (intermediates and end products) as well as volatile end products within the floral bouquet were analyzed for pool sizes and labeling kinetics by gas chromatography-mass spectrometry and liquid chromatography-mass spectrometry. Modeling of the benzenoid network revealed that both the CoA-dependent, beta-oxidative and CoA-independent, non-beta-oxidative pathways contribute to the formation of benzenoid compounds in petunia flowers. The flux through the CoA-independent, non-beta-oxidative pathway with benzaldehyde as a key intermediate was estimated to be about 2 times higher than the flux through the CoA-dependent, beta-oxidative pathway. Modeling of (2)H(5)-Phe labeling data predicted that in addition to benzaldehyde, benzylbenzoate is an intermediate between l-Phe and benzoic acid. Benzylbenzoate is the result of benzoylation of benzyl alcohol, for which activity was detected in petunia petals. A cDNA encoding a benzoyl-CoA:benzyl alcohol/phenylethanol benzoyltransferase was isolated from petunia cv Mitchell using a functional genomic approach. Biochemical characterization of a purified recombinant benzoyl-CoA:benzyl alcohol/phenylethanol benzoyltransferase protein showed that it can produce benzylbenzoate and phenylethyl benzoate, both present in petunia corollas, with similar catalytic efficiencies.     

A targeted metabolomics assay for cardiac metabolism and demonstration using a mouse model of dilated cardiomyopathy

Metabolomics can be performed either as an ‘open profiling’ tool where the aim is to measure, usually in a semi-quantitative manner, as many metabolites as possible or perform ‘closed’ or ‘targeted’ analyses where instead a pre-defined set of metabolites are measured. Targeted methods can be designed to be more sensitive and quantitative and so are particularly appropriate to systems biology for quantitative models of systems or when metabolomics is performed in a hypothesis driven manner to test whether a particular pathway is perturbed. We describe a targeted metabolomics assay that quantifies a broad range of over 130 metabolites relevant to cardiac metabolism including the pathways of the citric acid cycle, fatty acid oxidation, glycolysis, the pentose phosphate pathway, amino acid metabolism, the urea cycle, nucleotides and reactive oxygen species using tandem mass spectrometry to produce quantitative, sensitive and robust data. This assay is illustrated by profiling cardiac metabolism in a lamin A/C (Lmna) mouse model of dilated cardiomyopathy (DCM). The model of DCM was characterised by increases in concentrations of proline and methyl-histidine suggestive of increased myofibrillar and collagen degradation, as well as decreases in a number of citric acid cycle intermediates and carnitine derivatives indicating reduced energy metabolism in the dilated heart. These assays could be used for any other cardiac or cardiovascular disease in that they cover central core metabolism and key pathways involved in cardiac metabolism, and may provide a general start for many mammalian systems.     

Catabolism of 4-Hydroxyacids and 4-Hydroxynonenal via 4-Hydroxy-4-phosphoacyl-CoAs

4-Hydroxyacids are products of ubiquitously occurring lipid peroxidation (C₉, C₆) or drugs of abuse (C₄, C₅). We investigated the catabolism of these compounds using a combination of metabolomics and mass isotopomer analysis. Livers were perfused with various concentrations of unlabeled and labeled saturated 4-hydroxyacids (C₄ to C₁₁) or 4-hydroxynonenal. All the compounds tested form a new class of acyl-CoA esters, 4-hydroxy-4-phosphoacyl-CoAs, characterized by liquid chromatography-tandem mass spectrometry, accurate mass spectrometry, and ³¹P-NMR. All 4-hydroxyacids with five or more carbons are metabolized by two new pathways. The first and major pathway, which involves 4-hydroxy-4-phosphoacyl-CoAs, leads in six steps to the isomerization of 4-hydroxyacyl-CoA to 3-hydroxyacyl-CoAs. The latter are intermediates of physiological β-oxidation. The second and minor pathway involves a sequence of β-oxidation, α-oxidation, and β-oxidation steps. In mice deficient in succinic semialdehyde dehydrogenase, high plasma concentrations of 4-hydroxybutyrate result in high concentrations of 4-hydroxy-4-phospho-butyryl-CoA in brain and liver. The high concentration of 4-hydroxy-4-phospho-butyryl-CoA may be related to the cerebral dysfunction of subjects ingesting 4-hydroxybutyrate and to the mental retardation of patients with 4-hydroxybutyric aciduria. Our data illustrate the potential of the combination of metabolomics and mass isotopomer analysis for pathway discovery.     

GC/MS-based metabolomic studies reveal key roles of glycine in regulating silk synthesis in silkworm,Bombyx mori

Metabolic profiling of silkworm, especially the factors that affect silk synthesis at the metabolic level, is little known. Herein, metabolomic method based on gas chromatography-mass spectrometry was applied to identify key metabolic changes in silk synthesis deficient silkworms. Forty-six differential metabolites were identified in Nd group with the defect of silk synthesis. Significant changes in the levels of glycine and uric acid (up-regulation), carbohydrates and free fatty acids (down-regulation) were observed. The further metabolomics of silk synthesis deficient silkworms by decreasing silk proteins synthesis using knocking out fibroin heavy chain gene or extirpating silk glands operation showed that the changes of the metabolites were almost consistent with those of the Nd group. Furthermore, the increased silk yields by supplying more glycine or its related metabolite confirmed that glycine is a key metabolite to regulate silk synthesis. These findings provide important insights into the regulation between metabolic profiling and silk synthesis.     

Stimulation by glucose of gluconeogenesis in hepatocytes isolated from starved rats.

Control properties of the gluconeogenic pathway in hepatocytes isolated from starved rats were studied in the presence of glucose. The following observations were made. (1) Glucose stimulated the rate of glucose production from 20 mM-glycerol, from a mixture of 20 mM-lactate and 2 mM-pyruvate, or from pyruvate alone; no stimulation was observed with 20 mM-alanine or 20 mM-dihydroxyacetone. Maximal stimulation was obtained between 2 and 5 mM-glucose, depending on the conditions. At concentrations above 6 mM, gluconeogenesis declined again, so that at 10 mM-glucose the glucose production rate became equal to that in its absence. (2) With glycerol, stimulation of gluconeogenesis by glucose was accompanied by oxidation of cytosolic NADH and reduction of mitochondrial NAD+ and was insensitive to the transaminase inhibitor amino-oxyacetate; this indicated that glucose accelerated the rate of transport of cytosolic reducing equivalents to the mitochondria via the glycerol 1-phosphate shuttle. (3) With lactate plus pyruvate (10:1) as substrates, stimulation of gluconeogenesis by glucose was almost additive to that obtained with glucagon. From an analysis of the effect of glucose on the curves relating gluconeogenic flux and the steady-state intracellular concentrations of gluconeogenic intermediates under various conditions, in the absence and presence of glucagon, it was concluded that addition of glucose stimulated both phosphoenolpyruvate carboxykinase and pyruvate carboxylase activity.     

Development of a metabolomic approach based on urine samples and direct infusion mass spectrometry

The analysis of urine by direct infusion mass spectrometry suffers from ion suppression due to its high salt content and inter-sample variability caused by the differences in urine volume between persons. Thus, urine metabolomics requires a careful selection of the sample preparation procedure and a normalization strategy to deal with these problems. Several approaches were tested for metabolomic analysis of urine samples by direct infusion electrospray mass spectrometry (DI–ESI–MS), including solid phase extraction, liquid–liquid extraction, and sample dilution. In addition, normalization of results based on conductivity values and statistical treatment was performed to minimize sample variability. Both urine dilution and solid phase extraction with mixed mode sorbent considerably reduced the salt content in urine, providing comprehensive metabolomic fingerprints. Moreover, statistical data normalization enabled the correction of inter-sample physiological variability, improving the quality of results obtained. Therefore, high-throughput DI–ESI–MS fingerprinting of urine samples can be achieved with simple pretreatment procedures allowing the use of this noninvasive sampling in metabolomics. Finally, the optimized approach was tested in a pilot metabolomic investigation of urine samples from transgenic mice models of Alzheimer’s disease (APP/PS1) in order to illustrate the potential of the methodology.     

Evaluation of regression models in metabolic physiology: predicting fluxes from isotopic data without knowledge of the pathway

This study explores the ability of regression models, with no knowledge of the underlying physiology, to estimate physiological parameters relevant for metabolism and endocrinology. Four regression models were compared: multiple linear regression (MLR), principal component regression (PCR), partial least-squares regression (PLS) and regression using artificial neural networks (ANN). The pathway of mammalian gluconeogenesis was analyzed using [U−¹³C]glucose as tracer. A set of data was simulated by randomly selecting physiologically appropriate metabolic fluxes for the 9 steps of this pathway as independent variables. The isotope labeling patterns of key intermediates in the pathway were then calculated for each set of fluxes, yielding 29 dependent variables. Two thousand sets were created, allowing independent training and test data. Regression models were asked to predict the nine fluxes, given only the 29 isotopomers. For large training sets (>50) the artificial neural network model was superior, capturing 95% of the variability in the gluconeogenic flux, whereas the three linear models captured only 75%. This reflects the ability of neural networks to capture the inherent non-linearities of the metabolic system. The effect of error in the variables and the addition of random variables to the data set was considered. Model sensitivities were used to find the isotopomers that most influenced the predicted flux values. These studies provide the first test of multivariate regression models for the analysis of isotopomer flux data. They provide insight for metabolomics and the future of isotopic tracers in metabolic research where the underlying physiology is complex or unknown.We acknowledge the support of NIH Grant DK58533 and the DuPont-MIT Alliance.     

Characterization of oil-producing yeast Lipomyces starkeyi on glycerol carbon source based on metabolomics and 13C-labeling

Lipomyces starkeyi is an oil-producing yeast that can produce triacylglycerol (TAG) from glycerol as a carbon source. The TAG was mainly produced after nitrogen depletion alongside reduced cell proliferation. To obtain clues for enhancing the TAG production, cell metabolism during the TAG-producing phase was characterized by metabolomics with ¹³C labeling. The turnover analysis showed that the time constants of intermediates from glycerol to pyruvate (Pyr) were large, whereas those of tricarboxylic acid (TCA) cycle intermediates were much smaller than that of Pyr. Surprisingly, the time constants of intermediates in gluconeogenesis and the pentose phosphate (PP) pathway were large, suggesting that a large amount of the uptaken glycerol was metabolized via the PP pathway. To synthesize fatty acids that make up TAG from acetyl-CoA (AcCoA), 14 molecules of nicotinamide adenine dinucleotide phosphate (NADPH) per C16 fatty acid molecule are required. Because the oxidative PP pathway generates NADPH, this pathway would contribute to supply NADPH for fatty acid synthesis. To confirm that the oxidative PP pathway can supply the NADPH required for TAG production, flux analysis was conducted based on the measured specific rates and mass balances. Flux analysis revealed that the NADPH necessary for TAG production was supplied by metabolizing 48.2% of the uptaken glycerol through gluconeogenesis and the PP pathway. This result was consistent with the result of the ¹³C-labeling experiment. Furthermore, comparison of the actual flux distribution with the ideal flux distribution for TAG production suggested that it is necessary to flow more dihydroxyacetonephosphate (DHAP) through gluconeogenesis to improve TAG yield.

     

Interaction of oxamate with the gluconeogenic pathway in rat liver

Oxamate, a structural analog of pyruvate, known as a potent inhibitor of lactic dehydrogenase, lactic dehydrogenase, produces an inhibition of gluconeogenic flux in isolated perfused rat liver or hepatocyte suspensions from low concentrations of pyruvate (less than 0.5 mM) or substrates yielding pyruvate. The following observations indicate that oxamate inhibits flux through pyruvate carboxylase: accumulation of substrates and decreased concentration of all metabolic intermediates beyond pyruvate; decreased levels of aspartate, glutamate, and alanine; and enhanced ketone body production, which is a sensitive indicator of decreased mitochondrial free oxaloacetate levels. The decreased pyruvate carboxylase flux does not seem to be the result of a direct inhibitory action of oxamate on this enzyme but is secondary to a decreased rate of pyruvate entry into the mitochondria. This assumption is based on the following observations: Above 0.4 mM pyruvate, no significant inhibitory effect of oxamate on gluconeogenesis was observed. The competitive nature of oxamate inhibition is in conflict with its effect on isolated pyruvate carboxylase which is noncompetitive for pyruvate. Fatty acid oxidation was effective in stimulating gluconeogenesis in the presence of oxamate only at concentrations of pyruvate above 0.4 mM. Since only at low pyruvate concentrations its entry into the mitochondria occurs via the monocarboxylate translocator, from these observations it follows that pyruvate transport across the mitochondrial membrane, and not its carboxylation, is the first nonequilibrium step in the gluconeogenic pathway. In the presence of oxamate, fatty acid oxidation inhibited gluconeogenesis from lactate, alanine, and low pyruvate concentrations (less than 0.5 mM), and the rate of transfer of reducing equivalents to the cytosol was significantly decreased. Whether fatty acids stimulate or inhibit gluconeogenesis appears to correlate with the rate of flux through pyruvate carboxylase which ultimately seems to rely on pyruvate availability. Unless adequate rates of oxaloacetate formation are maintained, the shift of the mitochondrial NAD couple to a more reduced state during fatty acid oxidation seems to decrease mitochondrial oxaloacetate resulting in a decreased rate of transfer of carbon and reducing power to the cytosol.