- Open Access
13C-metabolic flux analysis of ethanol-assimilating Saccharomyces cerevisiae for S-adenosyl-l-methionine production
Microbial Cell Factories volume 17, Article number: 82 (2018)
Saccharomyces cerevisiae is a host for the industrial production of S-adenosyl-l-methionine (SAM), which has been widely used in pharmaceutical and nutritional supplement industries. It has been reported that the intracellular SAM content in S. cerevisiae can be improved by the addition of ethanol during cultivation. However, the metabolic state in ethanol-assimilating S. cerevisiae remains unclear. In this study, 13C-metabolic flux analysis (13C-MFA) was conducted to investigate the metabolic regulation responsible for the high SAM production from ethanol.
The comparison between the metabolic flux distributions of central carbon metabolism showed that the metabolic flux levels of the tricarboxylic acid cycle and glyoxylate shunt in the ethanol culture were significantly higher than that of glucose. Estimates of the ATP balance from the 13C-MFA data suggested that larger amounts of excess ATP was produced from ethanol via increased oxidative phosphorylation. The finding was confirmed by the intracellular ATP level under ethanol-assimilating condition being similarly higher than glucose.
These results suggest that the enhanced ATP regeneration due to ethanol assimilation was critical for the high SAM accumulation.
Saccharomyces cerevisiae has been used in several industrial processes such as for the production of S-adenosyl-l-methionine (SAM) [1, 2]. SAM is synthesized from l-methionine and ATP in S. cerevisiae, and acts as a biological methyl group donor involved in many metabolic reactions such as transmethylation of proteins. The metabolite has been widely used in pharmaceutical and nutritional supplement industries . For SAM production, S. cerevisiae strains used for Japanese sake brewing (Kyokai strains) are suitable owing to these higher intracellular contents among microorganisms . The metabolic analysis showed that enhanced energy (ATP) regeneration in Kyokai strains contribute to high SAM production .
It has been reported that the SAM production was improved by feeding with ethanol. For instance, the Kyokai strains produce 10.8 g/L of SAM in a 10-L fermenter under ethanol feed conditions . While the transcriptome analysis of S. cerevisiae grown on ethanol indicated the transcripts related to gluconeogenesis, the glyoxylate shunt, and the tricarboxylic acid (TCA) cycle were upregulated compared with those grown on glucose , the in vivo activity of the metabolic pathway involved in the SAM accumulation under ethanol-assimilating condition is still far from clear.
To further understand the mechanism in response to SAM production from ethanol, we performed a metabolic flux analysis for central carbon metabolism since SAM is biosynthesized from precursors to key intermediates including ribose 5-phosphate, ATP, oxaloacetate, and l-methionine using reductive power (NADPH) and chemical motive force (ATP). While a metabolic flux distribution in the ethanol-assimilating S. cerevisiae was estimated through metabolic flux balancing using a compartmented stoichiometric model in silico , the experimental determination by 13C-metabolic flux analysis (13C-MFA) is required to comprehensively quantify central carbon metabolism on ethanol in SAM-producing S. cerevisiae . To our knowledge, there have been no prior 13C-MFA studies on ethanol metabolism in S. cerevisiae.
In the present study, 13C-MFA was conducted on a high-SAM-producing S. cerevisiae strain with 13C-labeled ethanol as the sole carbon source. 13C-MFA is based on the cultivation of the cells in medium containing 13C-labeled carbon sources. A metabolic flux distribution was estimated from the 13C-labeling patterns of intracellular metabolites using mass spectrometry. The selection of a suitable 13C-labeled carbon source is critical in the design of a 13C-MFA experiment since the accuracy of the flux estimation depends on the labeling patterns of the carbon source. However, a suitable experimental design remains unclear for the 13C-MFA of the central carbon metabolism in S. cerevisiae cultured with ethanol as the sole carbon source. Thus the design of the 13C-MFA experiment was optimized in this study by a computer simulation and found that 100% [2-13C] ethanol was the best carbon source for the precise measurements of metabolic fluxes. Using the optimized experimental design, the intracellular metabolic flux distribution of central carbon metabolism was successfully determined for the high SAM-producing S. cerevisiae strain (Kyokai no. 6) in the ethanol limited-chemostat culture. The results showed that the metabolic flux levels through the glyoxylate shunt and the later TCA cycle were upregulated during growth on ethanol, and the resultant activation of the oxidative phosphorylation should contribute the high SAM accumulation in the ethanol-assimilating conditions.
Strain and growth conditions
The S. cerevisiae strain used for Japanese sake brewing, Kyokai no. 6 (NBRC2346), was purchased from the National Biological Resource Center (NBRC, Chiba, Japan). This S. cerevisiae strain was cultivated in an aerobic carbon-limited chemostat culture with a working volume of 100 mL in a 250-mL fermentor (ABLE Co., Tokyo, Japan). The pH was maintained at 5.5 by the automatic addition of 1.0 N NaOH. The temperature was maintained at 30 °C. The stirring speed was 1000 rpm. The aeration rate was 200 mL/min. The dilution rate was 0.06 h−1. The synthetic medium used for cultivation contained 3.8 g/L ethanol, 5.0 g/L (NH4)2SO4, 0.50 g/L K2HPO4, and 3.4 g/L of yeast nitrogen base without amino acids and ammonium sulfate (Difco Laboratories, Franklin Lakes, NJ, USA). Vitamins and minerals were added to the medium, with the final composition as follows: 2.0 g/L KH2PO4, 1.5 g/L MgSO4, 2.5 mg/L ZnSO4, 2.4 mg/L MnSO4, 0.27 mg/L CuSO4, 0.20 g/L CaCl2, 4.0 mg/L FeCl3, 0.40 mg/L NaMoO4∙2H2O, 1.0 mg/L H3BO3, 0.20 mg/L KI, 0.20 g/L NaCl, 34 µg/L biotin, 1.6 mg/L Ca-pantothenate, 13 mg/L inositol, 7.8 mg/L thiamine-HCl, 2.3 mg/L pyridoxine-HCl, 0.40 mg/L para-aminobenzoic acid, 0.40 mg/L riboflavin, 0.80 mg/L niacin, and 4.0 µg/L folic acid. For 13C-MFA, the carbon source in the medium was replaced with [2-13C]ethanol. Samples taken from the reactor were centrifuged at 18,800×g for 5 min at 4 °C. The supernatant was used for analysis of extracellular metabolites. The cell pellet was used for 13C-labeling analysis by gas chromatography-mass spectrometry (GC–MS).
The analysis of SAM, cell dry weight (CDW), glucose, ethanol, glycerol and organic acids such as acetate were performed as described by Hayakawa et al. . SAM concentration was measured by high-performance liquid chromatography (HPLC) LC2010A-HT (Shimadzu, Kyoto, Japan) after 10% HClO4 extraction. CDW was estimated using OD600 = 1 corresponding to 0.21 gCDW/L. Glucose and ethanol concentrations were measured enzymatically by using an analyzer (BF-7, Oji Scientific Instruments, Hyogo, Japan). An enzymatic kit was used for the quantification of glycerol (F-kit Glycerol, R-Biopharm, Washington, MO, USA) in the supernatant. The concentrations of acetate, pyruvate, succinate, fumarate, malate, citrate, and lactate in the supernatant were measured by HPLC (LC2010A-HT, Shimadzu).
The intracellular ATP, l-methionine, and SAM concentrations were measured using capillary electrophoresis time-of-flight mass spectrometry (CE-TOFMS) according to previously described methods . Cultivated cells were harvested by filtration and washed with Milli-Q water. After the membrane was soaked in methanol solution containing an internal standard (H3304-1002, Human Metabolome Technologies, Yamagata, Japan), ultra-sonication was performed. After removal of the membrane, chloroform and water were added to the solution. The aqueous portion collected after mixing with a vortex mixer was filtered by Ultrafree MC-PLHCC 250 (Human Metabolome Technologies) and dried. For CE-TOFMS, pellets were suspended in the second internal standard (H3304-1004, Human Metabolome Technologies) solution. Samples were analyzed using the Agilent 7100 CE system with Agilent 6224 TOF–MS (Agilent Technologies, Santa Clara, CA, USA).
Optimization of the 13C tracer for 13C-MFA by computational simulation
A computer simulation was conducted for the design of the 13C-MFA experiment as previously described . The intracellular flux distribution data was determined using the modified literature data through metabolic flux balancing . Data for the consumption and production rates were obtained in this study. A composition of [13C]ethanol was predetermined. The simulated 13C-enrichments were calculated for the 21 fragments of the following amino acids: (M−57)+ fragments for Ala, Asp, Glu, Gly, Phe, and Thr; (M−85)+ fragments for Ala, Asp, Glu, Gly, Ile, Leu, Pro, Thr, and Val; (M−159)+ fragments for Glu, Ile, Leu, Pro, and Val; and (M302)+ fragment for Asp; and then the Gaussian noise at 1% levels was added to produce hypothetical measurements for the 13C-enrichment data. A metabolic flux distribution was estimated by minimizing the residual between the hypothetical measurement and simulated 13C-enrichments of amino acid fragments. The 95% confidence intervals were determined using the grid search method [10, 11]. Precision scoring was calculated by the following equation to evaluate the range of the 95% confidence intervals:
where r i , r[10:0:0:0], i, and S i are the range of the 95% confidence intervals, the range of 95% confidence intervals in 13C-labeled ethanol composition of 10:0:0:0 (100% non-labeled ethanol), and precision score for the ith flux. The maximum net flux was 500, which was normalized to an ethanol consumption rate of 100. The metabolic network shown in the following was used.
The metabolic model of S. cerevisiae for the 13C-MFA was based on previously published models, including glycolysis, the pentose phosphate (PP) pathway, the anaplerotic pathways, the TCA cycle, gluconeogenesis, and the transport reactions between the cytosol and mitochondria (Additional file 1: Tables S1 and S2) [4, 12, 13]. In the metabolic model, the reaction from the mitochondrial oxaloacetate and acetyl-CoA to isocitrate via citrate was assumed to be reversible. The amino acid biosynthesis pathways, including two pathways for glycine and alanine biosynthesis, were employed (Additional file 1: Table S2) . The composition of the S. cerevisiae biomass was determined on basis of the literature (Additional file 1: Tables S3–S6) [15, 16]. Using these data and the stoichiometry of anabolic metabolism, the demands of precursor metabolites for cell growth were calculated (Additional file 1: Table S7) and used as constraint conditions in 13C-MFA (Additional file 1: Table S8).
GC–MS analysis of proteinogenic amino acids
Cell pellets were hydrolyzed in 6 N HCl at 105 °C for 18 h. After removal of the debris by filtration, the hydrolysate was dried and dissolved in acetonitrile. For GC–MS, the hydrolysate was mixed with an equal volume of N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide (MTBSTFA), and the mixture was incubated at 95 °C for 1 h. Next, 1 μL of the sample was injected into the GC–MS system (7890A GC and 5975C GC/MSD, Agilent Technologies, Santa Clara, CA, USA), according to a previously described method [17, 18]. The GC–MS data were corrected considering the natural abundance of C, H, N, O, and Si isotopes for 13C-MFA . The 13C-enrichments of the amino acid fragments at an isotopic steady state [χ(∞)] can be calculated from the following equation:
In Eq. (2), t represents the time that has elapsed since feeding 13C-labeled ethanol-containing medium, χ(t) and χ(0) are the 13C-enrichments of amino acid at t and start time since feeding the 13C-labeled ethanol-containing medium, respectively, and, µ is the dilution rate .
13C-metabolic flux analysis
The computational procedure for 13C-MFA was performed using a Python version of OpenMebius implemented in Python 2.7.9 [9, 21], by which 22 independent fluxes were iteratively tuned by minimizing the residual between the experimental and simulated 13C-enrichments of proteinogenic amino acid fragments. Furthermore, the 13C-labeling patterns of CO2 were independently optimized in the flux estimation. Amino acid fragments used for the flux optimization were shown above. Nonlinear optimization was performed using the SLSQP (sequential least squares programming) function implemented in PyOpt 1.2 . The 90% confidence intervals were determined using the grid search method [10, 11].
NADPH and ATP demands for cell growth and SAM production
The NADPH demands for cell growth in assimilating glucose and ethanol were calculated by using Additional file 1: Tables S3–S6 and the KEGG database (http://www.genome.jp/kegg/) (9.42 and 11.95 mmol/gCDW, respectively). The literature data were used for the ATP demands for cell growth on glucose and ethanol (39.78 and 105.56 mmol/gCDW, respectively) . A non-growth-associated maintenance ATP requirement of 1.0 mmol ATP/gCDW/h was also considered . The total ATP cost for SAM biosynthesis was found to be 13.0 mol ATP per mol SAM .
Physiological parameters of ethanol-assimilating S. cerevisiae
The S. cerevisiae strain used for Japanese sake brewing (Kyokai no. 6, NBRC2346) was cultivated in an ethanol-limited chemostat culture under aerobic condition at a dilution rate of 0.06 h−1. The physiological parameters in the cultivation are shown in Table 1. A comparison with the previous glucose-limited chemostat culture  showed that the SAM specific production rate and intracellular content in assimilating ethanol (0.27 μmol/gCDW/h and 1.8 mg/gCDW, respectively) were greater than that of a glucose-limited culture under the same dilution and carbon atom supply rate conditions (0.089 μmol/gCDW/h and 0.58 mg/gCDW, respectively). The biomass yield in assimilating ethanol was lower than that of glucose as with previous research [25, 26]. In contrast, the acetate yield in assimilating ethanol was higher than that of glucose. This might be related to the enhanced expression of acetaldehyde dehydrogenase (ALD4 and ALD6) under assimilating ethanol conditions .
Design of 13C-MFA using 13C-labeled ethanol
As mentioned in the background, 13C-MFA is based on the cultivation of cells in medium containing 13C-labeled carbon sources. Selection of a suitable 13C-labeled carbon source is critical in the design of a 13C-MFA experiment since the accuracy of the flux estimation depends on labeling patterns of the carbon source. In this study, a suitable experiment was designed for the 13C-MFA of the central carbon metabolism in S. cerevisiae cultured on ethanol as sole carbon source. A relationship between the composition of 13C-labeled ethanol including [1-13C], [2-13C], and [U-13C]ethanol, and 95% confidence intervals of the estimated flux levels were investigated by a computer simulation of 13C-MFA experiment.
The metabolic model of S. cerevisiae describing the stoichiometry equations and carbon atom transitions of each metabolic reaction was constructed from previously published models with modifications (Additional file 1: Tables S1 and S2) [4, 12, 13]. The computer simulation was conducted by the following procedure : using a flux distribution of ethanol-assimilating S. cerevisiae estimated by an in silico calculation method , a series of artificial mass spectra datasets of amino acid fragments were generated for various compositions of [13C]ethanol (all patterns of non-labeled, [1-13C], [2-13C], and [U-13C]ethanol with a 50% step size) by adding Gaussian noise at 1% levels. The metabolic flux distribution and 95% confidence intervals were determined using the artificial mass spectra datasets by the 13C-MFA procedure (Additional file 1: Table S9). A precision score S i was determined for each metabolic reaction i from the data of the 95% confidence intervals. Useful compositions of [13C]ethanol were investigated by comparing the precision score S i of estimated flux distributions.
For example, the 95% confidence intervals of isocitrate dehydrogenase (IDH) flux (isocitrate→α-ketoglutarate + CO2) with mixtures of non-labeled, [1-13C], [2-13C], [U-13C]ethanol at 10:0:0:0, 5:5:0:0, and 0:0:10:0 were estimated to be 2.5–45.1, 19.0–42.3, and 27.4–41.6, from which S IDH were estimated to be 1.0, 0.55, and 0.33, respectively (Additional file 1: Table S9). Since a smaller S IDH score shows a more precise flux estimation, the results indicated that the mixture of non-labeled, [1-13C], [2-13C], [U-13C]ethanol at 0:0:10:0 was a suitable carbon sources for the flux estimation of this reaction in the three mixtures. Figure 1 shows a heatmap in the precision score S i based on the range of 95% confidence interval for the estimated flux in each composition of 13C-labeled ethanol. The magenta and green colors represent narrower (better precision) and wider (poorer precision) 95% confidence intervals, respectively. This result indicates that the range of the 95% confidence interval depended greatly on the composition of 13C-labeled ethanol. The mixture ratio at 0:0:10:0 (100% [2-13C]ethanol) showed the smallest sum of the precision score Si in the TCA cycle and PP pathway among all 13C-labeling conditions (Fig. 1). Based on the results, [2-13C]ethanol was employed for the 13C-MFA of the central carbon metabolism in S. cerevisiae cultured on ethanol as the sole carbon source.
13C-Metabolic flux analysis
The Kyokai no. 6 strain was cultivated for 13C-MFA by an ethanol-limited chemostat culture under aerobic conditions at a dilution rate 0.06 h−1. At 98 h after the start of cultivation when the cell growth reached steady-state, the medium containing 100% [2-13C]ethanol was fed to the chemostat cultures. The biomass samples were harvested at 1.3, 2.3, 2.7, 3.7, 4.1, and 5.1 residential times after feeding 13C-labeled ethanol. Following acid hydrolysis and derivatization, 13C-enrichment of proteinogenic amino acids was measured using GC–MS analysis (Additional file 1: Fig. S1). The intracellular metabolism might reach an isotopic steady-state after an infinite amount of time. In this study, the 13C-enrichment data indicated by Eq. (2) were used to calculate metabolic flux distributions (Additional file 1: Table S10 and Fig. S1).
A metabolic flux distribution was estimated by minimizing the difference between the computationally simulated 13C-enrichments of proteinogenic amino acids and the experimentally obtained data by GC–MS. The differences (residual sum of squares) between experimental and simulated 13C-enrichment data for the best-fitted metabolic flux distributions that passed the χ2-test (α = 0.05) were small (69.5) , indicating that the estimated flux distributions could explain the experimentally obtained 13C-enrichments (Additional file 1: Table S10). Based on the best-estimate metabolic flux distribution, 90% confidence intervals of the metabolic flux levels were determined using the grid search method [10, 11].
The 90% confidence intervals of the flux levels on glucose and ethanol in the Kyokai no. 6 strain at a dilution rate 0.06 h−1 are shown in Fig. 2. The results for glucose have been previously published . The results revealed that the metabolic flux distribution was totally rewired in the ethanol culture from that under the glucose culture. The Embden–Meyerhof–Parnas (EMP) pathway was reversed from the glycolytic direction in the glucose culture to the direction for gluconeogenesis because the assimilated ethanol flowed into the EMP pathway via the glyoxylate shunt and then converted to glucose-6-phosphate. The metabolic flux in the glyoxylate shunt (IsoCit→Glxy + Suc, AcCOA_c + Glxy→Oxa_c) for the glucose and ethanol assimilations were 0.007–0.03 and 0.58–0.68 mmol/gCDW/h, respectively. This was in agreement with the gene transcription levels of ICL1, ICL2, MLS1, and MLS2 in the glyoxylate shunt that were upregulated in a transcriptome analysis of S. cerevisiae grown on ethanol described previously . Since the expression of malate synthase (AcCOA_c + Glxy→Oxa_c) was essential for growth on ethanol , the result confirmed that the pathway was active during the assimilation of ethanol.
To investigate the cofactor balance under the ethanol-assimilating conditions, metabolic flux levels responsible for the NADPH and ATP regeneration were compared (Fig. 3a, b). In S. cerevisiae, NADPH was mainly regenerated via the PP pathway (glucose-6-phosphate dehydrogenase and phosphogluconate dehydrogenase), malic enzyme, and isocitrate dehydrogenase. The results showed that the TCA cycle functioned as the main provider of NADPH on ethanol since the flux level of isocitrate dehydrogenase (IsoCit + NADP+→aKG + NADPH + CO2) in the TCA cycle was higher on ethanol (0.91–1.1 mmol/gCDW/h) compared to glucose (0.34–0.42 mmol/gCDW/h) (Figs. 2 and 3a). On the other hand, the glucose-6-phosphate dehydrogenase and phosphogluconate dehydrogenase flux (G6P + 2 NADP+→Ru5P + 2 NADPH + CO2) in the oxidative PP pathway was 0–0.25 mmol/gCDW/h on ethanol, which was lower compared to the 0.36–0.58 mmol/gCDW/h on glucose, indicating that the carbon fluxes of the central carbon metabolism on ethanol were redirected into the TCA cycle (Fig. 2). The malic enzyme flux (Mal + NADP+→Pyr_m + NADPH + CO2) was similar on glucose and ethanol (0.04–0.08 and 0–0.12 mmol/gCDW/h, respectively), which showed that this reaction did not contribute greatly in the regeneration of NADPH on both carbon sources (Fig. 2).
Figure 3b shows the ATP balance in the glucose- and ethanol-assimilating conditions. The data revealed that oxidative phosphorylation using the electron transport system during respiration was a major source of ATP regeneration. The metabolic flux distribution showed that the NADH supply for the electron transport system are from the EMP pathway, TCA cycle, glyoxylate shunt, alcohol dehydrogenase, and aldehyde dehydrogenase. The flux levels of these reactions were also largely different between glucose and ethanol. For instance, the malate dehydrogenase flux (Mal + NAD+→Oxa_m + NADH) in the TCA cycle on ethanol corresponded to 59–65% of ethanol consumption rate (1.5–1.6 mmol/gCDW/h), while the flux level on glucose corresponded to 35–48% of glucose consumption rate (0.23–0.31 mmol/gCDW/h) (Fig. 2). In contrast, the pyruvate dehydrogenase (Pyr_m + NAD+→AcCOA_m + CO2 + NADH) in the EMP pathway was inactive on ethanol because of a low flux (0–0.06 mmol/gCDW/h), indicating that the TCA cycle also functioned as the provider of NADH on ethanol (Fig. 2). Therefore, this result suggests that the high TCA activity in ethanol assimilation caused the lower cell yield on carbon source compared with glucose assimilation, since two CO2 molecules were released from the TCA cycle (Table 1).
Next, the carbon flows around cytosolic oxaloacetate, which is the precursor for SAM, were compared. The levels of SAM biosynthesis flux from cytosolic oxaloacetate on glucose- and ethanol-assimilating (0.09 and 0.3 µmol/gCDW/h, respectively) were only 0.01% or less of the glucose- and ethanol-consumption flux values (0.65 and 2.5 mmol/gCDW/h, respectively) (Fig. 2). In turn, the flux levels of pyruvate carboxylase (Pyr_c + CO2→Oxa_c) on glucose and ethanol were 0.14–0.18 and 0.42–1.4 mmol/gCDW/h, respectively, indicating that pyruvate carboxylase were used to supply cytosolic oxaloacetate on both carbon sources (Fig. 2). This result was confirmed by the detected enzyme activity under the ethanol-grown culture .
As mentioned above, 13C-MFA showed that the flux levels related to NADPH and ATP regeneration were largely different during growth on glucose and ethanol, while the flux of cytosolic oxaloacetate to SAM biosynthesis corresponded to 0.01% of the ethanol consumption flux value. Hence, the relationship between cofactor (NADPH and ATP) balances and SAM biosynthesis was investigated. To analyze the NADPH balance, the NADPH regeneration rates estimated from the metabolic flux data were compared with the NADPH consumption rates for cell growth. Since the composition of S. cerevisiae biomass differs among the carbon source used for the culture, NADPH demand levels were determined using Additional file 1: Tables S3–S6 and the KEGG database (see Materials and methods). As shown in Fig. 3a, NADPH regeneration rates were larger than the consumption rates for cell growth in both conditions. In assimilating glucose and ethanol, the differences between NADPH regeneration and consumption rates were 0.75 and 0.42 mmol/gCDW/h, respectively, indicating that the difference on glucose was greater than that of ethanol (Fig. 3a). Excess NADPH was assumed to be used for ATP production by the oxidation of NADPH as previously shown .
Since the high levels of SAM production by the Kyokai no. 6 strain could be attributed to enhanced ATP regeneration , ATP balances for regeneration and consumption were investigated (Fig. 3b). Total ATP regeneration rates were calculated based on the metabolic flux of reactions responsible for the respiration and substrate level phosphorylation. The biosynthetic ATP demand was estimated from the literature data (see “Methods”). The ATP consumption rate for cytosolic acetyl-CoA synthesis was calculated from acetyl-CoA synthase flux (Ac→AcCOA_c). The difference between ATP regeneration and consumption on ethanol was 5.1 mmol/gCDW/h, indicating that the difference on ethanol was 1.9-fold greater than that of glucose (2.8 mmol/gCDW/h) (Fig. 3b).
Intracellular ATP, l-methionine, and SAM concentration
The 13C-metabolic flux analysis revealed that SAM accumulation on ethanol should be derived from enhanced ATP regeneration. It was also expected that the intracellular ATP concentration on ethanol could be higher than that of glucose since the SAM content in assimilating ethanol at a dilution rate of 0.06 h−1 (1.8 mg/gCDW) (Table 1) was greater than that of glucose at a dilution rate of 0.1 h−1 (1.0 mg/gCDW) . Then, intracellular ATP concentrations were determined by CE-TOFMS analysis to further investigate the effects of ATP levels on SAM productivity. For the comparison of ATP levels, glucose- and ethanol-limited chemostat cultures were conducted at dilution rates of 0.1 and 0.06 h−1, respectively. After the batch phase (at 11 and 44 h after the initiation of cultivation for glucose and ethanol, respectively), the carbon-limited chemostat cultures were started. After reaching steady-state cell growth, the S. cerevisiae cells were repeatedly sampled (at 75, 82, and 99 h on glucose; and at 112, 119, and 136 h on ethanol; respectively) (Fig. 4a). In this study, three samples were considered to be experimental triplicate. The intracellular ATP levels in the cells grown on ethanol (24.8 ± 1.6 μmol/gCDW) was 1.4-fold higher than that of glucose (17.5 ± 2.6 μmol/gCDW) (Fig. 4b). This result confirmed that the S. cerevisiae cells grown on ethanol had the potential for high SAM content due to high intracellular ATP levels as shown in the metabolic flux distribution estimated using 13C-MFA. Similarly, the intracellular l-methionine and SAM levels on ethanol (0.10 ± 0.0019 μmol/gCDW and 1.2 ± 0.067 μmol/gCDW) were 1.7- and 1.2-fold higher than those of glucose (0.060 ± 0.0024 μmol/gCDW and 1.0 ± 0.12 μmol/gCDW), respectively, indicating that improvement of intracellular ATP and l-methionine levels enhanced SAM biosynthesis (Fig. 4c, d).
In the present study, 13C-MFA was conducted to investigate a metabolic regulation responsible for the high intracellular SAM content in ethanol-assimilating S. cerevisiae Kyokai no. 6 strain [5, 28]. The design of the 13C-MFA experiment was optimized by computer simulations, which showed that 100% [2-13C]ethanol was the best carbon source for the precise estimation of metabolic fluxes through the TCA cycle and PP pathway (Fig. 1 and Additional file 1: Table S9). An advantage of having 100% [2-13C]ethanol in 13C-MFA was that the second carbon of ethanol contributes to the generation of various metabolic intermediates with various 13C-labeling pattern, which in turn increases the sensitivity to fluxes. On the other hand, the first and fourth carbons of mitochondrial oxaloacetate derived from the first carbon of ethanol are released as CO2 in the TCA cycle . The 13C-MFA in Escherichia coli cultured on acetate was performed using the mixtures of non-labeled, [2-13C]acetate at a ratio of 8:2 because of the same reason .
A comparison of the metabolic flux distribution on glucose and ethanol revealed that the metabolic flux through the EMP pathway was reversed by the activation of the glyoxylate shunt and gluconeogenesis during growth on ethanol (Fig. 2). The metabolic redirection coincided with the regulation in mRNA and protein expression levels as reported in previous studies. It has been reported that S. cerevisiae grown on ethanol increased the transcript levels of genes involved in gluconeogenesis and the glyoxylate shunt compared to cells grown on glucose . Furthermore, comparable results have been observed in a proteome analysis under chemostat cultures limited for glucose and ethanol . The previous transcriptome and proteome studies during the diauxic shift in S. cerevisiae [32,33,34] also showed that the glyoxylate shunt, gluconeogenesis, and oxidative phosphorylation were activated after starting the assimilation of the produced ethanol.
On the other hand, the downregulation of the pyruvate dehydrogenase flux level (Pyr_m→AcCOA_m + CO2) at the entry point of the TCA cycle disagreed with a previous transcriptome analysis because the expression levels of PDA1, PDB1, PDX1, and LPD1 genes in the ethanol-assimilating conditions increased from that of glucose (Fig. 2) . This may have been due to pyruvate dehydrogenase being allosterically regulated by NADH . These results suggested that the fluxes through the central carbon metabolism in S. cerevisiae were regulated at the transcript, translational, and post-translational levels [6, 31, 36].
It has been reported that the ethanol-assimilating S. cerevisiae cells showed a higher oxygen uptake rate . The results of 13C-MFA confirmed that ATP regeneration via the oxidative phosphorylation significantly increased under the ethanol-assimilating conditions (Fig. 5). Furthermore, estimation of the ATP balance revealed that a higher amount of excess ATP was produced on ethanol (Fig. 3b), since upregulation of ATP regeneration via the oxidative phosphorylation was greater than the increase in ATP consumption for gluconeogenesis and cytosolic acetyl-CoA synthesis during ethanol assimilation . These results suggested that the excess ATP should activate SAM biosynthesis. However, the excess ATP regeneration flux (5.1 mmol/gCDW/h) was clearly greater than the ATP consumption flux for SAM production (3.5 µmol/gCDW/h) (Table 1 and Fig. 3b). The excess ATP may be utilized for the response of stress caused by acetaldehyde and acetate [37, 38].
The analysis of the intracellular metabolites (Fig. 4) revealed that the ATP level in the cells grown on ethanol was higher than that of glucose (Fig. 4b). It was reported that the increase of intracellular ATP level enhanced the SAM production in S. cerevisiae. The addition of sodium citrate improved the isocitric acid dehydrogenase activity and ATP level in the cell, which promoted the conversion of methionine into SAM . Optimization of the culture medium also revealed that the restriction of cell growth, and the enhancement of the intracellular ATP level and SAM production were achieved by reducing the supplemented yeast extract . Moreover, the result of this study that the intracellular l-methionine level in ethanol assimilating was higher than that of glucose (Fig. 4c). Since l-methionine biosynthesis consumes cytosolic acetyl-CoA in S. cerevisiae , the upregulation might be concerned with the higher flux level of acetyl-CoA synthase (Ac→AcCOA_c) on ethanol (2.5–2.6 mmol/gCDW/h) compared to glucose (0.13–0.43 mmol/gCDW/h) (Fig. 2). These results indicated that both ATP and l-methionine levels in cells had a positive effect on SAM biosynthesis (Table 1 and Fig. 4d).
In S. cerevisiae, cytosolic acetyl-CoA was synthesized from acetate and CoA by acetyl-CoA synthase with consuming ATP . The increased ATP regeneration in the ethanol assimilation conditions would also enhance the production of other useful chemicals, since the cytosolic acetyl-CoA is a common precursor for the biosynthesis of fatty acid ethyl ester (biodiesel), geraniol (flavor), and amorpha-4, 11-diene (precursor to artemisinin) [42,43,44]. The results of this study also suggested that the saving on the ATP consumption for cytosolic acetyl-CoA would further increase SAM production. For instance, it has been reported that introducing the metabolic reaction that produces acetyl-CoA from acetaldehyde using an ATP-independent enzyme could increase cell yield .
The 13C-MFA conduced in this study revealed that excess ATP regeneration via the activation of oxidative phosphorylation was a mechanism responsible for the SAM overproduction under the ethanol-assimilating conditions. A detailed mechanism for the metabolic redirection in the central carbon metabolism could be uncovered by a systems-level analysis of the metabolism in ethanol-assimilating S. cerevisiae using 13C-MFA [46, 47].
In this study, the design of the 13C-MFA experiment was optimized by computer simulations, which showed that 100% [2-13C]ethanol was the best carbon source for the precise estimation of metabolic fluxes through the TCA cycle and PP pathway. The 13C-MFA result revealed that the metabolic flux distribution was totally rewired in the ethanol culture from that of glucose. Estimates of the ATP balance from the 13C-MFA data suggested that larger amounts of excess ATP was produced from ethanol via increased oxidative phosphorylation. The finding was confirmed by the intracellular ATP level under ethanol-assimilating condition being similarly higher than glucose. These results suggest that the excess ATP regeneration via the activation of oxidative phosphorylation was a mechanism responsible for the SAM overproduction under the ethanol-assimilating conditions, which also provided general knowledge for the development of S. cerevisiae cell factories for the cytosolic acetyl-CoA-derived products.
Kanai M, Mizunuma M, Fujii T, Iefuji H. A genetic method to enhance the accumulation of S-adenosylmethionine in yeast. Appl Microbiol Biotechnol. 2017;101(4):1351–7.
Shiozaki S, Shimizu S, Yamada H. Unusual intracellular accumulation of S-Adenosyl-l-methionine by microorganisms. Agric Biol Chem. 1984;48(9):2293–300.
Shobayashi M, Mukai N, Iwashita K, Hiraga Y, Iefuji H. A new method for isolation of S-adenosylmethionine (SAM)-accumulating yeast. Appl Microbiol Biotechnol. 2006;69(6):704–10.
Hayakawa K, Kajihata S, Matsuda F, Shimizu H. 13C-metabolic flux analysis in S-adenosyl-l-methionine production by Saccharomyces cerevisiae. J Biosci Bioeng. 2015;120:532–8.
Shiozaki S, Shimizu S, Yamada H. Production of S-adenosyl-l-methionine by Saccharomyces sake. J Biotechnol. 1986;4(6):345–54.
Daran-Lapujade P, Jansen ML, Daran JM, van Pronk W, de Winde JH, Gulik JT. Role of transcriptional regulation in controlling fluxes in central carbon metabolism of Saccharomyces cerevisiae. A chemostat culture study. J Biol Chem. 2004;279(10):9125–38.
Toya Y, Shimizu H. Flux analysis and metabolomics for systematic metabolic engineering of microorganisms. Biotech Adv. 2013;31(6):818–26.
Hayakawa K, Matsuda F, Shimizu H. Metabolome analysis of Saccharomyces cerevisiae and optimization of culture medium for S-adenosyl-l-methionine production. AMB Express. 2016;6(1):38.
Maeda K, Okahashi N, Toya Y, Matsuda F, Shimizu H. Investigation of useful carbon tracers for 13C-metabolic flux analysis of Escherichia coli by considering five experimentally determined flux distributions. Metab Eng Commun. 2016;3:187–95.
van Winden WA, van Dam JC, Ras C, Kleijn RJ, Vinke JL, van Gulik WM, Heijnen JJ. Metabolic-flux analysis of Saccharomyces cerevisiae CEN.PK113-7D based on mass isotopomer measurements of 13C-labeled primary metabolites. FEMS Yeast Res. 2005;5(6–7):559–68.
Antoniewicz MR, Kelleher JK, Stephanopoulos G. Determination of confidence intervals of metabolic fluxes estimated from stable isotope measurements. Metab Eng. 2006;8(4):324–37.
Gombert AK, Moreira dos Santos M, Christensen B, Nielsen J. Network identification and flux quantification in the central metabolism of Saccharomyces cerevisiae under different conditions of glucose repression. J Bacteriol. 2001;183(4):1441–51.
Frick O, Wittmann C. Characterization of the metabolic shift between oxidative and fermentative growth in Saccharomyces cerevisiae by comparative 13C flux analysis. Microb Cell Fact. 2005;4:30.
Monschau N, Stahmann KP, Sahm H, McNeil JB, Bognar AL. Identification of Saccharomyces cerevisiae GLY1 as a threonine aldolase: a key enzyme in glycine biosynthesis. FEMS Microbiol Lett. 1997;150(1):55–60.
Stückrath I, Lange HC, Kötter P, van Gulik WM, Entian KD, Heijnen JJ. Characterization of null mutants of the glyoxylate cycle and gluconeogenic enzymes in S. cerevisiae through metabolic network modeling verified by chemostat cultivation. Biotechnol Bioeng. 2002;77(1):61–72.
Förster J, Famili I, Fu P, Palsson BØ, Nielsen J. Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. Genome Res. 2003;13(2):244–53.
Mori E, Furusawa C, Kajihata S, Shirai T, Shimizu H. Evaluating 13C enrichment data of free amino acids for precise metabolic flux analysis. Biotechnol J. 2011;6(11):1377–87.
Toya Y, Hirasawa T, Morimoto T, Masuda K, Kageyama Y, Ozaki K, Ogasawara N, Shimizu H. 13C-Metabolic flux analysis in heterologous cellulase production by Bacillus subtilis genome-reduced strain. J Biotechnol. 2014;179:42–9.
van Winden WA, Wittmann C, Heinzle E, Heijnen JJ. Correcting mass isotopomer distributions for naturally occurring isotopes. Biotechnol Bioeng. 2002;80(4):477–9.
van Winden W, Schipper D, Verheijen P, Heijnen J. Innovations in generation and analysis of 2D [13C, 1H] COSY NMR spectra for metabolic flux analysis purposes. Metab Eng. 2001;3(4):322–43.
Kajihata S, Furusawa C, Matsuda F, Shimizu H. OpenMebius: an open source software for isotopically nonstationary 13C-based metabolic flux analysis. Biomed Res Int. 2014;2014:627014.
Perez RE, Jansen PW, Martins JRRA. pyOpt: a Python-based object-oriented framework for nonlinear constrained optimization. Struct Multidiscip Optim. 2012;45:101–18.
Verduyn C. Physiology of yeasts in relation to biomass yields. Antonie Van Leeuwenhoek. 1991;60(3–4):325–53.
Verduyn C, Postma E, Scheffers WA, van Dijken JP. Energetics of Saccharomyces cerevisiae in anaerobic glucose-limited chemostat cultures. J Gen Microbiol. 1990;136(3):405–12.
Verduyn C, Stouthamer AH, Scheffers WA, van Dijken JP. A theoretical evaluation of growth yields of yeasts. Antonie Van Leeuwenhoek. 1991;59:49–63.
de Jong-Gubbels P, Vanrolleghem P, Heijnen S, van Dijken JP, Pronk JT. Regulation of carbon metabolism in chemostat cultures of Saccharomyces cerevisiae grown on mixtures of glucose and ethanol. Yeast. 1995;11(5):407–18.
Chen Y, Siewers V, Nielsen J. Profiling of cytosolic and peroxisomal acetyl-CoA metabolism in Saccharomyces cerevisiae. PLoS ONE. 2012;7(8):e42475.
Shiozaki S, Shimizu S, Yamada H. S-Adenosyl-l-methionine production by Saccharomyces sake: optimization of the culture conditions for the production of cells with a high S-Adenosyl-l-methionine content. Agric Biol Chem. 1989;53:3269–74.
Matsuoka Y, Shimizu K. The relationships between the metabolic fluxes and 13C-labeled isotopomer distribution for the flux analysis of the main metabolic pathways. Biochem Eng J. 2010;49(3):326–36.
Zhao J, Shimizu K. Metabolic flux analysis of Escherichia coli K12 grown on 13C-labeled acetate and glucose using GC-MS and powerful flux calculation method. J Biotechnol. 2003;101(2):101–17.
Kolkman A, Olsthoorn MM, Heeremans CE, Heck AJ, Slijper M. Comparative proteome analysis of Saccharomyces cerevisiae grown in chemostat cultures limited for glucose or ethanol. Mol Cell Proteomics. 2005;4(1):1–11.
DeRisi JL, Iyer VR, Brown PO. Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997;278(5338):680–6.
Murphy JP, Stepanova E, Everley RA, Paulo JA, Gygi SP. Comprehensive temporal protein dynamics during the diauxic shift in Saccharomyces cerevisiae. Mol Cell Proteomics. 2015;14(9):2454–65.
Zampar GG, Kümmel A, Ewald J, Jol S, Niebel B, Picotti P, Aebersold R, Sauer U, Zamboni N, Heinemann M. Temporal system-level organization of the switch from glycolytic to gluconeogenic operation in yeast. Mol Syst Biol. 2013;9:651.
Kim Y, Ingram LO, Shanmugam KT. Dihydrolipoamide dehydrogenase mutation alters the NADH sensitivity of pyruvate dehydrogenase complex of Escherichia coli K-12. J Bacteriol. 2008;190(11):3851–8.
Oliveira AP, Ludwig C, Picotti P, Kogadeeva M, Aebersold R, Sauer U. Regulation of yeast central metabolism by enzyme phosphorylation. Mol Syst Biol. 2012;8:623.
Tuma DJ, Casey CA. Dangerous byproducts of alcohol breakdown–focus on adducts. Alcohol Res Health. 2003;27(4):285–90.
Caspeta L, Castillo T, Nielsen J. Modifying yeast tolerance to inhibitory conditions of ethanol production processes. Front Bioeng Biotechnol. 2015;3:184.
Chen H, Wang Z, Wang Z, Dou J, Zhou C. Improving methionine and ATP availability by MET6 and SAM2 co-expression combined with sodium citrate feeding enhanced SAM accumulation in Saccharomyces cerevisiae. World J Microbiol Biotechnol. 2016;32(4):56.
Thomas D, Surdin-Kerjan Y. Metabolism of sulfur amino acids in Saccharomyces cerevisiae. Microbiol Mol Biol. 1997;61(4):503–32.
de Kok S, Kozak BU, Pronk JT, van Maris AJ. Energy coupling in Saccharomyces cerevisiae: selected opportunities for metabolic engineering. FEMS Yeast Res. 2012;12(4):387–97.
Westfall PJ, Pitera DJ, Lenihan JR, Eng D, Woolard FX, Regentin R, Horning T, Tsuruta H, Melis DJ, Owens A, Fickes S, Diola D, Benjamin KR, Keasling JD, Leavell MD, McPhee DJ, Renninger NS, Newman JD, Paddon CJ. Production of amorphadiene in yeast, and its conversion to dihydroartemisinic acid, precursor to the antimalarial agent artemisinin. Proc Natl Acad Sci USA. 2012;109(3):E111–8.
de Jong BW, Shi S, Siewers V, Nielsen J. Improved production of fatty acid ethyl esters in Saccharomyces cerevisiae through up-regulation of the ethanol degradation pathway and expression of the heterologous phosphoketolase pathway. Microb Cell Fact. 2014;13(1):39.
Jiang GZ, Yao MD, Wang Y, Zhou L, Song TQ, Liu H, Xiao WH, Yuan YJ. Manipulation of GES and ERG20 for geraniol overproduction in Saccharomyces cerevisiae. Metab Eng. 2017;41:57–66.
Kozak BU, van Rossum HM, Niemeijer MS, van Dijk M, Benjamin K, Wu L, Daran JM, Pronk JT, van Maris AJ. Replacement of the initial steps of ethanol metabolism in Saccharomyces cerevisiae by ATP-independent acetylating acetaldehyde dehydrogenase. FEMS Yeast Res. 2016;16(2):fow006.
Costenoble R, Picotti P, Reiter L, Stallmach R, Heinemann M, Sauer U, Aebersold R. Comprehensive quantitative analysis of central carbon and amino-acid metabolism in Saccharomyces cerevisiae under multiple conditions by targeted proteomics. Mol Syst Biol. 2011;7:464.
Hackett SR, Zanotelli VR, Xu W, Goya J, Park JO, Perlman DH, Gibney PA, Botstein D, Storey JD, Rabinowitz JD. Systems-level analysis of mechanisms regulating yeast metabolic flux. Science. 2016;354(6311):aaf2786.
KH designed and carried out all experiments, and drafted the manuscript. FM revised the manuscript and further strengthened the interpretation of the data. HS supervised and coordinated the research, and revised the manuscript. All authors read and approved the final manuscript.
We thank Dr. Yoshihiro Toya for his helpful comments regarding this manuscript, Dr. Nobuyuki Okahashi for his support to perform GC–MS analysis, Dr. Shuichi Kajihata and Mr. Kousuke Maeda for their support to perform 13C-MFA, and Dr. Yoshifumi Fukui and Dr. Shinji Ozawa for their helpful suggestions.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used in the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
This work was supported in-part by Grants-in-Aid for Scientific Research (grant no. 17H06303).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hayakawa, K., Matsuda, F. & Shimizu, H. 13C-metabolic flux analysis of ethanol-assimilating Saccharomyces cerevisiae for S-adenosyl-l-methionine production. Microb Cell Fact 17, 82 (2018). https://doi.org/10.1186/s12934-018-0935-6
- Saccharomyces cerevisiae
- 13C-based metabolic flux analysis
- Ethanol metabolism
- Central carbon metabolism
- Redox balance