Genome-scale reconstruction and in silico analysis of Klebsiella oxytoca for 2,3-butanediol production
© Park et al.; licensee BioMed Central Ltd. 2013
Received: 4 January 2013
Accepted: 17 February 2013
Published: 23 February 2013
Klebsiella oxytoca, a Gram-negative, rod-shaped, and facultative anaerobic bacterium, is one of the most promising 2,3-butanediol (2,3-BD) producers. In order to improve the metabolic performance of K. oxytoca as an efficient biofactory, it is necessary to assess its metabolic characteristics with a system-wide scope, and to optimize the metabolic pathways at a systems level. Provision of the complete genome sequence of K. oxytoca enabled the construction of genome-scale metabolic model of K. oxytoca and its in silico analyses.
The genome-scale metabolic model of K. oxytoca was constructed using the annotated genome with biochemical and physiological information. The stoichiometric model, KoxGSC1457, is composed of 1,457 reactions and 1,099 metabolites. The model was further refined by applying biomass composition equations and comparing in silico results with experimental data based on constraints-based flux analyses. Then, the model was applied to in silico analyses to understand the properties of K. oxytoca and also to improve its capabilities for 2,3-BD production according to genetic and environmental perturbations. Firstly, in silico analysis, which tested the effect of augmenting the metabolic flux pool of 2,3-BD precursors, elucidated that increasing the pyruvate pool is primarily important for 2,3-BD synthesis. Secondly, we performed in silico single gene knockout simulation for 2,3-BD overproduction, and investigated the changes of the in silico flux solution space of a ldhA gene knockout mutant in comparison with that of the wild-type strain. Finally, the KoxGSC1457 model was used to optimize the oxygen levels during fermentation for 2,3-BD production.
The genome-scale metabolic model, KoxGSC1457, constructed in this study successfully investigated metabolic characteristics of K. oxytoca at systems level. The KoxGSC1457 model could be employed as an useful tool to analyze its metabolic capabilities, to predict its physiological responses according to environmental and genetic perturbations, and to design metabolic engineering strategies to improve its metabolic performance.
There is growing interest in the production of 2,3-butanediol (2,3-BD) by microbial fermentation, as it can be easily converted to methyl ethyl ketone and tetramethyl ether, blending agents for gasoline, and 1,3-butadiene, an intermediate in synthetic rubber manufacture [1–4]. Klebsiella oxytoca is known as one of the most promising 2,3-BD producers [1–10], and its whole genome sequences have been reported recently. The genome of K. oxytoca KCTC1686 consists of a chromosome of 5,974,109 bp with a 56.05% GC content, including 5,488 coding genes . The genome of K. oxytoca E718 is composed of a chromosome of 6,097,032 bp and two plasmids of 324,906 bp and 110,781 bp with a 55.5% GC content, including 5,909 coding genes . More recently, we isolated K. oxytoca KCTC12133BP from a cattle farm and sequenced its whole genome, but the sequence information, which has not yet been published, consists of a chromosome of 5,903,932 bp and a plasmid of 109,773 bp with a 55.4% GC content, including 5,793 coding genes.
The most important characteristics of K. oxytoca is to produce large amounts of C3/C4 diols, 1,3-propanediol (1,3-PD) and 2,3-BD, using various carbon sources [10, 13–20]. For example, a lactate deficient mutant of K. oxytoca coproduced 83.6 and 60.1 g/L of 1,3-PD and 2,3-BD, respectively, in a fed-batch fermentation utilizing mixed substrates of glycerol and sucrose . K. oxytoca also could produce more than 95 g/L of 2,3-BD from glucose with a yield of 0.478 g/g (95.6% of the theoretical maximum yield) and a productivity of 1.71 g/L/h in a batch fermentation, in which the agitation speed was switched from 300 to 200 rpm during the fermentation . The availability of oxygen to K. oxytoca significantly affects its physiology and 2,3-BD production [21–23]. Another striking aspect of K. oxytoca is to readily metabolize glycerol, an inevitable by-product of biodiesel production, into biomass and products of value [13, 24, 25].
The aforementioned advantages make K. oxytoca an attractive host for industrial applications. Furthermore, the U.S. National Institute of Health (NIH, Guidelines for Research Involving Recombinant DNA Molecules, 2002) has reported that K. oxytoca belongs to risk group 1 (RG 1), recognizing it as a GRAS (Generally Regarded As Safe) organism. However, in order to use this organism on an industrial scale, the strain should be further developed. Systems metabolic engineering allows the rational design of metabolic networks for the overproduction of target compounds and the creation of industrially useful microorganisms [26–31]. Here, and in silico genome-scale metabolic model of K. oxytoca, KoxGSC1457, was constructed based on genome information, databases, and experimental data. The KoxGSC1457 model is composed of 1,457 reactions and 1,099 metabolites (Additional file 1 and Additional file 2). The model was carefully examined by in silico analyses for genetic and environmental perturbations. The in silico analysis using the model predicted that the pyruvate pool is mostly important for 2,3-BD synthesis, and this was verified by fermentation of the ldhA gene knockout mutant. Also, the model showed that the availability of oxygen strongly affected the production of 2,3-BD by K. oxytoca.
Results and discussion
Genome-scale reconstruction and general features of K. oxytoca metabolic network
Features of the in silico genome-scale metabolic model of K. oxytoca
Genome size (base pairs, bp)
No. of open reading frames (ORFs)
In silico metabolic model
No. of reactions (redundant) included in the model
No. of biochemical reactions
No. of transport reactions
No. of reactions (unique) included in the model
No. of metabolites
No. of ORFs assigned in metabolic network
ORF coverageb (%)
The metabolic reactions in KoxGSC1457 were classified into 8 different subsystems, comprising 97 sub-metabolisms: amino acid metabolism, carbohydrate metabolism, metabolism of cofactors and vitamins, nucleotide metabolism, lipid metabolism, energy metabolism, transporters, and others (Figure 2 and Additional file 1). Figure 2A details the numbers and percentages of reactions associated with genes for each subsystems. Amino acid metabolism ranks as the largest subsystem in KoxGSC1457, followed by carbohydrate metabolism. The sum of the three largest subsystems, amino acid metabolism, carbohydrate metabolism, and metabolism of cofactors and vitamins, account for more than a half of the total number of reactions. The percentage of reactions assigned to ORFs in KoxGSC1457 was 79.7%. The remaining 20.3% of reactions included in the model contained lumped multi-step reactions, spontaneous reactions, reactions added to fill missing links, or several transport reactions (Table 1). More than 90% of the reactions in 6 subsystems, except for transporters and energy metabolism are associated with genes (Figure 2A). Figure 2B and C describe the classification of genes and metabolites in KoxGSC1457. The smaller percentages of unique metabolites for each category compared with those of unique genes imply that several metabolites participate with reactions in different subsystems, while most genes function to a specific subsystem (Figure 2B).
The KoxGSC1457 model contains the reactions for 2,3-BD biosynthesis catalyzed by acetolactate synthase (E.C. 22.214.171.124), acetolactate decarboxylase (E.C. 126.96.36.199), and acetoin reductase (E.C. 188.8.131.52). The model also includes a spontaneous reaction that converts α-acetolactate into diacetyl (C4H6O2) and CO2 in the presence of oxygen. Diacetyl can then be catalyzed into acetoin (C4H8O2), which is the precursor of 2,3-BD, by diacetyl reductase (E.C. 184.108.40.2063). K. oxytoca has its specialty for glycerol utilization by oxidative and reductive routes. Thus, the model encompasses the reactions for glycerol utilization, which are catalyzed by glycerol kinase (E.C. 220.127.116.11) and glycerol-3-phosphate dehydrogenase (E.C. 18.104.22.168) of the oxidative route, and NAD-dependent glycerol dehydrogenase (E.C. 22.214.171.124) and dihydroxyacetone kinase (E.C. 126.96.36.199) of the reductive route. Finally, the refined KoxGSC1457 model was used to understand the metabolic characteristics of K oxytoca and was applied to analyze the strategies for 2,3-BD production.
Effects of increasing the metabolic flux pools of 2,3-BD precursors in K. oxytoca
Likewise, the KoxGSC1457 model was used to examine the importance of 2,3-BD precursors, and it gave better insights into augmenting the metabolic flux pools of 2,3-BD precursors for 2,3-BD overproduction (Figure 3B and C). In order to increase the metabolic flux pools of 2,3-BD precursors as well as to investigate its effects on 2,3-BD overproduction, all reactions participating in the consumption of target precursors were eliminated except for reactions involved in 2,3-BD biosynthetic pathways and essential reactions required for cell growth on the basis of the wild-type strain of KoxGSC1467, as described in Figure 3B. Accordingly, the flux solution spaces were scanned by increasing the flux pool of each precursor (Figure 3C). As a result, increasing the pools of most precursors had no effects on the production of 2,3-BD compared to that of the wild-type strain, as determined through comparison of the flux solution spaces. The model predicted that the surplus fluxes from artificially increased flux pools of most 2,3-BD precursors were not redirected to 2,3-BD synthesis but went toward other pathways for survival related with biomass synthesis, such as amino acids, nucleotides, lipids, and cofactors, prior to arriving at 2,3-BD biosynthetic reactions. In some cases, artificially increasing the flux pool of a 2,3-BD precursor (e.g. phosphoenolpyruvate) had negative influences on both cell growth and 2,3-BD synthesis by unbalancing the precursor requirements for biomass synthesis. However, if the pyruvate pool was augmented sufficiently, the flux solution space was altered so that KoxGSC1457 was enabled to overproduce 2,3-BD by lowering byproducts formation rates at an optimal cell growth rate (Figure 3C). The surplus of the pyruvate pool was sufficiently redirected to 2,3-BD synthesis, resulting in an approximate doubling of the 2,3-BD production rate. Correspondingly, byproducts formation rates decreased dramatically compared with those of the wild-type strain, and consequently, augmentation of the pyruvate pool is one of most important keys for 2,3-BD overproduction.
In silico single gene knockout simulation and changes of flux solution space by ldhA gene knockout in K. oxytoca
During fermentation by the K. oxytoca wild-type strain, several byproducts, including lactic acid, ethanol, and formic acid were produced (Figure 4). For the overproduction of 2,3-BD in K. oxytoca by reducing the formation of byproducts and redirecting the remaining metabolic fluxes towards the synthesis of 2,3-BD, in silico single gene knockout simulation was performed using flux balance analysis (FBA). As a result, the knockout of ldhA gene encoding lactate dehydrogenase, which converts pyruvate into lactic acid with NADH oxidation, was targeted with top priority for a single gene knockout strategy based on the criteria of the maximization of 2,3-BD production rate and minimization of byproducts formation rates. Then, the KoxGSC1457 model investigated the metabolic characteristics of the ldhA gene knockout mutant by examining the changes of flux solution space compared with those of the wild-type strain (Figure 4A and B). In the wild-type strain, the model predicted that K. oxytoca produces large amounts of byproducts along with 2,3-BD. In the ldhA gene knockout mutant, the optimal point shifted to an increase of 2,3-BD production but to a decreased formation of byproducts, including lactic acid. ldhA gene knockout significantly increased the pyruvate pool that redirected the metabolic flux to 2,3-BD synthesis. This prediction was validated by batch fermentations of K. oxytoca wild-type and ldhA gene knockout strains (Figure 4C and D). K. oxytoca wild-type strain produced 32 g/L of lactic acid and 17 g/L of 2,3-BD by consuming 90 g/L of glucose. However, the ldhA gene knockout mutant dramatically decreased lactic acid production to 1.9 g/L, which is about 6% of that produced by the wild-type strain (Figure 4D). Accordingly, 2,3-BD production in the mutant was increased to about 30 g/L, which is about 176% of that produced by the wild-type strain. However, cell growth and glucose uptake rates of the ldhA gene knockout mutant were maintained in comparison with those of the wild-type strain because the redox balance of NAD+/NADH was not disrupted (Figure 4C and D). The deletion of lactate dehydrogenase forming lactic acid with the oxidation of NADH seemed to be compensated with 2,3-BD dehydrogenase, producing 2,3-BD with the oxidation of NADH as well, which made 2,3-BD production increase without any retardation of growth. When we compared the in silico flux solution space of the ldhA gene knockout mutant in Figure 4B with that of the increasing highly pyruvate pool in Figure 3C, more improvements for 2,3-BD overproduction could be made by further genetic manipulation. In particular, the fermentation results of the ldhA gene knockout mutant in Figure 4D suggests the construction of a mutant preventing the formation of ethanol, formic acid, and succinic acid. For this, the KoxGSC1457 model will be applicable to identify suitable candidates for further genetic manipulations.
Effects of varying oxygen uptake rate for 2,3-BD production in K. oxytoca
This simulation was validated by batch fermentation experiments under different agitation speeds (rpm) and DO levels (Figure 5B and C). The agitation speeds were a150, 250, 350, 450, 550, 650, and 750 rpm. DO levels were maintained at 5% and 10% by automatically adjusting the agitation speed between 50 and 1,000 rpm. The aeration rate was set at 1 vvm (air volume∙ working volume-1∙ minute-1) during the whole period of fermentation. The fermentations were performed by using the ldhA gene knockout mutant of K. oxytoca, as the most productive 2,3-BD producer up to this point. As a result, the final 2,3-BD concentration (g/L) and yield (g/g glucose) increased through increasing the agitation speed to 350 rpm, but further increasing agitation speed caused a decrease of the final 2,3-BD concentration and yield. Meanwhile, the maximal value of 2,3-BD productivity of 2.1 g/L/h was observed at 450 rpm. The highest glucose consumption was observed at 550 rpm, however, most consumption of carbon source was used to synthesize not 2,3-BD but biomass at 550 rpm and agitation speeds greater than 550 rpm. In addition, fermentation at 10% DO level showed less 2,3-BD concentration, yield, and productivity than at 5% DO. These results are a strong indication that excessive oxygen exposure during fermentation impairs 2,3-BD overproduction. Since the optimal oxygen level exists beyond the scope of DO as represented in Figure 5, it is recommended to adjust the agitation speed depending on cellular status for 2,3-BD overproduction. These observations were also supported by a previous study, in which the optimal agitation speed existed on 2,3-BD production .
K. oxytoca is a promising microorganism which has great potential in the industrial production of 2,3-BD. To make the better use of this organism and to improve industrial applicability, we constructed a genome-scale metabolic model of K. oxytoca, KoxGSC1457, consisting of 1,457 reactions and 1,099 metabolites based on information obtained from genome annotation data, database, literature and validation experiments. Using the KoxGSC1457 model, in silico metabolic fluxes were analyzed to investigate the effects of environmental and genetic variation. First, the effects of increasing the metabolic flux pools of 2,3-BD precursors were examined to understand the most important precursor to be augmented for 2,3-BD overproduction. In silico single gene knockout simulation was then performed to identify gene knockout candidates for the reduction of byproduct formation, and the change of flux solution space was then analyzed by elimination of ldhA gene resulting in the K. oxytoca ldhA gene knockout mutant, which is one of organisms that shows best performance for 2,3-BD production so far. Finally, to design strategies for better 2,3-BD overproduction, in silico metabolic flux analyses of KoxGSC1457 were executed to examine the effects of varying oxygen uptake rate during 2,3-BD production, and to determine the optimal agitation speed. Through the successful use of the KoxGSC1457 genome-scale metabolic model like those, it gave us the confidence that the model was thoroughly validated in comparison with various experimental data. There are still many difficulties associated with the overall economics of producing 2,3-BD through K. oxytoca on an industrial scale. Hence, it is expected that the model constructed in this study will support us to understand cellular physiology on systems level, broaden our insight on this organism, and conquer problems more systematically.
Reconstruction of genome-scale metabolic model
Physiological data and biochemical reactions of K. oxytoca were collected from various sources, such as experiments, literatures, and public databases, in order to reconstruct the genome-scale metabolic model of K. oxytoca. First, information of genome annotation was gained from genome data in NCBI (http://www.ncbi.nlm.nih.gov/). Then, the biochemical reactions of K. oxytoca were assembled from various literatures, experiments, and public databases, including the Kyoto Encyclopedia of Genes and Genomes (KEGG, ), Biocyc , Metacyc , TCDB , and TransportDB . Second, the draft version of the metabolic model was manually curated and refined based on experimentation, literature, databases, gene-protein-reaction correlations, comparative genome analysis, error corrections, and gap filling. Third, the refined model was validated through comparison with experimental data, including fermentation data. The biomass composition data of K. oxytoca were then analyzed and applied to the model. Finally, the validated model was used to analyze the physiology of K. oxytoca and to apply to in silico simulations for strain improvement and fermentation optimization (Figure 1).
The composition of biomass components was experimentally measured to construct the biomass formation reactions of K. oxytoca more accurately. The cells of K. oxytoca were cultured in minimal culture medium containing no yeast extract, and the samples were collected in the mid-exponential growth phase. The carbohydrate compositions were determined for neutral and amino sugars (Korea Basic Science Institute, Daejeon, Korea). The neutral and amino sugars were extracted from cells by treating them with 2 M trifluoroacetic acid and 6 N HCl, respectively, at 100°C for 4 hours. The carbohydrates were analyzed by ICD-5000 (Dionex, Sunnyvale, CA, USA) equipped with a CarboPac PA10 column (4.5 × 250 mm, Dionex) and a CarboPac PA10 cartridge (4 × 50 mm). 16 mM NaOH was used as a mobile phase with a flow rate of 1.0 mL∙ min-1. Then, data were analyzed by Chromelon ver 6.8 software. Amino acid compositions were analyzed by a Hewlett Packard 1100 series HPLC systems equipped with Waters Nova-Pak C18 4 um column (3.9 × 300 mm) (Korea Basic Science Institute, Daejeon, Korea). The fatty acid composition of K. oxytoca was determined by Sherlock microbial identification system of Sherlock version 6.1 (Korea Research Institute of Bioscience and Biotechnology, Daejeon, Korea). DNA and RNA compositions were analyzed based on genome information of K. oxytoca. The DNA composition was calculated by using the GC content (53%) of K. oxytoca. The RNA composition was calculated based on nucleotide sequences for rRNA and tRNA derived from genome information. The compositions for other components were obtained from literature or assumed reasonably described in Additional file 3.
Constraints-based flux analysis
In order to perform in silico simulations, and to predict the metabolic characteristics of K. oxytoca, constraints-based flux analysis, including FBA, was carried out under the assumption of a pseudo-steady state [35, 40–43]. In order to simulate the in silico model more accurately, the limits of uptake and secretion rates for some metabolites, including amino acids and organic acids such as acetic acid, ethanol, formic acid, lactic acid, pyruvic acid, and succinic acid, were constrained by experimentally measured flux values. Likewise, the limits related with secretion for some metabolites, which were not produced during fermentation such amino acids, were constrained to zero.
The effect of 2,3-BD production rate in response to varying oxygen uptake rate was examined by flux response analysis [44–46]. The 2,3-BD production rate, as objective function, was maximized and minimized according to the changes of oxygen uptake rate from 0 to 30 mmol/gDCW/h with a 10 mmol/gDCW/h glucose uptake rate.
For seed preparation, suspended cells from single colonies on Luria-Bertani (LB) agar (Difco Laboratories, Detroit, MI) plates were precultured in 20 mL test tubes containing 5 mL culture medium at 37°C for 5 h. 1 mL aliquots of the preculture were then transferred to 500 mL Erlenmeyer flasks containing 300 mL culture medium, and the cells were cultivated to an optical density of 1.5 ~ 2.0 at 600 nm (OD600). The tube and flask cultivations were placed in a rotary shaker at 150 rpm and 37°C (JEIO Tech. Co. SI-900R). The culture medium contained per liter: yeast extract (Becton Dickinson, Le Pont de Claix, France), 5 g; FeSO4 · 7H2O, 0.05 g; ZnSO4 · 7H2O, 0.001 g, MnSO4 · H2O, 0.001 g; CaCl2∙ 2H2O, 0.001 g; MgSO4∙ 7H2O, 0.25 g; (NH4)2SO4, 6.6 g; K2HPO4, 8.7 g; KH2PO4, 6.8 g; trace metal solution, 10 mL. The trace metal solution contains per liter: FeSO4∙ 7H2O, 5 g; ZnSO4∙ 7H2O, 0.1 g; MnSO4∙ H2O, 0.1 g; CaCl2∙ 2H2O, 0.1 g; HCl 10 mL.
Batch fermentations were performed in a 5-L BIOFLO & CELLIGEN 310 bioreactor (New Brunswick. Scientific Co., Edison, NJ) with 3 L culture medium containing 90 g/L of D-glucose and 300 mL seed culture. The fermenter was continuously aerated through a 0.2 μm membrane filter at a flow rate of 1 vvm (air volume∙ working volume-1∙ minute-1). The temperature was maintained at 37°C. The pH was controlled at 6.5 ± 0.1 by the automatic feeding of 5 N NaOH. Foaming was controlled by the addition of Antifoam 289 (Sigma, St. Louis, MO). For biomass composition analysis, the cells were cultured in a medium in which only yeast extract was removed from the culture medium. All fermentations were performed at least three times independently, and the representative results are shown in Figures.
The concentrations of D-glucose and metabolites, including 2,3-BD, formic acid, ethanol, acetic acid, lactic acid, succinic acid, and acetoin, were determined by a high-performance liquid chromatography equipped with UV/VIS and RI detectors (Agilent 1260 series, Agilent Technologies, Waldbronn, Germany). An Aminiex HPX-87H column (300 mm × 7.8 mm, Bio-Rad, Hercules, CA) was isocratically eluted with 0.01 N H2SO4 at 80°C and a flow rate of 0.6 mL/min. The OD600 was measured using a UV–vis spectrophotometry (DR5000, Hach Company, CO) to monitor cell growth. Cell concentration, DCW per liter of culture broth, was calculated from the pre-determined standard curve relating OD600 to DCW (1 OD600 = 0.3877 ± 0.0136 g DCW∙ L-1). DCW was determined by filtering the culture broth through a 0.45 μm membrane and washing it twice with an equal volume of deionized distilled water. The filtered cells were then dried at 80 ± 5°C overnight and cooled to room temperature in a desiccator prior to weighting.
Open reading frame
Flux balance analysis
Kyoto encyclopedia of genes and genomes
Air volume• working volume-1• minute-1
Gram dry cell weight
Revolutions per minute
This work was supported by the Industrial Strategic Technology Development Program (No. 10035241) funded by the Ministry of Knowledge Economy (MKE, Korea).
- Volch M, Jansen NB, Ladisch MR, Tsao GT, Narayan R, Rodwell VW: 2, 3-butanediol.Comprehensive biotechnology. Volume 3. Edited by: Murray M-Y. 1985, 932-947. New York: Pergamon Press, 1
- Ji XJ, Huang H, Ouyang PK: Microbial 2, 3-butanediol production: a state-of-the-art review. Biotechnol Adv. 2011, 29: 351-364. 10.1016/j.biotechadv.2011.01.007.View Article
- Celinska E, Grajek W: Biotechnological production of 2, 3-butanediol–current state and prospects. Biotechnol Adv. 2009, 27: 715-725. 10.1016/j.biotechadv.2009.05.002.View Article
- Zeng AP, Sabra W: Microbial production of diols as platform chemicals: recent progresses. Curr Opin Biotechnol. 2011, 22: 749-757. 10.1016/j.copbio.2011.05.005.View Article
- Cho JH, Rathnasingh C, Song H, Chung BW, Lee HJ, Seung D: Fermentation and evaluation of Klebsiella pneumoniae and K. oxytoca on the production of 2, 3-butanediol. Bioprocess Biosyst Eng. 2012, 35: 1081-1088. 10.1007/s00449-012-0691-7.View Article
- Qureshi N, Cheryan M: Production of 2, 3-butanediol by Klebsiella oxytoca. Appl Microbiol Biotechnol. 1989, 30: 440-443.
- Ji XJ, Nie ZK, Huang H, Ren LJ, Peng C, Ouyang PK: Elimination of carbon catabolite repression in Klebsiella oxytoca for efficient 2, 3-butanediol production from glucose-xylose mixtures. Appl Microbiol Biotechnol. 2011, 89: 1119-1125. 10.1007/s00253-010-2940-5.View Article
- Ji XJ, Huang H, Li S, Du J, Lian M: Enhanced 2, 3-butanediol production by altering the mixed acid fermentation pathway in Klebsiella oxytoca. Biotechnol Lett. 2008, 30: 731-734. 10.1007/s10529-007-9599-8.View Article
- Ji XJ, Huang H, Zhu JG, Ren LJ, Nie ZK, Du J, Li S: Engineering Klebsiella oxytoca for efficient 2, 3-butanediol production through insertional inactivation of acetaldehyde dehydrogenase gene. Appl Microbiol Biotechnol. 2010, 85: 1751-1758. 10.1007/s00253-009-2222-2.View Article
- Ji XJ, Huang H, Du J, Zhu JG, Ren LJ, Hu N, Li S: Enhanced 2, 3-butanediol production by Klebsiella oxytoca using a two-stage agitation speed control strategy. Bioresour Technol. 2009, 100: 3410-3414. 10.1016/j.biortech.2009.02.031.View Article
- Shin SH, Kim S, Kim JY, Lee S, Um Y, Oh MK, Kim YR, Lee J, Yang KS: Complete genome sequence of Klebsiella oxytoca KCTC 1686, used in production of 2, 3-butanediol. J Bacteriol. 2012, 194: 2371-2372. 10.1128/JB.00026-12.View Article
- Liao TL, Lin AC, Chen E, Huang TW, Liu YM, Chang YH, Lai JF, Lauderdale TL, Wang JT, Chang SC: Complete genome sequence of Klebsiella oxytoca E718, a New Delhi metallo-beta-lactamase-1-producing nosocomial strain. J Bacteriol. 2012, 194: 5454- 10.1128/JB.01216-12.View Article
- Yang G, Tian J, Li J: Fermentation of 1, 3-propanediol by a lactate deficient mutant of Klebsiella oxytoca under microaerobic conditions. Appl Microbiol Biotechnol. 2007, 73: 1017-1024.View Article
- Champluvier B, Decallonne J, Rouxhet PG: Influence of sugar source (lactose, glucose, galactose) on 2, 3-butanediol production by Klebsiella oxytoca NRRL-B199. Arch Microbiol. 1989, 152: 411-414. 10.1007/BF00446920.View Article
- Jansen NB, Flickinger MC, Tsao GT: Production of 2, 3-butanediol from D-xylose by Klebsiella oxytoca ATCC 8724. Biotechnol Bioeng. 1984, 26: 362-369. 10.1002/bit.260260411.View Article
- Cao N, Xia Y, Gong CS, Tsao GT: Production of 2, 3-butanediol from pretreated corn cob by Klebsiella oxytoca in the presence of fungal cellulase. Appl Biochem Biotechnol. 1997, 63–65: 129-139.View Article
- Homman T, Tag C, Biebl H, Deckwer WD, Schink B: Fermentation of glycerol to 1, 3-propanediol by Klebsiella and Citrobacter strains. Appl Microbiol Biotechnol. 1990, 33: 121-126.
- Ji XJ, Huang H, Du J, Zhu JG, Ren LJ, Li S, Nie ZK: Development of an industrial medium for economical 2, 3-butanediol production through co-fermentation of glucose and xylose by Klebsiella oxytoca. Bioresour Technol. 2009, 100: 5214-5218. 10.1016/j.biortech.2009.05.036.View Article
- Nie ZK, Ji XJ, Huang H, Du J, Li ZY, Qu L, Zhang Q, Ouyang PK: An effective and simplified fed-batch strategy for improved 2, 3-butanediol production by Klebsiella oxytoca. Appl Biochem Biotechnol. 2011, 163: 946-953. 10.1007/s12010-010-9098-6.View Article
- Garg SK, Jain A: Fermentative production of 2, 3-butanediol: a review. Bioresour Technol. 1995, 51: 103-109. 10.1016/0960-8524(94)00136-O.View Article
- Beronio PB, Tsao GT: Optimization of 2, 3-butanediol production by Klebsiella oxytoca through oxygen transfer rate control. Biotechnol Bioeng. 1993, 42: 1263-1269. 10.1002/bit.260421102.View Article
- Beronio PB, Tsao GT: An energetic model for oxygen-limited metabolism. Biotechnol Bioeng. 1993, 42: 1270-1276. 10.1002/bit.260421103.View Article
- Qureshi N, Cheryan M: Effects of aeration on 2, 3-butanediol production from glucose by Klebsiella oxytoca. J Ferment Bioeng. 1989, 67: 415-418. 10.1016/0922-338X(89)90147-5.View Article
- Zhang G, Yang G, Wang X, Guo Q, Li Y, Li J: Influence of blocking of 2, 3-butanediol pathway on glycerol metabolism for 1, 3-propanediol production by Klebsiella oxytoca. Appl Biochem Biotechnol. 2012, 168: 116-128. 10.1007/s12010-011-9363-3.View Article
- Zhang Y, Huang Z, Du C, Li Y, Cao Z: Introduction of an NADH regeneration system into Klebsiella oxytoca leads to an enhanced oxidative and reductive metabolism of glycerol. Metab Eng. 2009, 11: 101-106. 10.1016/j.ymben.2008.11.001.View Article
- Lee JW, Na D, Park JM, Lee J, Choi S, Lee SY: Systems metabolic engineering of microorganisms for natural and non-natural chemicals. Nat Chem Biol. 2012, 8: 536-546. 10.1038/nchembio.970.View Article
- Lee JW, Kim TY, Jang YS, Choi S, Lee SY: Systems metabolic engineering for chemicals and materials. Trends Biotechnol. 2011, 29: 370-378. 10.1016/j.tibtech.2011.04.001.View Article
- Park JH, Lee SY: Towards systems metabolic engineering of microorganisms for amino acid production. Curr Opin Biotechnol. 2008, 19: 454-460. 10.1016/j.copbio.2008.08.007.View Article
- Furusawa C, Horinouchi T, Hirasawa T, Shimizu H: Systems metabolic engineering: the creation of microbial cell factories by rational metabolic design and evolution. Adv Biochem Eng Biotechnol. 2012, in press
- Becker J, Wittmann C: Systems and synthetic metabolic engineering for amino acid production - the heartbeat of industrial strain development. Curr Opin Biotechnol. 2012, 23: 718-726. 10.1016/j.copbio.2011.12.025.View Article
- Becker J, Zelder O, Hafner S, Schroder H, Wittmann C: From zero to hero–design-based systems metabolic engineering of Corynebacterium glutamicum for L-lysine production. Metab Eng. 2011, 13: 159-168. 10.1016/j.ymben.2011.01.003.View Article
- Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40: D109-D114. 10.1093/nar/gkr988.View Article
- Caspi R, Altman T, Dreher K, Fulcher CA, Subhraveti P, Keseler IM, Kothari A, Krummenacker M, Latendresse M, Mueller LA: The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res. 2012, 40: D742-D753. 10.1093/nar/gkr1014.View Article
- Blomqvist K, Nikkola M, Lehtovaara P, Suihko ML, Airaksinen U, Straby KB, Knowles JK, Penttila ME: Characterization of the genes of the 2, 3-butanediol operons from Klebsiella terrigena and Enterobacter aerogenes. J Bacteriol. 1993, 175: 1392-1404.
- Park JM, Kim TY, Lee SY: Constraints-based genome-scale metabolic simulation for systems metabolic engineering. Biotechnol Adv. 2009, 27: 979-988. 10.1016/j.biotechadv.2009.05.019.View Article
- Orth JD, Conrad TM, Na J, Lerman JA, Nam H, Feist AM, Palsson BO: A comprehensive genome-scale reconstruction of Escherichia coli metabolism–2011. Mol Syst Biol. 2011, 7: 535-View Article
- Schuetz R, Kuepfer L, Sauer U: Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli. Mol Syst Biol. 2007, 3: 119-View Article
- Saier MH, Tran CV, Barabote RD: TCDB: the transporter classification database for membrane transport protein analyses and information. Nucleic Acids Res. 2006, 34: D181-D186. 10.1093/nar/gkj001.View Article
- Ren Q, Chen K, Paulsen IT: TransportDB: a comprehensive database resource for cytoplasmic membrane transport systems and outer membrane channels. Nucleic Acids Res. 2007, 35: D274-D279. 10.1093/nar/gkl925.View Article
- Gombert AK, Nielsen J: Mathematical modelling of metabolism. Curr Opin Biotechnol. 2000, 11: 180-186. 10.1016/S0958-1669(00)00079-3.View Article
- Edwards JS, Ramakrishna R, Schilling CH, Palsson BO: Metabolic flux balance analysis. Metabolic engineering. Edited by: Lee SY, Papoutsakis ET. 1999, 13-57. New York: Marcel Dekker Inc
- Stephanopoulos G, Aristidou AA, Nielsen JH: Metabolic engineering: principles and methodologies. 1998, San Diego: Academic
- Varma A, Palsson BO: Metabolic flux balancing: basic concepts, scientific and practical use. Biotechnology. 1994, 12: 994-998. 10.1038/nbt1094-994.View Article
- Delgado J, Liao JC: Inverse flux analysis for reduction of acetate excretion in Escherichia coli. Biotechnol Prog. 1997, 13: 361-367. 10.1021/bp970047x.View Article
- Lee KH, Park JH, Kim TY, Kim HU, Lee SY: Systems metabolic engineering of Escherichia coli for L-threonine production. Mol Syst Biol. 2007, 3: 149-View Article
- Jung YK, Kim TY, Park SJ, Lee SY: Metabolic engineering of Escherichia coli for the production of polylactic acid and its copolymers. Biotechnol Bioeng. 2010, 105: 161-171. 10.1002/bit.22548.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.