#arnold9

FIX: H NADP NADPH Pext

#reactions
ATP_S:
 $pool = {1.0}ATP
  chloroplast * function_6(Vm, ADP, Pi, K1, K2)
F6P_TK:
 {1.0}GAP + {1.0}F6P = {1.0}X5P + {1.0}E4P
  chloroplast * (k1 * GAP * F6P - k2 * X5P * E4P)
FBP_A:
 {1.0}DHAP + {1.0}GAP = {1.0}FBP
  chloroplast * (k1 * DHAP * GAP - k2 * FBP)
FBPase:
 {1.0}FBP = {1.0}F6P
  chloroplast * function_2(Vm, FBP, K, F6P, KR1, Pi, KR2)
GAP_DH:
 {1.0}DPGA + {1.0}NADPH + {1.0}H = {1.0}GAP + {1.0}NADP
  chloroplast * (k1 * DPGA * NADPH * H - k2 * GAP * NADP * Pi)
PGA_K:
 {1.0}PGA + {1.0}ATP = {1.0}DPGA
  chloroplast * (k1 * PGA * ATP - k2 * DPGA * ADP)
PG_I:
 {1.0}F6P = {1.0}G6P
  chloroplast * (k1 * F6P - k2 * G6P)
PG_M:
 {1.0}G6P = {1.0}G1P
  chloroplast * (k1 * G6P - k2 * G1P)
R5P_I:
 {1.0}R5P = {1.0}Ru5P
  chloroplast * (k1 * R5P - k2 * Ru5P)
Ru5P_E:
 {1.0}X5P = {1.0}Ru5P
  chloroplast * (k1 * X5P - k2 * Ru5P)
Ru5P_K:
 {1.0}Ru5P + {1.0}ATP = {1.0}RuBP
  chloroplast * function_4(Vm, Ru5P, ATP, K1, PGA, KR1, RuBP, KR2, Pi, KR3, ADP, KR41, K2, KR42)
RuBisCO:
 {1.0}RuBP = {2.0}PGA
  chloroplast * function_1(Vm, RuBP, K, PGA, KR1, FBP, KR2, SBP, KR3, Pi, KR4, NADPH, KR5)
S7P_TK:
 {1.0}GAP + {1.0}S7P = {1.0}X5P + {1.0}R5P
  chloroplast * (k1 * GAP * S7P - k2 * X5P * R5P)
SBP_A:
 {1.0}DHAP + {1.0}E4P = {1.0}SBP
  chloroplast * (k1 * DHAP * E4P - k2 * SBP)
SBPase:
 {1.0}SBP = {1.0}S7P
  chloroplast * function_3(Vm, SBP, K, Pi, KR1)
Starch_P:
 $pool = {1.0}G1P
  chloroplast * function_3(Vm, Pi, K, G1P, KR1)
Starch_S:
 {1.0}G1P + {1.0}ATP = $pool
  chloroplast * function_5(Vm, G1P, ATP, K1, K2, ADP, KR1, Pi, KA1, PGA, KA2, F6P, KA3, FBP)
TPT_DHAP:
 {1.0}DHAP = $pool
  chloroplast * function_7(Vm, DHAP, KA, Pext, K, Pi, KR1, PGA, KR2, GAP, KR3)
TPT_GAP:
 {1.0}GAP = $pool
  chloroplast * function_7(Vm, GAP, KA, Pext, K, Pi, KR1, PGA, KR2, DHAP, KR3)
TPT_PGA:
 {1.0}PGA = $pool
  chloroplast * function_7(Vm, PGA, KA, Pext, K, Pi, KR1, GAP, KR2, DHAP, KR3)
TP_I:
 {1.0}GAP = {1.0}DHAP
  chloroplast * (k1 * GAP - k2 * DHAP)


#parameters
H = 1.25892541179417e-05
NADP = 0.29
NADPH = 0.21
Pext = 0.5
RuBisCO_Vm = 340.0
RuBisCO_K = 0.02
RuBisCO_KR1 = 0.84
RuBisCO_KR2 = 0.04
RuBisCO_KR3 = 0.0075
RuBisCO_KR4 = 0.9
RuBisCO_KR5 = 0.07
PGA_K_k1 = 500000000.0
PGA_K_k2 = 1612900000000.0
GAP_DH_k1 = 500000000.0
GAP_DH_k2 = 31.25
TP_I_k1 = 500000000.0
TP_I_k2 = 22727000.0
FBP_A_k1 = 500000000.0
FBP_A_k2 = 70423000.0
FBPase_Vm = 200.0
FBPase_K = 0.03
FBPase_KR1 = 0.7
FBPase_KR2 = 12.0
F6P_TK_k1 = 500000000.0
F6P_TK_k2 = 5952400000.0
SBP_A_k1 = 500000000.0
SBP_A_k2 = 38461500.0
SBPase_Vm = 40.0
SBPase_K = 0.02
SBPase_KR1 = 12.0
S7P_TK_k1 = 500000000.0
S7P_TK_k2 = 588240000.0
R5P_I_k1 = 500000000.0
R5P_I_k2 = 1250000000.0
Ru5P_E_k1 = 500000000.0
Ru5P_E_k2 = 746269000.0
Ru5P_K_Vm = 1000.0
Ru5P_K_K1 = 0.05
Ru5P_K_KR1 = 2.0
Ru5P_K_KR2 = 0.7
Ru5P_K_KR3 = 4.0
Ru5P_K_KR41 = 2.5
Ru5P_K_K2 = 0.05
Ru5P_K_KR42 = 0.4
PG_I_k1 = 500000000.0
PG_I_k2 = 217400000.0
PG_M_k1 = 500000000.0
PG_M_k2 = 8621000000.0
Starch_S_Vm = 40.0
Starch_S_K1 = 0.08
Starch_S_K2 = 0.08
Starch_S_KR1 = 10.0
Starch_S_KA1 = 0.1
Starch_S_KA2 = 0.02
Starch_S_KA3 = 0.02
Starch_P_Vm = 40.0
Starch_P_K = 0.1
Starch_P_KR1 = 0.05
ATP_S_Vm = 350.0
ATP_S_K1 = 0.014
ATP_S_K2 = 0.3
TPT_PGA_Vm = 250.0
TPT_PGA_KA = 0.74
TPT_PGA_K = 0.25
TPT_PGA_KR1 = 0.63
TPT_PGA_KR2 = 0.075
TPT_PGA_KR3 = 0.077
TPT_GAP_Vm = 250.0
TPT_GAP_KA = 0.74
TPT_GAP_K = 0.075
TPT_GAP_KR1 = 0.63
TPT_GAP_KR2 = 0.25
TPT_GAP_KR3 = 0.077
TPT_DHAP_Vm = 250.0
TPT_DHAP_KA = 0.74
TPT_DHAP_K = 0.077
TPT_DHAP_KR1 = 0.63
TPT_DHAP_KR2 = 0.25
TPT_DHAP_KR3 = 0.075
chloroplast = 1.0
cytosol = 1.0


#initial values
ATP = 0.68
DHAP = 0.48
DPGA = 0.0011
E4P = 0.05
F6P = 0.640764257004719
FBP = 0.669999999999999
G1P = 0.0854779518844294
G6P = 1.47375779111085
GAP = 0.02
PGA = 2.4
R5P = 0.125186846038864
Ru5P = 0.0500747384155456
RuBP = 2.0
S7P = 2.0
SBP = 0.3
X5P = 0.0747384155455904


#assignment rules
!F ADP = 1.5 - ATP
!F Pi = 15.0 - 2.0 * (RuBP + DPGA + FBP + SBP) - (PGA + GAP + DHAP + F6P + E4P + S7P + X5P + R5P + Ru5P + G6P + G1P + ATP)