Exploring the genetic control of glycolytic oscillations in Saccharomyces Cerevisiae
© Williamson et al.; licensee BioMed Central Ltd. 2012
Received: 3 April 2012
Accepted: 24 July 2012
Published: 24 August 2012
Skip to main content
© Williamson et al.; licensee BioMed Central Ltd. 2012
Received: 3 April 2012
Accepted: 24 July 2012
Published: 24 August 2012
A well known example of oscillatory phenomena is the transient oscillations of glycolytic intermediates in Saccharomyces cerevisiae, their regulation being predominantly investigated by mathematical modeling. To our knowledge there has not been a genetic approach to elucidate the regulatory role of the different enzymes of the glycolytic pathway.
We report that the laboratory strain BY4743 could also be used to investigate this oscillatory phenomenon, which traditionally has been studied using S. cerevisiae X2180. This has enabled us to employ existing isogenic deletion mutants and dissect the roles of isoforms, or subunits of key glycolytic enzymes in glycolytic oscillations. We demonstrate that deletion of TDH3 but not TDH2 and TDH1 (encoding glyceraldehyde-3-phosphate dehydrogenase: GAPDH) abolishes NADH oscillations. While deletion of each of the hexokinase (HK) encoding genes (HXK1 and HXK2) leads to oscillations that are longer lasting with lower amplitude, the effect of HXK2 deletion on the duration of the oscillations is stronger than that of HXK1. Most importantly our results show that the presence of beta (Pfk2) but not that of alpha subunits (Pfk1) of the hetero-octameric enzyme phosphofructokinase (PFK) is necessary to achieve these oscillations. Furthermore, we report that the cAMP-mediated PKA pathway (via some of its components responsible for feedback down-regulation) modulates the activity of glycoytic enzymes thus affecting oscillations. Deletion of both PDE2 (encoding a high affinity cAMP-phosphodiesterase) and IRA2 (encoding a GTPase activating protein- Ras-GAP, responsible for inactivating Ras-GTP) abolished glycolytic oscillations.
The genetic approach to characterising the glycolytic oscillations in yeast has demonstrated differential roles of the two types of subunits of PFK, and the isoforms of GAPDH and HK. Furthermore, it has shown that PDE2 and IRA2, encoding components of the cAMP pathway responsible for negative feedback regulation of PKA, are required for glycolytic oscillations, suggesting an enticing link between these cAMP pathway components and the glycolysis pathway enzymes shown to have the greatest role in glycolytic oscillation. This study suggests that a systematic genetic approach combined with mathematical modelling can advance the study of oscillatory phenomena.
Different factors have been shown to control the oscillations. These include temperature , cell density  and plasma membrane H+ ATPase and mitochondrial ATP-ase . It is believed that either acetaldehyde or ethanol acts as a synchronization mediator [15, 16] although the evidence for the role of ethanol had been disputed [15–17]. The controller of the oscillations (the so called “oscillophore”) has attracted a lot of attention but has yet to be fully established. The contributions of each element of the glycolytic pathway have been studied with mathematical models of yeast glycolytic oscillations [18–20]. Of the different glycolytic enzymes, phosphofructokinase (PFK), with its allosteric regulation (it is inhibited by ATP while AMP and fructose 1,6-bisphosphate reversing the inhibition) has been shown to have the greatest control over the oscillations. A role for the hexokinases has also been demonstrated . In contrast, there is evidence that each glycolytic enzyme plays a role in controlling the oscillatory phenomenon – this is described as “distributed control” . Both theories have been put forward as a result of analysis of the mathematical models that have been developed to explain the regulation of yeast glycolytic oscillations [17, 20].
To our knowledge there has not been previously an attempt to investigate the yeast glycolytic oscillations at genetic level. This could be due to the fact that until now, the oscillations have only been observed in the diploid S. cerevisiae strain X2180. This strain has never been used to generate targeted deletion mutants; therefore, the potential effects of deletion or over-expression of glycolytic genes have, until now, remained unknown. In contrast, almost complete collections of isogenic deletion mutants are available in the BY4743 sequenced standard laboratory genetic background .
In this study, we have demonstrated that glycolytic oscillations can be observed in the diploid S. cerevisiae strain BY4743. We have subsequently used the genetic resources available in this genetic background to investigate the effects of deletion of different glycolytic enzymes encoding genes on the NADH-mediated glycolytic oscillations in the respective mutants. We have observed differential roles of the two subunit types of the phosphofructokinase, as well as the different isoforms of the hexokinase and the glyceraldehyde-3-phosphate dehydrogenase. We have used this experimental data to evaluate via parameter sensitivity analysis and representative simulations the mathematical model of yeast glycolytic oscillations developed by Wolf et al. . Furthermore we have provided evidence for a role of the cAMP signal transduction pathway in modulating glycolytic oscillations.
Strains used in this study
MAT a/α SUC2 mal mel gal2 CUP1
MAT a/α; his3Δ1/his3Δ1; leu2Δ0/leu2Δ0; lys2Δ0/LYS2; MET15/met15Δ0; ura3Δ0/ura3Δ0
The strains were grown essentially as described by Poulsen et al. . For each strain, a single colony was used to inoculate minimal media containing appropriate supplements, and 100 mM potassium phthalate at pH 5. The cultures were incubated at 30°C overnight with shaking, and used to inoculate 200 ml of the same synthetic media. Strains were grown at 30°C with shaking until glucose was depleted (approximately 16–20 hours). The level of glucose in the media was tested with Clinistix glucose strips (Bayer). Cells were harvested by centrifugation at 5000 rpm, washed twice with buffer (50 mM K2HPO4, pH 6.8) and suspended to 10% wet weight in the same buffer. They were then incubated at 30C with shaking for three hours and kept on ice until use.
Oscillations were followed using a protocol adapted from Poulsen et al. . Following harvesting, 3 ml of yeast suspension was placed in a 4.5 ml PMMA cuvette (Fisher). The cuvette was placed in a Varian Carey Eclipse fluorescence spectrophotometer, and the temperature of the cell suspension was adjusted to 30°C. Cells were stirred at all times during the experiment using a magnetic stirrer. NADH fluorescence was followed with an excitation wavelength of 366 nm and an emission wavelength of 450 nm. Other settings were optimized so the measured intensity was always between 10 and 30 arbitrary units. During each run, the intensity was sampled 10 times every second. Oscillations were induced by the addition of glucose to a concentration of 30 mM after 60 seconds, followed by addition of KCN to a final concentration of 5 mM after 140 seconds. NADH levels were followed until the oscillations ceased (around 22 minutes for wild-type strains).
This equation gives the FFT of the data. The frequency is found by taking the maximum frequency value from a plot of frequency versus power spectrum.
The steady-state parameter sensitivity analysis was carried out using SBToolbox  in Matlab. The SBparameterestimation function was used to calculate normalized parameter sensitivities. The SBToolbox functions SBsensamplitude and SBsensperiod were used to find the sensitivities of the amplitude and periodicity of the oscillations.
where p j is the reference parameter value, δp j is the perturbation applied to the parameter, A i (p j ) is the amplitude for the reference parameter value and A i (p j + δp j ) is the amplitude for the perturbed parameter value.
The same definitions apply to the period sensitivity by replacing amplitudes with periods. These functions were used in this form by Stelling et al. .
Parameters and initial values of the Wolf model of glycolytic oscillations
Rate constant for lumped HK-PFK reaction
550 mM-1 · min-1
Inhibition constant for lumped HK-PFK reaction
Rate constant for aldolase reaction
Forward rate constant for GAPDH reaction
323.8 mM-1 · min-1
Reverse rate constant for GAPDH reaction
57823.1 mM-1 · min-1
Forward rate constant for PGK reaction
76411.1 mM-1 · min-1
Reverse rate constant for PGK reaction
23.7 mM-1 · min-1
Rate constant for PK reaction
80 mM-1 · min-1
Forward rate constant for PDC reaction
9.7 mM-1 · min-1
Forward rate constant for ADH reaction
2000 mM-1 · min-1
Rate constant for ATPase reaction
Rate constant for GLY reaction
85.7 mM-1 · min-1
Rate constant for CYA reaction
Rate constant for glucose import
50 mM min-1
Rate constant for acetaldehyde diffusion
The phosphofructokinase has long been considered to be the key regulator of yeast glycolytic oscillations [3, 29, 30]. It is a complex hetero-oligomeric enzyme composed of four alpha (encoded by PFK1) and four beta (encoded by PFK2) subunits. Unlike mathematical modelling analysis which does not discriminate the two subunit types, our genetic approach was able to test their individual contributions in glycolytic oscillations. Deletion of PFK1 resulted in a slight decrease in oscillation frequency. In contrast, deletion of PFK2 abolished oscillations completely (Figure 3b). Therefore, while the strains that have only one of the Pfk subunits can ferment glucose , the homomeric nature of the PFK enzyme in these strains results in differences in the glycolytic oscillations phenotype with the pfk2Δ but not the pfk1Δ deletion mutant being affected in NADH-mediated glycolytic oscillations. Furthermore, the characteristic “spike” of NADH that we observed in all other strains was absent in the pfk2Δ but not in the pfk1Δ strain. The pfk2Δ mutant was the only mutant displaying a drop in NADH levels after glucose addition.
Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) catalyses the reaction of glyceraldehyde-3-phosphate to 1,3 bis-phosphoglycerate. The enzyme has three isoforms encoded by three unlinked genes, TDH1, TDH2, and TDH3. All isoforms are catalytically active and have different specific enzymatic activities [32, 33]. While Tdh2p and Tdh3p are present in exponentially growing cells Tdh1p is the predominant isoform in stationary phase . As shown in Figure 3c, deletion of TDH1 and TDH2 had no effect on the glycolytic oscillations. In contrast, deletion of TDH3 abolished them almost completely.
Phosphoglycerate mutase mediates the conversion of 3-phosphoglycerate to 2-phosphoglycerate during glycolysis, and the reverse reaction during gluconeogenesis. Three genes: GPM1,2,3, code for phosphoglycerate mutase in the yeast genome. However, while GPM1 is functional and essential for viability , GPM2 and GPM3 are most-likely non-functional homologues of GPM1, since mutation in neither GPM2 nor GPM3 appears to affect glycolysis, or confer any obvious phenotype. In the current investigation, glycolytic oscillations in the gpm2Δ and gpm3Δ mutants were found to be wild type (data not shown). This observation gives further support to the theory that GPM2 and GPM3 were generated from GPM1 by a gene duplication event and then diverged from the parent copy by mutation .
There are two phosphopyruvate hydratase isoforms encoded by ENO1 and ENO2 respectively. Both function in dimeric complexes  and catalyze the conversion of 2-phosphoglycerate to phosphoenolpyruvate during glycolysis. The diploid homozygous eno1Δ was available in the EUROSCARF collection but eno2Δ was not. When tested in our experimental conditions the NADH profile of eno1Δ mutant was very similar to the wild type (data not shown).
The conversion of phosphoenolpyruvate to pyruvate (the last step in glycolysis) is catalysed by two pyruvate kinases encoded by PYK2 and CDC19, respectively. However, Cdc19p appears to be the main pyruvate kinase because the cdc19∆ mutant is not able to grow in the presence of glucose . Furthermore CDC19 and not PYK2 is tightly regulated and sensitive to levels of fructose-1,6-bisphospate. In agreement with these data, we found (data not shown) that glycolytic oscillations were present in the pyk2Δ and were comparable to the wild type.
The results described above provided a large amount of data that could on one hand be used to develop future models of yeast glycolytic oscillations, and on the other, to analyze published models. We have used the parameter sensitivity analysis tool in SBToolbox  as described in the Methods to analyze the sensitivities of the frequency and amplitude of NADH oscillations to model parameters using an existing model of glycolysis. Of the several published available models [18, 20, 22, 30, 39] we chose that of Wolf et al.  because it is comparatively easy to investigate and simulate (Table 2 and Figure 1); and because as shown by Brusch et al. , its dynamic properties are similar to those of the full scale glycolysis model .
Results of simulations of the Wolf model upon changing GAPDH levels
This study reports the results of the first genetic investigation of the pattern of oscillations of intracellular NADH of single deletion mutants in key components of the glycolysis pathway in S. cerevisiae. It has become possible because of the genetic resources available when using the sequenced genetic background strain BY4743, which as our study shows (Figure 2) displays an NADH oscillatory pattern identical to that of strain X2180 used in previously reported experiments.
One of the most significant results of this study is the genetic evidence of the differential roles of the two subunits of phosphofructokinase (PFK) in glycolytic oscillations. Previously, mathematical modelling which has demonstrated that the phosphofructokinase, with its allosteric regulation, has the greatest control over the oscillations has not discriminated between the two types of subunits that constitute this enzyme. Lack of PFK activity in vivo is only observed in double pfk1 pfk2 deletion mutants. In contrast, single deletion mutants (pfk1 and pfk2) can grow on glucose, although there is no in vitro detectable PFK activity in these mutants . Arvanitidis and Heinisch  demonstrated that each of the phosphofructokinase subunits can serve catalytic as well as regulatory functions. However, while the homomeric nature of the PFK enzyme supports growth on glucose, strains that have only one of the Pfk subunits have other different phenotypes. For example Rodicio et al.  have shown that only the pfk2 deletion mutant is temperature sensitive for growth at 37°C, sensitive to caffeine and Congo red. The pfk1 deletion mutant on the other hand is more resistant to caffeine compared to the wild type. And in the current study we show further phenotypic differences between the two types of deletion mutants. Firstly, the pfk2 but not pfk1 deletion mutant is deficient in NADH-mediated glycolytic oscillations. Secondly, pfk2 and not pfk1 mutant does not display the characteristic NADH “spike” after glucose addition (Figure 3b), an observation which will be the subject of further investigations.
The specific role of Pfk2 in glycolytic oscillations could be due to the fact that it is the subunit which contains binding sites for fructose 6-phosphate , as previous studies have shown that yeast extracts fed with the PFK substrate fructose 6-phosphate can exhibit oscillations, whereas extracts fed with fructose 1–6 bisphosphate cannot . On the other hand, the differential role of the two subunits of PFK could also be due to the differences in their associations with the cAMP pathway, which as this study shows is responsible for modulating the enzyme activities of the glycolytic enzymes in such a way that the system becomes unstable and starts to oscillate. Unlike PFK1, PFK2 is involved in genetic interactions with PDE2 and IRA2, the components of the cAMP pathway which affect NADH oscillatory activity (Figure 10). Pde2 and Ira2 are both responsible for feedback down-regulation of the cAMP-mediated activity of PKA ; hence deletion in either encoding gene causes constitutive activation of the cAMP pathway. The exact mechanism by which the cAMP-PKA pathway affects glycolytic oscillations is currently unknown, but it could be hypothesized that it might be via phosphorylation of Pfk2 by PKA, as Pfk2 (but not Pfk1) has been shown to physically interact with Tpk1. However, the fact that the pattern of glycolytic oscillations of the tpk1Δ mutant is wild type suggests that Tpk1-mediated phosphorylation of Pfk2 does not contribute to its role in glycolytic oscillations. However, it cannot be ruled out that in the absence of Tpk1, another PKA isoform (Tpk2 or Tpk3) is capable of phosphorylating Pfk2.
Deletion of each of the hexokinases encoding genes leads to oscillations that are longer lasting with lower amplitude (Figure 3a), which agrees with earlier reports which have shown that reduced hexokinase activity results in lower flow of glucose into the glycolytic pathway . Our study however, shows that the effect of HXK2 deletion on the duration of the oscillations is stronger than that of HXK1. Although the deletion mutant for GLK1 whose product is also known to phosphorylate glucose at the C6 position was not investigated, we predict that its deletion would yield results comparable to those of HXK1, as both GLK1 and HXK1 are repressed when grown on a fermentable medium using glucose as a carbon source . Interestingly, the transcript abundances of both HXK1 and GLK1 have been reported to increase, whilst those of HXK2 and PFK2 decrease respectively, upon the diauxic shift , the transition known to be regulated by the cAMP signal transduction pathway. And remarkably, HXK2 (but not HXK1), and PFK2 (but not PFK1) have been found to interact genetically with PDE2 and IRA2, deletions of which affect glycolytic oscillations (Figure 10). We therefore like to speculate that because of the role of the cAMP pathway in this oscillatory phenomenon, the glycolytic enzymes encoded by genes found to be genetically associated with PDE2 and IRA2 have the strongest regulatory role in glycolytic oscillations.
Deletion of TDH1 and TDH2 genes has little effect on glycolytic oscillations (Figure 3c). Deletion of TDH3 however leads to a loss of oscillations. These isoforms are growth-phase specific: Tdh1 is present in stationary phase cells, while Tdh2p and Tdh3p are present in exponentially growing cells . The observation that deletion of TDH3 and not that of TDH2, affects oscillations could now be explained by the fact that expression of TDH3 contributes around 60% of glyceraldehyde-3-phosphate dehydrogenase activity . As demonstrated by our simulations of the GAPDH reaction in the Wolf model of glycolysis (Figure 9), glycolytic oscillations are maintained but affected in amplitude and period upon 50% reduction of the GAPDH activity in the forward reaction, and further decrease to 30% abolished them almost completely.
The question regarding effects of deletion of ENO2 and CDC19 (PYK1) in glycolytic oscillations remains as the corresponding mutants were not available to test in the current study. It could be hypothesised that their transcripts’ down-regulation upon the diauxic shift  suggests a role in glycolytic oscillations since a similar correlation exists for PFK2. And although there might be an enticing link between transcript down-regulation upon the diauxic shift of genes such as HXK2 and PFK2 and their role in glycolytic oscillations, this relationship is not as straightforward. Our study shows that deletion of both HXK1 and HXK2 affects glycolytic oscillations, however their transcript levels respond in a contrasting manner to the diauxic transition: HXK1 is up-regulated, whilst HXK2 is down-regulated. Moreover, deletion of TDH3 abolishes glycolytic oscillations completely although neither this gene, nor the remaining two (TDH1,2) are subject to differential regulation of transcription upon the diauxic shift.
Previously, Reijenga et al.  and Madsen et al.  reported the control coefficients of enzymatic steps on the amplitude and frequency of oscillations, providing a measure of the effect that changes in enzyme activity produce on these characteristics. The deletion mutations in the present investigation represent large perturbations of the glycolysis pathway components whereas parameter sensitivities refer to changes in individual reaction parameters. This type of analysis has enabled us to distinguish between forward and reverse reaction parameters in reversible reactions, and between rate and inhibition constants for the lumped HK-PFK reaction. It has identified r1ki1, the inhibition constant for the combined HK-PFK reaction, as the most sensitive parameter. Moreover, as our experimental evidence demonstrated it is the deletion of PFK2, HXK1 and HXK2 that produce the greatest effect on glycolytic oscillations. Deletion of PFK1 on the other hand does not significantly affect glycolytic oscillations showing that the components of the lumped HK-PFK reaction have specific properties. More precise models will be needed that differentiate between the different isoforms or subunit components of glycolytic enzymes. Oscillatory parameter sensitivity analysis is a useful technique to identify crucial parameters for the control of biological rhythms and investigate their robustness. It has been applied to several other types of oscillating systems, such as circadian rhythms [28, 57], NF-κB signalling pathways  and calcium dynamics .
Interestingly, the forward constant for the lumped HK-PFK reaction (r1k1) is not as sensitive as the inhibition constant. This suggests that the reverse feedback by ATP on PFK is potentially the most important controller of the oscillations which is in agreement with results obtained in yeast extract experiments by Madsen et al. . So not only have we identified the role of the beta subunit (Pfk2) of phosphofructokinase in glycolytic oscillations (as deletion of PFK2 abolishes them completely), we have confirmed that the negative feedback of Pfk2 by ATP is a key regulatory factor in glycolytic oscillations.
The choice of the glycolysis model that we used in the current analysis might seem debatable. In light of the report by Madsen et al.  who argued firstly that there are two types of oscillations, namely relaxation type observed in cell free extracts and harmonic (sinusoidal) oscillations observed in intact cells, and secondly that the Wolf model is not applicable for oscillations in intact cells. However, it appears that the type of oscillations is not strictly coupled to the whole cells or cell extracts, but rather to how far the system is removed from the Hopf bifurcation. Close to this point the oscillations are harmonic, independent of whether they are observed in cell extracts or whole cells. Further from this point, the behaviour can move away from harmonic oscillations. This is for instance shown in Nielsen et al. , where cell free extract oscillations move from harmonic to complex and relaxation type of oscillations dependent on the distance from the Hopf bifurcation point. Furthermore, P. Richard  reports that the shape of the glycolytic oscillations varies from sinusoidal to relaxation depending on the rate of glucose addition.
This study demonstrated that r7k, the rate constant for the ATP-ase reaction, is also a highly sensitive parameter (Figures 6 and 7). Recently, Kloster and Olsen  investigated the role of intracellular ATPase activity using specific inhibitors for the three types of ATPases, namely F0F1 ATPase, Pma1 and vacuolar ATPase, and showed that both inhibiting and stimulating ATPase activity suppress glycolytic oscillations. Further support of the role of ATPase in glycolytic oscillations comes from the recent study of Ytting et al. (2012) . However, the latter study uses a GAL1 promoter replacement pma1 mutant and does not present data for its isogenic wild type parental strain. Instead the report shows as supplementary data the results of an alternative wild type strain (a derivative of W303) with a NADH oscillatory pattern which is very far from comparable to X2180 and that of BY4743 shown in our study. There are different ATP-ase encoding genes in the yeast genome; however none of the respective viable deletion mutants were included in our analysis. Given the well established role of the genetic background as a modifier of different phenotypes, we propose that in future studies it would be of interest to investigate the role of the different ATP-ases further by analysing the NADH-mediated oscillations in the standard BY4743 genetic background.
This study has for the first time provided genetic evidence for the differential roles of the two types of subunits of PFK, and the isoforms of GAPDH and HK in glycolytic oscillation in the yeast S. cerevisiae. Importantly it shows genetic evidence that PDE2 and IRA2, which encode components of the cAMP pathway responsible for negative feedback regulation of PKA, are required for glycolytic oscillations, suggesting an enticing link between these cAMP pathway components and the glycolysis pathway enzymes with a key role in glycolytic oscillation. We demonstrate that a systematic genetic approach combined with mathematical modelling can advance the study of oscillatory phenomena.
TW and DA are grateful recipients of Biotechnology and Biological Sciences Research Council (BBSRC UK) PhD studentship and a Doctoral Training Award, respectively. The authors are grateful to H. Westerhoff and J. Snoep for helpful discussions, and the anonymous reviewers for their comments and suggestions for improvement.
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.