Transhydrogenase promotes the robustness and evolvability of E. coli deficient in NADPH production

A citric acid cycle-deficient Escherichia coli as an efficient chassis for aerobic fermentations

Tricarboxylic acid cycle (TCA cycle) plays an important role for aerobic growth of heterotrophic bacteria. Theoretically, eliminating TCA cycle would decrease carbon dissipation and facilitate chemicals biosynthesis. Here, we construct an E. coli strain without a functional TCA cycle that can serve as a versatile chassis for chemicals biosynthesis. We first use adaptive laboratory evolution to recover aerobic growth in minimal medium of TCA cycle-deficient E. coli. Inactivation of succinate dehydrogenase is a key event in the evolutionary trajectory. Supply of succinyl-CoA is identified as the growth limiting factor. By replacing endogenous succinyl-CoA dependent enzymes, we obtain an optimized TCA cycle-deficient E. coli strain. As a proof of concept, the strain is engineered for high-yield production of four separate products. This work enhances our understanding of the role of the TCA cycle in E. coli metabolism and demonstrates the advantages of using TCA cycle-deficient E. coli strain for biotechnological applications.

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 (δ²H) and carbon (δ¹³C) compound specific isotopic values of amino acids, and whether diagnostic isotopic fingerprints are revealed by these paired analyses. We measured δ²H and δ¹³C 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 δ²H values, which collectively vary in composition by more than 400‰, are altered along with δ¹³C 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.

Ligand binding and conformational dynamics of the E. coli nicotinamide nucleotide transhydrogenase revealed by hydrogen/deuterium exchange mass spectrometry

Nicotinamide nucleotide transhydrogenases are integral membrane proteins that utilizes the proton motive force to reduce NADP⁺ to NADPH while converting NADH to NAD⁺. Atomic structures of various transhydrogenases in different ligand-bound states have become available, and it is clear that the molecular mechanism involves major conformational changes. Here we utilized hydrogen/deuterium exchange mass spectrometry (HDX-MS) to map ligand binding sites and analyzed the structural dynamics of E. coli transhydrogenase. We found different allosteric effects on the protein depending on the bound ligand (NAD⁺, NADH, NADP⁺, NADPH). The binding of either NADP⁺ or NADPH to domain III had pronounced effects on the transmembrane helices comprising the proton-conducting channel in domain II. We also made use of cyclic ion mobility separation mass spectrometry (cyclic IMS-MS) to maximize coverage and sensitivity in the transmembrane domain, showing for the first time that this technique can be used for HDX-MS studies. Using cyclic IMS-MS, we increased sequence coverage from 68 % to 73 % in the transmembrane segments. Taken together, our results provide important new insights into the transhydrogenase reaction cycle and demonstrate the benefit of this new technique for HDX-MS to study ligand binding and conformational dynamics in membrane proteins.

Sensitive and Specific Cadmium Biosensor Developed by Reconfiguring Metal Transport and Leveraging Natural Gene Repositories

Whole-cell biosensors are useful for monitoring heavy metal toxicity in public health and ecosystems, but their development has been hindered by intrinsic trade-offs between sensitivity and specificity. Here, we demonstrated an effective engineering solution by building a sensitive, specific, and high-response biosensor for carcinogenic cadmium ions. We genetically programmed the metal transport system of Escherichia coli to enrich intracellular cadmium ions and deprive interfering metal species. We then selected 16 cadmium-sensing transcription factors from the GenBank database and tested their reactivity to 14 metal ions in the engineered E. coli using the expression of the green fluorescent protein as the readout. The resulting cadmium biosensor was highly specific and showed a detection limit of 3 nM, a linear increase in fluorescent intensities from 0 to 200 nM, and a maximal 777-fold signal change. Using this whole-cell biosensor, a smartphone, and low-tech equipment, we developed a simple assay capable of measuring cadmium ions at the same concentration range in irrigation water and human urine. This method is user-friendly and cost-effective, making it affordable to screen large amounts of samples for cadmium toxicity in agriculture and medicine. Moreover, our work highlights natural gene repositories as a treasure chest for bioengineering.

Dynamic control over feedback regulatory mechanisms improves NADPH flux and xylitol biosynthesis in engineered E. coli

We report improved NADPH flux and xylitol biosynthesis in engineered E. coli. Xylitol is produced from xylose via an NADPH dependent reductase. We utilize 2-stage dynamic metabolic control to compare two approaches to optimize xylitol biosynthesis, a stoichiometric approach, wherein competitive fluxes are decreased, and a regulatory approach wherein the levels of key regulatory metabolites are reduced. The stoichiometric and regulatory approaches lead to a 20-fold and 90-fold improvement in xylitol production, respectively. Strains with reduced levels of enoyl-ACP reductase and glucose-6-phosphate dehydrogenase, led to altered metabolite pools resulting in the activation of the membrane bound transhydrogenase and an NADPH generation pathway, consisting of pyruvate ferredoxin oxidoreductase coupled with NADPH dependent ferredoxin reductase, leading to increased NADPH fluxes, despite a reduction in NADPH pools. These strains produced titers of 200 g/L of xylitol from xylose at 86% of theoretical yield in instrumented bioreactors. We expect dynamic control over the regulation of the membrane bound transhydrogenase as well as NADPH production through pyruvate ferredoxin oxidoreductase to broadly enable improved NADPH dependent bioconversions or production via NADPH dependent metabolic pathways.

Global fitness landscapes of the Shine-Dalgarno sequence

Shine-Dalgarno sequences (SD) in prokaryotic mRNA facilitate protein translation by pairing with rRNA in ribosomes. Although conventionally defined as AG-rich motifs, recent genomic surveys reveal great sequence diversity, questioning how SD functions. Here, we determined the molecular fitness (i.e., translation efficiency) of 4⁹ synthetic 9-nt SD genotypes in three distinct mRNA contexts in Escherichia coli. We uncovered generic principles governing the SD fitness landscapes: (1) Guanine contents, rather than canonical SD motifs, best predict the fitness of both synthetic and endogenous SD; (2) the genotype-fitness correlation of SD promotes its evolvability by steadily supplying beneficial mutations across fitness landscapes; and (3) the frequency and magnitude of deleterious mutations increase with background fitness, and adjacent nucleotides in SD show stronger epistasis. Epistasis results from disruption of the continuous base pairing between SD and rRNA. This “chain-breaking” epistasis creates sinkholes in SD fitness landscapes and may profoundly impact the evolution and function of prokaryotic translation initiation and other RNA-mediated processes. Collectively, our work yields functional insights into the SD sequence variation in prokaryotic genomes, identifies a simple design principle to guide bioengineering and bioinformatic analysis of SD, and illuminates the fundamentals of fitness landscapes and molecular evolution.

Evolutionary Repair Experiments as a Window to the Molecular Diversity of Life

Comparative genomics reveals an unexpected diversity in the molecular mechanisms underlying conserved cellular functions, such as DNA replication and cytokinesis. However, the genetic bases and evolutionary processes underlying this 'molecular diversity' remain to be explained. Here, we review a tool to generate alternative mechanisms for conserved cellular functions and test hypotheses concerning the generation of molecular diversity - evolutionary repair experiments, in which laboratory microbial populations adapt in response to a genetic perturbation. We summarize the insights gained from evolutionary repair experiments, the spectrum and dynamics of compensatory mutations, and the alternative molecular mechanisms used to repair perturbed cellular functions. We relate these experiments to the modifications of conserved functions that have occurred outside the laboratory. We end by proposing strategies to improve evolutionary repair experiments as a tool to explore the molecular diversity of life.

High-Efficient Production of (S)-1-[3,5-Bis(trifluoromethyl)phenyl]ethanol via Whole-Cell Catalyst in Deep-Eutectic Solvent-Containing Micro-Aerobic Medium System

The ratio of substrate to catalyst (S/C) is a prime target for the application of asymmetric production of enantiomerically enriched intermediates by whole-cell biocatalyst. In the present study, an attractive increase in S/C was achieved in a natural deep-eutectic solvent (NADES) containing reaction system under microaerobic condition for high production of (S)-1-[3,5-bis(trifluoromethyl)phenyl]ethanol ((S)-3,5-BTPE) with Candida tropicalis 104. In PBS buffer (0.2 M, pH 8.0) at 200 rpm and 30 °C, 79.5 g (Dry Cell Weight, DCW)/L C. tropicalis 104 maintained the same yield of 73.7% for the bioreduction of 3,5-bis(trifluoromethyl)acetophenone (BTAP) under an oxygen-deficient environment compared with oxygen-sufficient conditions, while substrate load increased 4.0-fold (from 50 mM to 200 mM). Furthermore, when choline chloride:trehalose (ChCl:T, 1:1 molar ratio) was introduced into the reaction system for its versatility of increasing cell membrane permeability and declining BTAP cytotoxicity to biocatalyst, the yields were further increased to 86.2% under 200 mM BTAP, or 72.9% at 300 mM BTAP. After the optimization of various reaction parameters involved in the bioreduction, and the amount of biocatalyst and maltose co-substrate remained 79.5 g (DCW)/L and 50 g/L, the S/C for the reduction elevated 6.3 times (3.8 mM/g versus 0.6 mM/g). By altering the respiratory pattern of the whole-cell biocatalyst and exploiting the ChCl:T-containing reaction system, the developed strategy exhibits an attractive potential for enhancing catalytic efficiency of whole-cell-mediated reduction, and provides valuable insight for the development of whole-cell catalysis.

Transcriptome-guided parsimonious flux analysis improves predictions with metabolic networks in complex environments

The metabolic responses of bacteria to dynamic extracellular conditions drives not only the behavior of single species, but also entire communities of microbes. Over the last decade, genome-scale metabolic network reconstructions have assisted in our appreciation of important metabolic determinants of bacterial physiology. These network models have been a powerful force in understanding the metabolic capacity that species may utilize in order to succeed in an environment. Increasingly, an understanding of context-specific metabolism is critical for elucidating metabolic drivers of larger phenotypes and disease. However, previous approaches to use network models in concert with omics data to better characterize experimental systems have met challenges due to assumptions necessary by the various integration platforms or due to large input data requirements. With these challenges in mind, we developed RIPTiDe (Reaction Inclusion by Parsimony and Transcript Distribution) which uses both transcriptomic abundances and parsimony of overall flux to identify the most cost-effective usage of metabolism that also best reflects the cell’s investments into transcription. Additionally, in biological samples where it is difficult to quantify specific growth conditions, it becomes critical to develop methods that require lower amounts of user intervention in order to generate accurate metabolic predictions. Utilizing a metabolic network reconstruction for the model organism Escherichia coli str. K-12 substr. MG1655 (iJO1366), we found that RIPTiDe correctly identifies context-specific metabolic pathway activity without supervision or knowledge of specific media conditions. We also assessed the application of RIPTiDe to in vivo metatranscriptomic data where E. coli was present at high abundances, and found that our approach also effectively predicts metabolic behaviors of host-associated bacteria. In the setting of human health, understanding metabolic changes within bacteria in environments where growth substrate availability is difficult to quantify can have large downstream impacts on our ability to elucidate molecular drivers of disease-associated dysbiosis across the microbiota. Our results indicate that RIPTiDe may have potential to provide understanding of context-specific metabolism of bacteria within complex communities.

Transcriptomic analyses of bacteria have become instrumental to our understanding of their responses to changes in their environment. While traditional analyses have been informative, leveraging these datasets within genome-scale metabolic network reconstructions (GENREs) can provide greatly improved context for shifts in pathway utilization and downstream/upstream ramifications for changes in metabolic regulation. Many previous techniques for GENRE transcript integration have focused on creating maximum consensus with input datasets, but these approaches were recently shown to generate less accurate metabolic predictions than a transcript-agnostic method of flux minimization (pFBA), which identifies the most efficient/economic patterns of metabolism given certain growth constraints. Despite this success, growth conditions are not always easily quantifiable and highlights the need for novel platforms that build from these findings. Our new method, RIPTiDe, combines these concepts and utilizes overall minimization of flux weighted by transcriptomic analysis to identify the most energy efficient pathways to achieve growth that include more highly transcribed enzymes, without previous insight into extracellular conditions. Utilizing a well-studied GENRE from Escherichia coli, we demonstrate that this new approach correctly predicts patterns of metabolism utilizing a variety of both in vitro and in vivo transcriptomes. This platform could be important for revealing context-specific bacterial phenotypes in line with governing principles of adaptive evolution, that drive disease manifestation or interactions between microbes.

RNA accumulation in Candida tropicalis based on cofactor engineering

Redox cofactors play an important role in biosynthetic and catabolic reactions and the transfer of energy for the cell. Therefore, studying the relationship between cofactor perturbation and metabolism is a useful approach to improve the yield of target products. To study RNA accumulation and metabolism when intracellular cofactor balance was impaired, the water-forming NADH oxidase (NoxE) from Lactococcus lactis and membrane-bound transhydrogenase (PntAB) from Escherichia coli were expressed in Candidatropicalis no. 121, respectively. Expression of noxE significantly decreased the intracellular NADH/NAD+ ratio, but the NADPH/NADP+ ratio did not differ significantly. PntAB increased the intracellular NADH pool, while the NADPH/NADP+ ratio decreased. The perturbation of the cofactors caused a large redistribution of metabolic fluxes. The biomass and RNA content decreased by 11.0% and 10.6% in pAUR-noxE strain, respectively, while the RNA content increased by 5.5% and the biomass showed no signification difference in pAUR-pntAB strain. Expression of noxE and pntAB led to decreases and increases in the ATP concentration and yield of RNA, respectively, which also indicated that ATP plays an important role in the RNA biosynthesis.

Metabolic perturbations in mutants of glucose transporters and their applications in metabolite production in Escherichia coli

BACKGROUND

Most microorganisms have evolved to maximize growth rate, with rapid consumption of carbon sources from the surroundings. However, fast growing phenotypes usually feature secretion of organic compounds. For example, E. coli mainly produced acetate in fast growing condition such as glucose rich and aerobic condition, which is troublesome for metabolic engineering because acetate causes acidification of surroundings, growth inhibition and decline of production yield. The overflow metabolism can be alleviated by reducing glucose uptake rate.

RESULTS

As glucose transporters or their subunits were knocked out in E. coli, the growth and glucose uptake rates decreased and biomass yield was improved. Alteration of intracellular metabolism caused by the mutations was investigated with transcriptome analysis and 13C metabolic flux analysis (13C MFA). Various transcriptional and metabolic perturbations were identified in the sugar transporter mutants. Transcription of genes related to glycolysis, chemotaxis, and flagella synthesis was downregulated, and that of gluconeogenesis, Krebs cycle, alternative transporters, quorum sensing, and stress induced proteins was upregulated in the sugar transporter mutants. The specific production yields of value-added compounds (enhanced green fluorescent protein, γ-aminobutyrate, lycopene) were improved significantly in the sugar transporter mutants.

CONCLUSIONS

The elimination of sugar transporter resulted in alteration of global gene expression and redirection of carbon flux distribution, which was purposed to increase energy yield and recycle carbon sources. When the pathways for several valuable compounds were introduced to mutant strains, specific yield of them were highly improved. These results showed that controlling the sugar uptake rate is a good strategy for ameliorating 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.

Equilibrium of the intracellular redox state for improving cell growth and L-lysine yield of Corynebacterium glutamicum by optimal cofactor swapping

Background

NAD(H/⁺) and NADP(H/⁺) are the most important redox cofactors in bacteria. However, the intracellular redox balance is in advantage of the cell growth and production of NAD(P)H-dependent products.

Results

In this paper, we rationally engineered glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and isocitrate dehydrogenase (IDH) to switch the nucleotide-cofactor specificity resulting in an increase in final titer [from 85.6 to 121.4 g L⁻¹] and carbon yield [from 0.33 to 0.46 g (g glucose)⁻¹] of l-lysine in strain RGI in fed-batch fermentation. To do this, we firstly analyzed the production performance of original strain JL-6, indicating that the imbalance of intracellular redox was the limiting factor for l-lysine production. Subsequently, we modified the native GAPDH and indicated that recombinant strain RG with nonnative NADP-GAPDH dramatically changed the intracellular levels of NADH and NADPH. However, l-lysine production did not significantly increase because cell growth was harmed at low NADH level. Lastly, the nonnative NAD-IDH was introduced in strain RG to increase the NADH availability and to equilibrate the intracellular redox. The resulted strain RGI showed the stable ratio of NADPH/NADH at about 1.00, which in turn improved cell growth (μ_max. = 0.31 h⁻¹) and l-lysine productivity (q_{Lys, max.} = 0.53 g g⁻¹ h⁻¹) as compared with strain RG (μ_max. = 0.14 h⁻¹ and q_{Lys, max.} = 0.42 g g⁻¹ h⁻¹).

Conclusions

This is the first report of balancing the intracellular redox state by switching the nucleotide-cofactor specificity of GAPDH and IDH, thereby improving cell growth and l-lysine production.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1114-0) contains supplementary material, which is available to authorized users.

Mutations in sigma 70 transcription factor improves expression of functional eukaryotic membrane proteins in Escherichia coli

Eukaryotic integral membrane proteins (IMPs) are difficult to study due to low functional expression levels. To investigate factors for efficient biogenesis of eukaryotic IMPs in the prokaryotic model organism Escherichia coli, important, e.g., for isotope-labeling for NMR, we selected for E. coli cells expressing high levels of functional G protein-coupled receptors (GPCRs) by FACS. Utilizing an E. coli strain library with all non-essential genes systematically deleted, we unexpectedly discovered upon whole-genome sequencing that the improved phenotype was not conferred by the deleted genes but by various subtle alterations in the “housekeeping” sigma 70 factor (RpoD). When analyzing effects of the rpoD mutations at the transcriptome level we found that toxic effects incurred on wild-type E. coli during receptor expression were diminished by two independent and synergistic effects: a slower but longer-lasting GPCR biosynthesis and an optimized transcriptional pattern, augmenting growth and expression at low temperature, setting the basis for further bacterial strain engineering.

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.

How adaptive evolution reshapes metabolism to improve fitness: recent advances and future outlook

Adaptive laboratory evolution (ALE) has emerged as a powerful tool in basic microbial research and strain development. In the context of metabolic science and engineering, it has been applied to study gene knockout responses, expand substrate ranges, improve tolerance to process conditions, and to improve productivity via designed growth coupling. In recent years, advancements in ALE methods and systems biology measurement technologies, particularly genome sequencing and ¹³C metabolic flux analysis (¹³C-MFA), have enabled detailed study of the mechanisms and dynamics of evolving metabolism. In this review, we discuss a range of studies that have applied flux analysis to adaptively evolved strains, as well as modeling frameworks developed to predict and interpret evolved fluxes. These efforts link mutations to fitness-enhanced phenotypes, identify bottlenecks and approaches to resolve them, and address systems concepts such as optimality.

Multiple Optimal Phenotypes Overcome Redox and Glycolytic Intermediate Metabolite Imbalances in Escherichia coli pgi Knockout Evolutions

A mechanistic understanding of how new phenotypes develop to overcome the loss of a gene product provides valuable insight on both the metabolic and regulatory functions of the lost gene. The pgi gene, whose product catalyzes the second step in glycolysis, was deleted in a growth-optimized Escherichia coli K-12 MG1655 strain. The initial knockout (KO) strain exhibited an 80% drop in growth rate that was largely recovered in eight replicate, but phenotypically distinct, cultures after undergoing adaptive laboratory evolution (ALE). Multi-omic data sets showed that the loss of pgi substantially shifted pathway usage, leading to a redox and sugar phosphate stress response. These stress responses were overcome by unique combinations of innovative mutations selected for by ALE. Thus, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.IMPORTANCE A mechanistic understanding of how microbes are able to overcome the loss of a gene through regulatory and metabolic changes is not well understood. Eight independent adaptive laboratory evolution (ALE) experiments with pgi knockout strains resulted in eight phenotypically distinct endpoints that were able to overcome the gene loss. Utilizing multi-omics analysis, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.

Growth Adaptation of gnd and sdhCB Escherichia coli Deletion Strains Diverges From a Similar Initial Perturbation of the Transcriptome

Adaptive laboratory evolution (ALE) has emerged as a new approach with which to pursue fundamental biological inquiries and, in particular, new insights into the systemic function of a gene product. Two E. coli knockout strains were constructed: one that blocked the Pentose Phosphate Pathway (gnd KO) and one that decoupled the TCA cycle from electron transport (sdhCDAB KO). Despite major perturbations in central metabolism, minimal growth rate changes were found in the two knockout strains. More surprisingly, many similarities were found in their initial transcriptomic states that could be traced to similarly perturbed metabolites despite the differences in the network location of the gene perturbations and concomitant re-routing of pathway fluxes around these perturbations. However, following ALE, distinct metabolomic and transcriptomic states were realized. These included divergent flux and gene expression profiles in the gnd and sdhCDAB KOs to overcome imbalances in NADPH production and nitrogen/sulfur assimilation, respectively, that were not obvious limitations of growth in the unevolved knockouts. Therefore, this work demonstrates that ALE provides a productive approach to reveal novel insights of gene function at a systems level that cannot be found by observing the fresh knockout alone.

Improving methyl ketone production in Escherichia coli by heterologous expression of NADH-dependent FabG

We previously engineered Escherichia coli to overproduce medium- to long-chain saturated and monounsaturated methyl ketones, which could potentially be applied as diesel fuel blending agents or in the flavor and fragrance industry. Recent efforts at strain optimization have focused on cofactor balance, as fatty acid-derived pathways face the systematic metabolic challenge of net NADPH consumption (in large part, resulting from the key fatty acid biosynthetic enzyme FabG [β-ketoacyl-ACP reductase]) and net NADH production. In this study, we attempted to mitigate cofactor imbalance by heterologously expressing NADH-dependent, rather than NADPH-dependent, versions of FabG identified in previous studies. Of the four NADH-dependent versions of FabG tested in our previously best-reported methyl ketone-producing strain (EGS1895), the version from Acholeplasma laidlawii (Al_FabG) showed the greatest increase in methyl ketone yield in shake flasks (35-75% higher than for an RFP negative-control strain, depending on sugar loading). An improved strain (EGS2920) attained methyl ketone titers during fed-batch fermentation of 5.4 ± 0.5 g/L, which were, on average, ca. 40% greater than those for the base strain (EGS1895) under fermentation conditions optimized in this study. Shotgun proteomic data for strains EGS2920 and EGS1895 during fed-batch fermentation were consistent with the goal of alleviating NADPH limitation through expression of Al_FabG. For example, relative to strain EGS1895, strain EGS2920 significantly upregulated glucose-6-phosphate isomerase (directing flux into glycolysis rather than the NADPH-producing pentose phosphate pathway) and downregulated MaeB (a NADP⁺ -dependent malate dehydrogenase). Overall, the results suggest that heterologous expression of NADH-dependent FabG in E. coli may improve sustained production of fatty acid-derived renewable fuels and chemicals.

NADPH metabolism: a survey of its theoretical characteristics and manipulation strategies in amino acid biosynthesis

Reduced nicotinamide adenine nucleotide phosphate (NADPH), which is one of the key cofactors in the metabolic network, plays an important role in the biochemical reactions, and physiological function of amino acid-producing strains. The manipulation of NADPH availability and form is an efficient and easy method of redirecting the carbon flux to the amino acid biosynthesis in industrial strains. In this review, we survey the metabolic mode of NADPH. Furthermore, we summarize the research developments in the understanding of the relationship between NADPH metabolism and amino acid biosynthesis. Detailed strategies to manipulate NADPH availability are addressed based on this knowledge. Finally, the uses of NADPH manipulation strategies to enhance the metabolic function of amino acid-producing strains are discussed.

Control of primary metabolism by a virulence regulatory network promotes robustness in a plant pathogen

Robustness is a key system-level property of living organisms to maintain their functions while tolerating perturbations. We investigate here how a regulatory network controlling multiple virulence factors impacts phenotypic robustness of a bacterial plant pathogen. We reconstruct a cell-scale model of Ralstonia solanacearum connecting a genome-scale metabolic network, a virulence macromolecule network, and a virulence regulatory network, which includes 63 regulatory components. We develop in silico methods to quantify phenotypic robustness under a broad set of conditions in high-throughput simulation analyses. This approach reveals that the virulence regulatory network exerts a control of the primary metabolism to promote robustness upon infection. The virulence regulatory network plugs into the primary metabolism mainly through the control of genes likely acquired via horizontal gene transfer, which results in a functional overlay with ancestral genes. These results support the view that robustness may be a selected trait that promotes pathogenic fitness upon infection.

Thermodynamic analysis of computed pathways integrated into the metabolic networks of E. coli and Synechocystis reveals contrasting expansion potential

Introducing biosynthetic pathways into an organism is both reliant on and challenged by endogenous biochemistry. Here we compared the expansion potential of the metabolic network in the photoautotroph Synechocystis with that of the heterotroph E. coli using the novel workflow POPPY (Prospecting Optimal Pathways with PYthon). First, E. coli and Synechocystis metabolomic and fluxomic data were combined with metabolic models to identify thermodynamic constraints on metabolite concentrations (NET analysis). Then, thousands of automatically constructed pathways were placed within each network and subjected to a network-embedded variant of the max-min driving force analysis (NEM). We found that the networks had different capabilities for imparting thermodynamic driving forces toward certain compounds. Key metabolites were constrained differently in Synechocystis due to opposing flux directions in glycolysis and carbon fixation, the forked tri-carboxylic acid cycle, and photorespiration. Furthermore, the lysine biosynthesis pathway in Synechocystis was identified as thermodynamically constrained, impacting both endogenous and heterologous reactions through low 2-oxoglutarate levels. Our study also identified important yet poorly covered areas in existing metabolomics data and provides a reference for future thermodynamics-based engineering in Synechocystis and beyond. The POPPY methodology represents a step in making optimal pathway-host matches, which is likely to become important as the practical range of host organisms is diversified.

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.

Improving metabolic efficiency of the reverse beta-oxidation cycle by balancing redox cofactor requirement

Previous studies have made many exciting achievements on pushing the functional reversal of beta-oxidation cycle (r-BOX) to more widespread adoption for synthesis of a wide variety of fuels and chemicals. However, the redox cofactor requirement for the efficient operation of r-BOX remains unclear. In this work, the metabolic efficiency of r-BOX for medium-chain fatty acid (C₆-C₁₀, MCFA) production was optimized by redox cofactor engineering. Stoichiometric analysis of the r-BOX pathway and further experimental examination identified NADH as a crucial determinant of r-BOX process yield. Furthermore, the introduction of formate dehydrogenase from Candida boidinii using fermentative inhibitor byproduct formate as a redox NADH sink improved MCFA titer from initial 1.2g/L to 3.1g/L. Moreover, coupling of increasing the supply of acetyl-CoA with NADH to achieve fermentative redox balance enabled product synthesis at maximum titers. To this end, the acetate re-assimilation pathway was further optimized to increase acetyl-CoA availability associated with the new supply of NADH. It was found that the acetyl-CoA synthetase activity and intracellular ATP levels constrained the activity of acetate re-assimilation pathway, and 4.7g/L of MCFA titer was finally achieved after alleviating these two limiting factors. To the best of our knowledge, this represented the highest titer reported to date. These results demonstrated that the key constraint of r-BOX was redox imbalance and redox engineering could further unleash the lipogenic potential of this cycle. The redox engineering strategies could be applied to acetyl-CoA-derived products or other bio-products requiring multiple redox cofactors for biosynthesis.

Genomic and Transcriptomic Analysis of Growth-Supporting Dehalogenation of Chlorinated Methanes in Methylobacterium

Bacterial adaptation to growth with toxic halogenated chemicals was explored in the context of methylotrophic metabolism of Methylobacterium extorquens, by comparing strains CM4 and DM4, which show robust growth with chloromethane and dichloromethane, respectively. Dehalogenation of chlorinated methanes initiates growth-supporting degradation, with intracellular release of protons and chloride ions in both cases. The core, variable and strain-specific genomes of strains CM4 and DM4 were defined by comparison with genomes of non-dechlorinating strains. In terms of gene content, adaptation toward dehalogenation appears limited, strains CM4 and DM4 sharing between 75 and 85% of their genome with other strains of M. extorquens. Transcript abundance in cultures of strain CM4 grown with chloromethane and of strain DM4 grown with dichloromethane was compared to growth with methanol as a reference C₁ growth substrate. Previously identified strain-specific dehalogenase-encoding genes were the most transcribed with chlorinated methanes, alongside other genes encoded by genomic islands (GEIs) and plasmids involved in growth with chlorinated compounds as carbon and energy source. None of the 163 genes shared by strains CM4 and DM4 but not by other strains of M. extorquens showed higher transcript abundance in cells grown with chlorinated methanes. Among the several thousand genes of the M. extorquens core genome, 12 genes were only differentially abundant in either strain CM4 or strain DM4. Of these, 2 genes of known function were detected, for the membrane-bound proton translocating pyrophosphatase HppA and the housekeeping molecular chaperone protein DegP. This indicates that the adaptive response common to chloromethane and dichloromethane is limited at the transcriptional level, and involves aspects of the general stress response as well as of a dehalogenation-specific response to intracellular hydrochloric acid production. Core genes only differentially abundant in either strain CM4 or strain DM4 total 13 and 58 CDS, respectively. Taken together, the obtained results suggest different transcriptional responses of chloromethane- and dichloromethane-degrading M. extorquens strains to dehalogenative metabolism, and substrate- and pathway-specific modes of growth optimization with chlorinated methanes.

Functional Trade-Offs in Promiscuous Enzymes Cannot Be Explained by Intrinsic Mutational Robustness of the Native Activity

The extent to which an emerging new function trades off with the original function is a key characteristic of the dynamics of enzyme evolution. Various cases of laboratory evolution have unveiled a characteristic trend; a large increase in a new, promiscuous activity is often accompanied by only a mild reduction of the native, original activity. A model that associates weak trade-offs with “evolvability” was put forward, which proposed that enzymes possess mutational robustness in the native activity and plasticity in promiscuous activities. This would enable the acquisition of a new function without compromising the original one, reducing the benefit of early gene duplication and therefore the selection pressure thereon. Yet, to date, no experimental study has examined this hypothesis directly. Here, we investigate the causes of weak trade-offs by systematically characterizing adaptive mutations that occurred in two cases of evolutionary transitions in enzyme function: (1) from phosphotriesterase to arylesterase, and (2) from atrazine chlorohydrolase to melamine deaminase. Mutational analyses in various genetic backgrounds revealed that, in contrast to the prevailing model, the native activity is less robust to mutations than the promiscuous activity. For example, in phosphotriesterase, the deleterious effect of individual mutations on the native phosphotriesterase activity is much larger than their positive effect on the promiscuous arylesterase activity. Our observations suggest a revision of the established model: weak trade-offs are not caused by an intrinsic robustness of the native activity and plasticity of the promiscuous activity. We propose that upon strong adaptive pressure for the new activity without selection against the original one, selected mutations will lead to the largest possible increases in the new function, but whether and to what extent they decrease the old function is irrelevant, creating a bias towards initially weak trade-offs and the emergence of generalist enzymes.

Understanding how enzymes evolve is a fundamental question that can help us decipher not only the mechanisms of evolution on a higher level, i.e., whole organisms, but also advances our knowledge of sequence-structure-function relationships as a guide to artificial evolution in the test tube. An important yet unexplained phenomenon occurs during the evolution of a new enzymatic function; it has been observed that new and ancestral functions often trade-off only weakly, meaning the original native activity is initially maintained at a high level despite drastic improvement of the new promiscuous activity. It has previously been proposed that weak trade-offs occur because the native activity is robust to mutations while the promiscuous activity is not. However, the present work contradicts this hypothesis, based on the detailed characterization of mutational effects on both activities in two examples of enzyme evolution. We propose an alternative explanation: the weak activity trade-off is consistent with being a by-product of strong selection for the new activity rather than an intrinsic property of the native activity.

Distinct transcriptional regulation of the two Escherichia coli transhydrogenases PntAB and UdhA

Transhydrogenases catalyse interconversion of the redox cofactors NADH and NADPH, thereby conveying metabolic flexibility to balance catabolic NADPH formation with anabolic or stress-based consumption of NADPH. Escherichia coli is one of the very few microbes that possesses two isoforms: the membrane-bound, proton-translocating transhydrogenase PntAB and the cytosolic, energy-independent transhydrogenase UdhA. Despite their physiological relevance, we have only fragmented information on their regulation and the signals coordinating their counteracting activities. Here we investigated PntAB and UdhA regulation by studying transcriptional responses to environmental and genetic perturbations. By testing pntAB and udhA GFP reporter constructs in the background of WT E. coli and 62 transcription factor mutants during growth on different carbon sources, we show distinct transcriptional regulation of the two transhydrogenase promoters. Surprisingly, transhydrogenase regulation was independent of the actual catabolic overproduction or underproduction of NADPH but responded to nutrient levels and growth rate in a fashion that matches the cellular need for the redox cofactors NADPH and/or NADH. Specifically, the identified transcription factors Lrp, ArgP and Crp link transhydrogenase expression to particular amino acids and intracellular concentrations of cAMP. The overall identified set of regulators establishes a primarily biosynthetic role for PntAB and link UdhA to respiration.

Pyridine nucleotide transhydrogenases enable redox balance of Pseudomonas putida during biodegradation of aromatic compounds

The metabolic versatility of the soil bacterium Pseudomonas putida is reflected by its ability to execute strong redox reactions (e.g., mono- and di-oxygenations) on aromatic substrates. Biodegradation of aromatics occurs via the pathway encoded in the archetypal TOL plasmid pWW0, yet the effect of running such oxidative route on redox balance against the background metabolism of P. putida remains unexplored. To answer this question, the activity of pyridine nucleotide transhydrogenases (that catalyze the reversible interconversion of NADH and NADPH) was inspected under various physiological and oxidative stress regimes. The genome of P. putida KT2440 encodes a soluble transhydrogenase (SthA) and a membrane-bound, proton-pumping counterpart (PntAB). Mutant strains, lacking sthA and/or pntAB, were subjected to a panoply of genetic, biochemical, phenomic and functional assays in cells grown on customary carbon sources (e.g., citrate) versus difficult-to-degrade aromatic substrates. The results consistently indicated that redox homeostasis is compromised in the transhydrogenases-defective variant, rendering the mutant sensitive to oxidants. This metabolic deficiency was, however, counteracted by an increase in the activity of NADP⁺ -dependent dehydrogenases in central carbon metabolism. Taken together, these observations demonstrate that transhydrogenases enable a redox-adjusting mechanism that comes into play when biodegradation reactions are executed to metabolize unusual carbon compounds.

13C-Metabolic Flux Analysis: An Accurate Approach to Demystify Microbial Metabolism for Biochemical Production

Metabolic engineering of various industrial microorganisms to produce chemicals, fuels, and drugs has raised interest since it is environmentally friendly, sustainable, and independent of nonrenewable resources. However, microbial metabolism is so complex that only a few metabolic engineering efforts have been able to achieve a satisfactory yield, titer or productivity of the target chemicals for industrial commercialization. In order to overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been continuously developed and widely applied to rigorously investigate cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, many ¹³C-MFA studies have been performed in academic labs and biotechnology industries to pinpoint key issues related to microbe-based chemical production. Insightful information about the metabolic rewiring has been provided to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this review, we will introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied via integration with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production for various host microorganisms

Protein Homeostasis Imposes a Barrier on Functional Integration of Horizontally Transferred Genes in Bacteria

Horizontal gene transfer (HGT) plays a central role in bacterial evolution, yet the molecular and cellular constraints on functional integration of the foreign genes are poorly understood. Here we performed inter-species replacement of the chromosomal folA gene, encoding an essential metabolic enzyme dihydrofolate reductase (DHFR), with orthologs from 35 other mesophilic bacteria. The orthologous inter-species replacements caused a marked drop (in the range 10–90%) in bacterial growth rate despite the fact that most orthologous DHFRs are as stable as E.coli DHFR at 37°C and are more catalytically active than E. coli DHFR. Although phylogenetic distance between E. coli and orthologous DHFRs as well as their individual molecular properties correlate poorly with growth rates, the product of the intracellular DHFR abundance and catalytic activity (k_cat/K_M), correlates strongly with growth rates, indicating that the drop in DHFR abundance constitutes the major fitness barrier to HGT. Serial propagation of the orthologous strains for ~600 generations dramatically improved growth rates by largely alleviating the fitness barriers. Whole genome sequencing and global proteome quantification revealed that the evolved strains with the largest fitness improvements have accumulated mutations that inactivated the ATP-dependent Lon protease, causing an increase in the intracellular DHFR abundance. In one case DHFR abundance increased further due to mutations accumulated in folA promoter, but only after the lon inactivating mutations were fixed in the population. Thus, by apparently distinguishing between self and non-self proteins, protein homeostasis imposes an immediate and global barrier to the functional integration of foreign genes by decreasing the intracellular abundance of their products. Once this barrier is alleviated, more fine-tuned evolution occurs to adjust the function/expression of the transferred proteins to the constraints imposed by the intracellular environment of the host organism.

Horizontal gene transfer (HGT) is central to bacterial evolution. The outcome of an HGT event (fixation in a population, elimination, or separation as a subdominant clone) depends not only on the availability of a new gene but crucially on the fitness cost or benefit of the genomic incorporation of the foreign gene and its expression in recipient bacteria. Here we studied the fitness landscape for inter-species chromosomal replacement of an essential protein, dihydrofolate reductase (DHFR) encoded by the folA gene, by its orthologs from other mesophilic bacteria. We purified and biochemically characterized 33 out of 35 orthologous DHFRs and found that most of them are stable and more catalytically active than E. coli DHFR. However, the inter-species replacement of DHFR caused significant fitness loss for most transgenic strains due to low abundance of orthologous DHFRs in E. coli cytoplasm. Laboratory evolution resulted in an increase in orthologous DHFR abundance leading to a dramatic fitness improvement. Genomic and proteomic analyses of “naive” and evolved strains suggest a new function of protein homeostasis to discriminate between “self” and “non-self” proteins, thus creating fitness barriers to HGT.

Review and Hypothesis. New insights into the reaction mechanism of transhydrogenase: Swivelling the dIII component may gate the proton channel

The membrane protein transhydrogenase in animal mitochondria and bacteria couples reduction of NADP⁺ by NADH to proton translocation. Recent X-ray data on Thermus thermophilus transhydrogenase indicate a significant difference in the orientations of the two dIII components of the enzyme dimer (Leung et al., 2015). The character of the orientation change, and a review of information on the kinetics and thermodynamics of transhydrogenase, indicate that dIII swivelling might assist in the control of proton gating by the redox state of bound NADP⁺/NADPH during enzyme turnover.