Dissecting the genetic and metabolic mechanisms of adaptation to the knockout of a major metabolic enzyme in Escherichia coli

Unraveling the mechanisms of microbial adaptive evolution following genetic or environmental challenges is of fundamental interest in biological science and engineering. When the challenge is the loss of a metabolic enzyme, adaptive responses can also shed significant insight into metabolic robustness, regulation, and areas of kinetic limitation. In this study, whole-genome sequencing and highresolution C-13-metabolic flux analysis were performed on 10 adaptively evolved pgi knockouts of Escherichia coli. Pgi catalyzes the first reaction in glycolysis, and its loss results in major physiological and carbon catabolism pathway changes, including an 80% reduction in growth rate. Following adaptive laboratory evolution (ALE), the knockouts increase their growth rate by up to 3.6-fold. Through combined genomic-fluxomic analysis, we characterized the mutations and resulting metabolic fluxes that enabled this fitness recovery. Large increases in pyridine cofactor transhydrogenase flux, correcting imbalanced production of NADPH and NADH, were enabled by direct mutations to the transhydrogenase genes sthA and pntAB. The phosphotransferase system component crr was also found to be frequently mutated, which corresponded to elevated flux from pyruvate to phosphoenolpyruvate. The overall energy metabolism was found to be strikingly robust, and what have been previously described as latently activated Entner-Doudoroff and glyoxylate shunt pathways are shown here to represent no real increases in absolute flux relative to the wild type. These results indicate that the dominant mechanism of adaptation was to relieve the rate-limiting steps in cofactor metabolism and substrate uptake and to modulate global transcriptional regulation from stress response to catabolism

Evolution of <i>E. coli</i> on [U^-13C] Glucose Reveals a Negligible Isotopic Influence on Metabolism and Physiology

13C-Metabolic flux analysis (13C-MFA) traditionally assumes that kinetic isotope effects from isotopically labeled compounds do not appreciably alter cellular growth or metabolism, despite indications that some biochemical reactions can be non-negligibly impacted. Here, populations of Escherichia coli were adaptively evolved for ~1000 generations on uniformly labeled 13C-glucose, a commonly used isotope for 13C-MFA. Phenotypic characterization of these evolved strains revealed ~40% increases in growth rate, with no significant difference in fitness when grown on either labeled (13C) or unlabeled (12C) glucose. The evolved strains displayed decreased biomass yields, increased glucose and oxygen uptake, and increased acetate production, mimicking what is observed after adaptive evolution on unlabeled glucose. Furthermore, full genome re-sequencing revealed that the key genetic changes underlying these phenotypic alterations were essentially the same as those acquired during adaptive evolution on unlabeled glucose. Additionally, glucose competition experiments demonstrated that the wild-type exhibits no isotopic preference for unlabeled glucose, and the evolved strains have no preference for labeled glucose. Overall, the results of this study indicate that there are no significant differences between 12C and 13C-glucose as a carbon source for E. coli growth

Fast growth phenotype of E. coli K-12 from adaptive laboratory evolution does not require intracellular flux rewiring

Adaptive laboratory evolution (ALE) is a widely-used method for improving the fitness of microorganisms in selected environmental conditions. It has been applied previously to Escherichia coli K-12 MG1655 during aerobic exponential growth on glucose minimal media, a frequently used model organism and growth condition, to probe the limits of E. coli growth rate and gain insights into fast growth phenotypes. Previous studies have described up to 1.6-fold increases in growth rate following ALE, and have identified key causal genetic mutations and changes in transcriptional patterns. Here, we report for the first time intracellular metabolic fluxes for six such adaptively evolved strains, as determined by high-resolution 13C-metabolic flux analysis. Interestingly, we found that intracellular metabolic pathway usage changed very little following adaptive evolution. Instead, at the level of central carbon metabolism the faster growth was facilitated by proportional increases in glucose uptake and all intracellular rates. Of the six evolved strains studied here, only one strain showed a small degree of flux rewiring, and this was also the strain with unique genetic mutations. A comparison of fluxes with two other wild-type (unevolved) E. coli strains, BW25113 and BL21, showed that inter-strain differences are greater than differences between the parental and evolved strains. Principal component analysis highlighted that nearly all flux differences (95%) between the nine strains were captured by only two principal components. The distance between measured and flux balance analysis predicted fluxes was also investigated. It suggested a relatively wide range of similar stoichiometric optima, which opens new questions about the path-dependency of adaptive evolution

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

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

13C-assisted metabolic flux analysis to investigate heterotrophic and mixotrophic metabolism in Cupriavidus necator H16

Introduction. Cupriavidus necator H16 is a gram-negative bacterium, capable of lithoautotrophic growth by utilizing hydrogen as an energy source and fixing carbon dioxide (CO2) through Calvin-Benson-Bassham (CBB) cycle. The potential to utilize synthesis gas (Syngas) and the prospects of rerouting carbon from polyhydroxybutyrate synthesis to value-added compounds makes C. necator an excellent chassis for industrial application. Objectives. In the context of lack of sufficient quantitative information of the metabolic pathways and to advance in rational metabolic engineering for optimized product synthesis in C. necator H16, we carried out a metabolic flux analysis based on steady-state 13C-labelling. Methods. In this study, steady-state carbon labelling experiments, using either D-[1-13C]fructose or [1,2-13C]glycerol, were undertaken to investigate the carbon flux through the central carbon metabolism in C. necator H16 under heterotrophic and mixotrophic growth conditions, respectively. Results. We found that the CBB cycle is active even under heterotrophic condition, and growth is indeed mixotrophic. While Entner-Doudoroff (ED) pathway is shown to be the major route for sugar degradation, tricarboxylic acid (TCA) cycle is highly active in mixotrophic condition. Enhanced flux is observed in reductive pentose phosphate pathway (redPPP) under the mixotrophic condition to supplement the precursor requirement for CBB cycle. The flux distribution was compared to the mRNA abundance of genes encoding enzymes involved in key enzymatic reactions of the central carbon metabolism. Conclusion. This study leads the way to establishing 13C-based quantitative fluxomics for rational pathway engineering in C. necator H16

Oxygen dependence of metabolic fluxes and energy generation of Saccharomyces cerevisiae CEN.PK113-1A

<p>Abstract</p> <p>Background</p> <p>The yeast <it>Saccharomyces cerevisiae </it>is able to adjust to external oxygen availability by utilizing both respirative and fermentative metabolic modes. Adjusting the metabolic mode involves alteration of the intracellular metabolic fluxes that are determined by the cell's multilevel regulatory network. Oxygen is a major determinant of the physiology of <it>S. cerevisiae </it>but understanding of the oxygen dependence of intracellular flux distributions is still scarce.</p> <p>Results</p> <p>Metabolic flux distributions of <it>S. cerevisiae </it>CEN.PK113-1A growing in glucose-limited chemostat cultures at a dilution rate of 0.1 h^-1with 20.9%, 2.8%, 1.0%, 0.5% or 0.0% O₂in the inlet gas were quantified by ¹³C-MFA. Metabolic flux ratios from fractional [U-¹³C]glucose labelling experiments were used to solve the underdetermined MFA system of central carbon metabolism of <it>S. cerevisiae</it>.</p> <p>While ethanol production was observed already in 2.8% oxygen, only minor differences in the flux distribution were observed, compared to fully aerobic conditions. However, in 1.0% and 0.5% oxygen the respiratory rate was severely restricted, resulting in progressively reduced fluxes through the TCA cycle and the direction of major fluxes to the fermentative pathway. A redistribution of fluxes was observed in all branching points of central carbon metabolism. Yet only when oxygen provision was reduced to 0.5%, was the biomass yield exceeded by the yields of ethanol and CO₂. Respirative ATP generation provided 59% of the ATP demand in fully aerobic conditions and still a substantial 25% in 0.5% oxygenation. An extensive redistribution of fluxes was observed in anaerobic conditions compared to all the aerobic conditions. Positive correlation between the transcriptional levels of metabolic enzymes and the corresponding fluxes in the different oxygenation conditions was found only in the respirative pathway.</p> <p>Conclusion</p> <p>¹³C-constrained MFA enabled quantitative determination of intracellular fluxes in conditions of different redox challenges without including redox cofactors in metabolite mass balances. A redistribution of fluxes was observed not only for respirative, respiro-fermentative and fermentative metabolisms, but also for cells grown with 2.8%, 1.0% and 0.5% oxygen. Although the cellular metabolism was respiro-fermentative in each of these low oxygen conditions, the actual amount of oxygen available resulted in different contributions through respirative and fermentative pathways.</p

13C labeling experiments at metabolic nonstationary conditions: An exploratory study

<p>Abstract</p> <p>Background</p> <p>Stimulus Response Experiments to unravel the regulatory properties of metabolic networks are becoming more and more popular. However, their ability to determine enzyme kinetic parameters has proven to be limited with the presently available data. In metabolic flux analysis, the use of ¹³C labeled substrates together with isotopomer modeling solved the problem of underdetermined networks and increased the accuracy of flux estimations significantly.</p> <p>Results</p> <p>In this contribution, the idea of increasing the information content of the dynamic experiment by adding ¹³C labeling is analyzed. For this purpose a small example network is studied by simulation and statistical methods. Different scenarios regarding available measurements are analyzed and compared to a non-labeled reference experiment. Sensitivity analysis revealed a specific influence of the kinetic parameters on the labeling measurements. Statistical methods based on parameter sensitivities and different measurement models are applied to assess the information gain of the labeled stimulus response experiment.</p> <p>Conclusion</p> <p>It was found that the use of a (specifically) labeled substrate will significantly increase the parameter estimation accuracy. An overall information gain of about a factor of six is observed for the example network. The information gain is achieved from the specific influence of the kinetic parameters towards the labeling measurements. This also leads to a significant decrease in correlation of the kinetic parameters compared to an experiment without ¹³C-labeled substrate.</p

Carbon Metabolism of Enterobacterial Human Pathogens Growing in Epithelial Colorectal Adenocarcinoma (Caco-2) Cells

Analysis of the genome sequences of the major human bacterial pathogens has provided a large amount of information concerning their metabolic potential. However, our knowledge of the actual metabolic pathways and metabolite fluxes occurring in these pathogens under infection conditions is still limited. In this study, we analysed the intracellular carbon metabolism of enteroinvasive Escherichia coli (EIEC HN280 and EIEC 4608-58) and Salmonella enterica Serovar Typhimurium (Stm 14028) replicating in epithelial colorectal adenocarcinoma cells (Caco-2). To this aim, we supplied [U-13C6]glucose to Caco-2 cells infected with the bacterial strains or mutants thereof impaired in the uptake of glucose, mannose and/or glucose 6-phosphate. The 13C-isotopologue patterns of protein-derived amino acids from the bacteria and the host cells were then determined by mass spectrometry. The data showed that EIEC HN280 growing in the cytosol of the host cells, as well as Stm 14028 replicating in the Salmonella-containing vacuole (SCV) utilised glucose, but not glucose 6-phosphate, other phosphorylated carbohydrates, gluconate or fatty acids as major carbon substrates. EIEC 4608-58 used C3-compound(s) in addition to glucose as carbon source. The labelling patterns reflected strain-dependent carbon flux via glycolysis and/or the Entner-Doudoroff pathway, the pentose phosphate pathway, the TCA cycle and anapleurotic reactions between PEP and oxaloacetate. Mutants of all three strains impaired in the uptake of glucose switched to C3-substrate(s) accompanied by an increased uptake of amino acids (and possibly also other anabolic monomers) from the host cell. Surprisingly, the metabolism of the host cells, as judged by the efficiency of 13C-incorporation into host cell amino acids, was not significantly affected by the infection with either of these intracellular pathogens

Influence of band width on the scattered ion yield spectra of a He + Ion by resonant or quasi-resonant charge exchange neutralization

The influence of the band structure, especially the bandwidth, on the scattered ion yield spectra of a He+ ion by the resonant or quasi-resonant neutralization was theoretically examined using quantum rate equations. When calculating the scattered ion yield spectra of He+ to simulate the experimental data, we observed that the band structure, especially the bandwidth, had a strong influence on the spectra at relatively low incident He+ ion energies of less than several hundred eV. Through many simulations, it was determined that theoretical calculations that include bandwidth calculation can simulate or reproduce the experimentally observed spectra of He+-In, He+-Ga, and He+-Sn systems. In contrast, simulations not including bandwidth simulation could neither reproduce nor account for such spectra. Furthermore, the calculated ion survival probability (ISP) at low incident ion energies tended to decrease with increasing bandwidth. This decrease in ISP probably corresponds to the relatively small scattered ion yield usually observed at low incident ion energies. Theoretically, such a decrease indicates that a He+ ion with a low incident energy can be easily neutralized on the surface when the bandwidth is large