#goldbeter7



#reactions
AK_dissoc:
 $pool = {1.0}A + {1.0}K
  theta * cytosol * (d1 * AK - a1 * A * K)
A_degradation:
 {1.0}A = $pool
  theta * cytosol * vdAx * A / (KdAx + A)
A_translation:
 $pool = {1.0}A
  theta * cytosol * ksAx * MAx
BP_dephosphorylation:
 {1.0}Bp = {1.0}B
  theta * cytosol * VMP * Bp / (K2 + Bp)
B_degradation:
 {1.0}B = $pool
  theta * cytosol * kd1 * B
B_phosphorylation:
 {1.0}B = {1.0}Bp
  theta * cytosol * VMK * KID / (KID + D) * B / (K1 + B) * AK / Kt
B_shuttling:
 {1.0}BN = {1.0}B
  theta * cytosol * (kt4 * BN - kt3 * B)
B_synth:
 $pool = {1.0}B
  theta * cytosol * vsB
Bp_degradation:
 {1.0}Bp = $pool
  theta * cytosol * kd2 * Bp
Dusp_degradation:
 {1.0}Dusp = $pool
  eta * cytosol * vdDusp * Dusp / (KdDusp + Dusp)
Dusp_translation:
 $pool = {1.0}Dusp
  eta * cytosol * ksDusp * MDusp
Erk_activation:
 $pool = {1.0}ERKa
  eta * cytosol * VMaErk * Rasa / Rast * ERKi / (KaErk + ERKi)
Erk_inactivation:
 {1.0}ERKa = $pool
  eta * cytosol * kcDusp * Dusp * ERKa / (KdErk + ERKa)
F_degradation:
 {1.0}F = $pool
  epsilon * cytosol * vdF * F / (KdF + F)
F_translation:
 $pool = {1.0}F
  epsilon * cytosol * ksF * MF
MAx_degradation:
 {1.0}MAx = $pool
  theta * cytosol * vmd * MAx / (Kmd + MAx)
MAx_trans_BN:
 $pool = {1.0}MAx
  theta * cytosol * (vMB * pow(BN, n) / (pow(KaB, n) + pow(BN, n)))
MAx_trans_Xa:
 $pool = {1.0}MAx
  theta * cytosol * (vMXa * pow(Xa, m) / (pow(KaXa, m) + pow(Xa, m)))
MAx_trans_basal:
 $pool = {1.0}MAx
  theta * cytosol * v0
MDusp_degradation:
 {1.0}MDusp = $pool
  eta * cytosol * VMdMDusp * MDusp / (KdMDusp + MDusp)
MDusp_transkription:
 $pool = {1.0}MDusp
  eta * cytosol * VMsMDusp * pow(Xa, q) / (pow(KaMDusp, q) + pow(Xa, q))
MF_degradation:
 {1.0}MF = $pool
  epsilon * cytosol * vmF * MF / (KdMF + MF)
MF_transkription:
 $pool = {1.0}MF
  epsilon * cytosol * vsFK * pow(Nan, p) / (pow(KA, p) + pow(Nan, p))
N_activation:
 {1.0}N = {1.0}Na
  epsilon * cytosol * kc * N * pow(KIF, j) / (pow(KIF, j) + pow(F, j))
N_degradation:
 {1.0}N = $pool
  epsilon * cytosol * vdN * N / (KdN + N)
Na_degradation:
 {1.0}Na = $pool
  epsilon * cytosol * VdNa * Na / (KdNa + Na)
Na_transport:
 {1.0}Na = {1.0}Nan
  epsilon * cytosol * (kt1 * Na - kt2 * Nan)
Nan_degradation:
 {1.0}Nan = $pool
  epsilon * cytosol * VdNan * Nan / (KdNan + Nan)
Ras_activation:
 $pool = {1.0}Rasa
  eta * cytosol * VMaRas * pow(Fgf, r) / (pow(KaFgf, r) + pow(Fgf, r)) * Rasi / (KaRas + Rasi)
Ras_inactivation:
 {1.0}Rasa = $pool
  eta * cytosol * VMdRas * Rasa / (KdRas + Rasa)
X_activation:
 $pool = {1.0}Xa
  eta * cytosol * VMaX * ERKa / ERKt * Xi / (KaX + Xi)
X_inactivation:
 {1.0}Xa = $pool
  eta * cytosol * VMdX * Xa / (KdX + Xa)
n_synth:
 $pool = {1.0}N
  cytosol * epsilon * vsN


#parameters
K1 = 0.28
K2 = 0.03
KA = 0.05
KID = 0.5
KIF = 0.5
KIG1 = 2.5
KaB = 0.7
KaErk = 0.05
KaFgf = 0.5
KaMDusp = 0.5
KaRas = 0.103
KaX = 0.05
KaXa = 0.05
KdAx = 0.63
KdDusp = 0.5
KdErk = 0.05
KdF = 0.37
KdMDusp = 0.5
KdMF = 0.768
KdN = 1.4
KdNa = 0.001
KdNan = 0.001
KdRas = 0.1
KdX = 0.05
Kmd = 0.48
VMK = 5.08
VMP = 1.0
VMaErk = 3.3
VMaRas = 4.968
VMaX = 1.6
VMdMDusp = 0.5
VMdRas = 0.41
VMdX = 0.5
VMsMDusp = 0.9
VdNa = 0.01
VdNan = 0.1
a1 = 1.8
d1 = 0.1
epsilon = 0.3
eta = 0.3
j = 2.0
kc = 3.45
kcDusp = 1.35
kd1 = 0.0
kd2 = 7.062
ksAx = 0.02
ksDusp = 0.5
ksF = 0.3
kt1 = 0.1
kt2 = 0.1
kt3 = 0.7
kt4 = 1.5
m = 2.0
n = 2.0
p = 2.0
q = 2.0
r = 2.0
theta = 1.5
v0 = 0.06
vMB = 1.64
vMXa = 0.5
vdAx = 0.6
vdDusp = 2.0
vdF = 0.39
vdN = 2.82
vmF = 1.92
vmd = 0.8
vsB = 0.087
vsF = 3.0
vsN = 0.23
cytosol = 1.0


#initial values
A = 0.1
B = 0.1
BN = 0.001
Bp = 0.1
D = 2.0
Dusp = 0.1
ERKa = 0.2
ERKt = 2.0
F = 0.001
Fgf = 1.0
K = 3.0
Kt = 3.0
MAx = 0.1
MDusp = 0.1
MF = 0.1
N = 0.5
Na = 0.2
Nan = 0.0
Rasa = 0.5
Rast = 2.0
Xa = 0.1
Xt = 2.0


#assignment rules
!F Xi = Xt - Xa
!F ERKi = ERKt - ERKa
!F Rasi = Rast - Rasa
!F AK = Kt - K
!F vsFK = vsF * (KIG1 / (KIG1 + K))