vanHeerden1

v1

v1

GLCi + Prb = G6P

v10

v10

P2G = PEP

v11

v11

PEP = Prb + PYR

v12

v12

PYR = ACE + CO2

v13

v13

{2.0}ACE + {3.0}NAD = SUCC + {3.0}NADH

v14

v14

GLCo = GLCi

v15

v15

ACE + NADH = ETOH + NAD

v16

v16

NADH + TRIO = PHOS + GLY + NAD

v17

v17

Prb = PHOS

v18

v18

X = PHOS

v2

v2

G6P = F6P

v3

v3

G6P + Prb = Glyc + {2.0}PHOS

v4

v4

Prb + {2.0}G6P = Trh + {3.0}PHOS

v5

v5

F6P + Prb = F16P

v6

v6

F16P = {2.0}TRIO

v7

v7

PHOS + TRIO + NAD = BPG + NADH

v8

v8

BPG = P3G + Prb

v9

v9

P3G = P2G

Global parameters

Assignment rules

ATP = (-SUMAXP + (1.0 - 4.0 * KEQAK) * Prb + pow(4.0 * (1.0 - 4.0 * KEQAK) * KEQAK * pow(Prb, 2.0) + pow(SUMAXP - (1.0 - 4.0 * KEQAK) * Prb, 2.0), 0.5)) / (2.0 * (1.0 - 4.0 * KEQAK))

Function definitions

Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.


Species:

Reactions:


Middle-click: pin/unpin nodes
Shift-click: pool/unpool species
Right-click: context menu

Apply alternate model layout to overlapping elements in current model:

log scales

y-axis min/max

x-axis min/max

Lost in transition: start-up of glycolysis yields subpopulations of nongrowing cells.

  • Johan H van Heerden
  • Meike T Wortel
  • Frank J Bruggeman
  • Joseph J Heijnen
  • Yves J M Bollen
  • Robert Planqué
  • Josephus Hulshof
  • Tom G O'Toole
  • S Aljoscha Wahl
  • Bas Teusink
Science 2014; 343 (6174): 1245114
Abstract
Cells need to adapt to dynamic environments. Yeast that fail to cope with dynamic changes in the abundance of glucose can undergo growth arrest. We show that this failure is caused by imbalanced reactions in glycolysis, the essential pathway in energy metabolism in most organisms. The imbalance arises largely from the fundamental design of glycolysis, making this state of glycolysis a generic risk. Cells with unbalanced glycolysis coexisted with vital cells. Spontaneous, nongenetic metabolic variability among individual cells determines which state is reached and, consequently, which cells survive. Transient ATP (adenosine 5'-triphosphate) hydrolysis through futile cycling reduces the probability of reaching the imbalanced state. Our results reveal dynamic behavior of glycolysis and indicate that cell fate can be determined by heterogeneity purely at the metabolic level.
This model reproduces Figure 1B of the paper.