#kongas1



#reactions
ADP_diffusion:
 {1.0}ADPi = {1.0}ADP
  IMS * k1_7 * ADPi - CYT * k1_7 * ADP
ATP_diffusion:
 {1.0}ATPi = {1.0}ATP
  IMS * k1_9 * ATPi - CYT * k1_9 * ATP
ATPase:
 {1.0}ATP = {1.0}ADP + {1.0}P
  CYT * v_4 * ATP
Cr_diffusion:
 {1.0}Cri = {1.0}Cr
  IMS * k1_6 * Cri - CYT * k1_6 * Cr
MMCK:
 {1.0}ATP + {1.0}Cr = {1.0}PCr + {1.0}ADP
  CYT * ((Vf_3 * ATP * Cr / (Kia_3 * Kb_3) - Vb_3 * ADP * PCr / (Kic_3 * Kd_3)) / (1 + Cr / Kib_3 + PCr / Kid_3 + ATP * (1 / Kia_3 + Cr / (Kia_3 * Kb_3)) + ADP * (1 / Kic_3 + Cr / (Kic_3 * Kib_3) + PCr / (Kid_3 * (Kic_3 * Kd_3 / Kid_3)))))
MiCK:
 {1.0}ATPi + {1.0}Cri = {1.0}ADPi + {1.0}PCri
  IMS * ((Vf_2 * ATPi * Cri / (Kia_2 * Kb_2) - Vb_2 * ADPi * PCri / (Kic_2 * Kd_2)) / (1 + Cri / Kib_2 + PCri / Kid_2 + ATPi * (1 / Kia_2 + Cri / (Kia_2 * Kb_2)) + ADPi * (1 / Kic_2 + Cri / (Kic_2 * Kib_2) + PCri / (Kid_2 * (Kic_2 * Kd_2 / Kid_2)))))
OxPhos:
 {1.0}ADPi + {1.0}Pi = {1.0}ATPi
  IMS * (V_1 * ADPi * Pi / (Ka_1 * Kb_1 * (1 + ADPi / Ka_1 + Pi / Kb_1 + ADPi * Pi / (Ka_1 * Kb_1))))
PCr_diffusion:
 {1.0}PCri = {1.0}PCr
  IMS * k1_8 * PCri - CYT * k1_8 * PCr
Pi_diffusion:
 {1.0}Pi = {1.0}P
  IMS * k2_5 * Pi - CYT * k2_5 * P


#parameters
OxPhos_V_1 = 4600.0
OxPhos_Ka_1 = 800.0
OxPhos_Kb_1 = 20.0
MiCK_Vf_2 = 2658.0
MiCK_Kia_2 = 750.0
MiCK_Kb_2 = 5200.0
MiCK_Vb_2 = 11160.0
MiCK_Kic_2 = 204.8
MiCK_Kd_2 = 500.0
MiCK_Kib_2 = 28800.0
MiCK_Kid_2 = 1600.0
MMCK_Vf_3 = 6966.0
MMCK_Kia_3 = 900.0
MMCK_Kb_3 = 15500.0
MMCK_Vb_3 = 29250.0
MMCK_Kic_3 = 222.4
MMCK_Kd_3 = 1670.0
MMCK_Kib_3 = 34900.0
MMCK_Kid_3 = 4730.0
ATPase_v_4 = 4600.0
Pi_diffusion_k2_5 = 18.4
Cr_diffusion_k1_6 = 14.6
ADP_diffusion_k1_7 = 8.16
PCr_diffusion_k1_8 = 14.6
ATP_diffusion_k1_9 = 8.16
CYT = 0.75
IMS = 0.0625


#initial values
ADP = 0.0
ADPi = 0.0
ATP = 9700.0
ATPi = 0.0
Cr = 26000.0
Cri = 0.0
P = 0.0
PCr = 0.0
PCri = 0.0
Pi = 32000.0


#assignment rules