begin name
liu1
end name

begin reactions
reaction1 {1.0}$EXTERNAL = {1.0}Auxin
reaction10 {1.0}$EXTERNAL = {1.0}Restar
reaction11 {1.0}Restar = {1.0}$EXTERNAL
reaction12 {1.0}$EXTERNAL = {1.0}ET
reaction13 {1.0}ET = {1.0}$EXTERNAL
reaction14 {1.0}$EXTERNAL = {1.0}CTR1star
reaction15 {1.0}CTR1star = {1.0}$EXTERNAL
reaction16 {1.0}$EXTERNAL = {1.0}X
reaction17 {1.0}X = {1.0}$EXTERNAL
reaction18 {1.0}$EXTERNAL = {1.0}CK
reaction19 {1.0}CK = {1.0}$EXTERNAL
reaction2 {1.0}$EXTERNAL = {1.0}Auxin
reaction3 {1.0}Auxin = {1.0}$EXTERNAL
reaction4 {1.0}$EXTERNAL = {1.0}Rastar
reaction5 {1.0}Rastar = {1.0}$EXTERNAL
reaction6 {1.0}$EXTERNAL = {1.0}PLSm
reaction7 {1.0}PLSm = {1.0}$EXTERNAL
reaction8 {1.0}$EXTERNAL = {1.0}PLSp
reaction9 {1.0}PLSp = {1.0}$EXTERNAL
vAuxin {1.0}$IAA = {1.0}Auxin
vCytokinin {1.0}$CKex = {1.0}CK
vEthylene {1.0}$ACC = {1.0}ET
end reactions

begin rate equations
reaction1 = compartment1*function2[k1a, X[t], k1]
reaction10 = compartment1*function10[k10, PLSp[t], k10a, Re]
reaction11 = compartment1*function12[Restar[t], ET[t], k11]
reaction12 = compartment1*function13[Auxin[t], CK[t], k12, k12a]
reaction13 = compartment1*function15[ET[t], k13]
reaction14 = compartment1*function16[Restar[t], k14, CTR1]
reaction15 = compartment1*function17[CTR1star[t], k15]
reaction16 = compartment1*function18[CTR1star[t], k16, k16a]
reaction17 = compartment1*function19[X[t], k17]
reaction18 = compartment1*function20[Auxin[t], k18a, k18]
reaction19 = compartment1*function21[CK[t], k19]
reaction2 = compartment1*function3[k2, k2a, ET[t], CK[t], k2b, PLSp[t], k2c]
reaction3 = compartment1*function1[k3, k3a, X[t], Auxin[t]]
reaction4 = compartment1*function4[k4, Auxin[t], Ra]
reaction5 = compartment1*k5*Rastar[t]
reaction6 = compartment1*function6[k6, Rastar[t], ET[t], k6a]
reaction7 = compartment1*function7[k7, PLSm[t]]
reaction8 = compartment1*function8[k8, PLSm[t]]
reaction9 = compartment1*function9[k9, PLSp[t]]
vAuxin = compartment1*IAA*kauxin
vCytokinin = CKex*compartment1*kcytokinin
vEthylene = ACC*compartment1*kethylene
end rate equations

begin parameters
k1 = 1.0
k10 = 0.0003
k10a = 0.5
k11 = 5.0
k12 = 0.1
k12a = 0.1
k13 = 1.0
k14 = 3.0
k15 = 0.085
k16 = 0.3
k16a = 1.0
k17 = 0.1
k18 = 0.1
k18a = 1.0
k19 = 1.0
k1a = 1.0
k2 = 0.2
k2a = 2.8
k2b = 1.0
k2c = 0.01
k3 = 2.0
k3a = 0.45
k4 = 1.0
k5 = 1.0
k6 = 0.3
k6a = 0.2
k7 = 1.0
k8 = 1.0
k9 = 1.0
kauxin = 70.0
kcytokinin = 10.0
kethylene = 0.5
ACC = 0.0
CKex = 0.0
CTR1T = 0.3
IAA = 0.0
RaT = 1.0
ReT = 0.3
compartment1 = 1.0
extra = 1.0
end parameters

begin initial conditions
Auxin[0] = Auxini
CK[0] = CKi
CTR1star[0] = CTR1stari
ET[0] = ETi
PLSm[0] = PLSmi
PLSp[0] = PLSpi
Rastar[0] = Rastari
Restar[0] = Restari
X[0] = Xi
end initial conditions

begin initial values
Auxini = 0.1
CKi = 0.1
CTR1stari = 0.3
ETi = 0.1
PLSmi = 0.1
PLSpi = 0.1
Rastari = 1.0
Restari = 0.3
Xi = 0.1
end initial values

begin assignment rules
CTR1 := CTR1T - CTR1star[t]
Re := ReT - Restar[t]
Ra := RaT - Rastar[t]
end assignment rules

begin function definitions
function16[Restar_,k14_,CTR1_] := CTR1*k14*Restar
function15[ET_,k13_] := ET*k13
function3[k2_,k2a_,ET_,CK_,k2b_,PLSp_,k2c_] := k2 + (ET*k2a*PLSp)/((1 + CK/k2b)*(k2c + PLSp))
function2[k1a_,X_,k1_] := k1a/(1 + X/k1)
function1[k3_,k3a_,X_,Auxin_] := Auxin*(k3 + k3a*X)
function21[CK_,k19_] := CK*k19
function19[X_,k17_] := k17*X
function18[CTR1star_,k16_,k16a_] := k16 - CTR1star*k16a
function20[Auxin_,k18a_,k18_] := k18a/(1 + Auxin/k18)
function17[CTR1star_,k15_] := CTR1star*k15
function13[Auxin_,CK_,k12_,k12a_] := k12 + Auxin*CK*k12a
function12[Restar_,ET_,k11_] := ET*k11*Restar
function10[k10_,PLSp_,k10a_,Re_] := (k10 + k10a*PLSp)*Re
function9[k9_,PLSp_] := k9*PLSp
function8[k8_,PLSm_] := k8*PLSm
function7[k7_,PLSm_] := k7*PLSm
function6[k6_,Rastar_,ET_,k6a_] := (k6*Rastar)/(1 + ET/k6a)
function4[k4_,Auxin_,Ra_] := Auxin*k4*Ra
end function definitions

begin events
end events

begin process annotations
reaction_8 http://identifiers.org/obo.go/GO:0006412
end process annotations

begin species annotations
ACC http://identifiers.org/obo.chebi/CHEBI:18053
ET http://identifiers.org/obo.chebi/CHEBI:18153
ET http://identifiers.org/kegg.compound/C06547
IAA http://identifiers.org/obo.chebi/CHEBI:16411
IAA http://identifiers.org/kegg.compound/C00954
end species annotations

begin units
end units