Feedback regulation and coordination of the main metabolism for bacterial growth and metabolic engineering for amino acid fermentation.

Living organisms such as bacteria are often exposed to continuous changes in the nutrient availability in nature. Therefore, bacteria must constantly monitor the environmental condition, and adjust the metabolism quickly adapting to the change in the growth condition. For this, bacteria must orchestrate (coordinate and integrate) the complex and dynamically changing information on the environmental condition. In particular, the central carbon metabolism (CCM), monomer synthesis, and macromolecular synthesis must be coordinately regulated for the efficient growth. It is a grand challenge in bioscience, biotechnology, and synthetic biology to understand how living organisms coordinate the metabolic regulation systems. Here, we consider the integrated sensing of carbon sources by the phosphotransferase system (PTS), and the feed-forward/feedback regulation systems incorporated in the CCM in relation to the pool sizes of flux-sensing metabolites and αketoacids. We also consider the metabolic regulation of amino acid biosynthesis (as well as purine and pyrimidine biosyntheses) paying attention to the feedback control systems consisting of (fast) enzyme level regulation with (slow) transcriptional regulation. The metabolic engineering for the efficient amino acid production by bacteria such as Escherichia coli and Corynebacterium glutamicum is also discussed (in relation to the regulation mechanisms). The amino acid synthesis is important for determining the rate of ribosome biosynthesis. Thus, the growth rate control (growth law) is further discussed on the relationship between (p)ppGpp level and the ribosomal protein synthesis.

Computer-Aided Rational Design of Efficient NADPH Production System by Escherichia coli pgi Mutant Using a Mixture of Glucose and Xylose.

Lignocellulosic biomass can be hydrolyzed into two major sugars of glucose and xylose, and thus the strategy for the efficient consumption of both sugars is highly desirable. NADPH is the essential molecule for the production of industrially important value-added chemicals, and thus its availability is quite important. Escherichia coli mutant lacking the pgi gene encoding phosphoglucose isomerase (Pgi) has been preferentially used to overproduce the NADPH. However, there exists a disadvantage that the cell growth rate becomes low for the mutant grown on glucose. This limits the efficient NADPH production, and therefore, it is quite important to investigate how addition of different carbon source such as xylose (other than glucose) effectively improves the NADPH production. In this study, we have developed a kinetic model to propose an efficient NADPH production system using E. coli pgi-knockout mutant with a mixture of glucose and xylose. The proposed system adds xylose to glucose medium to recover the suppressed growth of the pgi mutant, and determines the xylose content to maximize the NADPH productivity. Finally, we have designed a mevalonate (MVA) production system by implementing ArcA overexpression into the pgi-knockout mutant using a mixture of glucose and xylose. In addition to NADPH overproduction, the accumulation of acetyl-CoA (AcCoA) is necessary for the efficient MVA production. In the present study, therefore, we considered to overexpress ArcA, where ArcA overexpression suppresses the TCA cycle, causing the overflow of AcCoA, a precursor of MVA. We predicted the xylose content that maximizes the MVA production. This approach demonstrates the possibility of a great progress in the computer-aided rational design of the microbial cell factories for useful metabolite production.

Redox rebalance against genetic perturbations and modulation of central carbon metabolism by the oxidative stress regulation.

The micro-aerophilic organisms and aerobes as well as yeast and higher organisms have evolved to gain energy through respiration (via oxidative phosphorylation), thereby enabling them to grow much faster than anaerobes. However, during respiration, reactive oxygen species (ROSs) are inherently (inevitably) generated, and threaten the cell's survival. Therefore, living organisms (or cells) must furnish the potent defense systems to keep such ROSs at harmless level, where the cofactor balance plays crucial roles. Namely, NADH is the source of energy generation (catabolism) in the respiratory chain reactions, through which ROSs are generated, while NADPH plays important roles not only for the cell synthesis (anabolism) but also for detoxifying ROSs. Therefore, the cell must rebalance the redox ratio by modulating the fluxes of the central carbon metabolism (CCM) by regulating the multi-level regulation machinery upon genetic perturbations and the change in the growth conditions. Here, we discuss about how aerobes accomplish such cofactor homeostasis against redox perturbations. In particular, we consider how single-gene mutants (including pgi, pfk, zwf, gnd and pyk mutants) modulate their metabolisms in relation to cofactor rebalance (and also by adaptive laboratory evolution). We also discuss about how the overproduction of NADPH (by the pathway gene mutation) can be utilized for the efficient production of useful value-added chemicals such as medicinal compounds, polyhydroxyalkanoates, and amino acids, all of which require NADPH in their synthetic pathways. We then discuss about the metabolic responses against oxidative stress, where αketoacids play important roles not only for the coordination between catabolism and anabolism, but also for detoxifying ROSs by non-enzymatic reactions, as well as for reducing the production of ROSs by repressing the activities of the TCA cycle and respiration (via carbon catabolite repression). Thus, we discuss about the mechanisms (basic strategies) that modulate the metabolism from respiration to respiro-fermentative metabolism causing overflow, based on the role of Pyk activity, affecting the NADPH production at the oxidative pentose phosphate (PP) pathway, and the roles of αketoacids for the change in the source of energy generation from the oxidative phosphorylation to the substrate level phosphorylation.

Regulation of glycolytic flux and overflow metabolism depending on the source of energy generation for energy demand.

Overflow metabolism is a common phenomenon observed at higher glycolytic flux in many bacteria, yeast (known as Crabtree effect), and mammalian cells including cancer cells (known as Warburg effect). This phenomenon has recently been characterized as the trade-offs between protein costs and enzyme efficiencies based on coarse-graining approaches. Moreover, it has been recognized that the glycolytic flux increases as the source of energy generation changes from energetically efficient respiration to inefficient respiro-fermentative or fermentative metabolism causing overflow metabolism. It is highly desired to clarify the metabolic regulation mechanisms behind such phenomena. Metabolic fluxes are located on top of the hierarchical regulation systems, and represent the outcome of the integrated response of all levels of cellular regulation systems. In the present article, we discuss about the different levels of regulation systems for the modulation of fluxes depending on the growth rate, growth condition such as oxygen limitation that alters the metabolism towards fermentation, and genetic perturbation affecting the source of energy generation from respiration to respiro-fermentative metabolism in relation to overflow metabolism. The intracellular metabolite of the upper glycolysis such as fructose 1,6-bisphosphate (FBP) plays an important role not only for flux sensing, but also for the regulation of the respiratory activity either directly or indirectly (via transcription factors) at higher growth rate. The glycolytic flux regulation is backed up (enhanced) by unphosphorylated EIIA and HPr of the phosphotransferase system (PTS) components, together with the sugar-phosphate stress regulation, where the transcriptional regulation is further modulated by post-transcriptional regulation via the degradation of mRNA (stability of mRNA) in Escherichia coli. Moreover, the channeling may also play some role in modulating the glycolytic cascade reactions.

Modeling and simulation of the redox regulation of the metabolism in Escherichia coli at different oxygen concentrations.

BACKGROUND: Microbial production of biofuels and biochemicals from renewable feedstocks has received considerable recent attention from environmental protection and energy production perspectives. Many biofuels and biochemicals are produced by fermentation under oxygen-limited conditions following initiation of aerobic cultivation to enhance the cell growth rate. Thus, it is of significant interest to investigate the effect of dissolved oxygen concentration on redox regulation in Escherichia coli, a particularly popular cellular factory due to its high growth rate and well-characterized physiology. For this, the systems biology approach such as modeling is powerful for the analysis of the metabolism and for the design of microbial cellular factories. RESULTS: Here, we developed a kinetic model that describes the dynamics of fermentation by taking into account transcription factors such as ArcA/B and Fnr, respiratory chain reactions and fermentative pathways, and catabolite regulation. The hallmark of the kinetic model is its ability to predict the dynamics of metabolism at different dissolved oxygen levels and facilitate the rational design of cultivation methods. The kinetic model was verified based on the experimental data for a wild-type E. coli strain. The model reasonably predicted the metabolic characteristics and molecular mechanisms of fnr and arcA gene-knockout mutants. Moreover, an aerobic-microaerobic dual-phase cultivation method for lactate production in a pfl-knockout mutant exhibited promising yield and productivity. CONCLUSIONS: It is quite important to understand metabolic regulation mechanisms from both scientific and engineering points of view. In particular, redox regulation in response to oxygen limitation is critically important in the practical production of biofuel and biochemical compounds. The developed model can thus be used as a platform for designing microbial factories to produce a variety of biofuels and biochemicals.

Development of an accurate kinetic model for the central carbon metabolism of Escherichia coli.

BACKGROUND: A kinetic model provides insights into the dynamic response of biological systems and predicts how their complex metabolic and gene regulatory networks generate particular functions. Of many biological systems, Escherichia coli metabolic pathways have been modeled extensively at the enzymatic and genetic levels, but existing models cannot accurately reproduce experimental behaviors in a batch culture, due to the inadequate estimation of a specific cell growth rate and a large number of unmeasured parameters. RESULTS: In this study, we developed a detailed kinetic model for the central carbon metabolism of E. coli in a batch culture, which includes the glycolytic pathway, tricarboxylic acid cycle, pentose phosphate pathway, Entner-Doudoroff pathway, anaplerotic pathway, glyoxylate shunt, oxidative phosphorylation, phosphotransferase system (Pts), non-Pts and metabolic gene regulations by four protein transcription factors: cAMP receptor, catabolite repressor/activator, pyruvate dehydrogenase complex repressor and isocitrate lyase regulator. The kinetic parameters were estimated by a constrained optimization method on a supercomputer. The model estimated a specific growth rate based on reaction kinetics and accurately reproduced the dynamics of wild-type E. coli and multiple genetic mutants in a batch culture. CONCLUSIONS: This model overcame the intrinsic limitations of existing kinetic models in a batch culture, predicted the effects of multilayer regulations (allosteric effectors and gene expression) on central carbon metabolism and proposed rationally designed fast-growing cells based on understandings of molecular processes.

S-system-based analysis of the robust properties common to many biochemical network models.

Robustness is a key feature to characterize the adaptation of organisms to changes in their internal and external environments. A broad range of kinetic or dynamic models of biochemical systems have been developed. Robustness analyses are attractive for exploring some common properties of many biochemical models. To reveal such features, we transform different types of mathematical equations into a standard or intelligible formula and use the multiple parameter sensitivity (MPS) to identify some factors critically responsible for the total robustness to many perturbations. The MPS would be determined by the top quarter of the highly sensitive parameters rather than the single parameter with the maximum sensitivity. The MPS did not show any correlation to the network size. The MPS is closely related to the standard deviation of the sensitivity profile. A decrease in the standard deviation enhanced the total robustness, which shows the hallmark of distributed robustness that many factors (pathways) involve the total robustness.

Metabolic flux analysis for Escherichia coli by flux balance analysis.

Conventional metabolic flux analysis (MFA) of Escherichia coli wild type and of pathway gene knockout mutants cultivated under anaerobic condition is explained in detail in this chapter. To place the MFA results into the context of the literature, the regulation of central carbon metabolism in terms of catabolite regulation by the phosphotransferase system (PTS) and the response to oxygen limitations via global regulators is reviewed. The effects of gene deletions such as pflA, pta, ppc, pykF, adhE, and ldhA on the metabolic network are presented. Moreover, for the pflA mutant the effects of various carbon sources were quantified. The chapter thereby contributes to the discussion of metabolic network function and the design of microbial cell factories.

¹³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.

CADLIVE toolbox for MATLAB: automatic dynamic modeling of biochemical networks with comprehensive system analysis.

Mathematical modeling has become a standard technique to understand the dynamics of complex biochemical systems. To promote the modeling, we had developed the CADLIVE dynamic simulator that automatically converted a biochemical map into its associated mathematical model, simulated its dynamic behaviors and analyzed its robustness. To enhance the feasibility by CADLIVE and extend its functions, we propose the CADLIVE toolbox available for MATLAB, which implements not only the existing functions of the CADLIVE dynamic simulator, but also the latest tools including global parameter search methods with robustness analysis. The seamless, bottom-up processes consisting of biochemical network construction, automatic construction of its dynamic model, simulation, optimization, and S-system analysis greatly facilitate dynamic modeling, contributing to the research of systems biology and synthetic biology. This application can be freely downloaded from http://www.cadlive.jp/CADLIVE_MATLAB/ together with an instruction.

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.

Importance of understanding the main metabolic regulation in response to the specific pathway mutation for metabolic engineering of Escherichia coli.

Recent metabolic engineering practice was briefly reviewed in particular for the useful metabolite production such as natural products and biofuel productions. With the emphasis on systems biology approach, the metabolic regulation of the main metabolic pathways in E. coli was discussed from the points of view of enzyme level (allosteric and phosphorylation/ dephosphorylation) regulation, and gene level (transcriptional) regulation. Then the effects of the specific pathway gene knockout such as pts, pgi, zwf, gnd, pyk, ppc, pckA, lpdA, pfl gene knockout on the metabolism in E. coli were overviewed from the systems biology point of view with possible application for strain improvement point.

Metabolic regulation in Escherichia coli in response to culture environments via global regulators.

One of the ultimate goal of systems biology is to realize a virtual cell system in the computer. If this could be attained, it might be possible, for example, to quantitatively predict the effects of a culture environment and/or the removal/inactivation of specific genes on the metabolism without conducting many experiments. Thus, it may be possible to design cells, e.g., for more efficient production of a specific metabolite. To achieve this, it is important to properly understand the metabolic regulation mechanism and to develop a robust model by incorporating gene-level regulation into the enzymatic reaction model with the integration of different levels of information. However, the metabolic regulation in response to the change in culture environment is itself not well understood. Here, we overview how the culture environment affects cell metabolism via global regulators with sigma factors, considering the effects of carbon, nitrogen, and phosphate sources as well as oxygen, temperature, pH, and nutrient stress, etc., on transcriptional regulation. A variety of controlled strategies for the specific stimuli imposed on the cell appear to exist, and some of the regulations are interconnected by gene level regulation. Quantitative modeling for these regulation mechanisms is critical for efficient metabolic engineering of a cell.

Kinetics of dextran-independent α-(1→3)-glucan synthesis by Streptococcus sobrinus glucosyltransferase I.

Glucosyltransferase (GTF)-I from cariogenic Streptococcus sobrinus elongates the α-(1→3)-linked glucose polymer branches on the primer dextran bound to the C-terminal glucan-binding domain. We investigated the GTF-I-catalyzed glucan synthesis reaction in the absence of the primer dextran. The time course of saccharide production during dextran-independent glucan synthesis from sucrose was analyzed. Fructose and glucose were first produced by the sucrose hydrolysis. Leucrose was subsequently produced, followed by insoluble glucan [α-(1→3)-linked glucose polymers] after a lag phase. High levels of intermediate nigerooligosaccharide series accumulation were characteristically not observed during the lag phase. The results from the enzymatic activity of the acceptor reaction for the nigerooligosaccharide with a degree of polymerization of 2-6 and methyl α-D-glucopyranoside as a glucose analog indicate that the activity increased with an increase in the degree of polymerization. The production of insoluble glucan was numerically simulated using the fourth-order Runge-Kutta method with the kinetic parameters estimated from the enzyme assay. The simulated time course provided a profile similar to that of experimental data. These results define the relationship between the kinetic properties of GTF-I and the time course of saccharide production. These results are discussed with respect to a mechanism that underlies efficient glucan synthesis.