begin name
hong1
end name

begin reactions
r1 {1.0}$EXTERNAL = {1.0}M
r10 {1.0}$EXTERNAL = {1.0}TF + {1.0}CP2
r11 {1.0}CP2 = {1.0}$EXTERNAL
r12 {1.0}TF + {1.0}CP2 = {1.0}$EXTERNAL
r13 {1.0}CP2 = {1.0}$EXTERNAL
r14 {1.0}CP2 = {1.0}$EXTERNAL
r15 {1.0}$EXTERNAL = {1.0}TF
r16 {1.0}$EXTERNAL = {1.0}TF
r17 {1.0}$EXTERNAL = {1.0}TF
r2 {1.0}$EXTERNAL = {1.0}M
r3 {1.0}M = {1.0}$EXTERNAL
r4 {1.0}$EXTERNAL = {1.0}CP
r5 {1.0}CP = {1.0}$EXTERNAL
r6 {2.0}CP = {1.0}CP2
r7 {1.0}CP2 = {2.0}CP
r8 {1.0}CP = {1.0}$EXTERNAL
r9 {1.0}CP = {1.0}$EXTERNAL
end reactions

begin rate equations
r1 = function4r1*system
r10 = IC*kicd*system
r11 = kcp2d*system*CP2[t]
r12 = kica*system*CP2[t]*TF[t]
r13 = function4r13*system
r14 = chk2*system*CP2[t]
r15 = IC*kcp2d*system
r16 = chk2c*IC*system
r17 = function4r17*system
r2 = function4r2*system
r3 = kmd*system*M[t]
r4 = function1*system
r5 = kcpd*system*CP[t]
r6 = ka*system*CP[t]^2
r7 = kd*system*CP2[t]
r8 = function4r8*system
r9 = chk2*system*CP[t]
end rate equations

begin parameters
Dex = 0.0
ICtot = 1.0
IRpulseduration = 2.0
J = 0.3
Jp = 0.05
TFtot = 0.5
chk2ToggleOnOff = 1.0
chk2value = 0.25
chk2cvalue = 0.1
chk2pulseactive = 1.0
chk2pulseactiveelim = 1.0
chk2pulseactiveelimelim = 0.0
initialValueCP = 0.039
initialValueCP2 = 0.088
initialValueM = 1.35
initialValueTF = 0.06
ka = 100.0
kcp2d = 0.0525
kcpd = 0.525
kcps = 0.5
kd = 0.01
kica = 20.0
kicd = 0.01
kmd = 0.1
kms = 1.0
kp1 = 10.0
kp2 = 0.1
n = 2.0
trigONtime = 4.0
system = 1.0
end parameters

begin initial conditions
CP[0] = CPi
CP2[0] = CP2i
M[0] = Mi
TF[0] = TFi
end initial conditions

begin initial values
CPi = 0.039
CP2i = 0.088
Mi = 1.35
TFi = 0.06
end initial values

begin assignment rules
IC := TFtot - TF[t]
function1 := kcps*M[t]
function4r17 := (IC*kp2)/(system*(2*IC + Jp + CP[t] + 2*CP2[t]))
function4r2 := (kms*TF[t]^n)/(system*(J^n + TF[t]^n))
chk2pulseonoff := chk2pulseactive - chk2pulseactiveelim - chk2pulseactiveelimelim
function4r8 := (kp1*CP[t])/(system*(2*IC + Jp + CP[t] + 2*CP2[t]))
function4r13 := (kp2*CP2[t])/(system*(2*IC + Jp + CP[t] + 2*CP2[t]))
function4r1 := Dex/system
CPtot := 2*IC + CP[t] + 2*CP2[t]
trigOFFtime := IRpulseduration + trigONtime
chk2c := (chk2cvalue*chk2pulseonoff*chk2ToggleOnOff)/IRpulseduration
chk2 := (chk2pulseonoff*chk2ToggleOnOff*chk2value)/IRpulseduration
end assignment rules

begin function definitions
end function definitions

begin events
trigON = (t >= trigONtime, {chk2pulseactiveelim -> 0}, 0)
trigOFF = (t >= trigOFFtime, {chk2pulseactiveelimelim -> 1}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units