Glycolytic oscillation was observed first in yeast suspension by pulsing-addition of glucose to the system  , registering the transient oscillation of concentration of NAD (Nicotinamide-Adenine Dinucleotide [reduced form]). Cells were studied under anaerobic conditions and oxygen addition terminated the transient oscillation phenomenon. Soon, many biochemical oscillations were measured,   pushing the process into the front of the non-linear reaction-kinetics described by ordinary differential equations  . This activity accelerated when the first glycolytic oscillation was explained strictly on a mathematical basis called Selkov’s model  . This bifurcation-based phenomenon soon reached its direct biomedical applications too  , mainly centered on the pancreatic activity  . Mathematical models (like  ) had great support by the more and more precise experiments   . Presently, precise metabolic activity in cells is measured by showing entrainment of heterogeneous glycolytic oscillations in single cells   , in which work was one of the vast news of the year. Interaction of glycolysis with external electric field stimulation has been known for a long time    . Our objective is to clear the following question in detail: are there any similarities between the temperature stimuli and electric stimuli in transient glycolytic oscillation?
Temperature is a trivial driving parameter of glycolysis   . It was widely investigated experimentally  and explained theoretically too  as further-developed Selkov’s model. Two factors were introduced for this study: a constant α for positive feedback catalytic effect and a β(T) for temperature dependence.
where and w are constants, while β(T) satisfies the Arrhenius-dependence of temperature.
Solving these equations, the results were as follows:
1) The average concentration of the reagents decreases by time.
2) The frequency of oscillation depends on the temperature by Arrhenius function.
3) The form of oscillation changes by the temperature.
4) Oscillation can be modulated by the periodic changing of the temperature.
Let us suppose that the activation energy in Arrhenius law depends on the electric field which can have strong synergy with the temperature  .
Giving the perturbation of the temperature
after a Taylor-expansion and stop at the linear term:
When the changes have periodic part, then:
If P is the polarization and is the electric field strength, then and consequently
Comparing (4) and (6), the effects which we promote by temperature could be constructed on the same way by electric field, making substitution.
Comparing (4) and (6), electric field has a formally similar effect on the modulation of glycolytic oscillations to the temperature. Driving force is the change of
the field strength (instead of the relative temperature change),
and chemical activation energy is replaced simply by the polarization of the material. The abovementioned similarity is based on the Arrhenius equation which is a general expression for the temperature dependence of chemical reaction rates, so in the simple empirical comparison of the reaction kinetics it allows the interchange of the temperature and electric field effects. So, the effect is formally similar, but their real action is different, the electric field needs polarizability of the materials which is not the condition for temperature. The glycolysis is a special effect involving the membranes (cytoplasmic and mitochondrial) which have a definite and well distinguishable polarization vector. This means that this similarity makes high parallel effects nearby the membrane of the cells. In the other chemical mechanisms of the living material the difference between the temperature effect and the electric fields is simply the mechanism which these effects induce in the process.
The electric field in this meaning makes equivalent changes in the chemical effects of glycolytic oscillations like temperature does when the heat-exchange is concentrated on a certain chemical reaction and not spread all over the target.
Heat spreads in the target by various convection and conduction effects. Therefore temperature rise in this case cannot be equal to the precisely concentrated electric filed for chemical changes. This is indispensable for making hyperthermia treatment precise, where the accurate approximation of the chemical changes is blocked by the spreading of the heat energy. According to the Pennes Equation   :
The electric energy source has modified the Pennes Equation, and the changes are realized in reactions (introduced reaction coordinates ζi (i = 1, 2).
Moreover, two extra reaction-equations made the description valid for the system, where reactions could occur.
This interconnected (coupled) equation-system must be solved for describing the system. A more straightforward form of the conventional system is:
If the initial condition was zero product (best assumption) then the solution of the second equation is:
where the saturated (equilibrium) value and the time-constant could be temperature dependent functions. The generalized Pennes Equation with these mathematical treatments:
The Pennes Equation formally could be rewritten as follows;
where L is a linear differential-operator and g the inhomogeneity term, constructed from the blood-perfusion and the introduced electric-power. If chemical reactions (cellular disruptions) additionally happen, that makes non-linearity to the equation, and the original Equation (13) will be modified:
where N(.) term contains the non-linearities. In a general form Equation (14) could be:
where λ is an arbitrary parameter. (Naturally (14) and (15) are identical if ).
The solution of (15) could be approached by perturbation approximation having a power-series of λ:
Because L(.) is linear than:
However, N(.) is an analytical function of T; consequently, it could be presented as a series of λ powers.
Moreover, using these series in the Equation (15) grouping the terms by the power of λ as well as using the arbitrariness of λ, we get:
Solved and substituted to Equation (16) we have the solution of Equation (14) at :
The is the ordinary solution of the usual Pennes Equation, and the further terms are the consequences of the non-linear extension, so these are the essentials in our approach.
Let us study the one-dimensional case at first. The one-dimensional Pennes Equation with neglecting of the thermal conduction, using Equation (12):
which has a solution:
If the electromagnetic power is constant (not time-dependent), then:
Naturally, g could depend on the space-vector so as well.
The first term of the perturbation solution by (19) is the solution of the following equation:
Its solution (because the initial condition has been satisfied before) is:
Substitute N(.) from Equation (12) we have:
We only have numerical possibility to solve this, although it has no great importance. We must have the stationery solution, which is enough for hyperthermia conditions.
In stationery cases, is constant in time, and the first term of the integrands is zero, so:
where the only stationer solution is:
What did we get for the value of the temperature correction?
Nothing on the absolute value, because no experimental data exists. However, the relative change could be guessed by the experimental dose-function.
The correction at 43˚C for arbitrary T temperature, based on the Equation (28), is:
Using the results from the dose-idea:
where according to Sapareto and Dewey  :
Based on these equations:
For example, when we compare two temperature points: T41 = 273 + 41 = 314 K, and T45 = 273 + 45 = 318 K. In the first interval (41˚C - 43˚C) R41 = 0.25 and the second: (43˚C - 45˚C) R45 = 0.5. The blood-perfusion rate at T41wb41 = 20 ml/min/100g, and at T45 is wb45 5 ml/min/100g  .
Then the ratio of corrections according to Equation (32) is
This means if at temperature 45˚C we have an inaccuracy ΔT45C = 0.2˚C due to the chemical reactions, then the correction on T = 41˚C became ΔT41C = 3.2˚C. Due to the general hyperthermia practice, the temperature does not exceed the 41˚C in average, so the inaccuracy in these practical temperatures is large. The application of electric field enhances the accuracy to target the chemical bonds directly. Consequently, this energy-absorption has better efficacy at the same temperature when applied according to the targeted bonds. General electromagnetic radiation or capacitive heating affects all dipoles, not only the concentration for those which are devoted for cellular destruction. This is the reason why certain cell-destruction can occur at a lower temperature   , when precise impedance-matched electric field action is applied,  .
Modulation of the glycolysis both by temperature and external field is feasible. The values of the reaction change depend on the temperature and field changes by activation energy multiplied by change of temperature, and by the polarization multiplied by the change of the field. The electric field has a similar dynamic change to heating. Temperature, as the average energy of the included particles has various “waste” energy used to heat up the particles which do not participate in the desired changes. The various molecular changes in their spatiotemporal structure and the dynamical, geometrical, chemical changes could be promoted by electric field without increasing the temperature. The inaccuracy of the temperature effect could be high and must be corrected by the well-applied, electric field application.
This work was supported by the Hungarian Competitiveness and Excellence Program grant (NVKP_16-1-2016-0042).
 Gilon, P., Ravier, M.A., Jonas, J.C. and Henguin, J.C. (2002) Control Mechanisms of the Oscillations of Insulin Secretion in Vitro and in Vivo. Diabetes, 51, S144-S151.
 Westermark, P.O. and Lansner, A. (2003) A Model of Phosphofructokinase and Glycolytic Oscillations in the Pancreatic β-Cell. Biophysical Journal, 85,126-139.
 Merrins, M.J., Van Dyke, A.R., Mapp, A.K., Rizzo, M.A. and Satin, L.S. (2013) Direct Measurements of Oscillatory Glycolysis in Pancreatic Islet β-Cells Using Novel Fluorescence Resonance Energy Transfer (FRET) Biosensors for Pyruvate Kinase M2 Activity. Journal of Biological Chemistry, 288, 33312-33322.
 Olsen, L.F., Andersen, A.Z., Lunding, A., Brasen, J.C. and Poulsen, AK. (2009) Regulation of Glycolytic Oscillations by Mitochondrial and Plasma Membrane H+-ATPases. Biophysical Journal, 96, 3850-3861.
 Gustavsson, A.K., Adiels, C.B., Mehlig, B. and Goksör, M. (2015) Entrainment of Heterogeneous Glycolytic Oscillations in Single Cells. Scientific Reports, 5, Article No. 9404.
 Bertram, R., Sherman, A. and Satin, L.S. (2007) Metabolic and Electrical Oscillations: Partners in Controlling Pulsatile Insulin Secretion. American Journal of Physiology-Endocrinology and Metabolism, 293, E890-E900.
 Cruz, A.L.B., Hebly, M., Duong, G.H., Wahl, S.A., Pronk, J.T., Heijnen, J.J., Daran-Lapujade, P. and van Gulik, W.M. (2012) Similar Temperature Dependencies of Glycolytic Enzymes: An Evolutionary Adaptation to Temperature Dynamics? BMC Systems Biology, 6,151.
 Postnikov, E.B., Verveyko, D.V. and Verisokin, Y.A. (2011) Simple Model for Temperature Control of Glycolytic Oscillations. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 83, Article ID: 062901.
 Andocs, G., Renner, H., Balogh, L., Fonyad, L., Jakab, C. and Szasz, A. (2009) Strong Synergy of Heat and Modulated Electromagnetic Field in Tumor Cell Killing. Strahlentherapie und Onkologie, 185,120-126.
 Yang, K.L., Huang, C.C., Chi, M.S., Chiang, H.C., Wang, Y.S., Andocs, G., et al. (2016) In Vitro Comparison of Conventional Hyperthermia and Modulated Electro-Hyperthermia. Oncotarget, 7, 84082-84092.
 Andocs, G., Rehman, M.U., Zhao, Q.L., Tabuchi, Y., Kanamori, M. and Kondo, T. (2016) Comparison of Biological Effects of Modulated Electro-Hyperthermia and Conventional Heat Treatment in Human Lymphoma U937 Cell. Cell Death Discovery, 2, Article No. 16039.
 Szasz, A. (2015) Bioelectromagnetic Paradigm of Cancer Treatment Oncothermia. In: Rosch, P.J., Ed., Bioelectromagnetic and Subtle Energy Medicine, CRC Press, Taylor & Francis Group, Boca Raton, 323-336.