Yeast glycolytic oscillations that are not controlled by a single oscillophore: a new definition of oscillophore strength

Biochemical oscillations, such as glycolytic oscillations, are often believed to be caused by a single so-called ‘oscillophore’. The main characteristics of yeast glycolytic oscillations, such as frequency and amplitude, are however controlled by several enzymes. In this paper, we develop a method to quantify to which extent any enzyme determines the occurrence of oscillations. Principles extrapolated from metabolic control analysis are applied to calculate the control exerted by individual enzymes on the real and imaginary parts of the eigenvalues of the Jacobian matrix. We propose that the control exerted by an enzyme on the real part of the smallest eigenvalue, in terms of absolute value, quantifies to which extent that enzyme contributes to the emergence of instability. Likewise the control exerted by an enzyme on the imaginary part of complex eigenvalues may serve to quantify the extent to which that enzyme contributes to the tendency of the system to oscillate. The method was applied both to a core model and to a realistic model of yeast glycolytic oscillations. Both the control over stability and the control over oscillatory tendency were distributed among several enzymes, of which glucose transport, pyruvate decarboxylase and ATP utilization were the most important. The distributions of control were different for stability and oscillatory tendency, showing that control of instability does not imply control of oscillatory tendency nor vice versa. The control coefficients summed up to 1, suggesting the existence of a new summation theorem. These results constitute proof that glycolytic oscillations in yeast are not caused by a single oscillophore and provide a new, subtle, definition for the oscillophore strength of an enzyme.     

Control Analysis for Autonomously Oscillating Biochemical Networks

It has hitherto not been possible to analyze the control of oscillatory dynamic cellular processes in other than qualitative ways. The control coefficients, used in metabolic control analyses of steady states, cannot be applied directly to dynamic systems. We here illustrate a way out of this limitation that uses Fourier transforms to convert the time domain into the stationary frequency domain, and then analyses the control of limit cycle oscillations. In addition to the already known summation theorems for frequency and amplitude, we reveal summation theorems that apply to the control of average value, waveform, and phase differences of the oscillations. The approach is made fully operational in an analysis of yeast glycolytic oscillations. It follows an experimental approach, sampling from the model output and using discrete Fourier transforms of this data set. It quantifies the control of various aspects of the oscillations by the external glucose concentration and by various internal molecular processes. We show that the control of various oscillatory properties is distributed over the system enzymes in ways that differ among those properties. The models that are described in this paper can be accessed on http://jjj.biochem.sun.ac.za.     

Temperature dependence and temperature compensation of kinetics of chemical oscillations; Belousov-Zhabotinskii reaction, glycolysis and circadian rhythms

Based on the distribution of activation energies around the experimental mean and averaging of rate constants we propose a theoretical scheme to examine the temperature dependence and temperature compensation of time periods of chemical oscillations. The critical finite width of the distribution is characteristic of endogeneous oscillations for compensating kinetics as observed in circadian oscillations, while the vanishing width corresponds to Arrhenius temperature dependent kinetics of non-endogeneous chemical oscillation in Belousov-Zhabotinskii reaction in a CSTR or glycolysis in cell-free yeast extracts. Our theoretical analysis is corroborated with experimental data.     

Control of glycolytic oscillations by temperature

External control of oscillatory glycolysis in yeast extract has been performed by application of either homogeneous temperature oscillations or stationary, spatial temperature gradients. Entrainment of the glycolytic oscillations by the 1/2- and 1/3-harmonic, as well as the fundamental input frequency, could be observed. From the phase response curve to a single temperature pulse, a distinct sensitivity of NADH-oxidizing processes, compared with NAD-reducing processes, is visible. Determination of glycolytic intermediates shows that the feedback-regulated phosphofructokinase as well as the glyceraldehyde-3-phosphate dehydrogenase are the most temperature-sensitive steps of glycolysis. We also find strong concentration changes in ATP and AMP at varying temperatures and, accordingly, in the energy charge. Construction of a feedback loop for spatial control of temperature by means of a Peltier element allowed us to apply a temperature gradient to the yeast extract. With this setup it is possible to initiate traveling waves and to control the wave velocity.     

Long term oscillation in glycolysis

To increase the period of glycolytic oscillations in yeast extracts (Saccharomyces uvarum), the dependence of the period on pH, on concentrations of phosphate and enzymes, and on temperature has been studied. Stable oscillatory trans were obtained at a pH value of about 6.5. Increasing the phosphate and decreasing the enzyme concentrations as well as decreasing temperature lengthened the period. By dilution of the extract with buffer while maintaining the metabolite concentrations at their initial level the period could be successively prolonged from 20 min to about 6 h.     

How yeast cells synchronize their glycolytic oscillations: a perturbation analytic treatment

Of all the lifeforms that obtain their energy from glycolysis, yeast cells are among the most basic. Under certain conditions the concentrations of the glycolytic intermediates in yeast cells can oscillate. Individual yeast cells in a suspension can synchronize their oscillations to get in phase with each other. Although the glycolytic oscillations originate in the upper part of the glycolytic chain, the signaling agent in this synchronization appears to be acetaldehyde, a membrane-permeating metabolite at the bottom of the anaerobic part of the glycolytic chain. Here we address the issue of how a metabolite remote from the pacemaking origin of the oscillation may nevertheless control the synchronization. We present a quantitative model for glycolytic oscillations and their synchronization in terms of chemical kinetics. We show that, in essence, the common acetaldehyde concentration can be modeled as a small perturbation on the "pacemaker" whose effect on the period of the oscillations of cells in the same suspension is indeed such that a synchronization develops.     

Temperature compensation through systems biology

Temperature has a strong influence on most individual biochemical reactions. Despite this, many organisms have the remarkable ability to keep certain physiological fluxes approximately constant over an extended temperature range. In this study, we show how temperature compensation can be considered as a pathway phenomenon rather than the result of a single-enzyme property. Using metabolic control analysis, it is possible to identify reaction networks that exhibit temperature compensation. Because most activation enthalpies are positive, temperature compensation of a flux can occur when certain control coefficients are negative. This can be achieved in networks with branching reactions or if the first irreversible reaction is regulated by a feedback loop. Hierarchical control analysis shows that networks that are dynamic through regulated gene expression or signal transduction may offer additional possibilities to bring the apparent activation enthalpies close to zero and lead to temperature compensation. A calorimetric experiment with yeast provides evidence that such a dynamic temperature adaptation can actually occur.     

Exploring the genetic control of glycolytic oscillations in Saccharomyces Cerevisiae

ABSTRACT: BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.     

The rhythm of yeast

Although yeast are unicellular and comparatively simple organisms, they have a sense of time which is not related to reproduction cycles. The glycolytic pathway exhibits oscillatory behaviour, i.e. the metabolite concentrations oscillate around phosphofructokinase. The frequency of these oscillations is about 1 min when using intact cells. Also a yeast cell extract can oscillate, though with a lower frequency. With intact cells the macroscopic oscillations can only be observed when most of the cells oscillate in concert. Transient oscillations can be observed upon simultaneous induction; sustained oscillations require an active synchronisation mechanism. Such an active synchronisation mechanism, which involves acetaldehyde as a signalling compound, operates under certain conditions. How common these oscillations are in the absence of a synchronisation mechanism is an open question. Under aerobic conditions an oscillatory metabolism can also be observed, but with a much lower frequency than the glycolytic oscillations. The frequency is between one and several hours. These oscillations are partly related to the reproductive cycle, i.e. the budding index also oscillates; however, under some conditions they are unrelated to the reproductive cycle, i.e. the budding index is constant. These oscillations also have an active synchronisation mechanism, which involves hydrogen sulfide as a synchronising agent. Oscillations with a frequency of days can be observed with yeast colonies on plates. Here the oscillations have a synchronisation mechanism which uses ammonia as a synchronising agent.     

Calculating activation energies for temperature compensation in circadian rhythms

Many biological species possess a circadian clock, which helps them anticipate daily variations in the environment. In the absence of external stimuli, the rhythm persists autonomously with a period of approximately 24 h. However, single pulses of light, nutrients, chemicals or temperature can shift the clock phase. In the case of light- and temperature-cycles, this allows entrainment of the clock to cycles of exactly 24 h. Circadian clocks have the remarkable property of temperature compensation, that is, the period of the circadian rhythm remains relatively constant within a physiological range of temperatures. For several organisms, temperature-regulated processes within the circadian clock have been identified in recent years. However, how these processes contribute to temperature compensation is not fully understood. Here, we theoretically investigate temperature compensation in general oscillatory systems. It is known that every oscillator can be locally temperature compensated around a reference temperature, if reactions are appropriately balanced. A balancing is always possible if the control coefficient with respect to the oscillation period of at least one reaction in the oscillator network is positive. However, for global temperature compensation, the whole physiological temperature range is relevant. Here, we use an approach which leads to an optimization problem subject to the local balancing principle. We use this approach to analyse different circadian clock models proposed in the literature and calculate activation energies that lead to temperature compensation.     

Is a constant low-entropy process at the root of glycolytic oscillations?

We measured temporal oscillations in thermodynamic variables such as temperature, heat flux, and cellular volume in suspensions of non-dividing yeast cells which exhibit temporal glycolytic oscillations. Oscillations in these variables have the same frequency as oscillations in the activity of intracellular metabolites, suggesting strong coupling between them. These results can be interpreted in light of a recently proposed theoretical formalism in which isentropic thermodynamic systems can display coupled oscillations in all extensive and intensive variables, reminiscent of adiabatic waves. This interpretation suggests that oscillations may be a consequence of the requirement of living cells for a constant low-entropy state while simultaneously performing biochemical transformations, i.e., remaining metabolically active. This hypothesis, which is in line with the view of the cellular interior as a highly structured and near equilibrium system where energy inputs can be low and sustain regular oscillatory regimes, calls into question the notion that metabolic processes are essentially dissipative.     

On the mechanisms of glycolytic oscillations in yeast

This work concerns the cause of glycolytic oscillations in yeast. We analyse experimental data as well as models in two distinct cases: the relaxation-like oscillations seen in yeast extracts, and the sinusoidal Hopf oscillations seen in intact yeast cells. In the case of yeast extracts, we use flux-change plots and model analyses to establish that the oscillations are driven by on/off switching of phosphofructokinase. In the case of intact yeast cells, we find that the instability leading to the appearance of oscillations is caused by the stoichiometry of the ATP-ADP-AMP system and the allosteric regulation of phosphofructokinase, whereas frequency control is distributed over the reaction network. Notably, the NAD+/NADH ratio modulates the frequency of the oscillations without affecting the instability. This is important for understanding the mutual synchronization of oscillations in the individual yeast cells, as synchronization is believed to occur via acetaldehyde, which in turn affects the frequency of oscillations by changing this ratio.     

Glucose modulates [Ca2+]i oscillations in pancreatic islets via ionic and glycolytic mechanisms

Pancreatic islets of Langerhans display complex intracellular calcium changes in response to glucose that include fast (seconds), slow ( approximately 5 min), and mixed fast/slow oscillations; the slow and mixed oscillations are likely responsible for the pulses of plasma insulin observed in vivo. To better understand the mechanisms underlying these diverse patterns, we systematically analyzed the effects of glucose on period, amplitude, and plateau fraction (the fraction of time spent in the active phase) of the various regimes of calcium oscillations. We found that in both fast and slow islets, increasing glucose had limited effects on amplitude and period, but increased plateau fraction. In some islets, however, glucose caused a major shift in the amplitude and period of oscillations, which we attribute to a conversion between ionic and glycolytic modes (i.e., regime change). Raising glucose increased the plateau fraction equally in fast, slow, and regime-changing islets. A mathematical model of the pancreatic islet consisting of an ionic subsystem interacting with a slower metabolic oscillatory subsystem can account for these complex islet calcium oscillations by modifying the relative contributions of oscillatory metabolism and oscillatory ionic mechanisms to electrical activity, with coupling occurring via K(ATP) channels.     

Interaction of D-fructose and fructose 1-phosphate with yeast phosphofructokinase and its influence on glycolytic oscillations.

Fermentation of D-fructose- and D-glucose induced glycolytic oscillations of different period lengths in Saccharomyces carlsbergensis. Recent studies suggested, that D-fructose or one of its metabolites interacted with phosphofructokinase (ATP:D-fructo-6-phosphate 1-phosphofructokinase, EC 2.7.1.11), the core of the glycolytic 'oscillator'. In order to explore the kinetics of interaction, the influence of D-fructose and fructose 1-phosphate on purified yeast phosphofructokinase was studied. D-fructose concentrations up to 0.3 mM stimulated the enzyme, while a further increase led to competitive inhibition. The Hill coefficient for fructose 6-phosphate decreased from 2.8 to 1.0. Fructose 1-phosphate acted in a similar way, up to 1 mM activation and inhibition competitive to fructose 6-phosphate at higher concentration (2.0--3.5 mM) with the same effect on the Hill coefficient. The inhibition patterns obtained with D-fructose or fructose 1-phosphate suggest a sequential random reaction mechanism of yeast phosphofructokinase with fructose 6-phosphate and MgATP2-. The mode of interaction of phosphofructokinase with D-fructose and fructose 1-phosphate is discussed. The influence of both effectors resulted in altered enzyme kinetics, which may cause the different period lengths of glycolytic oscillations.     

The role of fructose 2,6-bisphosphate in glycolytic oscillations in extracts and cells of Saccharomyces cerevisiae

Fructose 2,6-bisphosphate is physiologically one of the most potent activators of yeast 6-phosphofructo-1-kinase. The glycolytic oscillation observed in cell-free cytoplasmic extracts of the yeast Saccharomyces cerevisiae responds to the addition of fructose 2,6-bisphosphate in micromolar concentrations by showing a pronounced decrease of both the amplitude and the period. The oscillations can be suppressed completely by 10 microM and above of this activator but recovers almost fully (95%) to the unperturbed state after 3 h. Fructose 2,6-bisphosphate shifts the phases of the oscillations by a maximal +/- 60 degrees. Oscillations in concentration of endogenous fructose 2,6-bisphosphate in the extract were also observed. Fructose 2,6-bisphosphate alters the dynamic properties of 6-phosphofructo-1-kinase which are vital for its role as the 'oscillophore'. However, the minute amount (approximately 0.3 microM) of endogenous fructose 2,6-bisphosphate and the phase relationship of its oscillations compared with other metabolites indicate that this activator is not an essential component of the oscillatory mechanism. Further support for this conclusion is the observation of sustained oscillations in both the extracts and a population of intact cells of a mutant strain (YFA) of S. cerevisiae with no detectable fructose 2,6-bisphosphate (less than 5 nM).