#nakakuki2



#reactions
reaction_1:
 $pool = {1.0}x1
  compartment * (-x1 / tau1 + K1 * L / tau1)
reaction_2:
 $pool = {1.0}x2
  compartment * (-x2 / tau2 + K2 * L / tau2)
reaction_3:
 $pool = {1.0}ppERKn
  compartment * (k1 * ppERKc - k2 * ppERKn - k3 * DUSP * ppERKn)
reaction_4:
 $pool = {1.0}DUSP
  compartment * k * ppERKn
reaction_5:
 $pool = {1.0}pRSKn
  compartment * (k4 * ppERKn - k5 * pRSKn)
reaction_6:
 $pool = {1.0}cFOSp
  compartment * (pow(ppERKn * pRSKn, n) / (pow(k6, n) + pow(ppERKn * pRSKn, n)) - k7 * cFOSp)
reaction_7:
 $pool = {1.0}cFOSm
  compartment * (k7 * cFOSp - k8 * cFOSm)
reaction_8:
 $pool = {1.0}cFOS
  compartment * (k9 * cFOSm - k10 * cFOS - k11 * cFOS * ppERKc + k13 * pcFOS)
reaction_9:
 $pool = {1.0}pcFOS
  compartment * (k11 * cFOS * ppERKc - k12 * pcFOS - k13 * pcFOS)


#parameters
K1 = 1.09
K2 = 2.89
L = 1.0
k11 = 0.11
k13 = 0.06
k7 = 0.5
tau1 = 3.07
tau2 = 472.0
reaction_3_k1 = 15.0
reaction_3_k2 = 50.0
reaction_3_k3 = 14.0
reaction_4_k = 1.0
reaction_5_k4 = 0.1
reaction_5_k5 = 0.15
reaction_6_k6 = 0.13
reaction_6_n = 1.1
reaction_7_k8 = 0.08
reaction_8_k10 = 0.3
reaction_8_k9 = 0.3
reaction_9_k12 = 0.001
compartment = 1.0


#initial values
DUSP = 0.0
cFOS = 0.0
cFOSm = 0.0
cFOSp = 0.0
pRSKn = 0.0
pcFOS = 0.0
ppERKn = 0.0
x1 = 0.0
x2 = 0.0


#assignment rules
!F ppERKc = x1 - x2