begin name
goldbeter7
end name

begin reactions
AKdissoc {1.0}$EXTERNAL = {1.0}A + {1.0}K
Adegradation {1.0}A = {1.0}$EXTERNAL
Atranslation {1.0}$EXTERNAL = {1.0}A
BPdephosphorylation {1.0}Bp = {1.0}B
Bdegradation {1.0}B = {1.0}$EXTERNAL
Bphosphorylation {1.0}B = {1.0}Bp
Bshuttling {1.0}BN = {1.0}B
Bsynth {1.0}$EXTERNAL = {1.0}B
Bpdegradation {1.0}Bp = {1.0}$EXTERNAL
Duspdegradation {1.0}Dusp = {1.0}$EXTERNAL
Dusptranslation {1.0}$EXTERNAL = {1.0}Dusp
Erkactivation {1.0}$EXTERNAL = {1.0}ERKa
Erkinactivation {1.0}ERKa = {1.0}$EXTERNAL
Fdegradation {1.0}F = {1.0}$EXTERNAL
Ftranslation {1.0}$EXTERNAL = {1.0}F
MAxdegradation {1.0}MAx = {1.0}$EXTERNAL
MAxtransBN {1.0}$EXTERNAL = {1.0}MAx
MAxtransXa {1.0}$EXTERNAL = {1.0}MAx
MAxtransbasal {1.0}$EXTERNAL = {1.0}MAx
MDuspdegradation {1.0}MDusp = {1.0}$EXTERNAL
MDusptranskription {1.0}$EXTERNAL = {1.0}MDusp
MFdegradation {1.0}MF = {1.0}$EXTERNAL
MFtranskription {1.0}$EXTERNAL = {1.0}MF
Nactivation {1.0}N = {1.0}Na
Ndegradation {1.0}N = {1.0}$EXTERNAL
Nadegradation {1.0}Na = {1.0}$EXTERNAL
Natransport {1.0}Na = {1.0}Nan
Nandegradation {1.0}Nan = {1.0}$EXTERNAL
Rasactivation {1.0}$EXTERNAL = {1.0}Rasa
Rasinactivation {1.0}Rasa = {1.0}$EXTERNAL
Xactivation {1.0}$EXTERNAL = {1.0}Xa
Xinactivation {1.0}Xa = {1.0}$EXTERNAL
nsynth {1.0}$EXTERNAL = {1.0}N
end reactions

begin rate equations
AKdissoc = cytosol*theta*(AK*d1 - a1*A[t]*K[t])
Adegradation = (cytosol*theta*vdAx*A[t])/(KdAx + A[t])
Atranslation = cytosol*ksAx*theta*MAx[t]
BPdephosphorylation = (cytosol*theta*VMP*Bp[t])/(K2 + Bp[t])
Bdegradation = cytosol*kd1*theta*B[t]
Bphosphorylation = (AK*cytosol*KID*theta*VMK*B[t])/((K1 + B[t])*(KID + D[t])*Kt[t])
Bshuttling = cytosol*theta*(-(kt3*B[t]) + kt4*BN[t])
Bsynth = cytosol*theta*vsB
Bpdegradation = cytosol*kd2*theta*Bp[t]
Duspdegradation = (cytosol*eta*vdDusp*Dusp[t])/(KdDusp + Dusp[t])
Dusptranslation = cytosol*eta*ksDusp*MDusp[t]
Erkactivation = (cytosol*ERKi*eta*VMaErk*Rasa[t])/((ERKi + KaErk)*Rast[t])
Erkinactivation = (cytosol*eta*kcDusp*Dusp[t]*ERKa[t])/(KdErk + ERKa[t])
Fdegradation = (cytosol*epsilon*vdF*F[t])/(KdF + F[t])
Ftranslation = cytosol*epsilon*ksF*MF[t]
MAxdegradation = (cytosol*theta*vmd*MAx[t])/(Kmd + MAx[t])
MAxtransBN = (cytosol*theta*vMB*BN[t]^n)/(KaB^n + BN[t]^n)
MAxtransXa = (cytosol*theta*vMXa*Xa[t]^m)/(KaXa^m + Xa[t]^m)
MAxtransbasal = cytosol*theta*v0
MDuspdegradation = (cytosol*eta*VMdMDusp*MDusp[t])/(KdMDusp + MDusp[t])
MDusptranskription = (cytosol*eta*VMsMDusp*Xa[t]^q)/(KaMDusp^q + Xa[t]^q)
MFdegradation = (cytosol*epsilon*vmF*MF[t])/(KdMF + MF[t])
MFtranskription = (cytosol*epsilon*vsFK*Nan[t]^p)/(KA^p + Nan[t]^p)
Nactivation = (cytosol*epsilon*kc*KIF^j*N[t])/(KIF^j + F[t]^j)
Ndegradation = (cytosol*epsilon*vdN*N[t])/(KdN + N[t])
Nadegradation = (cytosol*epsilon*VdNa*Na[t])/(KdNa + Na[t])
Natransport = cytosol*epsilon*(kt1*Na[t] - kt2*Nan[t])
Nandegradation = (cytosol*epsilon*VdNan*Nan[t])/(KdNan + Nan[t])
Rasactivation = (cytosol*eta*Rasi*VMaRas*Fgf[t]^r)/((KaRas + Rasi)*(KaFgf^r + Fgf[t]^r))
Rasinactivation = (cytosol*eta*VMdRas*Rasa[t])/(KdRas + Rasa[t])
Xactivation = (cytosol*eta*VMaX*Xi*ERKa[t])/((KaX + Xi)*ERKt[t])
Xinactivation = (cytosol*eta*VMdX*Xa[t])/(KdX + Xa[t])
nsynth = cytosol*epsilon*vsN
end rate equations

begin 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
end parameters

begin initial conditions
A[0] = Ai
B[0] = Bi
BN[0] = BNi
Bp[0] = Bpi
D[0] = Di
Dusp[0] = Duspi
ERKa[0] = ERKai
ERKt[0] = ERKti
F[0] = Fi
Fgf[0] = Fgfi
K[0] = Ki
Kt[0] = Kti
MAx[0] = MAxi
MDusp[0] = MDuspi
MF[0] = MFi
N[0] = Ni
Na[0] = Nai
Nan[0] = Nani
Rasa[0] = Rasai
Rast[0] = Rasti
Xa[0] = Xai
Xt[0] = Xti
end initial conditions

begin initial values
Ai = 0.1
Bi = 0.1
BNi = 0.001
Bpi = 0.1
Di = 2.0
Duspi = 0.1
ERKai = 0.2
ERKti = 2.0
Fi = 0.001
Fgfi = 1.0
Ki = 3.0
Kti = 3.0
MAxi = 0.1
MDuspi = 0.1
MFi = 0.1
Ni = 0.5
Nai = 0.2
Nani = 0.0
Rasai = 0.5
Rasti = 2.0
Xai = 0.1
Xti = 2.0
end initial values

begin assignment rules
Xi := -Xa[t] + Xt[t]
ERKi := -ERKa[t] + ERKt[t]
Rasi := -Rasa[t] + Rast[t]
AK := -K[t] + Kt[t]
vsFK := (KIG1*vsF)/(KIG1 + K[t])
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units