#shen1



#reactions
v1:
 {1.0}NAD + {1.0}XYL = {1.0}XLAC
  (Eb*VmXDH*NAD*XYL)/(1000*KmXDHNAD*KmXDHXYL*(1 + (Ntot - NAD)/KmXDHNADH + NAD/KmXDHNAD)*(1 + XYL/KmXDHXYL + XLAC/KmXDHXLAC))
v2:
 {1.0}XLAC = {1.0}XA
  kXLA*XLAC
v3:
 {1.0}XLAC = {1.0}XA
  (VmXLA*Piecewise(0, time < 90, Ec/1000)*XLAC)/(KmXLAXLAC*(1 + XLAC/KmXLAXLAC + XA/KmXLAXA))
v4:
 {1.0}XA = {1.0}KDX
  (VmXAD*Piecewise(0, time < 180, Ed/1000)*XA)/(KmXADXA*(1 + (Ntot - NAD)/KiXADNADH)*(1 + KDX/KmXADKDX + XA/KmXADXA))
v5:
 {1.0}KDX = {1.0}KGSA
  (VmKDXD*KDX*Piecewise(0, time < 270, Ex/1000))/(KmKDXDKDX*(1 + KGSA/KmKDXDKGSA + KDX/KmKDXDKDX)*(1 + KG/KiKDXDakg + LAC/KiKDXDLAC + PYR/KiKDXDPYR + XA/KiKDXDxylonate))
v6:
 {1.0}KGSA + {1.0}NAD = {1.0}KG
  Piecewise(0, lt(time, 360), Ea / 1000) * (VmKGSADH * KGSA * NAD / (KmKGSADHKGSA * KmKGSADHNAD)) / ((1 + KGSA / KmKGSADHKGSA + KG / KmKGSADHKG) * (1 + NAD / KmKGSADHNAD + (Ntot - NAD) / KmKGSADHNADH) + KDX / KmKGSADHKDX * (1 + NAD / KmKGSADHNAD + (Ntot - NAD) / KmKGSADHNADH2))
v7:
 {1.0}PYR = {1.0}NAD + {1.0}LAC
  kLDH*(Ntot - NAD)*PYR


#parameters
Ea = 6.6
Eb = 2.5
Ec = 0.8
Ed = 6.5
Ex = 3.3
KiKDXDLAC = 27.9762
KiKDXDPYR = 17.9168
KiKDXDakg = 14.7919
KiKDXDxylonate = 18.3001
KiXADNADH = 10.4629
KmKDXDKDX = 0.207656
KmKDXDKGSA = 0.288658
KmKGSADHKDX = 0.2136
KmKGSADHKG = 0.279059
KmKGSADHKGSA = 0.0217501
KmKGSADHNAD = 0.596286
KmKGSADHNADH = 0.2674
KmKGSADHNADH2 = 0.0241
KmXADKDX = 0.862428
KmXADXA = 0.793955
KmXDHNAD = 0.161239
KmXDHNADH = 0.0297091
KmXDHXLAC = 0.535734
KmXDHXYL = 0.19851
KmXLAXA = 0.0381004
KmXLAXLAC = 0.445269
Ntot = 10.0
VmKDXD = 107.075
VmKGSADH = 49.3374
VmXAD = 42.359
VmXDH = 119.743
VmXLA = 944.406
kLDH = 0.0
kXLA = 0.00718272
default = 1.0


#initial values
KDX = 0.0
KG = 0.0
KGSA = 0.0
LAC = 0.0
NAD = 10.0
PYR = 0.0
XA = 0.0
XLAC = 0.0
XYL = 5.0


#assignment rules