Dynamic metabolomics differentiates between carbon and energy starvation in recombinant Saccharomyces cerevisiae fermenting xylose
© Bergdahl et al.; licensee BioMed Central Ltd. 2012
Received: 15 November 2011
Accepted: 23 April 2012
Published: 15 May 2012
The concerted effects of changes in gene expression due to changes in the environment are ultimately reflected in the metabolome. Dynamics of metabolite concentrations under a certain condition can therefore give a description of the cellular state with a high degree of functional information. We used this potential to evaluate the metabolic status of two recombinant strains of Saccharomyces cerevisiae during anaerobic batch fermentation of a glucose/xylose mixture. Two isogenic strains were studied, differing only in the pathways used for xylose assimilation: the oxidoreductive pathway with xylose reductase (XR) and xylitol dehydrogenase (XDH) or the isomerization pathway with xylose isomerase (XI). The isogenic relationship between the two strains ascertains that the observed responses are a result of the particular xylose pathway and not due to unknown changes in regulatory systems. An increased understanding of the physiological state of these strains is important for further development of efficient pentose-utilizing strains for bioethanol production.
Using LC-MS/MS we determined the dynamics in the concentrations of intracellular metabolites in central carbon metabolism, nine amino acids, the purine nucleotides and redox cofactors. The general response to the transition from glucose to xylose was increased concentrations of amino acids and TCA-cycle intermediates, and decreased concentrations of sugar phosphates and redox cofactors. The two strains investigated had significantly different uptake rates of xylose which led to an enhanced response in the XI-strain. Despite the difference in xylose uptake rate, the adenylate energy charge remained high and stable around 0.8 in both strains. In contrast to the adenylate pool, large changes were observed in the guanylate pool.
The low uptake of xylose by the XI-strain led to several distinguished responses: depletion of key metabolites in glycolysis and NADPH, a reduced GTP/GDP ratio and accumulation of PEP and aromatic amino acids. These changes are strong indicators of carbon starvation. The XR/XDH-strain displayed few such traits. The coexistence of these traits and a stable adenylate charge indicates that xylose supplies energy to the cells but does not suppress a response similar to carbon starvation. Particular signals may play a role in the latter, of which the GTP/GMP ratio could be a candidate as it decreased significantly in both strains.
KeywordsMetabolomics Yeast metabolism Xylose fermentation Metabolic status Starvation Bioethanol
The yeast Saccharomyces cerevisiae has been the organism of choice in the food and beverage industry due to its excellent growth and fermentation capabilities under anaerobic conditions. These characteristics combined with high tolerance to low pH and high sugar and ethanol concentrations yields a production organism, which is very robust in industrial processes. Thus, S. cerevisiae has also become the preferred organism for the production of biofuels and fine chemicals . An important step towards the replacement of fossil raw materials is the efficient utilization of renewable lignocellulose. This raw material is generated in the forest and agricultural sectors, and contains carbohydrates which can be converted into fuel-grade ethanol by fermentation . The pentose sugar xylose constitutes a major fraction of the sugar monomers obtained after hydrolysis of certain lignocellulose materials [3, 4] and complete utilization of xylose is therefore necessary to obtain competitive process economics .
The rational design of industrially applicable microorganisms requires an understanding of the cellular processes that give rise to a certain phenotype. The functional information about a cellular phenotype is found in the biochemical processes that occur in response to environmental conditions and is commonly viewed as the connection between the three major ‘omes in the cell: transcriptome, proteome and metabolome [16, 17]. Both proteins and metabolites are directly involved in the cellular biochemistry and thereby closely dictate the function of the organism . The functional role of intracellular metabolites is eminent as they can influence the conversion rate by an enzyme, either as substrate, product or allosteric effector, act as signalling molecules and even influence the direction of enzymatic reactions through their thermodynamic properties . The concerted effect of changes in gene expression in response to environmental variation is thus ultimately reflected in the metabolome. The dynamic changes in metabolite concentrations under a certain perturbation can therefore provide a description of the phenotype and the cellular state with a higher degree of functional information than a snap-shot of either the transcriptome or the proteome. We used this potential to evaluate the metabolic status of two isogenic recombinant strains of S. cerevisiae harbouring either the XR/XDH or the XI pathway during anaerobic batch fermentation of a glucose/xylose mixture. A recently developed method for quantification of intracellular concentrations of metabolites by LC-MS/MS  was used to determine the dynamic response of targeted metabolites in these recombinant strains when the carbon and energy source changed from glucose to xylose. The difference in xylose uptake rate between the XI-strain and the XR/XDH-strain allows us to propose responses in metabolite concentrations expected to occur during starvation for carbon and energy. The results suggest that xylose is primarily used as an energy source as both strains maintained a high energy charge during the transition to xylose fermentation while varying signs of carbon starvation were observed.
In this study we determined the dynamics in the concentrations of metabolites in glycolysis, the pentose phosphate pathway (PPP), the TCA cycle, nine amino acids, the purine nucleotides and redox cofactors, in two recombinant strains of Saccharomyces cerevisiae. The two strains have identical genetic background but employ two different pathways for xylose utilization; strain TMB 3057  has the oxido-reductive fungal pathway consisting of XR and XDH from Pichia stipitis and strain TMB 3359  has the isomerisation pathway consisting of the XI from Piromyces sp. Samples for determination of metabolite concentrations were collected at specific time points during anaerobic batch fermentation with 20 g/L glucose and 40 g/L xylose as substrates. Additionally, for the XR/XDH-strain measured intracellular metabolite concentrations were validated by evaluating the simultaneous thermodynamic feasibility of 237 metabolic and transport reactions at each sampling point.
Validation of measured metabolite concentrations in the XR/XDH-strain
The metabolic response to a change in carbon source from glucose to xylose was determined by measuring intracellular metabolites at 11 time points during 24 hours of anaerobic batch cultivation. The first samples were collected at mid-exponential growth on glucose while the majority of samples were collected during the transition towards glucose depletion. The last samples were collected 18 hours after glucose exhaustion and represent a state when xylose is the sole carbon and energy source. Quadruplicate samples were collected at each time point and experiments were performed in duplicate for each strain. The final data set consisted of 185 samples and 44 measured metabolites.
The thermodynamic feasibility of 38 out of the 44 metabolite ranges measured in the XR/XDH-strain were validated using anNET toolbox . These 38 metabolites were quantified using external standards. The metabolites that could not be quantified, due to lack of standards, were sedoheptulose 7-phosphate (S7P), acetyl-CoA (ACCOA), reduced and oxidized glutathione (GTTred and GTTox), FAD and FMN. All metabolites evaluated, apart from 1,3-bisphosphoglycerate (BPG), glucose 1-phosphate (G1P) and glyceraldehyde 3-phosphate (GAP), were within the feasible ranges in all 11 time points (see Additional file 2 for further details). This means that the measurements of the remaining 35 metabolites are of good quality as they allowed 237 reactions to be thermodynamically feasible simultaneously.
Global analysis of the metabolite responses
Identification of key metabolites and phases during batch fermentation
Principal component analysis (PCA) was used to determine the natural grouping of the sampling points based on quantified and validated concentrations, and to identify which metabolites have the highest influence on the separation of the data. The concentrations of pyruvate were not included in this analysis due to large errors in the measurements from the XI-strain (Additional file 1: Figure S6E). The first and second principal components (PC1 and PC2, respectively) accounted for 76.2% and 10.4% of the variance in the data, respectively (Additional file 1: Figure S1). The data points were separated by sampling time in PC1 and by strain background in PC2. The grouping of the data closely reflects the experimental structure. Samples from the glucose-xylose phase cluster together, illustrating similar metabolite profiles in the two strains during glucose metabolism, which reflects the isogenicity of the strains with the only difference being the xylose pathway. This difference becomes evident when glucose is depleted. Samples collected from the xylose phase clearly form two separate clusters, one for each strain. The grouping of the data from the XI-strain also reveals that there are three distinct phases during the fermentation: a glucose phase between 14.5 h and 18.5 h, a transition phase between 18.5 h and 22 h and a xylose phase between 22 h and 38.5 h. Among the 20 most influential metabolites on PC1 and PC2 the distribution was as follows: nine amino acids, five sugar phosphates, three carboxylic acids, two nucleotides and one redox cofactor (Additional file 1: Figure S1). This distribution agrees well with the results of the cluster analysis (Figure 3). It is worth noting that several parts of metabolism are represented by these metabolites and most importantly that the adenine nucleotides are not among them.
A two-sided t-test was used to establish whether the average concentration of each quantified metabolite was different between the two strains at each sampling point. This analysis showed that the number of metabolites that were significantly different between the strains (at a significance level of α = 0.05) increased from between 16 and 23 in the mixed sugar phase to between 29 and 33 in the xylose phase (Additional file 1: Table S2), which confirms the general pattern of the PCA analysis. The majority of the quantified metabolites (81%) were significantly different in more than half of the time points, including all the metabolites in Figure S1 (Additional file 1). The remaining metabolites that showed few differences between the two strains were the nucleotides and intermediates of the TCA cycle which are in Group B of the cluster analysis (Figure 3).
Dynamic changes in sugar phosphates, precursor metabolites and amino acids
All the free amino acids measured in the current study, except glutamine (GLN-L), accumulated during the glucose and/or the transition phase (Figure 3; Additional file 1: Figure S7). The largest increase was seen in the pools of the aromatic amino acids phenylalanine (PHE-L), tryptophan (TRP-L) and tyrosine (TYR-L). The response was stronger in the XI-strain where the concentrations of PHE-L, TRP-L and TYR-L increased by 9-, 6- and 17-fold, respectively, while in the XR/XDH-strain the response was more moderate with 4-, 2- and 11-fold increases, respectively (Additional file 1: Figure S2). Concentrations of L-malate and fumarate also increased during the glucose and transition phases by ca. 2.3-fold and 1.4-fold in the XR/XDH-strain and XI-strain, respectively (Additional file 1: Figure S6 C and D). However, in the xylose phase these two carboxylic acids decreased in concentration in the XI-strain, whereas they continued to accumulate in the XR/XDH-strain, reaching around 5-fold higher concentrations compared with the glucose phase.
Dynamic changes in redox cofactors and nucleotides
For a culture of microorganisms to obtain balanced growth a careful coordination between nutrient assimilation, transformation of energy to generate ATP, redox balancing, biosynthesis and progression through the cell division cycle is required. In this study we determined the time-dependent changes of targeted metabolites in two recombinant strains of S. cerevisiae, harbouring either the XR/XDH- or the XI-pathway. This data allows us to evaluate the status of the five cellular processes required for balanced growth, from a metabolic perspective, when the cells shift from glucose to xylose fermentation.
Carbon assimilation is not the main limitation of growth on xylose
When the carbon source changed from glucose to xylose there was a significant decrease in sugar uptake rate (Figure 2). Both strains had the same specific uptake rate of glucose, but that of xylose was 60% and 94% lower in the XR/XDH-strain and the XI-strain, respectively. These results agree well with previously reported performance data for both the XR/XDH-strain  and the XI-strain . The results for the XI-strain also compare well with data from a related strain with the same genetic background, but with a different promoter controlling the expression of the XI-encoding gene . The XI-strain used in the current study does not perform as well as other yeast strains harbouring an XI but which have also undergone extensive evolutionary engineering . However, for the current investigation rationally designed strains were deliberately chosen to ensure isogenicity, except for the two initial xylose assimilating pathways, so that it can be reasonably assumed that the metabolism of the strains is controlled by the same regulatory system. The isogenic relationship between the two strains thus ensures that the observed responses are primarily a consequence of the metabolic impact of the adopted xylose assimilation pathway. When xylose is the sole carbon source mainly the high affinity glucose transporters are expressed  which have a Km ≈ 130 mM for xylose . Since the concentration of xylose during our experiments was on average double the Km the transporters would be operating at 2/3 of the maximum velocity, which corresponds to approximately 0.26 g/g CDW/h (assuming a single transporter with Vmax = 110 nmol/min/mg protein  and a cellular protein content of 40%). This value is probably an underestimation since more than one transporter is actually present, but the magnitude can be considered reasonable as it is close to the measured uptake rate in the XR/XDH-strain (0.3 g/g CDW/h). In the XI-strain on the other hand the measured uptake rate is several folds lower than the potential capacity, indicating that other factors downstream sugar transport is limiting the uptake of xylose. From this calculation it is difficult to evaluate whether transport is limiting in the XR/XDH-strain, but even though the uptake rate is 10-fold higher than in the XI-strain it is not able to grow on xylose. Attempts to improve the uptake of xylose in other XR/XDH-strains by introducing a heterologous glucose/xylose facilitator only had an effect at low xylose concentrations and the increased uptake did not improve biomass or ethanol yield from xylose . This suggests that also in the XR/XDH-strain there are other limitations in addition to sugar uptake.
A carbon starvation response is displayed in biosynthetic metabolites
In this study we measured 16 metabolites that are related to biosynthesis. Seven of these were precursor metabolites  required for biosynthesis of macromolecular building blocks (i.e. G6P, F6P, R5P, PEP, 3PG, AKG and ACCOA) (Figure 4; Figure 5) and the remaining nine were amino acids (Figure 3; Additional file 1: Figure S2 and S7). In total there are twelve essential precursor metabolites which are all produced in glycolysis, PPP and the TCA-cycle (Figure 1A). Hence, the rate of their formation is closely linked to the carbon uptake rate.
Concentrations of PPP intermediates reflect the xylose uptake rate
The measured concentrations of the metabolites in the non-oxidative PPP were markedly higher in the XR/XDH-strain compared to the XI-strain, even during growth on glucose (Figure 4A and B). This is most likely caused by the difference in xylose uptake. Between 14.5 h and 18.5 h the specific glucose uptake rate differed by less than 9% between the two strains (data not shown), however, when the uptake of xylose is also taken into account the difference in total sugar uptake rate is 45% (data not shown). The small difference with regard to glucose is consistent with the isogenic relationship between the two strains. Also the large difference in total sugar uptake rate, caused by a higher uptake of xylose by the XR/XDH-strain, agrees well with the higher concentrations of pentose sugar phosphates measured in this strain. The carbon channelled through the PPP during fermentative growth of natural S. cerevisiae is only around 15% of the glucose utilized  and the PEP generated by employing the PPP has been estimated to be less than 4% [46, 47]. The role of PPP has thus been suggested as exclusively anabolic and not catabolic . Native S. cerevisiae accumulates intracellular S7P during xylulose fermentation which was interpreted as insufficient activity of subsequent PPP enzymes . This hypothesis seemed to be confirmed when overexpression of the four genes of the non-oxidative PPP led to increased fermentation of xylulose  and xylose . Many recombinant xylose fermenting strains of S. cerevisiae, like the two strains used in the current investigation, currently carry this trait . Nevertheless, recent metabolomics studies continue to point out overexpression of the non-oxidative PPP enzymes as targets for improved pentose utilization [37, 38]. Based on the dynamics of the PPP intermediates in the current investigation with two isogenic strains, the limitation rather seems to have been shifted toward the initial xylose assimilating steps. This appears to be the case with the XI pathway in particular. One consequence of this limitation is that the XI-strain is unable to sustain a flux towards 6PG formation (Figure 4D). 6PG is formed and consumed in two irreversible reactions catalysed by 6-phosphogluconolactonase (PGL; EC 22.214.171.124) and 6-phosphogluconate dehydrogenase (GND; EC 126.96.36.199), respectively (Figure 1A). Hence, the concentration of 6PG depends on the ratio between the rates of formation and consumption. The depletion of 6PG in the XI-strain after glucose exhaustion indicates a much slower rate of 6PG formation than consumption and consequently a lower flux through the oxidative PPP. The flux through the oxidative PPP depends on the formation of G6P, which during xylose assimilation can only occur via the non-oxidative PPP and phosphoglucoisomerase (PGI; EC 188.8.131.52) (Figure 1A). The depletion of 6PG in the XI-strain is therefore most likely caused by a reduced flux through the non-oxidative PPP during xylose assimilation.
Significant decrease of hexose phosphates influences metabolite homeostasis in lower glycolysis
The concentration of the hexose phosphates in glycolysis decreased significantly in both strains during the transition from glucose to xylose fermentation (Figure 5; Additional file 1: Figure S5B). This was a general response in all glycolytic metabolites except for PEP. In bacteria, increased concentration of PEP is a generic response to carbon starvation . Elevated levels of PEP have also been reported in carbon limited or starved yeast cells under aerobic conditions [52–54] but not under anaerobic conditions . In mutant strains it was shown, using traditional biochemical methods, that the hexose phosphates in upper glycolysis play an important role for the activation of the enzymes in lower glycolysis [55, 56]. ATP generation and production of ethanol depends on a high activity of pyruvate kinase (PYK; EC 184.108.40.206) which catalyses the last reaction in glycolysis where PEP is converted to pyruvate (PYR) (Figure 1A). This enzyme requires both G6P and F6P for induction  and FBP for activation  (Figure 1B). FBP is formed from F6P and ATP in the reaction catalysed by 6-phosphofructo kinase (PFK; EC 220.127.116.11) whose activity significantly increases in the presence of fructose 2,6-bisphosphate (F26BP)  (Figure 1B). The formation of F26BP is activated through phosphorylation of 6-phosphofructo-2-kinase (PFK26; EC 18.104.22.168) by protein kinase A (PKA)  which is most active in the presence of glucose  (Figure 1B). This series of metabolic interactions could be the underlying mechanism of the stimulating effect low glucose concentration has on xylose uptake [61–63]. However, it does not explain why xylose is not taken by the engineered yeast cells in a medium with this compound in excess. The accumulation of PEP in the XI-strain suggests a reduced activity of PYK which is consistent with the depletion of FBP. This is most likely connected with the low uptake of xylose by this strain as FBP was not depleted in the XR/XDH-strain (Additional file 1: Figure S5B).
Accumulation of aromatic amino acids is a general response to carbon starvation
Increased concentrations of alanine (ALA-L), asparagine (ASN-L), glutamate (GLU-L) and the aromatic amino acids has been observed in yeast when starved for carbon under aerobic conditions . Similar results for GLU-L and the aromatic amino acids PHE-L and TRP-L have also been reported for anaerobic conditions . Thus the currently observed accumulation of amino acid pools (Figure 3) corresponds well with previously reported response to carbon starvation, both qualitatively and in many cases also quantitatively (Figure S2D). Aromatic amino acids are the most “expensive” amino acids to synthesize, in terms of ATP requirement . It is therefore unlikely that the observed accumulation is a result of continued synthesis. Instead it is more a reflection of protein degradation with subsequent accumulation of aromatic amino acids, whose aromatic ring is energetically “expensive” or unfavourable to break open. This hypothesis is supported by a measured decrease in total protein content from 0.5 ± 0.05 g protein/g CDW to 0.4 ± 0.03 g protein/g CDW in the XR/XDH-strain when the carbon source changed from glucose to xylose (data not shown).
The XR/XDH pathway has no apparent effect on the redox balance
Declining concentrations of redox cofactors have been observed in cells entering stationary phase  and during carbon starvation [37, 52]. The introduction of the XR/XDH pathway in S. cerevisiae increases the demand for NADPH and the need to re-oxidize NADH (Figure 1A). The increased demand for NADPH is met by increased expression of genes encoding enzymes which catalyse NADPH-yielding reactions, such as those of the oxidative pentose phosphate pathway [24, 66]. During xylose assimilation NADPH generation in the oxidative PPP depends on the flux from the non-oxidative PPP via F6P to G6P. The current study shows that the XR/XDH-strain, despite its increased demand for NADPH, is able to maintain a low but constant concentration throughout the xylose phase (Figure S3B). In contrast, the XI-strain is depleted of NADPH already at the beginning of the xylose phase which coincides with the depletion of 6PG after glucose exhaustion (Figure 4D). The depletion of NADPH supports the previous suggestion that the flux through the non-oxidative PPP is low in the XI-strain and it also indicates that this strain suffers from significant carbon limitation.
When the demand for oxidized NADH, i.e. NAD, is not fully met xylitol is secreted as was observed during the transition and xylose phases for the XR/XDH-strain only (Figure 2). The incapacity of both native and recombinant xylose-utilizing XR/XDH yeast strains to reconstitute the NAD/NADH balance under anaerobic conditions has frequently been suggested to be the most important bottleneck in the fermentation of xylose [13, 67] and has been the major incentive to develop XI-strains . However, the current targeted metabolite analysis recorded little significant difference between the XR/XDH-strain and the XI-strain in NADH concentration (Figure S3D) and NADH/NAD ratio (Figure 6B). Furthermore, no accumulation of NADH was observed in the XR/XDH-strain during the xylose phase. Although unexpected, it is most likely a result of mixing the cytosolic and mitochondrial pools of redox cofactors which is unavoidable with the extraction method used. If this is the case, it indicates that it is necessary to take into account the compartmentalization of metabolites when evaluating the balancing of redox cofactors in yeast. Another aspect of the redox balance is that several genes associated with respiratory metabolism become derepressed, regardless of oxygenation level, when xylose is the sole carbon source [24, 66, 68, 69]. This specifically leads to increased transcript levels of the three genes coding for malate dehydrogenase (MDH): MDH1 (mitochondrial), MDH2 (cytosolic) and MDH3 (peroxisomal). Under anaerobic conditions Mdh1p converts oxaloacetate to L-malate and is required for reoxidizing the NADH formed during amino acid synthesis inside the mitochondria  (Figure 1A). The physiological role of Mdh2p is to generate oxaloacetate for G6P-synthesis during growth on non-fermentable carbon sources and is thus repressed by glucose . The increased expression of both MDH1 and MDH2 on xylose gives the cells an enhanced NADH-oxidizing capability. This capability seems to be utilized by the XR/XDH-strain as the concentration of L-malate, the product of the MDH-catalysed reaction, accumulated intracellularly during xylose consumption (Figure 3; Additional file 1: Figure S6D). Even though this path can aid in the balancing of redox cofactors it does not appear to be sufficient as a significant amount of xylitol is still being formed (Figure 2). Instead this path might actually be detrimental for the cells, rather than beneficial, since oxaloacetate is an important precursor metabolite for synthesis of several amino acids, including methionine (Figure 1A). A diversion of flux from these reactions towards L-malate production could result in depletion of these amino acids. Unfortunately we cannot confirm this hypothesis due to lack of measurements of these amino acids and oxaloacetate. However, some support is given by the fact that yeast cells with high Mdh2p-activity only grow half as fast on glucose as the wild-type cells, thus illustrating its negative effect .
Both yeast strains are able to maintain high energy charge on xylose
The most commonly used measure to assess the energetic status of cells is the adenylate energy charge (AEC) . It is defined as the ratio between the total adenine nucleotide pool and the sum of the ATP pool and half the ADP pool: (ATP + ½ADP)/(ATP + ADP + AMP). The adenine nucleotides are balanced by adenylate kinase (AK) which catalyses the reaction ATP + AMP ↔ 2 ADP with the ADP pool serving as energy reserve to be used when ATP levels decrease. The AEC has been suggested as an important regulatory signal for controlling the energy balance by directly modulating enzyme activities [73–75]. In S. cerevisiae the AEC drops sharply from 0.8 to around 0.5 in response to glucose depletion and continues to fall during prolonged glucose starvation [36, 37]. It would therefore have been expected that the low uptake rate of xylose by the two strains (Figure 2) would lead to a decrease in the AEC and ATP concentration, especially in the XI-strain. However, the results show that there is a remarkable consistency between the adenine nucleotide pools in the two strains, clearly demonstrated by the very small confidence intervals in the energy charge (Figure 7A). This implies that the regulatory signal the AEC provides is the same during xylose fermentation as during glucose fermentation. Furthermore, it may reflect an ability of the cells to sense the surplus of external xylose even at very low uptake rates. The observed consistency could thus be the outcome of a mechanism to prevent depletion of ATP which is required to phosphorylate whatever small amount of carbon source is taken up. Therefore the cells do not respond with an energy starvation but only a carbon starvation phenotype. However, the absence of growth despite a high and stable AEC strongly indicates that it is not the required signal for continued growth on xylose.
Role of guanine nucleotides in cell division cycle control
Before a yeast cell commits to start the cell division cycle two criteria need to be fulfilled: accumulation of cell mass to reach a critical size  and establishment of a critical rate of protein synthesis . This occurs during the G1 phase of the cycle in which the cell evaluates the nutritional condition of the environment. The guanine nucleotides have been implicated as important regulatory factors in both these processes [78–80]. Rudoni et al. reported a value for the GTP/GDP ratio in exponentially growing yeast cells of ca. 4 which is reduced to below one in stationary phase and during glucose starvation. Furthermore, they found a correlation between the cytosolic GTP/GDP ratio and the level of GTP-bound Ras2, which is the active form of this protein. In this form it can bind to adenylate cyclase and trigger the production of cAMP  (Figure 1B). The active Ras/cAMP/PKA pathway initiates the biosynthetic machinery by increasing the expression of genes involved in ribosome biogenesis  and is thereby required for the cell to reach the size threshold and start the cell cycle . The synthesis of proteins is highly dependent on GTP availability, not only in the elongation step but also in the initiation of translation . The two initiation factors eIF2α and eIF2B have a higher affinity for GDP than GTP and it is therefore not surprising that an imbalance in the GTP/GDP ratio also can affect the rate of protein synthesis [79, 80].
The value of the GTP/GDP ratio during exponential growth on glucose obtained in the current study (Figure 7D) agrees well with the results presented by Rudoni et al.. The dynamic change in the XI-strain also corresponds with the reported response to carbon exhaustion. The very low GTP/GDP ratio in this strain on xylose could result in low Ras/cAMP activity and rate of protein synthesis leading to an arrest in the G0 phase of the cell cycle. This would support the idea that xylose is not sensed as a fermentable carbon source by the yeast cells [24, 66, 69] and does not activate the biosynthetic machinery efficiently. However, the GTP/GDP ratio did not change significantly in the XR/XDH-strain during the transition to xylose fermentation, suggesting that the GTP/GDP ratio rather responds to the rate of carbon supply than to what type of carbon source is actually taken up. Hence, other constituents of the nutrient sensing pathways must be involved, e.g. Snf1 and the Tor kinases which can control the initiation of translation by modulating the phosphorylation of eIF2α .
Other indicators of metabolic status
In contrast to the GTP/GDP ratio both the guanylate energy charge (GEC) and the GTP/GMP ratio decreased significantly in the two strains when shifting to xylose fermentation (Figure 7B and C). The fast and distinct response in the GEC indicates that this measure could serve as a more sensitive indicator of metabolic status than the AEC. However, yeast lacks an enzyme equivalent to adenylate kinase that can catalyse the conversion of two molecules of GDP to GTP and GMP. Thus, the GEC can be calculated analogous to the AEC, but it is not related to any biological process. Even so, the pools of the adenine and guanine nucleotides are directly interconnected as GTP is used as a phosphate donor in the de novo synthesis of ATP and vice versa . This could explain the similar values of the AEC and the GEC (0.82 and 0.86, respectively) during the glucose phase (Figure 7A and B). A 10-fold decrease in the GTP/GMP ratio, similar to the one reported here, as a response to a change from glucose to xylose fermentation has also been measured in a recent study . This highlights the sensitivity of this ratio and could therefore be a better indicator of metabolic status as it is directly related to metabolic processes. However, there have been no studies on how and to what extent the GTP/GMP ratio affects or controls metabolic processes. Even though the role of the guanine nucleotides in yeast metabolism has been less explored than that of the adenine nucleotides, it is clear that the biochemical processes in which these nucleotides participate are required for both growth and survival [86, 87]. Recent studies have shown that an imbalance in the guanine nucleotide pool can have detrimental effects on nitrogen metabolism, protein synthesis, viability and proliferation of yeast cells [78, 79, 88, 89]. Still, more systematic investigations of the response in the guanine nucleotide pools under different conditions are needed to understand the role they play in the physiology of yeast.
In this study we determined the time-dependent changes in a selected part of the metabolome of two recombinant strains of S. cerevisiae, harbouring either the XR/XDH or the XI pathway, in response to a transition from glucose to xylose fermentation. A fundamental difference between the two strains was the very low uptake rate of xylose by the XI-strain, which gave rise to larger changes in metabolite concentrations than in the XR/XDH-strain. However, even though the XR/XDH-strain had a 10-fold higher uptake rate of xylose it did not continue to grow after glucose depletion. This indicates that factors other than sugar transport limit the growth on xylose. Several distinct responses in the XI-strain were not observed in the XR/XDH-strain: depletion of NADPH and several metabolites in glycolysis, accumulation of PEP and a reduced GTP/GDP ratio. These changes in the metabolome typically direct towards carbon starvation, which appears to be the general response of the XI-strain to xylose. Still, both strains maintained a high and stable energy charge during the transition to xylose, indicating that neither strain is energy limited. In fact, the only response observed in both strains was the accumulation of aromatic amino acids. We suggest that this is related to a general decrease in protein concentration. We also suggest that the redox imbalance between XR and XDH, combined with the poor ability of xylose to repress genes associated with respiratory pathways, give rise to reactions unbeneficial for xylose fermentation.
In light of the stable energy charge, xylose fermentation seems to present a unique condition for S. cerevisiae where the energy metabolism is uncoupled from carbon metabolism. The observed discrimination between energy and carbon starvation during xylose fermentation has not been reported previously, but is now made visible through metabolomic analysis of isogenic strains. However, the maintenance of a high energy charge does evidently not correlate with the general metabolic status of the cells; hence, other internal signals could be better indicators. Potential candidates are the guanylate energy charge and the GTP/GMP ratio as they decreased significantly in both strains. The importance of these relationships as indicators of metabolic status merits further investigation as the guanine nucleotides have been shown to play an important regulatory role in cell biology and metabolism. Dynamic metabolomics is a suitable method for exposing such underlying metabolic interactions and to unravel the behaviour of rationally designed cell factories and wild type microorganisms.
Strains and cultivation conditions
In this study, two recombinant strains of Saccharomyces cerevisiae, with the same genetic background but with two different pathways for xylose utilization, were used, i.e. TMB 3057 employing the XR/XDH pathway  and TMB 3359 utilizing a xylose isomerase from Piromyces sp . Liquid cultures were grown in yeast nitrogen base medium (YNB) without amino acids (1.7 g/L; Difco, Becton, Dickinson and Company, Sparks, MD, USA) with 20 g/L glucose and 5 g/L ammonium sulphate. In anaerobic fermentation, double concentration of YNB (2x YNB) with 20 g/L glucose and 40 g/L xylose was used, supplemented with 5 g/L ammonium sulphate, 400 mg/L Tween 80 and 10 mg/L ergosterol.
Cells were pre-grown in YNB medium until late exponential phase and inoculated into the bioreactor at a concentration of 0.04 g CDW/L. Anaerobic fermentation was performed in a 3.1 L Benchtop Fermenter Type KLF2000 bioreactor (Bioengineering AG, Switzerland) with a working volume of 2 litres. The temperature was 30°C, stirring was set at 900 rpm and pH was controlled at 5.5 with 3 M KOH and 1.5 M H2SO4. Anaerobic conditions were maintained by sparging the culture broth with 0.5 L/min nitrogen gas (Grade 5.0, containing less than 3 ppm oxygen; PanGas AG, Switzerland) controlled by a gas flow controller (Get red-y, Vögtlin Intruments AG, Switzerland). The condenser was connected to a Thermostatic Circulator 2219 (LKB Produkter AB, Sweden) and cooled to 4°C. Foaming was prevented by adding three to four drops of sterile 50% PEG2000. Fermentation experiments were performed in duplicate.
Quantification of extracellular substrates and products
Concentrations of substrates and products in the culture broth were determined by HPLC (Waters, USA) with Aminex HPX-87 H ion exchange column (Bio-Rad, USA) and refractive index detector (RID-6A, Shimadzu, Japan). Mobile phase was 5 mM H2SO4, the temperature 45°C and the flow rate 0.6 mL/min. Apart from glucose, xylose, xylitol, glycerol and ethanol, no other compounds were detected in the analysis. Stated concentrations of ethanol are not corrected for losses due to evaporation.
Optical density was measured in triplicate at each time point at 600 nm. Cell dry weight was determined in triplicate by filtering 10 mL of culture through a dried and pre-weighed nitrocellulose filter with a 0.45 μm pore size. The filters were washed with an equal amount of water and dried at 80°C overnight. The filters were allowed to equilibrate to room temperature before weighed again.
Quenching, extraction of intracellular metabolites and sample preparation
The quenching procedure of cell samples was based on the method developed by de Koning et al.. Quenching was achieved by quickly adding 1 mL of cell culture to 4 mL of cold, buffered methanol (60% methanol, 10 mM ammonium acetate, pH 7.5, -40°C). The cells were cooled at −40°C in an ethylene glycol bath (FP40-MC, Julabo, Germany) for maximum 30 seconds and then collected by centrifugation at 5000 rpm and −9°C for 5 min. The supernatant was decanted and any residual liquid was removed with a pipette. The pellet was then immediately frozen in liquid nitrogen and stored at −80°C until extraction. Quadruplicate samples were quenched at each time point.
Intracellular metabolites were extracted using the boiling ethanol procedure . Frozen cell pellets were briefly thawed at −40°C before extraction. To extract the metabolites, 1 mL of boiling, buffered ethanol (75% ethanol, 10 mM ammonium acetate, pH 7.5, 85°C) was added to the pellet, directly followed by addition of 100 μL of an internal standard (metabolite extract of yeast grown on 13C-labeled glucose) and vortexed quickly. The mixture was incubated 3 min at 85°C with vigorous mixing every 30 seconds. The solution was then shortly cooled back to −40°C and cell debris was spun down by centrifugation at 5000 rpm and −9°C for 3 min. The supernatant was collected in a 1.5 mL tube, frozen in liquid nitrogen and stored at −80°C until further use.
Extracted metabolites were dried by evaporating the solvent at low pressure, 0.01-0.1 mbar, using a Christ Alpha 2–4 LD Plus rotational vacuum concentrator (Adolf Kühner AG, Switzerland) connected to a Chemistry Hybrid Pump RC 6 (Vacuubrand GmbH, Germany) and a Christ AVC 2–33 centrifuge (Adolf Kühner AG, Switzerland). Dried samples were stored at −80°C until further use.
Metabolite analysis by LC-MS/MS
Separation and detection of compounds was achieved on a Waters Acquity UPLC (Waters Corporation, Milford, MA, USA) using a Waters Acquity T3 end-capped reverse phase column with dimensions 150 mm × 2.1 mm × 1.8 mm (Waters Corporation) coupled to a Thermo TSQ Quantum Ultra triple quadrupole mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA) as previously described .
Data normalization, calibration and quantification
Detection and integration of chromatographic peaks was performed using Xcalibur software version 2.0.7 (Thermo Fisher Scientific Inc, USA). The automatic integrations were manually curated to obtain consistent peak boundaries and the results were exported to Matlab R2008a (The MathWorks, USA) where subsequent data processing was performed.
where X ijk and are the measured and normalized signals of analyte i in sample k of time point j, respectively, Z ijk is the signal of the 13C-labeled analyte i and is a robust mean of the internal standard signals of analyte i in all samples. The robust mean was calculated according to an iterative procedure as outlined in Miller & Miller  p. 174].
Seven external standards with known concentrations ranging from 0.14 μM to 100 μM were analyzed together with each batch of samples to be quantified. An initial robust linear regression with a logistic estimator was used to identify possible outliers in the calibration data. Data points with a z-score above three were regarded as outliers and removed from the data set . An ordinary least squares (OLS) regression was then performed with the remaining data points to obtain unbiased model and regression parameters . All calibration curves had at least five data points in the final OLS regression and the adjusted coefficients of determination (adjusted R2) were between 0.90 and 1, with the majority around 0.99. A term for the intercept was included in the equation only if the root mean squared error (RMSE) was lower than the RMSE obtained without including the intercept and if it was statistically significant according to a t-test at the 0.05 probability level.
The intracellular concentration corresponding to an estimated sample concentration was calculated from the sample volume and the measured cell dry weight. All intracellular concentrations stated in this paper are mean values of concentrations from duplicate biological experiments and errors are given as a 95% confidence interval of the calculated mean. Samples collected at a particular time point were analyzed for outliers before calculating the mean. This was done by calculating the ratio between |c ijk – median(c ij·)| and the median absolute deviation (MAD): median(|c ijk – median(c ij·)|). Measurements with a ratio above four were regarded as outliers and were removed. The total number of replicates at each time point was between four and ten, but the majority of the means were calculated from seven replicates.
Principal component analysis (PCA) is a statistical method commonly used to analyse data from metabolomics studies . PCA is an unsupervised method used to reduce the dimension of the data by calculating new components which are linear combinations of the original values . The first component is constructed so that it accounts for the highest fraction of the total variance in the data; the second component is constructed to be orthogonal (uncorrelated) to the first component and account for the highest fraction of the variance in this direction. Subsequent components are constructed similarly resulting in a number of components that are fewer than the original number of variables (dimension reduction) and uncorrelated and thus give a unique description of each sample in the new dimension. PCA is classified as an unsupervised method because no prior grouping of the data is required. The result can therefore give an insight into the natural grouping of the data based on the maximisation of the variance criterion. PCA was performed in MATLAB R2010b (Mathworks, USA) using the princomp function on quantified concentrations of 37 metabolites from 185 samples. After removing outliers according to the above mentioned criterion the data matrix contained 5.6% missing data points. If only one data point was missing in a certain sampling point and strain it was replaced by the average value of the measured concentrations. If more than one data point was missing, they were replaced by randomly generated numbers from a t-distribution with the number of degrees of freedom equal to the number of samples at that sampling point minus one. The random numbers were converted to concentrations by multiplication with the standard deviation and addition of the mean of the known concentrations. Each metabolite in the imputed data matrix was normalized to the median concentration measured in both strains at the first sampling point and then converted to log2 scale.
Hierarchical clustering was performed using the clustergram function in MATLAB R2010b with quantified concentrations and normalized chromatographic peak areas of 38 and six metabolites, respectively. The data for the XR/XDH-strain and the XI-strain contained 95 and 90 samples, respectively. Each set of samples was normalized to the corresponding mean concentration or peak area at the first sampling point before converted to log2 scale. Euclidean distance was used to calculate the distance between the metabolite profiles and average linkage was used to generate the dendrogram. Data used for the PCA and the hierarchical clustering are given as supplementary information in Additional file 3.
The thermodynamic analysis of measured metabolite concentrations in the XR/XDH-strain was performed using anNET version 1.1.06 . Detailed information regarding the calculations and discussion of the results is given in Additional file 2. All the data used in the calculations are given as supplementary information in Additional file 4.
2PG: 2-phosphoglycerate; 3PG: 3-phosphoglycerate; 6PG: 6- phosphogluconate; 6pgl: 6-phosphogluconolactonate; Acald: Acetaldehyde; ACCOA: Acetyl-Coenzyme A; ADP: Adenosine 5'-diphosphate; AKG: a-ketoglutarate; ALA-L: Alanine; AMP: Adenosine 5'-monophosphate; ARG-L: Arginine; ASN-L: Asparagine; ASP-L: Aspartate; ATP: Adenosine 5'-triphosphate; BPG, Bpg: 1,3-bisphosphoglycerate; cAMP: Cyclic AMP; CIT: Citrate; CoA: Coenzyme A; Cys: Cysteine; DHAP: Dihydroxy acetone phosphate; E4p: Erythrose 4-phosphate; F26BP: Fructose 2,6-bisphosphate; F6P: Fructose 6-phosphate; FAD: Flavin adenine dinucleotide; FBP: Fructose 1,6-bisphosphate; FMN: Riboflavin 5'-monophosphate; FUM: Fumarate; G1P: Glucose 1-phosphate; G6P: Glucose 6-phosphate; GAP, Gap: Glyceraldehyde 3-phosphate; GDP: Guanosine 5'-diphosphate; GLN-L: Glutamine; GLU-L: Glutamate; Gly: Glycine; GLYC3P: Glycerol 3-phosphate; GMP: Guanosine 5'-monophosphate; GTP: Guanosine 5'-triphosphate; GTTox: Oxidized glutathione; GTTred: Reduced glutathione; His: Histidine; ICIT: Isocitrate; Ile: Isoleucine; Leu: Leucine; L-MAL: L-malate; Lys: lysine; Met: Methionine; NAD: Oxidized nicotinamide adenine dinucleotide; NADH: Reduced nicotinamide adenine dinucleotide; NADP: Oxidized nicotinamide adenine dinucleotide phosphate; NADPH: Reduced nicotinamide adenine dinucleotide phosphate; Oaa: Oxaloacetate; PEP: Phoshoenolpyruvate; PHE-L: Phenylalanine; Pro: Proline; PYR, pyr: Pyruvate; R5P: Ribose 5-phosphate; RU5P: Ribulose 5-phosphate; S7P: Sedoheptulose 7-phosphate; Ser: Serine; SUCC: Succinate; Succoa: Succinyl-Coenzyme A; T6P: Trehalose 6-phosphate; Thr: Threonine; TRP-L: Tryptophane; TYR-L: Tyrosine; Val: Valine; xPG: Combined pool of 2PG and 3PG; XU5P: Xylulose 5-phosphate; Xylu: Xylulose;
AEC, Adenylate energy charge; GEC, Guanylate energy charge; PCA, Principal component analysis; PPP, Pentose phosphate pathway; TCA cycle, Tricarboxylic acid cycle; 2PG, 2-phosphoglycerate; 3PG, 3-phosphoglycerate; 6PG, 6-phosphogluconate; 6pgl, 6-phosphogluconolactonate; Acald, Acetaldehyde; ACCOA, Acetyl-Coenzyme A; ADP, Adenosine 5'-diphosphate; AKG, α-ketoglutarate; ALA-L, Alanine; AMP, Adenosine 5'-monophosphate; ARG-L, Arginine; ASN-L, Asparagine; ASP-L, Aspartate; ATP, Adenosine 5'-triphosphate; BPG, Bpg, 1,3-bisphosphoglycerate; cAMP, Cyclic AMP; CIT, Citrate; CoA, Coenzyme A; Cys, Cysteine; DHAP, Dihydroxy acetone phosphate; E4p, Erythrose 4-phosphate; F26BP, Fructose 2,6-bisphosphate; F6P, Fructose 6-phosphate; FAD, Flavin adenine dinucleotide; FBP, Fructose 1,6-bisphosphate; FMN, Riboflavin 5'-monophosphate; FUM, Fumarate; G1P, Glucose 1-phosphate; G6P, Glucose 6-phosphate; GAP, Gap, Glyceraldehyde 3-phosphate; GDP, Guanosine 5'-diphosphate; GLN-L, Glutamine; GLU-L, Glutamate; Gly, Glycine; GLYC3P, Glycerol 3-phosphate; GMP, Guanosine 5'-monophosphate; GTP, Guanosine 5'-triphosphate; GTTox, Oxidized glutathione; GTTred, Reduced glutathione; His, Histidine; ICIT, Isocitrate; Ile, Isoleucine; Leu, Leucine; L-MAL, L-malate; Lys, lysine; Met, Methionine; NAD, Oxidized nicotinamide adenine dinucleotide; NADH, Reduced nicotinamide adenine dinucleotide; NADP, Oxidized nicotinamide adenine dinucleotide phosphate; NADPH, Reduced nicotinamide adenine dinucleotide phosphate; Oaa, Oxaloacetate; PEP, Phoshoenolpyruvate; PHE-L, Phenylalanine; Pro, Proline; PYR, pyr, Pyruvate; R5P, Ribose 5-phosphate; RU5P, Ribulose 5-phosphate; S7P, Sedoheptulose 7-phosphate; Ser, Serine; SUCC, Succinate; Succoa, Succinyl-Coenzyme A; T6P, Trehalose 6-phosphate; Thr, Threonine; TRP-L, Tryptophane; TYR-L, Tyrosine; Val, Valine; xPG, Combined pool of 2PG and 3PG; XU5P, Xylulose 5-phosphate; Xylu, Xylulose; ACOm, Aconitase EC 22.214.171.124; ACS, Acetyl-CoA synthetase EC 126.96.36.199; ADH, Alcohol dehydrogenase EC 188.8.131.52; AK, Adenylate kinase EC 184.108.40.206; AKGDm, Alpha-ketoglutarate dehydrogenase EC 220.127.116.11; ALDH, Aldehyde dehydrogenase EC 18.104.22.168; CSm, Citrate synthase EC 22.214.171.124; CYR1, Adenylate cyclase EC 126.96.36.199; ENO, Enolase EC 188.8.131.52; FBA, Fructose 1,6-bisphosphate aldolase EC 184.108.40.206; FBP1, Fructose-1,6-bisphosphatase EC 220.127.116.11; FBP26, Fructose-2,6-bisphosphatase EC 18.104.22.168; FRD, FRDm, Fumarate reductase EC 22.214.171.124; FUMASE, FUMASEm, Fumarase EC 126.96.36.199; G6PDH, Glucose-6-phosphate dehydrogenase EC 188.8.131.52; GAPD, Glyceraldehyde-3-phosphate dehydrogenase EC 184.108.40.206; GND, 6-phosphogluconate dehydrogenase EC 220.127.116.11; GPD, NAD-dependent glycerol 3-phosphate dehydrogenase EC 18.104.22.168; GPP, DL-glycerol-3-phosphatase EC 22.214.171.124; HXK, Hexokinase EC 126.96.36.199; IDHm, Isocitrate dehydrogenase EC 188.8.131.52; MDH2, MDH1m, Malate dehydrogenase EC 184.108.40.206; NTH, Trehalase, EC 220.127.116.11; PC, Pyruvate carboxylase EC 18.104.22.168; PDC, Pyruvate decarboxylase EC 22.214.171.124; PDHm, Pyruvate dehydrogenase complex EC 126.96.36.199; PFK, 6-phosphofructo kinase EC 188.8.131.52; PFK26, 6-phosphofructo-2-kinase EC 184.108.40.206; PGAM, Phosphoglycerate mutase EC 220.127.116.11; PGI, Phosphoglucose isomerase EC 18.104.22.168; PGK, 3-phosphoglycerate kinase EC 22.214.171.124; PGL, 6-phosphogluconolactonase EC 126.96.36.199; PKA, cAMP-dependent protein kinase EC 188.8.131.52; PYK, Pyruvate kinase EC 184.108.40.206; RKI, Ribose-5-phosphate ketol-isomerase EC 220.127.116.11; RPE, D-ribulose-5-phosphate 3-epimerase EC 18.104.22.168; SUCOASm, Succinate-CoA ligase EC 22.214.171.124; TAL, Transaldolase EC 126.96.36.199; TKL1, TKL2, Transketolase EC 188.8.131.52; TPI, Triose phosphate isomerase EC 184.108.40.206; TPS1, Trehalose-6-phosphate synthase EC 220.127.116.11; TPS2, Trehalose-6-phosphate phosphatase EC 18.104.22.168; XDH, Xylitol dehydrogenase EC 22.214.171.124; XI, Xylose isomerase EC 126.96.36.199; XR, Xylose reductase EC 188.8.131.527.
This work was supported by funding from the EU project NILE within the 6th framework programme (Contract N° 019882). BB also received financial support from the Swedish Energy Agency. We wish to thank Professor Marie-Francoise Gorwa-Grauslund for critically reading the manuscript.
- Nevoigt E: Progress in metabolic engineering ofSaccharomyces cerevisiae. Microbiol Mol Biol Rev 2008, 72:379–412.View Article
- Lin Y, Tanaka S: Ethanol fermentation from biomass resources: Current state and prospects. Appl Microbiol Biotechnol 2006, 69:627–642.View Article
- Wiselogel AE, Agblevor FA, Johnson DK, Deutch S, Fennell JA, Sanderson MA: Compositional changes during storage of large round switchgrass bales. Bioresour Technol 1996, 56:103–109.View Article
- Perez DD, Guillemain A, Berthelot A, N'Guyen-The N, De Morogues F, Gomes C: Evaluation of forestry biomass quality for the production of second-generation biofuels. Cell Chem Technol 2010, 44:1–14.
- Sassner P, Galbe M, Zacchi G: Techno-economic evaluation of bioethanol production from three different lignocellulosic materials. Biomass Bioenerg 2008, 32:422–430.View Article
- Van Vleet JH, Jeffries TW: Yeast metabolic engineering for hemicellulosic ethanol production. Curr Opin Biotechnol 2009, 20:300–306.View Article
- Matsushika A, Inoue H, Kodaki T, Sawayama S: Ethanol production from xylose in engineeredSaccharomyces cerevisiaestrains: Current state and perspectives. Appl Microbiol Biotechnol 2009, 84:37–53.View Article
- Hahn-Hagerdal B, Karhumaa K, Fonseca C, Spencer-Martins I, Gorwa-Grauslund MF: Towards industrial pentose-fermenting yeast strains. Appl Microbiol Biotechnol 2007, 74:937–953.View Article
- Kern A, Tilley E, Hunter IS, Legisa M, Glieder A: Engineering primary metabolic pathways of industrial micro-organisms. J Biotechnol 2007, 129:6–29.View Article
- Karhumaa K, Sanchez RG, Hahn-Hagerdal B, Gorwa-Grauslund MF: Comparison of the xylose reductase-xylitol dehydrogenase and the xylose isomerase pathways for xylose fermentation by recombinantSaccharomyces cerevisiae. Microb Cell Fact 2007, 6:5.View Article
- Bengtsson O: Genetic traits beneficial for xylose utilization by recombinantSaccharomyces cerevisiae.Doctoral Thesis. Lund University, Division of Applied Microbiology 2008.
- Bruinenberg PM, de Bot PHM, van Dijken JP, Scheffers WA: NADH-linked aldose reductase - the key to anaerobic alcoholic fermentation of xylose by yeasts. Appl Microbiol Biotechnol 1984, 19:256–260.View Article
- Bruinenberg PM, de Bot PHM, van Dijken JP, Scheffers WA: The role of redox balances in the anaerobic fermentation of xylose by yeasts. Appl Microbiol Biotechnol 1983, 18:287–292.View Article
- Kotter P, Ciriacy M: Xylose fermentation bySaccharomyces cerevisiae. Appl Microbiol Biotechnol 1993, 38:776–783.View Article
- van Maris AJA, Winkler AA, Kuyper M, de Laat W, van Dijken JP, Pronk JT: Development of efficient xylose fermentation inSaccharomyces cerevisiae: Xylose lsomerase as a key component. In Biofuels. Edited by: Olsson L. Springer-Verlag Berlin, Berlin; 2007:179–204. [Scheper T (Series Editor): Advances in Biochemical Engineering/Biotechnology, vol 108.]View Article
- Smedsgaard J, Nielsen J: Metabolite profiling of fungi and yeast: from phenotype to metabolome by MS and informatics. J Exp Bot 2005, 56:273–286.View Article
- Delneri D, Brancia FL, Oliver SG: Towards a truly integrative biology through the functional genomics of yeast. Curr Opin Biotechnol 2001, 12:87–91.View Article
- Heinemann M, Sauer U: Systems biology of microbial metabolism. Curr Opin Microbiol 2010, 13:337–343.View Article
- Gerosa L, Sauer U: Regulation and control of metabolic fluxes in microbes. Curr Opin Biotechnol 2011, 22:566–575.View Article
- Buescher JM, Moco S, Sauer U, Zamboni N: Ultrahigh performance liquid chromatography-tandem mass spectrometry method for fast and robust quantification of anionic and aromatic metabolites. Anal Chem 2010, 82:4403–4412.View Article
- Karhumaa K, Hahn-Hägerdal B, Gorwa-Grauslund MF: Investigation of limiting metabolic steps in the utilization of xylose by recombinantSaccharomyces cerevisiaeusing metabolic engineering. Yeast 2005, 22:359–368.View Article
- Parachin NS, Bergdahl B, van Niel EWJ, Gorwa-Grauslund MF: Kinetic modeling reveals current limitations in the production of ethanol from xylose by recombinantSaccharomyces cerevisiae. Metab Eng 2011, 13:508–517.View Article
- Zamboni N, Kummel A, Heinemann M, anNET: A tool for network-embedded thermodynamic analysis of quantitative metabolome data. BMC Bioinformatics 2008, 9:199.View Article
- Runquist D, Hahn-Hagerdal B, Bettiga M: Increased expression of the oxidative pentose phosphate pathway and gluconeogenesis in anaerobically growing xylose-utilizingSaccharomyces cerevisiae. Microb Cell Fact 2009, 8:49.View Article
- de Koning W, van Dam K: A method for the determination of changes of glycolytic metabolites in yeast on a subsecond time scale using extraction at neutral pH. Anal Biochem 1992, 204:118–123.View Article
- Zamboni N, Kummel A, Heinemann M: anNET: A tool for networkembedded thermodynamic analysis of quantitative metabolome data. BMC Bioinformatics 2008, 9:199.View Article
- Mo M, Palsson B, Herrgard M: Connecting extracellular metabolomic measurements to intracellular flux states in yeast. BMC Syst Biol 2009, 3:37.View Article
- Becker SA, Feist AM, Mo ML, Hannum G, Palsson BO, Herrgard MJ: Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox. Nat Protocols 2007, 2:727–738.View Article
- Wahlbom CF, Eliasson A, Hahn-Hägerdal B: Intracellular fluxes in a recombinant xylose-utilizingSaccharomyces cerevisiaecultivated anaerobically at different dilution rates and feed concentrations. Biotechnol Bioeng 2001, 72:289–296.View Article
- Snitkin E, Dudley A, Janse D, Wong K, Church G, Segre D: Model-driven analysis of experimentally determined growth phenotypes for 465 yeast gene deletion mutants under 16 different conditions. Genome Biol 2008, 9:R140.View Article
- Jankowski MD, Henry CS, Broadbelt LJ, Hatzimanikatis V: Group contribution method for thermodynamic analysis of complex metabolic networks. Biophys J 2008, 95:1487–1499.View Article
- Klimacek M, Krahulec S, Sauer U, Nidetzky B: Limitations in xylosefermentingSaccharomyces cerevisiae, made evident through comprehensive metabolite profiling and thermodynamic analysis. Appl Environ Microbiol 2010, 76:7566–7574.View Article
- Cipollina C, ten Pierick A, Canelas AB, Seifar RM, van Maris AJA, van Dam JC, Heijnen JJ: A comprehensive method for the quantification of the nonoxidative pentose phosphate pathway intermediates inSaccharomyces cerevisiaeby GC-IDMS. J Chromatogr B 2009, 877:3231–3236.View Article
- Hynne F, Dano S, Sorensen PG: Full-scale model of glycolysis inSaccharomyces cerevisiae. Biophys Chem 2001, 94:121–163.View Article
- Altintas MM, Eddy CK, Zhang M, McMillan JD, Kompala DS: Kinetic modeling to optimize pentose fermentation inZymomonas mobilis. Biotechnol Bioeng 2006, 94:273–295.View Article
- Ball WJ, Atkinson DE: Adenylate energy charge inSaccharomyces cerevisiaeduring starvation. J Bacteriol 1975, 121:975–982.
- Klimacek M, Krahulec S, Sauer U, Nidetzky B: Limitations in xylose-fermentingSaccharomyces cerevisiae, made evident through comprehensive metabolite profiling and thermodynamic analysis. Appl Environ Microbiol 2010, 76:7566–7574.View Article
- Wisselink HW, Cipollina C, Oud B, Crimi B, Heijnen JJ, Pronk JT, van Maris AJA: Metabolome, transcriptome and metabolic flux analysis of arabinose fermentation by engineeredSaccharomyces cerevisiae. Metab Eng 2010, 12:537–551.View Article
- Karhumaa K, Fromanger R, Hahn-Hagerdal B, Gorwa-Grauslund MF: High activity of xylose reductase and xylitol dehydrogenase improves xylose fermentation by recombinantSaccharomyces cerevisiae. Appl Microbiol Biotechnol 2007, 73:1039–1046.View Article
- Kuyper M, Toirkens MJ, Diderich JA, Winkler AA, van Dijken JP, Pronk JT: Evolutionary engineering of mixed-sugar utilization by a xylose-fermenting Saccharomyces cerevisiae strain. Fems Yeast Research 2005, 5:925–934.View Article
- Sedlak M, Ho NWY: Characterization of the effectiveness of hexose transporters for transporting xylose during glucose and xylose co-fermentation by a recombinant Saccharomyces yeast. Yeast 2004, 21:671–684.View Article
- Saloheimo A, Rauta J, Stasyk OV, Sibirny AA, Penttila M, Ruohonen L: Xylose transport studies with xylose-utilizing Saccharomyces cerevisiae strains expressing heterologous and homologous permeases. Appl Microbiol Biotechnol 2007, 74:1041–1052.View Article
- Fonseca C, Olofsson K, Ferreira C, Runquist D, Fonseca LL, Hahn-Hagerdal B, Liden G: The glucose/xylose facilitator Gxf1 fromCandida intermediaexpressed in a xylose-fermenting industrial strain ofSaccharomyces cerevisiaeincreases xylose uptake in SSCF of wheat straw. Enzyme Microb Technol 2011, 48:518–525.View Article
- Stephanopoulos G, Aristidou A, Nielsen J: Metabolic Engineering: Principles and Methodologies. Elsevier Academic Press, San Diego, CA, USA; 1998.
- Frick O, Wittmann C: Characterization of the metabolic shift between oxidative and fermentative growth inSaccharomyces cerevisiaeby comparative C-13 flux analysis. Microb Cell Fact 2005, 4:30.View Article
- Gancedo JM, Lagunas R: Contribution of the pentose-phosphate pathway to glucose metabolism inSaccharomyces cerevisiae: A critical analysis on the use of labelled glucose. Plant Science Letters 1973, 1:193–200.View Article
- Maaheimo H, Fiaux J, Çakar ZP, Bailey JE, Sauer U, Szyperski T: Central carbon metabolism ofSaccharomyces cerevisiaeexplored by biosynthetic fractional13C labeling of common amino acids. Eur J Biochem 2001, 268:2464–2479.View Article
- Senac T, Hahn-Hägerdal B: Intermediary metabolite concentrations in xylulose-fermenting and glucose-fermentingSaccharomyces cerevisiaecells. Appl Environ Microbiol 1990, 56:120–126.
- Johansson B, Hahn-Hagerdal B: The non-oxidative pentose phosphate pathway controls the fermentation rate of xylulose but not of xylose inSaccharomyces cerevisiaeTMB3001. Fems Yeast Research 2002, 2:277–282.
- Chu BCH, Lee H: Genetic improvement ofSaccharomyces cerevisiaefor xylose fermentation. Biotechnol Adv 2007, 25:425–441.View Article
- Evans PR, Farrants GW, Hudson PJ: Phosphofructokinase - structure and control. Philos Trans R Soc Lond B Biol Sci 1981, 293:53–62.View Article
- Brauer MJ, Yuan J, Bennett BD, Lu W, Kimball E, Botstein D, Rabinowitz JD: Conservation of the metabolomic response to starvation across two divergent microbes. Proc Natl Acad Sci USA 2006, 103:19302–19307.View Article
- Walther T, Novo M, Roessger K, Letisse F, Loret M-O, Portais J-C, Francois J-M: Control of ATP homeostasis during the respiro-fermentative transition in yeast. Mol Syst Biol 2010, 6:344.View Article
- Kresnowati MTAP, van Winden WA, Almering MJH, ten Pierick A, Ras C, Knijnenburg TA, Daran-Lapujade P, Pronk JT, Heijnen JJ, Daran JM: When transcriptome meets metabolome: fast cellular responses of yeast to sudden relief of glucose limitation. Mol Syst Biol 2006, 2:49.View Article
- Boles E, Zimmermann FK, Heinisch J: Different signals control the activation of glycolysis in the yeastSaccharomyces cerevisiae. Yeast 1993, 9:761–770.View Article
- Mueller S, Boles E, May M, Zimmermann FK: Different internal metabolites trigger the induction of glycolytic gene expression in Saccharomyces cerevisiae. J Bacteriol 1995, 177:4517–4519.
- Barwell CJ, Woodward B, Brunt RV: Regulation of pyruvate kinase by fructose 1,6-diphosphate inSaccharomyces cerevisiae. Eur J Biochem 1971, 18:59–64.View Article
- Gad A: Stimulation of yeast phosphofructokinase activity by Fructose 2,6-bisphosphate. Biochem Biophys Res Commun 1981, 102:985–991.View Article
- Francois J, Vanschaftingen E, Hers HG: The mechanism by which glucose increases fructose 2,6-bisphosphate concentration inSaccharomyces cerevisiae- A cyclic-AMP-dependent activation of phosphofructokinase-2. Eur J Biochem 1984, 145:187–193.View Article
- Busti S, Coccetti P, Alberghina L, Vanoni M: Glucose signaling-mediated coordination of cell growth and cell cycle inSaccharomyces cerevisiae. Sensors 2010, 10:6195–6240.View Article
- Jeffries TW, Fady JH, Lightfoot EN: Effect of glucose supplements on the fermentation of xylose byPachysolen tannophilus. Biotechnology and Bioengineering 1985, 27:171–176.View Article
- Ohgren K, Bengtsson O, Gorwa-Grauslund MF, Galbe M, Hahn-Hagerdal B, Zacchi G: Simultaneous saccharification and co-fermentation of glucose and xylose in steam-pretreated corn stover at high fiber content withSaccharomyces cerevisiaeTMB3400. Journal of Biotechnology 2006, 126:488–498.View Article
- Krahulec S, Petschacher B, Wallner M, Longus K, Klimacek M, Nidetzky B: Fermentation of mixed glucose-xylose substrates by engineered strains ofSaccharomyces cerevisiae: Role of the coenzyme specificity of xylose reductase, and effect of glucose on xylose utilization. Microb Cell Fact 2010, 10:9–16.
- Braus GH: Aromatic amino acid biosynthesis in the yeastSaccharomyces-cerevisiae- a model system for the regulation of a eukaryotic biosynthetic pathway. Microbiol Rev 1991, 55:349–370.
- Mauricio JC, Pareja M, Ortega JM: Changes in the intracellular concentrations of the adenosine phosphates and nicotinamide adenine dinucleotides ofSaccharomyces cerevisiaeduring batch fermentation. World J Microbiol Biotechnol 1995, 11:196–201.View Article
- Salusjarvi L, Pitkanen JP, Aristidou A, Ruohonen L, Penttila M: Transcription analysis of recombinantSaccharomyces cerevisiaereveals novel responses to xylose. Appl Biochem Biotechnol 2006, 128:237–261.View Article
- Wahlbom CF, Hahn-Hagerdal B: Furfural, 5-hydroxymethyl furfural, and acetoin act as external electron acceptors during anaerobic fermentation of xylose in recombinantSaccharomyces cerevisiae. Biotechnol Bioeng 2002, 78:172–178.View Article
- Wahlbom CF, Otero RRC, van Zyl WH, Hahn-Hagerdal B, Jonsson LJ: Molecular analysis of aSaccharomyces cerevisiaemutant with improved ability to utilize xylose shows enhanced expression of proteins involved in transport, initial xylose metabolism, and the pentose phosphate pathway. Appl Environ Microbiol 2003, 69:740–746.View Article
- Jin YS, Laplaza JM, Jeffries TW: Saccharomyces cerevisiaeengineered for xylose metabolism exhibits a respiratory response. Appl Environ Microbiol 2004, 70:6816–6825.View Article
- Camarasa C, Grivet JP, Dequin S: Investigation by 13C-NMR and tricarboxylic acid (TCA) deletion mutant analysis of pathways for succinate formation inSaccharomyces cerevisiaeduring anaerobic fermentation. Microbiology 2003, 149:2669–2678.View Article
- Minard KI, McAlisterhenn L: Glucose-induced degradation of the MDH2 isozyme of malate dehydrogenase in yeast. J Biol Chem 1992, 267:17458–17464.
- Atkinson DE: Energy charge of adenylate pool as a regulatory parameter. Interaction with feedback modifiers. Biochemistry 1968, 7:4030–4034.View Article
- Dombek KM, Ingram LO: Intracellular accumulation of AMP as a cause for the decline in rate of ethanol-production bySaccharomyces cerevisiaeduring batch fermentation. Appl Environ Microbiol 1988, 54:98–104.
- Knowles CJ: Microbial metabolic regulation by adenine nucleotide pools. In Microbial Energetics: Twenty-Seventh Symposium of the Society for General Microbiology. Edited by: Haddock BA, Hamilton WA. Cambridge University Press, Cambridge; 1977:241–283.
- Chapman AG, Atkinson DE: Adenine nucleotide concentrations and turnover rates. Their correlation with biological activity in bacteria and yeast. In Advances in microbial physiology. Edited by: Rose AH, Tempest DW. Academic Press, London; 1977:253–306.View Article
- Johnston GC, Ehrhardt CW, Lorincz A, Carter BLA: Regulation of cell-size in the yeastSacchaormyces cerevisiae. J Bacteriol 1979, 137:1–5.
- Moore SA: Kinetic evidence for a critical rate of protein synthesis in theSaccharomyces cerevisiaeyeast cell cycle. J Biol Chem 1988, 263:9674–9681.
- Rudoni S, Colombo S, Coccetti P, Martegani E: Role of guanine nucleotides in the regulation of the Ras/cAMP pathway inSaccharomyces cerevisiae. BBA Mol Cell Res 2001, 1538:181–189.
- Iglesias-Gato D, Martin-Marcos P, Santos MA, Hinnebusch AG, Tamame M: Guanine nucleotide pool imbalance impairs multiple steps of protein synthesis and disrupts GCN4 translational control inSaccharomyces cerevisiae. Genetics 2011, 187:105–122.View Article
- Yalowitz JA, Jayaram HN: Molecular targets of guanine nucleotides in differentiation, proliferation and apoptosis. Anticancer Res 2000, 20:2329–2338.
- Toda T, Uno I, Ishikawa T, Powers S, Kataoka T, Broek D, Cameron S, Broach J, Matsumoto K, Wigler M: In yeast, Ras proteins are controlling elements of adenylate cyclase. Cell 1985, 40:27–36.View Article
- Thevelein JM: The Ras-adenylate cyclase pathway and cell cycle control inSaccharomyces cerevisiae. Antonie Van Leeuwenhoek 1992, 62:109–130.View Article
- Kapp LD, Lorsch JR: The molecular mechanics of eukaryotic translation. Annu Rev Biochem 2004, 73:657–704.View Article
- Cherkasova V, Qiu H, Hinnebusch AG: Snf1 promotes phosphorylation of the alpha subunit of eukaryotic translation initiation factor 2 by activating Gcn2 and inhibiting phosphatases Glc7 and Sit4. Molecular and Cellular Biology 2010, 30:2862–2873.View Article
- Henderson JF, Paterson ARP: Nucleotide metabolism: an introduction. Academic Press, New York; 1973.
- Thompson FM, Atkinson DE: Response of nucleoside diphosphate kinase to adenylate energy charge. Biochem Biophys Res Commun 1971, 45:1581.View Article
- Varma A, Freese EB, Freese E: Partial deprivation of GTP initiates meiosis and sporulation inSaccharomyces cerevisiae. Mol Gen Genet 1985, 201:1–6.View Article
- Breton A, Pinson B, Coulpier F, Giraud MF, Dautant A, Daignan-Fornier B: Lethal accumulation of guanylic nucleotides inSaccharomyces cerevisiaeHPT1-deregulated mutants. Genetics 2008, 178:815–824.View Article
- Saint-Marc C, Pinson B, Coulpier F, Jourdren L, Lisova O, Daignan-Fornier B: Phenotypic consequences of purine nucleotide imbalance inSaccharomyces cerevisiae. Genetics 2009, 183:529–538.View Article
- Gonzalez B, Francois J, Renaud M: A rapid and reliable method for metabolite extraction in yeast using boiling buffered ethanol. Yeast 1997, 13:1347–1355.View Article
- Miller JN, Miller JC: Statistics and Chemometrics for Analytical Chemistry. 5th edition. Pearson Education Limited, Harlow; 2005.
- Rousseeuw PJ: Tutorial to robust statistics. J Chemometr 1991, 5:1–20.View Article
- Ortiz MC, Sarabia LA, Herrero A: Robust regression techniques - A useful alternative for the detection of outlier data in chemical analysis. Talanta 2006, 70:499–512.View Article
- Scholz M, Selbig J: Visualization and analysis of molecular data. In Metabolomics. Humana Press, Totowa, USA; 2007:87–104. [Walker JM (Series Editor): Methods in Molecular Biology, vol 358.]View Article
- Jolliffe IT: Principal component analysis. 2nd edition. Springer-Verlag, New York; 2002.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.