begin name
kongas1
end name

begin reactions
ADPdiffusion {1.0}ADPi = {1.0}ADP
ATPdiffusion {1.0}ATPi = {1.0}ATP
ATPase {1.0}ATP = {1.0}ADP + {1.0}P
Crdiffusion {1.0}Cri = {1.0}Cr
MMCK {1.0}ATP + {1.0}Cr = {1.0}PCr + {1.0}ADP
MiCK {1.0}ATPi + {1.0}Cri = {1.0}ADPi + {1.0}PCri
OxPhos {1.0}ADPi + {1.0}Pi = {1.0}ATPi
PCrdiffusion {1.0}PCri = {1.0}PCr
Pidiffusion {1.0}Pi = {1.0}P
end reactions

begin rate equations
ADPdiffusion = -(ADPdiffusionk17*CYT*ADP[t]) + ADPdiffusionk17*IMS*ADPi[t]
ATPdiffusion = -(ATPdiffusionk19*CYT*ATP[t]) + ATPdiffusionk19*IMS*ATPi[t]
ATPase = ATPasev4*CYT*ATP[t]
Crdiffusion = -(Crdiffusionk16*CYT*Cr[t]) + Crdiffusionk16*IMS*Cri[t]
MMCK = (CYT*((MMCKVf3*ATP[t]*Cr[t])/(MMCKKb3*MMCKKia3) - (MMCKVb3*ADP[t]*PCr[t])/(MMCKKd3*MMCKKic3)))/(1 + Cr[t]/MMCKKib3 + ATP[t]*(MMCKKia3^(-1) + Cr[t]/(MMCKKb3*MMCKKia3)) + PCr[t]/MMCKKid3 + ADP[t]*(MMCKKic3^(-1) + Cr[t]/(MMCKKib3*MMCKKic3) + PCr[t]/(MMCKKd3*MMCKKic3)))
MiCK = (IMS*((MiCKVf2*ATPi[t]*Cri[t])/(MiCKKb2*MiCKKia2) - (MiCKVb2*ADPi[t]*PCri[t])/(MiCKKd2*MiCKKic2)))/(1 + Cri[t]/MiCKKib2 + ATPi[t]*(MiCKKia2^(-1) + Cri[t]/(MiCKKb2*MiCKKia2)) + PCri[t]/MiCKKid2 + ADPi[t]*(MiCKKic2^(-1) + Cri[t]/(MiCKKib2*MiCKKic2) + PCri[t]/(MiCKKd2*MiCKKic2)))
OxPhos = (IMS*OxPhosV1*ADPi[t]*Pi[t])/(OxPhosKa1*OxPhosKb1*(1 + ADPi[t]/OxPhosKa1 + Pi[t]/OxPhosKb1 + (ADPi[t]*Pi[t])/(OxPhosKa1*OxPhosKb1)))
PCrdiffusion = -(CYT*PCrdiffusionk18*PCr[t]) + IMS*PCrdiffusionk18*PCri[t]
Pidiffusion = -(CYT*Pidiffusionk25*P[t]) + IMS*Pidiffusionk25*Pi[t]
end rate equations

begin parameters
OxPhosV1 = 4600.0
OxPhosKa1 = 800.0
OxPhosKb1 = 20.0
MiCKVf2 = 2658.0
MiCKKia2 = 750.0
MiCKKb2 = 5200.0
MiCKVb2 = 11160.0
MiCKKic2 = 204.8
MiCKKd2 = 500.0
MiCKKib2 = 28800.0
MiCKKid2 = 1600.0
MMCKVf3 = 6966.0
MMCKKia3 = 900.0
MMCKKb3 = 15500.0
MMCKVb3 = 29250.0
MMCKKic3 = 222.4
MMCKKd3 = 1670.0
MMCKKib3 = 34900.0
MMCKKid3 = 4730.0
ATPasev4 = 4600.0
Pidiffusionk25 = 18.4
Crdiffusionk16 = 14.6
ADPdiffusionk17 = 8.16
PCrdiffusionk18 = 14.6
ATPdiffusionk19 = 8.16
CYT = 0.75
IMS = 0.0625
end parameters

begin initial conditions
ADP[0] = ADPi
ADPi[0] = ADPii
ATP[0] = ATPi
ATPi[0] = ATPii
Cr[0] = Cri
Cri[0] = Crii
P[0] = Pi
PCr[0] = PCri
PCri[0] = PCrii
Pi[0] = Pii
end initial conditions

begin initial values
ADPi = 0.0
ADPii = 0.0
ATPi = 9700.0
ATPii = 0.0
Cri = 26000.0
Crii = 0.0
Pi = 0.0
PCri = 0.0
PCrii = 0.0
Pii = 32000.0
end initial values

begin assignment rules
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
ADP http://identifiers.org/chebi/CHEBI:16761
ADP http://identifiers.org/kegg.compound/C00008
ADPi http://identifiers.org/chebi/CHEBI:16761
ADPi http://identifiers.org/kegg.compound/C00008
ATP http://identifiers.org/chebi/CHEBI:15422
ATP http://identifiers.org/kegg.compound/C00002
ATPi http://identifiers.org/chebi/CHEBI:15422
ATPi http://identifiers.org/kegg.compound/C00002
Cr http://identifiers.org/chebi/CHEBI:16919
Cr http://identifiers.org/kegg.compound/C00300
Cri http://identifiers.org/chebi/CHEBI:16919
Cri http://identifiers.org/kegg.compound/C00300
P http://identifiers.org/chebi/CHEBI:18367
P http://identifiers.org/kegg.compound/C00009
PCr http://identifiers.org/chebi/CHEBI:17287
PCr http://identifiers.org/kegg.compound/C02305
PCri http://identifiers.org/chebi/CHEBI:17287
PCri http://identifiers.org/kegg.compound/C02305
Pi http://identifiers.org/chebi/CHEBI:18367
Pi http://identifiers.org/kegg.compound/C00009
end species annotations

begin units
end units