#arnold10

FIX: ATPc H Hc NADP NADPH O2 PiPi PiPic Succ UDPc UTPc

#reactions
AGPase:
 {1.0}HeP + {1.0}ATP = {1.0}ADPG + {1.0}PiPi
  chloroplast * function_28(Vm, PGA, Pi, G1P, ATP, ADPG, PiPi, q, Ks1, Ks2, Kp1, Kp2)
ATP_S:
 $pool = {1.0}ATP
  chloroplast * function_24(ADP, Pi, ATP, q12, Ks1, Ks2, Kp12, Vm)
F26BPc_P:
 {1.0}F26BPc = {1.0}HePc
  cytosol * function_33(Vm, F26BPc, Ks1, TPc, PGAc, Kr1, Pic, Kr3, HePc, Kr4)
F6P_TK:
 {1.0}HeP + {1.0}TP = {1.0}E4P + {1.0}PeP
  chloroplast * Transketolase(Vm, q, F6P, GAP, E4P, X5P, K1, K2, K1s2, K2s1, S7P, K2r1, K2s2, R5P)
F6Pc_K:
 {1.0}HePc = {1.0}F26BPc
  cytosol * function_32(Vm, F6Pc, Ks1, Pic, Kr1, TPc, PGAc, Kr2)
FBP_A:
 {2.0}TP = {1.0}FBP
  chloroplast * function_24(GAP, DHAP, FBP, q, Ks1, Ks2, Kp1, Vm)
FBPase:
 {1.0}FBP = {1.0}HeP
  chloroplast * function_25(Vm, FBP, F6P, Pi, q, Ks1, Kp1, Kp2)
FBPc_A:
 {2.0}TPc = {1.0}FBPc
  cytosol * function_24(GAPc, DHAPc, FBPc, q, Ks1, Ks2, Kp1, Vm)
FBPcase:
 {1.0}FBPc = {1.0}HePc
  cytosol * function_30(Vm, FBPc, F6Pc, Pic, q, Ks1, F26BPc, Kr1, Kp1, Kp2)
PGA_K:
 {1.0}PGA + {1.0}ATP = {1.0}TP
  chloroplast * function_23(Vm, PGA, ATP, GAP, ADP, Pi, q1, Ks1, Ks2, Kp1, Kp2, Kp3)
Ru5P_K:
 {1.0}PeP + {1.0}ATP = {1.0}RuBP
  chloroplast * function_27(Ru5P, ATP, RuBP, ADP, q, Ks1, Ks2, Kp1, Kp2, Vm)
RuBisCO_1:
 {1.0}RuBP = {1.0}ER
  chloroplast * (k1 * RuBP * E - k2 * ER)
RuBisCO_2_CO2:
 {1.0}ER = {1.0}EPP
  chloroplast * k1 * ER * CO2
RuBisCO_4:
 {1.0}EPP = {1.0}PGA + {1.0}EP
  chloroplast * (k1 * EPP - k2 * PGA * EP)
RuBisCO_5_EOP:
 {1.0}EOP = $pool
  chloroplast * function_34(W4)
RuBisCO_5_EP:
 {1.0}EP = {1.0}PGA
  chloroplast * (k1 * EP - k2 * PGA * E)
RuBisCO_6_O2:
 {1.0}ER + {1.0}O2 + {0.5}ATP = {1.0}EPG + {0.5}PGA
  chloroplast * k1 * ER * O2
RuBisCO_7:
 {1.0}EPG = {1.0}EP
  chloroplast * k1 * EPG
S7P_TK:
 {1.0}S7P + {1.0}TP = {2.0}PeP
  chloroplast * Transketolase(Vm, q, S7P, GAP, R5P, X5P, K1, K2, K1s2, K2s1, F6P, K2r1, K2s2, E4P)
SBP_A:
 {1.0}E4P + {1.0}TP = {1.0}SBP
  chloroplast * function_24(E4P, DHAP, SBP, q, Ks1, Ks2, Kp1, Vm)
SBPase:
 {1.0}SBP = {1.0}S7P
  chloroplast * function_26(Vm, SBP, S7P, Pi, q, Ks1, Kp1, Kp2)
StS:
 {1.0}ADPG = $pool
  chloroplast * function_29(Vm, ADPG, ADP, q, Ks1, Kp1)
SucPc_P:
 {1.0}SucPc = {1.0}Succ
  cytosol * function_25(Vm, SucPc, Succ, Pic, q, Ks1, Kp1, Kp2)
SucPc_S:
 {1.0}HePc + {1.0}UDPGc = {1.0}UDPc + {1.0}SucPc + {1.0}Hc
  cytosol * function_31(Vm, F6Pc, UDPGc, UDPc, SucPc, Hc, q, Ks1, Pic, Kr11, Ks2, Kp1, Kp2, Kr12)
TPT_PGA:
 {1.0}PGA = {1.0}PGAc
  TPTout(Vm, PGA, Ks, TP, Kr1, Pi, Kr2, PGAc, Kp, TPc, Kr3, Pic, Kr4)
TPT_TP:
 {1.0}TP = {1.0}TPc
  TPTout(Vm, TP, Ks, PGA, Kr1, Pi, Kr2, TPc, Kp, PGAc, Kr3, Pic, Kr4)
UGPase:
 {1.0}HePc + {1.0}UTPc = {1.0}UDPGc + {1.0}PiPic
  cytosol * function_27(G1Pc, UTPc, UDPGc, PiPic, q, Ks1, Ks2, Kp1, Kp2, Vm)


#parameters
ADT = 0.0015
ADTc = 0.001
Et = 0.0028030303030303
Kp12 = 224014.808032967
NADPT = 0.0005
PiT = 0.0284090909090909
PiTc = 0.0170454545454545
UDTc = 0.001
V28 = 7.386364e-05
q1 = 0.129053067280279
q12 = 2227862541257.35
ATPc = 0.00036
H = 0.0891250931577478
Hc = 0.158489318357816
NADP = 0.00029
NADPH = 0.00021
O2 = 0.00026
PiPi = 1e-06
PiPic = 4e-05
Succ = 0.0
UDPc = 0.00064
UTPc = 0.00036
FBPcase_Vm = 0.00113636
RuBisCO_1_k1 = 50000.0
RuBisCO_1_k2 = 0.9
RuBisCO_2_CO2_k1 = 300000.0
RuBisCO_4_k1 = 6.0
RuBisCO_4_k2 = 0.0
RuBisCO_5_EP_k1 = 6.0
RuBisCO_5_EP_k2 = 70000.0
RuBisCO_6_O2_k1 = 3030.3
RuBisCO_7_k1 = 3.0
PGA_K_Vm = 0.0170455
PGA_K_Ks1 = 0.0011122
PGA_K_Ks2 = 0.0003307
PGA_K_Kp1 = 0.00027035
PGA_K_Kp2 = 0.00053013
PGA_K_Kp3 = 0.0027397
FBP_A_q = 1.18815
FBP_A_Ks1 = 0.00027035
FBP_A_Ks2 = 0.00036393
FBP_A_Kp1 = 2.0129e-05
FBP_A_Vm = 0.022727
FBPase_Vm = 0.011364
FBPase_q = 0.77294
FBPase_Ks1 = 3.2842e-05
FBPase_Kp1 = 6.3429e-05
FBPase_Kp2 = 0.0017914
F6P_TK_Vm = 0.170455
F6P_TK_q = 0.99943
F6P_TK_K1 = 0.00061349
F6P_TK_K2 = 0.00011438
F6P_TK_K1s2 = 0.00027035
F6P_TK_K2s1 = 0.0005407
F6P_TK_K2r1 = 0.00017677
F6P_TK_K2s2 = 9.0464e-05
SBP_A_q = 1.18815
SBP_A_Ks1 = 0.00017677
SBP_A_Ks2 = 0.00036393
SBP_A_Kp1 = 2.0129e-05
SBP_A_Vm = 0.011364
SBPase_Vm = 0.00568182
SBPase_q = 0.77294
SBPase_Ks1 = 1.2713e-05
SBPase_Kp1 = 1.5597e-05
SBPase_Kp2 = 0.006744
S7P_TK_Vm = 0.0821023
S7P_TK_q = 0.99996
S7P_TK_K1 = 0.00061349
S7P_TK_K2 = 0.00011438
S7P_TK_K1s2 = 0.00027035
S7P_TK_K2s1 = 0.00017677
S7P_TK_K2r1 = 0.0005407
S7P_TK_K2s2 = 9.0464e-05
Ru5P_K_q = 1.05289
Ru5P_K_Ks1 = 3.63934e-05
Ru5P_K_Ks2 = 0.00055117
Ru5P_K_Kp1 = 9.95868e-05
Ru5P_K_Kp2 = 9.11825e-05
Ru5P_K_Vm = 0.568182
ATP_S_Ks1 = 0.00031808
ATP_S_Ks2 = 0.00031612
ATP_S_Vm = 0.0284091
AGPase_Vm = 0.00113636
AGPase_q = 0.11059
AGPase_Ks1 = 0.0010398
AGPase_Ks2 = 0.00011023
AGPase_Kp1 = 0.00053013
AGPase_Kp2 = 0.01951
StS_Vm = 0.00284091
StS_q = 1.00326
StS_Ks1 = 0.000212052
StS_Kp1 = 0.000636157
TPT_TP_Vm = 0.0568182
TPT_TP_Ks = 9.3583e-05
TPT_TP_Kr1 = 0.00089213
TPT_TP_Kr2 = 9.8597e-05
TPT_TP_Kp = 9.6372e-05
TPT_TP_Kr3 = 0.00054107
TPT_TP_Kr4 = 9.4837e-05
TPT_PGA_Vm = 0.0568182
TPT_PGA_Ks = 0.00089213
TPT_PGA_Kr1 = 9.3583e-05
TPT_PGA_Kr2 = 9.8597e-05
TPT_PGA_Kp = 0.00054107
TPT_PGA_Kr3 = 9.6372e-05
TPT_PGA_Kr4 = 9.4837e-05
FBPc_A_q = 1.00224
FBPc_A_Ks1 = 0.000278407
FBPc_A_Ks2 = 0.000374778
FBPc_A_Kp1 = 2.10226e-05
FBPc_A_Vm = 0.00568182
FBPcase_q = 0.792367
FBPcase_Ks1 = 2.2129e-05
FBPcase_Kr1 = 1.1065e-06
FBPcase_Kp1 = 6.5319e-05
FBPcase_Kp2 = 0.0018624
UGPase_q = 1.6219
UGPase_Ks1 = 3.2124e-05
UGPase_Ks2 = 0.0002364
UGPase_Kp1 = 0.00014393
UGPase_Kp2 = 0.0013192
UGPase_Vm = 0.00410568
SucPc_S_Vm = 7.38636e-05
SucPc_S_q = 1.00012
SucPc_S_Ks1 = 0.000278407
SucPc_S_Kr11 = 0.00920241
SucPc_S_Ks2 = 0.000110717
SucPc_S_Kp1 = 0.000642157
SucPc_S_Kp2 = 0.000374778
SucPc_S_Kr12 = 0.00164329
SucPc_P_Vm = 0.0010267
SucPc_P_q = 1.35286
SucPc_P_Ks1 = 5.354e-05
SucPc_P_Kp1 = 0.01
SucPc_P_Kp2 = 0.002191
F6Pc_K_Vm = 1.02614e-07
F6Pc_K_Ks1 = 0.001
F6Pc_K_Kr1 = 0.001
F6Pc_K_Kr2 = 0.0015
F26BPc_P_Vm = 2.05284e-10
F26BPc_P_Ks1 = 1e-09
F26BPc_P_Kr1 = 0.002
F26BPc_P_Kr3 = 0.001
F26BPc_P_Kr4 = 0.0001
chloroplast = 1.0
cytosol = 1.0


#initial values
ADPG = 6.295e-06
ATP = 0.00068
E4P = 5e-05
EOP = 0.0002966
EP = 0.0007045
EPG = 0.0003206
EPP = 0.0002711
ER = 0.001201
F26BPc = 7.8e-09
FBP = 0.00067
FBPc = 0.002
HeP = 0.0022
HePc = 0.0058
PGA = 0.0024
PGAc = 0.0
PeP = 0.00025
RuBP = 0.002
S7P = 0.002
SBP = 0.0003
SucPc = 0.0
TP = 0.0005
TPc = 0.0023
UDPGc = 0.00057


#assignment rules
!F Pic = PiTc - 2.0 * (FBPc + UTPc + ATPc + PiPic) - (PGAc + TPc + HePc + SucPc + UDPGc + UDPc + ADPc)
!F F6P = HeP * 0.999837 / (1.0 + 0.999837 + 0.999308)
!F Ru5P = PeP * 0.99974 / (1.0 + 0.99974 + 0.99974 / 1.000053)
!F GAP = TP / (1.0 + 1.0007329)
!F ADPc = ADTc - ATPc
!F W4 = 6.0 * EOP - 70000.0 * Pi * E
!F E = Et - (ER + EPP + EPG + EP + EOP)
!F CO2 = (1.200000e-05 / (0.38 + 0.015) + 8.0 * 4.0 * 0.00055 * 3030.3 * ER * O2 / 2.0) / (1.0 / (0.38 + 0.015) + 8.0 * 4.0 * 0.00055 * 300000.0 * ER)
!F DHAPc = 1.0007329 * TPc / (1.0 + 1.0007329)
!F GAPc = TPc / (1.0 + 1.0007329)
!F Pi = PiT - 2.0 * (EPP + EPG + RuBP + FBP + SBP + ATP + PiPi) - (EP + PGA + TP + HeP + E4P + S7P + PeP + ADP + ADPG)
!F ADP = ADT - ATP
!F R5P = PeP / (1.0 + 0.99974 + 0.99974 / 1.000053)
!F F6Pc = HePc * 0.999837 / (1.0 + 0.999837 + 0.999308)
!F G6Pc = HePc / (1.0 + 0.999837 + 0.999308)
!F DHAP = 1.0007329 * TP / (1.0 + 1.0007329)
!F G1Pc = HePc * 0.999308 / (1.0 + 0.999837 + 0.999308)
!F G1P = HeP * 0.999308 / (1.0 + 0.999837 + 0.999308)
!F G6P = HeP / (1.0 + 0.999837 + 0.999308)
!F X5P = PeP * (0.99974 / 1.000053) / (1.0 + 0.99974 + 0.99974 / 1.000053)