begin name
krab
end name

begin reactions
v1 {1.0}R00 = {1.0}A00
v10 {1.0}F01 = {1.0}F02
v11 {1.0}F02 = {1.0}F03
v12 {1.0}F03 = {1.0}F04
v13 {1.0}F04 = {1.0}F05
v14 {1.0}F05 = {1.0}O00
v15 {1.0}O00 = {1.0}O01
v16 {1.0}O01 = {1.0}O02
v17 {1.0}O02 = {1.0}O03
v18 {1.0}O03 = {1.0}O04
v19 {1.0}O04 = {1.0}O05
v2 {1.0}A00 = {1.0}P00
v20 {1.0}O05 = {1.0}E00
v21 {1.0}E00 = {1.0}E01
v22 {1.0}E01 = {1.0}E02
v23 {1.0}E02 = {1.0}E03
v24 {1.0}E03 = {1.0}E04
v25 {1.0}E04 = {1.0}E05
v26 {1.0}E05 = {1.0}R00
v3 {1.0}P00 = {1.0}P01
v4 {1.0}P01 = {1.0}P02
v5 {1.0}P02 = {1.0}P03
v6 {1.0}P03 = {1.0}P04
v7 {1.0}P04 = {1.0}P05
v8 {1.0}P05 = {1.0}F00
v9 {1.0}F00 = {1.0}F01
end reactions

begin rate equations
v1 = -(kb01*A00[t]) + kf01*O2*R00[t]
v10 = (kf10*F01[t])/2.718281828^(19.46188812245315*deltapsi*q2) - 2.718281828^(19.46188812245315*deltapsi*q2)*kb10*F02[t]
v11 = (kf11*F02[t])/(2.718281828^(19.46188812245315*deltapsi*q3)*10^pHM) - 2.718281828^(19.46188812245315*deltapsi*q3)*kb11*F03[t]
v12 = kf12*F03[t] - kb12*F04[t]
v13 = (kf13*F04[t])/(2.718281828^(19.46188812245315*deltapsi*(1 - q2))*10^pHM) - 2.718281828^(19.46188812245315*deltapsi*(1 - q2))*kb13*F05[t]
v14 = (kf14*F05[t])/2.718281828^(19.46188812245315*deltapsi*(1 - q3)) - (2.718281828^(19.46188812245315*deltapsi*(1 - q3))*kb14*O00[t])/10^pHC
v15 = cred*kf15*O00[t] - (1 - cred)*kb15*O01[t]
v16 = (kf16*O01[t])/2.718281828^(19.46188812245315*deltapsi*q2) - 2.718281828^(19.46188812245315*deltapsi*q2)*kb16*O02[t]
v17 = (kf17*O02[t])/(2.718281828^(19.46188812245315*deltapsi*q3)*10^pHM) - 2.718281828^(19.46188812245315*deltapsi*q3)*kb17*O03[t]
v18 = kf18*O03[t] - kb18*O04[t]
v19 = (kf19*O04[t])/(2.718281828^(19.46188812245315*deltapsi*(1 - q2))*10^pHM) - 2.718281828^(19.46188812245315*deltapsi*(1 - q2))*kb19*O05[t]
v2 = kf02*A00[t] - kb02*P00[t]
v20 = -((2.718281828^(19.46188812245315*deltapsi*(1 - q3))*kb20*E00[t])/10^pHC) + (kf20*O05[t])/2.718281828^(19.46188812245315*deltapsi*(1 - q3))
v21 = cred*kf21*E00[t] - (1 - cred)*kb21*E01[t]
v22 = (kf22*E01[t])/2.718281828^(19.46188812245315*deltapsi*q2) - 2.718281828^(19.46188812245315*deltapsi*q2)*kb22*E02[t]
v23 = (kf23*E02[t])/(2.718281828^(19.46188812245315*deltapsi*q3)*10^pHM) - 2.718281828^(19.46188812245315*deltapsi*q3)*kb23*E03[t]
v24 = kf24*E03[t] - kb24*E04[t]
v25 = (kf25*E04[t])/(2.718281828^(19.46188812245315*deltapsi*(1 - q2))*10^pHM) - 2.718281828^(19.46188812245315*deltapsi*(1 - q2))*kb25*E05[t]
v26 = (kf26*E05[t])/2.718281828^(19.46188812245315*deltapsi*(1 - q3)) - (2.718281828^(19.46188812245315*deltapsi*(1 - q3))*kb26*R00[t])/10^pHC
v3 = cred*kf03*P00[t] - (1 - cred)*kb03*P01[t]
v4 = (kf04*P01[t])/2.718281828^(19.46188812245315*deltapsi*q2) - 2.718281828^(19.46188812245315*deltapsi*q2)*kb04*P02[t]
v5 = (kf05*P02[t])/(2.718281828^(19.46188812245315*deltapsi*q3)*10^pHM) - 2.718281828^(19.46188812245315*deltapsi*q3)*kb05*P03[t]
v6 = kf06*P03[t] - kb06*P04[t]
v7 = (kf07*P04[t])/(2.718281828^(19.46188812245315*deltapsi*(1 - q2))*10^pHM) - 2.718281828^(19.46188812245315*deltapsi*(1 - q2))*kb07*P05[t]
v8 = -((2.718281828^(19.46188812245315*deltapsi*(1 - q3))*kb08*F00[t])/10^pHC) + (kf08*P05[t])/2.718281828^(19.46188812245315*deltapsi*(1 - q3))
v9 = cred*kf09*F00[t] - (1 - cred)*kb09*F01[t]
end rate equations

begin parameters
EmaccE = 0.65
EmaccF = 0.35
EmaccP = 0.25
Emcyta = 0.27
O2 = 0.0001
cred = 0.04
deltapsi = 0.18
kb01 = 19200.0
kb02 = 64.0363912
kb15 = 111548.6393
kb16 = 31464.85383
kb17 = 55.55555556
kb18 = 400000000.0
kb19 = 0.078120117
kb20 = 72727272.73
kf01 = 192000000.0
kf02 = 28000.0
kf15 = 200000.0
kf16 = 68535.14617
kf17 = 66666666667.0
kf18 = 600000000.0
kf19 = 12499218799.0
kf20 = 400.0
pHC = 7.0
pHM = 7.0
pKaE = 6.0
pKaF = 7.0
phiE03 = 0.0
phiE04 = 0.0
phiF03 = 0.0
phiF04 = 0.0
phiP03 = 0.0
phiP04 = 0.0
q2 = 0.33
q3 = 0.67
rtotal = 0.002581067
defaultcompartment = 1.0
end parameters

begin initial conditions
A00[0] = A00i
E00[0] = E00i
E01[0] = E01i
E02[0] = E02i
E03[0] = E03i
E04[0] = E04i
E05[0] = E05i
F00[0] = F00i
F01[0] = F01i
F02[0] = F02i
F03[0] = F03i
F04[0] = F04i
F05[0] = F05i
O00[0] = O00i
O01[0] = O01i
O02[0] = O02i
O03[0] = O03i
O04[0] = O04i
O05[0] = O05i
P00[0] = P00i
P01[0] = P01i
P02[0] = P02i
P03[0] = P03i
P04[0] = P04i
P05[0] = P05i
R00[0] = R00i
end initial conditions

begin initial values
A00i = 0.0
E00i = 0.0
E01i = 0.0
E02i = 0.0
E03i = 0.0
E04i = 0.0
E05i = 0.0
F00i = 0.0
F01i = 0.0
F02i = 0.0
F03i = 0.0
F04i = 0.0
F05i = 0.0
O00i = 0.0
O01i = 0.0
O02i = 0.0
O03i = 0.0
O04i = 0.0
O05i = 0.0
P00i = 0.0
P01i = 0.0
P02i = 0.0
P03i = 0.0
P04i = 0.0
P05i = 0.0
R00i = 1.0
end initial values

begin assignment rules
kb24 := kb18*rE03^((-1 + phiE03)/2)
kf12 := kf18*rF03^((1 + phiF03)/2)
rF03 := (2.718281828^(38.9237762449063*(EmaccF - Emcyta))*kb18)/kf18
kb04 := kb16
kf06 := kf18*rP03^((1 + phiP03)/2)
kf13 := kf19*rF04^((1 + phiF04)/2)
kf24 := kf18*rE03^((1 + phiE03)/2)
kf04 := kf16
kf26 := kf20
kf22 := kf16
rP04 := rtotal/(rE03*rE04*rF03*rF04*rP03)
rP03 := (2.718281828^(38.9237762449063*(EmaccP - Emcyta))*kb18)/kf18
kb03 := kb15
rE03 := (2.718281828^(38.9237762449063*(EmaccE - Emcyta))*kb18)/kf18
kf14 := kf20
kb12 := kb18*rF03^((-1 + phiF03)/2)
kf07 := kf19*rP04^((1 + phiP04)/2)
kf10 := kf16
kb13 := kb19*rF04^((-1 + phiF04)/2)
kf08 := kf20
kf21 := kf15
kb11 := kb17
kb09 := kb15
kb10 := kb16
kb05 := kb17
rF04 := (10^pKaF*kb19)/kf19
kf11 := kf17
kb22 := kb16
kb21 := kb15
kf09 := kf15
kb26 := kb20
kb14 := kb20
kb08 := kb20
kb07 := kb19*rP04^((-1 + phiP04)/2)
kb23 := kb17
kf03 := kf15
kf25 := kf19*rE04^((1 + phiE04)/2)
kf05 := kf17
kf23 := kf17
rE04 := (10^pKaE*kb19)/kf19
kb06 := kb18*rP03^((-1 + phiP03)/2)
kb25 := kb19*rE04^((-1 + phiE04)/2)
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