A transient isotopic labeling methodology for 13C metabolic flux analysis of photoautotrophic microorganisms

Metabolic flux analysis is increasingly recognized as an integral component of systems biology. However, techniques for experimental measurement of system-wide metabolic fluxes in purely photoautotrophic systems (growing on CO(2) as the sole carbon source) have not yet been developed due to the unique problems posed by such systems. In this paper, we demonstrate that an approach that balances positional isotopic distributions transiently is the only route to obtaining system-wide metabolic flux maps for purely autotrophic metabolism. The outlined transient (13)C-MFA methodology enables measurement of fluxes at a metabolic steady-state, while following changes in (13)C-labeling patterns of metabolic intermediates as a function of time, in response to a step-change in (13)C-label input. We use mathematical modeling of the transient isotopic labeling patterns of central intermediates to assess various experimental requirements for photoautotrophic MFA. This includes the need for intracellular metabolite concentration measurements and isotopic labeling measurements as a function of time. We also discuss photobioreactor design and operation in order to measure fluxes under precise environmental conditions. The transient MFA technique can be used to measure and compare fluxes under different conditions of light intensity, nitrogen sources or compare strains with various mutations or gene deletions and additions.     

Detection of metabolic fluxes of o and h atoms into intracellular water in Mammalian cells

Metabolic processes result in the release and exchange of H and O atoms from organic material as well as some inorganic salts and gases. These fluxes of H and O atoms into intracellular water result in an isotopic gradient that can be measured experimentally. Using isotope ratio mass spectroscopy, we revealed that slightly over 50% of the H and O atoms in the intracellular water of exponentially-growing cultured Rat-1 fibroblasts were isotopically distinct from growth medium water. We then employed infrared spectromicroscopy to detect in real time the flux of H atoms in these same cells. Importantly, both of these techniques indicate that the H and O fluxes are dependent on metabolic processes; cells that are in lag phase or are quiescent exhibit a much smaller flux. In addition, water extracted from the muscle tissue of rats contained a population of H and O atoms that were isotopically distinct from body water, consistent with the results obtained using the cultured Rat-1 fibroblasts. Together these data demonstrate that metabolic processes produce fluxes of H and O atoms into intracellular water, and that these fluxes can be detected and measured in both cultured mammalian cells and in mammalian tissue.     

Growth Media Simulating Ileal and Colonic Environments Affect the Intracellular Proteome and Carbon Fluxes of Enterohemorrhagic Escherichia coli O157:H7 Strain EDL933

In this study, the intracellular proteome of Escherichia coli O157:H7 strain EDL933 was analyzed by two-dimensional gel electrophoresis and matrix-assisted laser desorption ionization–time-of-flight (MALDI-TOF) spectrometry after growth in simulated ileal environment media (SIEM) and simulated colonic environment media (SCEM) under aerobic and microaerobic conditions. Differentially expressed intracellular proteins were identified and allocated to functional protein groups. Moreover, metabolic fluxes were analyzed by isotopologue profiling with [U-¹³C₆]glucose as a tracer. The results of this study show that EDL933 responds with differential expression of a complex network of proteins and metabolic pathways, reflecting the high metabolic adaptability of the strain. Growth in SIEM and SCEM is obviously facilitated by the upregulation of nucleotide biosynthesis pathway proteins and could be impaired by exposition to 50 µM 6-mercaptopurine under aerobic conditions. Notably, various stress and virulence factors, including Shiga toxin, were expressed without having contact with a human host.     

Network identification and flux quantification of glucose metabolism in Rhodobacter sphaeroides under photoheterotrophic H(2)-producing conditions

The nonsulfur purple bacteria that exhibit unusual metabolic versatility can produce hydrogen gas (H(2)) using the electrons derived from metabolism of organic compounds during photoheterotrophic growth. Here, based on (13)C tracer experiments, we identified the network of glucose metabolism and quantified intracellular carbon fluxes in Rhodobacter sphaeroides KD131 grown under H(2)-producing conditions. Moreover, we investigated how the intracellular fluxes in R. sphaeroides responded to knockout mutations in hydrogenase and poly-β-hydroxybutyrate synthase genes, which led to increased H(2) yield. The relative contribution of the Entner-Doudoroff pathway and Calvin-Benson-Bassham cycle to glucose metabolism differed significantly in hydrogenase-deficient mutants, and this flux change contributed to the increased formation of the redox equivalent NADH. Disruption of hydrogenase and poly-β-hydroxybutyrate synthase resulted in a significantly increased flux through the phosphoenolpyruvate carboxykinase and a reduced flux through the malic enzyme. A remarkable increase in the flux through the tricarboxylic acid cycle, a major NADH producer, was observed for the mutant strains. The in vivo regulation of the tricarboxylic acid cycle flux in photoheterotrophic R. sphaeroides was discussed based on the measurements of in vitro enzyme activities and intracellular concentrations of NADH and NAD(+). Overall, our results provide quantitative insights into how photoheterotrophic cells manipulate the metabolic network and redistribute intracellular fluxes to generate more electrons for increased H(2) production.     

Metabolic fluxes and metabolic engineering

Metabolic engineering is the directed improvement of cellular properties through the modification of specific biochemical reactions or the introduction of new ones, with the use of recombinant DNA technology. As such, metabolic engineering emphasizes metabolic pathway integration and relies on metabolic fluxes as determinants of cell physiology and measures of metabolic control. The combination of analytical methods to quantify fluxes and their control with molecular biological techniques to implement genetic modifications is the essence of metabolic engineering. Strategies for metabolic flux determination are reviewed in this paper and it is shown how metabolic fluxes can be used in the systematic elucidation of metabolic control in the framework of reaction grouping and top-down metabolic control analysis.     

Defining hypoxia: a systems view of VO2, glycolysis, energetics, and intracellular PO2

The necessity for defining hypoxia as O2-limited energy flux rather than low partial pressure is explored from a systems perspective. Oxidative phosphorylation, the Krebs cycle, glycolysis, substrate supply, and cell energetics interact as subsystems; the set point is a match between ATP demand and aerobic ATP production. To this end the transport subsystem must match the transcapillary and mitochondrial O2 fluxes. High transcapillary O2 flux requires intracellular PO2 in the range 1-10 Torr. In this range the O2 drive on electron transport must be compensated by adaptive changes in the phosphorylation and redox drives. Thus the metabolic subsystem supports diffusive O2 transport by maintaining O2 flux at intracellular partial pressures required for O2 release from blood. Since responses to stress are distributed according to the state of the entire system, several simultaneous metabolic measurements, including intracellular PO2 (or a known direction of change in intracellular PO2) and the O2 dependence of a measurable function are required to judge the adequacy of O2 supply. ATP demand and aerobic capacity must also be evaluated, because the hypoxic threshold depends on the ratio of ATP demand to aerobic capacity. The application and limitation of commonly used criteria of hypoxia are discussed, and a more precise terminology is proposed.     

A single-cell technique for the measurement of membrane potential, membrane conductance, and the efflux of rapidly penetrating solutes in Amphiuma erythrocytes

We describe a single-cell technique for measuring membrane potential, membrane resistance, and the efflux of rapidly penetrating solutes such as Cl and H2O. Erythrocytes from Amphiuma means were aspirated into a Sylgard (Dow Corning Corp.)-coated capillary. The aspirated cell separated a solution within the capillary from a solution in the bath. Each of these two solutions was contiguous with approximately 5% of the total membrane surface. Microelectrodes placed concentrically within the capillary permit the measurement of intracellular voltage, specific membrane resistance, and the electrical seal between the two solutions. The intracellular voltage averaged -17.7 mV (pH 7.6) and changed as either intra- or extracellular chloride was varied. The average specific membrane resistance measured by passing current across the exposed membrane surface was 110 ohm-cm2. 36Cl and tritiated H2O fluxes (0.84 +/- 0.05 x 10(-6) M . cm-2 . min-1 and 6.4 +/- 1.5 x 10(-3) M . cm-2 . min-1, respectively) were determined by noting the rate at which isotope leaves the cell and crosses the membrane exposed to the bath. Our measured values for the flux of 36Cl and tritiated H2O approximate reported values for free-floating cells. 36Cl efflux, in addition, is inhibited by 4-acetamido-4'-isothiocyano-stilbene 2,2'-disulfonic acid (SITS) and furosemide, known inhibitors of the anion exchange mechanism responsible for the rapid anion fluxes of red blood cells. One can also demonstrate directly that > 89% of 36Cl efflux is "electrically silent" by analyzing the flux in the presence of an imposed transcellular voltage.     

Metabolic flux analysis of CHO cell metabolism in the late non‐growth phase

Chinese hamster ovary (CHO) cell cultures are commonly used for production of recombinant human therapeutic proteins. Often the goal of such a process is to separate the growth phase of the cells, from the non‐growth phase where ideally the cells are diverting resources to produce the protein of interest. Characterizing the way that the cells use nutrients in terms of metabolic fluxes as a function of culture conditions can provide a deeper understanding of the cell biology offering guidance for process improvements. To evaluate the fluxes, metabolic flux analysis of the CHO cell culture in the non‐growth phase was performed by a combination of steady‐state isotopomer balancing and stoichiometric modeling. Analysis of the glycolytic pathway and pentose phosphate pathway (PPP) indicated that almost all of the consumed glucose is diverted towards PPP with a high NADPH production; with even recycle from PPP to G6P in some cases. Almost all of the pyruvate produced from glycolysis entered the TCA cycle with little or no lactate production. Comparison of the non‐growth phase against previously reported fluxes from growth phase cultures indicated marked differences in the fluxes, in terms of the split between glycolysis and PPP, and also around the pyruvate node. Possible reasons for the high NADPH production are also discussed. Evaluation of the fluxes indicated that the medium strength, carbon dioxide level, and temperature with dissolved oxygen have statistically significant impacts on different nodes of the flux network. Biotechnol. Bioeng. 2011; 108:82–92. © 2010 Wiley Periodicals, Inc.     

Dynamic flux cartography of hairy roots primary metabolism

A dynamic model for plant cell and hairy root primary metabolism is presented. The model includes nutrient uptake (Pi, sugars, nitrogen sources), the glycolysis and pentose phosphate pathways, the TCA cycle, amino acid biosynthesis, respiratory chain, biosynthesis of cell building blocks (structural hexoses, organic acids, lipids, and organic phosphated molecules). The energy shuttles (ATP, ADP) and cofactors (NAD/H, NADP/H) are also included. The model describes the kinetics of 44 biochemical reactions (fluxes) of the primary metabolism of plant cells and includes 41 biochemical species (metabolites, nutrients, biomass components). Multiple Michaelis-Menten type kinetics are used to describe biochemical reaction rates. Known regulatory phenomena on metabolic pathways are included using sigmoid switch functions. A visualization framework showing fluxes and metabolite concentrations over time is presented. The visualization of fluxes and metabolites is used to analyze simulation results from Catharanthus roseus hairy root 50 d batch cultures. The visualization of the metabolic system allows analyzing split ratios between pathways and flux time-variations. For carbon metabolism, the cells were observed to have relatively high and stable fluxes for the central carbon metabolism and low and variable fluxes for anabolic pathways. For phosphate metabolism, a very high free intracellular Pi turnover rate was observed with higher flux variations than for the carbon metabolism. Nitrogen metabolism also exhibited large flux variations. The potential uses of the model are also discussed.     

Metabolic-flux analysis of continuously cultured hybridoma cells using (13)CO(2) mass spectrometry in combination with (13)C-lactate nuclear magnetic resonance spectroscopy and metabolite balancing

Protein production of mammalian-cell culture is limited due to accumulation of waste products such as lactate, CO(2), and ammonia. In this study, the intracellular fluxes of hybridoma cells are measured to determine the amount by which various metabolic pathways contribute to the secretion of waste products derived from glucose. Continuously cultured hybridoma cells are grown in medium containing either 1-(13)C-, 2-(13)C-, or 6-(13)C-glucose. The uptake and production rates of amino acids, glucose, ammonia, O(2), and CO(2) as well as the cellular composition are measured. In addition, the (13)C distribution of the lactate produced and alanine produced by the hybridomas is determined by (1)H-NMR spectroscopy, and the (13)CO(2)/(12)CO(2) ratio is measured by on-line mass spectrometry. These data are used to calculate the intracellular fluxes of the glycolysis, the pentose phosphate pathway, the TCA cycle, and fluxes involved in amino acid metabolism. It is shown that: (i) approximately 20% of the glucose consumed is channeled through the pentose shunt; (ii) the glycolysis pathway contributes the most to lactate production, and most of the CO(2) is produced by the TCA cycle; (iii) the pyruvate-carboxylase flux is negligibly small; and (iv) the malic-enzyme flux is estimated to be 10% of the glucose uptake rate. Based on these flux data suggestions are made to engineer a more efficient glucose metabolism in mammalian cells.     

Melatonin directly scavenges hydrogen peroxide: a potentially new metabolic pathway of melatonin biotransformation

A potential new metabolic pathway of melatonin biotransformation is described in this investigation. Melatonin was found to directly scavenge hydrogen peroxide (H(2)O(2)) to form N(1)-acetyl-N(2)-formyl-5-methoxykynuramine and, thereafter this compound could be enzymatically converted to N(1)-acetyl-5-methoxykynuramine by catalase. The structures of these kynuramines were identified using proton nuclear magnetic resonance, carbon nuclear magnetic resonance, and mass spectrometry. This is the first report to reveal a possible physiological association between melatonin, H(2)O(2), catalase, and kynuramines. Melatonin scavenges H(2)O(2) in a concentration-dependent manner. This reaction appears to exhibit two distinguishable phases. In the rapid reaction phase, the interaction between melatonin and H(2)O(2) reaches equilibrium rapidly (within 5 s). The rate constant for this phase was calculated to be 2.3 x 10(6) M(-1)s(-1). Thereafter, the relative equilibrium of melatonin and H(2)O(2) was sustained for roughly 1 h, at which time the content of H(2)O(2) decreased gradually over a several hour period, identified as the slow reaction phase. These observations suggest that melatonin, a ubiquitously distributed small nonenzymatic molecule, might serve to directly detoxify H(2)O(2) in living organisms. H(2)O(2) and melatonin are present in all subcellular compartments; thus, presumably, one important function of melatonin may be complementary in function to catalase and glutathione peroxidase in keeping intracellular H(2)O(2) concentrations at steady-state levels.     

Numerical bias estimation for mass spectrometric mass isotopomer analysis

Mass spectrometric (MS) isotopomer analysis has become a standard tool for investigating biological systems using stable isotopes. In particular, metabolic flux analysis uses mass isotopomers of metabolic products typically formed from ¹³C-labeled substrates to quantitate intracellular pathway fluxes. In the current work, we describe a model-driven method of numerical bias estimation regarding MS isotopomer analysis. Correct bias estimation is crucial for measuring statistical qualities of measurements and obtaining reliable fluxes. The model we developed for bias estimation corrects a priori unknown systematic errors unique for each individual mass isotopomer peak. For validation, we carried out both computational simulations and experimental measurements. From stochastic simulations, it was observed that carbon mass isotopomer distributions and measurement noise can be determined much more precisely only if signals are corrected for possible systematic errors. By removing the estimated background signals, the residuals resulting from experimental measurement and model expectation became consistent with normality, experimental variability was reduced, and data consistency was improved. The method is useful for obtaining systematic error-free data from ¹³C tracer experiments and can also be extended to other stable isotopes. As a result, the reliability of metabolic fluxes that are typically computed from mass isotopomer measurements is increased.     

Stable isotope compounds - production,detection, and application

Stable isotopes are used in wide fields of application from natural tracers in biology, geology and archeology through studies of metabolic fluxes to their application as tracers in quantitative proteomics and structural biology. We review the use of stable isotopes of biogenic elements (H, C, N, O, S, Mg, Se) with the emphasis on hydrogen and its heavy isotope deuterium. We will discuss the limitations of enriching various compounds in stable isotopes when produced in living organisms. Finally, we overview methods for measuring stable isotopes, focusing on methods for detection in single cells in situ and their exploitation in modern biotechnologies.     

Metabolic flux analysis of CHO cells at growth and non-growth phases using isotopic tracers and mass spectrometry

Chinese hamster ovary (CHO) cells are the main platform for production of biotherapeutics in the biopharmaceutical industry. However, relatively little is known about the metabolism of CHO cells in cell culture. In this work, metabolism of CHO cells was studied at the growth phase and early stationary phase using isotopic tracers and mass spectrometry. CHO cells were grown in fed-batch culture over a period of six days. On days 2 and 4, [1,2-¹³C] glucose was introduced and the labeling of intracellular metabolites was measured by gas chromatography-mass spectrometry (GC–MS) at 6, 12 and 24 h following the introduction of tracer. Intracellular metabolic fluxes were quantified from measured extracellular rates and ¹³C-labeling dynamics of intracellular metabolites using non-stationary ¹³C-metabolic flux analysis (¹³C-MFA). The flux results revealed significant rewiring of intracellular metabolic fluxes in the transition from growth to non-growth, including changes in energy metabolism, redox metabolism, oxidative pentose phosphate pathway and anaplerosis. At the exponential phase, CHO cell metabolism was characterized by a high flux of glycolysis from glucose to lactate, anaplerosis from pyruvate to oxaloacetate and from glutamate to α-ketoglutarate, and cataplerosis though malic enzyme. At the stationary phase, the flux map was characterized by a reduced flux of glycolysis, net lactate uptake, oxidative pentose phosphate pathway flux, and reduced rate of anaplerosis. The fluxes of pyruvate dehydrogenase and TCA cycle were similar at the exponential and stationary phase. The results presented here provide a solid foundation for future studies of CHO cell metabolism for applications such as cell line development and medium optimization for high-titer production of recombinant proteins.     

Metabolic flux analysis of hybridoma cells in different culture media using mass balances

The estimation of the intracellular fluxes of mammalian cells using only the mass balances of the relevant metabolites is not possible because the set of linear equations defined by these mass balances is underdetermined. Either additional experimental flux data or additional theoretical constraints are required to find one unique flux distribution out of the solution space that is bound by the mass balances. Here, a method is developed using the latter approach. The uptake and production rates of amino acids, glucose, lactate, O(2), CO(2), NH(4), MAB, and the intracellular amino acid pools have been determined for two different steady-states. The cellular composition {total protein and protein composition, total lipids and fatty acid distribution, total carbohydrates, DNA and RNA} has been measured to calculate the requirements for biosynthesis. It is shown to be essential to determine the uptake/production rates of ammonia and either carbon dioxide or oxygen. In mammalian cells these are cometabolites of cyclic metabolic pathways. The flux distribution that is found using the Euclidean minimum norm as the additional theoretical constraint and taking either the CO(2) or the NAD(P)H mass balance into account is shown to be in agreement with the measured O(2) and CO(2) metabolic rates.The metabolic fluxes in hybridoma cells in continuous culture at a specific growth rate of 0.83 day(-1) are estimated for a medium with (optimal medium) and without (suboptimal medium) Primatone RL, an enzymatic hydrolysate of animal tissue that causes a more than twofold increase in cell density. It is concluded that (i)The majority of the consumed glucose (>90%) is channeled through the pentose-phosphate pathway in rapidly proliferating cells.(ii)Pyruvate oxidation and tricarboxylic acid (TCA) cycle activity are relatively low, i.e., 8% of the glucose uptake in suboptimal and 14% in optimal medium, respectively. Under both conditions, only a small fraction of pyruvate is further oxidized to CO(2).(iii)The flux from glutamate to alpha-ketoglutarate (catalyzed by glutamate dehydrogenase) is almost zero in medium with and even slightly reversed in medium without Primatone RL. Almost all glutamate enters the TCA cycle due to the action of transaminases.(iv)Transhydrogenation plays a significant role in hybridoma cells under our experimental conditions. NADPH is produced at relatively high rates (11 x 10(-12) to 13 x 10(-12) mol . cell(-1) . day(-1)) compared to other fluxes in both culture media. (c) 1996 John Wiley & Sons, Inc.     

A genome-scale metabolic reconstruction of Pseudomonas putida KT2440: iJN746 as a cell factory

Background

The study of biological interaction networks is a central theme of systems biology. Here, we investigate the relationships between two distinct types of interaction networks: the metabolic pathway map and the protein-protein interaction network (PIN). It has long been established that successive enzymatic steps are often catalyzed by physically interacting proteins forming permanent or transient multi-enzymes complexes. Inspecting high-throughput PIN data, it was shown recently that, indeed, enzymes involved in successive reactions are generally more likely to interact than other protein pairs. In our study, we expanded this line of research to include comparisons of the underlying respective network topologies as well as to investigate whether the spatial organization of enzyme interactions correlates with metabolic efficiency.

Results

Analyzing yeast data, we detected long-range correlations between shortest paths between proteins in both network types suggesting a mutual correspondence of both network architectures. We discovered that the organizing principles of physical interactions between metabolic enzymes differ from the general PIN of all proteins. While physical interactions between proteins are generally dissortative, enzyme interactions were observed to be assortative. Thus, enzymes frequently interact with other enzymes of similar rather than different degree. Enzymes carrying high flux loads are more likely to physically interact than enzymes with lower metabolic throughput. In particular, enzymes associated with catabolic pathways as well as enzymes involved in the biosynthesis of complex molecules were found to exhibit high degrees of physical clustering. Single proteins were identified that connect major components of the cellular metabolism and may thus be essential for the structural integrity of several biosynthetic systems.

Conclusion

Our results reveal topological equivalences between the protein interaction network and the metabolic pathway network. Evolved protein interactions may contribute significantly towards increasing the efficiency of metabolic processes by permitting higher metabolic fluxes. Thus, our results shed further light on the unifying principles shaping the evolution of both the functional (metabolic) as well as the physical interaction network.     

Metabolic flux and metabolic network analysis of Penicillium chrysogenum using 2D [13C, 1H] COSY NMR measurements and cumulative bondomer simulation

At present two alternative methods are available for analyzing the fluxes in a metabolic network: (1) combining measurements of net conversion rates with a set of metabolite balances including the cofactor balances, or (2) leaving out the cofactor balances and fitting the resulting free fluxes to measured (13)C-labeling data. In this study these two approaches are applied to the fluxes in the glycolysis and pentose phosphate pathway of Penicillium chrysogenum growing on either ammonia or nitrate as the nitrogen source, which is expected to give different pentose phosphate pathway fluxes. The presented flux analyses are based on extensive sets of 2D [(13)C, (1)H] COSY data. A new concept is applied for simulation of this type of (13)C-labeling data: cumulative bondomer modeling. The outcomes of the (13)C-labeling based flux analysis substantially differ from those of the pure metabolite balancing approach. The fluxes that are determined using (13)C-labeling data are shown to be highly dependent on the chosen metabolic network. Extending the traditional nonoxidative pentose phosphate pathway with additional transketolase and transaldolase reactions, extending the glycolysis with a fructose 6-phosphate aldolase/dihydroxyacetone kinase reaction sequence or adding a phosphoenolpyruvate carboxykinase reaction to the model considerably improves the fit of the measured and the simulated NMR data. The results obtained using the extended version of the nonoxidative pentose phosphate pathway model show that the transketolase and transaldolase reactions need not be assumed reversible to get a good fit of the (13)C-labeling data. Strict statistical testing of the outcomes of (13)C-labeling based flux analysis using realistic measurement errors is demonstrated to be of prime importance for verifying the assumed metabolic model.     

in silico analyses of metabolic pathway and protein interaction network for identification of next gen therapeutic targets in Chlamydophila pneumoniae

Chlamydophila pneumoniae, the causative agent of chronic obstructive pulmonary disease (COPD), is presently the fifth mortality causing chronic disease in the world. The understanding of disease and treatment options are limited represents a severe concern and a need for better therapeutics. With the advancements in the field of complete genome sequencing and computational approaches development have lead to metabolic pathway analysis and protein-protein interaction network which provides vital evidence to the protein function and has been appropriate to the fields such as systems biology and drug discovery. Protein interaction network analysis allows us to predict the most potential drug targets among large number of the non-homologous proteins involved in the unique metabolic pathway. A computational comparative metabolic pathway analysis of the host H. sapiens and the pathogen C pneumoniae AR39 has been carried out at three level analyses. Firstly, metabolic pathway analysis was performed to identify unique metabolic pathways and non-homologous proteins were identified. Secondly, essentiality of the proteins was checked, where these proteins contribute to the growth and survival of the organism. Finally these proteins were further subjected to predict protein interaction networks. Among the total 65 pathways in the C pneumoniae AR39 genome 10 were identified as the unique metabolic pathways which were not found in the human host, 32 enzymes were predicted as essential and these proteins were considered for protein interaction analysis, later using various criteria''s we have narrowed down to prioritize ribonucleotide-diphosphate reductase subunit beta as a potential drug target which facilitate for the successful entry into drug designing.