(* Generated by JWS Online *)

(* This is an experimental feature of JWS Online. Please report any mistakes.*)

(* Note that the following notable SBML entities or features are not supported in notebook outputyet: *)
    (* Events *)
    (* Constraints *)
    (* Units and UnitDefinitions *)
    (* AlgebraicRules *)
    (* conversionFactors *)

variables = {
ACE[t], AMP[t], ATP[t], BPG[t], F16P[t], F6P[t], G3P[t], G6P[t], GLCi[t], NADH[t], P2G[t], P3G[t], PEP[t], PYR[t], TRIO[t]
};

initialValues = {
ACE[0] == 0.0320190578411128, AMP[0] == 0.169126793522056, ATP[0] == 2.88896546881255, BPG[0] == 0.000272378466415018, F16P[0] == 5.12303060629701, F6P[0] == 1.35316959717949, G3P[0] == 1.51464838664408, G6P[0] == 5.23300353550081, GLCi[0] == 0.0264527982731908, NADH[0] == 0.570243405517748, P2G[0] == 0.0219411145027765, P3G[0] == 0.182294273856502, PEP[0] == 0.0180777173496713, PYR[0] == 2.26719823411809, TRIO[0] == 2.65068578040418
};

rates = {
v\[LetterSpace]ADH, v\[LetterSpace]AK, v\[LetterSpace]ALD, v\[LetterSpace]ATP, v\[LetterSpace]ENO, v\[LetterSpace]G3PA, v\[LetterSpace]G3PDH, v\[LetterSpace]GAPDH, v\[LetterSpace]GLK, v\[LetterSpace]GLT, v\[LetterSpace]GLYCO, v\[LetterSpace]PDC, v\[LetterSpace]PFK, v\[LetterSpace]PGI, v\[LetterSpace]PGK, v\[LetterSpace]PGM, v\[LetterSpace]PYK, v\[LetterSpace]SUC, v\[LetterSpace]Treha
};

rateEquations = {
v\[LetterSpace]ADH -> -((VmADH*(ETOH*NAD - (ACE[t]*NADH[t])/KeqADH))/(KiADHNAD*KmADHETOH*(1 + (ETOH*KmADHNAD)/(KiADHNAD*KmADHETOH) + NAD/KiADHNAD + (ETOH*NAD)/(KiADHNAD*KmADHETOH) + (KmADHNADH*ACE[t])/(KiADHNADH*KmADHACE) + (ETOH*NAD*ACE[t])/(KiADHACE*KiADHNAD*KmADHETOH) + (KmADHNADH*NAD*ACE[t])/(KiADHNAD*KiADHNADH*KmADHACE) + NADH[t]/KiADHNADH + (ETOH*KmADHNAD*NADH[t])/(KiADHNAD*KiADHNADH*KmADHETOH) + (ACE[t]*NADH[t])/(KiADHNADH*KmADHACE) + (ETOH*ACE[t]*NADH[t])/(KiADHETOH*KiADHNADH*KmADHACE)))), v\[LetterSpace]AK -> 133.333*(ADP^2 - (AMP[t]*ATP[t])/KeqAK), v\[LetterSpace]ALD -> (VmALD*(F16P[t] - (KeqTPI*TRIO[t]^2)/(KeqALD*(1 + KeqTPI)^2)))/(KmALDF16P*(1 + F16P[t]/KmALDF16P + TRIO[t]/((1 + KeqTPI)*KmALDDHAP) + (KeqTPI*TRIO[t])/((1 + KeqTPI)*KmALDGAP) + (KeqTPI*F16P[t]*TRIO[t])/((1 + KeqTPI)*KmALDF16P*KmALDGAPi) + (KeqTPI*TRIO[t]^2)/((1 + KeqTPI)^2*KmALDDHAP*KmALDGAP))), v\[LetterSpace]ATP -> (KATPASE*ATP[t]^nATP)/(KmATP^nATP + ATP[t]^nATP), v\[LetterSpace]ENO -> (VmENO*(P2G[t] - PEP[t]/KeqENO))/(KmENOP2G*(1 + P2G[t]/KmENOP2G + PEP[t]/KmENOPEP)), v\[LetterSpace]G3PA -> (VmG3PA*G3P[t])/(KmG3PAG3P*(1 + Phi/KmG3PAPhi)*(1 + G3P[t]/KmG3PAG3P)), v\[LetterSpace]G3PDH -> (VmG3PDH*(-((NAD*G3P[t])/KeqG3PDH) + (NADH[t]*TRIO[t])/(1 + KeqTPI)))/(KmG3PDHDHAP*KmG3PDHNADH*(1 + ADP/KmG3PDHADP + ATP[t]/KmG3PDHATP + F16P[t]/KmG3PDHF16P)*(1 + NAD/KmG3PDHNAD + NADH[t]/KmG3PDHNADH)*(1 + G3P[t]/KmG3PDHG3P + TRIO[t]/((1 + KeqTPI)*KmG3PDHDHAP))), v\[LetterSpace]GAPDH -> (-((VmGAPDHf*BPG[t]*NADH[t])/(KeqGAPDH*KmGAPDHGAP*KmGAPDHNAD)) + (KeqTPI*NAD*VmGAPDHf*TRIO[t])/((1 + KeqTPI)*KmGAPDHGAP*KmGAPDHNAD))/((1 + NAD/KmGAPDHNAD + NADH[t]/KmGAPDHNADH)*(1 + BPG[t]/KmGAPDHBPG + (KeqTPI*TRIO[t])/((1 + KeqTPI)*KmGAPDHGAP))), v\[LetterSpace]GLK -> (VmGLK*(-((ADP*G6P[t])/KeqGLK) + ATP[t]*GLCi[t]))/(KmGLKATP*KmGLKGLCi*(1 + ADP/KmGLKADP + ATP[t]/KmGLKATP)*(1 + G6P[t]/KmGLKG6P + GLCi[t]/KmGLKGLCi)), v\[LetterSpace]GLT -> (VmGLT*(GLCo - GLCi[t]/KeqGLT))/(KmGLTGLCo*(1 + GLCo/KmGLTGLCo + GLCi[t]/KmGLTGLCi + (alpha*GLCo*GLCi[t])/(KmGLTGLCi*KmGLTGLCo))), v\[LetterSpace]GLYCO -> KGLYCOGEN*ATP[t]*G6P[t], v\[LetterSpace]PDC -> (VmPDC*PYR[t]^nPDC)/(KmPDCPYR^nPDC*(1 + PYR[t]^nPDC/KmPDCPYR^nPDC)), v\[LetterSpace]PFK -> (gR*VmPFK*ATP[t]*F6P[t]*(1 + ATP[t]/KmPFKATP + F6P[t]/KmPFKF6P + (gR*ATP[t]*F6P[t])/(KmPFKATP*KmPFKF6P)))/(KmPFKATP*KmPFKF6P*((L0*(1 + (CPFKAMP*AMP[t])/KPFKAMP)^2*(1 + (CiPFKATP*ATP[t])/KiPFKATP)^2*(1 + (CPFKATP*ATP[t])/KmPFKATP)^2*(1 + (CPFKF26BP*F26BP)/KPFKF26BP + (CPFKF16BP*F16P[t])/KPFKF16BP)^2)/((1 + AMP[t]/KPFKAMP)^2*(1 + ATP[t]/KiPFKATP)^2*(1 + F26BP/KPFKF26BP + F16P[t]/KPFKF16BP)^2) + (1 + ATP[t]/KmPFKATP + F6P[t]/KmPFKF6P + (gR*ATP[t]*F6P[t])/(KmPFKATP*KmPFKF6P))^2)), v\[LetterSpace]PGI -> (VmPGI*(-(F6P[t]/KeqPGI) + G6P[t]))/(KmPGIG6P*(1 + F6P[t]/KmPGIF6P + G6P[t]/KmPGIG6P)), v\[LetterSpace]PGK -> (VmPGK*(ADP*KeqPGK*BPG[t] - ATP[t]*P3G[t]))/(KmPGKATP*KmPGKP3G*(1 + ADP/KmPGKADP + ATP[t]/KmPGKATP)*(1 + BPG[t]/KmPGKBPG + P3G[t]/KmPGKP3G)), v\[LetterSpace]PGM -> (VmPGM*(-(P2G[t]/KeqPGM) + P3G[t]))/(KmPGMP3G*(1 + P2G[t]/KmPGMP2G + P3G[t]/KmPGMP3G)), v\[LetterSpace]PYK -> (VmPYK*(ADP*PEP[t] - (ATP[t]*PYR[t])/KeqPYK))/(KmPYKADP*KmPYKPEP*(1 + ADP/KmPYKADP + ATP[t]/KmPYKATP)*(1 + PEP[t]/KmPYKPEP + PYR[t]/KmPYKPYR)), v\[LetterSpace]SUC -> KSUCC*ACE[t], v\[LetterSpace]Treha -> KTREHALOSE*ATP[t]*G6P[t]
};

parameters = {
AXPsum -> 4.1, CPFKAMP -> 0.0845, CPFKATP -> 3.0, CPFKF16BP -> 0.397, CPFKF26BP -> 0.0174, CiPFKATP -> 100.0, ETOH -> 50.0, EXTERNAL -> 0.0, F26BP -> 0.02, GLCo -> 10.0, KATPASE -> 68.809631749526, KGLYCOGEN -> 1.68983106019077, KPFKAMP -> 0.0995, KPFKF16BP -> 0.111, KPFKF26BP -> 0.000682, KSUCC -> 19.6023811745556, KTREHALOSE -> 0.754128480342013, KeqADH -> 6.9*^-05, KeqAK -> 0.45, KeqALD -> 0.069, KeqENO -> 6.7, KeqG3PDH -> 4300.0, KeqGAPDH -> 0.00562639062770364, KeqGLK -> 3800.0, KeqGLT -> 1.0, KeqPGI -> 0.314, KeqPGK -> 3200.0, KeqPGM -> 0.19, KeqPYK -> 6500.0, KeqTPI -> 0.045, KiADHACE -> 1.1, KiADHETOH -> 90.0, KiADHNAD -> 0.92, KiADHNADH -> 0.031, KiPFKATP -> 0.65, KmADHACE -> 1.11, KmADHETOH -> 17.0, KmADHNAD -> 0.17, KmADHNADH -> 0.11, KmALDDHAP -> 2.4, KmALDF16P -> 0.3, KmALDGAP -> 2.0, KmALDGAPi -> 10.0, KmATP -> 0.263159506089159, KmENOP2G -> 0.04, KmENOPEP -> 0.5, KmG3PAG3P -> 4.21421727097407, KmG3PAPhi -> 0.798017147070638, KmG3PDHADP -> 1.60537927993171, KmG3PDHATP -> 0.568725492710765, KmG3PDHDHAP -> 0.4, KmG3PDHF16P -> 4.77110176421347, KmG3PDHG3P -> 1.088994646912, KmG3PDHNAD -> 0.93, KmG3PDHNADH -> 0.023, KmGAPDHBPG -> 0.0098, KmGAPDHGAP -> 0.21, KmGAPDHNAD -> 0.09, KmGAPDHNADH -> 0.06, KmGLKADP -> 0.23, KmGLKATP -> 0.15, KmGLKG6P -> 30.0, KmGLKGLCi -> 0.08, KmGLTGLCi -> 1.1918, KmGLTGLCo -> 1.1918, KmPDCPYR -> 4.33, KmPFKATP -> 0.71, KmPFKF6P -> 0.1, KmPGIF6P -> 0.3, KmPGIG6P -> 1.4, KmPGKADP -> 0.2, KmPGKATP -> 0.3, KmPGKBPG -> 0.003, KmPGKP3G -> 0.53, KmPGMP2G -> 0.08, KmPGMP3G -> 1.2, KmPYKADP -> 0.53, KmPYKATP -> 1.5, KmPYKPEP -> 0.14, KmPYKPYR -> 21.0, L0 -> 0.66, NADSUM -> 1.0, Phi -> 1.20470265921072, VmADH -> 834.684820558968, VmALD -> 232.895576807474, VmENO -> 462.02518044535, VmG3PA -> 542.31923582697, VmG3PDH -> 477.924460922096, VmGAPDHf -> 419.560093466267, VmGLK -> 329.958304618391, VmGLT -> 136.494419254188, VmPDC -> 223.978775020613, VmPFK -> 183.616866462272, VmPGI -> 467.754325004144, VmPGK -> 1419.49486841362, VmPGM -> 2762.32037661175, VmPYK -> 1751.9596196023, alpha -> 0.91, gR -> 5.12, nATP -> 1.0, nPDC -> 1.9, default\[LetterSpace]compartment -> 1.0
};

assignments = {
NAD -> NADSUM - NADH[t], ADP -> AXPsum - AMP[t] - ATP[t]
};

events = {

};

speciesAnnotations = {
ACE[t]->"http://identifiers.org/obo.chebi/CHEBI15343", ACE[t]->"http://identifiers.org/kegg.compound/C00084", AMP[t]->"http://identifiers.org/obo.chebi/CHEBI16027", AMP[t]->"http://identifiers.org/kegg.compound/C00020", ATP[t]->"http://identifiers.org/obo.chebi/CHEBI15422", ATP[t]->"http://identifiers.org/kegg.compound/C00002", BPG[t]->"http://identifiers.org/obo.chebi/28907", BPG[t]->"http://identifiers.org/kegg.compound/C00236", F16P[t]->"http://identifiers.org/obo.chebi/16905", F16P[t]->"http://identifiers.org/kegg.compound/C00354", F6P[t]->"http://identifiers.org/kegg.compound/C00085", F6P[t]->"http://identifiers.org/obo.chebi/15946", G3P[t]->"http://identifiers.org/kegg.compound/C00093", G3P[t]->"http://identifiers.org/obo.chebi/CHEBI15978", G6P[t]->"http://identifiers.org/kegg.compound/C00092", G6P[t]->"http://identifiers.org/obo.chebi/17665", GLCi[t]->"http://identifiers.org/kegg.compound/C00293", GLCi[t]->"http://identifiers.org/obo.chebi/CHEBI17234", NADH[t]->"http://identifiers.org/obo.chebi/16908", NADH[t]->"http://identifiers.org/kegg.compound/C00004", P2G[t]->"http://identifiers.org/kegg.compound/C00631", P2G[t]->"http://identifiers.org/obo.chebi/17835", P3G[t]->"http://identifiers.org/kegg.compound/C00197", P3G[t]->"http://identifiers.org/obo.chebi/17794", PEP[t]->"http://identifiers.org/obo.chebi/18021", PEP[t]->"http://identifiers.org/kegg.compound/C00074", PYR[t]->"http://identifiers.org/obo.chebi/15361", PYR[t]->"http://identifiers.org/kegg.compound/C00022", TRIO[t]->"http://identifiers.org/obo.chebi/CHEBI15978", TRIO[t]->"http://identifiers.org/obo.chebi/CHEBI17138", TRIO[t]->"http://identifiers.org/kegg.compound/C00093", TRIO[t]->"http://identifiers.org/kegg.compound/C00661"
};

reactionAnnotations = {
v\[LetterSpace]ADH->"http://identifiers.org/kegg.reaction/R00754", v\[LetterSpace]ALD->"http://identifiers.org/kegg.reaction/R01068", v\[LetterSpace]ATP->"http://identifiers.org/obo.go/0006200", v\[LetterSpace]ENO->"http://identifiers.org/kegg.reaction/R00658", v\[LetterSpace]GAPDH->"http://identifiers.org/kegg.reaction/R01061", v\[LetterSpace]GLK->"http://identifiers.org/kegg.reaction/R00299", v\[LetterSpace]GLT->"http://identifiers.org/obo.go/0015758", v\[LetterSpace]PDC->"http://identifiers.org/kegg.reaction/R00224", v\[LetterSpace]PFK->"http://identifiers.org/kegg.reaction/R04779", v\[LetterSpace]PGI->"http://identifiers.org/kegg.reaction/R02740", v\[LetterSpace]PGK->"http://identifiers.org/kegg.reaction/R01512", v\[LetterSpace]PGM->"http://identifiers.org/kegg.reaction/R01518"
};

units = {
{"time" -> "min", "metabolite" -> "mmol/L", "extent" -> "mM"}
};

(* Time evolution *)
odes = {
ACE'[t] == 1.0*v\[LetterSpace]PDC -1.0*v\[LetterSpace]ADH -2.0*v\[LetterSpace]SUC, AMP'[t] == 1.0*v\[LetterSpace]AK , ATP'[t] == 1.0*v\[LetterSpace]AK +1.0*v\[LetterSpace]PGK +1.0*v\[LetterSpace]PYK -1.0*v\[LetterSpace]ATP -1.0*v\[LetterSpace]GLK -1.0*v\[LetterSpace]GLYCO -1.0*v\[LetterSpace]PFK -4.0*v\[LetterSpace]SUC -1.0*v\[LetterSpace]Treha, BPG'[t] == 1.0*v\[LetterSpace]GAPDH -1.0*v\[LetterSpace]PGK, F16P'[t] == 1.0*v\[LetterSpace]PFK -1.0*v\[LetterSpace]ALD, F6P'[t] == 1.0*v\[LetterSpace]PGI -1.0*v\[LetterSpace]PFK, G3P'[t] == 1.0*v\[LetterSpace]G3PDH -1.0*v\[LetterSpace]G3PA, G6P'[t] == 1.0*v\[LetterSpace]GLK -1.0*v\[LetterSpace]GLYCO -1.0*v\[LetterSpace]PGI -2.0*v\[LetterSpace]Treha, GLCi'[t] == 1.0*v\[LetterSpace]GLT -1.0*v\[LetterSpace]GLK, NADH'[t] == 1.0*v\[LetterSpace]GAPDH +3.0*v\[LetterSpace]SUC -1.0*v\[LetterSpace]ADH -1.0*v\[LetterSpace]G3PDH, P2G'[t] == 1.0*v\[LetterSpace]PGM -1.0*v\[LetterSpace]ENO, P3G'[t] == 1.0*v\[LetterSpace]PGK -1.0*v\[LetterSpace]PGM, PEP'[t] == 1.0*v\[LetterSpace]ENO -1.0*v\[LetterSpace]PYK, PYR'[t] == 1.0*v\[LetterSpace]PYK -1.0*v\[LetterSpace]PDC, TRIO'[t] == 2.0*v\[LetterSpace]ALD -1.0*v\[LetterSpace]G3PDH -1.0*v\[LetterSpace]GAPDH
};
timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}];

(* Steady-state solution initialized with result of time evolution *)
findRootEquations = odes /.D[_[t],t]->0;
findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2];
steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100]
fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations

(* Plot the time evolution of the variables *)
plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]