begin name
muraro1
end name

begin reactions
r1 {1.0}$EXTERNAL = {1.0}IAAm
r10 {1.0}$EXTERNAL = {1.0}PINp
r11 {1.0}$EXTERNAL = {1.0}ARm
r12 {1.0}$EXTERNAL = {1.0}ARp
r13 {1.0}$EXTERNAL = {1.0}CRm
r14 {1.0}$EXTERNAL = {1.0}CRp
r15 {1.0}$EXTERNAL = {1.0}AHKph
r16 {1.0}$EXTERNAL = {1.0}Ck
r17 {1.0}$EXTERNAL = {1.0}ARRBph
r18 {1.0}$EXTERNAL = {1.0}ARRAph
r19 {1.0}$EXTERNAL = {1.0}ARRAm
r2 {1.0}$EXTERNAL = {1.0}IAAp
r20 {1.0}$EXTERNAL = {1.0}ARRAp
r3 {1.0}$EXTERNAL = {1.0}AuxTIR1
r4 {1.0}$EXTERNAL = {1.0}AuxTIAA
r5 {1.0}$EXTERNAL = {1.0}IAAs
r6 {1.0}$EXTERNAL = {1.0}ARFIAA
r7 {1.0}$EXTERNAL = {1.0}ARF2
r8 {1.0}$EXTERNAL = {1.0}Aux
r9 {1.0}$EXTERNAL = {1.0}PINm
end reactions

begin rate equations
r1 = (F2 + F1*lambda1 + F3*lambda3)*phiIAAp - IAAm[t]
r10 = (deltaPINp*PINm[t] - PINp[t])/eps
r11 = (F5b + F5a*lambda1)*phiARp - ARm[t]
r12 = (deltaARp*ARm[t] - ARp[t])/eps
r13 = F4*phiCRp - CRm[t]
r14 = (deltaCRp*CRm[t] - CRp[t])/eps
r15 = (CkAHKph*rd - ra*AHKph[t]*Ck[t])/eps
r16 = muCk*(alphaCk - Ck[t]) - (etaCkPh*(-(CkAHKph*rd) + ra*AHKph[t]*Ck[t]))/eps
r17 = (ARRBp*CkAHKph*ua - CkAHK*ud*ARRBph[t])/eps
r18 = (CkAHKph*sa*ARRAp[t] - CkAHK*sd*ARRAph[t])/eps
r19 = F6*phiARRAp - ARRAm[t]
r2 = etaARFIAA*(pd*ARFIAA[t] - ARF*pa*IAAp[t]) + (ld*AuxTIAA[t] + deltaIAAp*IAAm[t] - la*AuxTIR1[t]*IAAp[t])/eps
r20 = (deltaARRAp*ARRAm[t] - ARRAp[t] + etaAHKph*(-(CkAHKph*sa*ARRAp[t]) + CkAHK*sd*ARRAph[t]))/eps
r3 = (ka*TIR1*Aux[t] + (1 + ld)*AuxTIAA[t] - kd*AuxTIR1[t] - la*AuxTIR1[t]*IAAp[t])/eps
r4 = (-((1 + ld)*AuxTIAA[t]) + la*AuxTIR1[t]*IAAp[t])/eps
r5 = (AuxTIAA[t] - muIAAs*IAAs[t])/eps
r6 = -(pd*ARFIAA[t]) + ARF*pa*IAAp[t]
r7 = ARF^2*qa - qd*ARF2[t]
r8 = muAux*(alphaAux - Aux[t]) - (etaAuxTIR1*(ka*TIR1*Aux[t] - kd*AuxTIR1[t]))/eps
r9 = (F5b + F5a*lambda1)*phiPINp - PINm[t]
end rate equations

begin parameters
alphaAHK = 1.0
alphaARF = 1.0
alphaARRB = 2.0
alphaAux = 1.0
alphaCk = 1.0
alphaPH = 1.0
alphaTIR1 = 1.0
deltaARRAp = 1.0
deltaARp = 1.0
deltaCRp = 1.0
deltaIAAp = 1.0
deltaPINp = 1.0
eps = 0.01
etaAHKph = 1.0
etaARFIAA = 1.0
etaAuxTIR1 = 10.0
etaCkPh = 1.0
ka = 100.0
kd = 1.0
la = 0.5
lambda1 = 0.1
lambda3 = 0.02
ld = 0.1
muAux = 0.1
muCk = 0.1
muIAAs = 1.0
pa = 10.0
pd = 10.0
phiARRAp = 100.0
phiARp = 2.0
phiCRp = 2.0
phiIAAp = 100.0
phiPINp = 100.0
psiARF = 0.1
psiARFIAA = 0.1
qa = 1.0
qd = 1.0
ra = 1.0
rd = 1.0
sa = 1.0
sd = 1.0
thARFIAA = 0.1
thARRAph = 0.1
thARRBph = 0.1
thetaARF = 0.1
thetaARF2 = 0.01
thetaARp = 0.1
thetaIAAp = 0.1
ua = 1.0
ud = 1.0
cell = 1.0
end parameters

begin initial conditions
AHKph[0] = AHKphi
ARF2[0] = ARF2i
ARFIAA[0] = ARFIAAi
ARRAm[0] = ARRAmi
ARRAp[0] = ARRApi
ARRAph[0] = ARRAphi
ARRBph[0] = ARRBphi
ARm[0] = ARmi
ARp[0] = ARpi
Aux[0] = Auxi
AuxTIAA[0] = AuxTIAAi
AuxTIR1[0] = AuxTIR1i
CRm[0] = CRmi
CRp[0] = CRpi
Ck[0] = Cki
IAAm[0] = IAAmi
IAAp[0] = IAApi
IAAs[0] = IAAsi
PINm[0] = PINmi
PINp[0] = PINpi
end initial conditions

begin initial values
AHKphi = 1.0
ARF2i = 0.0
ARFIAAi = 0.0
ARRAmi = 0.0
ARRApi = 0.0
ARRAphi = 0.0
ARRBphi = 0.0
ARmi = 0.0
ARpi = 0.0
Auxi = 1.0
AuxTIAAi = 0.0
AuxTIR1i = 0.0
CRmi = 0.0
CRpi = 0.0
Cki = 1.0
IAAmi = 0.0
IAApi = 0.0
IAAsi = 0.0
PINmi = 0.0
PINpi = 0.0
end initial values

begin assignment rules
CkAHK := alphaAHK - etaAHKph*(CkAHKph + AHKph[t])
CkAHKph := alphaPH - AHKph[t] - ARRAph[t] - ARRBph[t]
ARRBp := alphaARRB - etaAHKph*ARRBph[t]
ARF := alphaARF - 2*ARF2[t] - ARFIAA[t]
TIR1 := alphaTIR1 - AuxTIAA[t] - AuxTIR1[t]
F6 := ARp[t]/(thetaARp*(1 + ARp[t]/thetaARp))
F5b := (ARF^2/psiARF + ARF2[t]/thetaARF2)/(1 + ARF^2/psiARF + ARF/thetaARF + ARF2[t]/thetaARF2 + ARFIAA[t]/thARFIAA + (ARF*IAAp[t])/psiARFIAA)
F5a := ARF/(thetaARF*(1 + ARF^2/psiARF + ARF/thetaARF + ARF2[t]/thetaARF2 + ARFIAA[t]/thARFIAA + (ARF*IAAp[t])/psiARFIAA))
F4 := ARRBph[t]/(thARRBph*(1 + ARRAph[t]/thARRAph + ARRBph[t]/thARRBph))
F3 := ARRBph[t]/(thARRBph*(1 + ARF^2/psiARF + ARF/thetaARF + ARF2[t]/thetaARF2 + ARFIAA[t]/thARFIAA + ARRBph[t]/thARRBph + (ARF*IAAp[t])/psiARFIAA))
F2 := (ARF^2/psiARF + ARF2[t]/thetaARF2)/(1 + ARF^2/psiARF + ARF/thetaARF + ARF2[t]/thetaARF2 + ARFIAA[t]/thARFIAA + ARRBph[t]/thARRBph + (ARF*IAAp[t])/psiARFIAA)
F1 := ARF/(thetaARF*(1 + ARF^2/psiARF + ARF/thetaARF + ARF2[t]/thetaARF2 + ARFIAA[t]/thARFIAA + ARRBph[t]/thARRBph + (ARF*IAAp[t])/psiARFIAA))
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