13C n.m.r. isotopomer and computer-simulation studies of the non-oxidative pentose phosphate pathway of human erythrocytes

Constructing and analysing dynamic models with modelbase v1.2.3: a software update

BACKGROUND

Computational mathematical models of biological and biomedical systems have been successfully applied to advance our understanding of various regulatory processes, metabolic fluxes, effects of drug therapies, and disease evolution and transmission. Unfortunately, despite community efforts leading to the development of SBML and the BioModels database, many published models have not been fully exploited, largely due to a lack of proper documentation or the dependence on proprietary software. To facilitate the reuse and further development of systems biology and systems medicine models, an open-source toolbox that makes the overall process of model construction more consistent, understandable, transparent, and reproducible is desired.

RESULTS AND DISCUSSION

We provide an update on the development of modelbase, a free, expandable Python package for constructing and analysing ordinary differential equation-based mathematical models of dynamic systems. It provides intuitive and unified methods to construct and solve these systems. Significantly expanded visualisation methods allow for convenient analysis of the structural and dynamic properties of models. After specifying reaction stoichiometries and rate equations modelbase can automatically assemble the associated system of differential equations. A newly provided library of common kinetic rate laws reduces the repetitiveness of the computer programming code. modelbase is also fully compatible with SBML. Previous versions provided functions for the automatic construction of networks for isotope labelling studies. Now, using user-provided label maps, modelbase v1.2.3 streamlines the expansion of classic models to their isotope-specific versions. Finally, the library of previously published models implemented in modelbase is growing continuously. Ranging from photosynthesis to tumour cell growth to viral infection evolution, all these models are now available in a transparent, reusable and unified format through modelbase.

CONCLUSION

With this new Python software package, which is written in currently one of the most popular programming languages, the user can develop new models and actively profit from the work of others. modelbase enables reproducing and replicating models in a consistent, tractable and expandable manner. Moreover, the expansion of models to their isotopic label-specific versions enables simulating label propagation, thus providing quantitative information regarding network topology and metabolic fluxes.

Metabolic incorporation of H218 O into specific glucose-6-phosphate oxygens by red-blood-cell lysates as observed by 13 C isotope-shifted NMR signals

Water enriched with oxygen-18 (H₂¹⁸ O) is a potential tracer for evaluating the sources of glucose and glycogen synthesis since it is incorporated into specific sites of glucose-6-phosphate via specific enzyme-mediated exchange/addition mechanisms. Unlike ² H, ¹⁸ O does not experience significant isotope effects for any of these processes. Therefore, H₂¹⁸ O might provide more precise estimates of endogenous carbohydrate synthesis compared with deuterated water provided that positional ¹⁸ O enrichments of glucose can be measured. As a proof of concept, H₂¹⁸ O was incorporated into a well characterized hemolysate model of sugar phosphate metabolism and ¹³ C NMR was applied to quantify positional ¹⁸ O enrichment of glucose-6-phosphate oxygens. Human erythrocyte hemolysate preparations were incubated overnight at 37 °C with a buffer containing sugar phosphate precursors and 20% (n = 5) and 80% (n = 1) H₂¹⁸ O. Enrichment of glucose-6-phosphate was analyzed by ¹³ C NMR analysis of ¹⁸ O-shifted versus unshifted signals following derivatization to monoacetone glucose (MAG). ¹³ C NMR MAG spectra from hemolysate revealed resolved ¹⁸ O-shifted signals in Positions 1-5. Mean ¹⁸ O enrichments were 16.4 ± 1.6% (Position 1), 13.3 ± 1.3% (Position 2), 4.1 ± 1.1% (Position 3), 12.6 ± 0.8% (Position 4), 10.7 ± 1.4% (Position 5), and no detectable enrichment of Position 6. No ¹⁸ O-shifted glucose-6-phosphate signals were detected in preparations containing sugar phosphate precursors only. H₂¹⁸ O is incorporated into Positions 1-5 of glucose-6-phosphate in accordance with spontaneous aldose hydration and specific enzymatic reaction mechanisms. This provides a basis for its deployment as a tracer for glucose and glycogen biosynthesis.

The importance of the photosynthetic Gibbs effect in the elucidation of the Calvin-Benson-Bassham cycle

The photosynthetic carbon reduction cycle, or Calvin-Benson-Bassham (CBB) cycle, is now contained in every standard biochemistry textbook. Although the cycle was already proposed in 1954, it is still the subject of intense research, and even the structure of the cycle, i.e. the exact series of reactions, is still under debate. The controversy about the cycle's structure was fuelled by the findings of Gibbs and Kandler in 1956 and 1957, when they observed that radioactive 14CO2 was dynamically incorporated in hexoses in a very atypical and asymmetrical way, a phenomenon later termed the 'photosynthetic Gibbs effect'. Now, it is widely accepted that the photosynthetic Gibbs effect is not in contradiction to the reaction scheme proposed by CBB, but the arguments given have been largely qualitative and hand-waving. To fully appreciate the controversy and to understand the difficulties in interpreting the Gibbs effect, it is illustrative to illuminate the history of the discovery of the CBB cycle. We here give an account of central scientific advances and discoveries, which were essential prerequisites for the elucidation of the cycle. Placing the historic discoveries in the context of the modern textbook pathway scheme illustrates the complexity of the cycle and demonstrates why especially dynamic labelling experiments are far from easy to interpret. We conclude by arguing that it requires sound theoretical approaches to resolve conflicting interpretations and to provide consistent quantitative explanations.

Removing the bottlenecks of cell culture metabolomics: fast normalization procedure, correlation of metabolites to cell number, and impact of the cell harvesting method

INTRODUCTION

Although cultured cells are nowadays regularly analyzed by metabolomics technologies, some issues in study setup and data processing are still not resolved to complete satisfaction: a suitable harvesting method for adherent cells, a fast and robust method for data normalization, and the proof that metabolite levels can be normalized to cell number.

OBJECTIVES

We intended to develop a fast method for normalization of cell culture metabolomics samples, to analyze how metabolite levels correlate with cell numbers, and to elucidate the impact of the kind of harvesting on measured metabolite profiles.

METHODS

We cultured four different human cell lines and used them to develop a fluorescence-based method for DNA quantification. Further, we assessed the correlation between metabolite levels and cell numbers and focused on the impact of the harvesting method (scraping or trypsinization) on the metabolite profile.

RESULTS

We developed a fast, sensitive and robust fluorescence-based method for DNA quantification showing excellent linear correlation between fluorescence intensities and cell numbers for all cell lines. Furthermore, 82-97 % of the measured intracellular metabolites displayed linear correlation between metabolite concentrations and cell numbers. We observed differences in amino acids, biogenic amines, and lipid levels between trypsinized and scraped cells.

CONCLUSION

We offer a fast, robust, and validated normalization method for cell culture metabolomics samples and demonstrate the eligibility of the normalization of metabolomics data to the cell number. We show a cell line and metabolite-specific impact of the harvesting method on metabolite concentrations.

Demonstration of glucose-6-phosphate hydrogen 5 enrichment from deuterated water by transaldolase-mediated exchange alone

PURPOSE

Enrichment of glucose position 5 (H5) from deuterated water ((2)H2O) is widely used for quantifying gluconeogenesis. Exchanges of hexose and triose phosphates mediated by transaldolase have been postulated to enrich H5 independently of gluconeogenesis, but to date this mechanism has not been proven. We determined the enrichment of glucose-6-phosphate (G6P), the immediate precursor of endogenously produced glucose, from (2)H2O in erythrocyte hemolysate preparations. Here, transaldolase exchange is active but gluconeogenesis is absent.

METHODS

Hemolysates were prepared from human erythrocytes and incubated with a buffer containing 5% [U-(13)C]G6P, unlabeled fructose 1,6-bisphosphate, and 10% (2)H2O. G6P (2)H-enrichment and (13)C-isotopomer distributions were analyzed by (2)H and (13)C NMR following derivatization to monoacetone glucose.

RESULTS

(2)H NMR analysis revealed high (2)H-enrichment of G6P hydrogens 2, 4, and 5; low enrichment of hydrogen 3, and residual enrichments of hydrogens 1, 6R, and 6S. (13)C NMR isotopomer analysis revealed that [U-(13)C]G6P was converted to [1,2,3-(13)C3]G6P, a predicted product of transaldolase-mediated exchange, as well as [1,2-(13)C2]G6P and [3-(13)C]G6P, predicted products of combined transaldolase and transketolase exchanges.

CONCLUSION

Hydrogen 5 of G6P was enriched from (2)H2O through exchanges mediated by transaldolase. These studies prove that G6P can be enriched in hydrogen 5 by (2)H2O independently of gluconeogenesis.

Pentose Phosphate Pathway in Disease and Therapy

Pentose phosphate (PP) pathway, which is ubiquitously present in all living organisms, is one of the major metabolic pathways associated with glucose metabolism. The most important functions of this pathway includes the generation of reducing equivalents in the form of NADPH for reductive biosynthesis, and production of ribose sugars for the biosynthesis of nucleotides, amino acids, and other macromolecules required by all living cells. Under normal conditions of growth, PP pathway is important for cell cycle progression, myelin formation, and the maintenance of the structure and function of brain, liver, cortex and other organs. Under diseased conditions, such as in cases of many metabolic, neurological or malignant diseases, pathological mechanisms augment due to defects in the PP pathway genes. Adoption of alternative metabolic pathways by cells that are metabolically abnormal, or malignant cells that are resistant to chemotherapeutic drugs often plays important roles in disease progression and severity. Accordingly, the PP pathway has been suggested to play critical roles in protecting cancer or abnormal cells by providing reduced environment, to protect cells from oxidative damage and generating structural components for nucleic acids biosynthesis. Novel drugs that targets one or more components of the PP pathway could potentially serve to overcome challenges associated with currently available therapeutic options for many metabolic and non-metabolic diseases. However, careful designing of drugs is critical that takes into the accounts of cell’s broader genomic, proteomic and metabolic contexts under consideration, in order to avoid undesirable side-effects. In this review, we discuss the role of PP pathway under normal and abnormal physiological conditions and the potential of the PP pathway as a target for new drug development to treat metabolic and non-metabolic diseases.

Determination of the fluxes in the central metabolism of Corynebacterium glutamicum by nuclear magnetic resonance spectroscopy combined with metabolite balancing

To determine the in vivo fluxes of the central metabolism we have developed a comprehensive approach exclusively based on the fundamental enzyme reactions known to be present, the fate of the carbon atoms of individual reactions, and the metabolite balance of the culture. No information on the energy balance is required, nor information on enzyme activities, or the directionalities of reactions. Our approach combines the power of (1)H-detected (13)C nuclear magnetic resonance spectroscopy to follow individual carbons with the simplicity of establishing carbon balances of bacterial cultures. We grew a lysine-producing strain of Corynebacterium glutamicum to the metabolic and isotopic steady state with [1-(13)C]glucose and determined the fractional enrichments in 27 carbon atoms of 11 amino acids isolated from the cell. Since precursor metabolites of the central metabolism are incorporated in an exactly defined manner in the carbon skeleton of amino acids, the fractional enrichments in carbons of precursor metabolites (oxaloacetate, glyceraldehyde 3-phosphate, erythrose 4-phosphate, etc.) became directly accessible. A concise and generally applicable mathematical model was established using matrix calculus to express all metabolite mass and carbon labeling balances. An appropriate all-purpose software for the iterative solution of the equations is supplied. Applying this comprehensive methodology to C. glutamicum, all major fluxes within the central metabolism were determined. The result is that the flux through the pentose phosphate pathway is 66.4% (relative to the glucose input flux of 1.49 mmol/g dry weight h), that of entry into the tricarboxylic acid cycle 62.2%, and the contribution of the succinylase pathway of lysine synthesis 13.7%. Due to the large amount and high quality of measured data in vivo exchange reactions could also be quantitated with particularly high exchange rates within the pentose phosphate pathway for the ribose 5-phosphate transketolase reaction. Moreover, the total net flux of the anaplerotic reactions was quantitated as 38.0%. Most importantly, we found that in vivo one component within these anaplerotic reactions is a back flux from the carbon 4 units of the tricarboxylic acid cycle to the carbon 3 units of glycolysis of 30.6%. (c) 1996 John Wiley & Sons, Inc.

The shikimate pathway and aromatic amino Acid biosynthesis in plants

L-tryptophan, L-phenylalanine, and L-tyrosine are aromatic amino acids (AAAs) that are used for the synthesis of proteins and that in plants also serve as precursors of numerous natural products, such as pigments, alkaloids, hormones, and cell wall components. All three AAAs are derived from the shikimate pathway, to which ≥30% of photosynthetically fixed carbon is directed in vascular plants. Because their biosynthetic pathways have been lost in animal lineages, the AAAs are essential components of the diets of humans, and the enzymes required for their synthesis have been targeted for the development of herbicides. This review highlights recent molecular identification of enzymes of the pathway and summarizes the pathway organization and the transcriptional/posttranscriptional regulation of the AAA biosynthetic network. It also identifies the current limited knowledge of the subcellular compartmentalization and the metabolite transport involved in the plant AAA pathways and discusses metabolic engineering efforts aimed at improving production of the AAA-derived plant natural products.

Riboneogenesis in yeast

Glucose is catabolized in yeast via two fundamental routes, glycolysis and the oxidative pentose phosphate pathway, which produces NADPH and the essential nucleotide component ribose-5-phosphate. Here, we describe riboneogenesis, a thermodynamically driven pathway that converts glycolytic intermediates into ribose-5-phosphate without production of NADPH. Riboneogenesis begins with synthesis, by the combined action of transketolase and aldolase, of the seven-carbon bisphosphorylated sugar sedoheptulose-1,7-bisphosphate. In the pathway's committed step, sedoheptulose bisphosphate is hydrolyzed to sedoheptulose-7-phosphate by the enzyme sedoheptulose-1,7-bisphosphatase (SHB17), whose activity we identified based on metabolomic analysis of the corresponding knockout strain. The crystal structure of Shb17 in complex with sedoheptulose-1,7-bisphosphate reveals that the substrate binds in the closed furan form in the active site. Sedoheptulose-7-phosphate is ultimately converted by known enzymes of the nonoxidative pentose phosphate pathway to ribose-5-phosphate. Flux through SHB17 increases when ribose demand is high relative to demand for NADPH, including during ribosome biogenesis in metabolically synchronized yeast cells.

Models of the human metabolic network: aiming to reconcile metabolomics and genomics

The metabolic syndrome, inborn errors of metabolism, and drug-induced changes to metabolic states all bring about a seemingly bewildering array of alterations in metabolite concentrations; these often occur in tissues and cells that are distant from those containing the primary biochemical lesion. How is it possible to collect sufficient biochemical information from a patient to enable us to work backwards and pinpoint the primary lesion, and possibly treat it in this whole human metabolic network? Potential analyses have benefited from modern methods such as ultra-high-pressure liquid chromatography, mass spectrometry, nuclear magnetic resonance spectroscopy, and more. A yet greater challenge is the prediction of outcomes of possible modern therapies using drugs and genetic engineering. This exposes the notion of viewing metabolism from a completely different perspective, with focus on the enzymes, regulators, and structural elements that are encoded by genes that specify the amino acid sequences, and hence encode the various interactions, be they regulatory or catalytic. The mainstream view of metabolism is being challenged, so we discuss here the reconciling of traditionally quantitative chemocentric metabolism with the seemingly 'parameter-free' genomic description, and vice versa.

Modeling isotopomer distributions in biochemical networks using isotopomer mapping matrices

Within the last decades NMR spectroscopy has undergone tremendous development and has become a powerful analytical tool for the investigation of intracellular flux distributions in biochemical networks using (13)C-labeled substrates. Not only are the experiments much easier to conduct than experiments employing radioactive tracer elements, but NMR spectroscopy also provides additional information on the labeling pattern of the metabolites. Whereas the maximum amount of information obtainable with (14)C-labeled substrates is the fractional enrichment in the individual carbon atom positions, NMR spectroscopy can also provide information on the degree of labeling at neighboring carbon atom positions by analyzing multiplet patterns in NMR spectra or using 2-dimensional NMR spectra. It is possible to quantify the mole fractions of molecules that show a specific labeling pattern, i.e., information of the isotopomer distribution in metabolite pools can be obtained. The isotopomer distribution is the maximum amount of information that in theory can be obtained from (13)C-tracer studies. The wealth of information contained in NMR spectra frequently leads to overdetermined algebraic systems. Consequently, fluxes must be estimated by nonlinear least squares analysis, in which experimental labeling data is compared with simulated steady state isotopomer distributions. Hence, mathematical models are required to compute the steady state isotopomer distribution as a function of a given set of steady state fluxes. Because 2(n) possible labeling patterns exist in a molecule of n carbon atoms, and each pattern corresponds to a separate state in the isotopomer model, these models are inherently complex. Model complexity, so far, has restricted usage of isotopomer information to relatively small metabolic networks. A general methodology for the formulation of isotopomer models is described. The model complexity of isotopomer models is reduced to that of classical metabolic models by expressing the 2(n) isotopomer mass balances of a metabolite pool in a single matrix equation. Using this approach an isotopomer model has been implemented that describes label distribution in primary carbon metabolism, i.e., in a metabolic network including the Embden-Meyerhof-Parnas and pentose phosphate pathway, the tricarboxylic acid cycle, and selected anaplerotic reaction sequences. The model calculates the steady state label distribution in all metabolite pools as a function of the steady state fluxes and is applied to demonstrate the effect of selected anaplerotic fluxes on the labeling pattern of the pathway intermediates.

Role of band 3 in regulating metabolic flux of red blood cells

Deoxygenation elevates glycolytic flux and lowers pentose phosphate pathway (PPP) activity in mammalian erythrocytes. The membrane anion transport protein (band 3 or AE1) is thought to facilitate this process by binding glycolytic enzymes (GEs) and inhibiting their activity in an oxygen-dependent manner. However, this regulatory mechanism has not been demonstrated under physiological conditions. In this study, we introduce a (1)H-(13)C NMR technique for measuring metabolic fluxes in intact cells. The role of band 3 in mediating the oxygenated/deoxygenated metabolic transition was examined by treating cells with pervanadate, a reagent that prevents the GE-band 3 complex from forming. We report that pervanadate suppresses oxygen-dependent changes in glycolytic and PPP fluxes. Moreover, these metabolic alterations were not attributable to modulation of bisphosphoglycerate mutase, direct inhibition of GEs by pervanadate, or oxidation, which are the major side effects of pervanadate treatment. These data provide direct evidence supporting the role of band 3 in mediating oxygen-regulated metabolic transitions.

Quantitative evaluation of respiration induced metabolic oscillations in erythrocytes

The changes in the partial pressures of oxygen and carbon dioxide (P(O(2)) and P(CO(2))) during blood circulation alter erythrocyte metabolism, hereby causing flux changes between oxygenated and deoxygenated blood. In the study we have modeled this effect by extending the comprehensive kinetic model by Mulquiney and Kuchel [P.J. Mulquiney, and P.W. Kuchel. Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: equations and parameter refinement, Biochem. J. 1999, 342, 581-596.] with a kinetic model of hemoglobin oxy-/deoxygenation transition based on an oxygen dissociation model developed by Dash and Bassingthwaighte [R. Dash, and J. Bassingthwaighte. Blood HbO(2) and HbCO(2) dissociation curves at varied O(2), CO(2), pH, 2,3-DPG and temperature levels, Ann. Biomed. Eng., 2004, 32(12), 1676-1693.]. The system has been studied during transitions from the arterial to the venous phases by simply forcing P(O(2)) and P(CO(2)) to follow the physiological values of venous and arterial blood. The investigations show that the system passively follows a limit cycle driven by the forced oscillations of P(O(2)) and is thus inadequately described solely by steady state consideration. The metabolic system exhibits a broad distribution of time scales. Relaxations of modes with hemoglobin and Mg(2+) binding reactions are very fast, while modes involving glycolytic, membrane transport and 2,3-BPG shunt reactions are much slower. Incomplete slow mode relaxations during the 60 s period of the forced transitions cause significant overshoots of important fluxes and metabolite concentrations - notably ATP, 2,3-BPG, and Mg(2+). The overshoot phenomenon arises in consequence of a periodical forcing and is likely to be widespread in nature - warranting a special consideration for relevant systems.

Isotopomer subspaces as indicators of metabolic-pathway structure

The relative abundances and rates of formation of particular isotopic isomers (isotopomers) of metabolic intermediates from (13)C-labelled substrates in living cells provide information on the routes taken by the initial (13)C-atoms. When a primary substrate such as [U,(13)C] d-glucose is added to human erythrocytes, the pattern of labels in terminal metabolites is determined by a set of carbon-group exchange reactions in both glycolysis and the pentose phosphate pathway (PPP). Of a given terminal metabolite, not all possible isotopomers will be produced from each possible primary substrate isotopomer. There are only 8 different (13)C-isotopomers of lactate but not all of these are produced when one of the 64 possible (13)C-isotopomers of glucose is used as the input substrate; thus a subset of all 63 glucose isotopomers x 8 lactate isotopomers+1 unlabelled glucose x 1 unlabelled lactate=505 pattern associations, would be produced if a complete experimental analysis were performed with all the glucose variants. The pattern of labelling in this isotopomer subspace reflects the nature of the re-ordering reactions that 'direct' the metabolism. Predicting the combinatorial rearrangements for particular sets of reactions and comparing these with real data should enable conclusions to be drawn about which enzymes are involved in the real metabolic system. An example of the glycolysis-PPP system is discussed in the context of a debate that occurred around the F- and L-type PPPs and which one actually operates in the human RBC. As part of this discussion we introduce the term 'combinatorial deficit' of all possible isotopomers and we show that this deficit is less for the F- than the L-type pathway.

Revisiting the 13C-label distribution of the non-oxidative branch of the pentose phosphate pathway based upon kinetic and genetic evidence

The currently applied reaction structure in stoichiometric flux balance models for the nonoxidative branch of the pentose phosphate pathway is not in accordance with the established ping-pong kinetic mechanism of the enzymes transketolase (EC 2.2.1.1) and transaldolase (EC 2.2.1.2). Based upon the ping-pong mechanism, the traditional reactions of the nonoxidative branch of the pentose phosphate pathway are replaced by metabolite specific, reversible, glycolaldehyde moiety (C(2)) and dihydroxyacetone moiety (C(3)) fragments producing and consuming half-reactions. It is shown that a stoichiometric model based upon these half-reactions is fundamentally different from the currently applied stoichiometric models with respect to the number of independent C(2) and C(3) fragment pools in the pentose phosphate pathway and can lead to different label distributions for (13)C-tracer experiments. To investigate the actual impact of the new reaction structure on the estimated flux patterns within a cell, mass isotopomer measurements from a previously published (13)C-based metabolic flux analysis of Saccharomyces cerevisiae were used. Different flux patterns were found. From a genetic point of view, it is well known that several micro-organisms, including Escherichia coli and S. cerevisiae, contain multiple genes encoding isoenzymes of transketolase and transaldolase. However, the extent to which these gene products are also actively expressed remains unknown. It is shown that the newly proposed stoichiometric model allows study of the effect of isoenzymes on the (13)C-label distribution in the nonoxidative branch of the pentose phosphate pathway by extending the half-reaction based stoichiometric model with two distinct transketolase enzymes instead of one. Results show that the inclusion of isoenzymes affects the ensuing flux estimates.

Current status and challenges in connecting models of erythrocyte metabolism to experimental reality

Detailed kinetic models of human erythrocyte metabolism have served to summarize the vast literature and to predict outcomes from laboratory and "Nature's" experiments on this simple cell. Mathematical methods for handling the large array of nonlinear ordinary differential equations that describe the time dependence of this system are well developed, but experimental methods that can guide the evolution of the models are in short supply. NMR spectroscopy is one method that is non-selective with respect to analyte detection but is highly specific with respect to their identification and quantification. Thus time courses of metabolism are readily recorded for easily changed experimental conditions. While the data can be simulated, the systems of equations are too complex to allow solutions of the inverse problem, namely parameter-value estimation for the large number of enzyme and membrane-transport reactions operating in situ as opposed to in vitro. Other complications with the modelling include the dependence of cell volume on time, and the rates of membrane transport processes are often dependent on the membrane potential. These matters are discussed in the light of new modelling strategies.

Carbohydrate cycling in micro-organisms: what can (13)C-NMR tell us?

The extension of (13)C-nuclear magnetic resonance (NMR) techniques to study cellular metabolism over recent years has provided valuable data supporting the occurrence, diversity and extent of carbon cycling in the carbohydrate metabolism of micro-organisms. The occurrence of such cycles, resulting from the simultaneous operation of different and sometimes opposite individual steps, is inherently related to the network organisation of cellular metabolism. These cycles are tentatively classified here as 'reversibility', 'metabolic' and 'substrate' cycles on the basis of their balance in carbon and cofactors. Current hypotheses concerning the physiological relevance of carbohydrate cycles are discussed in light of the (13)C-NMR data. They most likely represent system-level mechanisms for coherent and timely partitioning of carbon resources to fit with the various biosynthetic, energetic or redox needs of cells and/or additional strategies in the adaptive capacity of micro-organisms to face variation in environmental conditions.

Application of radiolabeled tracers to biocatalytic flux analysis

Radiolabeled tracers can provide valuable information about the structure of and flux distributions in biocatalytic reaction networks. This method derives from prior studies of glucose metabolism in mammalian systems and is implemented by pulsing a culture with a radiolabeled metabolite that can be transported into the cells and subsequently measuring the radioactivity of all network metabolites following separation by liquid chromatography. Intracellular fluxes can be directly determined from the transient radioactivity count data by tracking the depletion of the radiolabeled metabolite and/or the accompanying accumulation of any products formed. This technique differs from previous methods in that it is applied within a systems approach to the problem of flux determination. It has been used for the investigation of the indene bioconversion network expressed in Rhodococcus sp. KY1. Flux estimates obtained by radioactive tracers were confirmed by macroscopic metabolite balancing and showed that indene oxidation in steady state chemostat cultures proceeds primarily through a monooxygenase activity forming (1S,2R)-indan oxide, with no dehydrogenation of trans-(1R,2R)-indandiol. These results confirmed the significance of indan oxide formation and identified the hydrolysis of indan oxide as a key step in maximizing the production of (2R)-indandiol, a chiral precursor of the HIV protease inhibitor, Crixivan.

Possible pitfalls of flux calculations based on (13)C-labeling

Metabolic engineers have enthusiastically adopted the (13)C-labeling technique as a powerful tool for elucidating fluxes in metabolic networks. This tracer technique makes it possible to determine fluxes that are unobservable using only metabolite balances and allows the elimination of doubtful cofactor balances that are indispensable in flux analysis based on metabolite balancing alone. The (13)C-labeling technique, however, relies on a number of assumptions that are not free from uncertainties. Two possible errors in the models that are needed to determine the metabolic fluxes from labeling data are omitted reactions and ignored occurrence of channeling. By means of two representative examples it is shown that these modeling errors may lead to serious errors in the calculated flux distributions despite the use of labeling data. A complicating fact is that the model errors are not always easily detected as poor models may still yield good fits of experimental data. Results of (13)C-labeling experiments should therefore be interpreted with appropriate caution.

Metabolic characterization of Lactococcus lactis deficient in lactate dehydrogenase using in vivo 13C-NMR

The metabolism of glucose by nongrowing cells of Lactococcus lactis strain FI7851, constructed from the wild-type L. lactis strain MG1363 by disruption of the lactate dehydrogenase (ldh) gene [Gasson, M.J., Benson, K., Swindel, S. & Griffin, H. (1996) Lait 76, 33-40] was studied in a noninvasive manner by 13C-NMR. The kinetics of the build-up and consumption of the pools of intracellular intermediates mannitol 1-phosphate, fructose 1,6-bisphosphate, 3-phosphoglycerate, and phosphoenolpyruvate as well as the utilization of [1-13C]glucose and formation of products (lactate, acetate, mannitol, ethanol, acetoin, 2,3-butanediol) were monitored in vivo with a time resolution of 30 s. The metabolism of glucose by the parental wild-type strain was also examined for comparison. A clear shift from typical homolactic fermentation (parental strain) to a mixed acid fermentation (lactate dehdydrogenase deficient; LDHd strain) was observed. Furthermore, high levels of mannitol were transiently produced and metabolized once glucose was depleted. Mannitol 1-phosphate accumulated intracellularly up to 76 mM concentration. Mannitol was formed from fructose 6-phosphate by the combined action of mannitol-1-phosphate dehydrogenase and phosphatase. The results show that the formation of mannitol 1-phosphate by the LDHd strain during glucose catabolism is a consequence of impairment in NADH oxidation caused by a highly reduced LDH activity, the transient production of mannitol 1-phosphate serving as a regeneration pathway for NAD+ regeneration. Oxygen availability caused a drastic change in the pattern of intermediates and end-products, reinforcing the key-role of the fulfilment of the redox balance. The flux control coefficients for the step catalysed by mannitol-1-phosphate dehydrogenase were calculated and the implications in the design of metabolic engineering strategies are discussed.

Relevance and isotopic assessment of hexose-6-phosphate recycling in micro-organisms

Some pathways of hexose-6-phosphate recycling--those involving a breakdown of the hexose skeleton--through carbohydrate metabolism of micro-organisms were analyzed for both metabolic and isotopic effects. Two modes of recycling were proposed based on the degree of alteration of the hexose molecule through the catabolic part of the cycle. Simulated operation of most of these pathways resulted in increased synthesis of hexose-6-phosphate and NADPH, and reduced the NADH and moreover the ATP synthesis within the carbohydrate metabolism. A basic model for the quantitative assessment by means of isotopic studies of the processes of hexose-6-phosphate recycling is presented. The model was initially designed for the study of micro-organisms producing polysaccharides, but it can be extended to other situations.

Cyclic organization of the carbohydrate metabolism in Sinorhizobium meliloti

The pathways of polysaccharide biosynthesis were investigated in cells of Sinorhizobium meliloti (strain Su47) using a stable isotope approach. The isotopic labeling of the periplasmic beta-1,2-glucans synthesized from glucose labeled at various positions evidenced the involvement of catabolic pathways, namely the pentose-phosphate and Entner-Doudoroff pathways, into the early steps of polysaccharide synthesis. The exopolysaccharides produced at the same time had a labeling pattern similar to that of the beta-glucans, indicating similar early steps for both polysaccharides. The results emphasized a cyclic organization of the carbohydrate metabolism in S. meliloti, in which the carbons of the initial hexose were allowed to re-enter the catabolic pathways many times. The metabolic incidences of such metabolic topology are discussed.

Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: in vivo kinetic characterization of 2,3-bisphosphoglycerate synthase/phosphatase using 13C and 31P NMR

This is the first in a series of three papers [see also Mulquiney and Kuchel (1999) Biochem. J. 342, 579-594; Mulquiney and Kuchel (1999) Biochem. J. 342, 595-602] that present a detailed mathematical model of erythrocyte metabolism which explains the regulation and control of 2,3-bisphosphoglycerate (2,3-BPG) metabolism. 2,3-BPG is a modulator of haemoglobin oxygen affinity and hence plays an important role in blood oxygen transport and delivery. This paper presents an in vivo kinetic characterization of 2,3-BPG synthase/phosphatase (BPGS/P), the enzyme that catalyses both the synthesis and degradation of 2,3-BPG. Much previous work had indicated that the behaviour of this enzyme in vitro is markedly different from that in vivo. (13)C and (31)P NMR were used to monitor the time courses of selected metabolites when erythrocytes were incubated with or without [U-(13)C]glucose. Simulations of the experimental time courses were then made. By iteratively changing the parameters of the BPGS/P part of the model until a good match between the NMR-derived data and simulations were achieved, it was possible to characterize BPGS/P kinetically in vivo. This work revealed that: (1) the pH-dependence of the synthase activity results largely from a strong co-operative inhibition of the synthase activity by protons; (2) 3-phosphoglycerate and 2-phosphoglycerate are much weaker inhibitors of 2,3-BPG phosphatase in vivo than in vitro; (3) the K(m) of BPGS/P for 2,3-BPG is significantly higher than that measured in vitro; (4) the maximal activity of the phosphatase in vivo is approximately twice that in vitro, when P(i) is the sole activator (second substrate); and (5) 2-phosphoglycollate appears to play no role in the activation of the phosphatase in vivo. Using the newly determined kinetic parameters, the percentage of glycolytic carbon flux that passes through the 2, 3-BPG shunt in the normal in vivo steady state was estimated to be 19%.

Quantitative determination of the main glucose metabolic fluxes in human erythrocytes by 13C- and 1H-MR spectroscopy

Information displayed by homonuclear and heteronuclear spin-coupling patterns in 13C- and 1H-MR spectra allowed us to identify the major lactate isotopomers produced either from [1-(13)C]-glucose or from [2-(13)C]-glucose by human erythrocytes. Relative concentrations of detectable isotopomers were determined by integrating the corresponding MR signals. The interpretation of these data in terms of the fractional glucose metabolised through glycolysis and pentose phosphate pathway was performed by a computer simulation of the metabolism that took into account metabolic schemes pertaining to glycolysis and to the F-type of pentose phosphate pathway. The simulation was organised in a way to anticipate the populations of the isotopomers produced from any precursor at a priori established metabolic steady state. By the simulation, isotopomer populations were determined according to different values of pentose cycle, defined as the flux of glyceraldehyde 3-phosphate originating from pentose phosphate pathway at unitary glucose uptake. The populations of the isotopomers originating from [2-(13)C]-glucose were described by polynomials, and ratios between the polynomials were used in conjunction with 13C- and 1H-MR data to determine pentose cycle values. The knowledge of glucose uptake and of pentose cycle value allowed us to perform accurate measurement of the pentose phosphate pathway flux, of the hexokinase and phosphofructokinase fluxes as well as, indirectly, of the carbon dioxide production.

13C-NMR studies of transmembrane electron transfer to extracellular ferricyanide in human erythrocytes

Human erythrocytes are known to reduce ferricyanide (hexacyanoferrate) [Fe(CN)6]3- to ferrocyanide [Fe(CN)6]2- in an extracellular reaction that involves the transmembrane transfer of reducing equivalents; potentially these could be either electrons from NADH, formed in glycolysis inside the cells or transmembrane exchange of reduced solutes. The 13C-NMR resonance of [Fe(13CN)6]3- (which was synthesised in our laboratory) was seen to be very broad while that of ferrocyanide was narrow. This phenomenon formed the basis of a simple non-invasive procedure to study ferricyanide reduction in high-haematocrit suspensions of erythrocytes. The method should be directly applicable to other cell types. In a series of experiments, erythrocyte metabolism was studied in the presence of ferricyanide, using 1H, 13C, and 31P NMR spectroscopy. Incubating the cells with 13C-labelled glucose enabled the rate of ferricyanide reduction, glucose utilisation, and lactate and bicarbonate production to be measured simultaneously. Various metabolic states were imposed as follows: glycolysis was inhibited with F- and iodoacetate; glucose transport was inhibited with phloretin and cytochalasin B; and anion transport was inhibited with dinitrostilbene 2,2'-disulfonate and p-chloromercuriphenyl sulfonate. Earlier work was confirmed, showing that ascorbate is intimately involved in the reduction reaction; but its main action appears not to be mediated by membrane transport but in a membrane-associated redox-protein complex that is functionally linked to glycolysis. Also, large differences (factors of three) in the rate of the reduction reaction were recorded in erythrocytes from different, apparently healthy, donors.

Model of the pH-dependence of the concentrations of complexes involving metabolites, haemoglobin and magnesium ions in the human erythrocyte

The rate of glucose consumption and the concentrations of glycolytic intermediates in human erythrocytes have long been known to be pH sensitive. Despite the extensive literature on modelling erythrocyte metabolism, no model developed so far can adequately describe all of these pH-dependent changes. None of these models have included all the significant association reactions between metabolites, Hb and Mg2+ that will influence metabolism. As part of a larger enterprise to develop a detailed model of erythrocyte glycolysis, we present a sub-model which predicts, as a function of pH and oxygenation state, the concentrations of free and Mg(2+)-bound metabolites that are substrates, co-factors and effectors of glycolysis. This model shows that pH changes around physiological values can cause large changes in the distribution of metabolites between free, bound and Mg(2+)-complexed forms, based on binding interactions alone; in oxygenated cells, at pH 7.2-7.6, many glycolytic intermediates undergo changes in concentration of 50-100%. The model also predicts intracellular concentrations of free Mg2+ in erythrocytes to be 0.4 mM and 0.64 mM in oxygenated and deoxygenated cells, respectively, assuming a total magnesium concentration of 3 mM (approximately 88% of the total magnesium usually found in erythrocytes). This is in close agreement with the values found by Flatman [Flatman, P. W. (1980) J. Physiol. 300, 19-30] and the finding by Flatman and Lew [Flatman, P. & Lew, V. L. (1977) Nature 267, 360-362] that the main Mg2+ buffer systems bind approximately 90% of Mg2+ in the cell. Hexokinase has a high 'flux control coefficient' in human erythrocyte glycolysis, so the dependence of its rate on the pH and oxygenation state of haemoglobin is important. With a low oxygen tension and an intracellular pH of 7.34, the major inhibitor of its activity (2,3-bisphosphoglycerate) is 85% bound to either haemoglobin or Mg2+, and the maximum possible flux of substrate via it would be 2.05 mmol L erythrocytes-1 h-1. However, if the haemoglobin were saturated with oxygen, and the pH were 7.2, it was calculated that the maximum rate would be 1.48 mmol L erythrocytes-1 h-1; this is primarily due to a doubling of the free 2,3-bisphosphoglycerate concentration. However, the full extent of the inhibition is counteracted because the concentration of the Mg(2+)-2,3-bisphosphoglycerate would be approximately doubled. Many other similar comparisons are possible with this new model, which highlights the complex network of interactions between haemoglobin, Mg2+, H+ and the metabolites as substrates and effectors of the glycolytic reactions.

Nuclear magnetic resonance analysis of cell metabolism

Nuclear magnetic resonance (NMR) continues to be a useful tool for the study of cellular metabolism. A variety of NMR techniques have been developed or newly applied to the analysis of cell systems. Many of these techniques are particularly useful for the analysis of immobilized cell bioreactors. The use of several NMR techniques has been an integral part of recent comprehensive metabolic studies. Novel computer-based models and methods have been developed which may make NMR study of metabolism more accessible and powerful.