Systems biology approach reveals that overflow metabolism of acetate in Escherichia coli is triggered by carbon catabolite repression of acetyl-CoA synthetase

Acetate: friend or foe? Efficient production of a sweet protein in Escherichia coli BL21 using acetate as a carbon source

BACKGROUND

Escherichia coli is, to date, the most used microorganism for the production of recombinant proteins and biotechnologically relevant metabolites. High density cell cultures allow efficient biomass and protein yields. However, their main limitation is the accumulation of acetate as a by-product of unbalanced carbon metabolism. Increased concentrations of acetate can inhibit cellular growth and recombinant protein production, and many efforts have been made to overcome this problem. On the other hand, it is known that E. coli is able to grow on acetate as the sole carbon source, although this mechanism has never been employed for the production of recombinant proteins.

RESULTS

By optimization of the fermentation parameters, we have been able to develop a new acetate containing medium for the production of a recombinant protein in E. coli BL21(DE3). The medium is based on a buffering phosphate system supplemented with 0.5% yeast extract for essential nutrients and sodium acetate as additional carbon source, and it is compatible with lactose induction. We tested these culture conditions for the production of MNEI, a single chain derivative of the sweet plant protein monellin, with potential for food and beverage industries. We noticed that careful oxygenation and pH control were needed for efficient protein production. The expression method was also coupled to a faster and more efficient purification technique, which allowed us to obtain MNEI with a purity higher than 99%.

CONCLUSIONS

The method introduced represents a new strategy for the production of MNEI in E. coli BL21(DE3) with a simple and convenient process, and offers a new perspective on the capabilities of this microorganism as a biotechnological tool. The conditions employed are potentially scalable to industrial processes and require only low-priced reagents, thus dramatically lowering production costs on both laboratory and industrial scale. The yield of recombinant MNEI in these conditions was the highest to date from E. coli cultures, reaching on average ~180 mg/L of culture, versus typical LB/IPTG yields of about 30 mg/L.

Effect of Global Regulators RpoS and Cyclic-AMP/CRP on the Catabolome and Transcriptome of Escherichia coli K12 during Carbon- and Energy-Limited Growth

For heterotrophic microbes, limited availability of carbon and energy sources is one of the major nutritional factors restricting the rate of growth in most ecosystems. Physiological adaptation to this hunger state requires metabolic versatility which usually involves expression of a wide range of different catabolic pathways and of high-affinity carbon transporters; together, this allows for simultaneous utilization of mixtures of carbonaceous compounds at low concentrations. In Escherichia coli the stationary phase sigma factor RpoS and the signal molecule cAMP are the major players in the regulation of transcription under such conditions; however, their interaction is still not fully understood. Therefore, during growth of E. coli in carbon-limited chemostat culture at different dilution rates, the transcriptomes, expression of periplasmic proteins and catabolomes of strains lacking one of these global regulators, either rpoS or adenylate cyclase (cya), were compared to those of the wild-type strain. The inability to synthesize cAMP exerted a strong negative influence on the expression of alternative carbon source uptake and degradation systems. In contrast, absence of RpoS increased the transcription of genes belonging to high-affinity uptake systems and central metabolism, presumably due to reduced competition of σ^D with σ^S. Phenotypical analysis confirmed this observation: The ability to respire alternative carbon substrates and to express periplasmic high-affinity binding proteins was eliminated in cya and crp mutants, while these properties were not affected in the rpoS mutant. As expected, transcription of numerous stress defence genes was negatively affected by the rpoS knock-out mutation. Interestingly, several genes of the RpoS stress response regulon were also down-regulated in the cAMP-negative strain indicating a coordinated global regulation. The results demonstrate that cAMP is crucial for catabolic flexibility during slow, carbon-limited growth, whereas RpoS is primarily involved in the regulation of stress response systems necessary for the survival of this bacterium under hunger conditions.

Salmonella typhimurium and Escherichia coli dissimilarity: Closely related bacteria with distinct metabolic profiles

Live attenuated strains of Salmonella typhimurium have been extensively investigated as vaccines for a number of infectious diseases. However, there is still little information available concerning aspects of their metabolism. S. typhimurium and Escherichia coli show a high degree of similarity in terms of their genome contents and metabolic networks. However, this work presents experimental evidence showing that significant differences exist in their abilities to direct carbon fluxes to biomass and energy production. It is important to study the metabolism of Salmonella to elucidate the formation of acetate and other metabolites involved in optimizing the production of biomass, essential for the development of recombinant vaccines. The metabolism of Salmonella under aerobic conditions was assessed using continuous cultures performed at dilution rates ranging from 0.1 to 0.67 h(-1), with glucose as main substrate. Acetate assimilation and glucose metabolism under anaerobic conditions were also investigated using batch cultures. Chemostat cultivations showed deviation of carbon towards acetate formation, starting at dilution rates above 0.1 h(-1). This differed from previous findings for E. coli, where acetate accumulation was only detected at dilution rates exceeding 0.4 h(-1), and was due to the lower rate of acetate assimilation by S. typhimurium under aerobic conditions. Under anaerobic conditions, both microorganisms mainly produced ethanol, acetate, and formate. A genome-scale metabolic model, reconstructed for Salmonella based on an E. coli model, provided a poor description of the mixed fermentation pattern observed during Salmonella cultures, reinforcing the different patterns of carbon utilization exhibited by these closely related bacteria.

Protein costs do not explain evolution of metabolic strategies and regulation of ribosomal content: does protein investment explain an anaerobic bacterial Crabtree effect?

Protein investment costs are considered a major driver for the choice of alternative metabolic strategies. We tested this premise in Lactococcus lactis, a bacterium that exhibits a distinct, anaerobic version of the bacterial Crabtree/Warburg effect; with increasing growth rates it shifts from a high yield metabolic mode [mixed-acid fermentation; 3 adenosine triphosphate (ATP) per glucose] to a low yield metabolic mode (homolactic fermentation; 2 ATP per glucose). We studied growth rate-dependent relative transcription and protein ratios, enzyme activities, and fluxes of L. lactis in glucose-limited chemostats, providing a high-quality and comprehensive data set. A three- to fourfold higher growth rate rerouted metabolism from acetate to lactate as the main fermentation product. However, we observed hardly any changes in transcription, protein levels and enzyme activities. Even levels of ribosomal proteins, constituting a major investment in cellular machinery, changed only slightly. Thus, contrary to the original hypothesis, central metabolism in this organism appears to be hardly regulated at the level of gene expression, but rather at the metabolic level. We conclude that L. lactis is either poorly adapted to growth at low and constant glucose concentrations, or that protein costs play a less important role in fitness than hitherto assumed.

Optimization of Fermentation Conditions for Recombinant Human Interferon Beta Production by Escherichia coli Using the Response Surface Methodology

Background:

The periplasmic overexpression of recombinant human interferon beta (rhIFN-β)-1b using a synthetic gene in Escherichia coli BL21 (DE3) was optimized in shake flasks using Response Surface Methodology (RSM) based on the Box-Behnken Design (BBD).

Objectives:

This study aimed to predict and develop the optimal fermentation conditions for periplasmic expression of rhIFN-β-1b in shake flasks whilst keeping the acetate excretion as the lowest amount and exploit the best results condition for rhIFN-β in a bench top bioreactor.

Materials and Methods:

The process variables studied were the concentration of glucose as carbon source, cell density prior the induction (OD _{600 nm}) and induction temperature. Ultimately, a three-factor three-level BBD was employed during the optimization process. The rhIFN-β production and the acetate excretion served as the evaluated responses.

Results:

The proposed optimum fermentation condition consisted of 7.81 g L^-1 glucose, OD _{600 nm} prior induction 1.66 and induction temperature of 30.27°C. The model prediction of 0.267 g L^-1 of rhIFN-β and 0.961 g L^-1 of acetate at the optimum conditions was verified experimentally as 0.255 g L^-1 and 0.981 g L^-1 of acetate. This agreement between the predicted and observed values confirmed the precision of the applied method to predict the optimum conditions.

Conclusions:

It can be concluded that the RSM is an effective method for the optimization of recombinant protein expression using synthetic genes in E. coli.

Integrating -Omics: Systems Biology as Explored Through C. elegans Research

-Omics data have become indispensable to systems biology, which aims to describe the full complexity of functional cells, tissues, organs and organisms. Generating vast amounts of data via such methods, researchers have invested in ways of handling and interpreting these. From the large volumes of -omics data that have been gathered over the years, it is clear that the information derived from one -ome is usually far from complete. Now, individual techniques and methods for integration are maturing to the point that researchers can focus on network-based integration rather than simply interpreting single -ome studies. This review evaluates the application of integrated -omics approaches with a focus on Caenorhabditis elegans studies, intending to direct researchers in this field to useful databases and inspiring examples.

Metabolic Regulation and Coordination of the Metabolism in Bacteria in Response to a Variety of Growth Conditions

Living organisms have sophisticated but well-organized regulation system. It is important to understand the metabolic regulation mechanisms in relation to growth environment for the efficient design of cell factories for biofuels and biochemicals production. Here, an overview is given for carbon catabolite regulation, nitrogen regulation, ion, sulfur, and phosphate regulations, stringent response under nutrient starvation as well as oxidative stress regulation, redox state regulation, acid-shock, heat- and cold-shock regulations, solvent stress regulation, osmoregulation, and biofilm formation, and quorum sensing focusing on Escherichia coli metabolism and others. The coordinated regulation mechanisms are of particular interest in getting insight into the principle which governs the cell metabolism. The metabolism is controlled by both enzyme-level regulation and transcriptional regulation via transcription factors such as cAMP-Crp, Cra, Csr, Fis, P(II)(GlnB), NtrBC, CysB, PhoR/B, SoxR/S, Fur, MarR, ArcA/B, Fnr, NarX/L, RpoS, and (p)ppGpp for stringent response, where the timescales for enzyme-level and gene-level regulations are different. Moreover, multiple regulations are coordinated by the intracellular metabolites, where fructose 1,6-bisphosphate (FBP), phosphoenolpyruvate (PEP), and acetyl-CoA (AcCoA) play important roles for enzyme-level regulation as well as transcriptional control, while α-ketoacids such as α-ketoglutaric acid (αKG), pyruvate (PYR), and oxaloacetate (OAA) play important roles for the coordinated regulation between carbon source uptake rate and other nutrient uptake rate such as nitrogen or sulfur uptake rate by modulation of cAMP via Cya.

Protein acetylation affects acetate metabolism, motility and acid stress response in Escherichia coli

Although protein acetylation is widely observed, it has been associated with few specific regulatory functions making it poorly understood. To interrogate its functionality, we analyzed the acetylome in Escherichia coli knockout mutants of cobB, the only known sirtuin-like deacetylase, and patZ, the best-known protein acetyltransferase. For four growth conditions, more than 2,000 unique acetylated peptides, belonging to 809 proteins, were identified and differentially quantified. Nearly 65% of these proteins are related to metabolism. The global activity of CobB contributes to the deacetylation of a large number of substrates and has a major impact on physiology. Apart from the regulation of acetyl-CoA synthetase, we found that CobB-controlled acetylation of isocitrate lyase contributes to the fine-tuning of the glyoxylate shunt. Acetylation of the transcription factor RcsB prevents DNA binding, activating flagella biosynthesis and motility, and increases acid stress susceptibility. Surprisingly, deletion of patZ increased acetylation in acetate cultures, which suggests that it regulates the levels of acetylating agents. The results presented offer new insights into functional roles of protein acetylation in metabolic fitness and global cell regulation.

Improving the expression of recombinant proteins in E. coli BL21 (DE3) under acetate stress: an alkaline pH shift approach

Excess acetate has long been an issue for the production of recombinant proteins in E. coli cells. Recently, improvements in acetate tolerance have been achieved through the use of genetic strategies and medium supplementation with certain amino acids and pyrimidines. The aim of our study was to evaluate an alternative to improve the acetate tolerance of E. coli BL21 (DE3), a popular strain used to express recombinant proteins. In this work we reported the cultivation of BL21 (DE3) in complex media containing acetate at high concentrations. In the presence of 300 mM acetate, compared with pH 6.5, pH 7.5 improved cell growth by approximately 71%, reduced intracellular acetate by approximately 50%, and restored the expression of glutathione S-transferase (GST), green fluorescent protein (GFP) and cytochrome P450 monooxygenase (CYP). Further experiments showed that alkaline pHs up to 8.5 had little inhibition in the expression of GST, GFP and CYP. In addition, the detrimental effect of acetate on the reduction of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) by the cell membrane, an index of cellular metabolic capacity, was substantially alleviated by a shift to alkaline pH values of 7.5–8.0. Thus, we suggest an approach of cultivating E. coli BL21 (DE3) at pH 8.0±0.5 to minimize the effects caused by acetate stress. The proposed strategy of an alkaline pH shift is a simple approach to solving similar bioprocessing problems in the production of biofuels and biochemicals from sugars.

TreeEFM: calculating elementary flux modes using linear optimization in a tree-based algorithm

MOTIVATION

Elementary flux modes (EFMs) analysis constitutes a fundamental tool in systems biology. However, the efficient calculation of EFMs in genome-scale metabolic networks (GSMNs) is still a challenge. We present a novel algorithm that uses a linear programming-based tree search and efficiently enumerates a subset of EFMs in GSMNs.

RESULTS

Our approach is compared with the EFMEvolver approach, demonstrating a significant improvement in computation time. We also validate the usefulness of our new approach by studying the acetate overflow metabolism in the Escherichia coli bacteria. To do so, we computed 1 million EFMs for each energetic amino acid and then analysed the relevance of each energetic amino acid based on gene/protein expression data and the obtained EFMs. We found good agreement between previous experiments and the conclusions reached using EFMs. Finally, we also analysed the performance of our approach when applied to large GSMNs.

AVAILABILITY AND IMPLEMENTATION

The stand-alone software TreeEFM is implemented in C++ and interacts with the open-source linear solver COIN-OR Linear program Solver (CLP).

Basic regulatory principles of Escherichia coli's electron transport chain for varying oxygen conditions

For adaptation between anaerobic, micro-aerobic and aerobic conditions Escherichia coli's metabolism and in particular its electron transport chain (ETC) is highly regulated. Although it is known that the global transcriptional regulators FNR and ArcA are involved in oxygen response it is unclear how they interplay in the regulation of ETC enzymes under micro-aerobic chemostat conditions. Also, there are diverse results which and how quinones (oxidised/reduced, ubiquinone/other quinones) are controlling the ArcBA two-component system. In the following a mathematical model of the E. coli ETC linked to basic modules for substrate uptake, fermentation product excretion and biomass formation is introduced. The kinetic modelling focusses on regulatory principles of the ETC for varying oxygen conditions in glucose-limited continuous cultures. The model is based on the balance of electron donation (glucose) and acceptance (oxygen or other acceptors). Also, it is able to account for different chemostat conditions due to changed substrate concentrations and dilution rates. The parameter identification process is divided into an estimation and a validation step based on previously published and new experimental data. The model shows that experimentally observed, qualitatively different behaviour of the ubiquinone redox state and the ArcA activity profile in the micro-aerobic range for different experimental conditions can emerge from a single network structure. The network structure features a strong feed-forward effect from the FNR regulatory system to the ArcBA regulatory system via a common control of the dehydrogenases of the ETC. The model supports the hypothesis that ubiquinone but not ubiquinol plays a key role in determining the activity of ArcBA in a glucose-limited chemostat at micro-aerobic conditions.

The pyruvate-tricarboxylic acid cycle node: a focal point of virulence control in the enteric pathogen Yersinia pseudotuberculosis

Despite our increasing knowledge of the specific pathogenicity factors in bacteria, the contribution of metabolic processes to virulence is largely unknown. Here, we elucidate a tight connection between pathogenicity and core metabolism in the enteric pathogen Yersinia pseudotuberculosis by integrated transcriptome and [(13)C]fluxome analysis of the wild type and virulence-regulator mutants. During aerobic growth on glucose, Y. pseudotuberculosis reveals an unusual flux distribution with a high level of secreted pyruvate. The absence of the transcriptional and post-transcriptional regulators RovA, CsrA, and Crp strongly perturbs the fluxes of carbon core metabolism at the level of pyruvate metabolism and the tricarboxylic acid (TCA) cycle, and these perturbations are accompanied by transcriptional changes in the corresponding enzymes. Knock-outs of regulators of this metabolic branch point and of its central enzyme, pyruvate kinase (ΔpykF), result in mutants with significantly reduced virulence in an oral mouse infection model. In summary, our work identifies the pyruvate-TCA cycle node as a focal point for controlling the host colonization and virulence of Yersinia.

An individual-based modeling approach to simulate the effects of cellular nutrient competition on Escherichia coli K-12 MG1655 colony behavior and interactions in aerobic structured food systems

Traditional kinetic models in predictive microbiology reliably predict macroscopic dynamics of planktonically-growing cell cultures in homogeneous liquid food systems. However, most food products have a semi-solid structure, where microorganisms grow locally in colonies. Individual colony cells exhibit strongly different and non-normally distributed behavior due to local nutrient competition. As a result, traditional models considering average population behavior in a homogeneous system do not describe colony dynamics in full detail. To incorporate local resource competition and individual cell differences, an individual-based modeling approach has been applied to Escherichia coli K-12 MG1655 colonies, considering the microbial cell as modeling unit. The first contribution of this individual-based model is to describe single colony growth under nutrient-deprived conditions. More specifically, the linear and stationary phase in the evolution of the colony radius, the evolution from a disk-like to branching morphology, and the emergence of a starvation zone in the colony center are simulated and compared to available experimental data. These phenomena occur earlier at more severe nutrient depletion conditions, i.e., at lower nutrient diffusivity and initial nutrient concentration in the medium. Furthermore, intercolony interactions have been simulated. Higher inoculum densities lead to stronger intercolony interactions, such as colony merging and smaller colony sizes, due to nutrient competition. This individual-based model contributes to the elucidation of characteristic experimentally observed colony behavior from mechanistic information about cellular physiology and interactions.

A mathematical model of metabolism and regulation provides a systems-level view of how Escherichia coli responds to oxygen

The efficient redesign of bacteria for biotechnological purposes, such as biofuel production, waste disposal or specific biocatalytic functions, requires a quantitative systems-level understanding of energy supply, carbon, and redox metabolism. The measurement of transcript levels, metabolite concentrations and metabolic fluxes per se gives an incomplete picture. An appreciation of the interdependencies between the different measurement values is essential for systems-level understanding. Mathematical modeling has the potential to provide a coherent and quantitative description of the interplay between gene expression, metabolite concentrations, and metabolic fluxes. Escherichia coli undergoes major adaptations in central metabolism when the availability of oxygen changes. Thus, an integrated description of the oxygen response provides a benchmark of our understanding of carbon, energy, and redox metabolism. We present the first comprehensive model of the central metabolism of E. coli that describes steady-state metabolism at different levels of oxygen availability. Variables of the model are metabolite concentrations, gene expression levels, transcription factor activities, metabolic fluxes, and biomass concentration. We analyze the model with respect to the production capabilities of central metabolism of E. coli. In particular, we predict how precursor and biomass concentration are affected by product formation.

Coordination of microbial metabolism

Beyond fuelling cellular activities with building blocks and energy, metabolism also integrates environmental conditions into intracellular signals. The underlying regulatory network is complex and multifaceted: it ranges from slow interactions, such as changing gene expression, to rapid ones, such as the modulation of protein activity via post-translational modification or the allosteric binding of small molecules. In this Review, we outline the coordination of common metabolic tasks, including nutrient uptake, central metabolism, the generation of energy, the supply of amino acids and protein synthesis. Increasingly, a set of key metabolites is recognized to control individual regulatory circuits, which carry out specific functions of information input and regulatory output. Such a modular view of microbial metabolism facilitates an intuitive understanding of the molecular mechanisms that underlie cellular decision making.

The metabolic potential of Escherichia coli BL21 in defined and rich medium

BACKGROUND

The proteome reflects the available cellular machinery to deal with nutrients and environmental challenges. The most common E. coli strain BL21 growing in different, commonly employed media was evaluated using a detailed quantitative proteome analysis.

RESULTS

The presence of preformed biomass precursor molecules in rich media such as Luria Bertani supported rapid growth concomitant to acetate formation and apparently unbalanced abundances of central metabolic pathway enzymes, e.g. high levels of lower glycolytic pathway enzymes as well as pyruvate dehydrogenase, and low levels of TCA cycle and high levels of the acetate forming enzymes Pta and AckA. The proteome of cells growing exponentially in glucose-supplemented mineral salt medium was dominated by enzymes of amino acid synthesis pathways, contained more balanced abundances of central metabolic pathway enzymes, and a lower portion of ribosomal and other translational proteins. Entry into stationary phase led to a reconstruction of the bacterial proteome by increasing e.g. the portion of proteins required for scavenging rare nutrients and general cell protection. This proteomic reconstruction during entry into stationary phase was more noticeable in cells growing in rich medium as they have a greater reservoir of recyclable proteins from the translational machinery.

CONCLUSIONS

The proteomic comparison of cells growing exponentially in different media reflected the antagonistic and competitive regulation of central metabolic pathways through the global transcriptional regulators Cra, Crp, and ArcA. For example, the proteome of cells growing exponentially in rich medium was consistent with a dominating role of phosphorylated ArcA most likely a result from limitations in reoxidizing reduced quinones in the respiratory chain under these growth conditions. The proteomic alterations of exponentially growing cells into stationary phase cells were consistent with stringent-like and stationary phase responses and a dominating control through DksA-ppGpp and RpoS.

Escherichia coli achieves faster growth by increasing catalytic and translation rates of proteins

Regulation levels of the gene expression cascade controlling protein levels and metabolic fluxes for cells to achieve faster growth have not been elaborated in acceptable detail. Furthermore, there is need for specific growth rate (μ) dependent absolute quantitative transcriptome and proteome data to understand the molecular relationships for enabling cells to modify μ. We address these questions, for the first time, by presenting regulatory strategies for more efficient metabolism of Escherichia coli at higher μ by statistical covariance analysis of genome-wide intracellular mRNA and protein concentrations coupled to metabolic flux analysis in the steady state range of μ = 0.11-0.49 h(-1). Our analyses show dominating post-transcriptional control of protein abundances and post-translational control of flux rates. On average, E. coli achieved five-times faster growth through 3.7-fold increase of apparent catalytic rates of enzymes (kapp) and 2.5-fold increased translation rates, demonstrating the relevance of post-translational regulation for increasing flux throughput. Interestingly, pathways carrying the highest flux showed both high protein abundance and kapp values. Furthermore, co-regulation analysis of enzymatic capacities revealed tightly coupled regulatory dependencies of protein synthesis and RNA precursor synthesis, substrate utilization, biosynthetic and energy generation pathways carrying the highest flux. We also observed metabolic pathway and COG specific protein and metabolic flux control levels, protein expression costs and genome-wide principles for translation efficiency and transcription unit polarity. This work contributes to the much needed quantitative understanding of coordinated gene expression regulation and metabolic flux control. Our findings will also advance modeling and metabolic engineering of industrial strains.

Fed-batch operation in special microtiter plates: a new method for screening under production conditions

Batch and fed-batch operation result in completely different physiological conditions for cultivated microorganisms or cells. To close the gap between screening, which is hitherto exclusively performed in batch mode, and fed-batch production processes, a special microtiter plate was developed that allows screening in fed-batch mode. The fed-batch microtiter plate (FB-MTP) enables 44 parallel fed-batch experiments at small scale. A small channel filled with a hydrogel connects a reservoir well with a culture well. The nutrient compound diffuses from the reservoir well through the hydrogel into the culture well. Hence, the feed rate can easily be adjusted to the needs of the cultured microorganisms by changing the geometry of the hydrogel channel and the driving concentration gradient. Any desired compound including liquid nutrients like glycerol can be fed to the culture. In combination with an optical measuring device (BioLector), online monitoring of these 44 fed-batch cultures is possible. Two Escherichia coli strains and a Hansenula polymorpha strain were successfully cultivated in the new FB-MTP. As a positive impact of the fed-batch mode on the used strains, a fourfold increase in product formation was observed for E. coli. For H. polymorpha, the use of fed-batch mode resulted in a strong increase in product formation, whereas no measurable product formation was observed in batch mode. In conclusion, the newly developed fed-batch microtiter plate is a versatile, easy-to-use, disposable system to perform fed-batch cultivations at small scale. Screening cultures in high-throughput under online monitoring are possible similar to cultivations under production conditions.

Dual role of transcription and transcript stability in the regulation of gene expression in Escherichia coli cells cultured on glucose at different growth rates

Microorganisms extensively reorganize gene expression to adjust growth rate to changes in growth conditions. At the genomic scale, we measured the contribution of both transcription and transcript stability to regulating messenger RNA (mRNA) concentration in Escherichia coli. Transcriptional control was the dominant regulatory process. Between growth rates of 0.10 and 0.63 h(-1), there was a generic increase in the bulk mRNA transcription. However, many transcripts became less stable and the median mRNA half-life decreased from 4.2 to 2.8 min. This is the first evidence that mRNA turnover is slower at extremely low-growth rates. The destabilization of many, but not all, transcripts at high-growth rate correlated with transcriptional upregulation of genes encoding the mRNA degradation machinery. We identified five classes of growth-rate regulation ranging from mainly transcriptional to mainly degradational. In general, differential stability within polycistronic messages encoded by operons does not appear to be affected by growth rate. We show here that the substantial reorganization of gene expression involving downregulation of tricarboxylic acid cycle genes and acetyl-CoA synthetase at high-growth rates is controlled mainly by transcript stability. Overall, our results demonstrate that the control of transcript stability has an important role in fine-tuning mRNA concentration during changes in growth rate.

¹³C-metabolic flux analysis for Escherichia coli

(13)C-Metabolic flux analysis ((13)C-MFA) is used here to study the effects of the knockout of such genes as pgi, zwf, gnd, ppc, pck, pyk, and lpdA on the metabolic changes in Escherichia coli cultivated under aerobic condition. The metabolic regulation mechanisms were clarified by integrating such information as fermentation data, gene expression, enzyme activities, and metabolite concentrations as well the result of (13)C-MFA.

Regulation Systems of Bacteria such as Escherichia coli in Response to Nutrient Limitation and Environmental Stresses

An overview was made to understand the regulation system of a bacterial cell such as Escherichia coli in response to nutrient limitation such as carbon, nitrogen, phosphate, sulfur, ion sources, and environmental stresses such as oxidative stress, acid shock, heat shock, and solvent stresses. It is quite important to understand how the cell detects environmental signals, integrate such information, and how the cell system is regulated. As for catabolite regulation, F1,6B P (FDP), PEP, and PYR play important roles in enzyme level regulation together with transcriptional regulation by such transcription factors as Cra, Fis, CsrA, and cAMP-Crp. αKG plays an important role in the coordinated control between carbon (C)- and nitrogen (N)-limitations, where αKG inhibits enzyme I (EI) of phosphotransferase system (PTS), thus regulating the glucose uptake rate in accordance with N level. As such, multiple regulation systems are co-ordinated for the cell synthesis and energy generation against nutrient limitations and environmental stresses. As for oxidative stress, the TCA cycle both generates and scavenges the reactive oxygen species (ROSs), where NADPH produced at ICDH and the oxidative pentose phosphate pathways play an important role in coping with oxidative stress. Solvent resistant mechanism was also considered for the stresses caused by biofuels and biochemicals production in the cell.

Integrating gene and protein expression data with genome-scale metabolic networks to infer functional pathways

Background

The study of cellular metabolism in the context of high-throughput -omics data has allowed us to decipher novel mechanisms of importance in biotechnology and health. To continue with this progress, it is essential to efficiently integrate experimental data into metabolic modeling.

Results

We present here an in-silico framework to infer relevant metabolic pathways for a particular phenotype under study based on its gene/protein expression data. This framework is based on the Carbon Flux Path (CFP) approach, a mixed-integer linear program that expands classical path finding techniques by considering additional biophysical constraints. In particular, the objective function of the CFP approach is amended to account for gene/protein expression data and influence obtained paths. This approach is termed integrative Carbon Flux Path (iCFP). We show that gene/protein expression data also influences the stoichiometric balancing of CFPs, which provides a more accurate picture of active metabolic pathways. This is illustrated in both a theoretical and real scenario. Finally, we apply this approach to find novel pathways relevant in the regulation of acetate overflow metabolism in Escherichia coli. As a result, several targets which could be relevant for better understanding of the phenomenon leading to impaired acetate overflow are proposed.

Conclusions

A novel mathematical framework that determines functional pathways based on gene/protein expression data is presented and validated. We show that our approach is able to provide new insights into complex biological scenarios such as acetate overflow in Escherichia coli.

Analysis of fluorescent reporters indicates heterogeneity in glucose uptake and utilization in clonal bacterial populations

BACKGROUND

In this study, we aimed at investigating heterogeneity in the expression of metabolic genes in clonal populations of Escherichia coli growing on glucose as the sole carbon source. Different metabolic phenotypes can arise in these clonal populations through variation in the expression of glucose transporters and metabolic enzymes. First, we focused on the glucose transporters PtsG and MglBAC to analyze the diversity of glucose uptake strategies. Second, we analyzed phenotypic variation in the expression of genes involved in gluconeogenesis and acetate scavenging (as acetate is formed and excreted during bacterial growth on glucose), which can reveal, for instance, phenotypic subpopulations that cross-feed through the exchange of acetate. In these experiments, E. coli MG1655 strains containing different transcriptional GFP reporters were grown in chemostats and reporter expression was measured with flow cytometry.

RESULTS

Our results suggest heterogeneous expression of metabolic genes in bacterial clonal populations grown in glucose environments. The two glucose transport systems exhibited different level of heterogeneity. The majority of the bacterial cells expressed the reporters for both glucose transporters MglBAC and PtsG and a small fraction of cells only expressed the reporter for Mgl. At a low dilution rate, signals from transcriptional reporters for acetyl-CoA synthetase Acs and phosphoenolpyruvate carboxykinase Pck indicated that almost all cells expressed the genes that are part of acetate utilization and the gluconeogenesis pathway, respectively. Possible co-existence of two phenotypic subpopulations differing in acs expression occurred at the threshold of the switch to overflow metabolism. The overflow metabolism results in the production of acetate and has been previously reported to occur at intermediate dilution rates in chemostats with high concentration of glucose in the feed.

CONCLUSIONS

Analysis of the heterogeneous expression of reporters for genes involved in glucose and acetate metabolism raises new question whether different metabolic phenotypes are expressed in clonal populations growing in continuous cultures fed on glucose as the initially sole carbon source.

Physiological and Molecular Timing of the Glucose to Acetate Transition in Escherichia coli

The glucose-acetate transition in Escherichia coli is a classical model of metabolic adaptation. Here, we describe the dynamics of the molecular processes involved in this metabolic transition, with a particular focus on glucose exhaustion. Although changes in the metabolome were observed before glucose exhaustion, our results point to a massive reshuffling at both the transcriptome and metabolome levels in the very first min following glucose exhaustion. A new transcriptional pattern, involving a change in genome expression in one-sixth of the E. coli genome, was established within 10 min and remained stable until the acetate was completely consumed. Changes in the metabolome took longer and stabilized 40 min after glucose exhaustion. Integration of multi-omics data revealed different modifications and timescales between the transcriptome and metabolome, but both point to a rapid adaptation of less than an hour. This work provides detailed information on the order, timing and extent of the molecular and physiological events that occur during the glucose-acetate transition and that are of particular interest for the development of dynamic models of metabolism.

Catabolite regulation analysis of Escherichia coli for acetate overflow mechanism and co-consumption of multiple sugars based on systems biology approach using computer simulation

It is quite important to understand the basic principle embedded in the main metabolism for the interpretation of the fermentation data. For this, it may be useful to understand the regulation mechanism based on systems biology approach. In the present study, we considered the perturbation analysis together with computer simulation based on the models which include the effects of global regulators on the pathway activation for the main metabolism of Escherichia coli. Main focus is the acetate overflow metabolism and the co-fermentation of multiple carbon sources. The perturbation analysis was first made to understand the nature of the feed-forward loop formed by the activation of Pyk by FDP (F1,6BP), and the feed-back loop formed by the inhibition of Pfk by PEP in the glycolysis. Those together with the effect of transcription factor Cra caused by FDP level affected the glycolysis activity. The PTS (phosphotransferase system) acts as the feed-back system by repressing the glucose uptake rate for the increase in the glucose uptake rate. It was also shown that the increased PTS flux (or glucose consumption rate) causes PEP/PYR ratio to be decreased, and EIIA-P, Cya, cAMP-Crp decreased, where cAMP-Crp in turn repressed TCA cycle and more acetate is formed. This was further verified by the detailed computer simulation. In the case of multiple carbon sources such as glucose and xylose, it was shown that the sequential utilization of carbon sources was observed for wild type, while the co-consumption of multiple carbon sources with slow consumption rates were observed for the ptsG mutant by computer simulation, and this was verified by experiments. Moreover, the effect of a specific gene knockout such as Δpyk on the metabolic characteristics was also investigated based on the computer simulation.

Acetyl-phosphate is a critical determinant of lysine acetylation in E. coli

Lysine acetylation is a frequently occurring posttranslational modification in bacteria; however, little is known about its origin and regulation. Using the model bacterium Escherichia coli (E. coli), we found that most acetylation occurred at a low level and accumulated in growth-arrested cells in a manner that depended on the formation of acetyl-phosphate (AcP) through glycolysis. Mutant cells unable to produce AcP had significantly reduced acetylation levels, while mutant cells unable to convert AcP to acetate had significantly elevated acetylation levels. We showed that AcP can chemically acetylate lysine residues in vitro and that AcP levels are correlated with acetylation levels in vivo, suggesting that AcP may acetylate proteins nonenzymatically in cells. These results uncover a critical role for AcP in bacterial acetylation and indicate that most acetylation in E. coli occurs at a low level and is dynamically affected by metabolism and cell proliferation in a global, uniform manner.

Proteomic View of Interactions of Shiga Toxin-Producing Escherichia coli with the Intestinal Environment in Gnotobiotic Piglets

BACKGROUND

Shiga toxin (Stx)-producing Escherichia coli cause severe intestinal infections involving colonization of epithelial Peyer's patches and formation of attachment/effacement (A/E) lesions. These lesions trigger leukocyte infiltration followed by inflammation and intestinal hemorrhage. Systems biology, which explores the crosstalk of Stx-producing Escherichia coli with the in vivo host environment, may elucidate novel molecular pathogenesis aspects.

METHODOLOGY/PRINCIPAL FINDINGS

Enterohemorrhagic E. coli strain 86-24 produces Shiga toxin-2 and belongs to the serotype O157:H7. Bacterial cells were scrapped from stationary phase cultures (the in vitro condition) and used to infect gnotobiotic piglets via intestinal lavage. Bacterial cells isolated from the piglets' guts constituted the in vivo condition. Cell lysates were subjected to quantitative 2D gel and shotgun proteomic analyses, revealing metabolic shifts towards anaerobic energy generation, changes in carbon utilization, phosphate and ammonia starvation, and high activity of a glutamate decarboxylase acid resistance system in vivo. Increased abundance of pyridine nucleotide transhydrogenase (PntA and PntB) suggested in vivo shortage of intracellular NADPH. Abundance changes of proteins implicated in lipopolysaccharide biosynthesis (LpxC, ArnA, the predicted acyltransferase L7029) and outer membrane (OM) assembly (LptD, MlaA, MlaC) suggested bacterial cell surface modulation in response to activated host defenses. Indeed, there was evidence for interactions of innate immunity-associated proteins secreted into the intestines (GP340, REG3-γ, resistin, lithostathine, and trefoil factor 3) with the bacterial cell envelope.

SIGNIFICANCE

Proteomic analysis afforded insights into system-wide adaptations of strain 86-24 to a hostile intestinal milieu, including responses to limited nutrients and cofactor supplies, intracellular acidification, and reactive nitrogen and oxygen species-mediated stress. Protein and lipopolysaccharide compositions of the OM were altered. Enhanced expression of type III secretion system effectors correlated with a metabolic shift back to a more aerobic milieu in vivo. Apparent pathogen pattern recognition molecules from piglet intestinal secretions adhered strongly to the bacterial cell surface.

Coregulation of Terpenoid Pathway Genes and Prediction of Isoprene Production in Bacillus subtilis Using Transcriptomics

The isoprenoid pathway converts pyruvate to isoprene and related isoprenoid compounds in plants and some bacteria. Currently, this pathway is of great interest because of the critical role that isoprenoids play in basic cellular processes, as well as the industrial value of metabolites such as isoprene. Although the regulation of several pathway genes has been described, there is a paucity of information regarding system level regulation and control of the pathway. To address these limitations, we examined Bacillus subtilis grown under multiple conditions and determined the relationship between altered isoprene production and gene expression patterns. We found that with respect to the amount of isoprene produced, terpenoid genes fall into two distinct subsets with opposing correlations. The group whose expression levels positively correlated with isoprene production included dxs, which is responsible for the commitment step in the pathway, ispD, and two genes that participate in the mevalonate pathway, yhfS and pksG. The subset of terpenoid genes that inversely correlated with isoprene production included ispH, ispF, hepS, uppS, ispE, and dxr. A genome-wide partial least squares regression model was created to identify other genes or pathways that contribute to isoprene production. These analyses showed that a subset of 213 regulated genes was sufficient to create a predictive model of isoprene production under different conditions and showed correlations at the transcriptional level. We conclude that gene expression levels alone are sufficiently informative about the metabolic state of a cell that produces increased isoprene and can be used to build a model that accurately predicts production of this secondary metabolite across many simulated environmental conditions.

Metabolic engineering for high yielding L(-)-carnitine production in Escherichia coli

Background

L(-)-carnitine production has been widely studied because of its beneficial properties on various diseases and dysfunctions. Enterobacteria possess a specific biotransformation pathway which can be used for the enantioselective production of L(-)-carnitine. Although bioprocesses catalyzed by enzymes or whole cells can overcome the lack of enantioselectivity of chemical methods, current processes for L(−)-carnitine production still have severe disadvantages, such as the low yields, side reactions and the need of high catalyst concentrations and anaerobic conditions for proper expression of the biotransformation pathway. Additionally, genetically engineered strains so far constructed for L(-)-carnitine production are based on plasmids and, therefore, suffer from segregational unstability.

Results

In this work, a stable, high yielding strain for L(-)-carnitine production from low cost substrates was constructed. A metabolic engineering strategy was implemented in a multiple mutant for use in both growing and resting cells systems. The effect of mutations on gene expression and metabolism was analyzed to characterize the productivity constraints of the wild type and the overproducer strains. Precise deletion of genes which encode proteins of central and carnitine metabolisms were performed. Specifically, flux through the TCA cycle was increased by deletion of aceK (which encodes a bifunctional kinase/phosphatase which inhibits isocitrate dehydrogenase activity) and the synthesis of the by-product γ-butyrobetaine was prevented by deletion of caiA (which encodes a crotonobetainyl-CoA reductase). Both mutations led to improve the L(-)-carnitine production by 20 and 42%, respectively. Moreover, the highly regulated promoter of the cai operon was substituted by a constitutive artificial promoter increasing the biotransformation rate, even under aerobic conditions. Resting cells of the BW ΔaceK ΔcaiA p37cai strain produced 59.6 mmol l^-1 · h^-1 of L(−)-carnitine, doubling the productivity of the wild type strain. In addition, almost total conversion was attained in less than two hours without concomitant production of the side product γ–butyrobetaine.

Conclusions

L(-)-carnitine production has been enhanced by strain engineering. Metabolic engineering strategies herein implemented allowed obtaining a robust and high yielding E. coli strain. The new overproducer strain attained almost complete conversion of crotonobetaine into L(-)-carnitine with growing and resting cells, and even under aerobic conditions, overcoming the main environmental restriction to carnitine metabolism expression. So far, this is the best performing L(-)-carnitine production E. coli strain described.

Modification of glucose import capacity in Escherichia coli: physiologic consequences and utility for improving DNA vaccine production

Background

The bacterium Escherichia coli can be grown employing various carbohydrates as sole carbon and energy source. Among them, glucose affords the highest growth rate. This sugar is nowadays widely employed as raw material in industrial fermentations. When E. coli grows in a medium containing non-limiting concentrations of glucose, a metabolic imbalance occurs whose main consequence is acetate secretion. The production of this toxic organic acid reduces strain productivity and viability. Solutions to this problem include reducing glucose concentration by substrate feeding strategies or the generation of mutant strains with impaired glucose import capacity. In this work, a collection of E. coli strains with inactive genes encoding proteins involved in glucose transport where generated to determine the effects of reduced glucose import capacity on growth rate, biomass yield, acetate and production of an experimental plasmid DNA vaccine (pHN).

Results

A group of 15 isogenic derivatives of E. coli W3110 were generated with single and multiple deletions of genes encoding glucose, mannose, beta-glucoside, maltose and N-acetylglucosamine components of the phosphoenolpyruvate:sugar phosphotransferase system (PTS), as well as the galactose symporter and the Mgl galactose/glucose ABC transporter. These strains were characterized by growing them in mineral salts medium supplemented with 2.5 g/L glucose. Maximum specific rates of glucose consumption (q_s) spanning from 1.33 to 0.32 g/g h were displayed by the group of mutants and W3110, which resulted in specific growth rates ranging from 0.65-0.18 h^-1. Acetate accumulation was reduced or abolished in cultures with all mutant strains. W3110 and five selected mutant derivatives were transformed with pHN. A 3.2-fold increase in pHN yield on biomass was observed in cultures of a mutant strain with deletion of genes encoding the glucose and mannose PTS components, as well as Mgl.

Conclusions

The group of E. coli mutants generated in this study displayed a reduction or elimination of overflow metabolism and a linear correlation between q_s and the maximum specific growth rate as well as the acetate production rate. By comparing DNA vaccine production parameters among some of these mutants, it was possible to identify a near-optimal glucose import rate value for this particular application. The strains employed in this study should be a useful resource for studying the effects of different predefined q_s values on production capacity for various biotechnological products.

Bringing together models from bottom-up and top-down approaches: an application for growth of Escherichia coli on different carbohydrates

Modeling in systems biology follows two lines: a data driven top-down approach that integrates experimental data from various "omics" technologies and a model based bottom-up approach where the model structure is given and kinetic parameters are chosen in such a way that an experimental observation can be reproduced quantitatively or qualitatively. Mathematical models are frequently used to elucidate cellular design principles in order to understand complex biochemical networks better. To show that both approaches lead to a consistent description of cellular dynamics, mathematical models from both approaches are explored. On the level of transcription factor activities a sufficient qualitative agreement is observed. Experimental data for the classical growth experiment of Escherichia coli on two carbon sources, glucose and lactose is available to set up the data driven model and to support the theoretical findings from the bottom-up approach.

Analysis and Design of Stimulus Response Curves of E. coli

Metabolism and signalling are tightly coupled in bacteria. Combining several theoretical approaches, a core model is presented that describes transcriptional and allosteric control of glycolysis in Escherichia coli. Experimental data based on microarrays, signaling components and extracellular metabolites are used to estimate kinetic parameters. A newly designed strain was used that adjusts the incoming glucose flux into the system and allows a kinetic analysis. Based on the results, prediction for intracelluar metabolite concentrations over a broad range of the growth rate could be performed and compared with data from literature.

Influence of the RelA Activity on E. coli Metabolism by Metabolite Profiling of Glucose-Limited Chemostat Cultures

Metabolite profiling of E. coli W3110 and the isogenic ΔrelA mutant cells was used to characterize the RelA-dependent stringent control of metabolism under different growth conditions. Metabolic profiles were obtained by gas chromatography–mass spectrometry (GC-MS) analysis and revealed significant differences between E. coli strains grown at different conditions. Major differences between the two strains were assessed in the levels of amino acids and fatty acids and their precursor metabolites, especially when growing at the lower dilution rates, demonstrating differences in their metabolic behavior. Despite the fatty acid biosynthesis being the most affected due to the lack of the RelA activity, other metabolic pathways involving succinate, lactate and threonine were also affected. Overall, metabolite profiles indicate that under nutrient-limiting conditions the RelA-dependent stringent response may be elicited and promotes key changes in the E. coli metabolism.

Metabolic responses to recombinant bioprocesses in Escherichia coli

Escherichia coli has been widely used for the production of recombinant proteins. However, the unbalances between host metabolism and recombinant biosynthesis continue to hamper the efficiency of these recombinant bioprocesses. The additional drainage of biosynthetic precursors toward recombinant processes burdens severely the metabolism of cells that, ultimately, elicits a series of stress responses, reducing biomass growth and recombinant protein production. Several strategies to overcome these metabolic limitations have been implemented; however, in most cases, improvements in recombinant protein expression were achieved at the expense of biomass growth arrest, which significantly hampers the efficiency of recombinant bioprocesses. With the advent of high throughput techniques and modelling approaches that provide a system-level understanding of the cellular systems, it is now expected that new advances in recombinant bioprocesses are achieved. By providing means to deal with these systems, our understanding on the metabolic behaviour of recombinant cells will advance and can be further explored to the design of suitable hosts and more efficient and cost-effective bioprocesses. Here, we review the major metabolic responses associated with recombinant processes and the engineering strategies relevant to overcome these stresses. Moreover, the advantages of applying systems levels engineering strategies to enhance recombinant protein production in E. coli cells are discussed and future perspectives on the advances of mathematical modelling approaches to study these systems are exposed.

Extensive exometabolome analysis reveals extended overflow metabolism in various microorganisms

Overflow metabolism is well known for yeast, bacteria and mammalian cells. It typically occurs under glucose excess conditions and is characterized by excretions of by-products such as ethanol, acetate or lactate. This phenomenon, also denoted the short-term Crabtree effect, has been extensively studied over the past few decades, however, its basic regulatory mechanism and functional role in metabolism is still unknown. Here we present a comprehensive quantitative and time-dependent analysis of the exometabolome of Escherichia coli, Corynebacterium glutamicum, Bacillus licheniformis, and Saccharomyces cerevisiae during well-controlled bioreactor cultivations. Most surprisingly, in all cases a great diversity of central metabolic intermediates and amino acids is found in the culture medium with extracellular concentrations varying in the micromolar range. Different hypotheses for these observations are formulated and experimentally tested. As a result, the intermediates in the culture medium during batch growth must originate from passive or active transportation due to a new phenomenon termed “extended” overflow metabolism. Moreover, we provide broad evidence that this could be a common feature of all microorganism species when cultivated under conditions of carbon excess and non-inhibited carbon uptake. In turn, this finding has consequences for metabolite balancing and, particularly, for intracellular metabolite quantification and ¹³C-metabolic flux analysis.

Adaptation of glycolysis and growth to acetate in Sporolactobacillus sp. Y2-8

Exogenous addition of a low concentration of acetate (2 g/L sodium acetate) effectively decreased acetic acid excretion and lowered the ATP content in Sporolactobacillus sp. Y2-8 without any growth defect although the acetate could not be utilized at an initial glucose concentration of 150 g/L. This induced an enhanced glycolytic flux with increased specific activities of hexokinase and phosphofructokinase, probably to compensate for the lowered efficiency of ATP production. However, with increasing concentrations (5 g/L sodium acetate), acetate was utilized first before being produced again, causing a growth lag at the transition. Glucose consumption was also reduced at high acetate concentrations, resulting in decreased D-lactic acid production. These results demonstrate that acetate plays a significant role in regulating glycolysis and growth of Sporolactobacillus.

Comparison and applications of label-free absolute proteome quantification methods on Escherichia coli

Three different label-free proteome quantification methods--APEX, emPAI and iBAQ--were evaluated to measure proteome-wide protein concentrations in the cell. All the methods were applied to a sample from Escherichia coli chemostat culture. A Pearson squared correlation of approximately 0.6 among the three quantification methods was demonstrated. Importantly, the sum of quantified proteins by iBAQ and emPAI corresponded with the Lowry total protein quantification, demonstrating applicability of label-free methods for an accurate calculation of protein concentrations at the proteome level. The iBAQ method showed the best correlation between biological replicates, a normal distribution among all protein abundances, and the lowest variation among ribosomal protein abundances, which are expected to have equal amounts. Absolute quantitative proteome data enabled us to evaluate metabolic cost for protein synthesis and apparent catalytic activities of enzymes by integration with flux analysis. All the methods demonstrated similar ATP costs for protein synthesis for different cellular processes and that costs for expressing biomass synthesis related proteins were higher than those for energy generation. Importantly, catalytic activities of energy metabolism enzymes were an order or two higher than those of monomer synthesis. Interestingly, a staircase-like protein expression was demonstrated for most of the transcription units.

Aerobic glycolysis: meeting the metabolic requirements of cell proliferation

Warburg's observation that cancer cells exhibit a high rate of glycolysis even in the presence of oxygen (aerobic glycolysis) sparked debate over the role of glycolysis in normal and cancer cells. Although it has been established that defects in mitochondrial respiration are not the cause of cancer or aerobic glycolysis, the advantages of enhanced glycolysis in cancer remain controversial. Many cells ranging from microbes to lymphocytes use aerobic glycolysis during rapid proliferation, which suggests it may play a fundamental role in supporting cell growth. Here, we review how glycolysis contributes to the metabolic processes of dividing cells. We provide a detailed accounting of the biosynthetic requirements to construct a new cell and illustrate the importance of glycolysis in providing carbons to generate biomass. We argue that the major function of aerobic glycolysis is to maintain high levels of glycolytic intermediates to support anabolic reactions in cells, thus providing an explanation for why increased glucose metabolism is selected for in proliferating cells throughout nature.

Analyzing the metabolic stress response of recombinant Escherichia coli cultures expressing human interferon-beta in high cell density fed batch cultures using time course transcriptomic data

Fed batch cultures expressing recombinant interferon beta under the T7 promoter were run with different exponential feeding rates of a complex substrate and induced at varying cell densities. Post-induction profiles of the specific product formation rates showed a strong dependence on the specific growth rate with the maximum product yield obtained at 0.2 h(-1). A study of the relative transcriptomic profiles as a function of pre-induction μ was therefore done to provide insight into the role of cellular physiology in enhancing recombinant protein expression. Hierarchical clustering analysis of the significantly regulated genes allowed us to identify biologically important groups of genes which fall under specific master regulators. The groups were: rpoH, ArcB, CreB, Lrp, RelA, Fis and Hfq. The response of these regulators, which exert a feedback control on the growth and product formation rates correlated well with the expression levels obtained. Thus at the optimum pre-induction μ, the alternative sigma factors and ribosomal machinery genes did not get depressed till the 6th hour post-induction unlike at other specific growth rates, demonstrating a critical role for the genes in sustaining recombinant protein expression.

Acetate scavenging activity in Escherichia coli: interplay of acetyl-CoA synthetase and the PEP-glyoxylate cycle in chemostat cultures

Impairment of acetate production in Escherichia coli is crucial for the performance of many biotechnological processes. Aerobic production of acetate (or acetate overflow) results from changes in the expression of central metabolism genes. Acetyl-CoA synthetase scavenges extracellular acetate in glucose-limited cultures. Once converted to acetyl-CoA, it can be catabolized by the tricarboxylic acid cycle or the glyoxylate pathway. In this work, we assessed the significance of these pathways on acetate overflow during glucose excess and limitation. Gene expression, enzyme activities, and metabolic fluxes were studied in E. coli knock-out mutants related to the glyoxylate pathway operon and its regulators. The relevance of post-translational regulation by AceK-mediated phosphorylation of isocitrate dehydrogenase for pathway functionality was underlined. In chemostat cultures performed at increasing dilution rates, acetate overflow occurs when growing over a threshold glucose uptake rate. This threshold was not affected in a glyoxylate-pathway-deficient strain (lacking isocitrate lyase, the first enzyme of the pathway), indicating that it is not relevant for acetate overflow. In carbon-limited chemostat cultures, gluconeogenesis (maeB, sfcA, and pck), the glyoxylate operon and, especially, acetyl-CoA synthetase are upregulated. A mutant in acs (encoding acetyl-CoA synthetase) produced acetate at all dilution rates. This work demonstrates that, in E. coli, acetate production occurs at all dilution rates and that overflow is the result of unbalanced synthesis and scavenging activities. The over-expression of acetyl-CoA synthetase by cAMP-CRP-dependent induction limits this phenomenon in cultures consuming glucose at low rate, ensuring the recycling of the acetyl-CoA and acetyl-phosphate pools, although establishing an energy-dissipating substrate cycle.

Decrease of energy spilling in Escherichia coli continuous cultures with rising specific growth rate and carbon wasting

Background

Growth substrates, aerobic/anaerobic conditions, specific growth rate (μ) etc. strongly influence Escherichia coli cell physiology in terms of cell size, biomass composition, gene and protein expression. To understand the regulation behind these different phenotype properties, it is useful to know carbon flux patterns in the metabolic network which are generally calculated by metabolic flux analysis (MFA). However, rarely is biomass composition determined and carbon balance carefully measured in the same experiments which could possibly lead to distorted MFA results and questionable conclusions. Therefore, we carried out both detailed carbon balance and biomass composition analysis in the same experiments for more accurate quantitative analysis of metabolism and MFA.

Results

We applied advanced continuous cultivation methods (A-stat and D-stat) to continuously monitor E. coli K-12 MG1655 flux and energy metabolism dynamic responses to change of μ and glucose-acetate co-utilisation. Surprisingly, a 36% reduction of ATP spilling was detected with increasing μ and carbon wasting to non-CO₂ by-products under constant biomass yield. The apparent discrepancy between constant biomass yield and decline of ATP spilling could be explained by the rise of carbon wasting from 3 to 11% in the carbon balance which was revealed by the discovered novel excretion profile of E. coli pyrimidine pathway intermediates carbamoyl-phosphate, dihydroorotate and orotate. We found that carbon wasting patterns are dependent not only on μ, but also on glucose-acetate co-utilisation capability. Accumulation of these compounds was coupled to the two-phase acetate accumulation profile. Acetate overflow was observed in parallel with the reduction of TCA cycle and glycolysis fluxes, and induction of pentose phosphate pathway.

Conclusions

It can be concluded that acetate metabolism is one of the major regulating factors of central carbon metabolism. More importantly, our model calculations with actual biomass composition and detailed carbon balance analysis in steady state conditions with -omics data comparison demonstrate the importance of a comprehensive systems biology approach for more advanced understanding of metabolism and carbon re-routing mechanisms potentially leading to more successful metabolic engineering.

Quantification of metabolic limitations during recombinant protein production in Escherichia coli

Escherichia coli is one of the major microorganisms for recombinant protein production because it has been best characterized in terms of molecular genetics and physiology, and because of the availability of various expression vectors and strains. The synthesis of proteins is one of the most energy consuming processes in the cell, with the result that cellular energy supply may become critical. Indeed, the so called metabolic burden of recombinant protein synthesis was reported to cause alterations in the operation of the host's central carbon metabolism. To quantify these alterations in E. coli metabolism in dependence of the rate of recombinant protein production, (13)C-tracer-based metabolic flux analysis in differently induced cultures was used. To avoid dilution of the (13)C-tracer signal by the culture history, the recombinant protein produced was used as a flux probe, i.e., as a read out of intracellular flux distributions. In detail, an increase in the generation rate rising from 36 mmol(ATP)g(CDW)(-1)h(-1) for the reference strain to 45 mmol(ATP)g(CDW)(-1)h(-1) for the highest yielding strain was observed during batch cultivation. Notably, the flux through the TCA cycle was rather constant at 2.5±0.1 mmol g(CDW)(-1)h(-1), hence was independent of the induced strength for gene expression. E. coli compensated for the additional energy demand of recombinant protein synthesis by reducing the biomass formation to almost 60%, resulting in excess NADPH. Speculative, this excess NADPH was converted to NADH via the soluble transhydrogenase and subsequently used for ATP generation in the electron transport chain. In this study, the metabolic burden was quantified by the biomass yield on ATP, which constantly decreased from 11.7g(CDW)mmol(ATP)(-1) for the reference strain to 4.9g(CDW)mmol(ATP)(-1) for the highest yielding strain. The insights into the operation of the metabolism of E. coli during recombinant protein production might guide the optimization of microbial hosts and fermentation conditions.