Abstract
Isoprenoids are diverse natural compounds, which have various applications as pharmaceuticals, fragrances, and solvents. The low yield of isoprenoids in plants makes them difficult for cost-effective production, and chemical synthesis of complex isoprenoids is impractical. Microbial production of isoprenoids has been considered as a promising approach to increase the yield. In this study, we engineered the model cyanobacterium Synechocystis sp. PCC 6803 for sustainable production of a commercially valuable isoprenoid, limonene. Limonene synthases from the plants Mentha spicata and Citrus limon were expressed in cyanobacteria for limonene production. Production of limonene was two-fold higher with limonene synthase from M. spicata than that from C. limon. To enhance isoprenoid production, computational strain design was conducted by applying the OptForce strain design algorithm on Synechocystis 6803. Based on the metabolic interventions suggested by this algorithm, genes (ribose 5-phosphate isomerase and ribulose 5-phosphate 3-epimerase) in the pentose phosphate pathway were overexpressed, and a geranyl diphosphate synthase from the plant Abies grandis was expressed to optimize the limonene biosynthetic pathway. The optimized strain produced 6.7âmg/L of limonene, a 2.3-fold improvement in productivity. Thus, this study presents a feasible strategy to engineer cyanobacteria for photosynthetic production of isoprenoids.
Similar content being viewed by others
Introduction
Recent studies have demonstrated the potential of using cyanobacteria as biological platforms to produce fuels and high-value chemicals1,2. Harnessing solar energy using the photosynthetic apparatus, atmospheric CO2 is fixed into sugars, which can be further converted to desired products by engineered cyanobacteria. Due to the recent development of genetic tools for model cyanobacteria3, expression of heterologous genes and pathways has become more feasible, thus facilitating the construction of engineered cyanobacteria for biotechnological applications. In this study, we engineered the model cyanobacterium Synechocystis sp. PCC 6803 (hereafter, Synechocystis 6803) for production of a commercially valuable isoprenoid, limonene.
Isoprenoids are one of the most diverse groups of natural products, with more than 55,000 compounds4. Isoprenoids have multiple commercial applications, including natural pharmaceuticals, nutraceuticals, solvents, and perfume components5,6. To date, commercially-used isoprenoids are mainly extracted from plants, but the low quantities of these naturally-produced chemicals have become an impediment for cost-effective production. Successful microbial production of valuable isoprenoids by engineered yeast and E. coli have been demonstrated7,8, whereas fewer researchers have studied production of isoprenoids by cyanobacteria. To improve photosynthetic production of isoprenoids, optimization of isoprenoid biosynthetic pathways in cyanobacteria is needed using metabolic engineering coupled with computational approaches.
Limonene is a 10-carbon isoprenoid produced by plants. (R)-limonene has a characteristic fragrance of orange, and commonly exists in the rinds of citrus fruits. It is commercially used as a fragrance in perfumes or a solvent in cleaning products. (S)-Limonene is a precursor for the biosynthesis of (S)-menthol, which is the major component of mint. Recently, limonene has been evaluated as a âdrop-inâ replacement for diesel9 and jet fuels10. The fully hydrogenated form of limonene was used as a diesel additive, exhibiting similar chemical properties compared to diesel fuel9. Moreover, the physical properties of limonene, such as viscosity, freezing point, and boiling point, are highly comparable to aviation fuel Jet A-110.
Cyanobacteria use the methylerythritol 4-phosphate (MEP) pathway to produce isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP), which are the building blocks for isoprenoid biosynthesis. The MEP pathway is a seven-step pathway that starts with glyceraldehyde 3-phosphate (GAP) and pyruvate, and ends with IPP and DMAPP. Further, IPP and DMAPP undergo a series of head-to-tail condensations to produce diphosphate substrates, which are then converted to isoprenoids by isoprenoid synthases. To increase isoprenoid production, the amounts of IPP and DMAPP need to be enhanced by increasing the carbon flux toward the MEP pathway.
Attempts have been made to engineer the MEP pathway for improving cyanobacterial limonene production. However, production titers are extremely low compared to other compounds such as ethanol11, butanol12, and free fatty acid13. Genes involving in the bottlenecks of the MEP pathway were overexpressed in Synechocystis 680314. The recombinant strain showed a 1.4-fold increase of limonene, and the final titer reached 1âmg/L after 30-day cultivation14. In addition, researchers used similar strategies to engineer the MEP pathway in the nitrogen-fixing cyanobacterium Anabaena sp. PCC 7120 for production of limonene15. The limonene yield increased up to 6.8-fold. However, the final titer remained low (0.5âmg/L over 12-day incubation)15.
A previous in vitro study suggested that isoprenoid production in Synechocystis 6803 is stimulated by compounds in the pentose phosphate (PP) pathway but not by substrates in the MEP pathway16. Using Synechocystis 6803 cell extracts, isoprenoid biosynthesis was significantly improved by supplying xylulose 5-phosphate (X5P) in the PP pathway, whereas providing substrates (GAP, pyruvate, and MEP) in the MEP pathway showed lower stimulation of isoprenoid production16. These results indicated a connection between the PP pathway and isoprenoid production in Synechocystis 6803.
In addition to experimental engineering approaches, computational strain design techniques can be useful to develop non-intuitive genetic interventions to achieve the desired level of production of a particular bioproduct. To this end, the OptForce procedure17 first characterizes the wild-type strain in the form of reaction flux ranges by utilizing the 13C MFA (Metabolic Flux Analysis) flux estimations as additional regulations. OptForce then contrasts the wild-type flux ranges with those in the overproducing phenotype. As a result, the algorithm identifies a set of genetic interventions (i.e., up/down-regulations and deletions) that must happen in the metabolic reaction network for a desired level of yield. Finally, OptForce pinpoints the minimal interventions (from these changes) that are directly related to achieving the desired yield. These strategies can then be tested in an experimental setting.
In this work, we engineered Synechocystis 6803 for photosynthetic limonene production (Fig. 1). To construct limonene-producing strains, genes encoding limonene synthase (lims) from Mentha spicata and Citrus limon were introduced into Synechocystis 6803. For generating computation-driven non-intuitive strain engineering strategies, we applied the OptForce algorithm17 on the genome-scale Synechocystis 6803 model iSyn73118 and also utilized 13C MFA flux estimations19 under photosynthetic wild-type condition. OptForce predicted the up-regulation of two PP pathway genes, ribose 5-phosphate isomerase (rpi) and ribulose 5-phosphate 3-epimerase (rpe), in limonene-producing strains in order to divert the carbon flux toward limonene production. Furthermore, based on the prediction made by OptForce to further improve limonene production, a geranyl diphosphate synthase (gpps) from Abies grandis was expressed to optimize the limonene production pathway. The final recombinant strain led to a 2.3-fold improvement in yield, producing 6.7âmg/L of limonene in 7 days. The metabolic engineering strategies used in this study demonstrate the feasibility of increasing limonene production in Synechocystis 6803 and can be applied to phototrophic production of other high-value isoprenoids.
Results
Engineering Synechocystis 6803 for production of limonene
Limonene is a C10 cyclic isoprenoid converted from geranyl diphosphate (GPP). Due to the complex nature of carbocation rearrangement from GPP to limonene, limonene synthase produces not only limonene but also other monoterpenes such as bicyclic α-pinene and acyclic mycene20. To avoid the production of other unwanted byproducts, we chose limonene synthases which have the highest specificity for limonene production. Based on previous studies, limonene synthase from Citrus limon and Mentha spicata produce limonene of high purity. Expression of each of these limonene synthases in E. coli showed that the former produces 99% pure (R)-limonene21, and the latter generates 94% of (S)-limonene22. The coding sequences of lims were codon optimized for Synechocystis 6803, and the plastid targeting sequences were removed23,24. The truncated enzyme is known to have better catalytic activity than the native protein25. Genes were cloned into a pCC5.2 neutral-site-targeting plasmid and driven by the trc1O promoter for higher level expression of lims (Fig. 2A). Expression of an enhanced yellow fluorescent protein (EYFP) from the pCC5.2 endogenous plasmid is 8 to 14 times higher than that on the chromosome26.
When the lims was cloned into a suicide plasmid and transformed into E. coli, we found that the gene accumulated random mutations in the E. coli host, leading to changes in amino acid residues or truncated proteins. This was presumably because the lims product is toxic to E. coli cells. To introduce a lims without mutations into Synechocystis 6803, we circumvented the E. coli cloning step by first cloning the lims into the suicide plasmid via Gibson assembly, and used the assembled product as template for PCR to amplify the lims cassette flanked by upstream and downstream homologous sequences of the neutral site in pCC5.226. Subsequently, the PCR product was directly used for natural transformation into Synechocystis 6803. The lims was introduced into Synechocystis 6803 genome via double homologous recombination (Fig. 2A). DNA sequencing results showed that the lims has no mutation in Synechocystis 6803 (data not shown). Mutants were fully segregated after re-streaking the cells several times on BG-11 plates with antibiotics.
Limonene production by engineered Synechocystis 6803 was tested by incubating cultures for 7 days. Because of the volatility of limonene, a dodecane overlay was applied on cultures to collect limonene in the organic layer. It has been reported that over 99% of limonene escapes from the cyanobacterial cultures14, and covering an organic overlay on cultures had little influence on growth in cyanobacteria23. The limonene yield by the strain expressing lims from M. spicata was two-fold higher than that by the strain expressing lims from C. limon (Fig. 2B). These results suggest that the limonene synthase from M. spicata exhibited better catalytic activity in Synechocystis 6803, and hence, the strain was used for further engineering.
Computational modeling
The iSyn731 metabolic model of Synechocystis 680318 was used to perform the computational strain designs using the OptForce algorithm17 for overproduction of limonene. Based on the current understanding as reported in literature16,27, a connection between Calvin Benson Cycle (CBC)/PP pathway and MEP pathway (Fig. 3) was included in the iSyn731 model. By superimposing the photoautotrophic flux measurements19 of 31 reactions of central carbon metabolism including the CBC and PP pathways of Synechocystis 6803 onto the iSyn731 model, the phenotypic space of the base strain was defined. All simulations were performed for a basis of 100 millimoles of CO2 plus H2CO3 uptake and unlimited photon supply19. The uptake fluxes for the remaining metabolites present in the BG11 medium was set to -1,000 and the non-growth associated ATP maintenance was set at 8.39 mmole/gDW-h. In addition, the biomass flux was fixed at the optimal value subject to the experimental flux measurements19. The upper bound of the fluxes of the remaining reactions was set to 1,000 mmole/gDW-h, whereas the lower bound was set to zero and -1,000 mmole/gDW-h for irreversible and reversible reactions, respectively.
Similarly, the limonene overproducing phenotype was obtained by maximizing and minimizing each flux of the metabolic model iteratively subject to the network stoichiometry, uptake and medium conditions, regulatory constraints, and overproduction target. In this work, a minimum production yield of 85% of the theoretical maximum of limonene (i.e., 15.3 mmole/gDW-h) was set as the overproduction target, while the biomass flux was constrained to be at least 10% of its theoretical maximum (i.e., 0.021âhâ1) with the basis of 100 millimoles of carbon fixed (i.e., CO2 plus H2CO3). The remaining parameter values including medium conditions and regulatory constraints were the same as those in the wild-type. By contrasting the maximal range of flux variability between the wild-type strain and the over-producing strain to meet the pre-specified yield of limonene, OptForce was used to identify the minimal set of genetic interventions (i.e., deletions and up-/down-regulations). In order to first explore non-intuitive interventions, reactions from the MEP and isoprenoid biosynthesis pathways were not considered as the candidates for any form of intervention. Integer cuts were used to identify alternative optimal solutions (i.e., alternative genetic intervention choices) to achieve the minimum production yield of limonene as specified earlier. The termination criterion for the OptForce procedure was set as either meeting a production yield of at least 85% of the theoretical maximum for limonene or exceeding the maximum allowable number of reaction interventions (i.e., three). Note that the OptForce procedure works at the reaction level, which is why the set of genetic manipulations can subsequently be identified by using gene-protein-reaction (GPR) associations from the iSyn731 model. Thus, the OptForce procedure identified up-regulation of rpi and rpe as the best possible solution, which can lead up to limonene yield at 89% of its theoretical maximum (i.e., 16.02 mmole/gDW-h). By up-regulating these two genes, OptForce suggested to force more flux from the CBC/PP pathway toward MEP pathway that can ultimately increase the production yield of limonene (Fig. 3). Once the set of non-intuitive interventions was obtained, as a next step, it was logical to explore if their combination with any of the intuitive one(s) from the MEP and isoprenoid biosynthesis pathways could further improve the limonene production yield that was otherwise not possible to achieve individually (i.e., by the non-intuitive candidates or by the intuitive ones). With a target of a minimum production yield of 90% of the theoretical maximum of limonene, the OptForce procedure identified the up-regulation of gpps, rpe, and rpi that could lead the limonene production yield to 16.56 mmole/gDW.h (i.e., 92% of its theoretical maximum). Thus, the proposed interventions combined the amplification (i.e., push) of flux from the CBB/PP pathway to MEP pathway with a similar increase (i.e., pull) in the flux of the limonene synthesis. As reported in the literature28, this kind of push-and-pull strategy can achieve the desired level of production yield with minimal effects caused by feedback inhibition.
Genetic interventions of the PP pathway to improve limonene production
Based on the prediction of the OptForce procedure, up-regulation of rpi and rpe genes in the PP pathway increases the flux toward limonene production. To test this hypothesis, the rpi and rpe genes driven by the Synechocystis 6803 native rbcL promoter were expressed on a replicating plasmid in the limonene-producing strain, resulting in 1.3-fold increase in limonene yield (3.7âmg/L) after 7 days of cultivation (Fig. 4). Furthermore, we introduced a gene encoding a specific GPP synthase (GPPS) to optimize the limonene biosynthetic pathway. In Synechocystis 6803, formation of GPP is catalyzed by a farnesyl diphosphate (FPP) synthase, CrtE. It performs consecutive condensation of IPP with DMAPP, and only synthesizes GPP as an intermediate29. Although the PP pathway was engineered to stimulate the limonene yield, it is possible that the native isoprenoids pathway in Synechocystis 6803 provides insufficient GPP for limonene production since the flux is diverted toward FPP formation for pigment synthesis. In addition, OptForce also predicted an increase (i.e., from 89% to 92% of maximum theoretical limonene yield) when up-regulation of rpe and rpi was combined with the up-regulation of gpps. It was reported that the GPPS 2 from Abies grandis specifically produces GPP30. Expressing this specific gpps with lims, the limonene yield increased 1.4-fold (4.1âmg/L) (Fig. 4). Finally, coexpression of rpi, rpe, gpps and lims resulted in a remarkable (2.3-fold) enhancement in productivity (6.7âmg/L) (Fig. 4).
Pigment content in engineered Synechocystis 6803
Carotenoids and the phytol tail of chlorophyll, photosynthetic pigments, are derived from geranylgeranyl diphosphate (GGPP), a C20-intermediate for isoprenoid synthesis. Hence, production of limonene is expected to divert carbon flux away from pigment synthesis. To investigate the effect of limonene production on pigment content in engineered Synechocystis 6803, we extracted and quantified the chlorophyll and carotenoid contents. The chlorophyll content decreased over 30% in the gpps expression strains, whereas carotenoid levels were similar among the limonene-producing strains (Fig. 5). These results indicate that the specific GPPS diverts the carbon flux away from pigment synthesis.
Discussion
In this study, we combined metabolic engineering with model-driven strain design strategies to engineer Synechocystis 6803 for enhanced limonene production. To generate limonene-producing Synechocystis 6803, we first constructed a suicide plasmid26 to engineer the lims into the neutral site on the pCC5.2 endogenous plasmid via double homologous recombination. This is the first time that the endogenous plasmid of Synechocystis 6803 has been used for enhanced production for the purpose of metabolic engineering. Expression of a gene on the pCC5.2 plasmid leads to higher expression level than that on the chromosome as well as the RSF1010 self-replicating plasmid26. Furthermore, during the stationary phase of cell growth, the copy numbers of the endogenous plasmids (pCA2.4, pCB2.4, pCC5.2) in Synechocystis 6803 are 3 to 7 per chromosome31. Using the endogenous plasmid to express the lims gene driven by the constitutive promoter trc1O allows high expression level at the stationary phase, decoupling growth and production, and thus leading to higher levels of production of limonene.
The higher yield with limonene synthase from M. spicata than that from C. limon may be due to the difference in enzyme kinetics of LIMS. Unfortunately, the kinetic parameters (both K m and k cat) are only available for the enzyme from M. spicata 25. In addition, it may be attributed to different protein expression levels. Although the same promoter was used to control the lims from two plant species, protein expression may vary because of different mRNA sequences and codon usage. To date, the highest reported limonene productivity in cyanobacteria was achieved by engineered Synechococcus sp. PCC 700223. In their study, only a lims from M. spicata was overexpressed, and the yield was over 4âmg/L in 4 days23. Our results also suggested that the LIMS from M. spicata performed better in limonene production (Fig. 2B). The doubling time of Synechococcus 7002 is shorter than Synechocystis 680332. Thus, the higher limonene yield from Synechococcus 7002 may be due to its faster growth rate. A recent study engineered Synechococcus elongatus PCC 7942 to produce limonene, achieving a 100-fold improvement in productivity33. However, it should be noted that such significant increase is due to the low productivity of the original strain, which produced merely 8.5âμg/L/OD/d of limonene. The best producing strain in this study, with a lims (M. spicata) controlled by the pea plant psbA promoter, produced 2.5âmg/L limonene in 5 days33.
Previously, researchers have engineered Synechocystis 6803 for limonene production by overexpressing genes involved in the bottleneck steps of the MEP pathway14. It is known that enzymes 1-deoxy-D-xylulose-5-phosphate synthase (DXS) and isopentenyl diphosphate isomerase (IDI) catalyze the rate-limiting reactions in the MEP pathway34,35. With the introduction of an additional copy of endogenous dxs, idi, and gpps genes, the engineered Synechocystis 6803 produced 1.4-fold higher yield than that of the parent strain14. However, such improvement was less effective than the strategy used in the current study. As mentioned in the Results section, the endogenous gpps gene may not be suitable for enhancing the production of limonene. In addition, the MEP pathway is highly regulated at genetic and metabolic levels36. Expressing endogenous genes in the MEP pathway may be subject to native regulations, presenting a less effective engineering approach.
Instead of manipulating the MEP pathway, we took a systematic model-driven metabolic engineering approach for finding genetic interventions in order to increase the limonene production yield. As explained in the Materials and Methods section, the OptForce procedure finds the minimal interventions to reach a desired production target. To this end, we employed OptForce on our previously developed genome-scale model iSyn731 in order to âpushâ more flux to MEP pathway and also to create better âpullâ for limonene synthesis (Fig. 3). From this in silico analysis, by up-regulating rpe and rpi, the metabolite pool of X5âP was found to be increased that, eventually, led to increased flux through the connection between the CBC/PP pathway and the MEP pathway. In addition, up-regulation of gpps created an improved âpullâ for limonene synthesis. Thus, the combination of this push-and-pull mechanism was proposed to be the best strategy to improve limonene yield by circumventing additional regulations (e.g. feedback inhibition). Interestingly, the same rationale could be applied to engineer cyanobacteria to produce other isoprenoid compounds.
Expression of the specific gpps modestly increased the limonene titer (Fig. 4), whereas the cellular chlorophyll content was greatly influenced (Fig. 5). Synthesis of limonene and the phytol tail of chlorophyll requires the same precursors, IPP and DMAPP. Table 1 compares the changes in chlorophyll and limonene contents between the strain with lims only and the two gpps-expressing strains. Compared to the lims-expressing strain, additional expression of gpps resulted in similar level of decrease in chlorophyll content in both strains. However, expression of rpi and rpe genes further led to 1.5-fold higher limonene productivity in the gpps-expressing strains (122 vs. 81âμM/OD730 isoprene unit). This result suggests that up-regulation of rpi and rpe enhances the carbon flux toward limonene synthesis.
Our results showed that overexpressing the genes in the PP pathway led to improved limonene production, suggesting an unidentified connection between the PP pathway and isoprenoids biosynthesis (Fig. 4). Our observation is consistent with previous in vitro study using Synechocystis 6803 cell lysate16. However, the connection between the PP pathway and isoprenoids biosynthesis remains to be elucidated. It was first shown that in vitro isoprenoid production increased significantly by providing substrates in the PP pathway16, while a recent study showed that increased production of isoprenoids by PP pathway substrates does not occur through the MEP pathway37. By removing the terminal enzyme of the MEP pathway in Synechocystis 6803 cell lysate, isoprenoid synthesis still increased by substrates in the PP pathway37. Taken together, it is still unclear how the PP pathway and isoprenoid production are connected in Synechocystis 6803. While our results made a strong argument for this connection, further investigation needs to be conducted to explore the details in terms of chemical conversions and genes/enzymes associated. From the modeling context, these details sometimes do not make much of a difference if they only involve aggregating linear reaction steps.
Conclusions
In this study, we engineered the model cyanobacterium Synechocystis 6803 to produce the isoprenoid, limonene. We applied computational strain design by using the OptForce procedure to identify minimal genetic interventions for improving limonene yield. Based on the prediction, the rpi and rpe genes in the PP pathway were overexpressed, and a specific gpps was introduced to optimize the limonene biosynthetic pathway. The final engineered strain produced 6.7âmg/L of limonene, which is a 2.3-fold improvement in productivity. The approach that we demonstrated can be applied to engineer cyanobacteria to produce other valuable isoprenoids.
Materials and Methods
Chemicals and reagents
All chemicals were purchased from Sigma-Aldrich (St. Louis, MO, USA) unless otherwise specified. Phusion DNA polymerase were purchased from Thermo Fisher Scientific (Waltham, MA, USA).
Culture medium and condition
All strains were maintained in liquid BG-11 medium or on solid BG-11 plates with appropriate antibiotics at 30â°C continuous white light (50 μmoles photons mâ2sâ1).
DNA manipulations
Coding sequences of lims from Mentha spicata and Citrus limon were codon optimized for Synechocystis 6803 and synthesized by IDT (San Jose, CA, USA). The genes were cloned into a suicide plasmid, allowing gene insertion into the neutral site (NSP1) on the endogenous plasmid pCC5.2 in Synechocystis 680326. The constructed plasmids were directly used as templates for PCR to amplify a fragment which contains the lims and a kanamycin resistance cassette flanking by upstream and downstream homologous sequences of the NSP1 (Fig. 2A). The PCR product was then purified by DNA electrophoresis, and the linear DNA was transformed into Synechocystis 6803. The rpi, rpe, and gpps genes were cloned into a broad-host-range plasmid RSF1010 harboring a spectinomycin resistance cassette38. All the cloning works were done by Gibson isothermal DNA assembly method39.
Strains construction and transformation
The lims expression cassette was transformed into Synechocystis 6803 through homologous recombination. Cells at mid-log phase (OD730 of 0.4 to 0.6) were incubated with 600 ng of linear DNA overnight at 30â°C in the dark. Cells were then grown on BG-11 plates supplemented with 10âμg/mL of kanamycin for selection of transformants. Colonies were patched on BG-11 plates with 20âμg/mL of kanamycin for segregation. PCR was used to verify strain segregation. For the construction of rpi, rpe, and gpps expressing strains, self-replicating plasmids (600 ng per transformation) were transformed into the strain expressing lims. Transformants were selected by BG-11 plates with 2âμg/mL of spectinomycin and 5âμg/mL of kanamycin.
Limonene production by engineered cyanobacteria
Strains were inoculated in BG-11 medium with kanamycin (10âμg/mL) and spectinomycin (4âμg/mL) to mid-log phase at 30â°C with continuous white light (50 μmoles photons mâ2sâ1). Cells were collected by centrifugation at 7,000 x g, and washed by BG-11 medium to remove antibiotics. To test limonene production, the initial OD730 was adjusted to 0.34 (~0.5âg/L of biomass), and 50âmL of cell cultures were grown in 250-mL flasks at 30â°C with continuous white light (130 μmoles photons mâ2sâ1). A 10% (v/v) dodecane overlay was covered on top of cultures to trap evaporated limonene.
Quantification of limonene
Limonene samples were prepared by diluting 10âμL of dodecane overlay in 990âμL of ethyl acetate, and analyzed using a gas chromatography instrument with a flame ionization detector (Hewlett-Packard model 7890âA, Agilent Technologies, CA, USA) equipped with a 30âm DB5-MS column (J&W Scientific). The oven temperature program initiated at 60â°C, and increased at 12â°C/min to 300â°C. Limonene was quantified using a (R)-limonene standard.
Identification of engineering interventions via OptForce
We applied the OptForce algorithm17 on the genome-scale Synechocystis 6803 model iSyn73118. In order to characterize the wild-type phenotype, we utilized 13C MFA flux estimations19 under photosynthetic condition. Below is the step-by-step procedure that we followed:
Step 1
Identify the maximum biomass and limonene yields under photosynthetic condition.
Maximize v biomass or v ls
Subject to
Step 2
Characterize the wild-type phenotype
Maximize/Minimize \({v}_{j}\forall j\in \) reactions without experimental flux measurements
Subject to
Step 3
Characterize the limonene over-producing phenotype
Maximize/Minimizev j â j â 1, â¦â¦,m
Subject to
Step 4
Identify the MUST sets
In this step, fluxes ranging from step 2 and step 3 were compared to identify three different sets: reactions to be up-regulated (MUST U), down-regulated (MUST L), and deleted (MUST X).
Step 5
Identify the minimal engineering interventions
Maximize v j
(over MUST sets)
Subject to
Minimize v j
(over MUST sets)
Subject to
Here, S ij is the stoichiometric coefficient of metabolite i in reaction j and v j is the flux value of reaction j. Parameters v j,min and v j,max denote the minimum and maximum allowable fluxes for reaction j, respectively. V biomass and v ls represent biomass and limonene synthesis reactions under photosynthetic conditions, whereas v max biomass and v max ls represent the maximum theoretical yields of biomass and limonene under photosynthetic conditions. The minimal levels of biomass and the minimal target yield of limonene were set to be 10% of maximum biomass and 85% or 90% of maximum limonene yield, respectively. Finally, k represents the maximum number of interventions allowed.
Pigment content analysis
Cell cultures (1âmL) were collected by centrifugation at 16,000 x g for 7âmin, and the supernatants were removed. To extract pigments in Synechocystis, pre-cooled methanol (1âmL) was added to the pellets, and mixed thoroughly by pipetting and vortexing. Samples were incubated at 4â°C for 20âmins, and centrifuged at 16,000 x g for 7âmin. The supernatants were removed for a spectroph-otometer analysis to quantify the concentrations of carotenoids and chlorophyll. The following equations were used to calculate the pigment content: \(\text{chlorophyll}\,({\rm{\mu }}{\rm{g}}/{\rm{mL}})=(16.29\times {A}_{665})-(8.54\times {A}_{552})\) 40; \(\text{carotenoids}\,({\rm{\mu }}g/\text{mL})=[(1000\times {A}_{470})-\) \((2.86\times Ch{l}_{a}[\mu g/{mL}])]/221\) 41.
References
Oliver, J. W. & Atsumi, S. Metabolic design for cyanobacterial chemical synthesis. Photosynthesis research 120, 249â261 (2014).
Angermayr, S. A., Gorchs Rovira, A. & Hellingwerf, K. J. Metabolic engineering of cyanobacteria for the synthesis of commodity products. Trends in biotechnology 33, 352â361 (2015).
Berla, B. M. et al. Synthetic biology of cyanobacteria: unique challenges and opportunities. Frontiers in microbiology 4, 246 (2013).
Breitmaier, E. Terpenes: flavors, fragrances, pharmaca, pheromones. (John Wiley & Sons, 2006).
George, K. W., Alonso-Gutierrez, J., Keasling, J. D. & Lee, T. S. Isoprenoid drugs, biofuels, and chemicalsâartemisinin, farnesene, and beyond. Advances in biochemical engineering/biotechnology 148, 355â389 (2015).
Wang, G., Tang, W. & Bidigare, R. R. Terpenoids as therapeutic drugs and pharmaceutical agents. In Natural products pp. 197â227, Springer (2005).
Westfall, P. J. et al. Production of amorphadiene in yeast, and its conversion to dihydroartemisinic acid, precursor to the antimalarial agent artemisinin. Proceedings of the National Academy of Sciences of the United States of America 109, E111â118 (2012).
Martin, V. J., Pitera, D. J., Withers, S. T., Newman, J. D. & Keasling, J. D. Engineering a mevalonate pathway in Escherichia coli for production of terpenoids. Nat Biotechnol 21, 796â802 (2003).
Tracy, N. I., Chen, D., Crunkleton, D. W. & Price, G. L. Hydrogenated monoterpenes as diesel fuel additives. Fuel 88, 2238â2240 (2009).
Chuck, C. J. & Donnelly, J. The compatibility of potential bioderived fuels with Jet A-1 aviation kerosene. Applied Energy 118, 83â91 (2014).
Gao, Z., Zhao, H., Li, Z., Tan, X. & Lu, X. Photosynthetic production of ethanol from carbon dioxide in genetically engineered cyanobacteria. Energy & Environmental Science 5, 9857â9865 (2012).
Lan, E. I. & Liao, J. C. ATP drives direct photosynthetic production of 1-butanol in cyanobacteria. Proceedings of the National Academy of Sciences of the United States of America 109, 6018â6023 (2012).
Liu, X., Sheng, J. & Curtiss, R. 3rd Fatty acid production in genetically modified cyanobacteria. Proceedings of the National Academy of Sciences of the United States of America 108, 6899â6904 (2011).
Kiyota, H., Okuda, Y., Ito, M., Hirai, M. Y. & Ikeuchi, M. Engineering of cyanobacteria for the photosynthetic production of limonene from CO2. Journal of biotechnology 185, 1â7 (2014).
Halfmann, C., Gu, L. & Zhou, R. Engineering cyanobacteria for the production of a cyclic hydrocarbon fuel from CO2 and H2O. Green Chemistry 16, 3175â3185 (2014).
Ershov, Y. V., Gantt, R. R., Cunningham Jr, F. X. Jr. & Gantt, E. Isoprenoid biosynthesis in Synechocystis sp. strain PCC6803 is stimulated by compounds of the pentose phosphate cycle but not by pyruvate or deoxyxylulose-5-phosphate. Journal of bacteriology 184, 5045â5051 (2002).
Ranganathan, S., Suthers, P. F. & Maranas, C. D. OptForce: an optimization procedure for identifying all genetic manipulations leading to targeted overproductions. PLoS Comput Biol 6, e1000744 (2010).
Saha, R. et al. Reconstruction and comparison of the metabolic potential of cyanobacteria Cyanothece sp. ATCC 51142 and Synechocystis sp. PCC 6803. PloS one 7, e48285 (2012).
Young, J. D., Shastri, A. A., Stephanopoulos, G. & Morgan, J. A. Mapping photoautotrophic metabolism with isotopically nonstationary 13C flux analysis. Metabolic engineering 13, 656â665 (2011).
Bohlmann, J., Steele, C. L. & Croteau, R. Monoterpene synthases from grand fir (Abies grandis). cDNA isolation, characterization, and functional expression of myrcene synthase, (â)-(4S)-limonene synthase, and (-)-(1S,5S)-pinene synthase. The Journal of biological chemistry 272, 21784â21792 (1997).
Lucker, J. et al. Monoterpene biosynthesis in lemon (Citrus limon). cDNA isolation and functional analysis of four monoterpene synthases. European journal of biochemistry / FEBS 269, 3160â3171 (2002).
Colby, S. M., Alonso, W. R., Katahira, E. J., McGarvey, D. J. & Croteau, R. 4S-limonene synthase from the oil glands of spearmint (Mentha spicata). cDNA isolation, characterization, and bacterial expression of the catalytically active monoterpene cyclase. The Journal of biological chemistry 268, 23016â23024 (1993).
Davies, F. K., Work, V. H., Beliaev, A. S. & Posewitz, M. C. Engineering Limonene and Bisabolene Production in Wild Type and a Glycogen-Deficient Mutant of Synechococcus sp. PCC 7002. Frontiers in bioengineering and biotechnology 2, 21 (2014).
Alonso-Gutierrez, J. et al. Metabolic engineering of Escherichia coli for limonene and perillyl alcohol production. Metabolic engineering 19, 33â41 (2013).
Williams, D. C., McGarvey, D. J., Katahira, E. J. & Croteau, R. Truncation of limonene synthase preprotein provides a fully active âpseudomatureâ form of this monoterpene cyclase and reveals the function of the amino-terminal arginine pair. Biochemistry 37, 12213â12220 (1998).
Ng, A. H., Berla, B. M. & Pakrasi, H. B. Fine-Tuning of Photoautotrophic Protein Production by Combining Promoters and Neutral Sites in the Cyanobacterium Synechocystis sp. Strain PCC 6803. Applied and environmental microbiology 81, 6857â6863 (2015).
Poliquin, K. et al. Inactivation of sll1556 in Synechocystis strain PCC 6803 impairs isoprenoid biosynthesis from pentose phosphate cycle substrates in vitro. J Bacteriol 186, 4685-4693 (2004).
Tai, M. & Stephanopoulos, G. Engineering the push and pull of lipid biosynthesis in oleaginous yeast Yarrowia lipolytica for biofuel production. Metab Eng 15, 1â9 (2013).
Hosfield, D. J. et al. Structural basis for bisphosphonate-mediated inhibition of isoprenoid biosynthesis. The Journal of biological chemistry 279, 8526â8529 (2004).
Burke, C. & Croteau, R. Geranyl diphosphate synthase from Abies grandis: cDNA isolation, functional expression, and characterization. Archives of biochemistry and biophysics 405, 130â136 (2002).
Berla, B. M. & Pakrasi, H. B. Upregulation of plasmid genes during stationary phase in Synechocystis sp. strain PCC 6803, a cyanobacterium. Applied and environmental microbiology 78, 5448â5451 (2012).
Yu, J. et al. Synechococcus elongatus UTEX 2973, a fast growing cyanobacterial chassis for biosynthesis using light and CO(2). Scientific reports 5, 8132 (2015).
Wang, X. et al. Enhanced limonene production in cyanobacteria reveals photosynthesis limitations. Proceedings of the National Academy of Sciences of the United States of America 113, 14225â14230 (2016).
Estevez, J. M., Cantero, A., Reindl, A., Reichler, S. & Leon, P. 1-Deoxy-D-xylulose-5-phosphate synthase, a limiting enzyme for plastidic isoprenoid biosynthesis in plants. The Journal of biological chemistry 276, 22901â22909 (2001).
Wang, C. W., Oh, M. K. & Liao, J. C. Engineered isoprenoid pathway enhances astaxanthin production in Escherichia coli. Biotechnology and bioengineering 62, 235â241 (1999).
Banerjee, A. & Sharkey, T. D. Methylerythritol 4-phosphate (MEP) pathway metabolic regulation. Natural product reports 31, 1043â1055 (2014).
Poliquin, K., Cunningham, F. X., Gantt, R. R. & Gantt, E. In Isoprenoid Synthesis in Plants and Microorganisms: New Concepts and Experimental Approaches (eds Thomas J. Bach & Michel Rohmer) 51â63 (Springer New York, 2013).
Taton, A. et al. Broad-host-range vector system for synthetic biology and biotechnology in cyanobacteria. Nucleic Acids Res 42, e136 (2014).
Gibson, D. G. et al. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nature methods 6, 343â345 (2009).
Porra, R., Thompson, W. & Kriedemann, P. Determination of accurate extinction coefficients and simultaneous equations for assaying chlorophylls a and b extracted with four different solvents: verification of the concentration of chlorophyll standards by atomic absorption spectroscopy. Biochimica et Biophysica Acta (BBA)-Bioenergetics 975, 384â394 (1989).
Wellburn, A. R. The spectral determination of chlorophylls a and b, as well as total carotenoids, using various solvents with spectrophotometers of different resolution. Journal of plant physiology 144, 307â313 (1994).
Acknowledgements
Funding to support this work was provided by the Office of Science (BER), U. S. Department of Energy, to HBP. PCL was supported by a fellowship from the McDonnell International Scholars Academy at Washington University. We thank all members of the Pakrasi research group for collegial discussions.
Author information
Authors and Affiliations
Contributions
P.C.L. performed experiments; R.S. performed computational studies; All authors designed research, wrote and edited the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lin, PC., Saha, R., Zhang, F. et al. Metabolic engineering of the pentose phosphate pathway for enhanced limonene production in the cyanobacterium Synechocysti s sp. PCC 6803. Sci Rep 7, 17503 (2017). https://doi.org/10.1038/s41598-017-17831-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-17831-y
This article is cited by
-
Nutraceutical prospects of genetically engineered cyanobacteria- technological updates and significance
World Journal of Microbiology and Biotechnology (2024)
-
The current situations and limitations of genetic engineering in cyanobacteria: a mini review
Molecular Biology Reports (2023)
-
Exploring the potential of the model cyanobacterium Synechocystis PCC 6803 for the photosynthetic production of various high-value terpenes
Biotechnology for Biofuels and Bioproducts (2022)
-
Order-of-magnitude enhancement in photocurrent generation of Synechocystis sp. PCC 6803 by outer membrane deprivation
Nature Communications (2022)
-
Machine learning-informed and synthetic biology-enabled semi-continuous algal cultivation to unleash renewable fuel productivity
Nature Communications (2022)