Spatial desynchronization of glycolytic waves as revealed by Karhunen-Loeve analysis

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.

Nonequilibrium thermodynamics of glycolytic traveling wave: Benjamin-Feir instability

Evolution of the nonequilibrium thermodynamic entities corresponding to dynamics of the Hopf instabilities and traveling waves at a nonequilibrium steady state of a spatially extended glycolysis model is assessed here by implementing an analytically tractable scheme incorporating a complex Ginzburg-Landau equation (CGLE). In the presence of self and cross diffusion, a more general amplitude equation exploiting the multiscale Krylov-Bogoliubov averaging method serves as an essential tool to reveal the various dynamical instability criteria, especially Benjamin-Feir (BF) instability, to estimate the corresponding nonlinear dispersion relation of the traveling wave pattern. The critical control parameter, wave-number selection criteria, and magnitude of the complex amplitude for traveling waves are modified by self- and cross-diffusion coefficients within the oscillatory regime, and their variabilities are exhibited against the amplitude equation. Unlike the traveling waves, a low-amplitude broad region appears for the Hopf instability in the concentration dynamics as the system phase passes through minima during its variation with the control parameter. The total entropy production rate of the uniform Hopf oscillation and glycolysis wave not only qualitatively reflects the global dynamics of concentrations of intermediate species but almost quantitatively. Despite the crucial role of diffusion in generating and shaping the traveling waves, the diffusive part of the entropy production rate has a negligible contribution to the system's total entropy production rate. The Hopf instability shows a more complex and colossal change in the energy profile of the open nonlinear system than in the traveling waves. A detailed analysis of BF instability shows a contrary nature of the semigrand Gibbs free energy with discrete and continuous wave numbers for the traveling wave. We hope the Hopf and traveling wave pattern around the BF instability in terms of energetics and dissipation will open up new applications of such dynamical phenomena.

Partial synchronisation of glycolytic oscillations in yeast cell populations

The transition between synchronized and asynchronous behaviour of immobilized yeast cells of the strain Saccharomyces carlsbergensis was investigated by monitoring the autofluorescence of the coenzyme NADH. In populations of intermediate cell densities the individual cells remained oscillatory, whereas on the level of the cell population both a partially synchronized and an asynchronous state were accessible for experimental studies. In the partially synchronized state, the mean oscillatory frequency was larger than that of the cells in the asynchronous state. This suggests that synchronisation occurred due to entrainment by the cells that oscillated more rapidly. This is typical for synchronisation due to phase advancement. Furthermore, the synchronisation of the frequency of the glycolytic oscillations preceded the synchronisation of their phases. However, the cells did not synchronize completely, as the distribution of the oscillatory frequencies only narrowed but did not collapse to a unique frequency. Cells belonging to spatially denser clusters showed a slightly enhanced local synchronisation during the episode of partial synchronisation. Neither the clusters nor a transition from partially synchronized glycolytic oscillations to travelling glycolytic waves did substantially affect the degree of partial synchronisation. Chimera states, i.e., the coexistence of a synchronized and an asynchronous part of the population, could not be found.

Catalysis by Imaging: From Meso- to Nano-scale

In-situ imaging of catalytic reactions has provided insights into reaction front propagation, pattern formation and other spatio-temporal effects for decades. Most recently, analysis of the local image intensity opened a way towards evaluation of local reaction kinetics. Herein, our recent studies of catalytic CO oxidation on Pt(hkl) and Rh(hkl) via the kinetics by imaging approach, both on the meso- and nano-scale, are reviewed. Polycrystalline Pt and Rh foils and nanotips were used as µm- and nm-sized surface structure libraries as model systems for reactions in the 10–5–10–6 mbar pressure range. Isobaric light-off and isothermal kinetic transitions were visualized in-situ at µm-resolution by photoemission electron microscopy (PEEM), and at nm-resolution by field emission microscopy (FEM) and field ion microscopy (FIM). The local reaction kinetics of individual Pt(hkl) and Rh(hkl) domains and nanofacets of Pt and Rh nanotips were deduced from the local image intensity analysis. This revealed the structure-sensitivity of CO oxidation, both in the light-off and in the kinetic bistability: for different low-index Pt surfaces, differences of up to 60 K in the critical light-off temperatures and remarkable differences in the bistability ranges of differently oriented stepped Rh surfaces were observed. To prove the spatial coherence of light-off on nanotips, proper orthogonal decomposition (POD) as a spatial correlation analysis was applied to the FIM video-data. The influence of particular configurations of steps and kinks on kinetic transitions were analysed by using the average nearest neighbour number as a common descriptor. Perspectives of nanosized surface structure libraries for future model studies are discussed.

Reaction fronts of the autocatalytic hydrogenase reaction

We have built a model to describe the hydrogenase catalyzed, autocatalytic, reversible hydrogen oxidation reaction where one of the enzyme forms is the autocatalyst. The model not only reproduces the experimentally observed front properties, but also explains the found hydrogen ion dependence. Furthermore, by linear stability analysis, two different front types are found in good agreement with the experiments.

Spatially coupled catalytic ignition of CO oxidation on Pt: mesoscopic versus nano-scale

Spatial coupling during catalytic ignition of CO oxidation on μm-sized Pt(hkl) domains of a polycrystalline Pt foil has been studied in situ by PEEM (photoemission electron microscopy) in the 10(-5) mbar pressure range. The same reaction has been examined under similar conditions by FIM (field ion microscopy) on nm-sized Pt(hkl) facets of a Pt nanotip. Proper orthogonal decomposition (POD) of the digitized FIM images has been employed to analyze spatiotemporal dynamics of catalytic ignition. The results show the essential role of the sample size and of the morphology of the domain (facet) boundary in the spatial coupling in CO oxidation.

Diffusive coupling can discriminate between similar reaction mechanisms in an allosteric enzyme system

BACKGROUND

A central question for the understanding of biological reaction networks is how a particular dynamic behavior, such as bistability or oscillations, is realized at the molecular level. So far this question has been mainly addressed in well-mixed reaction systems which are conveniently described by ordinary differential equations. However, much less is known about how molecular details of a reaction mechanism can affect the dynamics in diffusively coupled systems because the resulting partial differential equations are much more difficult to analyze.

RESULTS

Motivated by recent experiments we compare two closely related mechanisms for the product activation of allosteric enzymes with respect to their ability to induce different types of reaction-diffusion waves and stationary Turing patterns. The analysis is facilitated by mapping each model to an associated complex Ginzburg-Landau equation. We show that a sequential activation mechanism, as implemented in the model of Monod, Wyman and Changeux (MWC), can generate inward rotating spiral waves which were recently observed as glycolytic activity waves in yeast extracts. In contrast, in the limiting case of a simple Hill activation, the formation of inward propagating waves is suppressed by a Turing instability. The occurrence of this unusual wave dynamics is not related to the magnitude of the enzyme cooperativity (as it is true for the occurrence of oscillations), but to the sensitivity with respect to changes of the activator concentration. Also, the MWC mechanism generates wave patterns that are more stable against long wave length perturbations.

CONCLUSIONS

This analysis demonstrates that amplitude equations, which describe the spatio-temporal dynamics near an instability, represent a valuable tool to investigate the molecular effects of reaction mechanisms on pattern formation in spatially extended systems. Using this approach we have shown that the occurrence of inward rotating spiral waves in glycolysis can be explained in terms of an MWC, but not with a Hill mechanism for the activation of the allosteric enzyme phosphofructokinase. Our results also highlight the importance of enzyme oligomerization for a possible experimental generation of Turing patterns in biological systems.

Inward rotating spiral waves in glycolysis

We report on the first observation of inward rotating spiral waves (antispirals) in a biochemical reaction-diffusion system. Experiments are performed with extracts from yeast cells in an open spatial reactor. By increasing the protein concentration of the extract we observe a transition from outward to inward propagating waves of glycolytic activity. Numerical simulations with an allosteric model for the phosphofructokinase can reproduce these inward propagating waves over a wide range of parameters if the octameric structure of yeast phosphofructokinase is taken into account.

Spatiotemporal dynamics of glycolytic waves provides new insights into the interactions between immobilized yeast cells and gels

The immobilization of cells or enzymes is a promising tool for the development of biosensors, yet the interactions between the fixative materials and the cells are not fully understood, especially with respect to their impact on both cell metabolism and cell-to-cell signaling. We show that the spatiotemporal dynamics of waves of metabolic synchronization of yeast cells provides a new criterion to distinguish the effect of different gels on the cellular metabolism, which otherwise could not be detected. Cells from the yeast Saccharomyces carlsbergensis were immobilized into agarose gel, silica gel (TMOS), or a mixture of TMOS and alginate. We compared these immobilized cells with respect to their ability to generate temporal, intracellular oscillations in glycolysis as well as propagating, extracellular synchronization waves. While the temporal dynamics, as measured by the period and the number of oscillatory cycles, was similar for all three immobilized cell populations, significant differences have been observed with respect to the shape of the waves, wave propagation direction and velocity in the three gel matrices used.

Phase reversal in the Selkov model with inhomogeneous influx

The dynamical reaction-diffusion Selkov system as a model describing the complex traveling wave behavior is presented. The approximate amplitude-phase solution allows us to extract the base properties of the biochemical distributed system, which determines such patterns. It is shown that this relatively simple model could describe qualitatively the main features of the glycolysis waves observed in the experiments.