Responses of the central metabolism in Escherichia coli to phosphoglucose isomerase and glucose-6-phosphate dehydrogenase knockouts

Identification of a novel NADPH generation reaction in the pentose phosphate pathway in Escherichia coli using mBFP

NADPH is a redox cofactor that drives the anabolic reactions. Although major NADPH generation reactions have been identified in Escherichia coli, some minor reactions have not been identified. In the present study, we explored novel NADPH generation reactions by monitoring the fluorescence dynamics after the addition of carbon sources to starved cells, using a metagenome-derived blue fluorescent protein (mBFP) as an intracellular NADPH reporter. Perturbation analyses were performed on a glucose-6-phosphate isomerase (PGI) deletion strain and its parental strain. Interestingly, mBFP fluorescence increased not only in the parental strain but also in the ΔPGI strain after the addition of xylose. Because the ΔPGI strain cannot metabolize xylose through the oxidative pentose phosphate pathway, this suggests that an unexpected NADPH generation reaction contributes to an increase in fluorescence. To unravel this mystery, we deleted the NADPH generation enzymes including transhydrogenase, isocitrate dehydrogenase, NADP+-dependent malic enzyme, glucose-6-phosphate dehydrogenase (G6PDH), and 6-phosphogluconate dehydrogenase (6PGDH) in the ΔPGI strain, and revealed that G6PDH and 6PGDH contribute to an increase in fluorescence under xylose conditions. In vitro assays using purified enzymes showed that G6PDH can produce NADPH using erythrose-4-phosphate (E4P) as a substitute for glucose-6-phosphate. Because the Km (0.65 mM) for E4P was much higher than the reported intracellular E4P concentrations in E. coli, little E4P must be metabolized through this bypass in the parental strain. However, the flux would increase when E4P accumulates in the cells owing to genetic modifications. This finding provides a metabolic engineering strategy for generating NADPH to produce useful compounds using xylose as a carbon source.IMPORTANCEBecause NADPH is consumed during the synthesis of various useful compounds, enhancing NADPH regeneration is highly desirable in metabolic engineering. In this study, we explored novel NADPH generation reactions in Escherichia coli using a fluorescent NADPH reporter and found that glucose-6-phosphate dehydrogenase can produce NADPH using erythrose-4-phosphate as a substrate under xylose conditions. Xylose is an abundant sugar in nature and is an attractive carbon source for bioproduction. Therefore, this finding contributes to novel pathway engineering strategies using a xylose carbon source in E. coli to produce useful compounds that consume NADPH for their synthesis.

Salmonella Typhimurium screen identifies shifts in mixed-acid fermentation during gut colonization

How enteric pathogens adapt their metabolism to a dynamic gut environment is not yet fully understood. To investigate how Salmonella enterica Typhimurium (S.Tm) colonizes the gut, we conducted an in vivo transposon mutagenesis screen in a gnotobiotic mouse model. Our data implicate mixed-acid fermentation in efficient gut-luminal growth and energy conservation throughout infection. During initial growth, the pathogen utilizes acetate fermentation and fumarate respiration. After the onset of gut inflammation, hexoses appear to become limiting, as indicated by carbohydrate analytics and the increased need for gluconeogenesis. In response, S.Tm adapts by ramping up ethanol fermentation for redox balancing and supplying the TCA cycle with α-ketoglutarate for additional energy. Our findings illustrate how S.Tm flexibly adapts mixed fermentation and its use of the TCA cycle to thrive in the changing gut environment. Similar metabolic wiring in other pathogenic Enterobacteriaceae may suggest a broadly conserved mechanism for gut colonization.

From metagenomes to metabolism: Systematically assessing the metabolic flux feasibilities for "Candidatus Accumulibacter" species during anaerobic substrate uptake

With the rapid growing availability of metagenome assembled genomes (MAGs) and associated metabolic models, the identification of metabolic potential in individual community members has become possible. However, the field still lacks an unbiassed systematic evaluation of the generated metagenomic information to uncover not only metabolic potential, but also feasibilities of these models under specific environmental conditions. In this study, we present a systematic analysis of the metabolic potential in species of "Candidatus Accumulibacter", a group of polyphosphate-accumulating organisms (PAOs). We constructed a metabolic model of the central carbon metabolism and compared the metabolic potential among available MAGs for "Ca. Accumulibacter" species. By combining Elementary Flux Modes Analysis (EFMA) with max-min driving force (MDF) optimization, we obtained all possible flux distributions of the metabolic network and calculated their individual thermodynamic feasibility. Our findings reveal significant variations in the metabolic potential among "Ca. Accumulibacter" MAGs, particularly in the presence of anaplerotic reactions. EFMA revealed 700 unique flux distributions in the complete metabolic model that enable the anaerobic uptake of acetate and its conversion into polyhydroxyalkanoates (PHAs), a well-known phenotype of "Ca. Accumulibacter". However, thermodynamic constraints narrowed down this solution space to 146 models that were stoichiometrically and thermodynamically feasible (MDF > 0 kJ/mol), of which only 8 were strongly feasible (MDF > 7 kJ/mol). Notably, several novel flux distributions for the metabolic model were identified, suggesting putative, yet unreported, functions within the PAO communities. Overall, this work provides valuable insights into the metabolic variability among "Ca. Accumulibacter" species and redefines the anaerobic metabolic potential in the context of phosphate removal. More generally, the integrated workflow presented in this paper can be applied to any metabolic model obtained from a MAG generated from microbial communities to objectively narrow the expected phenotypes from community members.

Metabolic engineering of Escherichia coli for enhanced production of 1,3-butanediol from glucose

1,3-Butanediol (1,3-BDO) finds versatile applications in the cosmetic, chemical, and food industries. This study focuses on the metabolic engineering of Escherichia coli K12 to achieve efficient production of 1,3-BDO from glucose via acetoacetyl-CoA, 3-hydroxybutyryl-CoA, and 3-hydroxybutyraldehyde. The accumulation of an intermediary metabolite (pyruvate) and a byproduct (3-hydroxybutyric acid) was reduced by disruption of the negative transcription factor (PdhR) for pyruvate dehydrogenase complex (PDHc) and employing an efficient alcohol dehydrogenase (YjgB), respectively. Additionally, to improve NADPH availability, carbon flux was redirected from the Embden-Meyerhof-Parnas (EMP) pathway to the Entner-Doudoroff (ED) pathway. One resulting strain achieved a record-high titer of 790 mM (∼71.1 g/L) with a yield of 0.65 mol/mol for optically pure (R)-1,3-BDO, with an enantiomeric excess (e.e.) value of 98.5 %. These findings are useful in the commercial production of 1,3-BDO and provide valuable insights into the development of an efficient cell factory for other acetyl-CoA derivatives.

Cytosolic phosphoglucose isomerase is essential for microsporogenesis and embryogenesis in Arabidopsis

Phosphoglucose isomerase (PGI) catalyzes the interconversion of fructose-6-phosphate and glucose-6-phosphate, which impacts cell carbon metabolic flow. Arabidopsis (Arabidopsis thaliana) contains two nuclear PGI genes respectively encoding plastidial PGI1 and cytosolic PGI (cPGI). The loss of PGI1 impairs the conversion of F6P of the Calvin-Benson cycle to G6P for the synthesis of transitory starch in leaf chloroplasts. Since cpgi knockout mutants have not yet been obtained, they are thought to be lethal. The cpgi lethality can be rescued by expressing CaMV 35S promoter (p35S)-driven cPGI; however, the complemented line is completely sterile due to pollen degeneration. Here, we generated a cpgi mutant expressing p35S::cPGI-YFP in which YFP fluorescence in developing anthers was undetectable specifically in the tapetum and in pollen, which could be associated with male sterility. We also generated RNAi-cPGI knockdown lines with strong cPGI repression in floral buds that exhibited reduced male fertility due to the degeneration of most pollen. Histological analyses indicated that the synthesis of intersporal callose walls was impaired, causing microsporocytes to fail to separate haploid daughter nuclei to form tetrads, which might be responsible for subsequent pollen degeneration. We successfully isolated cpgi knockout mutants in the progeny of a heterozygous cpgi mutant floral-dipped with sugar solutions. The rescued cpgi mutants exhibited diminished young vegetative growth, reduced female fertility, and impaired intersporal callose wall formation in a meiocyte, and, thus, male sterility. Collectively, our data suggest that cPGI plays a vital role in carbohydrate partitioning, which is indispensable for microsporogenesis and early embryogenesis.

Nutrient-Limited Operational Strategies for the Microbial Production of Biochemicals

Limiting an essential nutrient has a profound impact on microbial growth. The notion of growth under limited conditions was first described using simple Monod kinetics proposed in the 1940s. Different operational modes (chemostat, fed-batch processes) were soon developed to address questions related to microbial physiology and cell maintenance and to enhance product formation. With more recent developments of metabolic engineering and systems biology, as well as high-throughput approaches, the focus of current engineers and applied microbiologists has shifted from these fundamental biochemical processes. This review draws attention again to nutrient-limited processes. Indeed, the sophisticated gene editing tools not available to pioneers offer the prospect of metabolic engineering strategies which leverage nutrient limited processes. Thus, nutrient- limited processes continue to be very relevant to generate microbially derived biochemicals.

Central Metabolism and Growth Rate Impacts on Hydrogen and Carbon Isotope Fractionation During Amino Acid Synthesis in E. coli

Compound specific stable isotope analysis (CSIA) of amino acids from bacterial biomass is a newly emerging powerful tool for exploring central carbon metabolism pathways and fluxes. By comparing isotopic values and fractionations relative to water and growth substrate, the impact of variable flow path for metabolites through different central metabolic pathways, perturbations of these paths, and their resultant consequences on intracellular pools and resultant biomass may be elucidated. Here, we explore the effects that central carbon metabolism and growth rate can have on stable hydrogen (δ2H) and carbon (δ13C) compound specific isotopic values of amino acids, and whether diagnostic isotopic fingerprints are revealed by these paired analyses. We measured δ2H and δ13C in amino acids in the wild type Escherichia coli (MG1655) across a range of growth rates in chemostat cultures to address the unknown isotopic consequences as metabolic fluxes are shuffled between catabolic and anabolic metabolisms. Additionally, two E. coli knockout mutants, one with deficiency in glycolysis -pgi (LC1888) and another inhibiting the oxidative pentose phosphate pathway (OPPP) -zwf (LC1889), were grown on glucose and used as a comparison against the wild type E. coli (MG1655) to address the isotopic changes of amino acids produced in these perturbed metabolic pathways. Amino acid δ2H values, which collectively vary in composition by more than 400‰, are altered along with δ13C values demonstrating fundamental shifts in central metabolic pathways and/or fluxes. Within our linear discriminant analysis with a simple model organism to examine potential amino acid fingerprinting, our knockout strains and variable growth rate samples plot across a wider array of organism classification than merely within the boundaries of other bacterial data.

Exploring the Glucose Fluxotype of the E. coli y-ome Using High-Resolution Fluxomics

We have developed a robust workflow to measure high-resolution fluxotypes (metabolic flux phenotypes) for large strain libraries under fully controlled growth conditions. This was achieved by optimizing and automating the whole high-throughput fluxomics process and integrating all relevant software tools. This workflow allowed us to obtain highly detailed maps of carbon fluxes in the central carbon metabolism in a fully automated manner. It was applied to investigate the glucose fluxotypes of 180 Escherichia coli strains deleted for y-genes. Since the products of these y-genes potentially play a role in a variety of metabolic processes, the experiments were designed to be agnostic as to their potential metabolic impact. The obtained data highlight the robustness of E. coli’s central metabolism to y-gene deletion. For two y-genes, deletion resulted in significant changes in carbon and energy fluxes, demonstrating the involvement of the corresponding y-gene products in metabolic function or regulation. This work also introduces novel metrics to measure the actual scope and quality of high-throughput fluxomics investigations.

Evolutionary dynamics and structural consequences of de novo beneficial mutations and mutant lineages arising in a constant environment

BACKGROUND

Microbial evolution experiments can be used to study the tempo and dynamics of evolutionary change in asexual populations, founded from single clones and growing into large populations with multiple clonal lineages. High-throughput sequencing can be used to catalog de novo mutations as potential targets of selection, determine in which lineages they arise, and track the fates of those lineages. Here, we describe a long-term experimental evolution study to identify targets of selection and to determine when, where, and how often those targets are hit.

RESULTS

We experimentally evolved replicate Escherichia coli populations that originated from a mutator/nonsense suppressor ancestor under glucose limitation for between 300 and 500 generations. Whole-genome, whole-population sequencing enabled us to catalog 3346 de novo mutations that reached > 1% frequency. We sequenced the genomes of 96 clones from each population when allelic diversity was greatest in order to establish whether mutations were in the same or different lineages and to depict lineage dynamics. Operon-specific mutations that enhance glucose uptake were the first to rise to high frequency, followed by global regulatory mutations. Mutations related to energy conservation, membrane biogenesis, and mitigating the impact of nonsense mutations, both ancestral and derived, arose later. New alleles were confined to relatively few loci, with many instances of identical mutations arising independently in multiple lineages, among and within replicate populations. However, most never exceeded 10% in frequency and were at a lower frequency at the end of the experiment than at their maxima, indicating clonal interference. Many alleles mapped to key structures within the proteins that they mutated, providing insight into their functional consequences.

CONCLUSIONS

Overall, we find that when mutational input is increased by an ancestral defect in DNA repair, the spectrum of high-frequency beneficial mutations in a simple, constant resource-limited environment is narrow, resulting in extreme parallelism where many adaptive mutations arise but few ever go to fixation.

Metabolic flux configuration determination using information entropy

Constraint-based models use steady-state mass balances to define a solution space of flux configurations, which can be narrowed down by measuring as many fluxes as possible. Due to loops and redundant pathways, this process typically yields multiple alternative solutions. To address this ambiguity, flux sampling can estimate the probability distribution of each flux, or a flux configuration can be singled out by further minimizing the sum of fluxes according to the assumption that cellular metabolism favors states where enzyme-related costs are economized. However, flux sampling is susceptible to artifacts introduced by thermodynamically infeasible cycles and is it not clear if the economy of fluxes assumption (EFA) is universally valid. Here, we formulated a constraint-based approach, MaxEnt, based on the principle of maximum entropy, which in this context states that if more than one flux configuration is consistent with a set of experimentally measured fluxes, then the one with the minimum amount of unwarranted assumptions corresponds to the best estimation of the non-observed fluxes. We compared MaxEnt predictions to Escherichia coli and Saccharomyces cerevisiae publicly available flux data. We found that the mean square error (MSE) between experimental and predicted fluxes by MaxEnt and EFA-based methods are three orders of magnitude lower than the median of 1,350,000 MSE values obtained using flux sampling. However, only MaxEnt and flux sampling correctly predicted flux through E. coli’s glyoxylate cycle, whereas EFA-based methods, in general, predict no flux cycles. We also tested MaxEnt predictions at increasing levels of overflow metabolism. We found that MaxEnt accuracy is not affected by overflow metabolism levels, whereas the EFA-based methods show a decreasing performance. These results suggest that MaxEnt is less sensitive than flux sampling to artifacts introduced by thermodynamically infeasible cycles and that its predictions are less susceptible to overfitting than EFA-based methods.

Artificial consortium demonstrates emergent properties of enhanced cellulosic-sugar degradation and biofuel synthesis

Planktonic cultures, of a rationally designed consortium, demonstrated emergent properties that exceeded the sums of monoculture properties, including a >200% increase in cellobiose catabolism, a >100% increase in glycerol catabolism, a >800% increase in ethanol production, and a >120% increase in biomass productivity. The consortium was designed to have a primary and secondary-resource specialist that used crossfeeding with a positive feedback mechanism, division of labor, and nutrient and energy transfer via necromass catabolism. The primary resource specialist was Clostridium phytofermentans (a.k.a. Lachnoclostridium phytofermentans), a cellulolytic, obligate anaerobe. The secondary-resource specialist was Escherichia coli, a versatile, facultative anaerobe, which can ferment glycerol and byproducts of cellobiose catabolism. The consortium also demonstrated emergent properties of enhanced biomass accumulation when grown as biofilms, which created high cell density communities with gradients of species along the vertical axis. Consortium biofilms were robust to oxic perturbations with E. coli consuming O₂, creating an anoxic environment for C. phytofermentans. Anoxic/oxic cycling further enhanced biomass productivity of the biofilm consortium, increasing biomass accumulation ~250% over the sum of the monoculture biofilms. Consortium emergent properties were credited to several synergistic mechanisms. E. coli consumed inhibitory byproducts from cellobiose catabolism, driving higher C. phytofermentans growth and higher cellulolytic enzyme production, which in turn provided more substrate for E. coli. E. coli necromass enhanced C. phytofermentans growth while C. phytofermentans necromass aided E. coli growth via the release of peptides and amino acids, respectively. In aggregate, temporal cycling of necromass constituents increased flux of cellulose-derived resources through the consortium. The study establishes a consortia-based, bioprocessing strategy built on naturally occurring interactions for improved conversion of cellulose-derived sugars into bioproducts.

Pristane promotes anaerobic glycolysis to facilitate proinflammatory activation of macrophages and development of arthritis

Pristane-induced arthritis (PIA) could be adoptively transferred by splenic T cells in rats, and innate immunity should play critical roles in T cell activation. However, in pre-clinical stage, the activation mechanism of innate cells like macrophages remains unclear. Here we found that PIA was dependent on macrophages since cell depletion alleviated disease severity. Splenic macrophages of PIA rats showed M1 phenotypic shifting. The quantitative proteomics analysis suggested that macrophages initiated metabolic reprogramming with the conversion of aerobic oxidation to glycolysis in response to pristane in vivo. Notably, macrophages treated with pristane showed mitochondrial dysregulation and increased glycolysis flux and enzyme activity. Additionally, TNFα production, strongly associating with the glycolysis enzyme Ldha/Ldhb, could be reduced as glycolysis was inhibited or be enhanced as citrate cycle was blocked. This work provides detailed insights into the molecular mechanisms of pristane-mediated metabolic reprogramming in macrophages and suggests a new therapeutic strategy for arthritic disorders.

An oxidoreductase gene ZMO1116 enhances the p-benzoquinone biodegradation and chiral lactic acid fermentability of Pediococcus acidilactici

p-Benzoquinone (BQ) is a lignin-derived inhibitor to microbial strains. Unlike the furan inhibitors, p-benzoquinone is recalcitrant to traditional detoxification methods. This study shows a biological degradation of p-benzoquinone and a simultaneous D-lactic acid fermentation by an engineered Pediococcus acidilactici strain. The overexpression of an oxidoreductase gene ZMO1116 from Zymomonas mobilis encoding oxidoreductase was identified to improve the D-lactic acid fermentability of P. acidilactici against p-benzoquinone. The gene ZMO1116 was integrated into the genome of P. acidilactici and enabled the engineered P. acidilactici to convert p-benzoquinone into less toxic hydroquinone (HQ), resulting in the improved p-benzoquinone tolerance. Simultaneous saccharification and co-fermentation (SSCF) was conducted using the pretreated and biodetoxified corn stover containing p-benzoquinone, the D-lactic acid production of the engineered strain (123.8 g/L) was 21.4 % higher than the parental strain (102.0 g/L). This study provides a practical method on robust p-benzoquinone tolerance and efficient cellulosic chiral lactic acid fermentation from lignocellulose feedstock.

Control of the galactose-to-glucose consumption ratio in co-fermentation using engineered Escherichia coli strains

Marine biomasses capable of fixing carbon dioxide have attracted attention as an alternative to fossil resources for fuel and chemical production. Although a simple co-fermentation of fermentable sugars, such as glucose and galactose, has been reported from marine biomass, no previous report has discussed the fine-control of the galactose-to-glucose consumption ratio in this context. Here, we sought to finely control the galactose-to-glucose consumption ratio in the co-fermentation of these sugars using engineered Escherichia coli strains. Toward this end, we constructed E. coli strains GR2, GR2P, and GR2PZ by knocking out galRS, galRS-pfkA, and galRS-pfkA-zwf, respectively, in parent strain W3110. We found that strains W3110, GR2, GR2P, and GR2PZ achieved 0.03, 0.09, 0.12, and 0.17 galactose-to-glucose consumption ratio (specific galactose consumption rate per specific glucose consumption rate), respectively, during co-fermentation. The ratio was further extended to 0.67 by integration of a brief process optimization for initial sugar ratio using GR2P strain. The strategy reported in this study will be helpful to expand our knowledge on the galactose utilization under glucose conditions.

Escherichia coli metabolism under short-term repetitive substrate dynamics: adaptation and trade-offs

Background

Microbial metabolism is highly dependent on the environmental conditions. Especially, the substrate concentration, as well as oxygen availability, determine the metabolic rates. In large-scale bioreactors, microorganisms encounter dynamic conditions in substrate and oxygen availability (mixing limitations), which influence their metabolism and subsequently their physiology. Earlier, single substrate pulse experiments were not able to explain the observed physiological changes generated under large-scale industrial fermentation conditions.

Results

In this study we applied a repetitive feast–famine regime in an aerobic Escherichia coli culture in a time-scale of seconds. The regime was applied for several generations, allowing cells to adapt to the (repetitive) dynamic environment. The observed response was highly reproducible over the cycles, indicating that cells were indeed fully adapted to the regime. We observed an increase of the specific substrate and oxygen consumption (average) rates during the feast–famine regime, compared to a steady-state (chemostat) reference environment. The increased rates at same (average) growth rate led to a reduced biomass yield (30% lower). Interestingly, this drop was not followed by increased by-product formation, pointing to the existence of energy-spilling reactions. During the feast–famine cycle, the cells rapidly increased their uptake rate. Within 10 s after the beginning of the feeding, the substrate uptake rate was higher (4.68 μmol/g_CDW/s) than reported during batch growth (3.3 μmol/g_CDW/s). The high uptake led to an accumulation of several intracellular metabolites, during the feast phase, accounting for up to 34% of the carbon supplied. Although the metabolite concentrations changed rapidly, the cellular energy charge remained unaffected, suggesting well-controlled balance between ATP producing and ATP consuming reactions.

Conclusions

The adaptation of the physiology and metabolism of E. coli under substrate dynamics, representative for large-scale fermenters, revealed the existence of several cellular mechanisms coping with stress. Changes in the substrate uptake system, storage potential and energy-spilling processes resulted to be of great importance. These metabolic strategies consist a meaningful step to further tackle reduced microbial performance, observed under large-scale cultivations.

Effect of dissolved oxygen on L-methionine production from glycerol by Escherichia coli W3110BL using metabolic flux analysis method

L-Methionine is an essential amino acid in humans, which plays an important role in the synthesis of some important amino acids and proteins. In this work, metabolic flux of batch fermentation of L-methionine with recombinant Escherichia coli W3110BL was analyzed using the flux balance analysis method, which estimated the intracellular flux distributions under different dissolved oxygen conditions. The results revealed the producing L-methionine flux of 4.8 mmol/(g cell·h) [based on the glycerol uptake flux of 100 mmol/(g cell·h)] was obtained at 30% dissolved oxygen level which was higher than that of other dissolved oxygen levels. The carbon fluxes for synthesizing L-methionine were mainly obtained from the pathway of phosphoenolpyruvate to oxaloacetic acid [15.6 mmol/(g cell·h)] but not from the TCA cycle. Hence, increasing the flow from phosphoenolpyruvate to oxaloacetic acid by enhancing the enzyme activity of phosphoenolpyruvate carboxylase might be conducive to the production of L-methionine. Additionally, pentose phosphate pathway could provide a large amount of reducing power NADPH for the synthesis of amino acids and the flux could increase from 41 mmol/(g cell·h) to 51 mmol/(g cell·h) when changing the dissolved oxygen levels, thus meeting the requirement of NADPH for L-methionine production and biomass synthesis. Therefore, the following modification of the strains should based on the improvement of the key pathway and the NAD(P)/NAD(P)H metabolism.

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.

Lessons from Two Design-Build-Test-Learn Cycles of Dodecanol Production in Escherichia coli Aided by Machine Learning

The Design-Build-Test-Learn (DBTL) cycle, facilitated by exponentially improving capabilities in synthetic biology, is an increasingly adopted metabolic engineering framework that represents a more systematic and efficient approach to strain development than historical efforts in biofuels and biobased products. Here, we report on implementation of two DBTL cycles to optimize 1-dodecanol production from glucose using 60 engineered Escherichia coli MG1655 strains. The first DBTL cycle employed a simple strategy to learn efficiently from a relatively small number of strains (36), wherein only the choice of ribosome-binding sites and an acyl-ACP/acyl-CoA reductase were modulated in a single pathway operon including genes encoding a thioesterase (UcFatB1), an acyl-ACP/acyl-CoA reductase (Maqu_2507, Maqu_2220, or Acr1), and an acyl-CoA synthetase (FadD). Measured variables included concentrations of dodecanol and all proteins in the engineered pathway. We used the data produced in the first DBTL cycle to train several machine-learning algorithms and to suggest protein profiles for the second DBTL cycle that would increase production. These strategies resulted in a 21% increase in dodecanol titer in Cycle 2 (up to 0.83 g/L, which is more than 6-fold greater than previously reported batch values for minimal medium). Beyond specific lessons learned about optimizing dodecanol titer in E. coli, this study had findings of broader relevance across synthetic biology applications, such as the importance of sequencing checks on plasmids in production strains as well as in cloning strains, and the critical need for more accurate protein expression predictive tools.

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.

Dark fermentation metabolic models to study strategies for hydrogen consumers inhibition

A Flux Balance Analysis (FBA) metabolic model of dark fermentation was developed for anaerobic mixed cultures. In particular, the model was applied to evaluate the effect of a specific inoculum pre-treatment strategy, addition of waste frying oil (WFO) on H₂-producing and H₂-consuming metabolic pathways. Productions of volatile fatty acid (VFAs), CO₂, H₂ and CH₄ measured through triplicate batch experiments, were used as constraints for the FBA model, to compute fluxes trough different metabolic pathways. FBA model could estimate the effect of pre-treatment with WFO on major microbial populations present in the mixed community (H₂ producing bacteria, homoacetogen and methanogens). Results revealed that low concentrations of WFO did not completely inhibited hydrogenotrophic methanogens. FBA showed that acetoclastic methanogens were more sensitive to WFO, in comparison to hydrogenotrophic methanogens. The proposed model can be used to study H₂ production by any other mixed microbial culture with similar substrates.

Proteomic and metabolomic analysis of the cellular biomarkers related to inhibitors tolerance in Zymomonas mobilis ZM4

BACKGROUND

Toxic compounds present in both the hydrolysate and pyrolysate of lignocellulosic biomass severely hinder the further conversion of lignocellulose-derived fermentable sugars into useful chemicals by common biocatalysts like Zymomonas mobilis, which has remarkable advantages over yeast. Although the extra detoxification treatment prior to fermentation process can help biocatalysts to eliminate the inhibitory environment, it is not environment friendly and cost effective for industrial application. As also reported by previous studies, an ideal and holistic approach to solve this issue is to develop microbial strains with inhibitor tolerance. However, previously engineered strains had the limitation that they could not cope well with the synergistic effect of multiple inhibitors as they are resistant only to a single inhibitor. Hence, understanding the universal cellular responses of Z. mobilis to various inhibitors may guide the designing of rational strategies to obtain more robust engineered strains for biofuel production from lignocellulosic biomass.

RESULTS

Quantitative proteomics and metabolomics approaches were used to determine the cellular responses of Z. mobilis ZM4 to representative biomass-derived inhibitors like formic acid, acetic acid, furfural, 5-hydroxymethylfurfural, and phenol. The differentially expressed proteins identified under the challenge of single and combined inhibitors were involved in cell wall/membrane biogenesis, energy production, DNA replication, DNA recombination, DNA repair, DNA transcription, RNA translation, posttranslational modification, biosynthesis of amino acids, central carbon metabolism, etc. Metabolomics analysis showed that the up- or down-regulation pattern of metabolites was changed consistently with that of relevant proteins.

CONCLUSION

Fifteen up-regulated proteins (e.g., Isopropylmalate isomerase LeuC, transcription-repair-coupling factor Mfd, and phosphoglucose isomerase PGI) and thirteen down-regulated proteins (e.g., TonB-dependent transporter ZMO1522, transcription termination factor Rho, and S1/P1 nuclease ZMO0127) were identified as candidate proteins related to all the stress conditions, implying that these proteins are potential biomarkers for the improvement of Z. mobilis ZM4 to resist complex biomass-derived inhibitors. These data can be used to generate a database of inhibitor-tolerance biomarkers, which could provide a basis for engineering Z. mobilis that would be able to grow in the presence of multiple inhibitors and directly ferment the biomass-derived sugars into biofuels.

Engineering an in vivo EP-bifido pathway in Escherichia coli for high-yield acetyl-CoA generation with low CO2 emission

The low carbon yield from native metabolic machinery produces unfavorable process economics during the biological conversion of substrates to desirable bioproducts. To obtain higher carbon yields, we constructed a carbon conservation pathway named EP-bifido pathway in Escherichia coli by combining Embden-Meyerhof-Parnas Pathway, Pentose Phosphate Pathway and "bifid shunt", to generate high yield acetyl-CoA from glucose. ¹³C-Metabolic flux analysis confirmed the successful and appropriate employment of the EP-bifido pathway. The CO₂ release during fermentation significantly reduced compared with the control strains. Then we demonstrated the in vivo effectiveness of the EP-bifido pathway using poly-β-hydroxybutyrate (PHB), mevalonate and fatty acids as example products. The engineered EP-bifido strains showed greatly improved PHB yield (from 26.0 mol% to 63.7 mol%), fatty acid yield (from 9.17% to 14.36%), and the highest mevalonate yield yet reported (64.3 mol% without considering the substrates used for cell mass formation). The synthetic pathway can be employed in the production of chemicals that use acetyl-CoA as a precursor and can be extended to other microorganisms.

Re-engineering of an Escherichia coli K-12 strain for the efficient production of recombinant human Interferon Gamma

The Escherichia coli phosphoglucose isomerase (pgi) mutant strain GALG20 was developed previously from wild-type K12 strain MG1655 for increased plasmid yield. To investigate the potential effects of the pgi deletion/higher plasmid levels on recombinant human Interferon Gamma (IFN-γ) production, a detailed network of the central metabolic pathway (100 metabolites, 114 reactions) of GALG20 and MG1655 was constructed. Elementary mode analysis (EMA) was then performed to compare the phenotypic spaces of both the strains and to check the effect of the pgi deletion on flux efficiency of each metabolic reaction. The results suggested that pgi deletion increases amino acid biosynthesis and flux efficiency towards IFN-γ synthesis by 11%. To further confirm the qualitative prediction that the pgi mutation favours recombinant human IFN-γ expression, GALG20 and MG1655 were lysogenised, transformed with a plasmid coding for IFN-γ and tested alongside with BL21(DE3) for their expression capabilities in shake flask experiments using complex media. IFN-γ gene expression was analysed by quantifying plasmid and mRNA copy number per cell and IFN-γ protein production level. Specific IFN-γ yields confirmed the in silico metabolic network predictions, with GALG20(DE3) producing 3.0-fold and 1.5-fold more IFN-γ as compared to MG1655(DE3) and BL21(DE3), respectively. Most of the total IFN-γ was expressed as inclusion bodies across the three strains: 95% in GALG20(DE3), 97% in BL21(DE3) and 72% in MG1655(DE3). The copy number of mRNA coding for IFN-γ was found to be higher in GALG20(DE3) as compared to the other two strains. Overall, these findings show that GALG20(DE3) has the potential to become an excellent protein expression strain.

Engineering E. coli for simultaneous glucose-xylose utilization during methyl ketone production

Background

We previously developed an E. coli strain that overproduces medium-chain methyl ketones for potential use as diesel fuel blending agents or as flavors and fragrances. To date, the strain’s performance has been optimized during growth with glucose. However, lignocellulosic biomass hydrolysates also contain a substantial portion of hemicellulose-derived xylose, which is typically the second most abundant sugar after glucose. Commercialization of the methyl ketone-producing technology would benefit from the increased efficiency resulting from simultaneous, rather than the native sequential (diauxic), utilization of glucose and xylose.

Results

In this study, genetic manipulations were performed to alleviate carbon catabolite repression in our most efficient methyl ketone-producing strain. A strain engineered for constitutive expression of xylF and xylA (involved in xylose transport and metabolism) showed synchronized glucose and xylose consumption rates. However, this newly acquired capability came at the expense of methyl ketone titer, which decreased fivefold. Further efforts were made to improve methyl ketone production in this strain, and we found that two strategies were effective at enhancing methyl ketone titer: (1) chromosomal deletion of pgi (glucose-6-phosphate isomerase) to increase intracellular NADPH supply and (2) downregulation of CRP (cAMP receptor protein) expression by replacement of the native RBS with an RBS chosen based upon mutant library screening results. Combining these strategies resulted in the most favorable overall phenotypes for simultaneous glucose–xylose consumption without compromising methyl ketone titer at both 1 and 2% total sugar concentrations in shake flasks.

Conclusions

This work demonstrated a strategy for engineering simultaneous utilization of C₆ and C₅ sugars in E. coli without sacrificing production of fatty acid-derived compounds.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-0862-6) contains supplementary material, which is available to authorized users.

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 high-resolution ¹³C-metabolic flux analysis were performed on 10 adaptively evolved pgi knockouts of Escherichia coliPgi 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.

Expression of heterologous non-oxidative pentose phosphate pathway from Bacillus methanolicus and phosphoglucose isomerase deletion improves methanol assimilation and metabolite production by a synthetic Escherichia coli methylotroph

Synthetic methylotrophy aims to develop non-native methylotrophic microorganisms to utilize methane or methanol to produce chemicals and biofuels. We report two complimentary strategies to further engineer a previously engineered methylotrophic E. coli strain for improved methanol utilization. First, we demonstrate improved methanol assimilation in the presence of small amounts of yeast extract by expressing the non-oxidative pentose phosphate pathway (PPP) from Bacillus methanolicus. Second, we demonstrate improved co-utilization of methanol and glucose by deleting the phosphoglucose isomerase gene (pgi), which rerouted glucose carbon flux through the oxidative PPP. Both strategies led to significant improvements in methanol assimilation as determined by ¹³C-labeling in intracellular metabolites. Introduction of an acetone-formation pathway in the pgi-deficient methylotrophic E. coli strain led to improved methanol utilization and acetone titers during glucose fed-batch fermentation.

The genetic basis for the adaptation of E. coli to sugar synthesis from CO2

Understanding the evolution of a new metabolic capability in full mechanistic detail is challenging, as causative mutations may be masked by non-essential "hitchhiking" mutations accumulated during the evolutionary trajectory. We have previously used adaptive laboratory evolution of a rationally engineered ancestor to generate an Escherichia coli strain able to utilize CO₂ fixation for sugar synthesis. Here, we reveal the genetic basis underlying this metabolic transition. Five mutations are sufficient to enable robust growth when a non-native Calvin-Benson-Bassham cycle provides all the sugar-derived metabolic building blocks. These mutations are found either in enzymes that affect the efflux of intermediates from the autocatalytic CO₂ fixation cycle toward biomass (prs, serA, and pgi), or in key regulators of carbon metabolism (crp and ppsR). Using suppressor analysis, we show that a decrease in catalytic capacity is a common feature of all mutations found in enzymes. These findings highlight the enzymatic constraints that are essential to the metabolic stability of autocatalytic cycles and are relevant to future efforts in constructing non-native carbon fixation pathways.

Effect of NADPH availability on free fatty acid production in Escherichia coli

Microbial conversion of renewable carbon sources to free fatty acids has attracted significant attention in recent years. Accumulation of free fatty acids in Escherichia coli by overexpression of an acyl-ACP thioesterase which can break the fatty acid elongation has been well established. Various efforts have been made to increase fatty acid production in E. coli by enhancing the enzymes involved in the fatty acid synthesis cycle or host strain manipulations. The current study focused on the effect of NADPH availability on free fatty acids (FFAs) productivity. There are two reduction steps in the fatty acid elongation cycle which are catalyzed by beta keto-ACP reductase (FabG) and enoyl-ACP reductase (FabI), respectively. It is reported that FabI can use either NADH or NADPH as cofactor, while FabG only uses NADPH in E. coli. Fatty acid production dropped dramatically in the glucose-6-phosphate dehydrogenase (encoded by the zwf gene) deficient strain. Similarly, the pntB (which encodes one of the subunit of proton-translocating membrane bounded transhydrogenase PntAB) and udhA (which encodes the energy dependent cytoplasmic transhydrogenase UdhA) double mutant strain also showed an 88.8% decrease in free fatty acid production. Overexpression of PntAB and NadK restored the fatty acid production capability of these two mutant strains. These results indicated that the availability of NADPH played a very important role in fatty acid production.

Identification of GntR as regulator of the glucose metabolism in Pseudomonas aeruginosa

In contrast to Escherichia coli, glucose metabolism in pseudomonads occurs exclusively through the Entner-Doudoroff (ED) pathway. This pathway, as well as the three routes to generate the initial ED pathway substrate, 6-phosphogluconate, is regulated by the PtxS, HexR and GtrS/GltR systems. With GntR (PA2320) we report here the identification of an additional regulator in Pseudomonas aeruginosa PAO1. GntR repressed its own expression as well as that of the GntP gluconate permease. In contrast to PtxS and GtrS/GltR, GntR did not modulate expression of the toxA gene encoding the exotoxin A virulence factor. GntR was found to bind to promoters PgntR and PgntP and the consensus sequence of its operator was defined as 5'-AC-N-AAG-N-TAGCGCT-3'. Both operator sites overlapped with the RNA polymerase binding site and we show that GntR employs an effector mediated de-repression mechanism. The release of promoter bound GntR is induced by gluconate and 6-phosphogluconate that bind with similar apparent affinities to the GntR/DNA complex. GntR and PtxS are paralogous and may have evolved from a common ancestor. The concerted action of four regulatory systems in the regulation of glucose metabolism in Pseudomonas can be considered as a model to understand complex regulatory circuits in bacteria.

Co-utilization of hexoses by a microconsortium of sugar-specific E. coli strains

Escherichia coli is an important commercial species used for production of biofuels, biopolymers, organic acids, sugar alcohols, and natural compounds. Processed biomass and agroindustrial byproducts serve as low-cost nutrient sources and contain a variety of hexoses available for bioconversion. However, metabolism of hexose mixtures by E. coli is inefficient due to carbon catabolite repression (CCR), where the transport and catabolic activity of one or more carbon sources is repressed and/or inhibited by the transport and catabolism of another carbon source. In this work, we developed a microconsortium of different E. coli strains, each engineered to preferentially catabolize a different hexose-glucose, galactose, or mannose. We modified the specificity and preference of carbon source using a combination of rational strain design and adaptive evolution. The modifications ultimately resulted in strains that preferentially catabolized their specified sugar. Finally, comparative analysis in galactose- and mannose-rich sugar mixtures revealed that the consortium grew faster and to higher cell densities compared to the wild-type strain. Biotechnol. Bioeng. 2017;114: 2309-2318. © 2017 Wiley Periodicals, Inc.

Eliminating acetate formation improves citramalate production by metabolically engineered Escherichia coli

BACKGROUND

Citramalate, a chemical precursor to the industrially important methacrylic acid (MAA), can be synthesized using Escherichia coli overexpressing citramalate synthase (cimA gene). Deletion of gltA encoding citrate synthase and leuC encoding 3-isopropylmalate dehydratase were critical to achieving high citramalate yields. Acetate is an undesirable by-product potentially formed from pyruvate and acetyl-CoA, the precursors of citramalate during aerobic growth of E. coli. This study investigated strategies to minimize acetate and maximize citramalate production in E. coli mutants expressing the cimA gene.

RESULTS

Key knockouts that minimized acetate formation included acetate kinase (ackA), phosphotransacetylase (pta), and in particular pyruvate oxidase (poxB). Deletion of glucose 6-phosphate dehydrogenase (zwf) and ATP synthase (atpFH) aimed at improving glycolytic flux negatively impacted cell growth and citramalate accumulation in shake flasks. In a repetitive fed-batch process, E. coli gltA leuC ackA-pta poxB overexpressing cimA generated 54.1 g/L citramalate with a yield of 0.64 g/g glucose (78% of theoretical maximum yield), and only 1.4 g/L acetate in 87 h.

CONCLUSIONS

This study identified the gene deletions critical to reducing acetate accumulation during aerobic growth and citramalate production in metabolically engineered E. coli strains. The citramalate yield and final titer relative to acetate at the end of the fed-batch process are the highest reported to date (a mass ratio of citramalate to acetate of nearly 40) without being detrimental to citramalate productivity, significantly improving a potential process for the production of this five-carbon chemical.

Metabolic regulation is sufficient for global and robust coordination of glucose uptake, catabolism, energy production and growth in Escherichia coli

The metabolism of microorganisms is regulated through two main mechanisms: changes of enzyme capacities as a consequence of gene expression modulation (“hierarchical control”) and changes of enzyme activities through metabolite-enzyme interactions. An increasing body of evidence indicates that hierarchical control is insufficient to explain metabolic behaviors, but the system-wide impact of metabolic regulation remains largely uncharacterized. To clarify its role, we developed and validated a detailed kinetic model of Escherichia coli central metabolism that links growth to environment. Metabolic control analyses confirm that the control is widely distributed across the network and highlight strong interconnections between all the pathways. Exploration of the model solution space reveals that several robust properties emerge from metabolic regulation, from the molecular level (e.g. homeostasis of total metabolite pool) to the overall cellular physiology (e.g. coordination of carbon uptake, catabolism, energy and redox production, and growth), while allowing a large degree of flexibility at most individual metabolic steps. These properties have important physiological implications for E. coli and significantly expand the self-regulating capacities of its metabolism.

Metabolism is a fundamental biochemical process that enables cells to operate and grow by converting nutrients into ‘building blocks’ and energy. Metabolism happens through the work of enzymes, which are encoded by genes. Thus, genes and their regulation are often thought of controlling metabolism, somewhat at the top of a hierarchical control system. However, an increasing body of evidence indicates that metabolism plays an active role in the control of its own operation via a dense network of metabolite-enzyme interactions. The system-wide role of metabolic regulation is hard to dissect and so far remains largely uncharacterized. To better understand its role, we constructed a detailed kinetic model of the carbon and energy metabolism of the bacterium Escherichia coli, a model organism in Systems and Synthetic biology. Model simulations indicate that kinetic considerations of metabolism alone can explain data from hundreds of experiments, without needing to invoke regulation of gene expression. In particular, metabolic regulation is sufficient to coordinate carbon utilization, redox and energy production, and growth, while maintaining local flexibility at individual metabolic steps. These findings indicate that the self-regulating capacities of E. coli metabolism are far more significant than previously expected, and improve our understanding on how cells work.

Metabolic engineering of Escherichia coli for microbial production of L-methionine

L-methionine has attracted a great deal of attention for its nutritional, pharmaceutical, and clinical applications. In this study, Escherichia coli W3110 was engineered via deletion of a negative transcriptional regulator MetJ and over-expression of homoserine O-succinyltransferase MetA together with efflux transporter YjeH, resulting in L-methionine overproduction which is up to 413.16 mg/L. The partial inactivation of the L-methionine import system MetD via disruption of metI made the engineered E. coli ΔmetJ ΔmetI/pTrcA*H more tolerant to high L-ethionine concentration and accumulated L-methionine to a level 43.65% higher than that of E. coli W3110 ΔmetJ/pTrcA*H. Furthermore, deletion of lysA, which blocks the lysine biosynthesis pathway, led to a further 8.5-fold increase in L-methionine titer of E. coli ΔmetJ ΔmetI ΔlysA/pTrcA*H. Finally, addition of Na₂ S₂ O₃ to the media led to an increase of fermentation titer of 11.45%. After optimization, constructed E. coli ΔmetJ ΔmetI ΔlysA/pTrcA*H was able to produce 9.75 g/L L-methionine with productivity of 0.20 g/L/h in a 5 L bioreactor. This novel metabolically tailored strain of E. coli provides an efficient platform for microbial production of L-methionine. Biotechnol. Bioeng. 2017;114: 843-851. © 2016 Wiley Periodicals, Inc.

SUMOFLUX: A Generalized Method for Targeted 13C Metabolic Flux Ratio Analysis

Metabolic fluxes are a cornerstone of cellular physiology that emerge from a complex interplay of enzymes, carriers, and nutrients. The experimental assessment of in vivo intracellular fluxes using stable isotopic tracers is essential if we are to understand metabolic function and regulation. Flux estimation based on 13C or 2H labeling relies on complex simulation and iterative fitting; processes that necessitate a level of expertise that ordinarily preclude the non-expert user. To overcome this, we have developed SUMOFLUX, a methodology that is broadly applicable to the targeted analysis of 13C-metabolic fluxes. By combining surrogate modeling and machine learning, we trained a predictor to specialize in estimating flux ratios from measurable 13C-data. SUMOFLUX targets specific flux features individually, which makes it fast, user-friendly, applicable to experimental design and robust in terms of experimental noise and exchange flux magnitude. Collectively, we predict that SUMOFLUX's properties realistically pave the way to high-throughput flux analyses.

A simple strategy guides the complex metabolic regulation in Escherichia coli

A way to decipher the complexity of the cellular metabolism is to study the effect of different external perturbations. Through an analysis over a sufficiently large set of gene knockouts and growing conditions, one aims to find a unifying principle that governs the metabolic regulation. For instance, it is known that the cessation of the microorganism proliferation after a gene deletion is only transient. However, we do not know the guiding principle that determines the partial or complete recovery of the growth rate, the corresponding redistribution of the metabolic fluxes and the possible different phenotypes. In spite of this large variety in the observed metabolic adjustments, we show that responses of E. coli to several different perturbations can always be derived from a sequence of greedy and myopic resilencings. This simple mechanism provides a detailed explanation for the experimental dynamics both at cellular (proliferation rate) and molecular level (¹³C-determined fluxes), also in case of appearance of multiple phenotypes. As additional support, we identified an example of a simple network motif that is capable of implementing this myopic greediness in the regulation of the metabolism.

Reinventing the Wheel and Making It Round Again: Evolutionary Convergence in Buchnera-Serratia Symbiotic Consortia between the Distantly Related Lachninae Aphids Tuberolachnus salignus and Cinara cedri

Virtually all aphids (Aphididae) harbor Buchnera aphidicola as an obligate endosymbiont to compensate nutritional deficiencies arising from their phloem diet. Many species within the Lachninae subfamily seem to be consistently associated also with Serratia symbiotica We have previously shown that both Cinara (Cinara) cedri and Cinara (Cupressobium) tujafilina (Lachninae: Eulachnini tribe) have indeed established co-obligate associations with both Buchnera and S. symbiotica However, while Buchnera genomes of both Cinara species are similar, genome degradation differs greatly between the two S. symbiotica strains. To gain insight into the essentiality and degree of integration of S. symbiotica within the Lachninae, we sequenced the genome of both Buchnera and S. symbiotica endosymbionts from the distantly related aphid Tuberolachnus salignus (Lachninae: Tuberolachnini tribe). We found a striking level of similarity between the endosymbiotic system of this aphid and that of C. cedri In both aphid hosts, S. symbiotica possesses a highly reduced genome and is found exclusively intracellularly inside bacteriocytes. Interestingly, T. salignus' endosymbionts present the same tryptophan biosynthetic metabolic complementation as C. cedri's, which is not present in C. tujafilina's. Moreover, we corroborate the riboflavin-biosynthetic-role take-over/rescue by S. symbiotica in T. salignus, and therefore, provide further evidence for the previously proposed establishment of a secondary co-obligate endosymbiont in the common ancestor of the Lachninae aphids. Finally, we propose that the putative convergent split of the tryptophan biosynthetic role between Buchnera and S. symbiotica could be behind the establishment of S. symbiotica as an obligate intracellular symbiont and the triggering of further genome degradation.

Web application for genetic modification flux with database to estimate metabolic fluxes of genetic mutants

Computational analysis of metabolic fluxes is essential in understanding the structure and function of a metabolic network and in rationally designing genetically modified mutants for an engineering purpose. We had presented the genetic modification flux (GMF) that predicts the flux distribution of a broad range of genetically modified mutants. To enhance the feasibility and usability of GMF, we have developed a web application with a metabolic network database to predict a flux distribution of genetically modified mutants. One hundred and twelve data sets of Escherichia coli, Corynebacterium glutamicum, Saccharomyces cerevisiae, and Chinese hamster ovary were registered as standard models.

Systematic engineering of the central metabolism in Escherichia coli for effective production of n-butanol

BACKGROUND

Microbes have been extensively explored for production of environment-friendly fuels and chemicals. The microbial fermentation pathways leading to these commodities usually involve many redox reactions. This makes the fermentative production of highly reduced products challenging, because there is a limited NADH output from glucose catabolism. Microbial production of n-butanol apparently represents one typical example.

RESULTS

In this study, we addressed the issue by adjustment of the intracellular redox state in Escherichia coli. This was initiated with strain BuT-8 which carries the clostridial CoA-dependent synthetic pathway. Three metabolite nodes in the central metabolism of the strain were targeted for engineering. First, the pyruvate node was manipulated by enhancement of pyruvate decarboxylation in the oxidative pathway. Subsequently, the pentose phosphate (PP) pathway was amplified at the glucose-6-phosphate (G6P) node. The pathway for G6P isomerization was further blocked to force the glycolytic flux through the PP pathway. It resulted in a growth defect, and the cell growth was later recovered by limiting the tricarboxylic acid cycle at the acetyl-CoA node. Finally, the resulting strain exhibited a high NADH level and enabled production of 6.1 g/L n-butanol with a yield of 0.31 g/g-glucose and a productivity of 0.21 g/L/h.

CONCLUSIONS

The production efficiency of fermentative products in microbes strongly depends on the intracellular redox state. This work illustrates the flexibility of pyruvate, G6P, and acetyl-CoA nodes at the junction of the central metabolism for engineering. In principle, high production of reduced products of interest can be achieved by individual or coordinated modulation of these metabolite nodes.

Reconstruction and Use of Microbial Metabolic Networks: the Core Escherichia coli Metabolic Model as an Educational Guide

Biochemical network reconstructions have become popular tools in systems biology. Metabolicnetwork reconstructions are biochemically, genetically, and genomically (BiGG) structured databases of biochemical reactions and metabolites. They contain information such as exact reaction stoichiometry, reaction reversibility, and the relationships between genes, proteins, and reactions. Network reconstructions have been used extensively to study the phenotypic behavior of wild-type and mutant stains under a variety of conditions, linking genotypes with phenotypes. Such phenotypic simulations have allowed for the prediction of growth after genetic manipulations, prediction of growth phenotypes after adaptive evolution, and prediction of essential genes. Additionally, because network reconstructions are organism specific, they can be used to understand differences between organisms of species in a functional context.There are different types of reconstructions representing various types of biological networks (metabolic, regulatory, transcription/translation). This chapter serves as an introduction to metabolic and regulatory network reconstructions and models and gives a complete description of the core Escherichia coli metabolic model. This model can be analyzed in any computational format (such as MATLAB or Mathematica) based on the information given in this chapter. The core E. coli model is a small-scale model that can be used for educational purposes. It is meant to be used by senior undergraduate and first-year graduate students learning about constraint-based modeling and systems biology. This model has enough reactions and pathways to enable interesting and insightful calculations, but it is also simple enough that the results of such calculations can be understoodeasily.

13C metabolic flux analysis at a genome-scale

Metabolic models used in 13C metabolic flux analysis generally include a limited number of reactions primarily from central metabolism. They typically omit degradation pathways, complete cofactor balances, and atom transition contributions for reactions outside central metabolism. This study addresses the impact on prediction fidelity of scaling-up mapping models to a genome-scale. The core mapping model employed in this study accounts for (75 reactions and 65 metabolites) primarily from central metabolism. The genome-scale metabolic mapping model (GSMM) (697 reaction and 595 metabolites) is constructed using as a basis the iAF1260 model upon eliminating reactions guaranteed not to carry flux based on growth and fermentation data for a minimal glucose growth medium. Labeling data for 17 amino acid fragments obtained from cells fed with glucose labeled at the second carbon was used to obtain fluxes and ranges. Metabolic fluxes and confidence intervals are estimated, for both core and genome-scale mapping models, by minimizing the sum of square of differences between predicted and experimentally measured labeling patterns using the EMU decomposition algorithm. Overall, we find that both topology and estimated values of the metabolic fluxes remain largely consistent between core and GSM model. Stepping up to a genome-scale mapping model leads to wider flux inference ranges for 20 key reactions present in the core model. The glycolysis flux range doubles due to the possibility of active gluconeogenesis, the TCA flux range expanded by 80% due to the availability of a bypass through arginine consistent with labeling data, and the transhydrogenase reaction flux was essentially unresolved due to the presence of as many as five routes for the inter-conversion of NADPH to NADH afforded by the genome-scale model. By globally accounting for ATP demands in the GSMM model the unused ATP decreased drastically with the lower bound matching the maintenance ATP requirement. A non-zero flux for the arginine degradation pathway was identified to meet biomass precursor demands as detailed in the iAF1260 model. Inferred ranges for 81% of the reactions in the genome-scale metabolic (GSM) model varied less than one-tenth of the basis glucose uptake rate (95% confidence test). This is because as many as 411 reactions in the GSM are growth coupled meaning that the single measurement of biomass formation rate locks the reaction flux values. This implies that accurate biomass formation rate and composition are critical for resolving metabolic fluxes away from central metabolism and suggests the importance of biomass composition (re)assessment under different genetic and environmental backgrounds. In addition, the loss of information associated with mapping fluxes from MFA on a core model to a GSM model is quantified.

Physiological, biomass elemental composition and proteomic analyses of Escherichia coli ammonium-limited chemostat growth, and comparison with iron- and glucose-limited chemostat growth

Escherichia coli physiological, biomass elemental composition and proteome acclimations to ammonium-limited chemostat growth were measured at four levels of nutrient scarcity controlled via chemostat dilution rate. These data were compared with published iron- and glucose-limited growth data collected from the same strain and at the same dilution rates to quantify general and nutrient-specific responses. Severe nutrient scarcity resulted in an overflow metabolism with differing organic byproduct profiles based on limiting nutrient and dilution rate. Ammonium-limited cultures secreted up to 35% of the metabolized glucose carbon as organic byproducts with acetate representing the largest fraction; in comparison, iron-limited cultures secreted up to 70 % of the metabolized glucose carbon as lactate, and glucose-limited cultures secreted up to 4% of the metabolized glucose carbon as formate. Biomass elemental composition differed with nutrient limitation; biomass from ammonium-limited cultures had a lower nitrogen content than biomass from either iron- or glucose-limited cultures. Proteomic analysis of central metabolism enzymes revealed that ammonium- and iron-limited cultures had a lower abundance of key tricarboxylic acid (TCA) cycle enzymes and higher abundance of key glycolysis enzymes compared with glucose-limited cultures. The overall results are largely consistent with cellular economics concepts, including metabolic tradeoff theory where the limiting nutrient is invested into essential pathways such as glycolysis instead of higher ATP-yielding, but non-essential, pathways such as the TCA cycle. The data provide a detailed insight into ecologically competitive metabolic strategies selected by evolution, templates for controlling metabolism for bioprocesses and a comprehensive dataset for validating in silico representations of metabolism.

NADPH-generating systems in bacteria and archaea

Reduced nicotinamide adenine dinucleotide phosphate (NADPH) is an essential electron donor in all organisms. It provides the reducing power that drives numerous anabolic reactions, including those responsible for the biosynthesis of all major cell components and many products in biotechnology. The efficient synthesis of many of these products, however, is limited by the rate of NADPH regeneration. Hence, a thorough understanding of the reactions involved in the generation of NADPH is required to increase its turnover through rational strain improvement. Traditionally, the main engineering targets for increasing NADPH availability have included the dehydrogenase reactions of the oxidative pentose phosphate pathway and the isocitrate dehydrogenase step of the tricarboxylic acid (TCA) cycle. However, the importance of alternative NADPH-generating reactions has recently become evident. In the current review, the major canonical and non-canonical reactions involved in the production and regeneration of NADPH in prokaryotes are described, and their key enzymes are discussed. In addition, an overview of how different enzymes have been applied to increase NADPH availability and thereby enhance productivity is provided.

Deoxycytidine production by a metabolically engineered Escherichia coli strain

Background

Rational engineering studies for deoxycytidine production were initiated due to low intracellular levels and tight regulation. To achieve high-level production of deoxycytidine, a useful precursor of decitabine, genes related to feed-back inhibition as well as the biosynthetic pathway were engineered. Additionally, we predicted the impact of individual gene expression levels on a complex metabolic network by microarray analysis. Based on these findings, we demonstrated rational metabolic engineering strategies capable of producing deoxycytidine.

Results

To prepare the deoxycytidine producing strain, we first deleted 3 degradation enzymes in the salvage pathway (deoA, udp, and deoD) and 4 enzymes involved in the branching pathway (dcd, cdd, codA and thyA) to completely eliminate degradation of deoxycytidine. Second, purR, pepA and argR were knocked out to prevent feedback inhibition of CarAB. Third, to enhance influx to deoxycytidine, we investigated combinatorial expression of pyrG, T4 nrdCAB and yfbR. The best strain carried pETGY (pyrG-yfbR) from the possible combinatorial plasmids. The resulting strain showed high deoxycytidine yield (650 mg/L) but co-produced byproducts. To further improve deoxycytidine yield and reduce byproduct formation, pgi was disrupted to generate a sufficient supply of NADPH and ribose. Overall, in shake-flask cultures, the resulting strain produced 967 mg/L of dCyd with decreased byproducts.

Conclusions

We demonstrated that deoxycytidine could be readily achieved by recombineering with biosynthetic genes and regulatory genes, which appeared to enhance the supply of precursors for synthesis of carbamoyl phosphate, based on transcriptome analysis. In addition, we showed that carbon flux rerouting, by disrupting pgi, efficiently improved deoxycytidine yield and decreased byproduct content.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-015-0291-8) contains supplementary material, which is available to authorized users.

Synthetic methylotrophy: engineering the production of biofuels and chemicals based on the biology of aerobic methanol utilization

Synthetic methylotrophy is the development of non-native methylotrophs that can utilize methane and methanol as sole carbon and energy sources or as co-substrates with carbohydrates to produce metabolites as biofuels and chemicals. The availability of methane (from natural gas) and its oxidation product, methanol, has been increasing, while prices have been decreasing, thus rendering them as attractive fermentation substrates. As they are more reduced than most carbohydrates, methane and methanol, as co-substrates, can enhance the yields of biologically produced metabolites. Here we discuss synthetic biology and metabolic engineering strategies based on the native biology of aerobic methylotrophs for developing synthetic strains grown on methanol, with Escherichia coli as the prototype.