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

Dynamic fluctuations in a bacterial metabolic network

The operation of the central metabolism is typically assumed to be deterministic, but dynamics and high connectivity of the metabolic network make it potentially prone to generating fluctuations. However, time-resolved measurements of metabolite levels in individual cells that are required to characterize such fluctuations remained a challenge, particularly in small bacterial cells. Here we use single-cell metabolite measurements based on Förster resonance energy transfer, combined with computer simulations, to explore the real-time dynamics of the metabolic network of Escherichia coli. We observe that steplike exposure of starved E. coli to glycolytic carbon sources elicits large periodic fluctuations in the intracellular concentration of pyruvate in individual cells. These fluctuations are consistent with predicted oscillatory dynamics of E. coli metabolic network, and they are primarily controlled by biochemical reactions around the pyruvate node. Our results further indicate that fluctuations in glycolysis propagate to other cellular processes, possibly leading to temporal heterogeneity of cellular states within a population.

Synchronisation of glycolytic activity in yeast cells

Glycolysis is the central metabolic pathway of almost every cell and organism. Under appropriate conditions, glycolytic oscillations may occur in individual cells as well as in entire cell populations or tissues. In many biological systems, glycolytic oscillations drive coherent oscillations of other metabolites, for instance in cardiomyocytes near anorexia, or in pancreas where they lead to a pulsatile release of insulin. Oscillations at the population or tissue level require the cells to synchronize their metabolism. We review the progress achieved in studying a model organism for glycolytic oscillations, namely yeast. Oscillations may occur on the level of individual cells as well as on the level of the cell population. In yeast, the cell-to-cell interaction is realized by diffusion-mediated intercellular communication via a messenger molecule. The present mini-review focuses on the synchronisation of glycolytic oscillations in yeast. Synchronisation is a quorum-sensing phenomenon because the collective oscillatory behaviour of a yeast cell population ceases when the cell density falls below a threshold. We review the question, under which conditions individual cells in a sparse population continue or cease to oscillate. Furthermore, we provide an overview of the pathway leading to the onset of synchronized oscillations. We also address the effects of spatial inhomogeneities (e.g., the formation of spatial clusters) on the collective dynamics, and also review the emergence of travelling waves of glycolytic activity. Finally, we briefly review the approaches used in numerical modelling of synchronized cell populations.

Complex Stability and an Irrevertible Transition Reverted by Peptide and Fibroblasts in a Dynamic Model of Innate Immunity

We here apply a control analysis and various types of stability analysis to an in silico model of innate immunity that addresses the management of inflammation by a therapeutic peptide. Motivation is the observation, both in silico and in experiments, that this therapy is not robust. Our modeling results demonstrate how (1) the biological phenomena of acute and chronic modes of inflammation may reflect an inherently complex bistability with an irrevertible flip between the two modes, (2) the chronic mode of the model has stable, sometimes unique, steady states, while its acute-mode steady states are stable but not unique, (3) as witnessed by TNF levels, acute inflammation is controlled by multiple processes, whereas its chronic-mode inflammation is only controlled by TNF synthesis and washout, (4) only when the antigen load is close to the acute mode's flipping point, many processes impact very strongly on cells and cytokines, (5) there is no antigen exposure level below which reduction of the antigen load alone initiates a flip back to the acute mode, and (6) adding healthy fibroblasts makes the transition from acute to chronic inflammation revertible, although (7) there is a window of antigen load where such a therapy cannot be effective. This suggests that triple therapies may be essential to overcome chronic inflammation. These may comprise (1) anti-immunoglobulin light chain peptides, (2) a temporarily reduced antigen load, and (3a) fibroblast repopulation or (3b) stem cell strategies.

Synthetic biology and regulatory networks: where metabolic systems biology meets control engineering

Metabolic pathways can be engineered to maximize the synthesis of various products of interest. With the advent of computational systems biology, this endeavour is usually carried out through in silico theoretical studies with the aim to guide and complement further in vitro and in vivo experimental efforts. Clearly, what counts is the result in vivo, not only in terms of maximal productivity but also robustness against environmental perturbations. Engineering an organism towards an increased production flux, however, often compromises that robustness. In this contribution, we review and investigate how various analytical approaches used in metabolic engineering and synthetic biology are related to concepts developed by systems and control engineering. While trade-offs between production optimality and cellular robustness have already been studied diagnostically and statically, the dynamics also matter. Integration of the dynamic design aspects of control engineering with the more diagnostic aspects of metabolic, hierarchical control and regulation analysis is leading to the new, conceptual and operational framework required for the design of robust and productive dynamic pathways.

Allosteric regulation of phosphofructokinase controls the emergence of glycolytic oscillations in isolated yeast cells

Oscillations are widely distributed in nature and synchronization of oscillators has been described at the cellular level (e.g. heart cells) and at the population level (e.g. fireflies). Yeast glycolysis is the best known oscillatory system, although it has been studied almost exclusively at the population level (i.e. limited to observations of average behaviour in synchronized cultures). We studied individual yeast cells that were positioned with optical tweezers in a microfluidic chamber to determine the precise conditions for autonomous glycolytic oscillations. Hopf bifurcation points were determined experimentally in individual cells as a function of glucose and cyanide concentrations. The experiments were analyzed in a detailed mathematical model and could be interpreted in terms of an oscillatory manifold in a three-dimensional state-space; crossing the boundaries of the manifold coincides with the onset of oscillations and positioning along the longitudinal axis of the volume sets the period. The oscillatory manifold could be approximated by allosteric control values of phosphofructokinase for ATP and AMP.

DATABASE

The mathematical models described here have been submitted to the JWS Online Cellular Systems Modelling Database and can be accessed at http://jjj.mib.ac.uk/webMathematica/UItester.jsp?modelName=gustavsson5. [Database section added 14 May 2014 after original online publication].

Heterogeneity of glycolytic oscillatory behaviour in individual yeast cells

There are many examples of oscillations in biological systems and one of the most investigated is glycolytic oscillations in yeast. These oscillations have been studied since the 1950s in dense, synchronized populations and in cell-free extracts, but it has for long been unknown whether a high cell density is a requirement for oscillations to be induced, or if individual cells can oscillate also in isolation without synchronization. Here we present an experimental method and a detailed kinetic model for studying glycolytic oscillations in individual, isolated yeast cells and compare them to previously reported studies of single-cell oscillations. The importance of single-cell studies of this phenomenon and relevant future research questions are also discussed.

Effect of regulatory architecture on broad versus narrow sense heritability

Additive genetic variance (VA ) and total genetic variance (VG ) are core concepts in biomedical, evolutionary and production-biology genetics. What determines the large variation in reported VA /VG ratios from line-cross experiments is not well understood. Here we report how the VA /VG ratio, and thus the ratio between narrow and broad sense heritability (h(2) /H(2) ), varies as a function of the regulatory architecture underlying genotype-to-phenotype (GP) maps. We studied five dynamic models (of the cAMP pathway, the glycolysis, the circadian rhythms, the cell cycle, and heart cell dynamics). We assumed genetic variation to be reflected in model parameters and extracted phenotypes summarizing the system dynamics. Even when imposing purely linear genotype to parameter maps and no environmental variation, we observed quite low VA /VG ratios. In particular, systems with positive feedback and cyclic dynamics gave more non-monotone genotype-phenotype maps and much lower VA /VG ratios than those without. The results show that some regulatory architectures consistently maintain a transparent genotype-to-phenotype relationship, whereas other architectures generate more subtle patterns. Our approach can be used to elucidate these relationships across a whole range of biological systems in a systematic fashion.

Exploring the genetic control of glycolytic oscillations in Saccharomyces cerevisiae

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.

Understanding complexity in neurodegenerative diseases: in silico reconstruction of emergence

Healthy functioning is an emergent property of the network of interacting biomolecules that comprise an organism. It follows that disease (a network shift that causes malfunction) is also an emergent property, emerging from a perturbation of the network. On the one hand, the biomolecular network of every individual is unique and this is evident when similar disease-producing agents cause different individual pathologies. Consequently, a personalized model and approach for every patient may be required for therapies to become effective across mankind. On the other hand, diverse combinations of internal and external perturbation factors may cause a similar shift in network functioning. We offer this as an explanation for the multi-factorial nature of most diseases: they are "systems biology diseases," or "network diseases." Here we use neurodegenerative diseases, like Parkinson's disease (PD), as an example to show that due to the inherent complexity of these networks, it is difficult to understand multi-factorial diseases with simply our "naked brain." When describing interactions between biomolecules through mathematical equations and integrating those equations into a mathematical model, we try to reconstruct the emergent properties of the system in silico. The reconstruction of emergence from interactions between huge numbers of macromolecules is one of the aims of systems biology. Systems biology approaches enable us to break through the limitation of the human brain to perceive the extraordinarily large number of interactions, but this also means that we delegate the understanding of reality to the computer. We no longer recognize all those essences in the system's design crucial for important physiological behavior (the so-called "design principles" of the system). In this paper we review evidence that by using more abstract approaches and by experimenting in silico, one may still be able to discover and understand the design principles that govern behavioral emergence.

From steady-state to synchronized yeast glycolytic oscillations I: model construction

An existing detailed kinetic model for the steady-state behavior of yeast glycolysis was tested for its ability to simulate dynamic behavior. Using a small subset of experimental data, the original model was adapted by adjusting its parameter values in three optimization steps. Only small adaptations to the original model were required for realistic simulation of experimental data for limit-cycle oscillations. The greatest changes were required for parameter values for the phosphofructokinase reaction. The importance of ATP for the oscillatory mechanism and NAD(H) for inter-and intra-cellular communications and synchronization was evident in the optimization steps and simulation experiments. In an accompanying paper [du Preez F et al. (2012) FEBS J279, 2823-2836], we validate the model for a wide variety of experiments on oscillatory yeast cells. The results are important for re-use of detailed kinetic models in modular modeling approaches and for approaches such as that used in the Silicon Cell initiative.

DATABASE

The mathematical models described here have been submitted to the JWS Online Cellular Systems Modelling Database and can be accessed at http://jjj.biochem.sun.ac.za/database/dupreez/index.html.

From steady-state to synchronized yeast glycolytic oscillations II: model validation

In an accompanying paper [du Preez et al., (2012) FEBS J279, 2810-2822], we adapt an existing kinetic model for steady-state yeast glycolysis to simulate limit-cycle oscillations. Here we validate the model by testing its capacity to simulate a wide range of experiments on dynamics of yeast glycolysis. In addition to its description of the oscillations of glycolytic intermediates in intact cells and the rapid synchronization observed when mixing out-of-phase oscillatory cell populations (see accompanying paper), the model was able to predict the Hopf bifurcation diagram with glucose as the bifurcation parameter (and one of the bifurcation points with cyanide as the bifurcation parameter), the glucose- and acetaldehyde-driven forced oscillations, glucose and acetaldehyde quenching, and cell-free extract oscillations (including complex oscillations and mixed-mode oscillations). Thus, the model was compliant, at least qualitatively, with the majority of available experimental data for glycolytic oscillations in yeast. To our knowledge, this is the first time that a model for yeast glycolysis has been tested against such a wide variety of independent data sets.

DATABASE

The mathematical models described here have been submitted to the JWS Online Cellular Systems Modelling Database and can be accessed at http://jjj.biochem.sun.ac.za/database/dupreez/index.html.

Testing biochemistry revisited: how in vivo metabolism can be understood from in vitro enzyme kinetics

A decade ago, a team of biochemists including two of us, modeled yeast glycolysis and showed that one of the most studied biochemical pathways could not be quite understood in terms of the kinetic properties of the constituent enzymes as measured in cell extract. Moreover, when the same model was later applied to different experimental steady-state conditions, it often exhibited unrestrained metabolite accumulation.

Here we resolve this issue by showing that the results of such ab initio modeling are improved substantially by (i) including appropriate allosteric regulation and (ii) measuring the enzyme kinetic parameters under conditions that resemble the intracellular environment. The following modifications proved crucial: (i) implementation of allosteric regulation of hexokinase and pyruvate kinase, (ii) implementation of V_max values measured under conditions that resembled the yeast cytosol, and (iii) redetermination of the kinetic parameters of glyceraldehyde-3-phosphate dehydrogenase under physiological conditions.

Model predictions and experiments were compared under five different conditions of yeast growth and starvation. When either the original model was used (which lacked important allosteric regulation), or the enzyme parameters were measured under conditions that were, as usual, optimal for high enzyme activity, fructose 1,6-bisphosphate and some other glycolytic intermediates tended to accumulate to unrealistically high concentrations. Combining all adjustments yielded an accurate correspondence between model and experiments for all five steady-state and dynamic conditions. This enhances our understanding of in vivo metabolism in terms of in vitro biochemistry.

Baker's yeast is widely applied in modern biotechnology, for instance for production of heterologous protein or biofuel. For such applications a thorough understanding of the central energy metabolism of the bug is crucial. Nevertheless, even for this well-known organism, attempts to build models ab initio, based on independently measured characteristics of the catalysts (the enzymes), seldom gives reliable results. A key problem in this field is that enzyme characteristics are often studied under non-physiological conditions that do not resemble the environment inside the cell. In this study we measured the enzyme characteristics under physiological conditions and assembled the results into a computational model of yeast energy metabolism. We show that this simple trick greatly improves the predictive value of the computational model. This allowed us to predict correctly how yeast cells adapt to nitrogen starvation, an industrially relevant situation, in which remodeling of the proteome strongly affects cellular energy metabolism.

Systems biology from micro-organisms to human metabolic diseases: the role of detailed kinetic models

Human metabolic diseases are typically network diseases. This holds not only for multifactorial diseases, such as metabolic syndrome or Type 2 diabetes, but even when a single gene defect is the primary cause, where the adaptive response of the entire network determines the severity of disease. The latter may differ between individuals carrying the same mutation. Understanding the adaptive responses of human metabolism naturally requires a systems biology approach. Modelling of metabolic pathways in micro-organisms and some mammalian tissues has yielded many insights, qualitative as well as quantitative, into their control and regulation. Yet, even for a well-known pathway such as glycolysis, precise predictions of metabolite dynamics from experimentally determined enzyme kinetics have been only moderately successful. In the present review, we compare kinetic models of glycolysis in three cell types (African trypanosomes, yeast and skeletal muscle), evaluate their predictive power and identify limitations in our understanding. Although each of these models has its own merits and shortcomings, they also share common features. For example, in each case independently measured enzyme kinetic parameters were used as input. Based on these 'lessons from glycolysis', we will discuss how to make best use of kinetic computer models to advance our understanding of human metabolic diseases.

Mixed and diverse metabolic and gene-expression regulation of the glycolytic and fermentative pathways in response to a HXK2 deletion in Saccharomyces cerevisiae

In Saccharomyces cerevisiae the HXK2 gene, which encodes the glycolytic enzyme hexokinase II, is involved in the regulatory mechanism known as 'glucose repression'. Its deletion leads to fully respiratory growth at high glucose concentrations where the wild type ferments profusely. Here we describe that deletion of the HXK2 gene resulted in a 75% reduction in fermentative capacity. Using regulation analysis we found that the fluxes through most glycolytic and fermentative enzymes were regulated cooperatively by changes in their capacities (Vmax) and by changes in the way they interacted with the rest of the metabolism. Glucose transport and phosphofructokinase were regulated purely at the metabolic level. The reduction of fermentative capacity in the mutant was accompanied by a remarkable resilience of the remaining capacity to nutrient starvation. After starvation, the fermentative capacity of the hxk2Delta mutant was similar to that of the wild type. Based on our results and previous reports, we suggest an inverse correlation between glucose repression and the resilience of fermentative capacity towards nutrient starvation. Only a limited number of glycolytic enzyme activities changed upon starvation of the hxk2Delta mutant and we discuss to what extent this could explain the stability of the fermentative capacity.

Introduction to systems biology

The developments in the molecular biosciences have made possible a shift to combined molecular and system-level approaches to biological research under the name of Systems Biology. It integrates many types of molecular knowledge, which can best be achieved by the synergistic use of models and experimental data. Many different types of modeling approaches are useful depending on the amount and quality of the molecular data available and the purpose of the model. Analysis of such models and the structure of molecular networks have led to the discovery of principles of cell functioning overarching single species. Two main approaches of systems biology can be distinguished. Top-down systems biology is a method to characterize cells using system-wide data originating from the Omics in combination with modeling. Those models are often phenomenological but serve to discover new insights into the molecular network under study. Bottom-up systems biology does not start with data but with a detailed model of a molecular network on the basis of its molecular properties. In this approach, molecular networks can be quantitatively studied leading to predictive models that can be applied in drug design and optimization of product formation in bioengineering. In this chapter we introduce analysis of molecular network by use of models, the two approaches to systems biology, and we shall discuss a number of examples of recent successes in systems biology.

Rhythm of the beta-cell oscillator is not governed by a single regulator: multiple systems contribute to oscillatory behavior

Pulsatile insulin output, paralleled by oscillations in intracellular Ca(2+), reflect oscillating metabolism within beta-cells in response to secretory fuels. Here we question whether oscillatory periodicity is conserved or varied from stimulation to stimulation, whether glycolysis is essential for the manifestation of an oscillatory response, and if an environment of nutrient oversupply affects oscillatory regularity. We have determined that a beta-cell oscillatory Ca(2+) pattern is independent of the type of applied secretory fuel (glucose, methyl-pyruvate, or alpha-ketoisocaproate). In addition, single cells respond with the same pattern when repeatedly stimulated, regardless of the type of stimulatory fuel. Presence of substimulatory glucose is not necessary to obtain an oscillatory responses to methyl-pyruvate or alpha-ketoisocaproate. Glucose-6-phosphate, as a measure of glycolytic flux, is not detectable under these conditions. These data suggest that multiple systems, rather than a single enzyme component, can contribute to the beta-cell oscillatory behavior. Prolonged exposure to high levels of palmitate impaired oscillatory regularity in the individual beta-cells. This supports the hypothesis that a high-fat environment might contribute to loss of regular oscillatory pattern in diabetic subjects, acting, at least in part, at the level of the single beta-cell.

The nature of systems biology

The advent of functional genomics has enabled the molecular biosciences to come a long way towards characterizing the molecular constituents of life. Yet, the challenge for biology overall is to understand how organisms function. By discovering how function arises in dynamic interactions, systems biology addresses the missing links between molecules and physiology. Top-down systems biology identifies molecular interaction networks on the basis of correlated molecular behavior observed in genome-wide "omics" studies. Bottom-up systems biology examines the mechanisms through which functional properties arise in the interactions of known components. Here, we outline the challenges faced by systems biology and discuss limitations of the top-down and bottom-up approaches, which, despite these limitations, have already led to the discovery of mechanisms and principles that underlie cell function.

Approaches to biosimulation of cellular processes

Modelling and simulation are at the heart of the rapidly developing field of systems biology. This paper reviews various types of models, simulation methods, and theoretical approaches that are presently being used in the quantitative description of cellular processes. We first describe how molecular interaction networks can be represented by means of stoichiometric, topological and kinetic models. We briefly discuss the formulation of kinetic models using mesoscopic (stochastic) or macroscopic (continuous) approaches, and we go on to describe how detailed models of molecular interaction networks (silicon cells) can be constructed on the basis of experimentally determined kinetic parameters for cellular processes. We show how theory can help in analyzing models by applying control analysis to a recently published silicon cell model. Finally, we review some of the theoretical approaches available to analyse kinetic models and experimental data, respectively.

Training of yeast cell dynamics

In both industrial fermenters and in their natural habitats, microorganisms often experience an inhomogeneous and fluctuating environment. In this paper we mimicked one aspect of this nonideal behaviour by imposing a low and oscillating extracellular glucose concentration on nonoscillating suspensions of yeast cells. The extracellular dynamics changed the intracellular dynamics--which was monitored through NADH fluorescence--from steady to equally dynamic; the latter followed the extracellular dynamics at the frequency of glucose pulsing. Interestingly, the amplitude of the oscillation of the NADH fluorescence increased with time. This increase in amplitude was sensitive to inhibition of protein synthesis, and was due to a change in the cells rather than in the medium; the cell population was 'trained' to respond to the extracellular dynamics. To examine the mechanism behind this 'training', we subjected the cells to a low and constant extracellular glucose concentration. Seventy-five minutes of adaptation to a low and constant glucose concentration induced the same increase of the amplitude of the forced NADH oscillations as did the train of glucose pulses. Furthermore, 75 min of adaptation to a low (oscillating or continuous) glucose concentration decreased the K(M) of the glucose transporter from 26 mm to 3.5 mm. When subsequently the apparent K(M) was increased by addition of maltose, the amplitude of the forced oscillations dropped to its original value. This demonstrated that the increased affinity of glucose transport was essential for the training of the cells' dynamics.