Skip to main content

Engineering yeast mitochondrial metabolism for 3-hydroxypropionate production

Abstract

Background

With unique physiochemical environments in subcellular organelles, there has been growing interest in harnessing yeast organelles for bioproduct synthesis. Among these organelles, the yeast mitochondrion has been found to be an attractive compartment for production of terpenoids and branched-chain alcohols, which could be credited to the abundant supply of acetyl-CoA, ATP and cofactors. In this study we explored the mitochondrial potential for production of 3-hydroxypropionate (3-HP) and performed the cofactor engineering and flux control at the acetyl-CoA node to maximize 3-HP synthesis.

Results

Metabolic modeling suggested that the mitochondrion serves as a more suitable compartment for 3-HP synthesis via the malonyl-CoA pathway than the cytosol, due to the opportunity to obtain a higher maximum yield and a lower oxygen consumption. With the malonyl-CoA reductase (MCR) targeted into the mitochondria, the 3-HP production increased to 0.27 g/L compared with 0.09 g/L with MCR expressed in the cytosol. With enhanced expression of dissected MCR enzymes, the titer reached to 4.42 g/L, comparable to the highest titer achieved in the cytosol so far. Then, the mitochondrial NADPH supply was optimized by overexpressing POS5 and IDP1, which resulted in an increase in the 3-HP titer to 5.11 g/L. Furthermore, with induced expression of an ACC1 mutant in the mitochondria, the final 3-HP production reached 6.16 g/L in shake flask fermentations. The constructed strain was then evaluated in fed-batch fermentations, and produced 71.09 g/L 3-HP with a productivity of 0.71 g/L/h and a yield on glucose of 0.23 g/g.

Conclusions

In this study, the yeast mitochondrion is reported as an attractive compartment for 3-HP production. The final 3-HP titer of 71.09 g/L with a productivity of 0.71 g/L/h was achieved in fed-batch fermentations, representing the highest titer reported for Saccharomyces cerevisiae so far, that demonstrated the potential of recruiting the yeast mitochondria for further development of cell factories.

Background

The yeast Saccharomyces cerevisiae has been extensively engineered for synthesis of biofuels and other non-ethanol bioproducts [1, 2]. Most synthetic pathways are targeted to the cytosol due to its big volume, high coverage of metabolic enzymes and metabolites [3, 4]. Recent studies have revealed that the yeast mitochondrion could be attractive to produce terpenoids, fumarate, ornithine, and branched-chain alcohols [5,6,7]. With the mevalonate pathway and sequential synthetic genes targeted into the mitochondria, researchers have greatly improved the titers of isoprene, patchoulol, valencene, amorphadiene, linalool and santalene [8,9,10,11,12]. Moreover, mitochondrial compartmentalization of the 2-ketoacid elongation pathway has obtained 1.24 g/L isopentanol, the highest titer reported in yeast [13]. Therefore, there has been growing interest in exploring the potential of mitochondrial compartmentalization of synthetic pathways.

The improved production of subcellular compartmentation may be ascribed to physical separation, which concentrates overexpressed enzymes and substrates in a smaller volume to reduce the amount of the byproducts or the toxicity of some intermediate metabolites. Specially for mitochondrial compartmentation, there are more acetyl-CoA, ATP and NADH available compared with the cytosol. Furthermore, the mitochondrial matrix has a higher pH and redox potential than the cytosol, and this environment could be more suitable for bacterial enzymes [14]. Taking together, the mitochondrial compartment could be a promising alternative for production of 3-hydroxypropionate (3-HP), which could be synthesized from acetyl-CoA via the malonyl-CoA pathway catalyzed by acetyl-CoA carboxylase (ACC) and a bacterial enzyme malonyl-CoA reductase (MCR). 3-HP serves as a promising platform chemical to produce various C3-based chemicals, including acrylic acid and its derivatives, which are wildly used in the fields of plastics, pharmaceuticals and agriculture. Extensive studies of 3-HP bioproduction have been performed [15,16,17].

Previous studies on the malonyl-CoA pathway have been performed in yeast cytosol and revealed several crucial factors for 3-HP synthesis, including efficient supply of acetyl-CoA and malonyl-CoA, balanced supply of redox cofactor NADPH, and improved activities of MCR enzymes [18,19,20]. The native cytosolic acetyl-CoA synthesis is highly ATP demanding, and engineering efforts on acetyl-CoA biosynthetic pathways have been performed and enhanced 3-HP production by up to fivefold [21, 22]. Cytosolic malonyl-CoA is mainly involved in fatty acid synthesis, which is regulated stringently at multiple levels. The acetyl-CoA carboxylase Acc1 functions as one of the important nodes in fatty acid metabolism and is stringently regulated transcriptionally and post-translationally [23, 24]. With malonyl-CoA biosensor established in yeast, it has been found that overexpression, dynamic control and mutagenesis of dephosphorylation sites of Acc1 can direct more malonyl-CoA to 3-HP synthesis rather than fatty acid synthesis [25,26,27]. Combined with these strategies for enhanced precursor supply, engineering redox metabolism for optimized NADPH supply further improves 3-HP production [20,21,22]. Moreover, it has been reported that MCR engineering with domain dissection and balanced expression of the two dissected parts can significantly improve 3-HP production from 0.107 to 3.72 g/L in E. coli [28, 29]. With the dissected enzymes expressed in yeast cytosol, a significant improvement of 3-HP production is also achieved from 0.6 to 2.4 g/L, suggesting MCR as another limiting factor for 3-HP [18].

Compared with cytosolic precursor supply, the mitochondrial acetyl-CoA derived from pyruvate dehydrogenase (PDH) is used as fuel of the TCA cycle and is reported to be 10 times more concentrated than its cytosolic counterpart, which renders the mitochondria great potential for 3-HP synthesis. Mitochondrial acetyl-CoA can be converted to malonyl-CoA by Hfa1 [30, 31]. This mitochondrial isoform Hfa1 shares 72% sequence similarity with Acc1, and is essential for synthesis of mitochondrial fatty acids, especially lipoic acid cofactor [30]. Its initiation regulation is reported to be translated from non-canonical initial site with a mitochondrial targeting sequence at its N terminal [32]. So far, it remains unclear whether Hfa1 contains any regulation sites, and it would be interesting to investigate its effects on the synthesis of malonyl-CoA derived bioproduct.

While the cytosolic NADPH is mainly derived from the reactions catalyzed by glucose-6-phosphate dehydrogenase Zwf1, acetaldehyde dehydrogenase Ald6 and isocitrate dehydrogenase Idp2, the mitochondrial NADPH is relying on the reactions catalyzed by NAD+/NADH kinase Pos5 and NADP+ dependent enzymes, including acetaldehyde dehydrogenases Ald4 and Ald5, isocitrate dehydrogenase Idp1, and malic enzyme Mae1 [33, 34]. Pos5 is considered as the main source of mitochondrial NADPH and possesses higher NADH kinase activity than NAD kinase activity [35, 36], suggesting active NADPH synthesis from NADH with ATP consumed, both of which are abundant in the mitochondriona. Taking together, the yeast mitochondrion seems to be a suitable subcellular compartment for 3-HP production.

Here we explored yeast mitochondrion for its potential in 3-HP production via the malonyl-CoA pathway. Engineering efforts were performed to enhance both malonyl-CoA and NADPH supply (Fig. 1). The work presented here gains a better understanding for harnessing the yeast mitochondrion for synthesis of acetyl-CoA derived biochemicals, especially those with NADPH as a cofactor.

Fig. 1
figure 1

Schematic representation of rewiring mitochondrial metabolism for 3-HP production. The native enzymes were shown in blue, whereas heterologous enzymes in orange. Acc1 acetyl-CoA carboxylase, MCR malonyl-CoA reductase, PDH pyruvate dehydrogenase complex, Hfa1 mitochondrial acetyl-CoA carboxylase, Acc1** acetyl-CoA carboxylase mutant, Cit1 mitochondrial citrate synthase, Aco1/2 aconitase, Idh1/2 mitochondrial NAD+-dependent isocitrate dehydrogenase, Idp1 mitochondrial NADP-specific isocitrate dehydrogenase, Mdh1 mitochondrial malate dehydrogenase, Mae1 mitochondrial malic enzyme, EcMaeA* malic enzyme mutant from E. coli, Pos5 mitochondrial NADH kinase

Materials and methods

Construction of plasmids and strains

The E. coli DH5α strain was used for plasmid construction. The yeast strain CEN.PK 113-5D was used as the background strain. Genetic manipulations, including gene integration, knockout and promoter replacement, were performed using the GTR-CRISPR system [37]. Plasmids used in the study were constructed with Golden Gate methods if not specified [38]. The strains and plasmids used and constructed were listed in Additional file 1.

The intact MCR from Chloroflexus aurantiacus was cloned in pMCR1, while its dissected parts of MCR-N and MCR-C were cloned in pMCR2 in a previous study for mitochondrial localization [39]. The plasmid pMCR-mC was constructed to introduce the mutations N940V, K1106W, and S1114R into MCR-C by assembling 3 fragments amplified with primers MP1–4 from pMCR2 [39] into the vector pUGG1. The plasmid pMCR3 was constructed by assembling the fragments of MCR-N amplified with primers MP5&6 and MCR-mC amplified with primers MP7&8 into pUGG1. The plasmid pPOS5 was constructed to overexpress POS5 under the control of TDH3 promoter. The TDH3p fragment amplified with primers MP9&10 and the fragment containing POS5 and its terminator amplified with primers MP11&12 using yeast genomic DNA as the template were assembled into pUGG1. The plasmid pALD4 was constructed to overexpress ALD4 under the control of PGK1 promoter. The PGK1p fragment amplified with primers MP13&14 and the fragment containing ALD4 and its terminator amplified with primers MP15&16 were assembled into pUGG1. The plasmid pIDP1 was constructed to overexpress IDP1 under the control of TEF1 promoter. The TEF1p fragment amplified with primers MP17&18 and the fragment containing IDP1 its terminator amplified with primers MP19&20 were assembled into pUGG1. The plasmid pMaeA* was constructed to overexpress MaeA* from E. coli under the TEF1 promoter. The fragments containing TEF1p and the mitochondrial tag CoxIVm were amplified from pSSB02 [40] with primers MP21–24. The mutated MaeA fragments were amplified with primers MP25–28 from the E. coli genomic DNA. The ADH1t fragment was amplified with primers MP29&30 from yeast genomic DNA. Then, these fragments were assembled into pUGG1.

Mitochondrial overexpression of ACC1 was achieved using different promoters of varied strength and types, including TEF1p, CYC1p, CUP1p and GAL1p. The fragments containing ACC1S659A, S1157A were amplified from pCfB376 [27] with primers MP31–35 and MP29. The fragments containing TEF1p and CoxIVm were amplified from pSSB02 with primers MP21–24. Then, the two fragments were assembled into pUGG1 to generate pACC1**_TEF1p. Then, ACC1S659A, S1157A was mutated to the wild-type ACC1 by assembling two fragments amplified with primers MP36&37 and MP38&39, generating pACC1_TEF1p. Based on pACC1_TEF1p, TEF1p was replaced to CYC1p by assembling the fragment amplified with primers MP40–43 using the NEBuilder assembly to generate pACC1_CYC1p. Based on pACC1**_TEF1p, TEF1p and 2μ ori were replaced to CUP1pI and CEN4/6, respectively, by assembling the fragments amplified with primers MP44–51 using the NEBuilder assembly to generate pACC1**_CUP1p. Based on pACC1_TEF1p and pACC1**_TEF1p, TEF1p was replaced to GAL1p using NEBuilder assembly to generate pACC1_GAL1p and pACC1**_GAL1p. Based on pACC1**_CUP1p, CUP1p was replaced to GAL1p using the NEBuilder assembly to generate pACC1**_GAL1p_CEN.

For genome integration, reported integrative sites were used, including XI-3, X-4 and Int14 [41,42,43]. Intact MCR was integrated at XI-3 using pCas9_XI-3.3 and donor DNA amplified from pMCR1 with primers MP54&55. Dissected MCR was integrated at XI-3 site using pCas9_XI-3.3 and donor DNA fragments amplified from pMCR2 or pMCR3 with primers MP55&56. MCR-N and MCR-mC were integrated at X-4 and Int14 with donor DNA fragments amplified with MP57&58, and MP60&61, respectively. Moreover, several native promoters were replaced with strong constitutive promoters to enhance their expression. POS5p was replaced to TDH3p using pCas9_POS5p and donor DNA amplified with MP63&64. IDP1p was replaced to TEF1p using pCas9_IDP1p and donor DNA amplified with MP66&67. MAE1p was replaced to TDH3p using pCas9_MAE1p and donor DNA amplified with MP69&70. HFA1p was replaced to TEF1p using pCas9_HFA1p and donor DNA amplified with MP72&73. Meanwhile, IDH1 was deleted using pCas9_idh1 and donor DNA amplified with primers MP75&76. IDH2 was deleted using pCas9_idh2 and donor DNA amplified with primers MP78&79. CIT1 was deleted using pCas9_cit1 and donor DNA amplified with primers MP81&82. GAL80 was deleted using pCas9_gal80 and donor DNA amplified with primers MP85&86.

Media and cultivation condition

E. coli strains were cultivated in LB medium supplemented with 80 mg/L ampicillin when needed. Yeast strains were cultivated in YPD medium, or SC-URA medium to maintain URA3-marked plasmids, as previously described in [44]. Fermentation medium for 3-HP production in shake flasks was composed of 7.5 g/L (NH4)2SO4, 14.4 g/L KH2PO4, 0.5 g/L MgSO4·7H2O, 20 g/L glucose, trace metal solution and vitamin solution as described in [45] and 40 mg/L uracil was added when needed. For induction for the CUP1 promoter, 300 μM CuSO4 was supplemented at 24 h during the fermentation. For all shake flask experiments, 20 mL cultures were cultivated in 100 mL shake flasks with the initial OD600 of 0.1 at 30 °C and 200 rpm.

For fed-batch cultivations in 5 L bioreactors, the medium was composed of 5 g/L (NH4)2SO4, 3 g/L KH2PO4, 0.5 g/L MgSO4·7H2O, 20 g/L glucose, as well as the trace metal solution and vitamin solution as described in [45]. The seed was pre-cultured using the same medium in the 500 mL shake flasks, and then inoculated into the bioreactors with 2 L medium with the initial OD600 of 0.4. During the fermentation, the dissolved oxygen (DO) concentration was controlled above 15% by adjusting the stirring speed within 200–700 rpm and aeration with 0.5–1 vvm. pH was controlled at 5.0 by feeding NH3·H2O, and the temperature was controlled at 30 °C. The feed medium was added when glucose was almost exhausted, composing 5.6 g/L K2SO4, 15 g/L KH2PO4, 2.8 g/L MgSO4·7H2O, 700 g/L glucose, trace metal solution (7 ×) and vitamins solution (7 ×), and the feed rates were continuously adjusted to remain low glucose residue and ethanol residue in the culture.

Measurement of biomass and extracellular metabolites

The optical density was measured at 600 nm OD600 by GENESYS 30 Visible Spectrophotometer (Thermo Electron Scientific, Madison, USA).

Extracellular metabolites analysis was performed using the HPLC (Shimadzu LC-20AT, Japan) equipped with RID detector and PDA detector using Aminex HPX-87H column (Bio-Rad). The culture samples were centrifuged, and the supernatant was filtered by 0.22 μm membrane and stored in − 20 °C until analysis. The eluent of 0.5 mM H2SO4 was used at a flow rate of 0.5 mL/min with the column temperature of 65 °C for the measurement of glucose, ethanol and 3-HP. Glucose, ethanol and 3-HP were quantified using LabSolution Software (Shimadzu).

Flux balance analysis

The maximum theoretical 3-HP yields with the malonyl-CoA pathway expressed in the cytosol and mitochondria were calculated using flux balance analysis (FBA) [46] with the yeast genome-scale metabolic model Yeast8 [47], which was correspondingly modified by adding the malonyl-CoA pathway in the cytosol and mitochondria. The glucose uptake rate was fixed at 1 mmol/g DCW/h and the 3-HP secretion rate was maximized, and the simulated flux distribution can be then used to determine the active reactions, i.e., those whose flux is non-zero. The simulations were performed in MATLAB using the Cobra toolbox v.3.0 [48], and the simulated flux distributions were then used to calculate the maximum theoretical yields.

Results and discussion

Mitochondrial localization of MCR enhanced 3-HP synthesis

Flux balance analysis was performed based on the genome scale metabolic model of S. cerevisiae Yeast8 [47], to calculate the maximum theoretical yields of 3-HP for cytosolic and mitochondrial synthesis. A theoretical 3-HP yield of 1.71 mol/mol glucose was obtained with the oxygen requirement of 0.88 mol/mol glucose for mitochondrial synthesis, compared with those of 1.62 mol 3-HP/mol glucose and 1.15 mol O2/mol glucose for cytosolic synthesis (Additional file 2). The higher theoretical yield and lower oxygen requirement could be credited to the rich acetyl-CoA and ATP source in the mitochondria, as cytosolic acetyl-CoA synthesis was highly ATP demanding, rendering mitochondrial synthesis great potentials for future engineering in both basic and applied research.

In our previous research, we noticed that when the mitochondrial MCR was expressed in an episomal plasmid, the 3-HP titer reached 0.27 g/L [39], compared to the titer of 0.09 g/L when the cytosolic MCR was expressed [21], suggesting the yeast mitochondrion as an attractive subcellular compartment for 3-HP synthesis. With mitochondrial CaMCR integrated into the XI-3 site of the chromosome, the strain mMCR1 produced 0.11 g/L 3-HP, suggesting that CaMCR activity was crucial for 3-HP production. Therefore, dissected forms of CaMCR with higher enzymatic activities were then expressed in the mitochondria, including MCR-N&MCR-C, and MCR-N&MCR-mC with amino acid mutations. With the dissected genes integrated at XI-3 site, the constructed strains NC1 and NmC1 produced 3-HP at the titers of 0.43 g/L and 1.29 g/L, respectively (Fig. 2A). After a 2nd round integration of MCR-N & MCR-mC at the X-4 site in NmC1, the titer of the constructed strain NmC2 increased to 3.64 g/L, with its 70.6% accumulated during the ethanol phase, when the repressed mitochondrial activities were relieved. These results suggested that glucose derepression under glucose depletion conditions [49], contributed greatly to 3-HP synthesis in the presence of increased MCR activities. It also revealed that the MCR activity was a limiting factor for the mitochondrial 3-HP synthesis.

Fig. 2
figure 2

3-HP production in yeast mitochondria with the optimized expression of MCR genes. A Intact and dissected MCR gene expressed in the cytosol and mitochondria. cMCR_p, the strain harboring pYC01 with intact MCR targeted into the cytosol; mMCR_p, the strain harboring pMCR1 with intact MCR targeted into the mitochondria; mMCR1, the strain with intact MCR integrated at the site XI-3; NC1, the strain with MCR-N and MCR-C integrated at the site XI-3; NmC1, the strain with MCR-N and MCR-mC integrated at the site XI-3; NmC2, the strain with MCR-N and MCR-mC integrated at the sites XI-3 and X-4. B Varied expression of dissected MCR-N and MCR-mC in the mitochondria. NmC2 pUGG1, the strain NmC2 harboring an empty plasmid pUGG1; NmC2 pMCR_N, the strain NmC2 harboring pMCR-N; NmC2 pMCR_mC, the strain NmC2 harboring pMCR-mC; NmC3, the strain with MCR-N and MCR-mC integrated at the sites XI-3, X-4 and Int14. The cultivations were performed in biologically triplicate and error bars represent ± standard errors. Statistical analysis was performed using one-tailed Student’s t test (**p < 0.01)

A previous study in E. coli showed that functional balance between MCR-N and MCR-mC could greatly enhance 3-HP synthesis by expressing three copies of MCR-N with chromosomal integration and MCR-mC in a high copy number plasmid. We then investigated whether altering the expression of MCR-N and MCR-mC could benefit 3-HP production in yeast mitochondria. However, neither of the constructed strains exhibited significant improvements in cell growth or 3-HP production (Fig. 2B), suggesting that the current expression of MCR-N and MCR-mC might be well-balanced (Additional file 3: Fig. S1). With a third integration of MCR-N and MCR-mC at the Int14 site of NmC2, the constructed strain NmC3 reached a titer of 4.42 g/L at 60 h, which was comparable with the highest titer achieved with the production in the yeast cytosol with a substantially engineered yeast chassis [18].

Mitochondrial redox factor engineering for sufficient NADPH regeneration

Mitochondrial NADPH is generated from NADH by Pos5, whereas re-generation of NADPH from NADP+ is carried out by NADP+-dependent dehydrogenases including Ald4, Idp1 and Mae1 [33, 50, 51]. These enzymes were overexpressed under strong promoters to increase the availability of mitochondrial NADPH for 3-HP biosynthesis.

Indeed, when the plasmid pPOS5 for POS5 overexpression was transformed into the strain NmC2, the 3-HP titer increased by 23.4%, while the biomass decreased by 13.6% (Fig. 3A). Yet, when ALD4 was overexpressed in a high copy plasmid pALD4 into NmC2, no significant improvements were observed in 3-HP synthesis and cell growth. Similarly, when MAE1 was overexpressed in NmC2 with its native promoter switched to a strong promoter TDH3p, 3-HP titer decreased from 3.42 to 2.56 g/L, which might be resulted from the decrease in cell growth. Nevertheless, with the native promoter replaced with a strong promoter TEF1p, IDP1 overexpression in the strain NmC2 resulted in an increase in 3-HP production to 3.90 g/L, as well as a moderate improvement in cell growth. It has been reported that these dehydrogenases were transcriptionally regulated, and the varied effects of their overexpression on cell growth and 3-HP synthesis suggested that there might be other regulatory mechanisms involved in their activities.

Fig. 3
figure 3

3-HP production with increased NADPH synthesis. A Overexpression of NADPH generating enzymes in the strain NmC2. NmC2 pALD4, the strain NmC2 harboring pALD4; NmC2 MAE1, the strain NmC2 with MAE1 overexpressed; NmC2 pPOS5, the strain NmC2 harboring pPOS5; NmC2 IDP1, the strain NmC2 with IDP1 overexpressed; N2IP, the strain NmC2 with IDP1 and POS5 overexpressed. B Overexpression of NADPH generating enzymes in the strain NmC3. N3IP, the strain NmC3 with IDP1 and POS5 overexpressed; N3IP pMaeA*, the strain N3IP harboring pMaeA*. The cultivations were performed in biologically triplicate and error bars represent ± standard errors. Statistical analysis was performed using one-tailed Student’s t test (*p < 0.05, **p < 0.01)

With co-overexpression of POS5 and IDP1 in NmC2, the 3-HP titer of the constructed strain N2IP increased to 4.55 g/L, which was 33% higher than that of the strain NmC2. These results revealed that two routes for NADPH regeneration from NADH or NADP+ could synergistically boost 3-HP titer, with the one from NADH catalyzed by Pos5 as the main contributor. Then, both IDP1 and POS5 were overexpressed under strong promoters in the strain NmC3, generating the strain N3IP, with 3-HP titer increasing from 4.42 to 5.11 g/L (Fig. 3B).

However, overexpression of the native MAE1 significantly slowed down cell growth and reduced 3-HP titer (Fig. 3A). To test whether the adverse result of MAE1 overexpression for 3-HP production is caused by its cofactor preference of NAD+ over NADP+, a mutated malic enzyme MaeA* that prefers NADP+ as the cofactor [52], was expressed in the mitochondria of the strain N3IP pMaeA*. Yet, the cell biomass and 3-HP titer still reduced by 17% and 20% (Fig. 3B), respectively. The negative effect on cell growth of malic enzyme overexpression might be that it drained TCA cycle intermediates and slowed down the operation of the TCA cycle, which provided essential energy for cell growth and maintenance.

Limited flux of TCA cycle resulted in acetate overflow and impaired growth

Mitochondrial acetyl-CoA is dedicated to fuel the TCA cycle to provide energy for cell growth and maintenance. To redirect acetyl-CoA to 3-HP synthesis, attempts were made to limit the TCA cycle flux. The first reaction of the TCA cycle catalyzed by citrate synthase Cit1 is considered as the rate limiting step. Therefore, we first downregulated CIT1 transcription by transforming a plasmid harboring dCas9 protein and two gRNAs targeting the CIT1 promoter. However, CIT1 downregulation resulted in a decrease in 3-HP titer from 3.75 to 3.11 g/L (Fig. 4A). When CIT1 was deleted in the strain N3IP, the mutant showed significantly reduced growth and 3-HP production. The final biomass and 3-HP production were 56.0% and 24.4% of those in the strain N3IP, respectively. Meanwhile, acetate was accumulated during the cultivation of the mutant, which reached to 4.8 g/L (Additional file 3: Fig. S2). These results suggested that the disrupted TCA cycle resulted in acetate overflow, probably due to insufficient ATP production. The flux between 3-HP synthesis and the TCA cycle is, therefore, needed to be well-balanced.

Fig. 4
figure 4

Limited flux of TCA cycle impaired cell growth and 3-HP production. A CIT1 downregulation and deletion in strains NmC2 and N3IP. NmC2 pdCas9_CIT, the strain NmC2 harboring pdCas9_CIT; N3IP cit1, the strain N3IP with cit1 deletion. B IDH1/2 downregulation and deletion in NmC2 strains with and without IDP1 and POS5 overexpression. NmC2 pdCas_IDH, the strain NmC2 harboring pdCas_IDH; NmC2 IDP1 idh1, the strain NmC2 IDP1 with idh1 deletion; NmC2 IDP1 idh2, the strain NmC2 IDP1 with idh2 deletion; N2IP idh1, the strain N2IP with idh1 deletion; N2IP idh2, the strain N2IP with idh2 deletion. The cultivations were performed in biologically triplicate and error bars represent ± standard errors

Then another key enzyme of the TCA cycle, isocitrate dehydrogenase, that catalyzes the oxidation of isocitrate to alpha-ketoglutarate, was engineered to limit the TCA flux. Mitochondria contain both NAD+-dependent and NADPH+-dependent isozymes encoded by IDH1/2 and IDP1, respectively. These two isozymes are regulated independently as they played different roles [53, 54], and both of them could tune the mitochondrial TCA activity [55, 56]. When the strain NmC2 was transformed with a plasmid harboring the dCas9 protein and gRNAs targeting IDH1 and IDH2 promoters, the mutant exhibited impaired cell growth with 3-HP production decreased by 20.7% and 20.1%, respectively (Fig. 4B). Similarly, in the strain NmC2 with IDP1 overexpressed IDH1 and IDH2 deletion also decreased 3-HP production by 21.1%, and 16.2%, respectively (Fig. 4B). Meanwhile, a bigger biomass decrease observed in a CIT1 deletion strain may have resulted from acetate overflow, which accumulated to 1.65 g/L in the culture. In the strain NmC2 with IDP1 and POS5 overexpressed, IDH1 and IDH2 deletion resulted in decreases in 3-HP titer by 21.1%, and 17.9%, respectively, yet no decreases in cell growth (Fig. 4B). Although IDH1/2 downregulation and deletion resulted in impaired cell growth in the two strains, their deletion did not cause significant growth reduction in the strain with IDP1 and POS5 overexpression, probably because NADPH supply might still be the key factor for cell growth and 3-HP production in these strains.

Fine-tuned acetyl-CoA carboxylase expression provided sufficient malonyl-CoA

Besides sufficient acetyl-CoA supply, sufficient malonyl-CoA supply is required for production of malonyl-CoA derived product, including cytosolic fatty acids and 3-HP expression [22, 27], as ACC1 was stringently regulated in different levels. The regulation of the mitochondrial acetyl-CoA carboxylase Hfa1 has not been investigated in detail, and its impact on mitochondrial 3-HP production needs to be investigated. As shown in Fig. 5A, HFA1 overexpression in the strains mMCR1, NC1 and NmC1 boosted 3-HP production by 2.65-fold, 1.63-fold, and 1.79-fold, respectively. These results revealed that HFA1 overexpression greatly improved 3-HP synthesis during the ethanol growth phase (after 24 h). Yet, when HFA1 was overexpressed in the strain NmC2, it could not significantly improve the 3-HP titer. The results indicated that HFA1 might also be regulated to maintain the homeostasis of mitochondrial malonyl-CoA.

Fig. 5
figure 5

Overexpression of acetyl-CoA carboxylase in the mitochondria. A HFA1 overexpression in the 3-HP producing strains. mMCR1 HFA1, the strain mMCR1 with HFA1 overexpressed; NC1 HFA1, the strain NC1 with HFA1 overexpressed; NmC1 HFA1, the strain NmC1 with HFA1 overexpressed; NmC2 HFA1, the strain NmC2 with HFA1 overexpressed. B ACC1 and ACC1** overexpression in the stain NmC2. NmC2 pACC1_TEF1p, the strain NmC2 harboring pACC1_TEF1p; NmC2 pACC1_CYC1p, the strain NmC2 harboring pACC1_CYC1p; NmC2 pACC1**_TEF1p, the strain NmC2 harboring pACC1**_TEF1p. C Overexpression of ACC1** under the CUP1 promoter with and without Cu2+ induction. N3IP Cu2 + , the strain N3IP with 300 μM Cu2+ supplemented at 24 h; N3IP pACC1**_CUP1p, the strain N3IP harboring pACC1**_CUP1p; N3IP pACC1**_CUP1p Cu2 + , the strain N3IP pACC1**_CUP1p with 300 μM Cu2+ supplemented at 24 h. D Overexpression of ACC1 and ACC1** under the GAL1 promoter with galactose induction. NmC2 pUGG1 Gal, the strain NmC2 pUGG1 with 20 g/L galactose supplemented at 24 h; NmC2 pACC1_GAL1p Gal, the strain NmC2 harboring pACC1_GAL1p with 20 g/L galactose supplemented at 24 h; NmC2 pACC1**_GAL1p Gal, the strain NmC2 harboring pACC1**_GAL1p with 20 g/L galactose supplemented at 24 h; N3IP pUGG1 Gal, the strain N3IP harboring pUGG1 with 20 g/L galactose supplemented at 24 h; N3IP pACC1**_GAL1p Gal, the strain N3IP harboring pACC1**_GAL1p with 20 g/L galactose supplemented at 24 h. E Overexpression of ACC1** in the strain N3IP_1. N3IP_1, the strain N3IP with gal80 deletion; N3IP_2, the strain N3IP_1 harboring pACC1**_GAL1p_CEN. The cultivations were performed in biologically triplicate and error bars represent ± standard errors. Statistical analysis was performed using one-tailed Student’s t test (*p < 0.05, **p < 0.01)

We, therefore, overexpressed ACC1 in the mitochondria under the control of strong and weak promoters using TEF1p and CYC1p, respectively. Yet, impaired 3-HP synthesis was observed in both strains, with the overexpression under the strong promoter gave even lower titers compared with expression under the weak promoter (Fig. 5B), suggesting that ACC1 overexpression might impose great burden, or it may also be regulated post-translationally in the mitochondria. Therefore, we overexpressed an Acc1 mutant ACC1** abolishing posttranslational regulations, yet the constructed strain exhibited similar production as the strain with ACC1 overexpression but with impaired cell growth (Fig. 5B). Meanwhile, it was also observed that the strain with ACC1** overexpressed under TEF1p grew much slower, and the final 3-HP titer was achieved at 96 h compared with 60 h for other strains (Additional file 3: Fig. S3).

To release burden during early cell growth phase, ACC1** was then overexpressed in the mitochondria under the control of induced promoters CUP1p and Gal1p, respectively. In the strain N3IP, induced ACC1** expression with 300 μM Cu2+ supplemented at 24 h improved 3-HP titer from 5.01 to 5.50 g/L, while decreased 3-HP titer to 4.31 g/L without Cu2+ supplemented (Fig. 5C). The reduced biomass and 3-HP synthesis without Cu2+ supplemented at 24 h might be resulted from a basal expression as the medium contains 6 μM Cu2+ as a trace element. With the GAL1 promoter and 20 g/L galactose supplemented at 24 h, ACC1** overexpression in the strains NmC2 and N3IP improved 3-HP titer by 28.7% and 10.3%, respectively, while ACC1 overexpression failed to improve 3-HP synthesis (Fig. 5D). Then, GAL80 was deleted in N3IP strain to release glucose regulation on the GAL1 promoter [57]. In the constructed strain N3IP_2, the 3-HP titer increased by 18.6% to the final titer of 6.16 g/L (Fig. 5E), suggesting that fine-tuned ACC1* overexpression could coordinate malonyl-CoA requirements for cell growth and 3-HP synthesis, which thereof benefited 3-HP production in the mitochondria.

To further evaluate 3-HP production of the constructed strain, fed-batch cultivations were performed in 5 L bioreactors. To sustain the derepressed mitochondrial activity, continuous feeding was controlled under glucose limited conditions after 11 h, when glucose was almost exhausted (Fig. 6). The final 3-HP titer of 71.09 g/L was achieved at 100 h when the biomass reached 71.6 g/L, which was higher than the reported titer achieved through cytosol-based production so far. The 3-HP productivity reached 0.71 g/L/h, indicating its potential applications in industrial production. As shown in Additional file 3: Table S1, the titer and the productivity of N3IP_2 represented the highest values reported in S. cerevisiae via the malonyl-CoA pathway, suggesting that the yeast mitochondrion could be a suitable subcellular compartment for 3-HP production Yet, the yields on glucose of the strain N3IP_2 were 0.31 g/g and 0.23 g/g in the shake flask and fed-batch fermentation, respectively, compared with the maximum theoretical yield of 0.86 g/g. During the fed-batch cultivations, we observed that the 3-HP yield on glucose could reach as high as 0.40 g/g according to 3-HP production versus glucose consumption, suggesting that better feeding might improve 3-HP yields. Moreover, the yields could be further improved with better-tuned mitochondrial metabolism. As observed in this study, the flux distribution among the TCA cycle and 3-HP synthesis might require further investigation to balanced minimal growth requirements and maximum 3-HP production. In addition, it would be interesting to explore potential transporters involved in 3-HP transportation across the mitochondrial membranes and cytomembranes, as efficient transportation might relieve product inhibition and facilitate product synthesis.

Fig. 6
figure 6

Fed-batch fermentation profiles of N3IP_2 under glucose limited conditions. Blue circles, dry cell weight (DCW); red circles, glucose concentration; grey circles, ethanol concentration; green circles, 3-HP titer; yellow circles, consumed glucose. The cultivations were performed in biologically duplicate and error bars represent ± standard errors

Conclusions

Metabolic modeling revealed that the malonyl-CoA pathway targeted into the mitochondria could synthesize 3-HP with a higher yield and a lower oxygen consumption. Overexpression of the mitochondrial MCR enzymes obtained the comparable 3-HP titer and yield with the highest level reported in the cytosol, indicating great potentials in harnessing the mitochondrial for 3-HP production. With further engineering on the mitochondrial NADPH metabolism and flux control at the acetyl-CoA node, the constructed strain obtained increased 3-HP titer by 39.4% in shake flask fermentations, and a final 3-HP titer of 71.06 g/L with a productivity of 0.71 g/L/h in fed-batch fermentations, the highest levels reported in S. cerevisiae so far. The high 3-HP productivity could be greatly credited for the controlled glucose limited conditions during the process, as the derepressed mitochondrial activities not only provide sufficient acetyl-CoA and cofactor but also finely tune flux toward malonyl-CoA for 3-HP synthesis.

Availability of data and materials

The data sets supporting the conclusions of the article are available within this article and its additional files. All data generated or analyzed during this study are included in this published article.

Abbreviations

3-HP:

3-Hydroxypropionate

MCR:

Malonyl-CoA reductase

ACC:

Acetyl-CoA carboxylase

PDH:

Pyruvate dehydrogenase complex

Hfa1:

Mitochondrial acetyl-CoA carboxylase

Idp1:

Mitochondrial NADP-specific isocitrate dehydrogenase

Pos5:

Mitochondrial NADH kinase

Cit1:

Mitochondrial citrate synthase

Aco1/2:

Aconitase

Idh1/2:

Mitochondrial NAD+ -dependent isocitrate dehydrogenase

Mdh1:

Mitochondrial malate dehydrogenase

Mae1:

Mitochondrial malic enzyme

EcMaeA*:

Malic enzyme mutant from E. coli

References

  1. Gambacorta FV, Dietrich JJ, Yan Q, Pfleger BF. Rewiring yeast metabolism to synthesize products beyond ethanol. Curr Opin Chem Biol. 2020;59:182–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Liu H, Zhou P, Qi M, Guo L, Gao C, Hu G, et al. Enhancing biofuels production by engineering the actin cytoskeleton in Saccharomyces cerevisiae. Nat Commun. 2022;13(1):1886.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Zhang Y, Nielsen J, Liu Z. Yeast based biorefineries for oleochemical production. Curr Opin Biotechnol. 2021;67:26–34.

    Article  CAS  PubMed  Google Scholar 

  4. Liu Q, Liu Y, Chen Y, Nielsen J. Current state of aromatics production using yeast: achievements and challenges. Curr Opin Biotechnol. 2020;65:65–74.

    Article  CAS  PubMed  Google Scholar 

  5. Duran L, Lopez JM, Avalos JL. ¡Viva la mitochondria!: harnessing yeast mitochondria for chemical production. FEMS Yeast Res. 2020. https://doi.org/10.1093/femsyr/foaa037.

    Article  PubMed  Google Scholar 

  6. Yee DA, DeNicola AB, Billingsley JM, Creso JG, Subrahmanyam V, Tang Y. Engineered mitochondrial production of monoterpenes in Saccharomyces cerevisiae. Metab Eng. 2019;55:76–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Qin J, Zhou YJ, Krivoruchko A, Huang M, Liu L, Khoomrung S, et al. Modular pathway rewiring of Saccharomyces cerevisiae enables high-level production of L-ornithine. Nat Commun. 2015;6:8224.

    Article  PubMed  Google Scholar 

  8. Yao Z, Zhou P, Su B, Su S, Ye L, Yu H. Enhanced isoprene production by reconstruction of metabolic balance between strengthened precursor supply and improved isoprene synthase in Saccharomyces cerevisiae. ACS Synth Biol. 2018;7(9):2308–16.

    Article  CAS  PubMed  Google Scholar 

  9. Tao XY, Lin YC, Wang FQ, Liu QH, Ma YS, Liu M, et al. Production of sesquiterpene patchoulol in mitochondrion-engineered Saccharomyces cerevisiae. Biotechnol Lett. 2022;44(4):571–80.

    Article  CAS  PubMed  Google Scholar 

  10. Yuan J, Ching CB. Mitochondrial acetyl-CoA utilization pathway for terpenoid productions. Metab Eng. 2016;38:303–9.

    Article  CAS  PubMed  Google Scholar 

  11. Zhang Y, Wang J, Cao X, Liu W, Yu H, Ye L. High-level production of linalool by engineered Saccharomyces cerevisiae harboring dual mevalonate pathways in mitochondria and cytoplasm. Enzyme Microb Technol. 2020;134:109462.

    Article  CAS  PubMed  Google Scholar 

  12. Dong C, Shi Z, Huang L, Zhao H, Xu Z, Lian J. Cloning and characterization of a panel of mitochondrial targeting sequences for compartmentalization engineering in Saccharomyces cerevisiae. Biotechnol Bioeng. 2021;118(11):4269–77.

    Article  CAS  PubMed  Google Scholar 

  13. Hammer SK, Zhang Y, Avalos JL. Mitochondrial compartmentalization confers specificity to the 2-ketoacid recursive pathway: increasing isopentanol production in Saccharomyces cerevisiae. ACS Synth Biol. 2020;9(3):546–55.

    Article  CAS  PubMed  Google Scholar 

  14. Hammer SK, Avalos JL. Harnessing yeast organelles for metabolic engineering. Nat Chem Biol. 2017;13(8):823–32.

    Article  CAS  PubMed  Google Scholar 

  15. Liang B, Sun G, Zhang X, Nie Q, Zhao Y, Yang J. Recent advances, challenges and metabolic engineering strategies in the biosynthesis of 3-hydroxypropionic acid. Biotechnol Bioeng. 2022;119(10):2639–68.

    Article  CAS  PubMed  Google Scholar 

  16. Matsakas L, Hrůzová K, Rova U, Christakopoulos P. Biological production of 3-hydroxypropionic acid: an update on the current status. Fermentation. 2018;4(1):13.

    Article  Google Scholar 

  17. de Fouchecour F, Sanchez-Castaneda AK, Saulou-Berion C, Spinnler HE. Process engineering for microbial production of 3-hydroxypropionic acid. Biotechnol Adv. 2018;36(4):1207–22.

    Article  PubMed  Google Scholar 

  18. Yu W, Cao X, Gao J, Zhou YJ. Overproduction of 3-hydroxypropionate in a super yeast chassis. Bioresour Technol. 2022;361:127690.

    Article  CAS  PubMed  Google Scholar 

  19. Ji RY, Ding Y, Shi TQ, Lin L, Huang H, Gao Z, et al. Metabolic engineering of yeast for the production of 3-hydroxypropionic acid. Front Microbiol. 2018;9:2185.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Qin N, Li L, Ji X, Li X, Zhang Y, Larsson C, et al. Rewiring central carbon metabolism ensures increased provision of acetyl-CoA and NADPH required for 3-OH-propionic acid production. ACS Synth Biol. 2020;9(12):3236–44.

    Article  CAS  PubMed  Google Scholar 

  21. Chen Y, Bao J, Kim IK, Siewers V, Nielsen J. Coupled incremental precursor and co-factor supply improves 3-hydroxypropionic acid production in Saccharomyces cerevisiae. Metab Eng. 2014;22:104–9.

    Article  PubMed  Google Scholar 

  22. Kildegaard KR, Jensen NB, Schneider K, Czarnotta E, Ozdemir E, Klein T, et al. Engineering and systems-level analysis of Saccharomyces cerevisiae for production of 3-hydroxypropionic acid via malonyl-CoA reductase-dependent pathway. Microb Cell Fact. 2016;15:53.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Galdieri L, Vancura A. Acetyl-CoA carboxylase regulates global histone acetylation. J Biol Chem. 2012;287(28):23865–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Hasslacher M, Ivessa AS, Paltauf F, Kohlwein SD. Acetyl-coa carboxylase from yeast is an essential enzyme and is regulated by factors that control phospholipid-metabolism. J Biol Chem. 1993;268(15):10946–52.

    Article  CAS  PubMed  Google Scholar 

  25. Chen X, Yang X, Shen Y, Hou J, Bao X. Increasing malonyl-CoA derived product through controlling the transcription regulators of phospholipid synthesis in Saccharomyces cerevisiae. ACS Synth Biol. 2017;6(5):905–12.

    Article  CAS  PubMed  Google Scholar 

  26. David F, Nielsen J, Siewers V. Flux control at the malonyl-CoA node through hierarchical dynamic pathway regulation in Saccharomyces cerevisiae. ACS Synth Biol. 2016;5(3):224–33.

    Article  CAS  PubMed  Google Scholar 

  27. Shi S, Chen Y, Siewers V, Nielsen J. Improving production of malonyl coenzyme a-derived metabolites by abolishing snf1-dependent regulation of acc1. MBio. 2014. https://doi.org/10.1128/mBio.01130-14.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Liu CS, Ding YM, Zhang RB, Liu HZ, Xian M, Zhao G. Functional balance between enzymes in malonyl-CoA pathway for 3-hydroxypropionate biosynthesis. Metab Eng. 2016;34:104–11.

    Article  PubMed  Google Scholar 

  29. Liu C, Wang Q, Xian M, Ding Y, Zhao G. Dissection of malonyl-coenzyme a reductase of chloroflexus aurantiacus results in enzyme activity improvement. PLoS ONE. 2013;8(9):e75554.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Hoja U, Marthol S, Hofmann J, Stegner S, Schulz R, Meier S, et al. HFA1 encoding an organelle-specific acetyl-CoA carboxylase controls mitochondrial fatty acid synthesis in Saccharomyces cerevisiae. J Biol Chem. 2004;279(21):21779–86.

    Article  CAS  PubMed  Google Scholar 

  31. Al-Feel W, Chirala SS, Wakil SJ. Cloning of the yeast FAS3 gene and primary structure of yeast acetyl-CoA carboxylase. Proc Natl Acad Sci U S A. 1992;89(10):4534–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Suomi F, Menger KE, Monteuuis G, Naumann U, Kursu VA, Shvetsova A, et al. Expression and evolution of the non-canonically translated yeast mitochondrial acetyl-CoA carboxylase Hfa1p. PLoS ONE. 2014;9(12):e114738.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Miyagi H, Kawai S, Murata K. Two sources of mitochondrial NADPH in the yeast Saccharomyces cerevisiae. J Biol Chem. 2009;284(12):7553–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Minard KI, McAlister-Henn L. Sources of NADPH in yeast vary with carbon source. J Biol Chem. 2005;280(48):39890–6.

    Article  CAS  PubMed  Google Scholar 

  35. Strand MK, Stuart GR, Longley MJ, Graziewicz MA, Dominick OC, Copeland WC. POS5 gene of Saccharomyces cerevisiae encodes a mitochondrial NADH kinase required for stability of mitochondrial DNA. Eukaryot Cell. 2003;2(4):809–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Li YF, Shi F. Partial rescue of pos5 mutants by YEF1 and UTR1 genes in Saccharomyces cerevisiae. Acta Biochim Biophys Sin. 2006;38(5):293–8.

    Article  CAS  PubMed  Google Scholar 

  37. Zhang Y, Wang J, Wang Z, Zhang Y, Shi S, Nielsen J, et al. A gRNA-tRNA array for CRISPR-Cas9 based rapid multiplexed genome editing in Saccharomyces cerevisiae. Nat Commun. 2019;10(1):1053.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Engler C, Marillonnet S. Golden gate cloning. Methods Mol Biol. 2014;1116:119–31.

    Article  CAS  PubMed  Google Scholar 

  39. Zhang Y, Su M, Wang Z, Nielsen J, Liu Z. Rewiring regulation on respiro-fermentative metabolism relieved crabtree effects in Saccharomyces cerevisiae. Synth Syst Biotechnol. 2022;7(4):1034–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Shi S, Si T, Liu Z, Zhang H, Ang EL, Zhao H. Metabolic engineering of a synergistic pathway for n-butanol production in Saccharomyces cerevisiae. Sci Rep. 2016;6:25675.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhang Y, Su M, Qin N, Nielsen J, Liu Z. Expressing a cytosolic pyruvate dehydrogenase complex to increase free fatty acid production in Saccharomyces cerevisiae. Microb Cell Fact. 2020;19(1):226.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Mikkelsen MD, Buron LD, Salomonsen B, Olsen CE, Hansen BG, Mortensen UH, et al. Microbial production of indolylglucosinolate through engineering of a multi-gene pathway in a versatile yeast expression platform. Metab Eng. 2012;14(2):104–11.

    Article  CAS  PubMed  Google Scholar 

  43. Liu T, Gou Y, Zhang B, Gao R, Dong C, Qi M, et al. Construction of ajmalicine and sanguinarine de novo biosynthetic pathways using stable integration sites in yeast. Biotechnol Bioeng. 2022;119(5):1314–26.

    Article  CAS  PubMed  Google Scholar 

  44. Chen Y, Zhang Y, Siewers V, Nielsen J. Ach1 is involved in shuttling mitochondrial acetyl units for cytosolic C2 provision in Saccharomyces cerevisiae lacking pyruvate decarboxylase. FEMS Yeast Res. 2015. https://doi.org/10.1093/femsyr/fov015.

    Article  PubMed  Google Scholar 

  45. Verduyn C, Postma E, Scheffers WA, Vandijken JP. Effect of benzoic acid on metabolic fluxes in yeasts: a continuous-culture study on the regulation of respiration and alcoholic fermentation. Yeast. 1992;8(7):501–17.

    Article  CAS  PubMed  Google Scholar 

  46. Orth JD, Thiele I, Palsson BØ. What is flux balance analysis? Nat Biotechnol. 2010;28(3):245–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Lu H, Li F, Sanchez BJ, Zhu Z, Li G, Domenzain I, et al. A consensus S. cerevisiae metabolic model Yeast8 and its ecosystem for comprehensively probing cellular metabolism. Nat Commun. 2019;10(1):3586.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Heirendt L, Arreckx S, Pfau T, Mendoza SN, Richelle A, Heinken A, et al. Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0. Nat Protoc. 2019;14(3):639–702.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Kayikci O, Nielsen J. Glucose repression in Saccharomyces cerevisiae. FEMS Yeast Res. 2015. https://doi.org/10.1093/femsyr/fov068.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Outten CE, Culotta VC. A novel NADH kinase is the mitochondrial source of NADPH in Saccharomyces cerevisiae. EMBO J. 2003;22(9):2015–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhu J, Schwörer S, Berisa M, Kyung YJ, Ryu KW, Yi J, et al. Mitochondrial NADP(H) generation is essential for proline biosynthesis. Science. 2021;372(6545):968–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Zelle RM, Harrison JC, Pronk JT, van Maris AJ. Anaplerotic role for cytosolic malic enzyme in engineered Saccharomyces cerevisiae strains. Appl Environ Microbiol. 2011;77(3):732–8.

    Article  CAS  PubMed  Google Scholar 

  53. Zhao WN, McAlister-Henn L. Expression and gene disruption analysis of the isocitrate dehydrogenase family in yeast. Biochemistry. 1996;35(24):7873–8.

    Article  CAS  PubMed  Google Scholar 

  54. Haselbeck RJ, McAlister-Henn L. Function and expression of yeast mitochondrial NAD- and NADP-specific isocitrate dehydrogenases. J Biol Chem. 1993;268(16):12116–22.

    Article  CAS  PubMed  Google Scholar 

  55. Sazanov LA, Jackson JB. Proton-translocating transhydrogenase and NAD- and NADP-linked isocitrate dehydrogenases operate in a substrate cycle which contributes to fine regulation of the tricarboxylic acid cycle activity in mitochondria. FEBS Lett. 1994;344(2–3):109–16.

    Article  CAS  PubMed  Google Scholar 

  56. Contreras-Shannon V, Lin AP, McCammon MT, McAlister-Henn L. Kinetic properties and metabolic contributions of yeast mitochondrial and cytosolic NADP+-specific isocitrate dehydrogenases. J Biol Chem. 2005;280(6):4469–75.

    Article  CAS  PubMed  Google Scholar 

  57. Ricci-Tam C, Ben-Zion I, Wang J, Palme J, Li A, Savir Y, et al. Decoupling transcription factor expression and activity enables dimmer switch gene regulation. Science. 2021;372(6539):292–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank for Prof. Shuobo Shi and Prof. Irina Borodina for kindly providing the plasmids and the support from the Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, and the Knut and Alice Wallenberg Foundation.

Funding

The authors thank for the support from the National Key Research and Development Program of China (2018YFA0900201), National Natural Science Foundation of China (21808008 and 22211530047), Tianjin Synthetic Biotechnology Innovation Capacity Improvement Project (TSBICIP-KJGG-009).

Author information

Authors and Affiliations

Authors

Contributions

YZ contributed to conceptualization, investigation, formal analysis, and manuscript writing, review & editing. MS contributed to investigation, formal analysis, and manuscript review & editing. YC contributed to metabolic modeling analysis, and manuscript review & editing. ZW contributed to fed-batch fermentation, and manuscript review & editing. JN contributed to conceptualization, supervision, formal analysis, and manuscript review & editing. ZL contributed to conceptualization, supervision, formal analysis and manuscript review & editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Zihe Liu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. Plasmids, strains and primers constructed and used in the study.

Additional file 2

: Table S2. Flux distribution from flux balance analysis with cytosolic MCR and mitochondrial MCR expressed.

Additional file 3

: Fig. S1. Varied gene expression levels characterized by FPKM values. Fig. S2. Fermentation profiles of the strain N3IP cit1. Fig. S3. Fermentation profiles of the strain N3IP_pACC1**_TEF1p. Table S1. Summary of 3-HP bioproduction via malonyl-CoA pathway in S. cerevisiae.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Su, M., Chen, Y. et al. Engineering yeast mitochondrial metabolism for 3-hydroxypropionate production. Biotechnol Biofuels 16, 64 (2023). https://doi.org/10.1186/s13068-023-02309-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-023-02309-z

Keywords