begin name
smallbone19
end name

begin reactions
hxk {1.0}glc + {1.0}$atp = {1.0}g6p + {1.0}$adp + {1.0}$h
hxt {1.0}$glx = {1.0}glc
nth {1.0}trh + {1.0}$h2o = {2.0}glc
pgi {1.0}g6p = {1.0}$f6p
pgm {1.0}g6p = {1.0}g1p
tpp {1.0}t6p + {1.0}$h2o = {1.0}trh + {1.0}$pho
tps {1.0}g6p + {1.0}udg = {1.0}t6p + {1.0}$udp + {1.0}$h
ugp {1.0}g1p + {1.0}$utp + {1.0}$h = {1.0}udg + {1.0}$ppi
end reactions

begin rate equations
hxk = (cell*hxkshock^heat*hxkVmax*(-((adp*g6p[t])/hxkKeq) + atp*glc[t]))/((1 + adp/hxkKadp + atp/hxkKatp)*hxkKatp*hxkKglc*(1 + g6p[t]/hxkKg6p + glc[t]/hxkKglc + t6p[t]/hxkKit6p))
hxt = (cell*hxtshock^heat*hxtVmax*(glx - glc[t]))/(hxtKglc*(1 + (glx*hxtKi*glc[t])/hxtKglc^2 + (glx + glc[t])/hxtKglc))
nth = (cell*nthshock^heat*nthVmax*trh[t])/(nthKtrh*(1 + trh[t]/nthKtrh))
pgi = (cell*pgishock^heat*pgiVmax*(-(f6p/pgiKeq) + g6p[t]))/(pgiKg6p*(1 + f6p/pgiKf6p + g6p[t]/pgiKg6p))
pgm = (cell*pgmshock^heat*pgmVmax*(-(g1p[t]/pgmKeq) + g6p[t]))/(pgmKg6p*(1 + g1p[t]/pgmKg1p + g6p[t]/pgmKg6p))
tpp = (cell*tppshock^heat*tppVmax*t6p[t])/(tppKt6p*(1 + t6p[t]/tppKt6p))
tps = (cell*tpsactivity*tpsshock^heat*tpsVmax*g6p[t]*udg[t])/(tpsKg6p*tpsKudg*(1 + g6p[t]/tpsKg6p)*(1 + udg[t]/tpsKudg))
ugp = (cell*ugpshock^heat*ugpVmax*utp*g1p[t])/(ugpKg1p*ugpKutp*(ugpKiutp/ugpKutp + utp/ugpKutp + g1p[t]/ugpKg1p + (utp*g1p[t])/(ugpKg1p*ugpKutp) + (ugpKiutp*udg[t])/(ugpKiudg*ugpKutp) + (g1p[t]*udg[t])/(ugpKg1p*ugpKiudg)))
end rate equations

begin parameters
g1p0 = 0.1
g6p0 = 2.675
glc0 = 0.09765
heat = 0.0
t6p0 = 0.02
trh0 = 0.05
udg0 = 0.7
adp = 1.282
atp = 2.525
f6p = 0.625
glx = 100.0
h = 0.0
h2o = 0.0
pho = 0.0
ppi = 0.0
udp = 0.2815
utp = 0.6491
pgiVmax = 1071.0
pgiKg6p = 1.4
pgiKf6p = 0.29
pgiKeq = 0.3
pgishock = 1.0
hxtVmax = 97.24
hxtKglc = 1.1918
hxtKi = 0.91
hxtshock = 8.0
hxkVmax = 289.6
hxkKglc = 0.08
hxkKatp = 0.15
hxkKg6p = 30.0
hxkKadp = 0.23
hxkKeq = 2000.0
hxkKit6p = 0.04
hxkshock = 8.0
pgmVmax = 0.3545
pgmKg6p = 0.05
pgmKg1p = 0.023
pgmKeq = 0.1667
pgmshock = 16.0
tppVmax = 6.5
tppKt6p = 0.5
tppshock = 18.0
tpsVmax = 1.371
tpsKg6p = 3.8
tpsKudg = 0.886
tpsshock = 12.0
tpsactivity = 1.0
nthVmax = 15.2
nthKtrh = 2.99
nthshock = 6.0
ugpVmax = 36.82
ugpKutp = 0.11
ugpKiutp = 0.11
ugpKg1p = 0.32
ugpKiudg = 0.0035
ugpshock = 16.0
cell = 1.0
medium = 1.0
end parameters

begin initial conditions
g1p[0] = g1pi
g6p[0] = g6pi
glc[0] = glci
t6p[0] = t6pi
trh[0] = trhi
udg[0] = udgi
end initial conditions

begin initial values
g1pi = 0.1
g6pi = 2.675
glci = 0.09765
t6pi = 0.02
trhi = 0.05
udgi = 0.7
end initial values

begin assignment rules
udgchange := Log[udg[t]/udg0]/Log[10]
t6pchange := Log[t6p[t]/t6p0]/Log[10]
trhchange := Log[trh[t]/trh0]/Log[10]
g6pchange := Log[g6p[t]/g6p0]/Log[10]
g1pchange := Log[g1p[t]/g1p0]/Log[10]
glcchange := Log[glc[t]/glc0]/Log[10]
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
adp http://identifiers.org/obo.chebi/CHEBI:16761
atp http://identifiers.org/obo.chebi/CHEBI:15422
f6p http://identifiers.org/obo.chebi/CHEBI:16084
g1p http://identifiers.org/obo.chebi/CHEBI:16077
g6p http://identifiers.org/obo.chebi/CHEBI:17665
glc http://identifiers.org/obo.chebi/CHEBI:17925
glx http://identifiers.org/obo.chebi/CHEBI:17925
h http://identifiers.org/obo.chebi/CHEBI:24636
h2o http://identifiers.org/obo.chebi/CHEBI:15377
pho http://identifiers.org/obo.chebi/CHEBI:43474
ppi http://identifiers.org/obo.chebi/CHEBI:18361
t6p http://identifiers.org/obo.chebi/CHEBI:18283
trh http://identifiers.org/obo.chebi/CHEBI:16551
udg http://identifiers.org/obo.chebi/CHEBI:18066
udp http://identifiers.org/obo.chebi/CHEBI:17659
utp http://identifiers.org/obo.chebi/CHEBI:46398
end species annotations

begin units
end units