begin name
mosca1
end name

begin reactions
reaction1 {1.0}$species9 = {1.0}species1
reaction10 {1.0}species17 + {1.0}species16 = {1.0}species5 + {1.0}species7
reaction11 {1.0}species15 + {1.0}species4 = {1.0}$species20 + {1.0}$species21
reaction12 {1.0}species22 = {1.0}species2
reaction13 {1.0}$species24 + {1.0}$species23 = {1.0}$species24 + {1.0}species22
reaction14 {1.0}species6 = {1.0}species16 + {1.0}species27
reaction15 {1.0}species16 = {1.0}species27
reaction16 {1.0}species16 + {1.0}species19 + {1.0}$species23 = {1.0}species12
reaction17 {1.0}species12 = {1.0}species28 + {1.0}species4
reaction18 {1.0}$species24 + {1.0}$species23 = {1.0}$species24 + {1.0}species22
reaction19 {1.0}species22 + {1.0}species4 = {1.0}$species24 + {2.0}$species23
reaction2 {1.0}species1 + {1.0}species4 = {1.0}species2
reaction20 {1.0}species5 + {1.0}species4 = {1.0}species6
reaction21 {1.0}species28 = {1.0}species29
reaction22 {1.0}species29 = {1.0}species30
reaction23 {1.0}species30 = {1.0}species31 + {1.0}species4
reaction24 {1.0}species31 = {1.0}$species32 + {1.0}species19
reaction25 {1.0}$EXTERNAL = {1.0}species4 + {1.0}$species20
reaction26 {1.0}$EXTERNAL = {1.0}species19
reaction27 {1.0}species11 = {1.0}$EXTERNAL
reaction28 {0.08}species31 + {0.2}$species34 + {1.0}$species23 = {0.24}$species33 + {1.0}species4
reaction29 {1.0}species4 = {1.0}$species23
reaction3 {1.0}species2 = {1.0}species5
reaction4 {1.0}species2 = {1.0}species8 + {1.0}species11
reaction5 {1.0}species8 = {1.0}species13 + {1.0}species11
reaction6 {1.0}species13 = {1.0}species14
reaction7 {1.0}species13 = {1.0}species15
reaction8 {1.0}species15 + {1.0}species14 = {1.0}species16 + {1.0}species17
reaction9 {1.0}species14 + {1.0}species7 = {1.0}species16 + {1.0}species5
end reactions

begin rate equations
reaction1 = compartment1*GLUT[parameter30, species9, species1[t], reaction1keq, reaction1KGlce, reaction1KGlc]
reaction10 = compartment1*TAL[reaction10Vmax, species17[t], species16[t], species7[t], species5[t], reaction10KeqTAL, reaction10K1, reaction10K2, reaction10K6, reaction10K3, reaction10K5, reaction10K4, reaction10K7]
reaction11 = compartment1*PRPPS[reaction11Vmax, species15[t], species4[t], species21, species20, reaction11Kapp, reaction11KATP, reaction11KR5P]
reaction12 = compartment1*PGLM[parameter55, species22[t], parameter57, parameter15, species2[t], parameter56]
reaction13 = compartment1*GPa[parameter58, species24, species23, parameter61, parameter62, parameter4, species22[t], reaction13KGLYb, parameter60, reaction13KiPi, parameter59]
reaction14 = compartment1*FBA[parameter37, species6[t], parameter70, parameter8, species27[t], species16[t], parameter68, parameter69]
reaction15 = compartment1*TPI[parameter38, species16[t], parameter72, parameter9, species27[t], parameter71]
reaction16 = compartment1*GAPDH[parameter39, species19[t], species16[t], species23, parameter76, parameter75, parameter77, parameter10, species12[t], species18, parameter73, parameter74]
reaction17 = compartment1*PGK[parameter40, species12[t], species3, reaction17alfa, parameter51, parameter52, parameter11, species28[t], species4[t], reaction17beta, parameter49, parameter50]
reaction18 = compartment1*GPb[parameter63, species24, species23, parameter66, reaction18KPi, parameter7, species22[t], parameter64, parameter65, parameter67, reaction18KiG1P, parameter27, reaction18nH, reaction18Kamp]
reaction19 = compartment1*GS[parameter41, reaction19Kf, parameter17, species22[t], species4[t], species24, species23, species3, reaction19Keq, reaction19Kr]
reaction2 = compartment1*HK[parameter31, reaction2Ka, reaction2Kb, species1[t], species4[t], species2[t], species3, reaction2Kapp, reaction2Kp, reaction2Kq]
reaction20 = compartment1*PFK[parameter42, parameter17, species4[t], reaction20Katp, reaction20beta, species26, reaction20alfa, reaction20Kf26bp, species5[t], reaction20Kf6p, reaction20L, species25, reaction20Kcit, reaction20Kiatp, species3, species6[t], reaction20Kadp, reaction20Kfbp, reaction20Kapp]
reaction21 = compartment1*function1[parameter78, species28[t], parameter80, parameter22, species29[t], parameter79]
reaction22 = compartment1*function2[parameter43, species29[t], parameter54, parameter24, species30[t], parameter53]
reaction23 = compartment1*function3[parameter44, species3, reaction23Kadp, species30[t], reaction23Kpep, reaction23L, species4[t], reaction23Kiatp, species6[t], reaction23Kfbp, species31[t], reaction23Katp, reaction23Kpyr, reaction23Kapp, parameter17]
reaction24 = compartment1*function4[parameter45, species18, species31[t], reaction24alfa, parameter85, parameter86, parameter26, species32, species19[t], reaction24beta, parameter83, parameter84]
reaction25 = compartment1*function5[reaction25Vf, species3, species4[t], species20, reaction25Keq]
reaction26 = compartment1*function6[parameter46, species18, species19[t], reaction26Keq]
reaction27 = compartment1*function7[parameter47, species11[t], species10, reaction27Keq]
reaction28 = compartment1*function8[parameter48, species31[t], reaction28y, species23, species3, species34, species4[t], species33, reaction28Keq]
reaction29 = compartment1*reaction29k1*species4[t]
reaction3 = compartment1*PGI[parameter32, species2[t], parameter82, parameter13, species5[t], parameter81, species7[t], reaction3Kery4p, species6[t], reaction3Kfbp, species8[t], reaction3Kpg]
reaction4 = compartment1*G6PDH[parameter33, reaction4KG6P, reaction4KNADP, species2[t], species10, species8[t], species11[t], reaction4Kapp, species4[t], reaction4KATP, reaction4KNADPH, species12[t], reaction4KPGA23]
reaction5 = compartment1*PGDH[parameter34, reaction5K6PG1, reaction5KNADP, species8[t], species10, species13[t], species11[t], reaction5Kapp, species12[t], reaction5KPGA23, species4[t], reaction5KATP, reaction5K6PG2, reaction5KNADPH]
reaction6 = compartment1*RUPE[reaction6Vmax, species13[t], species14[t], reaction6KeqRUPE, reaction6KRu5P, reaction6KX5P]
reaction7 = compartment1*R5PI[reaction7Vmax, species13[t], species15[t], reaction7KeqR5PI, reaction7KRu5P, reaction7KR5P]
reaction8 = compartment1*TKL[parameter35, species15[t], species14[t], species16[t], species17[t], reaction8KeqTKL, reaction8K1, reaction8K2, reaction8K6, reaction8K3, reaction8K5, reaction8K4, reaction8K7]
reaction9 = compartment1*TKL2[parameter36, species7[t], species14[t], species16[t], species5[t], reaction9KeqTKL2, reaction9K1, reaction9K2, reaction9K6, reaction9K3, reaction9K5, reaction9K4, reaction9K7]
end rate equations

begin parameters
parameter1 = 0.0114
parameter12 = 0.0018
parameter14 = 0.0556
parameter16 = 17.2
parameter17 = 1000.0
parameter18 = 0.381
parameter19 = 0.3574
parameter2 = "1.932e-05"
parameter20 = 11.369
parameter21 = 1.6491
parameter23 = 1.4127
parameter25 = 3452.5
parameter28 = 1.0
parameter29 = 1.0
parameter3 = 0.001345
parameter49 = 0.00013
parameter5 = 0.42
parameter50 = 0.00027
parameter51 = "7.9e-05"
parameter52 = "4e-05"
parameter53 = "6e-05"
parameter54 = "3.8e-05"
parameter55 = 7.364
parameter56 = "3e-05"
parameter57 = "6.3e-05"
parameter58 = 0.03347
parameter59 = 0.00015
parameter6 = 16.62
parameter60 = 0.0101
parameter61 = 0.0017
parameter62 = 0.004
parameter63 = 0.01049
parameter64 = 0.0044
parameter65 = 0.0015
parameter66 = 0.015
parameter67 = 0.0046
parameter68 = "8e-05"
parameter69 = 0.00016
parameter70 = "9e-06"
parameter71 = 0.0016
parameter72 = 0.00051
parameter73 = "2.2e-05"
parameter74 = "1e-05"
parameter75 = 0.00019
parameter76 = "9e-05"
parameter77 = 0.029
parameter78 = 154.0
parameter79 = 0.00012
parameter80 = 0.00019
parameter81 = "5e-05"
parameter82 = 0.0004
parameter83 = 0.0047
parameter84 = "7e-05"
parameter85 = "2e-06"
parameter86 = 0.0003
species20 = 0.00311
species21 = 0.001
species23 = 0.02
species24 = 0.208403745497308
species25 = 0.00108
species26 = "3.67e-06"
species32 = 0.0155
species33 = 0.0214
species34 = "6.5e-05"
species9 = 0.01
reaction23Kadp = 0.4
reaction23Kpep = 0.014
reaction23L = 1.0
reaction1keq = 1.0
reaction1KGlce = 0.01
reaction1KGlc = 0.0093
reaction2Ka = 0.0001
reaction2Kb = 0.0011
reaction2Kapp = 651.0
reaction2Kp = "2e-05"
reaction2Kq = 0.0035
reaction3Kery4p = "1e-06"
reaction3Kfbp = "6e-05"
reaction3Kpg = "1.5e-05"
reaction4KG6P = "6.67e-08"
reaction4KNADP = "3.67e-09"
reaction4Kapp = 2000.0
reaction4KATP = "7.49e-07"
reaction4KNADPH = "3.12e-09"
reaction4KPGA23 = "2.289e-06"
reaction5K6PG1 = "1e-08"
reaction5KNADP = "1.8e-08"
reaction5Kapp = 141.7
reaction5KPGA23 = "1.2e-07"
reaction5KATP = "1.54e-07"
reaction5K6PG2 = "5.8e-08"
reaction5KNADPH = "4.5e-09"
reaction6Vmax = 1.471
reaction6KeqRUPE = 2.7
reaction6KRu5P = "1.9e-07"
reaction6KX5P = "5e-07"
reaction7Vmax = 0.7646
reaction7KeqR5PI = 3.0
reaction7KRu5P = "7.8e-07"
reaction7KR5P = "2.2e-06"
reaction8KeqTKL = 2.08
reaction8K1 = "4.177e-07"
reaction8K2 = "3.055e-07"
reaction8K6 = 0.00774
reaction8K3 = "1.2432e-05"
reaction8K5 = 0.41139
reaction8K4 = "4.96e-09"
reaction8K7 = 48.8
reaction9KeqTKL2 = 29.7
reaction9K1 = "1.84e-09"
reaction9K2 = "3.055e-07"
reaction9K6 = 0.122
reaction9K3 = "5.48e-08"
reaction9K5 = 0.0287
reaction9K4 = "3e-10"
reaction9K7 = 0.215
reaction10Vmax = 58.27
reaction10KeqTAL = 2.703
reaction10K1 = "8.23e-09"
reaction10K2 = "4.765e-08"
reaction10K6 = 0.4653
reaction10K3 = "1.733e-07"
reaction10K5 = 0.8683
reaction10K4 = "6.095e-09"
reaction10K7 = 2.524
reaction11Vmax = 0.5104
reaction11Kapp = 100000.0
reaction11KATP = "3e-08"
reaction11KR5P = "5.7e-07"
reaction13KGLYb = 0.00015
reaction13KiPi = 0.0047
reaction17alfa = 1.0
reaction17beta = 1.0
reaction18KPi = 0.0002
reaction18KiG1P = 0.0074
reaction18nH = 1.75
reaction18Kamp = "1.9e-12"
reaction19Kf = 17400.0
reaction19Keq = 267100.0
reaction19Kr = 158.0
reaction20Katp = "2.1e-05"
reaction20beta = 0.98
reaction20alfa = 0.32
reaction20Kf26bp = "8.4e-07"
reaction20Kf6p = 1.0
reaction20L = 4.1
reaction20Kcit = 6.8
reaction20Kiatp = 20.0
reaction20Kadp = 5.0
reaction20Kfbp = 5.0
reaction20Kapp = 247.0
reaction23Kiatp = 2.5
reaction23Kfbp = 0.0004
reaction23Katp = 0.86
reaction23Kpyr = 10.0
reaction23Kapp = 195172.0
reaction24alfa = 1.0
reaction24beta = 1.0
reaction25Vf = 141.2
reaction25Keq = 2.26
reaction26Keq = 300.0
reaction27Keq = 0.2
reaction28y = 12.5
reaction28Keq = 1000000.0
reaction29k1 = 6210.0
compartment1 = 1.0
end parameters

begin initial conditions
species1[0] = species1i
species11[0] = species11i
species12[0] = species12i
species13[0] = species13i
species14[0] = species14i
species15[0] = species15i
species16[0] = species16i
species17[0] = species17i
species19[0] = species19i
species2[0] = species2i
species22[0] = species22i
species27[0] = species27i
species28[0] = species28i
species29[0] = species29i
species30[0] = species30i
species31[0] = species31i
species4[0] = species4i
species5[0] = species5i
species6[0] = species6i
species7[0] = species7i
species8[0] = species8i
end initial conditions

begin initial values
species1i = 0.000897
species11i = 1.87082e-05
species12i = 6.29e-05
species13i = 0.000143
species14i = 0.000242
species15i = 2.74e-05
species16i = 0.000153
species17i = 8.58e-05
species19i = 0.00134
species2i = 0.00109
species22i = 3.41e-05
species27i = 0.000553
species28i = 3.07e-05
species29i = 4.98e-06
species30i = 5.79e-05
species31i = 0.00183
species4i = 0.0087
species5i = 3.62e-05
species6i = 0.000367
species7i = 0.00093
species8i = 0.0001
end initial values

begin assignment rules
parameter34 := 31.02*parameter28
parameter48 := 9801000*parameter29
parameter11 := (parameter40*parameter49*parameter50)/(parameter20*parameter51*parameter52)
parameter42 := 107.6*parameter28
parameter41 := 32040*parameter28
parameter10 := (parameter39*parameter73*parameter74)/(parameter19*parameter75*parameter76*parameter77)
parameter39 := 109.1*parameter28
parameter9 := (parameter38*parameter71)/(parameter18*parameter72)
parameter38 := 5.976*parameter28
parameter8 := (parameter37*parameter68*parameter69)/(parameter12*parameter70)
parameter37 := 14.63*parameter28
parameter7 := (parameter63*parameter64*parameter65)/(parameter6*parameter66*parameter67)
parameter4 := (parameter58*parameter59*parameter60)/(parameter5*parameter61*parameter62)
parameter26 := (parameter45*parameter83*parameter84)/(parameter25*parameter85*parameter86)
parameter45 := 340.3*parameter28
parameter27 := species20
parameter13 := (parameter32*parameter81)/(parameter14*parameter82)
parameter32 := 7778*parameter28
parameter33 := 1.008*parameter28
parameter31 := 86.85*parameter28
parameter22 := (parameter78*parameter79)/(parameter21*parameter80)
parameter46 := 4982000*parameter28
species18 := parameter3 - species19[t]
species3 := parameter1 - species4[t]
species10 := parameter2 - species11[t]
parameter47 := 127800*parameter28
parameter15 := (parameter55*parameter56)/(parameter16*parameter57)
parameter24 := (parameter43*parameter53)/(parameter23*parameter54)
parameter43 := 160.9*parameter28
parameter44 := 27.81*parameter28
parameter35 := 1056*parameter28
parameter36 := 0.1761*parameter28
parameter40 := 73.41*parameter28
parameter30 := 23.03*parameter28
end assignment rules

begin function definitions
PRPPS[Vmax_,R5P_,ATP_,PRPP_,AMP_,Kapp_,KATP_,KR5P_] := ((-((AMP*PRPP)/Kapp) + ATP*R5P)*Vmax)/((ATP + KATP)*(KR5P + R5P))
TKL2[Vmax_,E4P_,X5P_,GAP_,F6P_,KeqTKL2_,K1_,K2_,K6_,K3_,K5_,K4_,K7_] := (Vmax*(-((F6P*GAP)/KeqTKL2) + E4P*X5P))/(F6P*K4 + GAP*(K3 + F6P*K5) + E4P*(K2 + F6P*K6) + (E4P + K1)*X5P + GAP*K7*X5P)
RUPE[Vmax_,RU5P_,X5P_,KeqRUPE_,KRu5P_,KX5P_] := (Vmax*(RU5P - X5P/KeqRUPE))/(RU5P + KRu5P*(1 + X5P/KX5P))
PGI[Vmf_,A_,Kg6p_,Vmr_,P_,Kf6p_,E4P_,Kery4p_,F16P_,Kfbp_,PGN_,Kpg_] := ((A*Vmf)/Kg6p - (P*Vmr)/Kf6p)/(1 + E4P/Kery4p + F16P/Kfbp + A/Kg6p + P/Kf6p + PGN/Kpg)
R5PI[Vmax_,RU5P_,R5P_,KeqR5PI_,KRu5P_,KR5P_] := ((-(R5P/KeqR5PI) + RU5P)*Vmax)/(KRu5P*(1 + R5P/KR5P) + RU5P)
GLUT[Vmaxf_,GLCe_,GLC_,keq_,KGlce_,KGlc_] := ((GLCe - GLC/keq)*Vmaxf)/(GLCe + (1 + GLC/KGlc)*KGlce)
PGDH[Vmax_,K6PG1_,KNADP_,PGN_,NADP_,RU5P_,NADPH_,Kapp_,BPG_,KPGA23_,ATP_,KATP_,K6PG2_,KNADPH_] := ((NADP*PGN - (NADPH*RU5P)/Kapp)*Vmax)/(K6PG1*KNADP*(ATP/KATP + (1 + NADP/KNADP)*(1 + BPG/KPGA23 + PGN/K6PG1) + (NADPH*(1 + PGN/K6PG2))/KNADPH))
G6PDH[Vmax_,KG6P_,KNADP_,G6P_,NADP_,PGN_,NADPH_,Kapp_,ATP_,KATP_,KNADPH_,BPG_,KPGA23_] := ((G6P*NADP - (NADPH*PGN)/Kapp)*Vmax)/(KG6P*KNADP*(1 + ATP/KATP + BPG/KPGA23 + ((1 + G6P/KG6P)*NADP)/KNADP + NADPH/KNADPH))
TAL[Vmax_,S7P_,GAP_,E4P_,F6P_,KeqTAL_,K1_,K2_,K6_,K3_,K5_,K4_,K7_] := ((-((E4P*F6P)/KeqTAL) + GAP*S7P)*Vmax)/(F6P*K4 + E4P*(K3 + F6P*K5) + GAP*(K2 + F6P*K6) + (GAP + K1)*S7P + E4P*K7*S7P)
PGLM[Vmaxf_,G1P_,KG1P_,Vmaxr_,G6P_,KG6P_] := ((G1P*Vmaxf)/KG1P - (G6P*Vmaxr)/KG6P)/(1 + G1P/KG1P + G6P/KG6P)
GPa[Vmaxf_,GLY_,Pi_,KiGLYf_,KPi_,Vmaxr_,G1P_,KGLYb_,KiG1P_,KiPi_,KiGLYb_] := ((GLY*Pi*Vmaxf)/(KiGLYf*KPi) - (G1P*GLY*Vmaxr)/(KGLYb*KiG1P))/(1 + G1P/KiG1P + GLY/KiGLYb + (G1P*GLY)/(KiG1P*KiGLYb) + GLY/KiGLYf + Pi/KiPi + (GLY*Pi)/(KiGLYf*KiPi))
FBA[Vmf_,A_,Kfbp_,Vmr_,P_,Q_,Kdhap_,Kg3p_] := ((A*Vmf)/Kfbp - (P*Q*Vmr)/(Kdhap*Kg3p))/(1 + A/Kfbp + P/Kdhap + Q/Kg3p + (P*Q)/(Kdhap*Kg3p))
TPI[Vf_,GAP_,Kms_,Vr_,DHAP_,Kmp_] := ((GAP*Vf)/Kms - (DHAP*Vr)/Kmp)/(1 + DHAP/Kmp + GAP/Kms)
PGK[Vmf_,A_,B_,alfa_,Ka_,Kb_,Vmr_,P_,Q_,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))
function3[Vm_,B_,Kadp_,A_,Kpep_,L_,Q_,Kiatp_,F16P_,Kfbp_,P_,Katp_,Kpyr_,Kapp_,a_] := (-((a^2*P*Q)/(Kapp*Katp*Kpyr*(1 + (a*P)/Kpyr + (a*Q)/Katp + (a^2*P*Q)/(Katp*Kpyr)))) + (a^2*A*B*(1 + (a*A)/Kpep)^3)/((1 + (a*B)/Kadp)*Kadp*Kpep*((1 + (a*A)/Kpep)^4 + (L*(1 + (a*Q)/Kiatp)^4)/(1 + (a*F16P)/Kfbp)^4)))*Vm
function2[Vmf_,PG2_,Kms_,Vmr_,PEP_,Kmp_] := ((PG2*Vmf)/Kms - (PEP*Vmr)/Kmp)/(1 + PEP/Kmp + PG2/Kms)
GS[Vmaxf_,Kf_,a_,G1P_,ATP_,GLY_,Pi_,ADP_,Keq_,Kr_] := (a^3*ATP*G1P*GLY*(1 - (a*ADP*Pi^2)/(ATP*G1P*Keq))*Vmaxf)/(Kf*(1 + (a^3*ATP*G1P*GLY)/Kf + (a^4*ADP*GLY*Pi^2)/Kr))
function1[Vmf_,PG3_,Kms_,Vmr_,PG2_,Kmp_] := ((PG3*Vmf)/Kms - (PG2*Vmr)/Kmp)/(1 + PG2/Kmp + PG3/Kms)
PFK[Vm_,a_,B_,Katp_,beta_,F26P_,alfa_,Kf26bp_,A_,Kf6p_,L_,CIT_,Kcit_,Kiatp_,Q_,P_,Kadp_,Kfbp_,Kapp_] := (a*B*(1 + (a*beta*F26P)/(alfa*Kf26bp))*((a*A*(1 + (a*F26P)/(alfa*Kf26bp))*(1 + (a*A*(1 + (a*F26P)/(alfa*Kf26bp)))/((1 + (a*F26P)/Kf26bp)*Kf6p))^3)/((1 + (a*F26P)/Kf26bp)*Kf6p*((1 + (a*A*(1 + (a*F26P)/(alfa*Kf26bp)))/((1 + (a*F26P)/Kf26bp)*Kf6p))^4 + ((1 + (a*CIT)/Kcit)^4*(1 + (a*B)/Kiatp)^4*L)/(1 + (a*F26P)/Kf26bp)^4)) - (a^2*P*Q)/(Kadp*Kapp*Kfbp*(1 + (a*P)/Kfbp + (a*Q)/Kadp + (a^2*P*Q)/(Kadp*Kfbp))))*Vm)/((1 + (a*B)/Katp)*Katp*(1 + (a*F26P)/(alfa*Kf26bp)))
function8[Vmf_,PYR_,y_,Pi_,ADP_,O2_,ATP_,CO2_,Keq_] := (ADP*O2^(5/(2*y))*Pi*PYR^y^(-1)*(1 - (ATP*CO2^(3/y))/(ADP*Keq*O2^(5/(2*y))*Pi*PYR^y^(-1)))*Vmf)/(-1 + (1 + ATP)*(1 + CO2)^(3/y) + (1 + ADP)*(1 + O2)^(5/(2*y))*(1 + Pi)*(1 + PYR)^y^(-1))
function7[Vf_,NADPH_,NADP_,Keq_] := ((1 - NADP/(Keq*NADPH))*NADPH*Vf)/(1 + NADP + NADPH)
function6[Vf_,NADH_,NAD_,Keq_] := ((1 - NAD/(Keq*NADH))*NADH*Vf)/(1 + NAD + NADH)
function5[Vf_,ADP_,ATP_,AMP_,Keq_] := (ADP^2*(1 - (AMP*ATP)/Keq)*Vf)/(-1 + (1 + ADP)^2 + (1 + AMP)*(1 + ATP))
function4[Vmf_,A_,B_,alfa_,Ka_,Kb_,Vmr_,P_,Q_,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))
GAPDH[Vmf_,A_,B_,C_,Knad_,Kg3p_,Kp_,Vmr_,P_,Q_,Kdpg_,Knadh_] := ((A*B*C*Vmf)/(Kg3p*Knad*Kp) - (P*Q*Vmr)/(Kdpg*Knadh))/(1 + A/Knad + (A*B)/(Kg3p*Knad) + (A*B*C)/(Kg3p*Knad*Kp) + Q/Knadh + (P*Q)/(Kdpg*Knadh))
GPb[Vmaxf_,GLY_,Pi_,KiGLYf_,KPi_,Vmaxr_,G1P_,KiGLYb_,KG1P_,KiPi_,KiG1P_,AMP_,nH_,Kamp_] := (AMP^nH*((GLY*Pi*Vmaxf)/(KiGLYf*KPi) - (G1P*GLY*Vmaxr)/(KG1P*KiGLYb)))/((1 + AMP^nH/Kamp)*Kamp*(1 + G1P/KiG1P + GLY/KiGLYb + (G1P*GLY)/(KG1P*KiGLYb) + GLY/KiGLYf + Pi/KiPi + (GLY*Pi)/(KiGLYf*KPi)))
HK[Vmf_,Ka_,Kb_,A_,B_,P_,Q_,Kapp_,Kp_,Kq_] := ((A*B - (P*Q)/Kapp)*Vmf)/(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)))
TKL[Vmax_,R5P_,X5P_,GAP_,S7P_,KeqTKL_,K1_,K2_,K6_,K3_,K5_,K4_,K7_] := (Vmax*(-((GAP*S7P)/KeqTKL) + R5P*X5P))/(K4*S7P + GAP*(K3 + K5*S7P) + R5P*(K2 + K6*S7P) + GAP*K7*X5P + (K1 + R5P)*X5P)
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units