- Open Access
Designing intracellular metabolism for production of target compounds by introducing a heterologous metabolic reaction based on a Synechosystis sp. 6803 genome-scale model
Microbial Cell Factoriesvolume 15, Article number: 13 (2016)
Designing optimal intracellular metabolism is essential for using microorganisms to produce useful compounds. Computerized calculations for flux balance analysis utilizing a genome-scale model have been performed for such designs. Many genome-scale models have been developed for different microorganisms. However, optimal designs of intracellular metabolism aimed at producing a useful compound often utilize metabolic reactions of only the host microbial cells. In the present study, we added reactions other than the metabolic reactions with Synechosystis sp. 6803 as a host to its genome-scale model, and constructed a metabolic model of hybrid cells (SyHyMeP) using computerized analysis. Using this model provided a metabolic design that improves the theoretical yield of succinic acid, which is a useful compound.
Constructing the SyHyMeP model enabled new metabolic designs for producing useful compounds. In the present study, we developed a metabolic design that allowed for improved theoretical yield in the production of succinic acid during glycogen metabolism by Synechosystis sp. 6803. The theoretical yield of succinic acid production using a genome-scale model of these cells was 1.00 mol/mol-glucose, but use of the SyHyMeP model enabled a metabolic design with which a 33 % increase in theoretical yield is expected due to the introduction of isocitrate lyase, adding activations of endogenous tree reactions via D-glycerate in Synechosystis sp. 6803.
The SyHyMeP model developed in this study has provided a new metabolic design that is not restricted only to the metabolic reactions of individual microbial cells. The concept of construction of this model requires only replacement of the genome-scale model of the host microbial cells and can thus be applied to various useful microorganisms for metabolic design to produce compounds.
Designing optimal intracellular metabolism, as typified by metabolic engineering or synthetic biology, is essential when the aim is mass production of useful compounds by using microbial cells. One of the powerful tools for this design method is the use of genome-scale models (GSMs). In intracellular metabolic reactions, kinetics, of substrates and enzymes are involved, and various mathematical models that take these into consideration have been developed. However, the intracellular metabolism can be assumed to be in a steady state (pseudo-steady state) as long as the cells continue to be in the same environment. In other words, no change in the amount of intermediate metabolites occurs in the cell, and the metabolic flux in the cell can be predicted by mathematical calculation in the form of flux balance analysis (FBA). These computer-calculated prediction results successfully reflect the phenotype of the microbial cells in the actual experiment and have been reported to be highly accurate [1–3]. Many GSMs for different microbial cells have already been developed, and the number is now in excess of 100 (http://gcrg.ucsd.edu/InSilicoOrganisms/OtherOrganisms). Advances in the tools used to organize databases automatically and describe metabolic networks  have also led to a rapid increase in the number of new GSMs being developed. Metabolic design tools for various cells have been developed on the basis of these GSMs [5–8]. Using these tools has made it possible to identify the reaction pathways that need to be engineered and achieve high productivity of the target compound with high throughput. Indeed, metabolic prediction tools using GSMs from model microorganisms such as Escherichia coli, Saccharomyces cerevisiae, and Corynebacterium glutamicum have enabled high-efficiency production of compounds such as succinic acid , lactic acid , lycopene , valine , vanillin , and 1,4-butanediol . In addition, recent years have seen the introduction of omics data such as transcriptomics and metabolomics as parameters for the development of more accurate simulation tools [15–17]. However, metabolic reactions that can be used when performing calculations for the prediction of metabolic flux using GSMs are often confined mainly to reactions that occur in the host cells. For example, when GSMs for C. glutamicum or S. cerevisiae are used to simulate metabolic design for cells capable of high productivity of useful compounds, it is impossible to automatically use metabolic reactions of the Entner-Doudoroff pathway of E. coli or other cells . Moreover, when a simulation is run with E. coli as the host, the reaction pyruvic acid → oxaloacetic acid, which is catalyzed by pyruvate carboxylase in C. glutamicum, cannot be incorporated into the calculation . Thus, if GSMs only for host microbial cells in actual experiments are considered, it is possible that they end up being limitations to the repertoire of metabolic designs for high productivity of target compounds.
In the present study, we developed a tool with which metabolic reactions that do not belong to a host microorganism can be automatically added to a GSM of that host and simulations can be run. With this tool, individual candidate reactions can be selected from the metabolic reactions in the Kyoto Encyclopedia of Genes and Genomes (KEGG) to construct a computerized metabolic model of a hybrid cell to prepare a new metabolic design with FBA. This tool has been named the Hybrid Metabolic Pathway design tool (HyMeP).
Here, we propose an example of a metabolic design for constructing a HyMeP model for Synechosystis sp. 6803 (SyHyMeP) and maximizing the theoretical yield of succinic acid production during glycogen metabolism. Synechosystis sp. 6803, which is a model cyanobacterial microorganism, is able to use carbon dioxide as a carbon source during photosynthesis and store glycogen. The complete genome was read in 1996 , and many GSMs for it have been developed previously [21–27]. Research on the intracellular metabolism of cyanobacteria and on the production of useful compounds has also advanced in recent years [28–32], and studies have also been conducted on controlling their transcription factors to produce succinic acid during glycogen metabolism . Major advances are expected in the future in research on the production of useful compounds by cyanobacteria. Succinic acid is one of the bio-based 12 building blocks for useful chemical compounds selected by the U.S. Department of Energy (DOE) (http://www.energy.gov/), and, more importantly, it is a compound that is produced from non-fossil raw materials.
Results and discussion
Design of metabolic pathway for succinic acid production by extended GSM for Synechosystis sp. 6803 (SyHyMeP).
In the present study, we used the Synechosystis sp. 6803 GSM (SyGSM) developed by Nogales et al. . Metabolic reactions were randomly selected out of all the metabolic reactions obtained from KEGG and added to the SyGSM, and maximum production flux of succinic acid was calculated by FBA. Table 1 shows the numbers of exogenous reactions showing binding to the metabolites in the SyGSM at a particular time. More detailed information is provided in Additional file 1: Table S1. When these were added to the SyGSM and calculations were run to maximize the succinic acid production flux by FBA, we obtained a metabolic design with an improved theoretical yield compared to that with the design when calculations were done only with the SyGSM (Table 2). The focus was on adding reactions where there was a possibility of increase in succinic acid yield. Introduction of isocitrate lyase (KEGG reaction number R00479) was found to be desirable for enabling efficient succinic acid production in Synechosystis sp. 6803.
R00479: Isocitrate → Succinate + Glyoxylate
This prediction result is similar to that for metabolic design during succinic acid production in E. coli [33–36] or C. glutamicum [37, 38]. The existence of fumarate reductases (frdA, B, C, D) as found in E. coli and similar organisms has not been confirmed in Synechosystis sp. 6803. In fact, the reaction fumarate → succinate has not even been introduced in Synechosystis sp. 6803 [21–27]. To that end, succinic acid needs to be produced from citrate in the TCA cycle, and in the original SyGSM, CO2 necessarily needs to be discharged. Introducing isocitrate lyase allows carbon to be consumed more sparingly, leading to succinic acid production, and is regarded as a preferred choice in the SyHyMeP as well. Synechosystis sp. 6803 does not have a complete TCA cycle, and according to the SyGSM, succinic acid needs to be generated from isocitrate by a GABA-mediated reaction. Introducing isocitrate lyase is therefore regarded as effective because succinic acid can be produced by skipping this relatively long pathway. Another compound that is produced along with isocitrate lyase, that is, glyoxylate, is preferably returned to the glycolytic system by three reactions carried out by Synechosystis sp. 6803 (Fig. 1). The succinic acid production yield is potentially increased by up to 33 % due to the introduction of exogenous isocitrate lyase and by the enhanced activity of three reactions that return glyoxylate to the glycolytic system.
R00751 is listed as a candidate reaction for increasing succinic acid production yield.
R00751: l-Threonine → Glycine + Acetaldehyde
The C2 glycine that is generated is returned to the glycolytic pathway via C3 serine by tetrahydrofolate-mediated C1 metabolism, thus allowing for recovery of the carbon. The acetaldehyde that is generated can be returned to the central metabolic pathway, via conversion into acetyl-CoA from acetic acid. Succinic acid yield can also be increased by the reaction R01867.
R01867: (S)-Dihydroorotate + Fumarate → Orotate + Succinate
The SyGSM includes a reaction where orotate and hydrogen peroxide are combined as (S)-dihydroorotate; therefore, functionalizing both reactions would yield a reaction similar to that catalyzed by fumarate reductases of E. coli and other organisms. However, this may actually not be practical, because it requires that hydrogen peroxide, which is toxic to the cell, be generated within the cell.
In some GSMs, adding two or more reactions increased succinic acid yield, compared to that in GSMs where only one reaction was added, but these would encompass any of the three reactions mentioned above (Table 2). According to the SyHyMeP, if the theoretical yield of succinic acid production was maximum when R00751 and R00479 are combined, then, it was 155 if a third reaction R10179 was also included. However, it was necessary to activate 19 reactions in Synechosystis sp. 6803 besides the three reactions mentioned above, and thus, this design for succinate production is not thought to be realistic (Additional file 2: Table S2). In SyHyMeP, no combination that provided improved yield of succinic acid production was found even when four or five reactions were added.
There have previously been limitations for succinic acid production yield with SyGSM alone. However, introducing the SyHyMeP model has made it possible to automatically add exogenous reactions from other species and to obtain new metabolic designs for improving the production yield of target compounds. It is difficult to produce succinic acid metabolism as designed as shown in Fig. 1, and thus, after introducing the SyHyMeP model, use of algorithms, that take into account both growth and the flux of the target compound by employing tools such as OptForce , MOMA , and OP-Synthetic , is necessary to accomplish a more efficient metabolic design. Moreover, in the actual production of a strain, artificial metabolic switching systems are considered necessary to induce changes in the metabolic flux for growth to produce metabolic flux geared towards the production of the target compound [41–43].
In the present study, we designed the SyHyMeP, which automatically adds other metabolic reactions from the KEGG to the metabolic reaction model of Synechosystis sp. 6803 (SyGSM) and enables new metabolic designs. Thus, new metabolic designs pertaining to succinic acid production during glycogen metabolism in Synechosystis sp. 6803 have become possible. The optimal combination of metabolic reactions that could increase the yield of succinic acid production was when isocitrate lyase, represented by KEGG reaction number R00479, was introduced; in this case, three reactions were highly activated via intracellular glyceric acid. The theoretical yield for succinic acid at this time was 133. The HyMeP model proposed here was developed for Synechosystis sp. 6803, but the concept of construction of this model requires only replacing the GSM of host microbial cells, and thus can be applied to various microorganisms for developing metabolic designs to produce useful compounds. In other words, it is now possible to enable new metabolic designs that could not be designed with GSMs previously constructed for the metabolic reactions of individual microbial cells.
The present study is based on the Synechosystis sp. 6803 GSM (SyGSM) developed by Nogales et al. . To calculate the succinic acid production flux during glycogen metabolism, the oxygen uptake, and the uptake of CO2 from photosynthesis into the cells were set as 0, respectively. SyHyMeP was constructed by using the following procedure: Reactions included in the SyGSM were extracted from the metabolic reaction list provided by the KEGG (http://www.genome.jp/kegg/), and then the other reactions were set as a exSy_list. Inorganic compounds such as H2O and ATP and metabolites such as co-factors were excluded (Additional file 3: Table S3, Sy-metabolites) from the metabolites in the SyGSM. These data were used to implement the SyHyMeP as follows: In each metabolic reaction formula in the exSy_list, we determined if metabolites in Sy_metabolites were present in either the Reactants or Products; if any were present, we searched if there were any identical metabolites in the remaining reaction formulae in the exSy_list, with respect to the metabolites on the other side of the reaction. We searched for identical metabolites in either the reactants or products, and if the metabolic reactions were linked, we focused on the metabolite on the other side of the reaction. This procedure was repeated to search and link a specified number of reactions. The SyHyMeP implementation flow diagram is shown in Fig. 2 where, for example, two reactions have been shown to be newly added. At a few instances, the total number of reactions was two; in some cases, single-reaction increments were independently linked to the SyGSM (Fig. 2a), and still in others, each of the reactions was linked (Fig. 2b). Additional file 4: Table S4 (nmlist) summarizes such combinations of up to five reactions. Next, we created a list of reactions that lead to metabolites in the SyGSM in exSy_list (Addditional file 1: Table 1), for each reaction number (maximum reaction number 5). If the total reaction number was set to two, then two reactions (each being n1 in the Additional file 4: Table S4) were combined and added to the SyGSM to prepare the SyHyMeP, and the production yield of succinic acid can be calculated as an FBA (Fig. 2, pattern 1). In addition, each reaction was added to the SyGSM from the list for a reaction number of two (n2 in the Additional file 4: Table S4) to make a separate SyHyMeP, which could be calculated with FBA (Fig. 2, pattern 2).
For metabolic simulation, cytosolic glycogen was used as the sole carbon source, and the uptake rate was set to 100. Other external metabolites such as NH3 and CO2 could be transported freely through the cell membrane. All simulations were performed using the Java language. SyHyMeP is a web-based application implemented in Java and runs on the Apache Tomcat web server (ver. 7.0.64., http://tomcat.apache.org/). FBA was performed with our original command line tool using OptFlux libraries, which use GLPK (GNU Linear Programming Kit) as a linear programming solver .
genome scale model
flux balance analysis
Hybrid Metabolic Pathway Design Tool
genome scale model of Synechosystis sp. 6803
Hybrid Metabolic Pathway design tool of Synechosystis sp. 6803
Varma A, Palsson BO. Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110. Appl Environ Microbiol. 1994;60:3724–31.
Edwards JS, Ibarra RU, Palsson BO. In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data. Nat Biotechnol. 2001;19:125–30.
Lee J, Yun H, Feist AM, Palsson BO, Lee SY. Genome-scale reconstruction and in silico analysis of the Clostridium acetobutylicum ATCC 824 metabolic network. Appl Microbiol Biotechnol. 2008;80:849–62.
Thiele I, Palsson BO. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010;5:93–121.
Segre D, Vitkup D, Church GM. Analysis of optimality in natural and perturbed metabolic networks. Proc Natl Acad Sci USA. 2002;99:15112–7.
Burgard AP, Pharkya P, Maranas CD. Optknock: a bilevel programming framework for identifying gene knockout strategies for microbial strain optimization. Biotechnol Bioeng. 2003;84:647–57.
Pharkya P, Burgard AP, Maranas CD. OptStrain: a computational framework for redesign of microbial production systems. Genome Res. 2004;14:2367–76.
Rocha I, Maia P, Evangelista P, Vilaca P, Soares S, Pinto JP, et al. OptFlux: an open-source software platform for in silico metabolic engineering. BMC Syst Biol. 2010;4:45.
Kim TY, Kim HU, Park JM, Song H, Kim JS, Lee SY. Genome-scale analysis of Mannheimia succiniciproducens metabolism. Biotechnol Bioeng. 2007;97:657–71.
Fong SS, Burgard AP, Herring CD, Knight EM, Blattner FR, Maranas CD, et al. In silico design and adaptive evolution of Escherichia coli for production of lactic acid. Biotechnol Bioeng. 2005;91:643–8.
Choi HS, Lee SY, Kim TY, Woo HM. In silico identification of gene amplification targets for improvement of lycopene production. Appl Environ Microbiol. 2010;76:3097–105.
Park JH, Lee KH, Kim TY, Lee SY. Metabolic engineering of Escherichia coli for the production of l-valine based on transcriptome analysis and in silico gene knockout simulation. Proc Natl Acad Sci USA. 2007;104:7797–802.
Brochado AR, Matos C, Moller BL, Hansen J, Mortensen UH, Patil KR. Improved vanillin production in baker’s yeast through in silico design. Microb Cell Fact. 2010;9:84.
Yim H, Haselbeck R, Niu W, Pujol-Baxley C, Burgard A, Boldt J, et al. Metabolic engineering of Escherichia coli for direct production of 1,4-butanediol. Nat Chem Biol. 2011;7:445–52.
Becker SA, Palsson BO. Context-specific metabolic networks are consistent with experiments. PLoS Comput Biol. 2008;4:e1000082.
Bordbar A, Mo ML, Nakayasu ES, Schrimpe-Rutledge AC, Kim YM, Metz TO, et al. Model-driven multi-omic data analysis elucidates metabolic immunomodulators of macrophage activation. Mol Syst Biol. 2012;8:558.
Schmidt BJ, Ebrahim A, Metz TO, Adkins JN, Palsson BO, Hyduke DR. GIM3E: condition-specific models of cellular metabolism developed from metabolomics and expression data. Bioinformatics. 2013;29:2900–8.
Eisenberg RC, Dobrogosz WJ. Gluconate metabolism in Escherichia coli. J Bacteriol. 1967;93:941–9.
Peters-Wendisch PG, Kreutzer C, Kalinowski J, Patek M, Sahm H, Eikmanns BJ. Pyruvate carboxylase from Corynebacterium glutamicum: characterization, expression and inactivation of the pyc gene. Microbiology. 1998;144:915–27.
Kaneko T, Sato S, Kotani H, Tanaka A, Asamizu E, Nakamura Y, et al. Sequence analysis of the genome of the unicellular cyanobacterium Synechocystis sp. strain PCC6803. II. Sequence determination of the entire genome and assignment of potential protein-coding regions. DNA Res. 1996;3:109–36.
Shastri AA, Morgan JA. Flux balance analysis of photoautotrophic metabolism. Biotechnol Prog. 2005;21:1617–26.
Fu PC. Genome-scale modeling of Synechocystis sp. PCC 6803 and prediction of pathway insertion. J Chem Technol Biotechnol. 2009;84:473–83.
Knoop H, Zilliges Y, Lockau W, Steuer R. The metabolic network of Synechocystis sp. PCC 6803: systemic properties of autotrophic growth. Plant Physiol. 2010;154:410–22.
Montagud A, Navarro E, de Cordoba FP, Urchueguia JF, Patil KR. Reconstruction and analysis of genome-scale metabolic model of a photosynthetic bacterium. BMC Syst Biol. 2010;4:156.
Montagud A, Zelezniak A, Navarro E, de Cordoba PF, Urchueguia JF, Patil KR. Flux coupling and transcriptional regulation within the metabolic network of the photosynthetic bacterium Synechocystis sp. PCC6803. Biotechnol J. 2011;6:330–42.
Yoshikawa K, Kojima Y, Nakajima T, Furusawa C, Hirasawa T, Shimizu H. Reconstruction and verification of a genome-scale metabolic model for Synechocystis sp. PCC6803. Appl Microbiol Biotechnol. 2011;92:347–58.
Nogales J, Gudmundsson S, Knight EM, Palsson BO, Thiele I. Detailing the optimality of photosynthesis in cyanobacteria through systems biology analysis. Proc Natl Acad Sci USA. 2012;109:2678–83.
McNeely K, Xu Y, Bennette N, Bryant DA, Dismukes GC. Redirecting reductant flux into hydrogen production via metabolic engineering of fermentative carbon metabolism in a cyanobacterium. Appl Environ Microbiol. 2010;76:5032–8.
McNeely K, Kumaraswamy GK, Guerra T, Bennette N, Ananyev G, Dismukes GC. Metabolic switching of central carbon metabolism in response to nitrate: application to autofermentative hydrogen production in cyanobacteria. J Biotechnol. 2014;182–183:83–91.
Baebprasert W, Jantaro S, Khetkorn W, Lindblad P, Incharoensakdi A. Increased H2 production in the cyanobacterium Synechocystis sp. strain PCC 6803 by redirecting the electron supply via genetic engineering of the nitrate assimilation pathway. Metab Eng. 2011;13:610–6.
Osanai T, Shirai T, Iijima H, Kuwahara A, Suzuki I, Kondo A, et al. Alteration of cyanobacterial sugar and amino acid metabolism by overexpression hik8, encoding a KaiC-associated histidine kinase. Environ Microbiol. 2015;17:2430–40.
Osanai T, Oikawa A, Numata K, Kuwahara A, Doi Y, et al. Pathway-level acceleration of glycogen catabolism by a response regulator in the cyanobacterium Synechocystis sp. PCC 6803. Plant Physiol. 1831;164:1841.
Yang J, Wang Z, Zhu N, Wang B, Chen T, Zhao X. Metabolic engineering of Escherichia coli and in silico comparing of carboxylation pathways for high succinate productivity under aerobic conditions. Microbiol Res. 2014;169:432–40.
Sanchez AM, Bennett GN, San KY. Novel pathway engineering design of the anaerobic central metabolic pathway in Escherichia coli to increase succinate yield and productivity. Metab Eng. 2005;7:229–39.
Lin H, Bennett GN, San KY. Metabolic engineering of aerobic succinate production systems in Escherichia coli to improve process productivity and achieve the maximum theoretical succinate yield. Metab Eng. 2005;7:116–27.
Zhu LW, Li XH, Zhang L, Li HM, Liu JH, Yuan ZP, et al. Activation of glyoxylate pathway without the activation of its related gene in succinate-producing engineered Escherichia coli. Metab Eng. 2013;20:9–19.
Zhu N, Xia H, Wang Z, Zhao X, Chen T. Engineering of acetate recycling and citrate synthase to improve aerobic succinate production in Corynebacterium glutamicum. PLoS One. 2013;8:e60659.
Litsanov B, Kabus A, Brocker M, Bott M. Efficient aerobic succinate production from glucose in minimal medium with Corynebacterium glutamicum. Microb Biotechnol. 2012;5:116–28.
Ranganathan S, Suthers PF, Maranas CD. OptForce: an optimization procedure for identifying all genetic manipulations leading to targeted overproductions. PLoS Comput Biol. 2010;6:e1000744.
Liu H, Li Y, Wang X. OP-synthetic: identification of optimal genetic manipulations for the overproduction of native and non-native metabolites. Quant Biol. 2014;2:100–9.
Farmer WR, Liao JC. Improving lycopene production in Escherichia coli by engineering metabolic control. Nat Biotechnol. 2000;18:533–7.
Zhang F, Carothers JM, Keasling JD. Design of a dynamic sensor-regulator system for production of chemicals and fuels derived from fatty acids. Nat Biotechnol. 2012;30:354–9.
Soma Y, Tsuruno K, Wada M, Yokota A, Hanai T. Metabolic flux redirection from a central metabolic pathway toward a synthetic pathway using a metabolic toggle switch. Metab Eng. 2014;23:175–84.
Orth JD, Thiele I, Palsson BO. What is flux balance analysis? Nat Biotechnol. 2010;28:245–8.
Oberhardt MA, Palsson BO, Papin JA. Applications of genome-scale metabolic reconstructions. Mol Syst Biol. 2009;5:320.
TS performed the experiment; TS, TO, and AK designed the research; TS wrote the paper. All authors read and approved the final manuscript.
This work was supported by the Ministry of Education, Culture, Sports, Science, and Technology, Japan, by a grant to T.O. from ALCA (Project name “Production of cyanobacterial succinate by the genetic engineering of transcriptional regulators and circadian clocks”) from the Japan Science and Technology Agency.
The authors declare that they have no competing interests.