#santolini1

FIX: NADPH NADPplus NO NO3 O2 citrulline

#reactions
r1:
 {1.0}FeIII + {1.0}NADPH = {1.0}FeII + {1.0}NADPplus
  cytosol * k1 * FeIII
r10:
 {1.0}FeII_NO + {1.0}O2 = {1.0}FeIII + {1.0}NO3
  cytosol * k10 * FeII_NO * O2
r2:
 {1.0}FeII + {1.0}O2 = {1.0}FeII_O2
  cytosol * k2 * FeII * O2
r3:
 {1.0}FeII_O2 = {1.0}FeIII_star
  cytosol * k3 * FeII_O2
r4:
 {1.0}FeIII_star + {0.5}NADPH = {1.0}FeII_star + {0.5}NADPplus
  cytosol * k4 * FeIII_star
r5:
 {1.0}FeII_star + {1.0}O2 = {1.0}FeII_star_O2
  cytosol * k5 * FeII_star * O2
r6:
 {1.0}FeII_star_O2 = {1.0}FeIII_NO + {1.0}citrulline
  cytosol * k6 * FeII_star_O2
r9:
 {1.0}FeII_NO = {1.0}FeII + {1.0}NO
  cytosol * k9 * FeII_NO
rF:
 {1.0}FeIII_NO = {1.0}FeIII + {1.0}NO
  cytosol * k7 * FeIII_NO
rG:
 {1.0}FeIII_NO + {0.5}NADPH = {1.0}FeII_NO + {0.5}NADPplus
  cytosol * k8 * FeIII_NO


#parameters
k1 = 2.6
k10 = 0.0013
k2 = 0.9
k3 = 26.0
k4 = 2.6
k5 = 0.9
k6 = 26.0
k7 = 5.0
k8 = 2.6
k9 = 0.0001
NADPH = 40.0
NADPplus = 0.0
NO = 0.0
NO3 = 0.0
O2 = 140.0
citrulline = 0.0
cytosol = 1.0


#initial values
FeII = 0.0
FeIII = 1.0
FeIII_NO = 0.0
FeIII_star = 0.0
FeII_NO = 0.0
FeII_O2 = 0.0
FeII_star = 0.0
FeII_star_O2 = 0.0


#assignment rules
!F FeIII_t = FeIII + FeIII_star