Previous experiments demonstrated that a population of HeLa cells starved of glucose or both glucose and serum exhibited a strong heterogeneity in the glycolytic oscillations in terms of the number of oscillatory cells, periods of oscillations, and duration of oscillations. Here, we report numerical simulations of this heterogeneous oscillatory behavior in HeLa cells by using a newly developed mathematical model. It is simple enough that we can apply a mathematical analysis, but capture the core of the glycolytic pathway and the activity of the glucose transporter (GLUT). Lognormal distributions of the values of the four rate constants in the model were obtained from the experimental distributions in the periods of oscillations. Thus, the heterogeneity in the periods of oscillations can be attributed to the difference in the rate constants of the enzymatic reactions. The activity of GLUT is found to determine whether the HeLa cells were oscillatory or non-oscillatory under the same experimental conditions. Simulation with the log-normal distribution of the maximum uptake velocity of glucose and the four randomized rate constants based on the log-normal distributions successfully reproduced the time-dependent number of oscillatory cells (oscillatory ratios) under the two starving conditions. The difference in the initial values of the metabolites has little effect on the simulated results.

1
Alberts
,
E.
,
Larsson
,
C.
,
Andlid
,
T.
,
Walsh
,
M. C.
, and
Gustafsson
,
L.
, “
Effect of nutrient starvation on the cellular composition and metabolic capacity of Saccharomyces cerevisiae
,”
Appl. Environ. Microbiol.
73
(
15
),
4839
4848
(
2007
).
2
Amemiya
,
T.
,
Obase
,
K.
,
Hiramatsu
,
N.
,
Itoh
,
K.
,
Shibata
,
K.
,
Takinoue
,
M.
,
Yamamoto
,
T.
, and
Yamaguchi
,
T.
, “
Collective and individual glycolytic oscillations in yeast cells encapsulated in alginate microparticles
,”
Chaos
25
,
064606
(
2015
).
3
Amemiya
,
T.
,
Shibata
,
K.
,
Itoh
,
Y.
,
Itoh
,
K.
,
Watanabe
,
M.
, and
Yamaguchi
,
T.
, “
Primordial oscillations in life: Direct observation of glycolytic oscillations in individual HeLa cervical cancer cells
,”
Chaos
27
,
104602
(
2017
).
4
Bruni
,
P.
,
Farnararo
,
M.
,
Vasta
,
V.
, and
D’Alessandro
,
A.
, “
Increase of the glycolytic rate in human resting fibroblasts following serum stimulation
,”
FEBS Lett.
159
,
39
42
(
1983
).
5
Chance
,
B.
and
Hess
,
B.
, “
Metabolic control mechanisms: IV. The effect of glucose upon the steady state of respiratory enzymes in the ascites cell
,”
J. Biol. Chem.
234
,
2421
2427
(
1959
).
6
Chandra
,
F. A.
,
Buzi
,
G.
, and
Doyle
,
J. C.
, “
Glycolytic oscillations and limits on robust efficiency
,”
Science
333
,
187
192
(
2011
).
7
Crabtree
,
H. G.
, “
Observations on the carbohydrate metabolism of tumors
,”
Biochem. J.
23
,
536
545
(
1929
).
8
Danø
,
S.
,
Sørensen
,
P.
, and
Hynne
,
F.
, “
Sustained oscillations in living cells
,”
Nature
402
,
320
322
(
1999
).
9
De Monte
,
S.
,
d’Ovidio
,
F.
,
Danø
,
S.
, and
Sørensen
,
P.
, “
Dynamical quorum sensing: Population density encoded in cellular dynamics
,”
Proc. Natl. Acad. Sci. U.S.A.
20
,
18377
18381
(
2007
).
10
Dodd
,
B. J. T.
and
Kralj
,
J. M.
, “
Live cell imaging reveals pH oscillations in Saccharomyces cerevisiae during metabolic transitions
,”
Sci. Rep.
7
,
13922
(
2017
).
11
Diaz-Ruiz
,
R.
,
Rigoulet
,
M.
, and
Devin
,
A.
, “
The Warburg and Crabtree effects: On the origin of cancer cell energy metabolism and of yeast glucose repression
,”
Biochim. Biophys. Acta
1807
,
568
576
(
2011
).
12
du Preez
,
F. B.
,
van Niekerk
,
D. D.
, and
Snoep
,
J. L.
, “
From steady-state to synchronized yeast glycolytic oscillations II: Model validation
,”
FEBS J.
279
,
2823
2836
(
2012
).
13
Gerlinger
,
M.
 et al., “
Intratumor heterogeneity and branched evolution revealed by multiregion sequencing
,”
N. Engl. J. Med.
366
,
883
892
(
2012
).
14
Goldbeter
,
A.
,
Biochemical Oscillations and Cellular Rhythms
(
Cambridge University Press
,
1996
).
15
Goldbeter
,
A.
and
Lefever
,
R.
, “
Dissipative structure for an allosteric model
,”
Biophys. J.
12
,
1302
1315
(
1972
).
16
Gustavsson
,
A.-K.
,
van Niekerk
,
D. D.
,
Adiels
,
C. B.
,
du Preez
,
F. B.
,
Goksör
,
M.
, and
Snoep
,
J. L.
, “
Sustained glycolytic oscillations in individual isolated yeast cells
,”
FEBS J.
279
,
2837
2847
(
2012
).
17
Hald
,
B. O.
and
Sørensen
,
P. G.
, “
Modeling diauxic glycolytic oscillations in yeast
,”
Biophys. J.
99
,
3191
3199
(
2010
).
18
Hamanaka
,
R. B.
and
Chandel
,
N. S.
, “
Targeting glucose metabolism for cancer therapy
,”
J. Exp. Med.
209
,
211
215
(
2012
).
19
Heiden
,
M. G. V.
,
Cantley
,
L. C.
, and
Thompson
,
C. B.
, “
Understanding the Warburg effect: The metabolic requirements of cell proliferation
,”
Science
324
,
1029
1033
(
2009
).
20
Ibsen
,
K. H.
and
Schiller
,
K. W.
, “
Oscillations of nucleotides and glycolytic intermediates in aerobic suspensions of Ehrlich ascites tumor cells
,”
Biochim. Biophys. Acta
131
,
405
407
(
1967
).
21
Khazaei
,
T.
,
McGuigan
,
A.
, and
Mahadevan
,
R.
, “
Ensemble modeling of cancer metabolism
,”
Front. Physiol.
3
,
135
(
2012
).
22
Marín-Hernández
,
A.
,
Gallardo-Pérez
,
J. C.
,
Rodríguez-Enríquez
,
S.
,
Encalada
,
R.
,
Moreno-Sánchez
,
R.
, and
Saavedra
,
E.
, “
Modeling cancer glycolysis
,”
Biochem. Biophys. Acta
1807
,
755
767
(
2011
).
23
Marín.-Hernández
,
A.
,
López-Ramirez
,
S. Y.
,
Mazo.-Monsalvo
,
I. D.
,
Gallardo-Perez
,
J. C.
,
Rodríguez-Enriquez
,
S.
,
Moreno-Sanchez
,
R.
, and
Saavera
,
E.
, “
Modeling cancer glycolysis under hypoglycemia, and the role played by the differential expression of glycolytic isoforms
,”
FEBS J.
281
,
3325
3345
(
2014
).
24
Martin
,
R. R.
,
Montero
,
S.
,
Silva
,
E.
,
Bizzarri
,
M.
,
Cocho
,
G.
,
Mansilla
,
R.
, and
Nieto-Villar
,
J. M.
, “
Phase transitions in tumor growth: V what can be expected from cancer glycolytic oscillations?
,”
Phys. A
486
,
762
771
(
2017
).
25
Nascimento
,
R. A.
,
Özel
,
R. E.
,
Mak
,
W. H.
,
Mulato
,
M.
,
Singaram
,
B.
, and
Pourmand
,
N.
, “
Single cell glucose nanosensor verifies elevated glucose levels in individual cancer cells
,”
Nano Lett.
16
,
1194
1200
(
2016
).
26
Olsen
,
L. F.
,
Andersen
,
A. Z.
,
Lunding
,
A.
,
Brasen
,
J. C.
, and
Poulsen
,
A. K.
, “
Regulation of glycolytic oscillations by mitochondrial and plasma membrane H+-ATPases
,”
Biophys. J.
96
,
3850
3861
(
2009
).
27
Pirkmajer
,
S.
and
Chibalin
,
A. V.
, “
Serum starvation: Caveat emptor
,”
Am. J. Physiol. Cell Physiol.
301
,
C272
C279
(
2011
).
28
Ramalingam
,
V.
,
Krishnamoorthy
,
G.
, and
Govindarajulu
,
P.
, “
Enzymes of carbohydrate metabolism in human breast carcinoma: Relationship with serum hormones
,”
Neoplasma
41
,
49
54
(
1994
).
29
Richard
,
P.
,
Bakker
,
B. M.
,
Teusink
,
B.
,
van Dam
,
K.
, and
Westerhoff
,
H. V.
, “
Acetaldehyde mediates the synchronization of sustained glycolytic oscillations in populations of yeast cells
,”
Eur. J. Biochem.
235
,
238
241
(
1996
).
30
Schütze
,
J.
and
Wolf
,
J.
, “
Spatio-temporal dynamics of glycolysis in cell layers. A mathematical model
,”
Biosystems
99
,
104
108
(
2010
).
31
Sel’kov
,
E. E.
, “
Self-oscillations in glycolysis. 1. A simple kinetic model
,”
Eur. J. Biochem.
4
,
79
86
(
1968
).
32
Shaw
,
S. N.
and
Amos
,
H.
, “
Insulin stimulation of glucose entry in chick fibroblasts and HeLa cells
,”
Biochem. Biophys. Res. Commun.
53
,
357
365
(
1973
).
33
Shibata
,
K.
,
Amemiya
,
T.
,
Kawakita
,
Y.
,
Obase
,
K.
,
Itoh
,
K.
,
Takinoue
,
M.
,
Nakata
,
S.
, and
Yamaguchi
,
T.
, “
Promotion and inhibition of synchronous glycolytic oscillations in yeast by chitosan
,”
FEBS J.
285
,
2679
2690
(
2018
).
34
Shinomoto
,
S.
and
Kuramoto
,
Y.
, “
Phase transition in active rotator systems
,”
Prog. Theor. Phys.
75
,
1105
1110
(
1986
).
35
Taylor
,
A. F.
,
Tinsley
,
M. R.
,
Wang
,
F.
,
Huang
,
Z.
, and
Showalter
,
K.
, “
Dynamical quorum sensing and synchronization in large populations of chemical oscillators
,”
Science
323
,
614
617
(
2009
).
36
Termonia
,
Y.
and
Ross
,
J.
, “
Oscillations and control features in glycolysis: Numerical analysis of a comprehensive model
,”
Proc. Natl. Acad. Sci. U.S.A.
78
,
2952
2956
(
1981
).
37
Wang
,
L.
,
Schumann
,
U.
,
Liu
,
Y.
,
Prokopchuk
,
O.
, and
Steinacker
,
J. M.
, “
Heat shock protein 70 (Hsp 70) inhibits oxidative phosphorylation and compensates ATP balance through enhanced glycolytic activity
,”
J. Appl. Physiol.
113
,
1669
1676
(
2012
).
38
Warburg
,
O.
, “
On the origin of cancer cells
,”
Science
123
,
309
314
(
1956
).
39
Weber
,
A.
,
Prokazov
,
Y.
,
Zuschratter
,
W.
, and
Hauser
,
M. J. B.
, “
Desynchronisation of glycolytic oscillations in yeast cell populations
,”
PLoS One
7
,
e43276
(
2012
).
40
Wolf
,
J.
and
Heinrich
,
R.
, “
Effect of cellular interaction on glycolytic oscillations in yeast: A theoretical investigation
,”
Biochem. J.
345
,
321
334
(
2000
).
41
Wolf
,
J.
,
Passarge
,
J.
,
Somsen
,
O. J. G.
,
Snoep
,
J. L.
,
Heinrich
,
R.
, and
Westerhoff
,
H. V.
, “
Transduction of intracellular and intercellular dynamics in yeast glycolytic oscillations
,”
Biophys. J.
78
,
1145
1153
(
2000
).
42
Wu
,
R.
and
Racker
,
E.
, “
Regulatory mechanisms in carbohydrate metabolism. IV. Pasteur effect and Crabtree effect in ascites tumor cells
,”
J. Biol. Chem.
234
,
1036
1041
(
1959
).
43
Zheng
,
N.
,
Wang
,
K.
,
He
,
J.
,
Qiu
,
Y.
,
Xie
,
G.
,
Su
,
M.
,
Jia
,
W.
, and
Li
,
H.
, “
Effects of ADMA on gene expression and metabolism in serum-starved LoVo cells
,”
Sci. Rep.
6
,
25892
(
2016
).
You do not currently have access to this content.