begin name
marinhernandez3
end name

begin reactions
AK {1.0}ATP + {1.0}AMP = {2.0}ADP
ALDO {1.0}FBP = {1.0}DHAP + {1.0}G3P
ATPases {1.0}ATP = {1.0}ADP + {1.0}$Pi
DHases {1.0}NADH = {1.0}NAD
ENO {1.0}2PG = {1.0}PEP
GAPDH {1.0}NAD + {1.0}G3P + {1.0}$Pi = {1.0}13BPG + {1.0}NADH
GLUT {1.0}$Gluout = {1.0}Gluin
Glycogendegradation {1.0}$glycogen + {1.0}$Pi = {1.0}G6P
Glycogensynthesis {1.0}G6P + {1.0}ATP = {1.0}$glycogen + {1.0}ADP + {2.0}$Pi
HK {1.0}Gluin + {1.0}ATP = {1.0}G6P + {1.0}ADP
HPI {1.0}G6P = {1.0}F6P
LDH {1.0}NADH + {1.0}Pyr = {1.0}Lacin + {1.0}NAD
MCT1 {1.0}Lacin = {1.0}$Lacout
MPM {1.0}Pyr + {13.0}ADP + {13.0}$Pi = {13.0}ATP
OxPhos {1.0}ADP + {1.0}$Pi = {1.0}ATP
PFK1 {1.0}F6P + {1.0}ATP = {1.0}FBP + {1.0}ADP
PGAM {1.0}3PG = {1.0}2PG
PGK {1.0}13BPG + {1.0}ADP = {1.0}3PG + {1.0}ATP
PPP {1.0}G6P = {1.0}$6PG
PYK {1.0}PEP + {1.0}ADP = {1.0}Pyr + {1.0}ATP
TK {1.0}$Xy5P + {1.0}$Ery4P = {1.0}G3P + {1.0}F6P
TPI {1.0}DHAP = {1.0}G3P
end reactions

begin rate equations
AK = compartment*(-(AKk2*ADP[t]^2) + AKk1*AMP[t]*ATP[t])
ALDO = compartment*ALDOrateequation[FBP[t], DHAP[t], G3P[t], ALDOVmf, ALDOKfbp, ALDOVmr, ALDOKdhap, ALDOKg3p]
ATPases = ATPasesk1*compartment*ATP[t]
DHases = compartment*(-(DHasesk2*NAD[t]) + DHasesk1*NADH[t])
ENO = compartment*ReversibleMichaelisMenten[2PG[t], PEP[t], ENOKms, ENOKmp, ENOVf, ENOVr]
GAPDH = compartment*GAPDH0[NAD[t], G3P[t], Pi, 13BPG[t], NADH[t], GAPDHVmf, GAPDHKnad, GAPDHKg3p, GAPDHKp, GAPDHVmr, GAPDHKdpg, GAPDHKnadh]
GLUT = compartment*GLUTisoformsrateequation[Gluout, Gluin[t], GLUTVmf, GLUTf1, GLUTKeq, GLUTKgluout, GLUTKgluin, GLUTf2, GLUTKeq1, GLUTKgluout1, GLUTKgluin1]
Glycogendegradation = compartment*Constantfluxirreversible[Glycogendegradationv]
Glycogensynthesis = compartment*Constantfluxirreversible[Glycogensynthesisv]
HK = compartment*HKisoformsrateequation[Gluin[t], ATP[t], G6P[t], ADP[t], HKVm, HKf1, HKKa, HKKb, HKKeq, HKKp, HKKq, HKf2, HKKa2]
HPI = compartment*HPIrateequation[G6P[t], F6P[t], Ery4P, FBP[t], 6PG, HPIVmf, HPIKg6p, HPIVmr, HPIKf6p, HPIKery4p, HPIKfbp, HPIKpg]
LDH = compartment*RandomBiBireversibleMichaelisMenten[NADH[t], Pyr[t], Lacin[t], NAD[t], LDHVmf, LDHalfa, LDHKa, LDHKb, LDHVmr, LDHbeta, LDHKp, LDHKq]
MCT1 = compartment*MCT10[Lacin[t], Lacout, MCT1Vmf, MCT1Keq, MCT1Klacin, MCT1Klacout]
MPM = compartment*Constantfluxirreversible[MPMv]
OxPhos = compartment*Constantfluxirreversible[OxPhosv]
PFK1 = compartment*PFK1rateequation[F6P[t], ATP[t], FBP[t], ADP[t], F26BP, Cit, PFK1Vm, PFK1Katp, PFK1beta, PFK1alfa, PFK1Kf26bp, PFK1Kf6p, PFK1L, PFK1Kcit, PFK1Kiatp, PFK1Kadp, PFK1Kfbp, PFK1Keq]
PGAM = compartment*ReversibleMichaelisMenten[3PG[t], 2PG[t], PGAMKms, PGAMKmp, PGAMVf, PGAMVr]
PGK = compartment*RandomBiBireversibleMichaelisMenten[13BPG[t], ADP[t], 3PG[t], ATP[t], PGKVmf, PGKalfa, PGKKa, PGKKb, PGKVmr, PGKbeta, PGKKp, PGKKq]
PPP = compartment*Constantfluxirreversible[PPPv]
PYK = compartment*PYKkinetics[PEP[t], ADP[t], Pyr[t], ATP[t], PYKVmax, PYKKpep, PYKKadp, PYKKeq, PYKKpyr, PYKKatp]
TK = compartment*Constantfluxirreversible[TKv]
TPI = compartment*ReversibleMichaelisMenten[DHAP[t], G3P[t], TPIKms, TPIKmp, TPIVf, TPIVr]
end rate equations

begin parameters
Cit = 1.7
Ery4P = 0.22
F26BP = 0.0042
Gluout = 5.0
Lacout = 1.9
Pi = 4.0
Xy5P = 0.016
6PG = 1.5
glycogen = 57.0
GLUTVmf = 0.023
GLUTf1 = 0.14
GLUTKeq = 1.0
GLUTKgluout = 1.8
GLUTKgluin = 10.0
GLUTf2 = 0.86
GLUTKeq1 = 1.0
GLUTKgluout1 = 9.3
GLUTKgluin1 = 10.0
HKVm = 0.036
HKf1 = 0.01
HKKa = 0.03
HKKb = 1.1
HKKeq = 651.0
HKKp = 0.02
HKKq = 3.5
HKf2 = 0.99
HKKa2 = 0.3
HPIVmf = 0.24
HPIKg6p = 0.4
HPIVmr = 0.54
HPIKf6p = 0.05
HPIKery4p = 0.001
HPIKfbp = 0.06
HPIKpg = 0.015
PFK1Vm = 0.022
PFK1Katp = 0.0292
PFK1beta = 1.18
PFK1alfa = 0.75
PFK1Kf26bp = 0.00099
PFK1Kf6p = 1.1
PFK1L = 6.6
PFK1Kcit = 6.7
PFK1Kiatp = 1.1
PFK1Kadp = 5.0
PFK1Kfbp = 5.0
PFK1Keq = 247.0
ALDOVmf = 0.08
ALDOKfbp = 0.009
ALDOVmr = 0.063
ALDOKdhap = 0.08
ALDOKg3p = 0.16
TPIKms = 1.6
TPIKmp = 0.51
TPIVf = 3.4
TPIVr = 28.0
GAPDHVmf = 0.28
GAPDHKnad = 0.09
GAPDHKg3p = 0.19
GAPDHKp = 11.0
GAPDHVmr = 0.35
GAPDHKdpg = 0.022
GAPDHKnadh = 0.01
PGKVmf = 8.7
PGKalfa = 1.0
PGKKa = 0.079
PGKKb = 0.04
PGKVmr = 2.5
PGKbeta = 1.0
PGKKp = 0.13
PGKKq = 0.27
PGAMKms = 0.19
PGAMKmp = 0.12
PGAMVf = 0.94
PGAMVr = 0.36
ENOKms = 0.038
ENOKmp = 0.06
ENOVf = 0.34
ENOVr = 0.38
PYKVmax = 0.072
PYKKpep = 0.05
PYKKadp = 0.4
PYKKeq = 195172.4
PYKKpyr = 10.0
PYKKatp = 0.86
LDHVmf = 0.44
LDHalfa = 1.0
LDHKa = 0.002
LDHKb = 0.3
LDHVmr = 0.07
LDHbeta = 1.0
LDHKp = 4.7
LDHKq = 0.07
Glycogendegradationv = 0.0045
ATPasesk1 = 0.00265
AKk1 = 1.0
AKk2 = 2.26
DHasesk1 = 250.0
DHasesk2 = 1.0
PPPv = "9.5e-05"
Glycogensynthesisv = 0.001
MPMv = 0.0001
TKv = "9.5e-05"
MCT1Vmf = 0.03
MCT1Keq = 1.0
MCT1Klacin = 8.5
MCT1Klacout = 0.5
OxPhosv = 0.01875
compartment = 1.0
end parameters

begin initial conditions
ADP[0] = ADPi
AMP[0] = AMPi
ATP[0] = ATPi
DHAP[0] = DHAPi
F6P[0] = F6Pi
FBP[0] = FBPi
G3P[0] = G3Pi
G6P[0] = G6Pi
Gluin[0] = Gluini
Lacin[0] = Lacini
NAD[0] = NADi
NADH[0] = NADHi
PEP[0] = PEPi
Pyr[0] = Pyri
13BPG[0] = 13BPGi
2PG[0] = 2PGi
3PG[0] = 3PGi
end initial conditions

begin initial values
ADPi = 2.1
AMPi = 1.6
ATPi = 11.0
DHAPi = 0.8
F6Pi = 2.7
FBPi = 0.5
G3Pi = 0.5
G6Pi = 4.0
Gluini = 0.001
Lacini = 33.0
NADi = 1.3
NADHi = 0.05
PEPi = 0.5
Pyri = 2.7
13BPGi = 0.001
2PGi = 0.001
3PGi = 0.001
end initial values

begin assignment rules
end assignment rules

begin function definitions
MCT10[Lacin_,Lacout_,Vmf_,Keq_,Klacin_,Klacout_] := ((Lacin - Lacout/Keq)*Vmf)/(Lacin + Klacin*(1 + Lacout/Klacout))
GLUTisoformsrateequation[Gluout_,Gluin_,Vmf_,f1_,Keq_,Kgluout_,Kgluin_,f2_,Keq1_,Kgluout1_,Kgluin1_] := ((f1*(Gluout - Gluin/Keq))/(Gluout + (1 + Gluin/Kgluin)*Kgluout) + (f2*(Gluout - Gluin/Keq1))/(Gluout + (1 + Gluin/Kgluin1)*Kgluout1))*Vmf
RandomBiBireversibleMichaelisMenten[A_,B_,P_,Q_,Vmf_,alfa_,Ka_,Kb_,Vmr_,beta_,Kp_,Kq_] := ((A*B*Vmf)/(alfa*Ka*Kb) - (P*Q*Vmr)/(beta*Kp*Kq))/(1 + A/Ka + B/Kb + (A*B)/(alfa*Ka*Kb) + P/Kp + Q/Kq + (P*Q)/(beta*Kp*Kq))
ReversibleMichaelisMenten[substrate_,product_,Kms_,Kmp_,Vf_,Vr_] := ((substrate*Vf)/Kms - (product*Vr)/Kmp)/(1 + product/Kmp + substrate/Kms)
Constantfluxirreversible[v_] := v
PYKkinetics[A_,B_,P_,Q_,Vmax_,Kpep_,Kadp_,Keq_,Kpyr_,Katp_] := (((A*B)/(Kadp*Kpep) - (P*Q)/(Kadp*Keq*Kpep))*Vmax)/((1 + A/Kpep + P/Kpyr)*(1 + B/Kadp + Q/Katp))
HPIrateequation[G6P_,F6P_,ERY4P_,FBP_,PG_,Vmf_,Kg6p_,Vmr_,Kf6p_,Kery4p_,Kfbp_,Kpg_] := ((G6P*Vmf)/Kg6p - (F6P*Vmr)/Kf6p)/(1 + ERY4P/Kery4p + F6P/Kf6p + FBP/Kfbp + G6P/Kg6p + PG/Kpg)
GAPDH0[NAD_,G3P_,P_,DPG_,NADH_,Vmf_,Knad_,Kg3p_,Kp_,Vmr_,Kdpg_,Knadh_] := ((G3P*NAD*P*Vmf)/(Kg3p*Knad*Kp) - (DPG*NADH*Vmr)/(Kdpg*Knadh))/(1 + NAD/Knad + (G3P*NAD)/(Kg3p*Knad) + NADH/Knadh + (DPG*NADH)/(Kdpg*Knadh) + (G3P*NAD*P)/(Kg3p*Knad*Kp))
PFK1rateequation[ATP_,F6P_,ADP_,FBP_,F26BP_,CIT_,Vm_,Katp_,beta_,alfa_,Kf26bp_,Kf6p_,L_,Kcit_,Kiatp_,Kadp_,Kfbp_,Keq_] := (ATP*(1 + (beta*F26BP)/(alfa*Kf26bp))*(-((ADP*FBP)/(Kadp*Keq*(1 + ADP/Kadp + FBP/Kfbp + (ADP*FBP)/(Kadp*Kfbp))*Kfbp)) + (F6P*(1 + F26BP/(alfa*Kf26bp))*(1 + (F6P*(1 + F26BP/(alfa*Kf26bp)))/((1 + F26BP/Kf26bp)*Kf6p))^3)/((1 + F26BP/Kf26bp)*Kf6p*((1 + (F6P*(1 + F26BP/(alfa*Kf26bp)))/((1 + F26BP/Kf26bp)*Kf6p))^4 + ((1 + CIT/Kcit)^4*(1 + ATP/Kiatp)^4*L)/(1 + F26BP/Kf26bp)^4)))*Vm)/((1 + ATP/Katp)*Katp*(1 + F26BP/(alfa*Kf26bp)))
ALDOrateequation[FBP_,DHAP_,G3P_,Vmf_,Kfbp_,Vmr_,Kdhap_,Kg3p_] := ((FBP*Vmf)/Kfbp - (DHAP*G3P*Vmr)/(Kdhap*Kg3p))/(1 + DHAP/Kdhap + FBP/Kfbp + G3P/Kg3p + (DHAP*G3P)/(Kdhap*Kg3p))
HKisoformsrateequation[A_,B_,P_,Q_,Vm_,f1_,Ka_,Kb_,Keq_,Kp_,Kq_,f2_,Ka2_] := ((f1*(A*B - (P*Q)/Keq))/(Ka*Kb*(1 + A/Ka + B/Kb + (A*B)/(Ka*Kb) + P/Kp + (B*P)/(Kb*Kp) + Q/Kq + (A*Q)/(Ka*Kq) + (P*Q)/(Kp*Kq))) + (f2*(A*B - (P*Q)/Keq))/(Ka2*Kb*(1 + A/Ka2 + B/Kb + (A*B)/(Ka2*Kb) + P/Kp + (B*P)/(Kb*Kp) + Q/Kq + (A*Q)/(Ka2*Kq) + (P*Q)/(Kp*Kq))))*Vm
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units