Kinetic modelling of the Zymomonas mobilis Entner-Doudoroff pathway: insights into control and functionality

A new Zymomonas mobilis platform strain for the efficient production of chemicals

BACKGROUND

Zymomonas mobilis is well known for its outstanding ability to produce ethanol with both high specific productivity and with high yield close to the theoretical maximum. The key enzyme in the ethanol production pathway is the pyruvate decarboxylase (PDC) which is converting pyruvate to acetaldehyde. Since it is widely considered that its gene pdc is essential, metabolic engineering strategies aiming to produce other compounds derived from pyruvate need to find ways to reduce PDC activity.

RESULTS

Here, we present a new platform strain (sGB027) of Z. mobilis in which the native promoter of pdc was replaced with the IPTG-inducible PT7A1, allowing for a controllable expression of pdc. Expression of lactate dehydrogenase from E. coli in sGB027 allowed the production of D-lactate with, to the best of our knowledge, the highest reported specific productivity of any microbial lactate producer as well as with the highest reported lactate yield for Z. mobilis so far. Additionally, by expressing the L-alanine dehydrogenase of Geobacillus stearothermophilus in sGB027 we produced L-alanine, further demonstrating the potential of sGB027 as a base for the production of compounds other than ethanol.

CONCLUSION

We demonstrated that our new platform strain can be an excellent starting point for the efficient production of various compounds derived from pyruvate with Z. mobilis and can thus enhance the establishment of this organism as a workhorse for biotechnological production processes.

Improving pathway prediction accuracy of constraints-based metabolic network models by treating enzymes as microcompartments

Metabolic network models have become increasingly precise and accurate as the most widespread and practical digital representations of living cells. The prediction functions were significantly expanded by integrating cellular resources and abiotic constraints in recent years. However, if unreasonable modeling methods were adopted due to a lack of consideration of biological knowledge, the conflicts between stoichiometric and other constraints, such as thermodynamic feasibility and enzyme resource availability, would lead to distorted predictions. In this work, we investigated a prediction anomaly of EcoETM, a constraints-based metabolic network model, and introduced the idea of enzyme compartmentalization into the analysis process. Through rational combination of reactions, we avoid the false prediction of pathway feasibility caused by the unrealistic assumption of free intermediate metabolites. This allowed us to correct the pathway structures of l-serine and l-tryptophan. A specific analysis explains the application method of the EcoETM-like model and demonstrates its potential and value in correcting the prediction results in pathway structure by resolving the conflict between different constraints and incorporating the evolved roles of enzymes as reaction compartments. Notably, this work also reveals the trade-off between product yield and thermodynamic feasibility. Our work is of great value for the structural improvement of constraints-based models.

Zymo-Parts: A Golden Gate Modular Cloning Toolbox for Heterologous Gene Expression in Zymomonas mobilis

Zymomonas mobilis is a microorganism with extremely high sugar consumption and ethanol production rates and is generally considered to hold great potential for biotechnological applications. However, its genetic engineering is still difficult, hampering the efficient construction of genetically modified strains. In this work, we present Zymo-Parts, a modular toolbox based on Golden-Gate cloning offering a collection of promoters (including native, inducible, and synthetic constitutive promoters of varying strength), an array of terminators and several synthetic ribosomal binding sites and reporter genes. All these parts can be combined in an efficient and flexible way to achieve a desired level of gene expression, either from plasmids or via genome integration. Use of the GoldenBraid-based system also enables an assembly of operons consisting of up to five genes. We present the basic structure of the Zymo-Parts cloning system, characterize several constitutive and inducible promoters, and exemplify the construction of an operon and of chromosomal integration of a reporter gene. Finally, we demonstrate the power and utility of the Zymo-Parts toolbox for metabolic engineering applications by overexpressing a heterologous gene encoding for the lactate dehydrogenase of Escherichia coli to achieve different levels of lactate production in Z. mobilis.

The potential use of Zymomonas mobilis for the food industry

Zymomonas mobilis is a gram-negative facultative anaerobic spore, which is generally recognized as a safe. As a promising ethanologenic organism for large-scale bio-ethanol production, Z. mobilis has also shown a good application prospect in food processing and food additive synthesis for its unique physiological characteristics and excellent industrial characteristics. It not only has obvious advantages in food processing and becomes the biorefinery chassis cell for food additives, but also has a certain healthcare effect on human health. Until to now, most of the research is still in theory and laboratory scale, and further research is also needed to achieve industrial production. This review summarized the physiological characteristics and advantages of Z. mobilis in food industry for the first time and further expounds its research status in food industry from three aspects of food additive synthesis, fermentation applications, and prebiotic efficacy, it will provide a theoretical basis for its development and applications in food industry. This review also discussed the shortcomings of its practical applications in the current food industry, and explored other ways to broaden the applications of Z. mobilis in the food industry, to promote its applications in food processing.

Evolution of diversity in metabolic strategies

Understanding the origin and maintenance of biodiversity is a fundamental problem. Many theoretical approaches have been investigating ecological interactions, such as competition, as potential drivers of diversification. Classical consumer-resource models predict that the number of coexisting species should not exceed the number of distinct resources, a phenomenon known as the competitive exclusion principle. It has recently been argued that including physiological tradeoffs in consumer-resource models can lead to violations of this principle and to ecological coexistence of very high numbers of species. Here, we show that these results crucially depend on the functional form of the tradeoff. We investigate the evolutionary dynamics of resource use constrained by tradeoffs and show that if the tradeoffs are non-linear, the system either does not diversify or diversifies into a number of coexisting species that do not exceed the number of resources. In particular, very high diversity can only be observed for linear tradeoffs.

Thermodynamic and Kinetic Modeling of Co-utilization of Glucose and Xylose for 2,3-BDO Production by Zymomonas mobilis

Prior engineering of the ethanologen Zymomonas mobilis has enabled it to metabolize xylose and to produce 2,3-butanediol (2,3-BDO) as a dominant fermentation product. When co-fermenting with xylose, glucose is preferentially utilized, even though xylose metabolism generates ATP more efficiently during 2,3-BDO production on a BDO-mol basis. To gain a deeper understanding of Z. mobilis metabolism, we first estimated the kinetic parameters of the glucose facilitator protein of Z. mobilis by fitting a kinetic uptake model, which shows that the maximum transport capacity of glucose is seven times higher than that of xylose, and glucose is six times more affinitive to the transporter than xylose. With these estimated kinetic parameters, we further compared the thermodynamic driving force and enzyme protein cost of glucose and xylose metabolism. It is found that, although 20% more ATP can be yielded stoichiometrically during xylose utilization, glucose metabolism is thermodynamically more favorable with 6% greater cumulative Gibbs free energy change, more economical with 37% less enzyme cost required at the initial stage and sustains the advantage of the thermodynamic driving force and protein cost through the fermentation process until glucose is exhausted. Glucose-6-phosphate dehydrogenase (g6pdh), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and phosphoglycerate mutase (pgm) are identified as thermodynamic bottlenecks in glucose utilization pathway, as well as two more enzymes of xylose isomerase and ribulose-5-phosphate epimerase in xylose metabolism. Acetolactate synthase is found as potential engineering target for optimized protein cost supporting unit metabolic flux. Pathway analysis was then extended to the core stoichiometric matrix of Z. mobilis metabolism. Growth was simulated by dynamic flux balance analysis and the model was validated showing good agreement with experimental data. Dynamic FBA simulations suggest that a high agitation is preferable to increase 2,3-BDO productivity while a moderate agitation will benefit the 2,3-BDO titer. Taken together, this work provides thermodynamic and kinetic insights of Z. mobilis metabolism on dual substrates, and guidance of bioengineering efforts to increase hydrocarbon fuel production.

Ethanologenesis and respiration in a pyruvate decarboxylase-deficient Zymomonas mobilis

Objective

Zymomonas mobilis is an alpha-proteobacterium with a rapid ethanologenic pathway, involving Entner–Doudoroff (E–D) glycolysis, pyruvate decarboxylase (Pdc) and two alcohol dehydrogenase (ADH) isoenzymes. Pyruvate is the end-product of the E–D pathway and the substrate for Pdc. Construction and study of Pdc-deficient strains is of key importance for Z. mobilis metabolic engineering, because the pyruvate node represents the central branching point, most novel pathways divert from ethanol synthesis. In the present work, we examined the aerobic metabolism of a strain with partly inactivated Pdc.

Results

Relative to its parent strain the mutant produced more pyruvate. Yet, it also yielded more acetaldehyde, the product of the Pdc reaction and the substrate for ADH, although the bulk ADH activity was similar in both strains, while the Pdc activity in the mutant was reduced by half. Simulations with the kinetic model of Z. mobilis E-D pathway indicated that, for the observed acetaldehyde to ethanol production ratio in the mutant, the ratio between its respiratory NADH oxidase and ADH activities should be significantly higher, than the measured values. Implications of this finding for the directionality of the ADH isoenzyme operation in vivo and interactions between ADH and Pdc are discussed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05625-5.

Model-driven analysis of mutant fitness experiments improves genome-scale metabolic models of Zymomonas mobilis ZM4

Genome-scale metabolic models have been utilized extensively in the study and engineering of the organisms they describe. Here we present the analysis of a published dataset from pooled transposon mutant fitness experiments as an approach for improving the accuracy and gene-reaction associations of a metabolic model for Zymomonas mobilis ZM4, an industrially relevant ethanologenic organism with extremely high glycolytic flux and low biomass yield. Gene essentiality predictions made by the draft model were compared to data from individual pooled mutant experiments to identify areas of the model requiring deeper validation. Subsequent experiments showed that some of the discrepancies between the model and dataset were caused by polar effects, mis-mapped barcodes, or mutants carrying both wild-type and transposon disrupted gene copies-highlighting potential limitations inherent to data from individual mutants in these high-throughput datasets. Therefore, we analyzed correlations in fitness scores across all 492 experiments in the dataset in the context of functionally related metabolic reaction modules identified within the model via flux coupling analysis. These correlations were used to identify candidate genes for a reaction in histidine biosynthesis lacking an annotated gene and highlight metabolic modules with poorly correlated gene fitness scores. Additional genes for reactions involved in biotin, ubiquinone, and pyridoxine biosynthesis in Z. mobilis were identified and confirmed using mutant complementation experiments. These discovered genes, were incorporated into the final model, iZM4_478, which contains 747 metabolic and transport reactions (of which 612 have gene-protein-reaction associations), 478 genes, and 616 unique metabolites, making it one of the most complete models of Z. mobilis ZM4 to date. The methods of analysis that we applied here with the Z. mobilis transposon mutant dataset, could easily be utilized to improve future genome-scale metabolic reconstructions for organisms where these, or similar, high-throughput datasets are available.

Characterization and repurposing of the endogenous Type I-F CRISPR-Cas system of Zymomonas mobilis for genome engineering

Application of CRISPR-based technologies in non-model microorganisms is currently very limited. Here, we reported efficient genome engineering of an important industrial microorganism, Zymomonas mobilis, by repurposing the endogenous Type I-F CRISPR–Cas system upon its functional characterization. This toolkit included a series of genome engineering plasmids, each carrying an artificial self-targeting CRISPR and a donor DNA for the recovery of recombinants. Through this toolkit, various genome engineering purposes were efficiently achieved, including knockout of ZMO0038 (100% efficiency), cas2/3 (100%), and a genomic fragment of >10 kb (50%), replacement of cas2/3 with mCherry gene (100%), in situ nucleotide substitution (100%) and His-tagging of ZMO0038 (100%), and multiplex gene deletion (18.75%) upon optimal donor size determination. Additionally, the Type I-F system was further applied for CRISPRi upon Cas2/3 depletion, which has been demonstrated to successfully silence the chromosomally integrated mCherry gene with its fluorescence intensity reduced by up to 88%. Moreover, we demonstrated that genome engineering efficiency could be improved under a restriction–modification (R–M) deficient background, suggesting the perturbance of genome editing by other co-existing DNA targeting modules such as the R–M system. This study might shed light on exploiting and improving CRISPR–Cas systems in other microorganisms for genome editing and metabolic engineering practices.

In Vivo Thermodynamic Analysis of Glycolysis in Clostridium thermocellum and Thermoanaerobacterium saccharolyticum Using 13C and 2H Tracers

Thermodynamics constitutes a key determinant of flux and enzyme efficiency in metabolic networks. Here, we provide new insights into the divergent thermodynamics of the glycolytic pathways of C. thermocellum and T. saccharolyticum, two industrially relevant thermophilic bacteria whose metabolism still is not well understood. We report that while the glycolytic pathway in T. saccharolyticum is as thermodynamically favorable as that found in model organisms, such as E. coli or Saccharomyces cerevisiae, the glycolytic pathway of C. thermocellum operates near equilibrium. The use of a near-equilibrium glycolytic pathway, with potentially increased ATP yield, by this cellulolytic microbe may represent an evolutionary adaptation to growth on cellulose, but it has the drawback of being highly susceptible to product feedback inhibition. The results of this study will facilitate future engineering of high-performance strains capable of transforming cellulosic biomass to biofuels at high yields and titers.

Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are thermophilic anaerobic bacteria with complementary metabolic capabilities that utilize distinct glycolytic pathways for the conversion of cellulosic sugars to biofuels. We integrated quantitative metabolomics with ²H and ¹³C metabolic flux analysis to investigate the in vivo reversibility and thermodynamics of the central metabolic networks of these two microbes. We found that the glycolytic pathway in C. thermocellum operates remarkably close to thermodynamic equilibrium, with an overall drop in Gibbs free energy 5-fold lower than that of T. saccharolyticum or anaerobically grown Escherichia coli. The limited thermodynamic driving force of glycolysis in C. thermocellum could be attributed in large part to the small free energy of the phosphofructokinase reaction producing fructose bisphosphate. The ethanol fermentation pathway was also substantially more reversible in C. thermocellum than in T. saccharolyticum. These observations help explain the comparatively low ethanol titers of C. thermocellum and suggest engineering interventions that can be used to increase its ethanol productivity and glycolytic rate. In addition to thermodynamic analysis, we used our isotope tracer data to reconstruct the T. saccharolyticum central metabolic network, revealing exclusive use of the Embden-Meyerhof-Parnas (EMP) pathway for glycolysis, a bifurcated tricarboxylic acid (TCA) cycle, and a sedoheptulose bisphosphate bypass active within the pentose phosphate pathway.

IMPORTANCE Thermodynamics constitutes a key determinant of flux and enzyme efficiency in metabolic networks. Here, we provide new insights into the divergent thermodynamics of the glycolytic pathways of C. thermocellum and T. saccharolyticum, two industrially relevant thermophilic bacteria whose metabolism still is not well understood. We report that while the glycolytic pathway in T. saccharolyticum is as thermodynamically favorable as that found in model organisms, such as E. coli or Saccharomyces cerevisiae, the glycolytic pathway of C. thermocellum operates near equilibrium. The use of a near-equilibrium glycolytic pathway, with potentially increased ATP yield, by this cellulolytic microbe may represent an evolutionary adaptation to growth on cellulose, but it has the drawback of being highly susceptible to product feedback inhibition. The results of this study will facilitate future engineering of high-performance strains capable of transforming cellulosic biomass to biofuels at high yields and titers.

Metabolic Profiling of Glucose-Fed Metabolically Active Resting Zymomonas mobilis Strains

Zymomonas mobilis is the most efficient bacterial ethanol producer and its physiology is potentially applicable to industrial-scale bioethanol production. However, compared to other industrially important microorganisms, the Z. mobilis metabolome and adaptation to various nutritional and genetic perturbations have been poorly characterized. For rational metabolic engineering, it is essential to understand how central metabolism and intracellular redox balance are maintained in Z. mobilis under various conditions. In this study, we applied quantitative mass spectrometry-based metabolomics to explore how glucose-fed non-growing Z. mobilis Zm6 cells metabolically adapt to change of oxygen availability. Mutants partially impaired in ethanol synthesis (Zm6 adhB) or oxidative stress response (Zm6 cat) were also examined. Distinct patterns of adaptation of central metabolite pools due to the change in cultivation condition and between the mutants and Zm6 reference strain were observed. Decreased NADH/NAD ratio under aerobic incubation corresponded to higher concentrations of the phosphorylated glycolytic intermediates, in accordance with predictions of the kinetic model of Entner–Doudoroff pathway. The effects on the metabolite pools of aerobic to anaerobic transition were similar in the mutants, yet less pronounced. The present data on metabolic plasticity of non-growing Z. mobilis cells will facilitate the further metabolic engineering of the respective strains and their application as biocatalysts.

Zymomonas mobilis metabolism: Novel tools and targets for its rational engineering

Zymomonas mobilis is an α-proteobacterium that interests the biofuel industry due to its perfect ethanol fermentation yields. From its first description as a bacterial isolate in fermented alcoholic beverages to date, Z. mobilis has been rigorously studied in directions basic and applied. The Z. mobilis powerful Entner-Doudoroff glycolytic pathway has been the center of rigorous biochemical studies and, aside from ethanol, it has attracted interest in terms of high-added-value chemical manufacturing. Energetic balances and the effects of respiration have been explored in fundamental directions as also in applications pursuing strain enhancement and the utilization of alternative carbon sources. Metabolic modeling has addressed the optimization of the biochemical circuitry at various conditions of growth and/or substrate utilization; it has been also critical in predicting desirable end-product yields via flux redirection. Lastly, stress tolerance has received particular attention, since it directly determines biocatalytical performance at challenging bioreactor conditions. At a genetic level, advances in the genetic engineering of the organism have brought forth beneficial manipulations in the Z. mobilis gene pool, e.g., knock-outs, knock-ins and gene stacking, aiming to broaden the metabolic repertoire and increase robustness. Recent omic and expressional studies shed light on the genomic content of the most applied strains and reveal landscapes of activity manifested at ambient or reactor-based conditions. Studies such as those reviewed in this work, contribute to the understanding of the biology of Z. mobilis, enable insightful strain development, and pave the way for the transformation of Z. mobilis into a consummate organism for biomass conversion.

New technologies provide more metabolic engineering strategies for bioethanol production in Zymomonas mobilis

Bioethanol has been considered as a potentially renewable energy source, and metabolic engineering plays an important role in the production of biofuels. As an efficient ethanol-producing bacterium, Zymomonas mobilis has garnered special attention due to its high sugar uptake, ethanol yield, and tolerance. Different metabolic engineering strategies have been used to establish new metabolic pathways for Z. mobilis to broaden its substrate range, remove competing pathways, and enhance its tolerance to ethanol and lignocellulosic hydrolysate inhibitors. Recent advances in omics technology, computational modeling and simulation, system biology, and synthetic biology contribute to the efficient re-design and manipulation of microbes via metabolic engineering at the whole-cell level. In this review, we summarize the progress of some new technologies used for metabolic engineering to improve bioethanol production and tolerance in Z. mobilis. Some successful examples of metabolic engineering used to develop strains for ethanol production are described in detail. Lastly, some important strategies for future metabolic engineering efforts are also highlighted.

Model-based metabolism design: constraints for kinetic and stoichiometric models

The implementation of model-based designs in metabolic engineering and synthetic biology may fail. One of the reasons for this failure is that only a part of the real-world complexity is included in models. Still, some knowledge can be simplified and taken into account in the form of optimization constraints to improve the feasibility of model-based designs of metabolic pathways in organisms. Some constraints (mass balance, energy balance, and steady-state assumption) serve as a basis for many modelling approaches. There are others (total enzyme activity constraint and homeostatic constraint) proposed decades ago, but which are frequently ignored in design development. Several new approaches of cellular analysis have made possible the application of constraints like cell size, surface, and resource balance. Constraints for kinetic and stoichiometric models are grouped according to their applicability preconditions in (1) general constraints, (2) organism-level constraints, and (3) experiment-level constraints. General constraints are universal and are applicable for any system. Organism-level constraints are applicable for biological systems and usually are organism-specific, but these constraints can be applied without information about experimental conditions. To apply experimental-level constraints, peculiarities of the organism and the experimental set-up have to be taken into account to calculate the values of constraints. The limitations of applicability of particular constraints for kinetic and stoichiometric models are addressed.

Advances and prospects in metabolic engineering of Zymomonas mobilis

Biorefinery of biomass-based biofuels and biochemicals by microorganisms is a competitive alternative of traditional petroleum refineries. Zymomonas mobilis is a natural ethanologen with many desirable characteristics, which makes it an ideal industrial microbial biocatalyst for commercial production of desirable bioproducts through metabolic engineering. In this review, we summarize the metabolic engineering progress achieved in Z. mobilis to expand its substrate and product ranges as well as to enhance its robustness against stressful conditions such as inhibitory compounds within the lignocellulosic hydrolysates and slurries. We also discuss a few metabolic engineering strategies that can be applied in Z. mobilis to further develop it as a robust workhorse for economic lignocellulosic bioproducts. In addition, we briefly review the progress of metabolic engineering in Z. mobilis related to the classical synthetic biology cycle of "Design-Build-Test-Learn", as well as the progress and potential to develop Z. mobilis as a model chassis for biorefinery practices in the synthetic biology era.

Complete genome sequence and the expression pattern of plasmids of the model ethanologen Zymomonas mobilis ZM4 and its xylose-utilizing derivatives 8b and 2032

BACKGROUND

Zymomonas mobilis is a natural ethanologen being developed and deployed as an industrial biofuel producer. To date, eight Z. mobilis strains have been completely sequenced and found to contain 2-8 native plasmids. However, systematic verification of predicted Z. mobilis plasmid genes and their contribution to cell fitness has not been hitherto addressed. Moreover, the precise number and identities of plasmids in Z. mobilis model strain ZM4 have been unclear. The lack of functional information about plasmid genes in ZM4 impedes ongoing studies for this model biofuel-producing strain.

RESULTS

In this study, we determined the complete chromosome and plasmid sequences of ZM4 and its engineered xylose-utilizing derivatives 2032 and 8b. Compared to previously published and revised ZM4 chromosome sequences, the ZM4 chromosome sequence reported here contains 65 nucleotide sequence variations as well as a 2400-bp insertion. Four plasmids were identified in all three strains, with 150 plasmid genes predicted in strain ZM4 and 2032, and 153 plasmid genes predicted in strain 8b due to the insertion of heterologous DNA for expanded substrate utilization. Plasmid genes were then annotated using Blast2GO, InterProScan, and systems biology data analyses, and most genes were found to have apparent orthologs in other organisms or identifiable conserved domains. To verify plasmid gene prediction, RNA-Seq was used to map transcripts and also compare relative gene expression under various growth conditions, including anaerobic and aerobic conditions, or growth in different concentrations of biomass hydrolysates. Overall, plasmid genes were more responsive to varying hydrolysate concentrations than to oxygen availability. Additionally, our results indicated that although all plasmids were present in low copy number (about 1-2 per cell), the copy number of some plasmids varied under specific growth conditions or due to heterologous gene insertion.

CONCLUSIONS

The complete genome of ZM4 and two xylose-utilizing derivatives is reported in this study, with an emphasis on identifying and characterizing plasmid genes. Plasmid gene annotation, validation, expression levels at growth conditions of interest, and contribution to host fitness are reported for the first time.

Zymomonas mobilis as a model system for production of biofuels and biochemicals

Zymomonas mobilis is a natural ethanologen with many desirable industrial biocatalyst characteristics. In this review, we will discuss work to develop Z. mobilis as a model system for biofuel production from the perspectives of substrate utilization, development for industrial robustness, potential product spectrum, strain evaluation and fermentation strategies. This review also encompasses perspectives related to classical genetic tools and emerging technologies in this context.

The Low Energy-Coupling Respiration in Zymomonas mobilis Accelerates Flux in the Entner-Doudoroff Pathway

Performing oxidative phosphorylation is the primary role of respiratory chain both in bacteria and eukaryotes. Yet, the branched respiratory chains of prokaryotes contain alternative, low energy-coupling electron pathways, which serve for functions other than oxidative ATP generation (like those of respiratory protection, adaptation to low-oxygen media, redox balancing, etc.), some of which are still poorly understood. We here demonstrate that withdrawal of reducing equivalents by the energetically uncoupled respiratory chain of the bacterium Zymomonas mobilis accelerates its fermentative catabolism, increasing the glucose consumption rate. This is in contrast to what has been observed in other respiring bacteria and yeast. This effect takes place after air is introduced to glucose-consuming anaerobic cell suspension, and can be simulated using a kinetic model of the Entner-Doudoroff pathway in combination with a simple net reaction of NADH oxidation that does not involve oxidative phosphorylation. Although aeration hampers batch growth of respiring Z. mobilis culture due to accumulation of toxic byproducts, nevertheless under non-growing conditions respiration is shown to confer an adaptive advantage for the wild type over the non-respiring Ndh knock-out mutant. If cells get occasional access to limited amount of glucose for short periods of time, the elevated glucose uptake rate selectively improves survival of the respiring Z. mobilis phenotype.

Systematic development of biomass overproducing Scheffersomyces stipitis for high-cell-density fermentations

The development of economically feasible bio-based process requires efficient cell factories capable of producing the desired product at high titer under high-cell-density fermentation. Herein we present a combinatorial approach based on systems metabolic engineering and metabolic evolution for the development of efficient biomass-producing strain. Systems metabolic engineering guided by flux balance analysis (FBA) was first employed to rationally design mutant strains of Scheffersomyces stipitis with high biomass yield. By experimentally implementing these mutations, the biomass yield was improved by 30% in GPD1, 25% in TKL1, 30% in CIT1, and 44% in ZWF1 overexpressed mutants compared to wild-type. These designed mutants were further fine-tuned through metabolic evolution resulting in the maximal biomass yield of 0.49 g-cdw/g-glucose, which matches well with predicted yield phenotype. The constructed mutants are beneficial for biotechnology applications dealing with high cell titer cultivations. This work demonstrates a solid confirmation of systems metabolic engineering in combination with metabolic evolution approach for efficient strain development, which could assist in rapid optimization of cell factory for an economically viable and sustainable bio-based process.

Rites of passage: requirements and standards for building kinetic models of metabolic phenotypes

The overarching ambition of kinetic metabolic modeling is to capture the dynamic behavior of metabolism to such an extent that systems and synthetic biology strategies can reliably be tested in silico. The lack of kinetic data hampers the development of kinetic models, and most of the current models use ad hoc reduced stoichiometry or oversimplified kinetic rate expressions, which may limit their predictive strength. There is a need to introduce the community-level standards that will organize and accelerate the future developments in this area. We introduce here a set of requirements that will ensure the model quality, we examine the current kinetic models with respect to these requirements, and we propose a general workflow for constructing models that satisfy these requirements.

Improving prediction fidelity of cellular metabolism with kinetic descriptions

Several modeling frameworks for describing and redirecting cellular metabolism have been developed keeping pace with the rapid development in high-throughput data generation and advances in metabolic engineering techniques. The incorporation of kinetic information within stoichiometry-only modeling techniques offers potential advantages for improved phenotype prediction and consequently more precise computational strain design. In addition to substrate-level kinetic regulatory information, the integration of a number of additional layers of regulation at the transcription, translation, and post-translation levels is sought after by many research groups. However, the practical integration of these complex biological processes into a unified framework amenable to design remains an ongoing challenge.

CceR and AkgR regulate central carbon and energy metabolism in alphaproteobacteria

Many pathways of carbon and energy metabolism are conserved across the phylogeny, but the networks that regulate their expression or activity often vary considerably among organisms. In this work, we show that two previously uncharacterized transcription factors (TFs) are direct regulators of genes encoding enzymes of central carbon and energy metabolism in the alphaproteobacterium Rhodobacter sphaeroides. The LacI family member CceR (RSP_1663) directly represses genes encoding enzymes in the Entner-Doudoroff pathway, while activating those encoding the F1F0 ATPase and enzymes of the tricarboxylic acid (TCA) cycle and gluconeogenesis, providing a direct transcriptional network connection between carbon and energy metabolism. We identified bases that are important for CceR DNA binding and showed that DNA binding by this TF is inhibited by 6-phosphogluconate. We also showed that the GntR family TF AkgR (RSP_0981) directly activates genes encoding several TCA cycle enzymes, and we identified conditions where its activity is increased. The properties of single and double ΔCceR and ΔAkgR mutants illustrate that these 2 TFs cooperatively regulate carbon and energy metabolism. Comparative genomic analysis indicates that CceR and AkgR orthologs are found in other alphaproteobacteria, where they are predicted to have a conserved function in regulating central carbon metabolism. Our characterization of CceR and AkgR has provided important new insight into the networks that control central carbon and energy metabolism in alphaproteobacteria that can be exploited to modify or engineer new traits in these widespread and versatile bacteria.

IMPORTANCE

To extract and conserve energy from nutrients, cells coordinate a set of metabolic pathways into integrated networks. Many pathways that conserve energy or interconvert metabolites are conserved across cells, but the networks regulating these processes are often highly variable. In this study, we characterize two previously unknown transcriptional regulators of carbon and energy metabolism that are conserved in alphaproteobacteria, a group of abundant, environmentally and biotechnologically important organisms. We identify the genes they regulate, the DNA sequences they recognize, the metabolite that controls the activity of one of the regulators, and conditions where they are required for growth. We provide important new insight into conserved cellular networks that can also be used to improve a variety of hosts for converting feedstock into valuable products.

Elucidation of Zymomonas mobilis physiology and stress responses by quantitative proteomics and transcriptomics

Zymomonas mobilis is an excellent ethanologenic bacterium. Biomass pretreatment and saccharification provides access to simple sugars, but also produces inhibitors such as acetate and furfural. Our previous work has identified and confirmed the genetic change of a 1.5-kb deletion in the sodium acetate tolerant Z. mobilis mutant (AcR) leading to constitutively elevated expression of a sodium proton antiporter encoding gene nhaA, which contributes to the sodium acetate tolerance of AcR mutant. In this study, we further investigated the responses of AcR and wild-type ZM4 to sodium acetate stress in minimum media using both transcriptomics and a metabolic labeling approach for quantitative proteomics the first time. Proteomic measurements at two time points identified about eight hundreds proteins, or about half of the predicted proteome. Extracellular metabolite analysis indicated AcR overcame the acetate stress quicker than ZM4 with a concomitant earlier ethanol production in AcR mutant, although the final ethanol yields and cell densities were similar between two strains. Transcriptomic samples were analyzed for four time points and revealed that the response of Z. mobilis to sodium acetate stress is dynamic, complex, and involved about one-fifth of the total predicted genes from all different functional categories. The modest correlations between proteomic and transcriptomic data may suggest the involvement of posttranscriptional control. In addition, the transcriptomic data of forty-four microarrays from four experiments for ZM4 and AcR under different conditions were combined to identify strain-specific, media-responsive, growth phase-dependent, and treatment-responsive gene expression profiles. Together this study indicates that minimal medium has the most dramatic effect on gene expression compared to rich medium followed by growth phase, inhibitor, and strain background. Genes involved in protein biosynthesis, glycolysis and fermentation as well as ATP synthesis and stress response play key roles in Z. mobilis metabolism with consistently strong expression levels under different conditions.

The inefficient aerobic energetics of Zymomonas mobilis: identifying the bottleneck

To investigate the mechanisms of Zymomonas mobilis uncoupled aerobic metabolism, growth properties of the wild-type strain Zm6 were compared to those of its respiratory mutants cytB and cydB, and the effects of the ATPase inhibitor DCCD on growth and intracellular ATP concentration were studied. The effects of the ATPase inhibitor DCCD on growth and intracellular ATP concentration strongly indicated that the apparent lack of oxidative phosphorylation in aerobically growing Z. mobilis culture might be caused by the ATP hydrolyzing activity of the H(+) -dependent ATPase in all analyzed strains. Aerobic growth yields of the mutants, and their capacity of oxidative ATP synthesis with ethanol were closely similar, not supporting presence of one major, yet energetically inefficient electron transport branch causing the observed poor aerobic growth and lack of oxidative phosphorylation in Z. mobilis. We suggest that rapidly operating Entner-Doudoroff pathway generates too high phosphorylation potential for the weakly coupled respiratory system to shift the H(+) -dependent ATPase toward ATP synthesis.

Modeling of Zymomonas mobilis central metabolism for novel metabolic engineering strategies

Mathematical modeling of metabolism is essential for rational metabolic engineering. The present work focuses on several types of modeling approach to quantitative understanding of central metabolic network and energetics in the bioethanol-producing bacterium Zymomonas mobilis. Combined use of Flux Balance, Elementary Flux Mode, and thermodynamic analysis of its central metabolism, together with dynamic modeling of the core catabolic pathways, can help to design novel substrate and product pathways by systematically analyzing the solution space for metabolic engineering, and yields insights into the function of metabolic network, hardly achievable without applying modeling tools.