begin name
middleton1
end name

begin reactions
re1 {1.0}s1 + {1.0}s2 = {1.0}s65
re10 {1.0}s24 + {1.0}s27 = {1.0}s31
re11 {1.0}s31 = {1.0}s24 + {1.0}s27
re12 {1.0}s31 = {1.0}s27 + {1.0}s25
re13 {1.0}s25 + {1.0}s27 = {1.0}s30
re14 {1.0}s30 = {1.0}s25 + {1.0}s27
re15 {1.0}s30 = {1.0}s27 + {1.0}s26
re16 {1.0}s26 + {1.0}s28 = {1.0}s29
re17 {1.0}s29 = {1.0}s26 + {1.0}s28
re18 {1.0}s29 = {1.0}s28 + {1.0}s1
re19 {1.0}s65 = {1.0}s62
re2 {1.0}s65 = {1.0}s2 + {1.0}s1
re20 {1.0}s62 = {1.0}s65
re21 {1.0}$s7 = {1.0}s16
re22 {1.0}$s6 = {1.0}s27
re23 {1.0}$s5 = {1.0}s28
re24 {1.0}$s33 = {1.0}s2
re27 {1.0}$s11 = {1.0}s39
re28 {1.0}$s15 = {1.0}s42
re29 {1.0}$s34 = {1.0}s40
re3 {1.0}s62 + {1.0}s16 = {1.0}s45
re30 {1.0}$s35 = {1.0}s41
re31 {1.0}s40 = {1.0}$s34
re32 {1.0}s42 = {1.0}$s15
re33 {1.0}s39 = {1.0}$s11
re34 {1.0}s41 = {1.0}$s35
re35 {1.0}s27 = {1.0}$s6
re37 {1.0}s2 = {1.0}$s33
re38 {1.0}s28 = {1.0}$s5
re39 {1.0}s62 + {1.0}s16 = {1.0}s36
re4 {1.0}s45 = {1.0}s62 + {1.0}s16
re40 {1.0}s36 = {1.0}s62 + {1.0}s16
re41 {1.0}$s66 = {1.0}s1
re42 {1.0}s23 = {1.0}$s67
re43 {1.0}s24 = {1.0}$s68
re44 {1.0}s25 = {1.0}$s69
re45 {1.0}s26 = {1.0}$s70
re46 {1.0}s1 = {1.0}$s71
re5 {1.0}s45 = {1.0}s62 + {1.0}s22
re6 {1.0}$s3 = {1.0}s23
re7 {1.0}s23 + {1.0}s27 = {1.0}s32
re8 {1.0}s32 = {1.0}s23 + {1.0}s27
re9 {1.0}s32 = {1.0}s27 + {1.0}s24
end reactions

begin rate equations
re1 = re1la*s1[t]*s2[t]
re10 = re10ka15*s24[t]*s27[t]
re11 = re11kd15*s31[t]
re12 = re12km15*s31[t]
re13 = re13ka24*s25[t]*s27[t]
re14 = re14kd24*s30[t]
re15 = re15km24*s30[t]
re16 = re16ka9*s26[t]*s28[t]
re17 = re17kd9*s29[t]
re18 = re18km9*s29[t]
re19 = re19q*s65[t]
re2 = re2ld*s65[t]
re20 = re20p*s62[t]
re21 = re21deltaDELLA*s40[t]
re22 = re22deltaGA20ox*s39[t]
re23 = re23deltaGA3ox*s41[t]
re24 = re24deltaGID*s42[t]
re27 = (muGA20ox*s16[t])/(re27thetaGA20ox + s16[t])
re28 = (muGID*s16[t])/(re28thetaGID + s16[t])
re29 = (muDELLA*re29thetaDELLA)/(re29thetaDELLA + s16[t])
re3 = re3ua1*s16[t]*s62[t]
re30 = (muGA3ox*s16[t])/(re30thetaGA3ox + s16[t])
re31 = muDELLA*s40[t]
re32 = muGID*s42[t]
re33 = muGA20ox*s39[t]
re34 = muGA3ox*s41[t]
re35 = re35gammaGA20ox*s27[t]
re37 = re37gammaGID*s2[t]
re38 = re38gammaGA3ox*s28[t]
re39 = re39ua2*s16[t]*s62[t]
re4 = re4ud1*s45[t]
re40 = re40ud2*s36[t]
re41 = omegaGA4*Pmem*re41A1
re42 = muGA*s23[t]
re43 = muGA*s24[t]
re44 = muGA*s25[t]
re45 = muGA*s26[t]
re46 = (muGA + Pmem*re46B1)*s1[t]
re5 = re5um*s45[t]
re6 = omegaGA12
re7 = re7ka12*s23[t]*s27[t]
re8 = re8kd12*s32[t]
re9 = re9km12*s32[t]
end rate equations

begin parameters
Pmem = 2.66664
R = 1.0
appliedGA4 = 2.0
muDELLA = 0.070794578438414
muGA = 0.290804218727464
muGA20ox = 0.047770755070625
muGA3ox = 0.102600014140148
muGID = 0.045708818961487
omegaGA12ga13 = 0.006602803853512
omegaGA4 = 0.0
tGA4off = 620.0
tGA4on = 500.0
s11 = 0.0
s15 = 0.0
s3 = 0.0
s33 = 0.0
s34 = 0.0
s35 = 0.0
s5 = 0.0
s6 = 0.0
s66 = 0.0
s67 = 0.0
s68 = 0.0
s69 = 0.0
s7 = 0.0
s70 = 0.0
s71 = 0.0
re1la = 1.35
re2ld = 2.84315148627376
re3ua1 = 10.0
re4ud1 = 0.133045441797809
re5um = 6.92208879449283
re7ka12 = 2904.11853677638
re8kd12 = 2.67298621993027
re9km12 = 198.80427707769
re10ka15 = 2073.22402517968
re11kd15 = 0.008827838388125
re12km15 = 763.777072066507
re13ka24 = 3099.18915892587
re14kd24 = 0.01588492846351
re15km24 = 2.58846077319221
re16ka9 = 2073.22402517968
re17kd9 = 0.008827838388125
re18km9 = 763.777072066507
re19q = 0.025118864315096
re20p = 0.077624711662869
re21deltaDELLA = 0.000527749140286577
re22deltaGA20ox = 0.192990314378105
re23deltaGA3ox = 0.019299031437811
re24deltaGID = 19.2990314378105
re27thetaGA20ox = 0.6383
re28thetaGID = 0.00055995
re29thetaDELLA = 0.01
re30thetaGA3ox = 0.0082
re35gammaGA20ox = 3.514
re37gammaGID = 3.514
re38gammaGA3ox = 3.514
re39ua2 = 316.2278
re40ud2 = 2.8184
re41A1 = 0.0307
re46B1 = 0.00039795
default = 1.0
end parameters

begin initial conditions
s1[0] = s1i
s16[0] = s16i
s2[0] = s2i
s22[0] = s22i
s23[0] = s23i
s24[0] = s24i
s25[0] = s25i
s26[0] = s26i
s27[0] = s27i
s28[0] = s28i
s29[0] = s29i
s30[0] = s30i
s31[0] = s31i
s32[0] = s32i
s36[0] = s36i
s39[0] = s39i
s40[0] = s40i
s41[0] = s41i
s42[0] = s42i
s45[0] = s45i
s62[0] = s62i
s65[0] = s65i
end initial conditions

begin initial values
s1i = 0.01
s16i = 0.01
s2i = 0.01
s22i = 0.01
s23i = 0.01
s24i = 0.01
s25i = 0.01
s26i = 0.01
s27i = 0.01
s28i = 0.01
s29i = 0.01
s30i = 0.01
s31i = 0.01
s32i = 0.01
s36i = 0.01
s39i = 0.01
s40i = 0.01
s41i = 0.01
s42i = 0.01
s45i = 0.01
s62i = 0.01
s65i = 0.01
end initial values

begin assignment rules
omegaGA12 := omegaGA12ga13*R
end assignment rules

begin function definitions
end function definitions

begin events
AddGA4 = (t > tGA4on, {omegaGA4 -> appliedGA4}, 0)
RemoveGA4 = (t > tGA4off, {omegaGA4 -> 0}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units