From Metabolomics to Fluxomics: A Computational Procedure to Translate Metabolite Profiles into Metabolic Fluxes

We describe a believed-novel procedure for translating metabolite profiles (metabolome) into the set of metabolic fluxes (fluxome) from which they originated. Methodologically, computational modeling is integrated with an analytical platform comprising linear optimization, continuation and dynamic analyses, and metabolic control. The procedure was tested with metabolite profiles obtained from ex vivo mice Langendorff-heart preparations perfused with glucose. The metabolic profiles were analyzed using a detailed kinetic model of the glucose catabolic pathways including glycolysis, pentose phosphate (PP), glycogenolysis, and polyols to translate the glucose metabolome of the heart into the fluxome. After optimization, the ability of the model to simulate the initial metabolite profile was confirmed, and metabolic fluxes as well as the structure of control and regulation of the glucose catabolic network could be calculated. We show that the step catalyzed by phosphofructokinase together with ATP demand and glycogenolysis exert the highest control on the glycolytic flux. The negative flux control exerted by phosphofructokinase on the PP and polyol pathways revealed that the extent of glycolytic flux directly affects flux redirection through these pathways, i.e., the higher the glycolytic flux the lower the PP and polyols. This believed-novel methodological approach represents a step forward that may help in designing therapeutic strategies targeted to diagnose, prevent, and treat metabolic diseases.     

Glycolysis Is Governed by Growth Regime and Simple Enzyme Regulation in Adherent MDCK Cells

Due to its vital importance in the supply of cellular pathways with energy and precursors, glycolysis has been studied for several decades regarding its capacity and regulation. For a systems-level understanding of the Madin-Darby canine kidney (MDCK) cell metabolism, we couple a segregated cell growth model published earlier with a structured model of glycolysis, which is based on relatively simple kinetics for enzymatic reactions of glycolysis, to explain the pathway dynamics under various cultivation conditions. The structured model takes into account in vitro enzyme activities, and links glycolysis with pentose phosphate pathway and glycogenesis. Using a single parameterization, metabolite pool dynamics during cell cultivation, glucose limitation and glucose pulse experiments can be consistently reproduced by considering the cultivation history of the cells. Growth phase-dependent glucose uptake together with cell-specific volume changes generate high intracellular metabolite pools and flux rates to satisfy the cellular demand during growth. Under glucose limitation, the coordinated control of glycolytic enzymes re-adjusts the glycolytic flux to prevent the depletion of glycolytic intermediates. Finally, the model''s predictive power supports the design of more efficient bioprocesses.     

Microtubular protein in its polymerized or nonpolymerized states differentially modulates in vitro and intracellular fluxes catalyzed by enzymes of carbon metabolism

The fluxes through HK/G6PDH and PK/LDH coupled-enzymatic reactions were quantified in the presence of physiological concentrations (1–15 μM) of polymerized or non-polymerized microtubular protein (MTP) from rat brain and in a permeabilized yeast cell system. In vitro enzymatic fluxes were increased by either polymerized or nonpolymerized brain MTP mainly in the lower range of MTP concentration. At fixed MTP concentrations in the flux stimulatory range of HK/G6PDH (1 mg/ml MTP) or PK/LDH (0.4 mg/ml MTP), a hyperbolic and sigmoidal response to NADP and PEP, respectively, was detected. That dependence varied according to the polymeric status of MTP. The specificity of the phenomenon observed in vitro, was tested for the PK/LDH and HK/G6PDH enzymatic couples in the presence of neutral polymers such as glycogen (? 10 mg/ml), poly(ethylene glycol) (up to 10% w/w) or G-actin (? 1 mg/ml). In permeabilized Saccharomyces cerevisiae cells, the PK-catalyzed flux was sensitive to microtubule disruption by nocodazole (15 μg/ml). The HK/G6PDH system was not affected by nocodazole showing values of kinetic parameters close to those obtained in vitro in the presence of polymerized brain MTP. Indirect immunofluorescence with specific antibodies against tubulin allowed to confirm the microtubules disruption in the presence of nocodazole in permeabilized yeast cells under the same conditions in which enzymes were assayed intracellularly. The experimental evidence is in agreement with the observed phenomenon of increase in fluxes in the enzymatic reactions assayed to be specifically induced by MTP either in vitro or in situ. The results presented are discussed in terms of the assembly of large supramolecular structures as a supraregulatory mechanism of synchronization of systemic cellular processes such as metabolic fluxes. © 1994 Wiley-Liss, Inc.     

Retinoblastoma protein promotes oxidative phosphorylation through upregulation of glycolytic genes in oncogene‐induced senescent cells

Metabolism is closely linked with cellular state and biological processes, but the mechanisms controlling metabolic properties in different contexts remain unclear. Cellular senescence is an irreversible growth arrest induced by various stresses, which exhibits active secretory and metabolic phenotypes. Here, we show that retinoblastoma protein (RB) plays a critical role in promoting the metabolic flow by activating both glycolysis and mitochondrial oxidative phosphorylation (OXPHOS) in cells that have undergone oncogene‐induced senescence (OIS). A combination of real‐time metabolic monitoring, and metabolome and gene expression analyses showed that OIS‐induced fibroblasts developed an accelerated metabolic flow. The loss of RB downregulated a series of glycolytic genes and simultaneously reduced metabolites produced from the glycolytic pathway, indicating that RB upregulates glycolytic genes in OIS cells. Importantly, both mitochondrial OXPHOS and glycolytic activities were abolished in RB‐depleted or downstream glycolytic enzyme‐depleted OIS cells, suggesting that RB‐mediated glycolytic activation induces a metabolic flux into the OXPHOS pathway. Collectively, our findings reveal that RB essentially functions in metabolic remodeling and the maintenance of the active energy production in OIS cells.     

Improving lycopene production in Escherichia coli by engineering metabolic control

Metabolic engineering has achieved encouraging success in producing foreign metabolites in a variety of hosts. However, common strategies for engineering metabolic pathways focus on amplifying the desired enzymes and deregulating cellular controls. As a result, uncontrolled or deregulated metabolic pathways lead to metabolic imbalance and suboptimal productivity. Here we have demonstrated the second stage of metabolic engineering effort by designing and engineering a regulatory circuit to control gene expression in response to intracellular metabolic states. Specifically, we recruited and altered one of the global regulatory systems in Escherichia coli, the Ntr regulon, to control the engineered lycopene biosynthesis pathway. The artificially engineered regulon, stimulated by excess glycolytic flux through sensing of an intracellular metabolite, acetyl phosphate, controls the expression of two key enzymes in lycopene synthesis in response to flux dynamics. This intracellular control loop significantly enhanced lycopene production while reducing the negative impact caused by metabolic imbalance. Although we demonstrated this strategy for metabolite production, it can be extended into other fields where gene expression must be closely controlled by intracellular physiology, such as gene therapy.     

Overexpression of a cytosolic pyrophosphatase (TgPPase) reveals a regulatory role of PP(i) in glycolysis for Toxoplasma gondii

PP(i) is a critical element of cellular metabolism as both an energy donor and as an allosteric regulator of several metabolic pathways. The apicomplexan parasite Toxoplasma gondii uses PP(i) in place of ATP as an energy donor in at least two reactions: the glycolytic PP(i)-dependent PFK (phosphofructokinase) and V-H(+)-PPase [vacuolar H(+)-translocating PPase (pyrophosphatase)]. In the present study, we report the cloning, expression and characterization of cytosolic TgPPase (T. gondii soluble PPase). Amino acid sequence alignment and phylogenetic analysis indicates that the gene encodes a family I soluble PPase. Overexpression of the enzyme in extracellular tachyzoites led to a 6-fold decrease in the cytosolic concentration of PP(i) relative to wild-type strain RH tachyzoites. Unexpectedly, this subsequent reduction in PP(i) was associated with a higher glycolytic flux in the overexpressing mutants, as evidenced by higher rates of proton and lactate extrusion. In addition to elevated glycolytic flux, TgPPase-overexpressing tachyzoites also possessed higher ATP concentrations relative to wild-type RH parasites. These results implicate PP(i) as having a significant regulatory role in glycolysis and, potentially, other downstream processes that regulate growth and cell division.     

Experimental validation of metabolic pathway modeling

In the search for new drug targets in the human parasite Entamoeba histolytica, metabolic control analysis was applied to determine, experimentally, flux control distribution of amebal glycolysis. The first (hexokinase, hexose-6-phosphate isomerase, pyrophosphate-dependent phosphofructokinase (PP(i)-PFK), aldolase and triose-phosphate isomerase) and final (3-phosphoglycerate mutase, enolase and pyruvate phosphate dikinase) glycolytic segments were reconstituted in vitro with recombinant enzymes under near-physiological conditions of pH, temperature and enzyme proportion. Flux control was determined by titrating flux with each enzyme component. In parallel, both glycolytic segments were also modeled by using the rate equations and kinetic parameters previously determined. Because the flux control distribution predicted by modeling and that determined by reconstitution were not similar, kinetic interactions among all the reconstituted components were experimentally revised to unravel the causes of the discrepancy. For the final segment, it was found that 3-phosphoglycerate was a weakly competitive inhibitor of enolase, whereas PP(i) was a moderate inhibitor of 3-phosphoglycerate mutase and enolase. For the first segment, PP(i) was both a strong inhibitor of aldolase and a nonessential mixed-type activator of amebal hexokinase; in addition, lower V(max) values for hexose-6-phosphate isomerase, PP(i)-PFK and aldolase were induced by PP(i) or ATP inhibition. It should be noted that PP(i) and other metabolites were absent from the 3-phosphoglycerate mutase and enolase or aldolase and hexokinase kinetics experiments, but present in reconstitution experiments. Only by incorporating these modifications in the rate equations, modeling predicted values of flux control distribution, flux rate and metabolite concentrations similar to those experimentally determined. The experimentally validated segment models allowed 'in silico experimentation' to be carried out, which is not easy to achieve in in vivo or in vitro systems. The results predicted a nonsignificant effect on flux rate and flux control distribution by adding parallel routes (pyruvate kinase for the final segment and ATP-dependent PFK for the first segment), because of the much lower activity of these enzymes in the ameba. Furthermore, modeling predicted full flux-control by 3-phosphoglycerate mutase and hexokinase, in the presence of low physiological substrate and product concentrations. It is concluded that the combination of in vitro pathway reconstitution with modeling and enzyme kinetics experimentation permits a more comprehensive understanding of the pathway behavior and control properties.     

Flux balance analysis of Chlorella sp. FC2 IITG under photoautotrophic and heterotrophic growth conditions

Quantification of carbon flux distribution in the metabolic network of microalgae remains important to understand the complex interplay between energy metabolism, carbon fixation, and assimilation pathways. This is even more relevant with respect to cyclic metabolism of microalgae under light–dark cycle. In the present study, flux balance analysis (FBA) was carried out for an indigenous isolate Chlorella sp. FC2 IITG under photoautotrophic and heterotrophic growth conditions. A shift in intracellular flux distribution was predicted during transition from nutrient sufficient phase to nutrient starvation phase of growth. Further, dynamic flux analysis (dFBA) was carried out to capture light–dark metabolism over discretized pseudo steady state time intervals. Our key findings include the following: (i) unlike heterotrophic condition, oxidative pentose phosphate (PP) pathway, and Krebs cycle were relatively inactive under photoautotrophic growth; (ii) in both growth conditions, while transhydrogenation reaction was highly active, glyoxalate shunt was found to be nonoperative; (iii) flux distribution during transition period was marked with up regulation of carbon flux toward nongrowth associated (NGA) maintenance energy, oxidative phosphorylation, and photophosphorylation; (iv) redirection of carbon flux from polysaccharide and neutral lipid resulted in up regulation of Krebs cycle flux in the dark phase; (v) elevated glycolytic and acetyl-CoA flux were coupled with induction of neutral lipid during light cycle of the growth; (vi) significantly active photophosphorylation in the light phase was able to satisfy cellular energy requirement without need of oxidative PP pathway; and (vi) unlike static FBA, dFBA predicted an unaltered NGA maintenance energy of 1.5 mmol g^?1 DCW h^?1.     

Metabolic control analysis of the Warburg-effect in proliferating vascular smooth muscle cells

Summary The accumulation and proliferation of vascular smooth muscle cells (VSMC) within the vessel wall is an important pathogenic feature in the development of atherosclerosis. Glucose metabolism has been implicated to play an important role in this cellular mechanism. To further elucidate the role of glucose metabolism in atherogenesis, glycolysis and its regulation have been investigated in proliferating VSMC. Platelet derived growth factor (PDGF BB)-induced proliferation of VSMCs significantly stimulated glucose flux through glycolysis. Further evaluating the enzymatic regulation of this pathway, the analysis of flux:metabolite co-responses revealed that anaerobic glycolytic flux is controlled at different sites of gycolysis in proliferating VSMCs, being consistent with the concept of multisite modulation. These findings indicate that regulation of glycolytic flux in proliferating VSMCs differs from traditional concepts of metabolic control of the Embden–Meyerhof pathway.     

Kinetic modeling can describe in vivo glycolysis in Entamoeba histolytica

Glycolysis in the human parasite Entamoeba histolytica is characterized by the absence of cooperative modulation and the prevalence of pyrophosphate-dependent (over ATP-dependent) enzymes. To determine the flux-control distribution of glycolysis and understand its underlying control mechanisms, a kinetic model of the pathway was constructed by using the software gepasi. The model was based on the kinetic parameters determined in the purified recombinant enzymes, and the enzyme activities, and steady-state fluxes and metabolite concentrations determined in amoebal trophozoites. The model predicted, with a high degree of accuracy, the flux and metabolite concentrations found in trophozoites, but only when the pyrophosphate concentration was held constant; at variable pyrophosphate, the model was not able to completely account for the ATP production/consumption balance, indicating the importance of the pyrophosphate homeostasis for amoebal glycolysis. Control analysis by the model revealed that hexokinase exerted the highest flux control (73%), as a result of its low cellular activity and strong AMP inhibition. 3-Phosphoglycerate mutase also exhibited significant flux control (65%) whereas the other pathway enzymes showed little or no control. The control of the ATP concentration was also mainly exerted by ATP consuming processes and 3-phosphoglycerate mutase and hexokinase (in the producing block). The model also indicated that, in order to diminish the amoebal glycolytic flux by 50%, it was required to decrease hexokinase or 3-phosphoglycerate mutase by 24% and 55%, respectively, or by 18% for both enzymes. By contrast, to attain the same reduction in flux by inhibiting the pyrophosphate-dependent enzymes pyrophosphate-phosphofructokinase and pyruvate phosphate dikinase, they should be decreased > 70%. On the basis of metabolic control analysis, steps whose inhibition would have stronger negative effects on the energy metabolism of this parasite were identified, thus becoming alternative targets for drug design.     

Exploring the Altered Dynamics of Mammalian Central Carbon Metabolic Pathway in Cancer Cells: A Classical Control Theoretic Approach

Background

In contrast with normal cells, most of the cancer cells depend on aerobic glycolysis for energy production in the form of adenosine triphosphate (ATP) bypassing mitochondrial oxidative phosphorylation. Moreover, compared to normal cells, cancer cells exhibit higher consumption of glucose with higher production of lactate. Again, higher rate of glycolysis provides the necessary glycolytic intermediary precursors for DNA, protein and lipid synthesis to maintain high active proliferation of the tumor cells. In this scenario, classical control theory based approach may be useful to explore the altered dynamics of the cancer cells. Since the dynamics of the cancer cells is different from that of the normal cells, understanding their dynamics may lead to development of novel therapeutic strategies.

Method

We have developed a model based on the state space equations of classical control theory along with an order reduction technique to mimic the actual dynamic behavior of mammalian central carbon metabolic (CCM) pathway in normal cells. Here, we have modified Michaelis Menten kinetic equation to incorporate feedback mechanism along with perturbations and cross talks associated with a metabolic pathway. Furthermore, we have perturbed the proposed model to reduce the mitochondrial oxidative phosphorylation. Thereafter, we have connected proportional-integral (PI) controller(s) with the model for tuning it to behave like the CCM pathway of a cancer cell. This methodology allows one to track the altered dynamics mediated by different enzymes.

Results and Discussions

The proposed model successfully mimics all the probable dynamics of the CCM pathway in normal cells. Moreover, experimental results demonstrate that in cancer cells, a coordination among enzymes catalyzing pentose phosphate pathway and intermediate glycolytic enzymes along with switching of pyruvate kinase (M2 isoform) plays an important role to maintain their altered dynamics.     

Control analysis as a tool to understand the formation of the las operon in Lactococcus lactis

In Lactococcus lactis the enzymes phosphofructokinase (PFK), pyruvate kinase (PK) and lactate dehydrogenase (LDH) are uniquely encoded in the las operon. We used metabolic control analysis to study the role of this organization. Earlier studies have shown that, at wild-type levels, LDH has no control over glycolysis and growth rate, but high negative control over formate production (C(Jformate)LDH=-1.3). We found that PFK and PK exert no control over glycolysis and growth rate at wild-type enzyme levels but both enzymes exert strong positive control on the glycolytic flux at reduced activities. PK exerts high positive control over formate (C(Jformate)PK=0.9-1.1) and acetate production (C(Jacetate)PK=0.8-1.0), whereas PFK exerts no control over these fluxes at increased expression. Decreased expression of the entire las operon resulted in a strong decrease in the growth rate and glycolytic flux; at 53% expression of the las operon glycolytic flux was reduced to 44% and the flux control coefficient increased towards 3. Increased las expression resulted in a slight decrease in the glycolytic flux. At wild-type levels, control was close to zero on both glycolysis and the pyruvate branches. The sum of control coefficients for the three enzymes individually was comparable with the control coefficient found for the entire operon; the strong positive control exerted by PK almost cancels out the negative control exerted by LDH on formate production. Our analysis suggests that coregulation of PFK and PK provides a very efficient way to regulate glycolysis, and coregulating PK and LDH allows cells to maintain homolactic fermentation during glycolysis regulation.     

Mathematical modelling of metabolic pathways affected by an enzyme deficiency. A mathematical model of glycolysis in normal and pyruvate-kinase-deficient red blood cells

A mathematical model of glycolysis in human erythrocytes is proposed to study the influence of a pyruvate kinase deficiency on the energy metabolism. The model takes into account the main regulatory properties of the non-equilibrium enzymes and the magnesium-complex formation by the adenine nucleotides and by 2,3-bisphosphoglycerate. In the normal case (no enzyme defect) the calculated flux rates and metabolite concentrations are in a good agreement with experimental data. It is shown that a severe pyruvate kinase deficiency manifested in a tenfold diminished activity of that enzyme leads to a remarkable decrease of the glycolytic flux and the ATP concentration of about 50% of the normal values. On the other hand a lowering of the pyruvate kinase activity to half of the normal value, characteristic for the heterozygotes, gives no significant alterations of the metabolite concentrations and the flux rates compared with the normal case which is in accordance with the lack of clinical symptoms for a metabolic disease of these probands. For three patients with known alterations of their pyruvate kinase mutants the calculated metabolite concentrations and the control characteristics permit estimation of the degree of disorder of the glycolytic pathway. The resulting classification corresponds well to other independent experimental and clinical findings. In particular, the calculation demonstrates that there is no simple correlation between the lowered enzyme activity and the reduced flux rate through the affected pathway.     

Response of fluxome and metabolome to temperature-induced recombinant protein synthesis in Escherichia coli

The response of the central carbon metabolism of Escherichia coli to temperature-induced recombinant production of human fibroblast growth factor was studied on the level of metabolic fluxes and intracellular metabolite levels. During production, E. coli TG1:plambdaFGFB, carrying a plasmid encoded gene for the recombinant product, revealed stress related characteristics such as decreased growth rate and biomass yield and enhanced by-product excretion (acetate, pyruvate, lactate). With the onset of production, the adenylate energy charge dropped from 0.85 to 0.60, indicating the occurrence of a severe energy limitation. This triggered an increase of the glycolytic flux which, however, was not sufficient to compensate for the increased ATP demand. The activation of the glycolytic flux was also indicated by the readjustment of glycolytic pool sizes leading to an increased driving force for the reaction catalyzed by phosphofructokinase. Moreover, fluxes through the TCA cycle, into the pentose phosphate pathway and into anabolic pathways decreased significantly. The strong increase of flux into overflow pathways, especially towards acetate was most likely caused by a flux redirection from pyruvate dehydrogenase to pyruvate oxidase. The glyoxylate shunt, not active during growth, was the dominating anaplerotic pathway during production. Together with pyruvate oxidase and acetyl CoA synthase this pathway could function as a metabolic by-pass to overcome the limitation in the junction between glycolysis and TCA cycle and partly recycle the acetate formed back into the metabolism.     

The effect of orotic acid treatment on the energy and carbohydrate metabolism of the hypertrophying rat heart

• 1.1. Adenine nucleotide concentrations in normal and one day hypertrophied hearts of untreated, orotic acid (OA), uridine, uracil, dihydroorotate and reserpine pretreated rats were measured. OA treatment increased the ADP concentration 5-fold in one day hypertrophied hearts. Neither uracil, uridine, dihydroorotate nor reserpine treatments changed ADP or total adenylate concentrations at one day of hypertrophy.
• 2.2. The adenine nucleotide ratio (ANR) at 0.263 × 10³ M⁻¹ and the energy charge (0.66) were at their lowest values in OA and in reserpine treated one day hypertrophying hearts. The temporal decline in the indices of energy metabolism corresponded with the OA induced maximum stimulation of contractility and maximum rates of protein, RNA and glycogen synthesis.
• 3.3. The phosphorylation state of the adenine nucleotides (PSAN) was both the most sensitive and the best predictive index of the cellular energy status in normal and hypertrophying hearts. The pronounced ability of OA treatment to energize myocyte cytoplasm was shown by the 9- and 6-fold greater values of PSAN over ANR in one and three day hypertrophied hearts. The enhanced PSAN may be the key factor in the mechanism of OA induced enhancement of contractile and synthetic functions of the heart in compensatory hypertrophy.
• 4.4. The development of myocardial hypertrophy in untreated rats resulted in a 36% reduction in the cytoplasmic NAD/NADH ratio. In rats treated with OA this redox couple of the hypertrophying heart was more oxidized and was increased by 30% to restore it to the value range of normal heart.
• 5.5. The regulatory status of the glycolytic pathway in untreated and OA treated hypertrophying hearts was assessed by comparisons of the mass action ratio (MAR) and equilibrium constants for each of the individual glycolytic reactions. There was an OA induced 2.7-fold increase in glycogen, UDP-glucose and total uridine nucleotides in hypertrophied hearts. The concentrations of seven out of ten glycolytic intermediates, including pyruvate and lactate were increased as a consequence of OA treated hypertrophy. Glycolytic flux was not stalled, rather the pathway was “more open” permitting greater throughput of intermediates with individually increased levels of selected metabolites. OA stimulated hypertrophy did not change the canonical control of glycolysis by the activities and individual MAR values of phosphofructokinase and pyruvic kinase.
• 6.6. Elevated levels of Glu 6-P, Fru 6-P and DHAP can force glycolytic intermediate entry into the non-oxidative reaction segment of the pentose pathway (PP), thereby elevating Rib 5-P concentration by reversal of the conventional flux direction of PP. Rib 5-P is rate limiting for PRPP and nucleotide synthesis and increases in its concentration in OA treated hypertrophying hearts can inter alia explain the elevation in adenylate concentrations.
• 7.7. OA does not act directly on the isolated normal or hypertrophying heart neither as an inotropic agent nor as a metabolic substrate. Its myocardial action requires whole-body integration and its principal metabolic fate involves the liver and the activation of the salvage pathway of pyrimidine biosynthesis. Preformed bases and nucleosides, mostly as uridine, enter and enhance the domains of pyrimidine, purine nucleotide and RNA metabolism in the OA stimulated hypertrophying myocardium.

     

Robustness analysis of the Escherichia coli metabolic network

Genomic, biochemical, and strain-specific data can be assembled to define an in silico representation of the metabolic network for a select group of single cellular organisms. Flux-balance analysis and phenotypic phase planes derived therefrom have been developed and applied to analyze the metabolic capabilities and characteristics of Escherichia coli K-12. These analyses have shown the existence of seven essential reactions in the central metabolic pathways (glycolysis, pentose phosphate pathway, tricarboxylic acid cycle) for the growth in glucose minimal media. The corresponding seven gene products can be grouped into three categories: (1) pentose phosphate pathway genes, (2) three-carbon glycolytic genes, and (3) tricarboxylic acid cycle genes. Here we develop a procedure that calculates the sensitivity of optimal cellular growth to altered flux levels of these essential gene products. The results indicate that the E. coli metabolic network is robust with respect to the flux levels of these enzymes. The metabolic flux in the transketolase and the tricarboxylic acid cycle reactions can be reduced to 15% and 19%, respectively, of the optimal value without significantly influencing the optimal growth flux. The metabolic network also exhibited robustness with respect to the ribose-5-phosphate isomerase, and the ribose-5-phosephate isomerase flux was reduced to 28% of the optimal value without significantly effecting the optimal growth flux. The metabolic network exhibited limited robustness to the three-carbon glycolytic fluxes both increased and decreased. The development presented another dimension to the use of FBA to study the capabilities of metabolic networks.     

In vivo interpretation of model predicted inhibition in acrylate pathway engineered Lactococcus lactis

To maximize the productivity of engineered metabolic pathway, in silico model is an established means to provide features of enzyme reaction dynamics. In our previous study, Escherichia coli engineered with acrylate pathway yielded low propionic acid titer. To understand the bottleneck behind this low productivity, a kinetic model was developed that incorporates the enzymatic reactions of the acrylate pathway. The resulting model was capable of simulating the fluxes reported under in vitro studies with good agreement, suggesting repression of propionyl-CoA transferase (Pct) by carboxylate metabolites as the main limiting factor for propionate production. Furthermore, the predicted flux control coefficients of the pathway enzymes under steady state conditions revealed that the control of flux is shared between Pct and lactoyl-CoA dehydratase. Increase in lactate concentration showed gradual decrease in flux control coefficients of Pct that in turn confirmed the control exerted by the carboxylate substrate. To interpret these in silico predictions under in vivo system, an organized study was conducted with a lactic acid bacteria strain engineered with acrylate pathway. Analysis reported a decreased product formation rate on attainment of inhibitory titer by suspected metabolites and supported the model.