Jong Myoung Park

Systems metabolic engineering of microorganisms for natural and non-natural chemicals

Growing concerns over limited fossil resources and associated environmental problems are motivating the development of sustainable processes for the production of chemicals, fuels and materials from renewable resources. Metabolic engineering is a key enabling technology for transforming microorganisms into efficient cell factories for these compounds. Systems metabolic engineering, which incorporates the concepts and techniques of systems biology, synthetic biology and evolutionary engineering at the systems level, offers a conceptual and technological framework to speed the creation of new metabolic enzymes and pathways or the modification of existing pathways for the optimal production of desired products. Here we discuss the general strategies of systems metabolic engineering and examples of its application and offer insights as to when and how each of the different strategies should be used. Finally, we highlight the limitations and challenges to be overcome for the systems metabolic engineering of microorganisms at more advanced levels.

Engineering of microorganisms for the production of biofuels and perspectives based on systems metabolic engineering approaches

The increasing oil price and environmental concerns caused by the use of fossil fuel have renewed our interest in utilizing biomass as a sustainable resource for the production of biofuel. It is however essential to develop high performance microbes that are capable of producing biofuels with very high efficiency in order to compete with the fossil fuel. Recently, the strategies for developing microbial strains by systems metabolic engineering, which can be considered as metabolic engineering integrated with systems biology and synthetic biology, have been developed. Systems metabolic engineering allows successful development of microbes that are capable of producing several different biofuels including bioethanol, biobutanol, alkane, and biodiesel, and even hydrogen. In this review, the approaches employed to develop efficient biofuel producers by metabolic engineering and systems metabolic engineering approaches are reviewed with relevant example cases. It is expected that systems metabolic engineering will be employed as an essential strategy for the development of microbial strains for industrial applications.

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.

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.

The genome-scale metabolic network analysis of Zymomonas mobilis ZM4 explains physiological features and suggests ethanol and succinic acid production strategies.

BackgroundZymomonas mobilis ZM4 is a Gram-negative bacterium that can efficiently produce ethanol from various carbon substrates, including glucose, fructose, and sucrose, via the Entner-Doudoroff pathway. However, systems metabolic engineering is required to further enhance its metabolic performance for industrial application. As an important step towards this goal, the genome-scale metabolic model of Z. mobilis is required to systematically analyze in silico the metabolic characteristics of this bacterium under a wide range of genotypic and environmental conditions.ResultsThe genome-scale metabolic model of Z. mobilis ZM4, ZmoMBEL601, was reconstructed based on its annotated genes, literature, physiological and biochemical databases. The metabolic model comprises 579 metabolites and 601 metabolic reactions (571 biochemical conversion and 30 transport reactions), built upon extensive search of existing knowledge. Physiological features of Z. mobilis were then examined using constraints-based flux analysis in detail as follows. First, the physiological changes of Z. mobilis as it shifts from anaerobic to aerobic environments (i.e. aerobic shift) were investigated. Then the intensities of flux-sum, which is the cluster of either all ingoing or outgoing fluxes through a metabolite, and the maximum in silico yields of ethanol for Z. mobilis and Escherichia coli were compared and analyzed. Furthermore, the substrate utilization range of Z. mobilis was expanded to include pentose sugar metabolism by introducing metabolic pathways to allow Z. mobilis to utilize pentose sugars. Finally, double gene knock-out simulations were performed to design a strategy for efficiently producing succinic acid as another example of application of the genome-scale metabolic model of Z. mobilis.ConclusionThe genome-scale metabolic model reconstructed in this study was able to successfully represent the metabolic characteristics of Z. mobilis under various conditions as validated by experiments and literature information. This reconstructed metabolic model will allow better understanding of Z. mobilis metabolism and consequently designing metabolic engineering strategies for various biotechnological applications.

Prediction of metabolic fluxes by incorporating genomic context and flux-converging pattern analyses

Flux balance analysis (FBA) of a genome-scale metabolic model allows calculation of intracellular fluxes by optimizing an objective function, such as maximization of cell growth, under given constraints, and has found numerous applications in the field of systems biology and biotechnology. Due to the underdetermined nature of the system, however, it has limitations such as inaccurate prediction of fluxes and existence of multiple solutions for an optimal objective value. Here, we report a strategy for accurate prediction of metabolic fluxes by FBA combined with systematic and condition-independent constraints that restrict the achievable flux ranges of grouped reactions by genomic context and flux-converging pattern analyses. Analyses of three types of genomic contexts, conserved genomic neighborhood, gene fusion events, and co-occurrence of genes across multiple organisms, were performed to suggest a group of fluxes that are likely on or off simultaneously. The flux ranges of these grouped reactions were constrained by flux-converging pattern analysis. FBA of the Escherichia coli genome-scale metabolic model was carried out under several different genotypic (pykF, zwf, ppc, and sucA mutants) and environmental (altered carbon source) conditions by applying these constraints, which resulted in flux values that were in good agreement with the experimentally measured 13C-based fluxes. Thus, this strategy will be useful for accurately predicting the intracellular fluxes of large metabolic networks when their experimental determination is difficult.

Flux variability scanning based on enforced objective flux for identifying gene amplification targets.

BackgroundIn order to reduce time and efforts to develop microbial strains with better capability of producing desired bioproducts, genome-scale metabolic simulations have proven useful in identifying gene knockout and amplification targets. Constraints-based flux analysis has successfully been employed for such simulation, but is limited in its ability to properly describe the complex nature of biological systems. Gene knockout simulations are relatively straightforward to implement, simply by constraining the flux values of the target reaction to zero, but the identification of reliable gene amplification targets is rather difficult. Here, we report a new algorithm which incorporates physiological data into a model to improve the model’s prediction capabilities and to capitalize on the relationships between genes and metabolic fluxes.ResultsWe developed an algorithm, flux variability scanning based on enforced objective flux (FVSEOF) with grouping reaction (GR) constraints, in an effort to identify gene amplification targets by considering reactions that co-carry flux values based on physiological omics data via “GR constraints”. This method scans changes in the variabilities of metabolic fluxes in response to an artificially enforced objective flux of product formation. The gene amplification targets predicted using this method were validated by comparing the predicted effects with the previous experimental results obtained for the production of shikimic acid and putrescine in Escherichia coli. Moreover, new gene amplification targets for further enhancing putrescine production were validated through experiments involving the overexpression of each identified targeted gene under condition-controlled batch cultivation.ConclusionsFVSEOF with GR constraints allows identification of gene amplification targets for metabolic engineering of microbial strains in order to enhance the production of desired bioproducts. The algorithm was validated through the experiments on the enhanced production of putrescine in E. coli, in addition to the comparison with the previously reported experimental data. The FVSEOF strategy with GR constraints will be generally useful for developing industrially important microbial strains having enhanced capabilities of producing chemicals of interest.

Genome-scale reconstruction and in silico analysis of the Ralstonia eutropha H16 for polyhydroxyalkanoate synthesis, lithoautotrophic growth, and 2-methyl citric acid production

BackgroundRalstonia eutropha H16, found in both soil and water, is a Gram-negative lithoautotrophic bacterium that can utillize CO2 and H2 as its sources of carbon and energy in the absence of organic substrates. R. eutropha H16 can reach high cell densities either under lithoautotrophic or heterotrophic conditions, which makes it suitable for a number of biotechnological applications. It is the best known and most promising producer of polyhydroxyalkanoates (PHAs) from various carbon substrates and is an environmentally important bacterium that can degrade aromatic compounds. In order to make R. eutropha H16 a more efficient and robust biofactory, system-wide metabolic engineering to improve its metabolic performance is essential. Thus, it is necessary to analyze its metabolic characteristics systematically and optimize the entire metabolic network at systems level.ResultsWe present the lithoautotrophic genome-scale metabolic model of R. eutropha H16 based on the annotated genome with biochemical and physiological information. The stoichiometic model, RehMBEL1391, is composed of 1391 reactions including 229 transport reactions and 1171 metabolites. Constraints-based flux analyses were performed to refine and validate the genome-scale metabolic model under environmental and genetic perturbations. First, the lithoautotrophic growth characteristics of R. eutropha H16 were investigated under varying feeding ratios of gas mixture. Second, the genome-scale metabolic model was used to design the strategies for the production of poly[R-(-)-3hydroxybutyrate] (PHB) under different pH values and carbon/nitrogen source uptake ratios. It was also used to analyze the metabolic characteristics of R. eutropha when the phosphofructokinase gene was expressed. Finally, in silico gene knockout simulations were performed to identify targets for metabolic engineering essential for the production of 2-methylcitric acid in R. eutropha H16.ConclusionThe genome-scale metabolic model, RehMBEL1391, successfully represented metabolic characteristics of R. eutropha H16 at systems level. The reconstructed genome-scale metabolic model can be employed as an useful tool for understanding its metabolic capabilities, predicting its physiological consequences in response to various environmental and genetic changes, and developing strategies for systems metabolic engineering to improve its metabolic performance.

Design of homo-organic acid producing strains using multi-objective optimization.

Production of homo-organic acids without byproducts is an important challenge in bioprocess engineering to minimize operation cost for separation processes. In this study, we used multi-objective optimization to design Escherichia coli strains with the goals of maximally producing target organic acids, while maintaining sufficiently high growth rate and minimizing the secretion of undesired byproducts. Homo-productions of acetic, lactic and succinic acids were targeted as examples. Engineered E. coli strains capable of producing homo-acetic and homo-lactic acids could be developed by taking this systems approach for the minimal identification of gene knockout targets. Also, failure to predict effective gene knockout targets for the homo-succinic acid production suggests that the multi-objective optimization is useful in assessing the suitability of a microorganism as a host strain for the production of a homo-organic acid. The systems metabolic engineering-based approach reported here should be applicable to the production of other industrially important organic acids.

Development of a gene knockout system for Ralstonia eutropha H16 based on the broad-host-range vector expressing a mobile group II intron

Ralstonia eutropha H16 is a Gram-negative lithoautotrophic bacterium and is one of the best biopolymer-producing bacteria. It can grow to high cell densities either under lithoautotrophic or under heterotrophic conditions, which makes it suitable for a number of biotechnological applications. Also, R. eutropha H16 can degrade various aromatic compounds for environmental applications. The mobile group II intron can be used for the rapid and specific disruption of various bacterial genes by insertion into any desired target genes. Here, we applied the mobile group II intron to R. eutropha H16 and developed a markerless gene knockout system for R. eutropha: RalsTron. As a demonstration of the system, the phaC1 gene encoding polyhydroxyalkanoate synthase was successfully knocked out in R. eutropha H16. Furthermore, this knockout system would be useful for knocking out genes in other bacteria as well because it is based on a broad-host-range vector and the mobile group II intron that minimally depends on the bacterial hosts.