Metabolic engineering of Synechococcus elongatus PCC 7942 for improvement of 1,3-propanediol and glycerol production based on in silico simulation of metabolic flux distribution
Microbial Cell Factories volume 16, Article number: 212 (2017)
Production directly from carbon dioxide by engineered cyanobacteria is one of the promising technologies for sustainable future. Previously, we have successfully achieved 1,3-propanediol (1,3-PDO) production using Synechococcus elongatus PCC 7942 with a synthetic metabolic pathway. The strain into which the synthetic metabolic pathway was introduced produced 3.48 mM (0.265 g/L) 1,3-PDO and 14.3 mM (1.32 g/L) glycerol during 20 days of incubation. In this study, the productivities of 1,3-PDO were improved by gene disruption selected by screening with in silico simulation.
First, a stoichiometric metabolic model was applied to prediction of cellular metabolic flux distribution in a 1,3-PDO-producing strain of S. elongatus PCC 7942. A genome-scale model of S. elongatus PCC 7942 constructed by Knoop was modified by the addition of a synthetic metabolic pathway for 1,3-PDO production. Next, the metabolic flux distribution predicted by metabolic flux balance analysis (FBA) was used for in silico simulation of gene disruption. As a result of gene disruption simulation, NADPH dehydrogenase 1 (NDH-1) complexes were found by screening to be the most promising candidates for disruption to improve 1,3-PDO production. The effect of disruption of the gene encoding a subunit of the NDH-1 complex was evaluated in the 1,3-PDO-producing strain.
Results and Conclusions
During 20 days of incubation, the ndhF1-null 1,3-PDO-producing strain showed the highest titers: 4.44 mM (0.338 g/L) 1,3-PDO and 30.3 mM (2.79 g/L) glycerol. In this study, we successfully improved 1,3-PDO productivity on the basis of in silico simulation of gene disruption.
Direct production of chemicals and fuels from carbon dioxide by cyanobacteria is one of the promising technologies to reduce carbon dioxide emissions and the dependence on fossil fuels [1, 2]. Compared with biological production using heterotrophic microorganisms, biological production using cyanobacteria can avoid various complicated processes for harvesting and decomposition of biomass. The introduction of a synthetic metabolic pathway composed of genes from other organisms into cyanobacteria is a key technology for production of various chemicals directly from carbon dioxide. Some of these valuable chemicals include isobutyraldehyde, isobutanol [3, 4], ethanol , 2,3-butanediol [6, 7], 2-methyl-1-butanol , 1-butanol [9, 10], acetone [11, 12], 3-hydroxybutyrate , ethylene [14, 15], isoprene [16, 17], 1,2-propanediol (1,2-PDO) , and glycerol [19, 20]. In addition, we successfully developed Synechococcus elongatus PCC 7942 to produce 1,3-propanediol (1,3-PDO) by introduction of a synthetic metabolic pathway . 1,3-PDO is a valuable chemical widely used for manufacture of polymers, paint, solvents, and antifreeze agents . In the introduced synthetic metabolic pathway, dihydroxyacetone phosphate (DHAP) in the Calvin cycle is converted to 1,3-PDO via glycerol (Fig. 1). The strain into which the synthetic pathway was introduced was found to produce 3.79 mM (0.29 g/L) 1,3-PDO and 12.6 mM (1.16 g/L) glycerol during 14 days of incubation . Although chemical production directly from carbon dioxide is a key advantage of cyanobacteria compared with other heterotrophic organisms, even the ethanol production by cyanobacteria showing the highest titer (5.5 g/L)  among various chemicals has not been increased over 10 g/L. Therefore, the improvements of productivity are key issues to make production by engineered cyanobacteria practical and cost effective.
In silico simulation is an powerful methodology to improve productivity with low experimental costs. The stoichiometric metabolic model not requiring kinetic parameters of metabolic reactions is comparatively applicable to analysis of whole-cell metabolism. Flux balance analysis (FBA) can predict metabolic flux distribution under constraint conditions . FBA using the genome-scale metabolic model (GSM) has been applied to improvement, for example, of succinic acid , threonine , valine , and lycopene  production by Escherichia coli. For screening of effective candidates for gene manipulation, various strategies have been established for prediction of the effects of a gene knockout, knockdown, and overexpression . Shastri and Morgan  were the first to construct a stoichiometric metabolic model of cyanobacteria composed of the central metabolism and photosystem of Synechocystis sp. PCC 6803. So far, the metabolic models of some cyanobacteria used for bioproduction (Synechocystis sp. PCC 6803, Synechococcus sp. PCC 7002, Synechococcus elongatus PCC 7942, Arthrospira platensis, Cyanothece sp. ATCC 51142) have been constructed on the basis of an annotated genome sequence [30,31,32]. Anfelt et al.  reported that introduction of the phosphoketolase pathway increases cellular acetyl-CoA level and butanol productivity in Synechocystis sp. PCC 6803. FBA was used for validation of this genetic modification to increase theoretical productivity of butanol. Furthermore, some simulations based on FBA have predicted various metabolic modifications (knockout and overexpression) to increase theoretical production by cyanobacteria, for example, of biofuels (1-butanol, 1-octanol, limonene, ethanol) and succinic acid by Synechocystis sp. PCC 6803 [33,34,35] and biofuels (butanol, ethanol) by Synechococcus sp. PCC 7002 . Nonetheless, there are only a few reports that these predicted strategies for productivity improvement applied experimentally .
In this study, we applied FBA to prediction of metabolic flux distribution in a 1,3-PDO-producing strain of S. elongatus PCC 7942 constructed in our previous work . The candidates for the gene knockout to increase 1,3-PDO productivity were estimated by FBA with the minimization of metabolic adjustment (MOMA) algorithm . The positive effect of the gene knockout on 1,3-PDO productivity after screening by in silico simulation was experimentally validated in the 1,3-PDO-producing strain.
Chemicals and reagents
All the chemicals were purchased from Wako Pure Chemical Industries, Ltd. (Osaka, Japan) unless specified otherwise. Restriction enzymes, phosphatase (New England Biolabs; Ipswich, MA, USA), ligase (Rapid DNA Ligation Kit, Roche; Mannheim, Germany), and DNA polymerase (KOD Plus Neo DNA polymerase, TOYOBO Co., Ltd.; Osaka, Japan) were used for cloning. Oligonucleotides were synthesized by Life Technologies Japan, Ltd. (Tokyo, Japan).
A modified BG11 medium supplemented with 20 mM HEPES–NaOH (pH 7.5) for pH stabilization was used for cultivation of S. elongatus PCC 7942 (Life Technologies Corporation; Carlsbad, CA, USA). Hereafter, this medium is referred to as the BG11 medium. The composition of the BG11 medium used in this study was the same as previously reported . Antibiotics (10 μg/mL kanamycin, 20 µg/mL spectinomycin, and 5 µg/mL chloramphenicol) were added to the BG11 medium when appropriate.
The incubation conditions for 1,3-PDO production were the same as those described previously . Cell density (OD730) was measured using Infinite 200 PRO (TECAN; Männedorf, Switzerland). For determination of dry cell weight (DCW), cells in culture were harvested on the Ominipore™ membrane filter (JHWP04700; Merck Millipore; Darmstadt, Germany) by filtration under reduced pressure. Cell harvest filters were dried at 65 °C for 12 h. The weight differences before and after cell harvesting were determined as DCW in culture.
Plasmids and primers used in this study are listed in Additional file 1: Tables S1 and S2, respectively. To construct a plasmid for gene disruption, the upper and lower regions of a target gene were amplified for homologous recombination by PCR. The primer sets for the upstream region of ndhD1, ndhD2, ndhD3, ndhD4, and ndhF1 were T2711–T2712, T3189–T3190, T2723–T2724, T2719–T2720, and T3208–3209, respectively. The primer sets for the downstream region of ndhD1, ndhD2, ndhD3, ndhD4, and ndhF1 were T3187–T3188, T3191–T3192, T2787–T2726, T2721–2722, and T2713–T2714, respectively. The amplified upstream and downstream regions except for ndhF1 were digested by KpnI–HindIII and HindIII–BamHI, respectively, and ligated into KpnI–BamHI sites of pTA703. Only for ndhF1, the upstream and downstream regions were digested with SalI–HindIII and HindIII–BamHI, respectively, and ligated into SalI–BamHI sites of pTA703. A chloramphenicol resistance cassette was amplified from pZA31-luc (Expressys; Ruelzheim, Germany) by PCR using primers T2128–T2129. The amplified fragment digested with HindIII was ligated into the HindIII site of each plasmid harboring upstream and downstream regions of a target gene, thus creating plasmids pTA1811 (ndhD1 disruption), pTA1812 (ndhD2 disruption), pTA1641 (ndhD3 disruption), pTA1640 (ndhD4 disruption), and pTA1828 (ndhF1 disruption).
The constructed plasmids (pTA1811, 1812, 1641, 1640, and 1828) were used for gene disruption in TA1297 (wild-type strain), TA2984 (1,3-PDO-producing strain), and TA3800 (glycerol-producing strain). TA2984 and TA3800 were constructed in our previous studies [21, 38]. The strains used in this study are listed in Additional file 1: Table S1.
Culture supernatants were centrifuged (20,000×g, 10 min) and filtered using Minisart RC4 (Sartorius; Goettingen, Germany). Concentrations of glycerol and 1,3-PDO in filtered samples were determined using a high-performance liquid chromatograph (LC-20AD, Shimadzu; Kyoto, Japan) equipped with an autosampler (AOC-20; Shimadzu), a SUGAR SP-G guard column (Shodex; Tokyo, Japan), a SUGAR SP0810 column (Shodex), and a refractive index detector (RID-10A; Shimadzu). Operating conditions were as follows: injection volume, 5 µL; mobile phase, MilliQ water; flow rate, 1.0 mL/min in isocratic mode; column temperature, 80 °C.
In silico simulation of metabolic flux distribution
FBA for TA2984 (1,3-PDO-producing strain)
The stoichiometric metabolic model based on the annotated genome sequence of S. elongatus PCC 7942 published by Knoop  (hereafter called GSM_7942 in this study) served as a basic metabolic model of the wild-type strain in this study. GSM_7942 is composed of 595 genes, 666 reactions, and 533 metabolites. The information on the introduced synthetic metabolic pathway for 1,3-PDO production (8 genes, 6 reactions, and 3 metabolites) was added to GSM_7942, resulting in a metabolic model of TA2984 (hereafter called GSM_7942_PDO in this study). GSM_7942_PDO is composed of 603 genes, 672 reactions, and 536 metabolites. GSM_7942_PDO was used to perform FBA using the COBRA toolbox 2.05  on Matlab (Mathworks Inc., Natick, MA).
Gene disruption simulation in TA2984 with the MOMA algorithm
This simulation was performed on GSM_7942_PDO with scripts provided in COBRA toolbox 2.05 using Gurobi 6.03 (Gurobi optimization Inc., USA) as a solver.
Respiratory and photosynthetic activities
Chlorophyll concentrations were measured by extraction in methanol . Oxygen consumption (respiration) and evolution (photosynthesis) from whole cells were measured with a Clark-type electrode (OXYT-1, Hansatech Instruments, Norfolk, UK) at a chlorophyll concentration of 10 µg/mL. The temperature of the electrode chamber was maintained at 30°C, and the chamber contents were continuously stirred. Saturated light (1000 μE/m2/s) and growth light (100 μE/m2/s) provided by a digital LED source (PD3-5024-4-EI, CCS, Kyoto, Japan) were used for measurement of photosynthetic activity. Photon flux density was measured with IKS-27 (Koito; Tokyo, Japan). Respiratory activity was similarly measured in complete darkness.
Measurement of ATP content
ATP extraction and measurement were basically performed by the method of Wang et al. . The appropriate number of cells at growth phase (corresponding to 1-mL culture at OD730 of 1.0) was collected, and ATP was extracted from pelleted cells by resuspension in a 1% (w/v) trichloroacetic acid solution. After centrifugation at 20,000×g for 10 min at 4 °C, the supernatant was neutralized with 1 M Tris–acetate buffer (pH 7.8). The concentration of ATP in the neutralized sample was determined using the ATP determination kit (A22066, Molecular Probes Inc., Oregon, USA).
Measurement of intracellular metabolites
Intracellular metabolites of TA2984 (1,3-PDO-producing strain) and TA4058 (ndhF1-null 1,3-PDO-producing strain) were measured by LC-QqQ-MS and GC-Q-MS analyses. Sample preparation and analysis were performed by the method described by Soma and Hanai .
Measurement of glycogen content
Glycogen extraction was basically performed by the method of Suzuki et al. . Cells at early stationary phase (days 10) were used for glycogen measurement. The extracted glycogen was mixed with glucoamylase (Oriental Yeast Co., Ltd.; Tokyo, Japan) in 2.5 mM Na-acetate (pH 5.0) and incubated at 40 °C for 1 h. After centrifugation, glucose generated from glycogen was measured using BF-7 (Oji Scientific Instruments, Hyogo, Japan) equipping glucose electrode (ED07-0003, Oji Scientific Instruments).
Results and discussion
FBA of the 1,3-PDO-producing strain
We previously constructed a 1,3-PDO-producing strain, designated as TA2984, by introduction of a synthetic metabolic pathway. A four-step reaction of the introduced pathway converts DHAP, one of the metabolites in the Calvin cycle, into 1,3-PDO via glycerol (Fig. 1). Construction of the stoichiometric metabolic model based on the annotated genome sequence of S. elongatus PCC 7942 was first reported by Triana et al. . In the same year, the comparable scale of a metabolic model of this strain was published by Knoop et al. . We used Knoop’s model composed of 595 genes, 666 reactions, and 533 metabolites as a basic metabolic model of the wild-type strain of S. elongatus PCC 7942 (hereafter called GSM_7942 in this study). The additional information on the introduced synthetic metabolic pathway for 1,3-PDO production (8 genes, 6 reactions, and 3 metabolites) was added to GSM_7942, resulting in a metabolic model of the 1,3-PDO-producing strain composed of 603 genes, 672 reactions, and 536 metabolites (hereafter called GSM_7942_PDO in this study).
GSM_7942_PDO was used for prediction of metabolic flux distribution of TA2984 into which the synthetic metabolic pathway for 1,3-PDO was introduced . To predict metabolic flux distribution using a stoichiometric metabolic model, input fluxes to the cell and output fluxes from the cell were needed. For TA2984, the former represents the carbon fixation rate and photon uptake rate, and the latter represents the biomass production rate (corresponding to the specific growth rate) and specific production rates of 1,3-PDO and glycerol. To determine output fluxes, TA2984 was incubated under the same conditions as in our previous report . As a result, TA2984 produced 3.48 mM (0.265 g/L) 1,3-PDO and 14.3 mM (1.32 g/L) glycerol directly from carbon dioxide during 20 days of incubation (Fig. 2). Although the output fluxes were observed in the experimental data shown in Fig. 2, it was difficult to experimentally determine the input fluxes (carbon fixation rate and photon uptake rate). Therefore, GSM_7942_PDO was used for calculation of the input fluxes capable of satisfying the output fluxes observed experimentally. Because the assumption of a steady state of metabolism is one of the key constraints in FBA , the data from the exponential growth phase (days 2–4) were applied to the prediction. First, the specific growth rate and the specific production rates of 1,3-PDO and glycerol in the exponential growth phase were calculated to be 0.0250/h, 0.0273 mmol/(g DCW)/h, and 0.0515 mmol/(g DCW)/h, respectively. Next, the input fluxes of carbon fixation and photon use were calculated and found to be 1.35 and 19.0 mmol/(g DCW)/h, respectively. These input fluxes were calculated by the method of two-step optimization strategy of Shastri and Morgan  with some modifications. These calculated input fluxes and experimentally obtained output fluxes were used for prediction of the metabolic flux distribution of TA2984 in the exponential growth phase. Under the conditions restricted by input fluxes of carbon fixation and photon use, the maximized specific growth rate corresponded to 0.0290/h without 1,3-PDO and glycerol production. A two-dimensional solution space between the specific growth rate and the specific production rate of 1,3-PDO was then described (Fig. 3a). The experimental data shown in Fig. 2 were plotted in the described solution space.
Simulation of gene disruption in TA2984 (1,3-PDO-producing strain)
The metabolic flux distribution of TA2984 using FBA was applied to the gene disruption simulation based on the MOMA algorithm . This is one of the algorithms used for simulation of sequential gene deletion assuming that the cell adapts to a mutation to minimize the changes of metabolic flux distribution by the gene deletion. The estimated candidates for gene deletion based on this algorithm successfully increased lycopene production in E. coli . Because the objective of our simulation was to screen promising candidates for gene disruption for improvement of 1,3-PDO productivity, we set the specific production rate of 1,3-PDO obtained experimentally [0.0250 mmol/(g DCW)/h] as the lower limit. The metabolic flux distribution corresponding to each of the 595 gene disruptions in GSM_7942_PDO except for the 8 genes of 1,3-PDO producing synthetic metabolic pathway was simulated and the simulated results on 581 strains with single-gene disruption were plotted according to the estimated specific growth rate and estimated specific production rate of 1,3-PDO (Fig. 3b). The simulated results on the 581 strains with single-gene disruption were roughly subdivided into four groups. These four groups are presented as different symbols in Fig. 3b.
The first group (crosses in the figure) was composed of 74 strains that could not satisfy the constraint conditions and lower limit of 1,3-PDO production. The second group (triangles) included 254 strains showing a higher specific production rate of 1,3-PDO, but a null value of the specific growth rate. This may also mean that each gene disrupted in the 254 strains was essential for growth. It was likely that the disruption of any of the 328 genes classified into these two groups was not applicable to improvement of 1,3-PDO productivity. The metabolic flux distributions of the other 267 strains with single-gene disruption could show a certain specific growth rate and higher specific production rate of 1,3-PDO as compared with the experimental data. We screened the 267 candidates for more desirable ones using the product of the specific growth rate and specific production rate as an indicator. Out of the 267 candidates, 231 showed a lower product of the specific growth rate and specific production rate and were classified into the third group (squares).
Finally, 36 genes were classified into the fourth group (circles) as promising candidates to be disrupted for increasing 1,3-PDO productivity (Additional file 1: Table S3). Only 10 patterns of metabolic flux distribution were obtained from 36 strains in the fourth group. First, in GSM_7942, the annotated-gene information was defined by a gene–protein reaction (GPR) association . Under the definition of the GPR association, each disruption of a subunit of a protein complex indicates the same flux distribution. In this study, 22 genes involved in NADPH dehydrogenase-1 (NDH-1) complexes and four genes involved in cytochrome c oxidase were classified into the fourth group. Second, two cases indicating the same metabolic flux distribution were observed in different two-gene disruptions (Synpcc7942_0191 and 1501, and Synpcc7942_2079 and 2080). These genes were annotated as genes encoding an enzyme catalyzing a serial reaction without any blanching pathway. The other six metabolic flux distributions were examined by single-gene disruption. Based on the value of the product of the specific growth rate and specific production, the disruption of 22 genes involved in the NDH-1 complex was found by screening as the most promising candidate for improvement of 1,3-PDO production among the 10 patterns of the fourth group.
Disruption of genes involved in NDH-1 complexes
Simulated metabolic flux distribution for each strain by gene disruption indicated that the elimination of the NDH-1 complex is the most effective candidate for improvement of 1,3-PDO productivity. The estimated metabolic flux distributions with and without NDH-1 complex elimination were compared (Fig. 4). Each value with an arrow represents the flux [mmol/(g DCW)/h] calculated from experimental data (Fig. 4a) and predicted one of the NDH-1 complex-null strains using the MOMA algorithm (Fig. 4b). In all simulations of gene disruption, the flux of carbon fixation [1.35 mmol/(g DCW)/h] was assumed to not be changed by any gene disruptions. The stoichiometry of ribulose-bisphosphate carboxylase/oxygenase (RubisCO) reaction in GSM_7942 means that 1.97 molecules of 3-phosphoglycerate (3-PG) were produced from 1 molecule of ribulose-1,5-bisphosphate (RuBP), considering photorespiration . Therefore, the flux of 3-phosphoglycerate (3-PG) production corresponded to 2.66 mmol/(g DCW)/h. By means of elimination of the NDH-1 complex, some metabolic fluxes around the 1,3-PDO-producing pathway were predicted to be changed. First, the blanching ratio of the 3-PG-consuming flux into the glycolysis (2-phosphoglycerate; 2-PG) and the Calvin cycle (1,3-bisphosphoglycerate; BPG) was changed. NDH-1 complex elimination increased the flux of fixed carbon into the Calvin cycle. Second, each flux of DHAP production from GAP (glyceraldehyde 3-phosphate) and that of DHAP consumption for S7P (sedoheptulose 7-phosphate) or FBP (fructose-1,6-bisphosphate) synthesis was increased and decreased, respectively, resulting in the increase in the flux into the 1,3-PDO-producing pathway.
NDH-1 complexes in cyanobacteria are involved in several cellular functions (respiration, cyclic electron flow, and carbon dioxide acquisition) and are subdivided into four variants (NDH-1L, L′, MS, and MS′) by function and subunit composition according to reverse-genetic studies . The phenotypic analyses of gene knockout mutants revealed the functional differences among NDH-1 complexes, indicating that NDH-1L and L′ are involved in respiration and cyclic electron flow, whereas NDH-1MS and MS′ participate in carbon dioxide acquisition [48, 49]. Because NdhB is one of the common subunits for all four variants, it was expected that ndhB disruption would be desirable to represent the metabolic flux distribution observed in the gene disruption simulation (Fig. 4b, Additional file 1: Table S3). Nonetheless, all functions of NDH-1 complexes were drastically decreased in the ndhB-null mutant, in particular, this mutant could not grow under ambient air conditions [50,51,52]. Indeed, ndhB could not be completely disrupted in the wild-type strain and TA2984 (1,3-PDO-producing strain) in our experiment. On the other hand, it was reported that four mutants in which each gene encoding the unique subunit in four NDH-1 complexes is disrupted can grow under photoautotrophic and ambient air conditions . Therefore, we tried to disrupt the genes encoding these unique subunits in TA2984.
In contrast to the strain with ndhB disruption, the gene knockout strains of ndhD1 (Synpcc7942_1976, unique subunit gene for NDH-1L), ndhD2 (Synpcc7942_1439, that for NDH-1L′), ndhD3 (Synpcc7942_2092, that for NDH-MS), and ndhD4 (Synpcc7942_0609, that for NDH-MS′) based on the 1,3-PDO-producing strain (TA2984) were successfully constructed (Additional file 1: Table S1). These four strains were applied to 1,3-PDO production (Additional file 2: Figure S1). All the strains showed a growth defect as compared with TA2984, and the growth defect was more severe in three strains except for TA3672 (ndhD3-null strain; Additional file 2: Figure S1A). Among the four strains, TA4021 (ndhD1-null strain) showed the highest titers of 1,3-PDO (4.04 mM, 0.307 g/L) and glycerol (19.6 mM, 1.81 g/L) during 20 days of incubation (Additional file 2: Figure S1B, C). These titers of 1,3-PDO and glycerol shown by TA4021 were 1.16- and 1.38-fold higher than those of TA2984 (Fig. 2). Table 1 shows output fluxes (specific production rates or specific growth rates) for the wild type and the knockout strains observed in the exponential growth phase (days 2–4) and those calculated by the simulation. This strain revealed a higher specific production rate of 1,3-PDO in the exponential growth phase than TA2984 did (Table 1). The titer of 1,3-PDO produced by TA4022 (ndhD2-null strain) did not increase (Additional file 2: Figure S1B) while a significantly improved production rate of 1,3-PDO was observed in this strain (Table 1). The specific production rate of 1,3-PDO by TA3672 (ndhD3-null strain) was the lowest among the four gene knockout strains and TA2984 (Table 1). On the other hand, TA3671 (ndhD4-null strain) showed the highest production rate of 1,3-PDO (Table 1), but the titer of 1,3-PDO produced by this strain was lower only than the titer of TA2984 (Additional file 2: Figure S1).
NdhF1 is a unique subunit commonly present in two NDH-1 complexes: NDH-1L and L′ . It was expected that ndhF1 disruption would lead to the additive phenotype of ndhD1 and ndhD2 knockouts. Other studies on gene knockout mutants of ndhD1 and ndhF1 showed that the NDH-1L complex is missing in each mutant [53, 54]. Just as ndhD1-D4, the ndhF1-null 1,3-PDO-producing strain designated as TA4058 was successfully constructed (Additional file 1: Table S1). The growth of TA4058 was more defective than that of TA4021 and TA4022, indicating that the ndhF1 disruption affects cellular function more as compared with each single-gene disruption of ndhD1 or ndhD2 (Fig. 5a and Additional file 2: Figure S1A). During 20 days of incubation, the titers of 1,3-PDO and glycerol produced by TA4058 corresponded to 4.44 mM (0.338 g/L) and 30.3 mM (2.79 g/L), respectively (Fig. 5b, c). Both titers were the highest among the other gene knockout strains and the parental strain for 1,3-PDO production. These results indicated that the disruption of NDH-1 complexes (NDH-1L and L′), which have a function in respiration and cyclic electron flow, positively affected 1,3-PDO and glycerol production. The specific production rate of 1,3-PDO by TA4058 in the exponential growth phase (days 2–4) was significantly enhanced and the highest among all the constructed strains (Table 1). Unfortunately, this value [0.0497 mmol/(g DCW)/h] was inferior to that of simulation of the NDH-1 complex elimination [0.0852 mmol/(g DCW)/h]. Nonetheless, the sum of the specific production rates of 1,3-PDO and glycerol meaning the utilization rate of DHAP for the 1,3-PDO-producing pathway in TA4058 [0.213 mmol/(g DCW)/h] was surprisingly higher than that in simulation [0.162 mmol/(g DCW)/h]. The strain with higher total DHAP utilization showed a lower ratio of 1,3-PDO production to total DHAP utilization, indicating that the later part of the 1,3-PDO-producing pathway converting glycerol to 1,3-PDO is the bottleneck of effective production (Table 1). This result is consistent with findings of our recent study on 1,3-PDO production improvement by promoter exchange in the later part of the 1,3-PDO producing pathway . It is expected that the strategy involving a combination of promoter exchange and NDH-1 complex disruption is a promising method for further improvement of 1,3-PDO productivity.
Previously, we have constructed a glycerol-producing strain (TA3800) into which only the pathway converting DHAP to glycerol was introduced . The effect of ndhF1 disruption was also evaluated in TA3800. As in the 1,3-PDO producing strain, the ndhF1 disruption resulted in a growth defect and improvement of the production titer (Additional file 2: Figure S2). The growth defect caused by the ndhF1 disruption in the 1,3-PDO-producing strain was more significant than that in the wild-type strain and glycerol-producing strain (Additional file 2: Figure S2A–C). The highest titer of glycerol produced by the ndhF1-null glycerol-producing strain designated as TA4059 reached 36.7 mM (3.38 g/L; Additional file 2: Figure S2D). This value was almost equal to the sum of the titers of 1,3-PDO and glycerol (34.7 mM) shown by TA4058 (Fig. 5) and about twofold higher than that of TA3800 (Additional file 2: Figure S2E). Although the reason for the significant growth defect caused by ndhF1 disruption in the 1,3-PDO-producing strain was not clear, glycerol production was successfully stimulated by the ndhF1 disruption, as was 1,3-PDO production. The titer of glycerol produced by TA4059 was the highest among other chemicals (glycerol, 1,3-PDO, and 1,2-PDO) produced from DHAP by engineered cyanobacteria [18,19,20, 38] and this titer corresponded to high titer compared to other chemicals produced not from DHAP by engineered cyanobacteria, such as ethanol , sucrose , 2,3-butanediol , and lactate . Although various chemicals have been produced by engineered cyanobacteria as summarized in some reviews [1, 2], the research resulting from the titer of product over 3 g/L directly from carbon dioxide seems to be still rare.
Metabolism changes caused by ndhF1 disruption
The disruption of ndhF1 successfully increased metabolic flux of DHAP into the 1,3-PDO-producing pathway as expected on the basis of the in silico simulation (Table 1) and improved the titers of 1,3-PDO and glycerol production (Fig. 5 and Additional file 2: Figure S2). This result suggests that in silico simulation is an effective method for metabolic improvement of bioproduction in cyanobacteria. To identify the causes of improved productivity, the metabolic changes caused by the ndhF1 disruption were analyzed. It was expected that the ndhF1 disruption affected to respiration and cyclic electron flow . Therefore, at first, the respiration activities of TA4058 (ndhF1-null 1,3-PDO-producing strain) and TA2984 (parental strain) were measured. As the result, TA4058 showed lower respiratory activity in the dark and almost the same photosynthetic activity under saturated light conditions (1000 μE/m2/s) as compared with TA2984 (Additional file 2: Figure S3). This finding was consistent with the phenotype of a ndhF1-null strain of Synechococcus sp. PCC 7002 observed in another study . Under growth light conditions (100 μE/m2/s), the photosynthetic activity of TA4058 was lower than that of TA2984, being consistent with the growth defect of TA4058 (Fig. 5).
The physiological role of cyclic electron flow is known to increase the ATP/NADPH ratio produced by photosynthetic electron flow . So, next, ATP contents in TA4058 and TA2984 were measured. TA4058 showed lower cellular ATP content as compared with TA2984 (Additional file 2: Figure S4). In the respiratory chain, the reducing power (NADH and FADH2) mainly produced by glycolysis and the TCA cycle is used for generating a proton gradient and H2O from O2. The generated proton gradient functions as a driving force of ATP synthase for ATP production. Whole or partial loss of cyclic electron flow and respiratory activity under the influence of the ndhF1 disruption may decrease cellular ATP content.
The intracellular metabolites involved in a reaction with ATP seemed to reflect the ATP deficiency in TA4058. For example, ATP is consumed in reactions producing RuBP and ADP-glucose, and is produced in the reaction producing pyruvate. In ATP-consuming reactions, substrates (ribulose-5-phosphate: Ru5P and glucose-1-phosphate: G1P) and products (RuBP and ADP-glucose) were detected at high and low concentrations, respectively. Conversely, lower concentration of a substrate (phosphoenolpyruvate: PEP) and higher concentration of a product (pyruvate) were observed in an ATP-producing reaction. In particular, the metabolites around the glycogen synthesis pathway seemed to be affected by the ATP deficiency. Fructose-6-phosphate (F6P) is a metabolite of the branching point for glycogen synthesis from the Calvin cycle and is sequentially converted to glucose-6-phosphate (G6P), glucose-1-phosphate (G1P), and ADP-glucose. In this glycogen synthesis pathway, significant accumulation (more than fivefold) was observed for metabolites except for ADP-glucose produced with ATP consumption. This result suggests that glycogen synthesis is repressed in TA4058 caused by ATP deficiency. Because other studies showed that higher accumulation of glycogen is observed in stationary phase [44, 59], the glycogen concentrations of TA2984 and TA4058 were measured at 10 days (Additional file 2: Figure S6). As expected, the intracellular glycogen concentration of TA4058 was about half of that of TA2984 (Additional file 2: Figure S6A).
It is well known that the respiratory activity is strongly related to catabolic metabolism (glycolysis and TCA cycle) in heterotrophic organisms. We hypothesized that the repressed respiratory activity can decrease the metabolic flux in catabolic metabolism generating reducing power in photosynthetic organisms, too. The decreasing of 3-PG-consuming flux into glycolysis by the ndhF1 disruption in the simulation result (Fig. 4) was consistent with the hypothesis. Quantification of intracellular metabolites also supported this hypothesis (Additional file 2: Figure S5). Most of the metabolites in the Calvin cycle showed higher intracellular concentration in TA4058 than in TA2984. On the other hand, the metabolites in the glycolysis and TCA cycle were mostly lower in TA4058 compared with TA2984.
Not only greater flux of 3-PG into the Calvin cycle but also repression of glycogen synthesis may increase the productivity values of 1,3-PDO and glycerol. These results indicate that the productivity improvement was caused by the metabolism adapted to energy imbalance induced by the ndhF1 disruption. Hasunuma et al.  successfully improved the carbon fixation rate in Synechocystis sp. PCC 6803 by overexpression of the flv3 gene involved in NADPH consumption. The strategy of acceleration of NADPH consumption also showed the energy balance to change cellular metabolism of cyanobacteria. In the present study, the metabolism changes intended to adapt to the energy imbalance (ATP deficiency) were applied to chemical production (1,3-PDO and glycerol). Taking into consideration of the points mentioned above, ndhF1 disruption is a good strategy for enhancement of production of other chemicals (for instance, 1,2-PDO and sucrose) from metabolites in the Calvin cycle in cyanobacteria [18, 38, 55].
The simulation result of the metabolic flux distribution of the ndhF1-null and wild type strain (Fig. 4) showed the 3-PG-consuming flux to 2-PG or TCA cycle was decreased by the ndhF1 disruption. This result suggested that the ndhF1-null strain was not suitable for the chemical productions from pyruvate or acetyl-CoA. Yoshikawa et al. , however, showed the ndhF1-null PCC6803 strain with ethanol producing pathway was produced 145% more ethanol (0.132 g/L) from pyruvate than the wild type strain with the same pathway (0.191 g/L). It is difficult to explain clearly the reason why the suggestion based on our simulation result is different from their result, but one possible reason would be the driving force by ATP generation by the reaction from PEP to pyruvate to compensate the ATP deficiency in the ndhF1-null strain.
Recently, Broddrick et al.  constructed GEM of metabolism in PCC7942 composed of 785 genes. This model was constructed by using the information about essential genes revealed by the random barcode transposon site sequencing and physiological data for photoautotrophic metabolism. Using this model or the larger models than that we used, there are possibilities to screen out the other candidate genes for knockout to improve 1,3-PDO productivity.
In this study, the productivity of 1,3-PDO and glycerol in S. elongatus PCC 7942 were successfully improved by gene manipulation based on in silico simulation of gene disruptions. The best candidate for improvement of 1,3-PDO productivity was found to be the elimination of NDH-1 complexes involved in several cellular functions: respiration, cyclic electron flow, and carbon dioxide uptake (Additional file 1: Table S3). Genes encoding NDH-1 complexes were actually disrupted in TA2984 (parental strain for 1,3-PDO production). The effective improvement of productivity was achieved by disruption of the NDH-1 complexes that are related to respiration and cyclic electron flow (Additional file 2: Figure S1). Finally, TA4058 (ndhF1-null 1,3-PDO-producing strain) was found to produce 1.28-fold more 1,3-PDO (4.44 mM, 0.338 g/L) and 2.13-fold more glycerol (30.3 mM, 2.79 g/L) than TA2984 did (Fig. 5). The metabolism adapted to the energy imbalance caused by the ndhF1 disruption may enhance the 3-PG-consuming flux into the Calvin cycle and can be suitable for 1,3-PDO and glycerol production (Additional file 2: Figure S5). This result opens up exciting opportunities for applying in silico simulation to metabolic improvement of cyanobacteria and is suggestive of the feasibility of metabolic changes focused on the energy balance.
Angermayr SA, Rovira AG, Hellingwerf KJ. Metabolic engineering of cyanobacteria for the synthesis of commodity products. Trends Microbiol. 2015;33:352–61.
Oliver NJ, Rabinovitch-Deere CA, Carroll AL, Nozzi NE, Case AE, Atsumi S. Cyanobacterial metabolic engineering for biofuel and chemical production. Curr Opin Chem Biol. 2016;35:43–50.
Atsumi S, Higashide W, Liao JC. Direct photosynthetic recycling of carbon dioxide to isobutyraldehyde. Nat Biotechnol. 2009;27:1177–80.
Li X, Shen CR, Liao JC. Isobutanol production as an alternative metabolic sink to rescue the growth deficiency of the glycogen mutant of Synechococcus elongatus PCC 7942. Photosynth Res. 2014;120:301–10.
Gao Z, Zhao H, Li Z, Tan X, Lu X. Photosynthetic production of ethanol from carbon dioxide in genetically engineered cyanobacteria. Energy Environ Sci. 2012;5:9857–65.
Oliver JWK, Machado IM, Yoneda H, Atsumi S. Cyanobacterial conversion of carbon dioxide to 2,3-butanediol. Proc Natl Acad Sci USA. 2013;110:1249–54.
Oliver JWK, Machado IM, Yoneda H, Atsumi S. Combinatorial optimization of cyanobacterial 2,3-butanediol production. Metab Eng. 2014;22:76–82.
Shen CR, Liao JC. Photosynthetic production of 2-methyl-1-butanol from CO2 in cyanobacterium Synechococcus elongatus PCC 7942 and characterization of the native acetohydroxyacid synthase. Energy Environ Sci. 2012;5:9574–83.
Anfelt J, Kaczmarzyk D, Shabestary K, Renberg B, Rockberg J, Nielsen J, Uhlén M, Hudson EP. Genetic and nutrient modulation of acetyl-CoA levels in Synechocystis for n-butanol production. Microb Cell Fact. 2015;14:167–78.
Lan EI, Ro SY, Liao JC. Oxygen-tolerant coenzyme A-acylating aldehyde dehydrogenase facilitates efficient photosynthetic n-butanol biosynthesis in cyanobacteria. Energy Environ Sci. 2013;6:2672–81.
Chwa JW, Kim WJ, Sim SJ, Um Y, Woo HM. Engineering of a modular and synthetic phosphoketolase pathway for photosynthetic production of acetone from CO2 in Synechococcus elongatus PCC 7942 under light and aerobic condition. Plant Biotechnol J. 2016;14:1768–76.
Zhou J, Zhang H, Zhang Y, Li Y, Ma Y. Designing and creating a modularized synthetic pathway in cyanobacterium Synechocystis enables production of acetone from carbon dioxide. Metab Eng. 2012;14:394–400.
Wang B, Pugh S, Nielsen DR, Zhang W, Meldrum DR. Engineering cyanobacteria for photosynthetic production of 3-hydroxybutyrate directly from CO2. Metab Eng. 2013;16:68–77.
Lee T-C, Xiong W, Paddock T, Carrieri D, Chang I-F, Chiu H-F, Ungerer J, Juo S-HH, Maness P-C. Engineered xylose utilization enhances bio-products productivity in the cyanobacterium Synechocystis sp. PCC 6803. Metab Eng. 2015;30:179–89.
Ungerer J, Tao L, Davis M, Ghirardi M, Maness PC, Yu J. Sustained photosynthetic conversion of CO2 to ethylene in recombinant cyanobacterium Synechocystis 6803. Energy Environ Sci. 2012;5:8998–9006.
Bentley FC, Melis Anastasios. Diffusion-based process for carbon dioxide uptake and isoprene emission in gaseous/aqueous two-phase photobioreactors by photosynthetic microorganisms. Biotechnol Bioeng. 2012;109:100–9.
Gao X, Gao F, Liu D, Zhang H, Nie X, Yang C. Engineering the methylerythritol phosphate pathway in cyanobacteria for photosynthetic isoprene production from CO2. Energy Environ Sci. 2016;9:1400.
Li H, Liao JC. Engineering a cyanobacterium as the catalyst for the photosynthetic conversion of CO2 to 1,2-propanediol. Microb Cell Fact. 2013;12:4.
Savakis P, Tan X, Du W, Branco dos Santos F, Lu X, Hellingwerf KJ. Photosynthetic production of glycerol by a recombinant cyanobacterium. J Bacteriol. 2015;195:46–51.
Wang Y, Tao F, Ni J, Li C, Xu P. Production of C3 platform chemicals from CO2 by genetically engineered cyanobacteria. Green Chem. 2015;17:3100–10.
Hirokawa Y, Maki Y, Hanai T. Cyanobacterial production of 1,3-propanediol directly from carbon dioxide using a synthetic metabolic pathway. Metab Eng. 2016;34:97–103.
Saxena RK, Anand P, Saran S, Isar J. Microbial production of 1,3-propanediol: recent developments and emerging opportunities. Biotechnol Adv. 2009;27:895–913.
Orth JD, Thiele I, Palsson BØ. What is flux balance analysis? Nat Biotechnol. 2010;28:245–8.
Hong SH, Lee SY. Metabolic flux analysis for succinic acid production by recombinant Escherichia coli with amplified malic enzyme activity. Biotechnol Bioeng. 2001;74:89–95.
Lee KH, Park JH, Kim TY, Kim HU, Lee SY. Systems metabolic engineering of Escherichia coli for l-threonine production. Biol: Mol Syst; 2007. p. 3.
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.
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.
Yen JY, Tanniche I, Fischer AK, Gillaspy GE, Bevan DR, Senger RS. Designing metabolic engineering strategies with genome-scale metabolic flux modeling. Adv Genomics Genet. 2015;5:93–105.
Shastri AA, Morgan JA. Flux balance analysis of photoautotrophic metabolism. Biotechnol Prog. 2005;21:1617–26.
Baroukh C, Muñoz-Tamayo R, Bernard O, Steyer J-P. Mathematical modeling of unicellular microalgae and cyanobacteria metabolism for biofuel production. Curr Opin Biotechnol. 2015;33:198–205.
Triana J, Montagud A, Siurana M, Fuente D, Urchueguía A, Gamermann D, Torres J, Tena J, Fernández de Córdoba P, Urchueguía JF. Generation and evaluation of a genome-scale metabolic network model of Synechococcus elongatus PCC 7942. Metabolites. 2014;4:680–98.
Yoshikawa K, Aikawa S, Kojima Y, Toya Y, Furusawa C, Kondo A, Shimizu H. Construction of a genome-scale metabolic model of Arthrospira platensis NIES-39 and metabolic design for cyanobacterial bioproduction. PLoS ONE. 2015;7:e0144430.
Shabestary K, Hudson EP. Computational metabolic engineering strategies for growth-coupled biofuel production by Synechocystis. Metab Eng Commun. 2016;3:216–26.
Shirai T, Osanai T, Kondo A. Designing intracellular metabolism for production of target compounds by introducing a heterologous metabolic reaction based on a Synechocystis sp. PCC 6803 genome-scale model. Microb Cell Fact. 2016;15:13–9.
Yoshikawa K, Toya Y, Shimizu H. Metabolic engineering of Synechocystis sp. PCC6803 for enhanced ethanol production based on flux balance analysis. Bioprocess Biosyst Eng. 2017. https://doi.org/10.1007/s00449-017-1744-8.
Hendry JI, Prasannan CB, Joshi A, Dasgupta S, Wangikar PP. Metabolic model of Synechococcus sp. PCC 7002: prediction of flux distribution and network modification for enhanced biofuel production. Bioresour Technol. 2016;213:190–7.
Segrè D, Vitkup D, Church GM. Analysis of optimality in natural and perturbed metabolic networks. Proc Natl Acad Sci USA. 2002;99:15112–7.
Hirokawa Y, Maki Y, Hanai T. Improvement of 1,3-propanediol production using an engineered cyanobacterium, Synechococcus elongatus by optimization of the gene expression level of a synthetic metabolic pathway and production conditions. Metab Eng. 2017;39:192–9.
Knoop H. Integrative analyse des cyanobakteriellen stoffwechsels. PhD thesis. 2014. https://doi.org/10.18452/17055.
Schellenberger J, Que R, Fleming RM, Thiele I, Orth JD, Feist AM, Zielinski DC, Bordbar A, Lewis NE, Rahmanian S, Kang J, Hyduke DR, Palsson BØ. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011;6:1290–307.
Grimme LH, Boardman NK. Photochemical activities of a particle fraction P1 obtained from the green algae Chlorella fucsa. Biochem Biophys Res Commun. 1972;49:1617–23.
Wang X, Liu W, Xin C, Zheng Y, Cheng Y, Sun S, Li R, Zhu X-G, Dai SY, Rentzepis PM, Yuan JS. Enhanced limonene production in cyanobacteria reveals photosynthesis limitations. Proc Natl Acad Sci USA. 2016;113:14225–30.
Soma Y, Hanai T. Self-induced metabolic state switching by a tunable cell density sensor for microbial isopropanol production. Metab Eng. 2015;30:7–15.
Suzuki E, Ohkawa H, Moriya K, Matsubara T, Nagaike Y, Iwasaki I, Fujiwara S, Tsuzuki M, Nakamura Y. Carbohydrate metabolism in mutants of the cyanobacterium Synechococcus elongatus PCC 7942 defective in glycogen synthesis. Appl Environ Microbiol. 2010;76:3153–9.
Alper H, Jin Y-S, Moxley JF, Stephanopoulos G. Identifying gene tragets for the metabolic engineering of lycopene biosynthesis in Escherichia coli. Metab Eng. 2005;7:155–64.
Thiele I, Palsson BØ. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010;5:93–121.
Battchikova N, Eisenhut M, Aro E-M. Cyanobacterial NDH-1 complexes: novel insights and remaining puzzles. Biochim Biophys Acta. 2011;1807:935–44.
Ohkawa H, Pakrasi HB, Ogawa T. Two types of functionally distinct NAD(P)H dehydrogenases in Synechocystis sp. strain PCC 6803. J Biol Chem. 2000;275:31630–4.
Shibata M, Ohkawa H, Kaneko T, Fukuzawa H, Tabata S, Kaplan A, Ogawa T. Distinct constitutive and low-CO2-induced CO2 uptake systems in cyanobacteria: genes involved and their phylogenetic relationship with homologous genes in other organisms. Proc Natl Acad Sci USA. 2001;98:11789–94.
Marco E, Ohad N, Schwarz R, Lieman-Hurwitz J, Gabay C, Kaplan A. High CO2 concentration alleviates the block in photosynthetic electron transport in an ndhB-inactivated mutant of Synechococcus sp. PCC 7942. Plant Physiol. 1993;101:1047–53.
Mi H, Endo T, Schreiber U, Ogawa T, Asada K. Electron donation from cyclic and respiratory flows to the photosynthetic intersystem chain is mediated by pyridine nucleotide dehydrogenase in the cyanobacterium Synechocystis PCC 6803. Plant Cell Physiol. 1992;33:1233–7.
Ogawa T. A gene homologous to the subunit-2 gene of NADH dehydrogenase is essential to inorganic carbon transport of Synechocystis PCC 6803. Proc Natl Acad Sci USA. 1991;88:4275–9.
Zhang P, Battchikova N, Jansen T, Appel J, Ogawa T, Aro E-M. Expression and functional roles of the two distinct NDH-1 complexes and the carbon acquisition complex NdhD3/F3/CupA/Sll1735 in Synechocystis sp. PCC 6803. Plant Cell. 2004;16:3326–40.
Wulfhorst H, Franken LE, Wessinghage T, Boekema EJ, Nowaczyk MM. The 5 kDa protein NdhP is essential for stable NDH-1L assembly in Thermosynechococcus elongatus. PLoS ONE. 2014;9:e103584.
Ducat DC, Avelar-Rivas JA, Way JC, Silver PA. Rerouting carbon flux to enhance photosynthetic productivity. Appl Environ Microbiol. 2012;78:2660–8.
Angermayr SA, Hellingwerf KJ. On the use of metabolic control analysis in the optimization of cyanobacterial biosolar cell factories. J Phys Chem B. 2013;117:11169–75.
Schluchter WM, Zhao J, Bryant DA. Isolation and characterization of the ndhF gene of Synechococcus sp. strain PCC 7002 and initial characterization of an interposon mutant. J Bacteriol. 1993;175:3343–52.
Rochaix J-D. Reprint of: regulation of photosynthetic electron transport. Biochem Biophys Acta. 2011;1807:876–86.
Song K, Tan X, Liang Y, Lu X. The potential of Synechococcus elongatus UTEX 2973 for sugar feedstock production. Appl Microbiol Biotechnol. 2016;100:7865–75.
Hasunuma T, Matsuda M, Senga Y, Aikawa S, Toyoshima M, Shimakawa G, Miyake C, Kondo A. Overexpression of flv3 improves photosynthesis in the cyanobacterium Synechocystis sp. PCC 6803 by enhancement of alternative electron flow. Biotechnol Biofuels. 2014;7:493–503.
Broddrick JT, Rubin RE, Welkie DG, Du N, Mih N, Diamond S, Lee JJ, Golden SS, Palsson BO. Unique attributes of cyanobacterial metabolism revealed by improved genome-scale metabolic modeling and essential gene analysis. Proc Natl Acad Sci USA. 2016;113:E8344–53.
YH, SM, and TH designed the experiments. YH and SM carried out all experiments and SM, HH, and FM carried out theoretical analysis and YH drafted the figures and manuscript. TH edited the manuscript. All authors read and approved the final manuscript.
This study was financially supported by the Core Research of Evolutional Science and Technology program (CREST) from the Japan Science and Technology Agency (JST).
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used in the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
The Core Research of Evolutional Science and Technology program (CREST) from the Japan Science and Technology Agency (JST).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hirokawa, Y., Matsuo, S., Hamada, H. et al. Metabolic engineering of Synechococcus elongatus PCC 7942 for improvement of 1,3-propanediol and glycerol production based on in silico simulation of metabolic flux distribution. Microb Cell Fact 16, 212 (2017). https://doi.org/10.1186/s12934-017-0824-4
- Flux balance analysis
- Synthetic metabolic pathway