The heterologous production of terpenes by the thermophile Parageobacillus thermoglucosidasius in a consolidated bioprocess using waste bread

Parageobacillus thermoglucosidasius is a genetically tractable thermophile that grows rapidly at elevated temperatures, with a doubling time at 65 °C comparable to the shortest doubling times of Escherichia coli. It is capable of using a wide variety of substrates, including carbohydrate oligomers, and has been developed for the industrial production of ethanol. In this study, P. thermoglucosidasius NCIMB11955 has been engineered to produce the sesquiterpene τ-muurolol by introduction of a heterologous mevalonate pathway constructed using genes from several thermophilic archaea together with a recently characterised thermostable τ-muurolol synthase. P. thermoglucosidasius naturally uses the methylerythritol phosphate pathway for production of the terpene precursor, isopentenyl pyrophosphate, while archaea use a version of the mevalonate pathway. By introducing the orthogonal archaeal pathway it was possible to increase the flux through to sesquiterpene biosynthesis. Construction of such a large metabolic pathway created problems with genetic vector introduction and stability, so recombinant plasmids were introduced by conjugation, and a thermostable serine integrase system was developed for integration of large pathways onto the chromosome. Finally, by making the heterologous pathway maltose-inducible we demonstrate that the new strain is capable of using waste bread directly as an autoinduction carbon source for the production of terpenes in a consolidated bioprocess.     

Constraints-based genome-scale metabolic simulation for systems metabolic engineering

Random mutagenesis and selection approaches used traditionally for the development of industrial strains have largely been complemented by metabolic engineering, which allows purposeful modification of metabolic and cellular characteristics by using recombinant DNA and other molecular biological techniques. As systems biology advances as a new paradigm of research thanks to the development of genome-scale computational tools and high-throughput experimental technologies including omics, systems metabolic engineering allowing modification of metabolic, regulatory and signaling networks of the cell at the systems-level is becoming possible. In silico genome-scale metabolic model and its simulation play increasingly important role in providing systematic strategies for metabolic engineering. The in silico genome-scale metabolic model is developed using genomic annotation, metabolic reactions, literature information, and experimental data. The advent of in silico genome-scale metabolic model brought about the development of various algorithms to simulate the metabolic status of the cell as a whole. In this paper, we review the algorithms developed for the system-wide simulation and perturbation of cellular metabolism, discuss the characteristics of these algorithms, and suggest future research direction.     

Semi-automated Curation of Metabolic Models via Flux Balance Analysis: A Case Study with Mycoplasma gallisepticum

Primarily used for metabolic engineering and synthetic biology, genome-scale metabolic modeling shows tremendous potential as a tool for fundamental research and curation of metabolism. Through a novel integration of flux balance analysis and genetic algorithms, a strategy to curate metabolic networks and facilitate identification of metabolic pathways that may not be directly inferable solely from genome annotation was developed. Specifically, metabolites involved in unknown reactions can be determined, and potentially erroneous pathways can be identified. The procedure developed allows for new fundamental insight into metabolism, as well as acting as a semi-automated curation methodology for genome-scale metabolic modeling. To validate the methodology, a genome-scale metabolic model for the bacterium Mycoplasma gallisepticum was created. Several reactions not predicted by the genome annotation were postulated and validated via the literature. The model predicted an average growth rate of 0.358±0.12, closely matching the experimentally determined growth rate of M. gallisepticum of 0.244±0.03. This work presents a powerful algorithm for facilitating the identification and curation of previously known and new metabolic pathways, as well as presenting the first genome-scale reconstruction of M. gallisepticum.     

In silico metabolic engineering of Clostridium ljungdahlii for synthesis gas fermentation

Synthesis gas fermentation is one of the most promising routes to convert synthesis gas (syngas; mainly comprised of H₂ and CO) to renewable liquid fuels and chemicals by specialized bacteria. The most commonly studied syngas fermenting bacterium is Clostridium ljungdahlii, which produces acetate and ethanol as its primary metabolic byproducts. Engineering of C. ljungdahlii metabolism to overproduce ethanol, enhance the synthesize of the native byproducts lactate and 2,3-butanediol, and introduce the synthesis of non-native products such as butanol and butyrate has substantial commercial value. We performed in silico metabolic engineering studies using a genome-scale reconstruction of C. ljungdahlii metabolism and the OptKnock computational framework to identify gene knockouts that were predicted to enhance the synthesis of these native products and non-native products, introduced through insertion of the necessary heterologous pathways. The OptKnock derived strategies were often difficult to assess because increase product synthesis was invariably accompanied by decreased growth. Therefore, the OptKnock strategies were further evaluated using a spatiotemporal metabolic model of a syngas bubble column reactor, a popular technology for large-scale gas fermentation. Unlike flux balance analysis, the bubble column model accounted for the complex tradeoffs between increased product synthesis and reduced growth rates of engineered mutants within the spatially varying column environment. The two-stage methodology for deriving and evaluating metabolic engineering strategies was shown to yield new C. ljungdahlii gene targets that offer the potential for increased product synthesis under realistic syngas fermentation conditions.     

In silico genome-scale metabolic analysis of Pseudomonas putida KT2440 for polyhydroxyalkanoate synthesis,degradation of aromatics and anaerobic survival

Genome-scale metabolic models have been appearing with increasing frequency and have been employed in a wide range of biotechnological applications as well as in biological studies. With the metabolic model as a platform, engineering strategies have become more systematic and focused, unlike the random shotgun approach used in the past. Here we present the genome-scale metabolic model of the versatile Gram-negative bacterium Pseudomonas putida, which has gained widespread interest for various biotechnological applications. With the construction of the genome-scale metabolic model of P. putida KT2440, PpuMBEL1071, we investigated various characteristics of P. putida, such as its capacity for synthesizing polyhydroxyalkanoates (PHA) and degrading aromatics. Although P. putida has been characterized as a strict aerobic bacterium, the physiological characteristics required to achieve anaerobic survival were investigated. Through analysis of PpuMBEL1071, extended survival of P. putida under anaerobic stress was achieved by introducing the ackA gene from Pseudomonas aeruginosa and Escherichia coli.     

Engineering Bacteroides thetaiotaomicron to produce non-native butyrate based on a genome-scale metabolic model-guided design

Bacteroides thetaiotaomicron represents a major symbiont of the human gut microbiome that is increasingly viewed as a promising candidate strain for microbial therapeutics. Here, we engineer B. thetaiotaomicron for heterologous production of non-native butyrate as a proof-of-concept biochemical at therapeutically relevant concentrations. Since B. thetaiotaomicron is not a natural producer of butyrate, we heterologously expressed a butyrate biosynthetic pathway in the strain, which led to the production of butyrate at the final concentration of 12 mg/L in a rich medium. Further optimization of butyrate production was achieved by a round of metabolic engineering guided by an expanded genome-scale metabolic model (GEM) of B. thetaiotaomicron. The in silico knock-out simulation of the expanded model showed that pta and ldhD were the potent knock-out targets to enhance butyrate production. The maximum titer and specific productivity of butyrate in the pta-ldhD double knockout mutant increased by nearly 3.4 and 4.8 folds, respectively. To our knowledge, this is the first engineering attempt that enabled butyrate production from a non-butyrate producing commensal B. thetaiotaomicron. The study also highlights that B. thetaiotaomicron can serve as an effective strain for live microbial therapeutics in human.     

Flux analysis and metabolomics for systematic metabolic engineering of microorganisms

Rational engineering of metabolism is important for bio-production using microorganisms. Metabolic design based on in silico simulations and experimental validation of the metabolic state in the engineered strain helps in accomplishing systematic metabolic engineering. Flux balance analysis (FBA) is a method for the prediction of metabolic phenotype, and many applications have been developed using FBA to design metabolic networks. Elementary mode analysis (EMA) and ensemble modeling techniques are also useful tools for in silico strain design. The metabolome and flux distribution of the metabolic pathways enable us to evaluate the metabolic state and provide useful clues to improve target productivity. Here, we reviewed several computational applications for metabolic engineering by using genome-scale metabolic models of microorganisms. We also discussed the recent progress made in the field of metabolomics and ¹³C-metabolic flux analysis techniques, and reviewed these applications pertaining to bio-production development. Because these in silico or experimental approaches have their respective advantages and disadvantages, the combined usage of these methods is complementary and effective for metabolic engineering.     

Metabolic impact assessment for heterologous protein production in Streptomyces lividans based on genome-scale metabolic network modeling

The metabolic impact exerted on a microorganism due to heterologous protein production is still poorly understood in Streptomyces lividans. In this present paper, based on exometabolomic data, a proposed genome-scale metabolic network model is used to assess this metabolic impact in S. lividans. Constraint-based modeling results obtained in this work revealed that the metabolic impact due to heterologous protein production is widely distributed in the genome of S. lividans, causing both slow substrate assimilation and a shift in active pathways. Exchange fluxes that are critical for model performance have been identified for metabolites of mouse tumor necrosis factor, histidine, valine and lysine, as well as biomass. Our results unravel the interaction of heterologous protein production with intracellular metabolism of S. lividans, thus, a possible basis for further studies in relieving the metabolic burden via metabolic or bioprocess engineering.     

乳酸乳球菌NZ9000基因组规模代谢网络模型的构建与验证

随着后基因组时代的到来,工业微生物的代谢工程改造在工业生产上发挥着越来越重要的作用。而基因组规模代谢网络模型(Genome-scalemetabolicmodel,GSMM)将生物体体内所有已知代谢信息进行整合,为全局理解生物体的代谢状态、理性指导代谢工程改造提供了最佳的平台。乳酸乳球菌NZ9000(Lactococcuslactis NZ9000)作为工业发酵领域的重要菌株之一,由于其遗传背景清晰且几乎不分泌蛋白,是基因工程改造和外源蛋白表达的理想模式菌株。文中基于基因组功能注释和比较基因组学构建了L.lactisNZ9000的首个基因组规模代谢网络模型iWK557,包含557个基因、668个代谢物、840个反应,并进一步在定性和定量两个层次验证了iWK557的准确性,以期为理性指导L. lactis NZ9000代谢工程改造提供良好工具。     

Genome-scale metabolic model for Lactococcus lactis MG1363 and its application to the analysis of flavor formation

Lactococcus lactis subsp. cremoris MG1363 is a paradigm strain for lactococci used in industrial dairy fermentations. However, despite of its importance for process development, no genome-scale metabolic model has been reported thus far. Moreover, current models for other lactococci only focus on growth and sugar degradation. A metabolic model that includes nitrogen metabolism and flavor-forming pathways is instrumental for the understanding and designing new industrial applications of these lactic acid bacteria. A genome-scale, constraint-based model of the metabolism and transport in L. lactis MG1363, accounting for 518 genes, 754 reactions, and 650 metabolites, was developed and experimentally validated. Fifty-nine reactions are directly or indirectly involved in flavor formation. Flux Balance Analysis and Flux Variability Analysis were used to investigate flux distributions within the whole metabolic network. Anaerobic carbon-limited continuous cultures were used for estimating the energetic parameters. A thorough model-driven analysis showing a highly flexible nitrogen metabolism, e.g., branched-chain amino acid catabolism which coupled with the redox balance, is pivotal for the prediction of the formation of different flavor compounds. Furthermore, the model predicted the formation of volatile sulfur compounds as a result of the fermentation. These products were subsequently identified in the experimental fermentations carried out. Thus, the genome-scale metabolic model couples the carbon and nitrogen metabolism in L. lactis MG1363 with complete known catabolic pathways leading to flavor formation. The model provided valuable insights into the metabolic networks underlying flavor formation and has the potential to contribute to new developments in dairy industries and cheese-flavor research.     

Flux balance analysis of the ammonia-oxidizing bacterium Nitrosomonas europaea ATCC19718 unravels specific metabolic activities while degrading toxic compounds

The ammonia-oxidizing bacterium Nitrosomonas europaea has been widely recognized as an important player in the nitrogen cycle as well as one of the most abundant members in microbial communities for the treatment of industrial or sewage wastewater. Its natural metabolic versatility and extraordinary ability to degrade environmental pollutants (e.g., aromatic hydrocarbons such as benzene and toluene) enable it to thrive under various harsh environmental conditions. Constraint-based metabolic models constructed from genome sequences enable quantitative insight into the central and specialized metabolism within a target organism. These genome-scale models have been utilized to understand, optimize, and design new strategies for improved bioprocesses. Reduced modeling approaches have been used to elucidate Nitrosomonas europaea metabolism at a pathway level. However, genome-scale knowledge about the simultaneous oxidation of ammonia and pollutant metabolism of N. europaea remains limited. Here, we describe the reconstruction, manual curation, and validation of the genome-scale metabolic model for N. europaea, iGC535. This reconstruction is the most accurate metabolic model for a nitrifying organism to date, reaching an average prediction accuracy of over 90% under several growth conditions. The manually curated model can predict phenotypes under chemolithotrophic and chemolithoorganotrophic conditions while oxidating methane and wastewater pollutants. Calculated flux distributions under different trophic conditions show that several key pathways are affected by the type of carbon source available, including central carbon metabolism and energy production.     

Model-aided atpE gene knockout strategy in Escherichia coli for enhanced succinic acid production from glycerol

Succinic acid is an important platform chemical with a variety of applications. Model-guided metabolic engineering strategies in Escherichia coli for strain improvement to increase succinic acid production using glucose and glycerol remain largely unexplored. Herein, we report what are, to our knowledge, the first metabolic knockout of the atpE gene to have increased succinic acid production using both glucose and alternative glycerol carbon sources in E. coli. Guided by a genome-scale metabolic model, we engineered the E. coli host to enhance anaerobic production of succinic acid by deleting the atpE gene, thereby generating additional reducing equivalents by blocking H⁺ conduction across the mutant cell membrane. This strategy produced 1.58 and .49 g l^?1 of succinic acid from glycerol and glucose substrate, respectively. This work further elucidates a model-guided and/or system-based metabolic engineering, involving only a single-gene deletion strategy for enhanced succinic acid production in E. coli.     

Comprehensive analysis of glucose and xylose metabolism in Escherichia coli under aerobic and anaerobic conditions by 13C metabolic flux analysis

Glucose and xylose are the two most abundant sugars derived from the breakdown of lignocellulosic biomass. While aerobic glucose metabolism is relatively well understood in E. coli, until now there have been only a handful of studies focused on anaerobic glucose metabolism and no ¹³C-flux studies on xylose metabolism. In the absence of experimentally validated flux maps, constraint-based approaches such as MOMA and RELATCH cannot be used to guide new metabolic engineering designs. In this work, we have addressed this critical gap in current understanding by performing comprehensive characterizations of glucose and xylose metabolism under aerobic and anaerobic conditions, using recent state-of-the-art techniques in ¹³C metabolic flux analysis (¹³C-MFA). Specifically, we quantified precise metabolic fluxes for each condition by performing parallel labeling experiments and analyzing the data through integrated ¹³C-MFA using the optimal tracers [1,2-¹³C]glucose, [1,6-¹³C]glucose, [1,2-¹³C]xylose and [5-¹³C]xylose. We also quantified changes in biomass composition and confirmed turnover of macromolecules by applying [U-¹³C]glucose and [U-¹³C]xylose tracers. We demonstrated that under anaerobic growth conditions there is significant turnover of lipids and that a significant portion of CO₂ originates from biomass turnover. Using knockout strains, we also demonstrated that β-oxidation is critical for anaerobic growth on xylose. Quantitative analysis of co-factor balances (NADH/FADH₂, NADPH, and ATP) for different growth conditions provided new insights regarding the interplay of energy and redox metabolism and the impact on E. coli cell physiology.     

Pan-cancer analysis of the metabolic reaction network

Metabolic reprogramming is considered a hallmark of malignant transformation. However, it is not clear whether the network of metabolic reactions expressed by cancers of different origin differ from each other or from normal human tissues. In this study, we reconstructed functional and connected genome-scale metabolic models for 917 primary tumor samples across 13 types based on the probability of expression for 3765 reference metabolic genes in the sample. This network-centric approach revealed that tumor metabolic networks are largely similar in terms of accounted reactions, despite diversity in the expression of the associated genes. On average, each network contained 4721 reactions, of which 74% were core reactions (present in >95% of all models). Whilst 99.3% of the core reactions were classified as housekeeping also in normal tissues, we identified reactions catalyzed by ARG2, RHAG, SLC6 and SLC16 family gene members, and PTGS1 and PTGS2 as core exclusively in cancer. These findings were subsequently replicated in an independent validation set of 3388 genome-scale metabolic models. The remaining 26% of the reactions were contextual reactions. Their inclusion was dependent in one case (GLS2) on the absence of TP53 mutations and in 94.6% of cases on differences in cancer types. This dependency largely resembled differences in expression patterns in the corresponding normal tissues, with some exceptions like the presence of the NANP-encoded reaction in tumors not from the female reproductive system or of the SLC5A9-encoded reaction in kidney-pancreatic-colorectal tumors. In conclusion, tumors expressed a metabolic network virtually overlapping the matched normal tissues, raising the possibility that metabolic reprogramming simply reflects cancer cell plasticity to adapt to varying conditions thanks to redundancy and complexity of the underlying metabolic networks. At the same time, the here uncovered exceptions represent a resource to identify selective liabilities of tumor metabolism.     

基因组尺度代谢网络研究进展

基因组尺度代谢网络从基因组序列出发,结合基因、蛋白质、代谢数据库和实验数据,从系统的角度定量研究生命体的代谢过程,了解各个组分之间的相互作用关系。这类网络模型对于生命活动理论研究和优良工程菌的构建都具有重要的理论和实践意义。以下结合作者的实际研究经验,对基因组尺度代谢网络从重构到模拟直至应用进行了较为详细的介绍,并讨论了一些目前存在的难题和未来的研究方向。     

Genome-scale metabolic models as platforms for strain design and biological discovery

Genome-scale metabolic models (GEMs) have been developed and used in guiding systems’ metabolic engineering strategies for strain design and development. This strategy has been used in fermentative production of bio-based industrial chemicals and fuels from alternative carbon sources. However, computer-aided hypotheses building using established algorithms and software platforms for biological discovery can be integrated into the pipeline for strain design strategy to create superior strains of microorganisms for targeted biosynthetic goals. Here, I described an integrated workflow strategy using GEMs for strain design and biological discovery. Specific case studies of strain design and biological discovery using Escherichia coli genome-scale model are presented and discussed. The integrated workflow presented herein, when applied carefully would help guide future design strategies for high-performance microbial strains that have existing and forthcoming genome-scale metabolic models.     

Extension of the yeast metabolic model to include iron metabolism and its use to estimate global levels of iron-recruiting enzyme abundance from cofactor requirements

Metabolic networks adapt to changes in their environment by modulating the activity of their enzymes and transporters; often by changing their abundance. Understanding such quantitative changes can shed light onto how metabolic adaptation works, or how it can fail and lead to a metabolically dysfunctional state. We propose a strategy to quantify metabolic protein requirements for cofactor-utilising enzymes and transporters through constraint-based modelling. The first eukaryotic genome-scale metabolic model to comprehensively represent iron metabolism was constructed, extending the most recent community model of the Saccharomyces cerevisiae metabolic network. Partial functional impairment of the genes involved in the maturation of iron-sulphur (Fe-S) proteins was investigated employing the model and the in silico analysis revealed extensive rewiring of the fluxes in response to this functional impairment, despite its marginal phenotypic effect. The optimal turnover rate of enzymes bearing ion cofactors can be determined via this novel approach; yeast metabolism, at steady state, was determined to employ a constant turnover of its iron-recruiting enzyme at a rate of 3.02 × 10 ⁻¹¹ mmol·(g biomass) ⁻¹·h ⁻¹.     

Genomic and metabolic analysis of Komagataeibacter xylinus DSM 2325 producing bacterial cellulose nanofiber

Bacterial cellulose nanofiber (CNF) is a polymer with a wide range of potential industrial applications. Several Komagataeibacter species, including Komagataeibacter xylinus as a model organism, produce CNF. However, the industrial application of CNF has been hampered by inefficient CNF production, necessitating metabolic engineering for the enhanced CNF production. Here, we present complete genome sequence and a genome-scale metabolic model KxyMBEL1810 of K. xylinus DSM 2325 for metabolic engineering applications. Genome analysis of this bacterium revealed that a set of genes associated with CNF biosynthesis and regulation were present in this bacterium, which were also conserved in another six representative Komagataeibacter species having complete genome information. To better understand the metabolic characteristics of K. xylinus DSM 2325, KxyMBEL1810 was reconstructed using genome annotation data, relevant computational resources and experimental growth data generated in this study. Random sampling and correlation analysis of the KxyMBEL1810 predicted pgi and gnd genes as novel overexpression targets for the enhanced CNF production. Among engineered K. xylinus strains individually overexpressing heterologous pgi and gnd genes, either from Escherichia coli or Corynebacterium glutamicum, batch fermentation of a strain overexpressing the E. coli pgi gene produced 3.15 g/L of CNF in a complex medium containing glucose, which was the best CNF concentration achieved in this study, and 115.8% higher than that (1.46 g/L) obtained from the control strain. Genome sequence data and KxyMBEL1810 generated in this study should be useful resources for metabolic engineering of K. xylinus for the enhanced CNF production.