begin name
komarova1
end name

begin reactions
R1 {1.0}$EXTERNAL = {1.0}x1
R2 {1.0}x1 = {1.0}$EXTERNAL
R3 {1.0}$EXTERNAL = {1.0}x2
R4 {1.0}x2 = {1.0}$EXTERNAL
R5 {1.0}z = {1.0}$EXTERNAL
R6 {1.0}$EXTERNAL = {1.0}z
end reactions

begin rate equations
R1 = Ratex1appearance + alpha1*(1 - Fig7AR1alphaOnOff + Fig7AR1alphaOnOff*x1[t]^g11*x2[t]^g21)
R2 = beta1*x1[t]
R3 = alpha2*x1[t]^g12*x2[t]^g22
R4 = beta2*x2[t]
R5 = k1*y1
R6 = k2*y2
end rate equations

begin parameters
Fig7AR1alphaOnOff = 1.0
Ratex1appearanceOnOff = 0.0
Ratex1appearanceOnOffPiecewise = 0.0
Ratex1appearanceduration = 1.0
Ratex1appearancemagnitude = 7.0
Ratex1appearancetimefirst = 10.0
Ratex1appearancetimesecond = 150.0
alpha1 = 7.0
alpha2 = 7.0
beta1basalvalue = 0.2
beta1increaseOnOff = 1.0
beta1increaseOntime = 500.0
beta1increasevalue = 0.0
beta1increasevalue2 = 0.03
beta2 = 0.02
flagformation = 1.0
flagresorption = 1.0
g11 = 1.09
g12 = 1.0
g21 = -0.5
g22 = 0.09
initialvaluex1 = 5.08098
initialvaluex2 = 1397.73
initialvaluez = 90.0
k1 = 0.26
k2 = 0.0008
plotzyaxisdenominator = 113.749
default = 1.0
end parameters

begin initial conditions
x1[0] = x1i
x2[0] = x2i
z[0] = zi
end initial conditions

begin initial values
x1i = 5.08098
x2i = 1397.73
zi = 90.0
end initial values

begin assignment rules
y1 := flagresorption*(-x1bar + x1[t])
x2bar := (beta1/alpha1)^(g12/gamma)*(beta2/alpha2)^((1 - g11)/gamma)
x1variable := x1[t]
Ratex1appearance := Ratex1appearancemagnitude*Ratex1appearanceOnOff*Ratex1appearanceOnOffPiecewise
y2 := flagformation*(-x2bar + x2[t])
x2variable := x2[t]
gamma := g12*g21 - (1 - g11)*(1 - g22)
beta1 := beta1basalvalue + beta1increaseOnOff*beta1increasevalue
plotzyaxis := (100*z[t])/plotzyaxisdenominator
x1bar := (beta1/alpha1)^((1 - g22)/gamma)*(beta2/alpha2)^(g21/gamma)
plotx2yaxis := -initialvaluex2 + y2 + x2[t]
end assignment rules

begin function definitions
end function definitions

begin events
y2_gt_trigger = (x2[t] > x2bar, {flagformation -> 1}, 0)
Rx1_geq_2nd_trig = (t >= Ratex1appearancetimesecond, {Ratex1appearanceOnOffPiecewise -> 1}, 0)
y1_gt_trigger = (x1[t] > x1bar, {flagresorption -> 1}, 0)
y1_leq_trigger = (x1[t] <= x1bar, {flagresorption -> 0}, 0)
Rx1_gt_1stdur_trig = (t > Ratex1appearanceduration + Ratex1appearancetimefirst, {Ratex1appearanceOnOffPiecewise -> 0}, 0)
Rx1_lt_1st_trigger = (t < Ratex1appearancetimefirst, {Ratex1appearanceOnOffPiecewise -> 0}, 0)
beta1_geq_on_trigger = (t >= beta1increaseOntime, {beta1increasevalue -> beta1increasevalue2, flagresorption -> 1}, 0)
Rx1_geq_1st_trigger = (t >= Ratex1appearancetimefirst, {Ratex1appearanceOnOffPiecewise -> 1}, 0)
Rx1_gt_2nddur_trig = (t > Ratex1appearanceduration + Ratex1appearancetimesecond, {Ratex1appearanceOnOffPiecewise -> 0}, 0)
y2_leq_trigger = (x2[t] <= x2bar, {flagformation -> 0}, 0)
end events

begin process annotations
R5 http://identifiers.org/go/GO:0045453
end process annotations

begin species annotations
end species annotations

begin units
end units