Phase reversal in the Selkov model with inhomogeneous influx

Glycolytic Wave Patterns in a Simple Reaction-diffusion System with Inhomogeneous Influx: Dynamic Transitions

An inhomogeneous profile of chemostatted species generates a rich variety of patterns in glycolytic waves depicted in a Selkov reaction-diffusion framework here. A key role played by diffusion amplitude and symmetry in the chemostatted species profile in dictating the fate of local spatial dynamics involving periodic, quasiperiodic, and chaotic patterns and transitions among them are investigated systematically. More importantly, various dynamic transitions, including wave propagation direction changes, are illustrated in interesting situations. Besides numerical results, our analytical formulation of the amplitude equation connecting complex Ginzburg-Landau and Lambda-omega representation shed light on the phase dynamics of the system. This systematic study of the glycolytic reaction-diffusion wave is in line with previous experimental results in open spatial reactor and will provide a knowledge about the dynamics that shape and control biological information processing and related phenomena.

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.

Energetic and entropic cost due to overlapping of Turing-Hopf instabilities in the presence of cross diffusion

A systematic introduction to nonequilibrium thermodynamics of dynamical instabilities are considered for an open nonlinear system beyond conventional Turing pattern in presence of cross diffusion. An altered condition of Turing instability in presence of cross diffusion is best reflected through a critical control parameter and wave number containing both the self- and cross-diffusion coefficients. Our main focus is on entropic and energetic cost of Turing-Hopf interplay in stationary pattern formation. Depending on the relative dispositions of Turing-Hopf codimensional instabilities from the reaction-diffusion equation it clarifies two aspects: energy cost of pattern formation, especially how Hopf instability can be utilized to dictate a stationary concentration profile, and the possibility of revealing nonequilibrium phase transition. In the Brusselator model, to understand these phenomena, we have analyzed through the relevant complex Ginzberg-Landau equation using multiscale Krylov-Bogolyubov averaging method. Due to Hopf instability it is observed that the cross-diffusion parameters can be a source of huge change in free-energy and concentration profiles.

Mathematical model of chaotic oscillations and oscillatory entrainment in glycolysis originated from periodic substrate supply

We study the influence of periodic influx on a character of glycolytic oscillations within the forced Selkov system. We demonstrate that such a simple system demonstrates a rich variety of dynamical regimes (domains of entrainment of different order (Arnold tongues), quasiperiodic oscillations, and chaos), which can be qualitatively collated with the known experimental data. We determine detailed dynamical regimes exploring the map of Lyapunov characteristic exponents obtained in numerical simulations of the Selkov system with periodic influx. In addition, a special study of the chaotic regime and the scenario of its origin in this system was evaluated and discussed.

Traveling glycolytic waves induced by a temperature gradient and determination of diffusivities for dense media

Here we consider the spatially extended model incorporating the temperature-dependent autocatalytic coefficient into the Merkin-Needham-Scott version of the Selkov system and show that this model with temperature gradient quite reasonably explains the experimentally detected traveling glycolytic nonstationary waves, which can be attributed as kinematic ones. Additionally, we analyze the influence of possibly incorporating diffusion terms into the equations. It is shown that the value of diffusivity influences the timetable for the birth of new wave and their further evolution. This result could be used as a method for the determination of diffusivity.

Simple model for temperature control of glycolytic oscillations

We introduce the temperature-dependent autocatalytic coefficient into the Merkin-Needham-Scott version of the Selkov system and consider the resulting equations as a model for temperature-controlled, self-sustained glycolytic oscillations in a closed reactor. It has been shown that this simple model reproduces key features observed in the experiments with temperature growth: (i) exponentially decreasing period of oscillations; (ii) reversal of relative duration leading and tail fronts. The applied model also reproduces the modulations of oscillations induced by the periodic temperature change.

Self-sustained biochemical oscillations and waves with a feedback determined only by boundary conditions

We discuss the biochemical three-dimensional reaction-diffusion model, which does not provide temporal self-sustained oscillations via reaction terms. However, the self-sustained oscillations and waves could be obtained using the proper boundary conditions for systems with a finite thickness. We have carried out in our numerical simulation the results quite corresponding to the experimental ones. We discuss the range of models for which our approach is applicable.