(* 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 = {
ADPc[t], ADPg[t], AMPc[t], AMPg[t], ATPc[t], ATPg[t], B13PGAg[t], DHAPc[t], DHAPg[t], Fru16BPg[t], Fru6Pg[t], GA3Pg[t], Glc6Pc[t], Glc6Pg[t], Glcc[t], Glcg[t], Gly3Pc[t], Gly3Pg[t], NADHg[t], NADPHc[t], NADPHg[t], NADPc[t], NADPg[t], NADg[t], P2GAc[t], P3GAc[t], P3GAg[t], P6GLc[t], P6GLg[t], P6Gc[t], P6Gg[t], PEPc[t], Pyrc[t], Rib5Pg[t], Rul5Pc[t], Rul5Pg[t], TS2c[t], TSH2c[t]
};

initialValues = {
ADPc[0] == 1.3165, ADPg[0] == 1.519, AMPc[0] == 2.2418, AMPg[0] == 4.2405, ATPc[0] == 0.3417, ATPg[0] == 0.2405, B13PGAg[0] == 0.5, DHAPc[0] == 2.23132912, DHAPg[0] == 8.483130623, Fru16BPg[0] == 10.0, Fru6Pg[0] == 0.5, GA3Pg[0] == 2.5, Glc6Pc[0] == 0.5, Glc6Pg[0] == 0.5, Glcc[0] == 0.1, Glcg[0] == 0.1, Gly3Pc[0] == 2.76867088, Gly3Pg[0] == 10.51686938, NADHg[0] == 2.0, NADPHc[0] == 3.9, NADPHg[0] == 3.9, NADPc[0] == 0.1, NADPg[0] == 0.1, NADg[0] == 2.0, P2GAc[0] == 0.1, P3GAc[0] == 0.1, P3GAg[0] == 0.1, P6GLc[0] == 0.0795278, P6GLg[0] == 0.0795278, P6Gc[0] == 0.0841958, P6Gg[0] == 0.0841958, PEPc[0] == 1.0, Pyrc[0] == 10.0, Rib5Pg[0] == 0.01, Rul5Pc[0] == 0.41282, Rul5Pg[0] == 0.41282, TS2c[0] == 0.37, TSH2c[0] == 0.01
};

rates = {
AKc, AKg, ALDg, ATPuc, ENOc, G3PDHg, G6PDHc, G6PDHg, G6PPc, GAPDHg, GDAg, GKg, GPOc, GlcTc, GlcTg, HXKc, HXKg, NADPHuc, NADPHug, P3GATg, P6GDHc, P6GDHg, PFKg, PGAMc, PGIg, PGKg, PGLc, PGLg, PPIc, PPIg, PYKc, PyrTc, RKg, TOXc, TPIg, TRc
};

rateEquations = {
AKc -> AKcAKck1*ADPc[t]^2 - AKcAKck2*AMPc[t]*ATPc[t], AKg -> AKgAKgk1*ADPg[t]^2 - AKgAKgk2*AMPg[t]*ATPg[t], ALDg -> (ALDgALDgVmax*Fru16BPg[t]*(1 - (DHAPg[t]*GA3Pg[t])/(ALDgALDgKeq*Fru16BPg[t])))/(ALDgALDgKmFru16BP*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP)*(1 + DHAPg[t]/ALDgALDgKmDHAP + Fru16BPg[t]/(ALDgALDgKmFru16BP*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP)) + GA3Pg[t]/ALDgALDgKmGA3P + (DHAPg[t]*GA3Pg[t])/(ALDgALDgKmDHAP*ALDgALDgKmGA3P) + (Fru16BPg[t]*GA3Pg[t])/(ALDgALDgKiGA3P*ALDgALDgKmFru16BP*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP)))), ATPuc -> (ATPucATPuck*ATPc[t])/ADPc[t], ENOc -> (ENOcENOcVmax*P2GAc[t]*(1 - PEPc[t]/(ENOcENOcKeq*P2GAc[t])))/(ENOcENOcKm2PGA*(1 + P2GAc[t]/ENOcENOcKm2PGA + PEPc[t]/ENOcENOcKmPEP)), G3PDHg -> (G3PDHgG3PDHgVmax*DHAPg[t]*(1 - (Gly3Pg[t]*NADg[t])/(G3PDHgG3PDHgKeq*DHAPg[t]*NADHg[t]))*NADHg[t])/(G3PDHgG3PDHgKmDHAP*G3PDHgG3PDHgKmNADH*(1 + DHAPg[t]/G3PDHgG3PDHgKmDHAP + Gly3Pg[t]/G3PDHgG3PDHgKmGly3P)*(1 + NADg[t]/G3PDHgG3PDHgKmNAD + NADHg[t]/G3PDHgG3PDHgKmNADH)), G6PDHc -> (G6PDHcG6PDHcVmax*Glc6Pc[t]*NADPc[t]*(1 - (NADPHc[t]*P6GLc[t])/(G6PDHcG6PDHcKeq*Glc6Pc[t]*NADPc[t])))/(G6PDHcG6PDHcKmGlc6P*G6PDHcG6PDHcKmNADP*(1 + NADPc[t]/G6PDHcG6PDHcKmNADP + NADPHc[t]/G6PDHcG6PDHcKmNADPH)*(1 + Glc6Pc[t]/G6PDHcG6PDHcKmGlc6P + P6GLc[t]/G6PDHcG6PDHcKm6PGL)), G6PDHg -> (G6PDHgG6PDHgVmax*Glc6Pg[t]*NADPg[t]*(1 - (NADPHg[t]*P6GLg[t])/(G6PDHgG6PDHgKeq*Glc6Pg[t]*NADPg[t])))/(G6PDHgG6PDHgKmGlc6P*G6PDHgG6PDHgKmNADP*(1 + NADPg[t]/G6PDHgG6PDHgKmNADP + NADPHg[t]/G6PDHgG6PDHgKmNADPH)*(1 + Glc6Pg[t]/G6PDHgG6PDHgKmGlc6P + P6GLg[t]/G6PDHgG6PDHgKm6PGL)), G6PPc -> (G6PPcG6PPcVmax*Glc6Pc[t]*(1 - Glcc[t]/(G6PPcG6PPcKeq*Glc6Pc[t])))/(G6PPcG6PPcKmGlc6P*(1 + Glc6Pc[t]/G6PPcG6PPcKmGlc6P + Glcc[t]/G6PPcG6PPcKmGlc)), GAPDHg -> (GAPDHgGAPDHgVmax*GA3Pg[t]*NADg[t]*(1 - (B13PGAg[t]*NADHg[t])/(GAPDHgGAPDHgKeq*GA3Pg[t]*NADg[t])))/(GAPDHgGAPDHgKmGA3P*GAPDHgGAPDHgKmNAD*(1 + B13PGAg[t]/GAPDHgGAPDHgKm13BPGA + GA3Pg[t]/GAPDHgGAPDHgKmGA3P)*(1 + NADg[t]/GAPDHgGAPDHgKmNAD + NADHg[t]/GAPDHgGAPDHgKmNADH)), GDAg -> -(GDAgGDAgk*DHAPg[t]*Gly3Pc[t]) + GDAgGDAgk*DHAPc[t]*Gly3Pg[t], GKg -> (GKgGKgVmax*ADPg[t]*(1 - (Glye*ATPg[t])/(GKgGKgKeq*ADPg[t]*Gly3Pg[t]))*Gly3Pg[t])/(GKgGKgKmADP*GKgGKgKmGly3P*(1 + ADPg[t]/GKgGKgKmADP + ATPg[t]/GKgGKgKmATP)*(1 + Glye/GKgGKgKmGly + Gly3Pg[t]/GKgGKgKmGly3P)), GPOc -> (GPOcGPOcVmax*Gly3Pc[t])/(GPOcGPOcKmGly3P*(1 + Gly3Pc[t]/GPOcGPOcKmGly3P)), GlcTc -> (GlcTcGlcTcVmax*(Glce - Glcc[t]))/(Glce + GlcTcGlcTcKmGlc + Glcc[t] + (Glce*GlcTcGlcTcalpha*Glcc[t])/GlcTcGlcTcKmGlc), GlcTg -> GlcTgGlcTgk1*Glcc[t] - GlcTgGlcTgk2*Glcg[t], HXKc -> (HXKcHXKcVmax*ATPc[t]*(1 - (ADPc[t]*Glc6Pc[t])/(HXKcHXKcKeq*ATPc[t]*Glcc[t]))*Glcc[t])/(HXKcHXKcKmATP*HXKcHXKcKmGlc*(1 + ADPc[t]/HXKcHXKcKmADP + ATPc[t]/HXKcHXKcKmATP)*(1 + Glc6Pc[t]/HXKcHXKcKmGlc6P + Glcc[t]/HXKcHXKcKmGlc)), HXKg -> (HXKgHXKgVmax*ATPg[t]*(1 - (ADPg[t]*Glc6Pg[t])/(HXKgHXKgKeq*ATPg[t]*Glcg[t]))*Glcg[t])/(HXKgHXKgKmATP*HXKgHXKgKmGlc*(1 + ADPg[t]/HXKgHXKgKmADP + ATPg[t]/HXKgHXKgKmATP)*(1 + Glc6Pg[t]/HXKgHXKgKmGlc6P + Glcg[t]/HXKgHXKgKmGlc)), NADPHuc -> NADPHucNADPHuck*NADPHc[t], NADPHug -> NADPHugNADPHugk*NADPHg[t], P3GATg -> -(P3GATgP3GATgk*P3GAc[t]) + P3GATgP3GATgk*P3GAg[t], P6GDHc -> (P6GDHcP6GDHcVmax*NADPc[t]*P6Gc[t]*(1 - (NADPHc[t]*Rul5Pc[t])/(P6GDHcP6GDHcKeq*NADPc[t]*P6Gc[t])))/(P6GDHcP6GDHcKm6PG*P6GDHcP6GDHcKmNADP*(1 + NADPc[t]/P6GDHcP6GDHcKmNADP + NADPHc[t]/P6GDHcP6GDHcKmNADPH)*(1 + P6Gc[t]/P6GDHcP6GDHcKm6PG + Rul5Pc[t]/P6GDHcP6GDHcKmRul5P)), P6GDHg -> (P6GDHgP6GDHgVmax*NADPg[t]*P6Gg[t]*(1 - (NADPHg[t]*Rul5Pg[t])/(P6GDHgP6GDHgKeq*NADPg[t]*P6Gg[t])))/(P6GDHgP6GDHgKm6PG*P6GDHgP6GDHgKmNADP*(1 + NADPg[t]/P6GDHgP6GDHgKmNADP + NADPHg[t]/P6GDHgP6GDHgKmNADPH)*(1 + P6Gg[t]/P6GDHgP6GDHgKm6PG + Rul5Pg[t]/P6GDHgP6GDHgKmRul5P)), PFKg -> (PFKgPFKgKi1*PFKgPFKgVmax*ATPg[t]*(1 - (ADPg[t]*Fru16BPg[t])/(PFKgPFKgKeq*ATPg[t]*Fru6Pg[t]))*Fru6Pg[t])/(PFKgPFKgKmATP*PFKgPFKgKmFru6P*(PFKgPFKgKi1 + Fru16BPg[t])*(PFKgPFKgKsATP/PFKgPFKgKmATP + ADPg[t]/PFKgPFKgKmADP + ATPg[t]/PFKgPFKgKmATP + (ADPg[t]*Fru16BPg[t])/(PFKgPFKgKi2*PFKgPFKgKmADP) + Fru6Pg[t]/PFKgPFKgKmFru6P + (ATPg[t]*Fru6Pg[t])/(PFKgPFKgKmATP*PFKgPFKgKmFru6P))), PGAMc -> (PGAMcPGAMcVmax*(1 - P2GAc[t]/(PGAMcPGAMcKeq*P3GAc[t]))*P3GAc[t])/(PGAMcPGAMcKm3PGA*(1 + P2GAc[t]/PGAMcPGAMcKm2PGA + P3GAc[t]/PGAMcPGAMcKm3PGA)), PGIg -> (PGIgPGIgVmax*(1 - Fru6Pg[t]/(PGIgPGIgKeq*Glc6Pg[t]))*Glc6Pg[t])/(PGIgPGIgKmGlc6P*(1 + Fru6Pg[t]/PGIgPGIgKmFru6P + Glc6Pg[t]/PGIgPGIgKmGlc6P + P6Gg[t]/PGIgPGIgKi6PG)), PGKg -> (PGKgPGKgVmax*ADPg[t]*B13PGAg[t]*(1 - (ATPg[t]*P3GAg[t])/(PGKgPGKgKeq*ADPg[t]*B13PGAg[t])))/(PGKgPGKgKm13BPGA*PGKgPGKgKmADP*(1 + ADPg[t]/PGKgPGKgKmADP + ATPg[t]/PGKgPGKgKmATP)*(1 + B13PGAg[t]/PGKgPGKgKm13BPGA + P3GAg[t]/PGKgPGKgKm3PGA)), PGLc -> cytosol*PGLcPGLck*(-(P6Gc[t]/PGLcPGLcKeq) + P6GLc[t]) + (PGLcPGLcVmax*(1 - P6Gc[t]/(PGLcPGLcKeq*P6GLc[t]))*P6GLc[t])/(PGLcPGLcKm6PGL*(1 + P6Gc[t]/PGLcPGLcKm6PG + P6GLc[t]/PGLcPGLcKm6PGL)), PGLg -> glycosome*PGLgPGLgk*(-(P6Gg[t]/PGLgPGLgKeq) + P6GLg[t]) + (PGLgPGLgVmax*(1 - P6Gg[t]/(PGLgPGLgKeq*P6GLg[t]))*P6GLg[t])/(PGLgPGLgKm6PGL*(1 + P6Gg[t]/PGLgPGLgKm6PG + P6GLg[t]/PGLgPGLgKm6PGL)), PPIc -> (PPIcPPIcVmax*(1 - Rib5Pc/(PPIcPPIcKeq*Rul5Pc[t]))*Rul5Pc[t])/(PPIcPPIcKmRul5P*(1 + Rib5Pc/PPIcPPIcKmRib5P + Rul5Pc[t]/PPIcPPIcKmRul5P)), PPIg -> (PPIgPPIgVmax*(1 - Rib5Pg[t]/(PPIgPPIgKeq*Rul5Pg[t]))*Rul5Pg[t])/(PPIgPPIgKmRul5P*(1 + Rib5Pg[t]/PPIgPPIgKmRib5P + Rul5Pg[t]/PPIgPPIgKmRul5P)), PYKc -> (PYKcPYKcVmax*ADPc[t]*(PEPc[t]/(PYKcPYKcKmPEP*(1 + ADPc[t]/PYKcPYKcKiADP + ATPc[t]/PYKcPYKcKiATP)))^PYKcPYKcn*(1 - (ATPc[t]*Pyrc[t])/(PYKcPYKcKeq*ADPc[t]*PEPc[t])))/(PYKcPYKcKmADP*(1 + ADPc[t]/PYKcPYKcKmADP + ATPc[t]/PYKcPYKcKmATP)*(1 + (PEPc[t]/(PYKcPYKcKmPEP*(1 + ADPc[t]/PYKcPYKcKiADP + ATPc[t]/PYKcPYKcKiATP)))^PYKcPYKcn + Pyrc[t]/PYKcPYKcKmPyr)), PyrTc -> (PyrTcPyrTcVmax*Pyrc[t])/(PyrTcPyrTcKmPyr*(1 + Pyrc[t]/PyrTcPyrTcKmPyr)), RKg -> (RKgRKgVmax*ADPg[t]*(1 - (Ribg*ATPg[t])/(RKgRKgKeq*ADPg[t]*Rib5Pg[t]))*Rib5Pg[t])/(RKgRKgKmADP*RKgRKgKmRib5P*(1 + ADPg[t]/RKgRKgKmADP + ATPg[t]/RKgRKgKmATP)*(1 + Ribg/RKgRKgKmRib + Rib5Pg[t]/RKgRKgKmRib5P)), TOXc -> TOXcTOXck*TSH2c[t], TPIg -> (TPIgTPIgVmax*DHAPg[t]*(1 - GA3Pg[t]/(TPIgTPIgKeq*DHAPg[t])))/(TPIgTPIgKmDHAP*(1 + DHAPg[t]/TPIgTPIgKmDHAP + GA3Pg[t]/TPIgTPIgKmGA3P)), TRc -> (TRcTRcVmax*NADPHc[t]*TS2c[t]*(1 - (NADPc[t]*TSH2c[t])/(TRcTRcKeq*NADPHc[t]*TS2c[t])))/(TRcTRcKmNADPH*TRcTRcKmTS2*(1 + NADPc[t]/TRcTRcKmNADP + NADPHc[t]/TRcTRcKmNADPH)*(1 + TS2c[t]/TRcTRcKmTS2 + TSH2c[t]/TRcTRcKmTSH2))
};

parameters = {
AKcAKck1 -> 480.0, AKcAKck2 -> 1000.0, AKgAKgk1 -> 480.0, AKgAKgk2 -> 1000.0, ALDgALDgKeq -> 0.084, ALDgALDgKiADP -> 1.51, ALDgALDgKiAMP -> 3.65, ALDgALDgKiATP -> 0.68, ALDgALDgKiGA3P -> 0.098, ALDgALDgKmDHAP -> 0.015, ALDgALDgKmFru16BP -> 0.009, ALDgALDgKmGA3P -> 0.067, ALDgALDgVmax -> 560.0, ATPucATPuck -> 50.0, CO2c -> 0.0, CO2g -> 0.0, ENOcENOcKeq -> 4.17, ENOcENOcKm2PGA -> 0.054, ENOcENOcKmPEP -> 0.24, ENOcENOcVmax -> 598.0, G3PDHgG3PDHgKeq -> 17085.0, G3PDHgG3PDHgKmDHAP -> 0.1, G3PDHgG3PDHgKmGly3P -> 2.0, G3PDHgG3PDHgKmNAD -> 0.4, G3PDHgG3PDHgKmNADH -> 0.01, G3PDHgG3PDHgVmax -> 465.0, G6PDHcG6PDHcKeq -> 5.02, G6PDHcG6PDHcKm6PGL -> 0.04, G6PDHcG6PDHcKmGlc6P -> 0.058, G6PDHcG6PDHcKmNADP -> 0.0094, G6PDHcG6PDHcKmNADPH -> 0.0001, G6PDHcG6PDHcVmax -> 8.4, G6PDHgG6PDHgKeq -> 5.02, G6PDHgG6PDHgKm6PGL -> 0.04, G6PDHgG6PDHgKmGlc6P -> 0.058, G6PDHgG6PDHgKmNADP -> 0.0094, G6PDHgG6PDHgKmNADPH -> 0.0001, G6PDHgG6PDHgVmax -> 8.4, G6PPcG6PPcKeq -> 263.0, G6PPcG6PPcKmGlc -> 5.6, G6PPcG6PPcKmGlc6P -> 5.6, G6PPcG6PPcVmax -> 28.0, GAPDHgGAPDHgKeq -> 0.066, GAPDHgGAPDHgKm13BPGA -> 0.1, GAPDHgGAPDHgKmGA3P -> 0.15, GAPDHgGAPDHgKmNAD -> 0.45, GAPDHgGAPDHgKmNADH -> 0.02, GAPDHgGAPDHgVmax -> 720.9, GDAgGDAgk -> 600.0, GKgGKgKeq -> 0.000837, GKgGKgKmADP -> 0.56, GKgGKgKmATP -> 0.24, GKgGKgKmGly -> 0.44, GKgGKgKmGly3P -> 3.83, GKgGKgVmax -> 200.0, GPOcGPOcKmGly3P -> 1.7, GPOcGPOcVmax -> 368.0, GlcTcGlcTcKmGlc -> 1.0, GlcTcGlcTcVmax -> 111.7, GlcTcGlcTcalpha -> 0.75, GlcTgGlcTgk1 -> 250000.0, GlcTgGlcTgk2 -> 250000.0, Glce -> 5.0, Glye -> 0.0, HXKcHXKcKeq -> 759.0, HXKcHXKcKmADP -> 0.126, HXKcHXKcKmATP -> 0.116, HXKcHXKcKmGlc -> 0.1, HXKcHXKcKmGlc6P -> 2.7, HXKcHXKcVmax -> 154.32, HXKgHXKgKeq -> 759.0, HXKgHXKgKmADP -> 0.126, HXKgHXKgKmATP -> 0.116, HXKgHXKgKmGlc -> 0.1, HXKgHXKgKmGlc6P -> 2.7, HXKgHXKgVmax -> 1774.68, NADPHucNADPHuck -> 2.0, NADPHugNADPHugk -> 2.0, P3GATgP3GATgk -> 250.0, P6GDHcP6GDHcKeq -> 47.0, P6GDHcP6GDHcKm6PG -> 0.0035, P6GDHcP6GDHcKmNADP -> 0.001, P6GDHcP6GDHcKmNADPH -> 0.0006, P6GDHcP6GDHcKmRul5P -> 0.03, P6GDHcP6GDHcVmax -> 10.6, P6GDHgP6GDHgKeq -> 47.0, P6GDHgP6GDHgKm6PG -> 0.0035, P6GDHgP6GDHgKmNADP -> 0.001, P6GDHgP6GDHgKmNADPH -> 0.0006, P6GDHgP6GDHgKmRul5P -> 0.03, P6GDHgP6GDHgVmax -> 10.6, PFKgPFKgKeq -> 1035.0, PFKgPFKgKi1 -> 15.8, PFKgPFKgKi2 -> 10.7, PFKgPFKgKmADP -> 1.0, PFKgPFKgKmATP -> 0.0648, PFKgPFKgKmFru6P -> 0.999, PFKgPFKgKsATP -> 0.0393, PFKgPFKgVmax -> 1708.0, PGAMcPGAMcKeq -> 0.17, PGAMcPGAMcKm2PGA -> 0.16, PGAMcPGAMcKm3PGA -> 0.15, PGAMcPGAMcVmax -> 225.0, PGIgPGIgKeq -> 0.457, PGIgPGIgKi6PG -> 0.14, PGIgPGIgKmFru6P -> 0.12, PGIgPGIgKmGlc6P -> 0.4, PGIgPGIgVmax -> 1305.0, PGKgPGKgKeq -> 3377.0, PGKgPGKgKm13BPGA -> 0.003, PGKgPGKgKm3PGA -> 1.62, PGKgPGKgKmADP -> 0.1, PGKgPGKgKmATP -> 0.29, PGKgPGKgVmax -> 2862.0, PGLcPGLcKeq -> 20000.0, PGLcPGLcKm6PG -> 0.05, PGLcPGLcKm6PGL -> 0.05, PGLcPGLcVmax -> 28.0, PGLcPGLck -> 0.055, PGLgPGLgKeq -> 20000.0, PGLgPGLgKm6PG -> 0.05, PGLgPGLgKm6PGL -> 0.05, PGLgPGLgVmax -> 5.0, PGLgPGLgk -> 0.055, PPIcPPIcKeq -> 5.6, PPIcPPIcKmRib5P -> 4.0, PPIcPPIcKmRul5P -> 1.4, PPIcPPIcVmax -> 72.0, PPIgPPIgKeq -> 5.6, PPIgPPIgKmRib5P -> 4.0, PPIgPPIgKmRul5P -> 1.4, PPIgPPIgVmax -> 72.0, PYKcPYKcKeq -> 10800.0, PYKcPYKcKiADP -> 0.64, PYKcPYKcKiATP -> 0.57, PYKcPYKcKmADP -> 0.114, PYKcPYKcKmATP -> 15.0, PYKcPYKcKmPEP -> 0.34, PYKcPYKcKmPyr -> 50.0, PYKcPYKcVmax -> 1020.0, PYKcPYKcn -> 2.5, Pig -> 0.0, PyrTcPyrTcKmPyr -> 1.96, PyrTcPyrTcVmax -> 230.0, Pyre -> 0.0, RKgRKgKeq -> 0.0036, RKgRKgKmADP -> 0.25, RKgRKgKmATP -> 0.24, RKgRKgKmRib -> 0.51, RKgRKgKmRib5P -> 0.39, RKgRKgVmax -> 10.0, Rib5Pc -> 0.01, Ribg -> 0.01, TOXcTOXck -> 2.0, TPIgTPIgKeq -> 0.046, TPIgTPIgKmDHAP -> 1.2, TPIgTPIgKmGA3P -> 0.25, TPIgTPIgVmax -> 999.3, TRcTRcKeq -> 434.0, TRcTRcKmNADP -> 0.081, TRcTRcKmNADPH -> 0.00077, TRcTRcKmTS2 -> 0.0069, TRcTRcKmTSH2 -> 0.0018, TRcTRcVmax -> 252.0, cytosol -> 5.4549, glycosome -> 0.2451, default -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "", "metabolite" -> "", "extent" -> ""}
};

(* Time evolution *)
odes = {
ADPc'[t] == 0.183321*ATPuc +0.183321*HXKc -0.366643*AKc -0.183321*PYKc, ADPg'[t] == 4.07997*PFKg +4.07997*HXKg -4.07997*GKg -8.15993*AKg -4.07997*PGKg -4.07997*RKg, AMPc'[t] == 0.183321*AKc , AMPg'[t] == 4.07997*AKg , ATPc'[t] == 0.183321*AKc +0.183321*PYKc -0.183321*ATPuc -0.183321*HXKc, ATPg'[t] == 4.07997*GKg +4.07997*AKg +4.07997*PGKg +4.07997*RKg -4.07997*PFKg -4.07997*HXKg, B13PGAg'[t] == 4.07997*GAPDHg -4.07997*PGKg, DHAPc'[t] == 0.183321*GPOc -0.183321*GDAg, DHAPg'[t] == 4.07997*GDAg +4.07997*ALDg -4.07997*G3PDHg -4.07997*TPIg, Fru16BPg'[t] == 4.07997*PFKg -4.07997*ALDg, Fru6Pg'[t] == 4.07997*PGIg -4.07997*PFKg, GA3Pg'[t] == 4.07997*TPIg +4.07997*ALDg -4.07997*GAPDHg, Glc6Pc'[t] == 0.183321*HXKc -0.183321*G6PDHc -0.183321*G6PPc, Glc6Pg'[t] == 4.07997*HXKg -4.07997*G6PDHg -4.07997*PGIg, Glcc'[t] == 0.183321*GlcTc +0.183321*G6PPc -0.183321*GlcTg -0.183321*HXKc, Glcg'[t] == 4.07997*GlcTg -4.07997*HXKg, Gly3Pc'[t] == 0.183321*GDAg -0.183321*GPOc, Gly3Pg'[t] == 4.07997*G3PDHg -4.07997*GKg -4.07997*GDAg, NADHg'[t] == 4.07997*GAPDHg -4.07997*G3PDHg, NADPHc'[t] == 0.183321*G6PDHc +0.183321*P6GDHc -0.183321*NADPHuc -0.183321*TRc, NADPHg'[t] == 4.07997*G6PDHg +4.07997*P6GDHg -4.07997*NADPHug, NADPc'[t] == 0.183321*NADPHuc +0.183321*TRc -0.183321*G6PDHc -0.183321*P6GDHc, NADPg'[t] == 4.07997*NADPHug -4.07997*G6PDHg -4.07997*P6GDHg, NADg'[t] == 4.07997*G3PDHg -4.07997*GAPDHg, P2GAc'[t] == 0.183321*PGAMc -0.183321*ENOc, P3GAc'[t] == 0.183321*P3GATg -0.183321*PGAMc, P3GAg'[t] == 4.07997*PGKg -4.07997*P3GATg, P6GLc'[t] == 0.183321*G6PDHc -0.183321*PGLc, P6GLg'[t] == 4.07997*G6PDHg -4.07997*PGLg, P6Gc'[t] == 0.183321*PGLc -0.183321*P6GDHc, P6Gg'[t] == 4.07997*PGLg -4.07997*P6GDHg, PEPc'[t] == 0.183321*ENOc -0.183321*PYKc, Pyrc'[t] == 0.183321*PYKc -0.183321*PyrTc, Rib5Pg'[t] == 4.07997*PPIg -4.07997*RKg, Rul5Pc'[t] == 0.183321*P6GDHc -0.183321*PPIc, Rul5Pg'[t] == 4.07997*P6GDHg -4.07997*PPIg, TS2c'[t] == 0.183321*TOXc -0.183321*TRc, TSH2c'[t] == 0.183321*TRc -0.183321*TOXc
};
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]}]