Similar temperature dependencies of glycolytic enzymes: an evolutionary adaptation to temperature dynamics?
© Cruz et al.; licensee BioMed Central Ltd. 2012
Received: 15 June 2012
Accepted: 6 November 2012
Published: 7 December 2012
Temperature strongly affects microbial growth, and many microorganisms have to deal with temperature fluctuations in their natural environment. To understand regulation strategies that underlie microbial temperature responses and adaptation, we studied glycolytic pathway kinetics in Saccharomyces cerevisiae during temperature changes.
Saccharomyces cerevisiae was grown under different temperature regimes and glucose availability conditions. These included glucose-excess batch cultures at different temperatures and glucose-limited chemostat cultures, subjected to fast linear temperature shifts and circadian sinoidal temperature cycles. An observed temperature-independent relation between intracellular levels of glycolytic metabolites and residual glucose concentration for all experimental conditions revealed that it is the substrate availability rather than temperature that determines intracellular metabolite profiles. This observation corresponded with predictions generated in silico with a kinetic model of yeast glycolysis, when the catalytic capacities of all glycolytic enzymes were set to share the same normalized temperature dependency.
From an evolutionary perspective, such similar temperature dependencies allow cells to adapt more rapidly to temperature changes, because they result in minimal perturbations of intracellular metabolite levels, thus circumventing the need for extensive modification of enzyme levels.
Growth and survival of microorganisms is strongly affected by environmental variables such as temperature, nutrient and oxygen availability, pH and osmolarity. Since, in natural environments, these parameters are highly dynamic, microorganisms have to cope with fluctuating, often non-optimal growth conditions. Suboptimal growth temperatures have major impacts on cell physiology including decreasing membrane fluidity and a reduced efficiency of protein synthesis and folding [1–3]. In addition, the catalytic capacity of each enzyme in the cell decreases when the temperature is lowered. This temperature impact can, in many cases, be described by an Arrhenius equation .
In the past decade, the response of the mesophilic yeast Saccharomyces cerevisiae to suboptimal temperatures has been the focus of several studies [5–7]. Interest in this subject is motivated by the biotechnological applications of S. cerevisiae. In particular, brewing and winemaking are two processes in which yeast is subjected to suboptimal temperatures (typically 12 to 15°C) to obtain specific desired flavour compounds [5, 8]. Moreover, its experimental accessibility to genome-scale analysis makes S. cerevisiae an attractive model organism for systems biology studies on temperature responses. With a few exceptions [2, 9], studies on low temperature responses of S. cerevisiae have focussed on so-called cold shock experiments. In such experiments, instantaneous exposure to low temperatures triggers a general environmental stress response in addition to temperature-specific responses [2, 9–12]. To investigate long-term acclimation rather than rapid adaptation to low temperature, thereby preventing a cold shock effect, growth of S. cerevisiae has been studied at 30 and 12°C in anaerobic glucose-limited chemostat cultures . Since the maximum specific growth rate of S. cerevisiae at 12°C is circa sevenfold lower than at 30°C [6, 13], a low dilution rate of 0.03 h-1 was used for both temperatures in this chemostat study [5, 6, 14]. In anaerobic cultures, substrate-level phosphorylation in glycolysis is the sole mechanism for ATP synthesis. Tai et al. observed that, despite substantially lower specific catalytic capacities of the glycolytic enzymes at 12°C as compared to 30°C, chemostat cultures maintained the same glycolytic flux at these two temperatures. Because in these chemostat cultivations the growth rate was set to be the same for both conditions (0.03 h-1), this indicated that the biomass yield was the same for both temperatures. Moreover, yeast did not compensate for the lower temperature by increased synthesis of glycolytic enzymes. Instead, metabolic regulation, i.e., regulation by changes in the concentrations of substrates, products and effectors  was identified as the main strategy for temperature compensation. Especially for highly expressed pathways such as glycolysis, whose enzymes can account for up to 20% of the protein content of S. cerevisiae, repeated cycles of protein degradation and synthesis would represent a substantial burden. It was therefore hypothesized that the observed dominant role of metabolic regulation represents an evolutionary adaptation to environments with frequent (e.g. circadian) temperature fluctuations .
Nevertheless, physiological studies carried out at tightly controlled constant temperatures or during very fast temperature changes (e.g. cold or heat shock experiments) represent artificial conditions, considering that many microorganisms are exposed to circadian and seasonal temperature fluctuations in their natural habitats. Evolution in such habitats is likely to have resulted in regulatory strategies to optimize performance under dynamic temperature regimes. Analysing and understanding such strategies is a typical systems biology challenge, and requires integration of biological experiments with mathematical modelling [16–18].
A minimal model to describe temperature dependency of metabolic fluxes
In this rate equation, the specific format of the mechanistic kinetic function f(K P , K S , K E , S, P, E) depends on the catalytic mechanism of the reaction. The mass-action ratio (Γ) only depends on the stoichiometry and thermodynamic properties of the reaction .
Which kinetic parameters are temperature dependent?
How can temperature dependency be described for those who are?
Can the same mechanism be applied to all enzymes?
It is well known that each enzyme has an optimum temperature at which its catalytic capacity (k cat ) is maximal, while k cat decreases at values below and above the optimum temperature (for a review see ). Quantitative relations, describing the impact of temperature on k cat , can be established by in vitro activity measurements at different temperatures and at saturating reactant concentrations. There are, however, only few data available on in vitro enzyme activity measurements at different temperatures for yeast glycolytic enzymes. A recent study focused on only two different temperatures (12 and 30°C), thus precluding the derivation of a proper quantitative relation for the temperature dependency of k cat .
Even when information about in vitro temperature-dependencies of enzymes was available, this would not necessarily provide an accurate reflection of the in vivo situation . Analysis of the temperature impact on enzyme levels (e) needs to be done carefully, since these are intrinsically context dependent. For example, glycolytic enzyme levels in S. cerevisiae are strongly influenced by specific growth rate , nutrient limitation regimes  and number of generations . Tai et al. minimized these influences by growing S. cerevisiae in anaerobic chemostat cultures at 12 and 30°C under otherwise identical conditions. Under these conditions, the glycolytic flux was the same for both temperatures and the levels of glycolytic enzymes were found to be very similar (average absolute fold change 1.5 ± 0.1) at both temperatures. In the same experiments, k cat , estimated from in vitro enzyme assays, was 3.9 ± 0.5 fold lower at 12°C than at 30°C. It was therefore inferred that changes in enzyme levels are not the primary regulation mechanism used by cells to, at equal flux, compensate for the loss of catalytic capacity at lower temperatures. Similar conclusions were drawn by Postmus et al., when studying the impact of supra-optimal temperatures on the regulation of glycolytic flux in aerobically grown S. cerevisiae.
Thermodynamic equilibrium constants (K eq ) of selected glycolytic reactions at 12 and 30°C
When it is furthermore assumed that the pathway is irreversible (as, for instance, anaerobic yeast glycolysis), the concentration of P plays no role. This set of algebraic relations can then be solved to obtain the intracellular metabolite levels at the reference temperature, for given values of the extracellular concentration of S, thus yielding relations for the intracellular levels of A and B as a function of the extracellular concentration of substrate S for the reference temperature. It should be noted here that this conclusion is also valid for pathways containing reversible reactions, as long as the steady-state assumption is fulfilled.
If each enzyme has its own temperature dependent function R i (T), the ratio between A, B and S will be different for each temperature and dependent on the parameters of the corresponding temperature function. If, on the other hand, the temperature function of the single enzymes follow the same mechanism with the same parameters, e.g. R 1 (T) = R 2 (T) = R 3 (T), then eq. (6) will become equal to eq. (5). This would then imply that the relations describing the intracellular concentrations of A and B as a function of the extracellular substrate concentration are temperature independent. The differences between these two scenarios can be tested in vivo to provide insight into the temperature dependency of enzymes.
Considering the simple network described above, we now assume a dynamic temperature situation in which the concentration of S is high and non-limiting and thus fluxes v 1 to v 3 are at their maximum values. If the catalytic capacities of the enzymes in a pathway or network share the same temperature dependency, changes in temperature will cause all enzyme activities to change with the same factor. In such a scenario, no changes in intracellular metabolite levels will occur and thus the cells maintain complete homeostasis during temperature changes, without the need to adjust enzyme levels. If, however, for a substrate-limited system (i.e. the flux through the pathway is limited by the supply of S) it is experimentally imposed that the fluxes remain the same even at lower temperature (e.g. in a chemostat cultures at constant dilution rate), S as well as the metabolite concentrations A and B will change to compensate for the lower enzymatic capacities.
In silico evaluation of different k cat -temperature relationships in yeast glycolysis
To evaluate the impact of different or identical temperature dependencies of the k cat of the enzymes in yeast glycolysis, simulations were carried out with a detailed kinetic model of this pathway developed by Teusink et al.. Several modifications were made to the published version of the model to, amongst others, account for biomass growth and to include temperature dependency of enzyme capacities (see Materials and Methods). Using this model different substrate feeding regimes and temperature profiles were simulated to distinguish between the effects of temperature and of extracellular glucose concentration on glycolytic flux. This distinction is especially important for glucose-limited conditions, where the glycolytic flux is highly correlated with the extracellular glucose concentration because, as shown by  and , glucose transport is a key growth-rate-controlling process. Therefore, in glucose-limited chemostat cultures, temperature dynamics will affect the extracellular glucose concentration via changes in glycolytic capacity. These changes in the extracellular glucose concentration will then propagate through all intracellular metabolite levels, yielding temperature-independent relations between intracellular metabolite and extracellular glucose levels. If, on the other hand, a different relationship between intracellular metabolite levels and extracellular glucose concentrations is observed for different temperature profiles, this would indicate that the temperature impact on k cat differs for the different enzymes in yeast glycolysis (see simple example in the preceding paragraph).
Two distinct rounds of simulations with different substrate feeding regimes were performed to evaluate the impact of the parameters of the temperature-dependent function R i (T) on the intracellular metabolite profiles. In the first round, the temperature sensitivity of one of the glycolytic enzymes was set to be different from all others. In the second round, temperature sensitivities were set to be identical for all enzymes, implying that a decrease in temperature resulted in the same relative decrease of k cat for all enzymes. In both rounds, simulations were performed for circadian temperature cycles, linear temperature shifts starting from 30°C steady-state conditions and batch cultures grown at 12, 18, 24 and 30°C (temperature profiles are illustrated in Figure 1).
When temperature sensitivities were assumed to be identical for all glycolytic enzymes, the intracellular concentrations of glycolytic intermediates under glucose-excess conditions were predicted to be temperature independent (Figure 2B). Conversely, simulations of glucose-limited conditions revealed strong changes of intracellular metabolite concentrations during the temperature shifts. The range of intracellular levels was broader than in the previous simulation (Figure 2A). Nevertheless, because the temperature sensitivity was set to be identical for all enzymes, intracellular metabolite concentrations showed highly similar correlations with the extracellular glucose concentration, irrespective of the simulated temperature profile. The results of these simulations agree with the results obtained with the simple example network (equation 3), thus showing that model predictions are independent of the assumed kinetic mechanisms of the enzymes involved, as long as the assumption of temperature independent K i holds true. Therefore, experimental analysis of the relations between the extracellular glucose concentration and intracellular metabolite levels for different cultivation conditions and temperature regimes should resolve the question whether or not all enzymes in a pathway share the same k cat -temperature relationship. Although based on model simulations the major differences between the two hypotheses are expected in glucose-excess conditions, in vivo there might be more than one enzyme with significantly different temperature dependency. Therefore, the deviating trends of intracellular metabolite levels simulated in glucose-limited conditions compared to glucose-excess conditions might be augmented. Experimental analysis under both conditions is important to decipher and quantify the temperature dependency of the glycolytic enzymes.
Experimental evaluation of model predictions
To experimentally investigate k cat -temperature relationships in yeast glycolysis, anaerobic cultures were grown under different substrate feeding regimes and dynamic temperature conditions. The experimental setups used included batch cultures at 12, 18, 24 and 30°C, to evaluate the temperature impact on intracellular metabolite levels under glucose excess conditions (sequential batch experiments (SBR)) and different sets of dynamic temperature shifts between 12 and 30°C performed in glucose-limited chemostat cultures. These dynamic experiments comprised short term (3 h) linear temperature shifts (12→30→12°C) and (30→12→30°C) applied to steady-state cultures at 12 and 30°C (LTS12 and LTS30) and sinoidal circadian temperature cycles (CTC) with 24 h frequency (Figure 1). Experiments LTS12 and LTS30 aimed at understanding if different pre-cultivation temperatures would affect the metabolic response to temperature shifts. The different timescale of the temperature perturbations between LTS and CTC allowed to better discriminate between impact of temperature and substrate concentration (see previous paragraphs). The CO2 production rate was monitored on-line during all experiments since, in anaerobic cultures, it provides an accurate measure for (changes in) glycolytic flux.
In the SBR experiments, glucose was present in excess, and therefore the glycolytic flux depended exclusively on the cultivation temperature, through its impact on k cat . Even though the glycolytic flux in batch cultures was observed to be six-fold lower at 12°C than at 30°C, intracellular metabolite levels were independent of the growth temperature (Figure 4D).
Considering that chemical reactions in aqueous solutions have activation energies between 42 and 125 kJ.mol-1, it can be estimated that reaction rates change by two fold for every 10°C temperature change . However, temperature effects on individual enzyme-catalysed reactions can deviate substantially from this ‘rule of thumb’, for example due to the impact of temperature on protein structure. As illustrated with a simple linear metabolic pathway (equation 3), different temperature dependencies of enzyme reactions inevitably result in (large) changes in concentrations of pathway metabolites. Similar changes in metabolite levels are observed when expression levels of individual enzymes in a pathway are modified . Hierarchical regulation, i.e. regulation at the level of enzyme synthesis , would then be needed to restore and maintain homeostasis.
Previous in vitro assays in cell extracts, performed at 12 and 30°C, suggested large differences in temperature dependency for the ten glycolytic enzymes in S. cerevisiae. However, in vitro studies with optimized assays, in which parameters such as pH and concentrations of salts, cofactors and effectors differ for each enzyme, are not representative for the intracellular environment . The differences in, among others, protein content, osmotic pressure, substrate diffusion between the in vitro and in vivo conditions might lead to different flux versus temperature relationships.
Experimental results obtained under conditions of glucose excess, were consistent with this scenario and could not be reconciled with the alternative scenario in which temperature dependencies of the k cat values of individual glycolytic enzymes were different. Although the experimental results of metabolite levels versus residual glucose for the glucose limited chemostat experiments (LTS and CTC) showed similar trends as the model simulations, the experimental errors of the metabolite data did not allow discriminating between the two simulated scenarios (i.e. the same temperature dependency for all enzymes versus a single enzyme different). However, if several glycolytic enzymes have different temperature sensitivities, the differences should be larger. Unfortunately, simulations in which several glycolytic enzymes have different temperature sensitivities did not result in a stable model.
All simulations were carried out using the model parameters published in . It should be realized, however, that different initial conditions, in our case anaerobic glucose-limited chemostat cultures at 12 and 30°C, could have resulted in different expression levels of glycolytic isoenzymes, which could have led to differences in enzyme kinetic properties in both conditions. A genome wide transcriptome analysis of cells cultivated under both conditions  did not provide indications for significant differences in the expression of glycolytic isoenzymes. Nevertheless, changes were observed in the expression levels of the different hexose transporters between the 12 and 30°C chemostat cultivations. Determination of the glucose transport kinetics for both conditions indeed showed that the glucose transport capacity of cells grown in 12°C chemostats was higher than that of cells grown at 30°C . Changes in the kinetics of glucose transport could have had a repercussion on the metabolite profiles as a function of the residual glucose concentration between shifts LTS12 and LTS30. Nevertheless, it was measured that the increase in the k cat of glucose transport at 12°C was coupled to an increase in the glucose saturation constant (K S ), in such a way that the ratio K S /k cat is not significantly different for the two conditions. Moreover, the increase in glucose transport capacity observed for the 12°C chemostat (1.59) is significantly lower than the decrease in glucose transport capacity caused by the temperature decrease (6.4 fold change). Considering the similar K S /k cat ratio and the limited difference in k cat at 30°C for both conditions, the effect must have been small in our experiments.
Evolutionary ‘synchronization’ of the temperature dependencies of the catalytic capacities of enzymes in a metabolic pathway, may confer selective advantages to microorganisms that are subjected to frequent temperature changes. Firstly, changes of metabolite levels in response to temperature changes are minimized, thereby avoiding the need for energetically costly cycles of de novo synthesis and degradation of enzymes. This preservation of homeostasis may be particularly important for central metabolism, of which intermediates serve as precursors for several biosynthetic pathways. Secondly, based on metabolic control analysis, it can be inferred that this minimization of changes in metabolite concentrations causes pathway fluxes to be minimally sensitive to temperature variations (see Additional file 4).
Recent studies on temperature as an evolutionary pressure indicate that temperature optima for growth are a key factor for survival of sympatric strains [34–36]. Our results indicate that, in addition to temperature optima, microbial strategies for kinetic adaptation to temperature dynamics (e.g. circadian temperature cycles) may play a key role in determining microbial fitness in natural environments. Further research should reveal whether microorganisms evolved in environments with a steady temperature (e.g. obligate commensals of homoeothermic vertebrates) show a different kinetic regulatory strategy than found for glycolysis in S. cerevisiae.
In view of their size, microorganisms are intrinsically unable to control their temperature. Temperature dynamics in natural environments therefore force microorganisms to continually adapt to temperature changes. In the present study, the glycolytic pathway of S. cerevisiae was used as a model to study kinetic regulation strategies of organisms in response to temperature dynamics. By combining both modeling and experimental results we obtained strong indications that the temperature dependencies of the catalytic capacities of yeast glycolytic enzymes are highly similar under in vivo conditions. These similarities allow the cells to maintain their homeostasis during circadian temperature cycles, thereby avoiding extensive changes in enzyme levels and minimizing the temperature impact on the growth rate. These results provide new insights in the robustness of yeast central carbon metabolism and open the gate for new studies on temperature as a key parameter for evolution of sympatric strains.
The published kinetic model for yeast glycolysis  was developed to simulate a buffered cell environment where growth is absent and extracellular concentrations are constant. These constraints and the fact that this model was solely based on kinetic parameters that were estimated from in vitro experiments, render it unstable in simulations of dynamic conditions. To avoid stability problems when incorporating temperature changes and to better mimic the experimental setups applied in this study, several adjustments were made to the original model.
Equation (7) indicates how the amount of component i in the broth (N i (t)) will change based on its biomass specific production or consumption rate (q i ) and the amount of biomass present in the broth (N X (t)). The terms F in C i,in and F out C i,out refer to transport of the component via, respectively, the inlet and outlet streams of the fermenter. F is the flow rate (L.h-1) whereas C i (mM) is the concentration of component i in each stream. The parameters were set to be the same as in the experimental setup. The ethanol evaporation rate was taken into account in the mass balance of this compound following the same approach as Cruz et al..
The model was implemented and run in gPROMS (Process Systems Enterprise).
Strain and growth conditions
Saccharomyces cerevisiae CEN.PK113-7D Mata[38, 39] was grown anaerobically in medium containing 0.3 g.L-1 of (NH4)2SO4, 0.3 g.L-1 of K2H2PO4, 3.0 g.L-1 NH4H2PO4, 0.5 g.L-1 MgSO4.7H2O, 0.4 g.L -1 of Tween 80, 10 mg.L-1 ergosterol and glucose (25 g.L-1). The medium was supplemented with 1 ml.L-1 each of a trace element solution and a vitamin solution  as well as with 0.15 g.L-1 of antifoam (Silcolapse 5020, Bluestar Silicones, St. Fons, France).
The sequential batch experiments were performed in 1L fermenters with a working volume of 750 ml, whereas 2L fermenters were used for chemostat cultivation (Applikon, Schiedam, The Netherlands). The working volumes for linear-temperature-shifts (LTS) and circadian-temperature-cycles (CTC) experiments were 1.0 and 1.4 L, respectively. The stirring speed was set at 600 rpm for all cultures and the pH was controlled to 5.0 through automatic addition of 2.0 M KOH or 2.0M of H2SO4 using a Biostat Bplus controller (Sartorius BBI Systems, Melsungen, Germany). The impact of temperature on the pH measurement was taken into account by sensor calibration. In order to avoid oxygen diffusion into the cultures, Norprene tubing was used for all connections and both medium vessel and fermenter were continuously sparged with pure nitrogen gas (N2) at a flow rate of 0.35 vvm via an Ion Science Saga digital flow meter (Cambridge, UK).
An overview of the different temperature studies made is given in Figure 1. For the experiments in glucose excess conditions (SBR), two sequential batch cycles were run at each temperature (30, 24, 18 and 12°C) before samples were taken in the third cycle, to assure that cells were fully adapted to the new temperature . Duplicate samples were taken at three different time points of the exponential phase for each temperature.
The linear-temperature-shifts (LTS) and circadian-temperature-cycles (CTC) experiments were performed in glucose-limited cultures grown at a dilution rate of 0.03 h-1. After 5 residence times at constant temperature (30 or 12°C), the temperature in LTS experiments was linearly increased or decreased at a rate of 0.2°C.min-1. Samples were taken when temperature reached 30, 24, 18 and 12°C and up to 3 h after the temperature returned to the initial set-point. The temperature profiles in the SBR and LTS experiments were defined and controlled via the MFCS/win 2.1 software (Sartorius BBI Systems, Melsungen, Germany).
For CTC-experiments, a pre-programmed sinoid temperature profile (temperature (°C) = 21+9sin(π/12·t (h) +1.57) was started after 3 residence times at a constant temperature of 30°C. This profile was designed to mimic a circadian temperature cycle. Low-temperature thermostats (Lauda RE304, Lauda-Königshofen, Germany) ensured that the temperature was precisely controlled throughout the experiments. Samples were taken during the 5th and the 6th temperature cycle, by which time carbon dioxide profiles and metabolite concentrations during consecutive cycles were highly similar. To minimize disturbance, sampling volumes did not exceed 5% of the reactor volume during a single temperature cycle and minimum intervals of 3 h were maintained between sampling points.
Extracellular glucose was measured in 2 ml of broth samples, rapidly taken with syringes containing steel beads at −20°C . The number of beads was adjusted for each initial broth temperature, such that the temperature of all samples would decrease instantaneously to 1°C. Residual glucose concentrations were measured via high-performance liquid chromatography with a Bio-Rad Aminex HPX-87H column at 60°C. The column was eluted with 5 mM phosphoric acid at a flow rate of 0.6 ml.min–1. Glucose was detected with a Waters 2410 refractive index detector. Biomass dry weight was measured in duplicate samples as described in .
Intracellular metabolite samples were taken by withdrawing 1.2 ml of broth directly to 6 ml of 100% methanol at −40°C via a rapid sampling setup. Samples were washed with cold methanol and extracted with boiling ethanol as described in . The concentrations of glucose-6-phosphate (G6P), fructose-6-phosphate (F6P), fructose-1,6-bisphosphate (FBP), phospho-enol-pyruvate (PEP), glycerol-3-phosphate (G3P), malate (MAL), fumarate (FUM), succinate (SUC), α-ketoglutarate (αKG), citrate (CIT), glucose-1-phosphate (G1P), UDP-glucose (UDPgluc), trehalose-6-phosphate (T6P), mannose-6-phosphate (M6P), 6-phosphogluconate (P6G) and sedoheptulose-7-phosphate (S7P) were measured by LC-MS according to the protocol developed by . The nucleotide concentrations (ATP, ADP and AMP) were measured according to . Uniformly labeled 13C-cell-extract was applied in both analytical platforms as internal standard .
In vitro enzyme assays of the glycolytic enzymes were performed with freshly prepared cell extracts on a Hitachi model 100–60 spectrophotometer at 30°C and 340 nm (ε340 of reduced pyridine-dinucleotide cofactors 6.3 mM-1). All enzymes were assayed as described previously , with the exception of phosphofructokinase (PFK; EC 126.96.36.199), which was assayed according to , with minor modifications. The assay mixture contained: imidazole/HCl (pH 7.0) 50 mM, MgCl2 5 mM, NADH 0.15 mM, fructose 2,6-diphosphate 0.10 mM, aldolase (EC 188.8.131.52) (Sigma) 2.1 U ml-1, α-glycerophosphate dehydrogenase-triosephosphate isomerase, 1.2 U ml-1 and 12.4 U ml-1, respectively (Sigma) and cell extract. After recording background activity with 0.5 mM fructose 6-phosphate, the reaction was started with 1.0 mM ATP. All assays were performed at two concentrations of cell extract. Specific activities in duplicate experiments differed by less than 12%. Enzyme activities are expressed as μmol substrate converted per min per mg protein [U (mg protein)-1. Protein concentrations in cell extracts were determined according to , using dried bovine serum albumin (fatty-acid free; Sigma) as the standard.
We acknowledge Cor Ras and Reza Seifar for the LC-MS measurements, Marijke Luttik for help with enzyme assays, Erik de Hulster and Dirk Geerts for their assistance with the dynamic temperature set-ups and Annelie Verbon for carrying out part of the SBR experiments. This project was carried out within the research programmes of the Kluyver Centre for Genomics of Industrial Fermentation and the Netherlands Consortium for Systems Biology, which are both sponsored by the Netherlands Genomics Initiative.
- Thieringer HA, Jones PG, Inouye M: Cold shock and adaptation. Bioessays. 1998, 20: 49-57. 10.1002/(SICI)1521-1878(199801)20:1<49::AID-BIES8>3.0.CO;2-N.View ArticleGoogle Scholar
- Sahara T, Goda T, Ohgiya S: Comprehensive expression analysis of time-dependent genetic responses in yeast cells to low temperature. J Biol Chem. 2002, 277: 50015-50021. 10.1074/jbc.M209258200.View ArticleGoogle Scholar
- Feller G, Gerday C: Psychrophilic enzymes: hot topics in cold adaptation. Nat Rev Microbiol. 2003, 1: 200-208. 10.1038/nrmicro773.View ArticleGoogle Scholar
- Russell NJ, Harrisson P, Johnston IA, Jaenicke R, Zuber M, Franks F, Wynn-Williams D: Cold adaptation of microorganisms [and discussion]. Philos Trans R Soc Lond B Biol Sci. 1990, 326: 595-611. 10.1098/rstb.1990.0034.View ArticleGoogle Scholar
- Pizarro FJ, Jewett MC, Nielsen J, Agosin E: Growth Temperature Exerts Differential Physiological and Transcriptional Responses in Laboratory and Wine Strains of Saccharomyces cerevisiae. Appl Environ Microbiol. 2008, 74: 6358-6368. 10.1128/AEM.00602-08.View ArticleGoogle Scholar
- Tai SL, Daran-Lapujade P, Luttik MAH, Walsh MC, Diderich JA, Krijger GC, Van Gulik WM, Pronk JT, Daran JM: Control of the glycolytic flux in Saccharomyces cerevisiae grown at low temperature: a multi-level analysis in anaerobic chemostat cultures. J Biol Chem. 2007, 282: 10243-10251. 10.1074/jbc.M610845200.View ArticleGoogle Scholar
- Beltran G, Novo M, Guillamon JM, Mas A, Rozes N: Effect of fermentation temperature and culture media on the yeast lipid composition and wine volatile compounds. Int J Food Microbiol. 2008, 121: 169-177. 10.1016/j.ijfoodmicro.2007.11.030.View ArticleGoogle Scholar
- Redon M, Guillamon JM, Mas A, Rozes N: Effect of growth temperature on yeast lipid composition and alcoholic fermentation at low temperature. Eur Food Res Technol. 2011, 232: 517-527. 10.1007/s00217-010-1415-3.View ArticleGoogle Scholar
- Schade B, Jansen G, Whiteway M, Entian KD, Thomas DY: Cold adaptation in budding yeast. Mol Biol Cell. 2004, 15: 5492-5502. 10.1091/mbc.E04-03-0167.View ArticleGoogle Scholar
- Homma T, Iwahashi H, Komatsu Y: Yeast gene expression during growth at low temperature. Cryobiology. 2003, 46: 230-237. 10.1016/S0011-2240(03)00028-2.View ArticleGoogle Scholar
- Murata Y, Homma T, Kitagawa E, Momose Y, Sato MS, Odani M, Shimizu H, Hasegawa-Mizusawa M, Matsumoto R, Mizukami S, et al.: Genome-wide expression analysis of yeast response during exposure to 4°C. Extremophiles. 2006, 10: 117-128. 10.1007/s00792-005-0480-1.View ArticleGoogle Scholar
- Al-Fageeh MB, Smales CM: Control and regulation of the cellular responses to cold shock: the responses in yeast and mammalian systems. Biochem J Rev. 2006, 397: 247-259. 10.1042/BJ20060166.View ArticleGoogle Scholar
- Cruz ALB, Verbon AJ, Geurink LJ, Verheijen PJT, Heijnen JJ, Van Gulik WM: Use of sequential batch fermentations to characterize the impact of mild hypothermic temperatures on the anaerobic stoichiometry and kinetics of Saccharomyces cerevisiae. Biotechnol Bioeng. 2012, 109: 1735-1744. 10.1002/bit.24454.View ArticleGoogle Scholar
- Postmus J, Canelas AB, Bouwman J, Bakker BM, Van Gulik W, De Mattos MJT, Brul S, Smits GJ: Quantitative analysis of the high temperature-induced glycolytic flux increase in Saccharomyces cerevisiae reveals dominant metabolic regulation. J Biol Chem. 2008, 283: 23524-23532. 10.1074/jbc.M802908200.View ArticleGoogle Scholar
- Ter Kuile BH, Westerhoff HV: Transcriptome meets metabolome: hierarchical and metabolic regulation of the glycolytic pathway. FEBS Lett. 2001, 500: 169-171. 10.1016/S0014-5793(01)02613-8.View ArticleGoogle Scholar
- Bruggeman FJ, Westerhoff HV: The nature of systems biology. Trends Microbiol. 2007, 15: 45-50. 10.1016/j.tim.2006.11.003.View ArticleGoogle Scholar
- Westerhoff HV, Winder C, Messiha H, Simeonidis E, Adamczyk M, Verma M, Bruggeman FJ, Dunn W: Systems biology: the elements and principles of life. FEBS Lett. 2009, 583: 3882-3890. 10.1016/j.febslet.2009.11.018.View ArticleGoogle Scholar
- Heinemann M, Sauer U: Systems biology of microbial metabolism. Curr Opin Microbiol. 2010, 13: 337-343. 10.1016/j.mib.2010.02.005.View ArticleGoogle Scholar
- Hess B, Brand K: Enzymes and metabolite profiles. Control of energy metabolism. Edited by: Chance B, Estabrook RK, Williamson JR. 1965, Academic Press, New York, USAGoogle Scholar
- Teusink B, Passarge J, Reijenga CA, Esgalhado E, Van der Weijden CC, Schepper M, Walsh MC, Bakker BM, Van Dam K, Westerhoff HV, Snoep JL: Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry. Eur J Biochem. 2000, 267: 5313-5329. 10.1046/j.1432-1327.2000.01527.x.View ArticleGoogle Scholar
- Van Hoek P, Van Dijken JP, Pronk JT: Regulation of fermentative capacity and levels of glycolytic enzymes in chemostat cultures of Saccharomyces cerevisiae. Enzyme Microb Technol. 2000, 26: 724-736. 10.1016/S0141-0229(00)00164-2.View ArticleGoogle Scholar
- Van den Brink J, Canelas AB, Van Gulik WM, Pronk JT, Heijnen JJ, De Winde JH, Daran-Lapujade P: Dynamics of glycolytic regulation during adaptation of Saccharomyces cerevisiae to fermentative metabolism. Appl Environ Microbiol. 2008, 74: 5710-5723. 10.1128/AEM.01121-08.View ArticleGoogle Scholar
- Mashego MR, Jansen ML, Vinke JL, Van Gulik WM, Heijnen JJ: Changes in the metabolome of Saccharomyces cerevisiae associated with evolution in aerobic glucose-limited chemostats. FEMS Yeast Res. 2005, 5: 419-430. 10.1016/j.femsyr.2004.11.008.View ArticleGoogle Scholar
- Alberty RA: Thermodynamics of Biochemical Reactions. 2003, John Wiley & Sons, New Jersey, USAView ArticleGoogle Scholar
- Feller G, Gerday C: Psychrophilic enzymes: molecular basis of cold adaptation. Cell Mol Life Sci. 1997, 53: 830-841. 10.1007/s000180050103.View ArticleGoogle Scholar
- Hazel JR, Prosser CL: Molecular mechanisms of temperature compensation in poikilotherms. Physiol Rev. 1974, 54: 620-677.Google Scholar
- Feller G: Protein stability and enzyme activity at extreme biological temperatures. J Phys Condens Matter. 2010, 22: 323101-10.1088/0953-8984/22/32/323101.View ArticleGoogle Scholar
- Diderich JA, Teusink B, Valkier J, Anjos J, Spencer-Martins I, Van Dam K, Walsh MC: Strategies to determine the extent of control exerted by glucose transport on glycolytic flux in the yeast Saccharomyces bayanus. Microbiology. 1999, 145: 3447-3454.View ArticleGoogle Scholar
- Schaaff I, Heinisch J, Zimmermann FK: Overproduction of glycolytic enzymes in yeast. Yeast. 1989, 5: 285-290. 10.1002/yea.320050408.View ArticleGoogle Scholar
- Fendt SM, Buescher JM, Rudroff F, Picotti P, Zamboni N, Sauer U: Tradeoff between enzyme and metabolite efficiency maintains metabolic homeostasis upon perturbations in enzyme capacity. Mol Syst Biol. 2010, 6: 365-Google Scholar
- Connors K: Introduction to chemical kinetics. Chemical Kinetics: the study of reactioon rates in solution. Edited by: Connors K. 1990, VCH Publishers, New York, USAGoogle Scholar
- Van Eunen K, Bouwman J, Daran-Lapujade P, Postmus J, Canelas AB, Mensonides FIC, Orij R, Tuzun I, Van den Brink J, Smits GJ, et al.: Measuring enzyme activities under standardized in vivo-like conditions for systems biology. FEBS J. 2010, 277: 749-760. 10.1111/j.1742-4658.2009.07524.x.View ArticleGoogle Scholar
- Nissen TL, Schulze U, Nielsen J, Villadsen J: Flux distributions in anaerobic, glucose-limited continuous cultures of Saccharomyces cerevisiae. Microbiology. 1997, 143: 203-218. 10.1099/00221287-143-1-203.View ArticleGoogle Scholar
- Gonçalves P, Valerio E, Correia C, de Almeida JM, Sampaio JP: Evidence for divergent evolution of growth temperature preference in sympatric Saccharomyces species. PLoS One. 2011, 6: e20739-10.1371/journal.pone.0020739.View ArticleGoogle Scholar
- Libkind D, Hittinger CT, Valerio E, Goncalves C, Dover J, Johnston M, Goncalves P, Sampaio JP: Microbe domestication and the identification of the wild genetic stock of lager-brewing yeast. Proc Natl Acad Sci U S A. 2011, 108: 14539-14544. 10.1073/pnas.1105430108.View ArticleGoogle Scholar
- Salvado Z, Arroyo-Lopez FN, Guillamon JM, Salazar G, Querol A, Barrio E: Temperature adaptation markedly determines evolution within the genus Saccharomyces. Appl Environ Microbiol. 2011, 77: 2292-2302. 10.1128/AEM.01861-10.View ArticleGoogle Scholar
- Ratkowsky DA, Olley J, McMeekin TA, Ball A: Relationship between temperature and growth rate of bacterial cultures. J Bacteriol Virol. 1982, 149: 1-5.Google Scholar
- Entian K-D, Kötter P: 25 Yeast Genetic Strain and Plasmid Collections. Methods Microbiol. Volume 36. Edited by: Ian S, Michael JRS. 2007, Academic Press, London, UK, 629-666.Google Scholar
- Nijkamp J, Van den Broek M, Datema E, De Kok S, Bosman L, Luttik M, Daran-Lapujade P, Vongsangnak W, Nielsen J, Heijne W, et al.: De novo sequencing, assembly and analysis of the genome of the laboratory strain Saccharomyces cerevisiae CEN.PK113-7D, a model for modern industrial biotechnology. Microb Cell Fact. 2012, 11: 36-10.1186/1475-2859-11-36.View ArticleGoogle Scholar
- Verduyn C, Postma E, Scheffers WA, Van Dijken JP: Physiology of Saccharomyces cerevisiae in anaerobic glucose-limited chemostat cultures. J Gen Microbiol. 1990, 136: 395-403. 10.1099/00221287-136-3-395.View ArticleGoogle Scholar
- Canelas A, Ras C, ten Pierick A, Van Dam J, Heijnen J, Van Gulik W: Leakage-free rapid quenching technique for yeast metabolomics. Metabolomics. 2008, 4: 226-239. 10.1007/s11306-008-0116-4.View ArticleGoogle Scholar
- Van Dam JC, Eman MR, Frank J, Lange HC, Van Dedem GWK, Heijnen SJ: Analysis of glycolytic intermediates in Saccharomyces cerevisiae using anion exchange chromatography and electrospray ionization with tandem mass spectrometric detection. Anal Chim Acta. 2002, 460: 209-218. 10.1016/S0003-2670(02)00240-4.View ArticleGoogle Scholar
- Seifar RM, Ras C, Van Dam JC, Van Gulik WM, Heijnen JJ, Van Winden WA: Simultaneous quantification of free nucleotides in complex biological samples using ion pair reversed phase liquid chromatography isotope dilution tandem mass spectrometry. AnalBiochem. 2009, 388: 213-219.Google Scholar
- Wu L, Mashego MR, van Dam JC, Proell AM, Vinke JL, Ras C, van Winden WA, van Gulik WM, Heijnen JJ: Quantitative analysis of the microbial metabolome by isotope dilution mass spectrometry using uniformly 13C-labeled cell extracts as internal standards. Anal Biochem. 2005, 336: 164-171. 10.1016/j.ab.2004.09.001.View ArticleGoogle Scholar
- Jansen MLA, Diderich JA, Mashego M, Hassane A, De Winde JH, Daran-Lapujade P, Pronk JT: Prolonged selection in aerobic, glucose-limited chemostat cultures of Saccharomyces cerevisiae causes a partial loss of glycolytic capacity. Microbiology. 2005, 151: 1657-1669. 10.1099/mic.0.27577-0.View ArticleGoogle Scholar
- De Jong-Gubbels P, Van Rolleghem P, Heijnen S, Vandijken JP, Pronk JT: Regulation of carbon metabolism in chemostat cultures of Saccharomyces cerevisiae grown on mixtures of glucose and ethanol. Yeast. 1995, 11: 407-418. 10.1002/yea.320110503.View ArticleGoogle Scholar
- Lowry OH, Rosebrough NJ, Farr AL, Randall RJ: Protein measurement with the folin phenol reagent. J Biol Chem. 1951, 193: 265-275.Google Scholar
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.