#koo4



#reactions
re10:
 {1.0}s11 = {1.0}s3
  0.5 * dot_q_inpass
re102:
 {1.0}s95 = {1.0}s93
  exp(unity - pow(s39 / tf, 0.35)) * normal * pow((s39 + tiny_num) / unimol, -0.65) * (unity - pow(s39 / tf, 0.35))
re103:
 {1.0}s94 = {1.0}s92
  exp(unity - pow(s39 / tf, 1.3)) * normal * pow(s39 / unimol, 0.3) * (unity - pow(s39 / tf, 1.3))
re104:
 {1.0}s91 + {1.0}s93 = {1.0}s110
  k105 * s91 * s93 - k_105 * s110
re105:
 {1.0}s110 = {1.0}s112
  kcat_src * s110 * s92 - k_6 * s112
re106:
 {1.0}s112 + {1.0}s111 = {1.0}s113
  k107 * s111 * s112 - k_107 * s113
re107:
 {1.0}s113 = {1.0}s114 + {1.0}s93
  k108 * s113 - k_108 * s93 * s114
re108:
 {1.0}s98 = {1.0}s91
  V10 * s98 / (K10 + s98)
re109:
 {1.0}s114 = {1.0}s98 + {1.0}s111
  k9 * s114
re11:
 {1.0}s13 = {1.0}s12
  unisec * 0.5
re110:
 {1.0}s115 = {1.0}s106
  s115 * kT
re111:
 {1.0}s108 = {1.0}s115
  ktr1 * s107 + ktr1k2 * s49
re112:
 {1.0}s103 = {1.0}s102
  k111 * s114 * s103 / (s103 + K111)
re113:
 {1.0}s102 = {1.0}s103
  V12 * s102 / (K12 + s102)
re114:
 {1.0}s97 = {1.0}s101
  k113 * s102 * s97 / (K113 + s97)
re115:
 {1.0}s101 = {1.0}s97
  V18 * s101 / (K18 + s101)
re116:
 {1.0}s96 = {1.0}s100
  k19 * s101 * s96 / (K19 + s96)
re117:
 {1.0}s100 = {1.0}s96
  V20 * s100 / (s100 + K20)
re118:
 {1.0}s99 = {1.0}s104
  k21 * s99 * s43 / (s99 + K21)
re119:
 {1.0}s104 = {1.0}s99
  V22 * s104 / (s104 + K22)
re12:
 {1.0}s3 = {1.0}s9
  0.5 * dot_Vp * pow(s3, 2) / (pow(K4, 2) + pow(s3, 2))
re120:
 {1.0}s106 = {1.0}s35
  kD * s106
re121:
 {1.0}s37 = {1.0}s51
  kP * s106 / (s106 + K30)
re122:
 {1.0}s38 = {1.0}s39
  unitime
re123:
 {1.0}s43 = {1.0}s100
  V37 * s43 / (s43 + K37)
re124:
 {1.0}s100 = {1.0}s43
  k38 * s100 * s101 / (K38 + s100)
re125:
 {1.0}s44 = {1.0}s104
  V39 * s44 / (K39 + s44)
re126:
 {1.0}s104 = {1.0}s44
  k40 * s104 * s43 / (K40 + s104)
re127:
 {1.0}s42 = {1.0}s107
  s44 * s42 * k43 / (K43 + s42)
re128:
 {1.0}s107 = {1.0}s42
  V44 * s107 / (s107 + K44)
re129:
 {1.0}s105 = {1.0}s49
  exp(tau * (delay - s39 / tc)) / (1 + 2 * exp(tau * (delay - s39 / tc)) + exp(2 * tau * (delay - s39 / tc))) / uc * 29.256
re131:
 {1.0}s51 = {1.0}s116
  kDD * s51
re132:
 {1.0}s45 = {1.0}s117 + {1.0}s47
  kDD * s45
re133:
 {1.0}s50 = {1.0}s117 + {1.0}s48
  kDD * s50
re134:
 {1.0}s66 = {1.0}s57 + {1.0}s118
  kDD * s66
re135:
 {1.0}s65 = {1.0}s118 + {1.0}s57
  kDD * s65
re136:
 {1.0}s61 = {1.0}s57 + {1.0}s48
  kDD * s61
re137:
 {1.0}s62 = {1.0}s57 + {1.0}s48
  kDD * s62
re138:
 {1.0}s58 = {1.0}s57 + {1.0}s47
  kDD * s58
re139:
 {1.0}s60 = {1.0}s57 + {1.0}s47
  kDD * s60
re3:
 {1.0}s4 = {1.0}s3
  -k6 * s3 * (B_T - s4) + k7 * s4
re37:
 {1.0}s51 + {1.0}s47 = {1.0}s45
  k15 * s47 * s51
re38:
 {1.0}s48 + {1.0}s51 = {1.0}s50
  k17 * s51 * s48 - k18 * s50
re4:
 {1.0}s5 = {1.0}s2
  k_CCE * (fracK * Cao / (K3 + Cao) - s2) * (s1 - s2)
re41:
 {1.0}s52 = {1.0}s48
  k11 * s3 * s52 - k12 * s48
re42:
 {1.0}s48 = {1.0}s47
  k13 * s3 * s48 - k14 * s47
re5:
 {1.0}s7 = {1.0}s6
  0.5 * k1 * (R_T - half * R_T * (exp(-s12 / tau_I) + exp(-s12 / tau_II) + (exp(-s12 / tau_I) - exp(-s12 / tau_II)) * (tau_I + tau_II) / (tau_I - tau_II))) * s3 / (K1 + s3)
re50:
 {1.0}s45 = {1.0}s50
  gam * k14 * s45 - k13 * s3 * s50
re51:
 {1.0}s45 + {1.0}s57 = {1.0}s58
  k90 * s45 * s57
re52:
 {1.0}s58 = {1.0}s61
  gam * k14 * s58 - k13 * s3 * s61
re53:
 {1.0}s61 = {1.0}s50 + {1.0}s57
  kr90 * s61
re54:
 {1.0}s60 = {1.0}s62
  gam * k14 * s60 - k13 * s3 * s62
re55:
 {1.0}s58 = {1.0}s60
  kp * s58 * s27 / (s58 + Kmp) - Vdp * s60 / (s60 + Kmdp)
re56:
 {1.0}s61 = {1.0}s62
  kp * s61 * s27 / (s61 + Kmp) - Vdp * s62 / (s62 + Kmdp)
re57:
 {1.0}s22 = {1.0}s20
  exp(unity - pow(s23 / tf, 1.8)) * normal * pow(s23 / unimol, 0.8) * (unity - pow(s23 / tf, 1.8))
re58:
 {1.0}s19 = {1.0}s17
  k58 * s19 * s20 / (Km58 + s19)
re59:
 {1.0}s17 = {1.0}s19
  k59 * s17 * s18 / (Km59 + s17)
re6:
 {1.0}s6 = {1.0}s8
  0.5 * k2 * s6
re60:
 {1.0}s17 + {1.0}s16 = {1.0}s28
  k60 * s17 * s16 - kr60 * s28
re61:
 {1.0}s28 = {1.0}s26
  k61 * s28 * s14 / (Km61 + s28)
re62:
 {1.0}s26 = {1.0}s27
  k62 * s26 * s24 / (Km62 + s26)
re63:
 {1.0}s26 = {1.0}s28
  k63 * s26 * s15 / (Km63 + s26)
re64:
 {1.0}s27 = {1.0}s26
  k64 * s27 * s15 / (Km64 + s27)
re65:
 {1.0}s27 = {1.0}s17 + {1.0}s16
  k65 * s27 * s15 / (Km65 + s27)
re66:
 {1.0}s25 = {1.0}s14
  k66 * s17 * s25
re67:
 {1.0}s14 = {1.0}s25
  s14 * k67
re68:
 {1.0}s21 = {1.0}s23
  unitime
re69:
 {1.0}s63 = {1.0}s64
  K_cam_no * (s45 + s58) + K_p_no * (s62 + s65) + K_pcam_no * s60
re7:
 {1.0}s2 = {1.0}s3
  0.5 * k3 * (k_CICR * s3 / (K_CICR + s3)) * pow(s6 / (K2 + s6), 3) * s2 - 0.5 * k4 * pow(s3 / (K3 + s3), 2) + 0.5 * k5 * s2 * s2
re70:
 {1.0}s65 = {1.0}s66
  Vdp * s65 / (s65 + Kmdp)
re71:
 {1.0}s66 = {1.0}s51 + {1.0}s57
  kr90 * s66
re72:
 {1.0}s62 = {1.0}s65 + {1.0}s48
  k18 * s62 - k17 * s65 * s48
re8:
 {1.0}s3 = {1.0}s9
  0.5 * dot_Vhi * pow(s3, 4) / (pow(K_hi, 4) + pow(s3, 4))
re9:
 {1.0}s3 = {1.0}s10
  0.5 * dot_Vex * s3 / (K5 + s3)


#parameters
B_T = 120000.0
Cao = 100.0
K1 = 0.0
K2 = 200.0
K3 = 150.0
K4 = 80.0
K5 = 321.0
K_CICR = 0.0
K_hi = 380.0
Kmdp = 20.0
Kmp = 5.0
R_T = 44000.0
Vc_Vs = 3.5
Vdp = 4.0
alp = 10.0
dot_Vex = 18330.0
dot_Vhi = 4760.0
dot_Vp = 1630.0
dot_q_inpass = 6000.0
dot_q_instim = 2500.0
fracK = 7071067.81
gam = 0.1
half = 0.5
k1 = 0.0012
k11 = 0.004
k12 = 10.3
k13 = 0.08
k14 = 1152.0
k15 = 0.015
k16 = 0.0
k17 = 0.00015
k18 = 1.5
k2 = 2.0
k3 = 6.64
k4 = 5000.0
k5 = 1e-10
k6 = 0.1
k7 = 300.0
k8 = 7.5e-05
k90 = 0.002
kDD = 9.45e-05
k_CCE = 0.0
k_CICR = 1.0
kp = 0.1
kr90 = 1.5
tau_I = 33.0
tau_II = 0.005
re11_unisec = 1.0
re57_normal = 0.907
re57_unity = 1.0
re57_unimol = 1.0
re57_tf = 15.0
re58_k58 = 0.2
re58_Km58 = 6170.0
re59_k59 = 7.5
re59_Km59 = 80.9
re60_k60 = 0.045
re60_kr60 = 0.089
re61_k61 = 20.0
re61_Km61 = 80000.0
re62_k62 = 20.0
re62_Km62 = 80000.0
re63_k63 = 0.037
re63_Km63 = 8800.0
re64_k64 = 0.04
re64_Km64 = 48000.0
re65_k65 = 0.163
re65_Km65 = 48000.0
re66_k66 = 0.0007
re67_k67 = 0.98
re68_unitime = 1.0
re69_K_cam_no = 17.0
re69_K_p_no = 5.0
re69_K_pcam_no = 17.0
re102_normal = 4.0
re102_tf = 60.0
re102_unity = 1.0
re102_unimol = 1.0
re102_tiny_num = 1e-06
re103_normal = 0.026
re103_tf = 540.0
re103_unity = 1.0
re103_unimol = 1.0
re104_k105 = 0.1
re104_k_105 = 1.0
re105_kcat_src = 8.33
re105_k_6 = 5.0
re106_k107 = 60.0
re106_k_107 = 546.0
re107_k108 = 2040.0
re107_k_108 = 15700.0
re108_V10 = 154.0
re108_K10 = 340.0
re109_k9 = 40.8
re110_kT = 0.001
re111_ktr1 = 0.00012
re111_ktr1k2 = 9e-06
re111_tr2 = 3e-06
re112_k111 = 0.222
re112_K111 = 0.181
re113_V12 = 0.289
re113_K12 = 0.0571
re114_k113 = 0.035
re114_K113 = 10.0
re115_K18 = 8.0
re115_V18 = 0.125
re116_k19 = 0.005
re116_K19 = 15.0
re117_K20 = 15.0
re117_V20 = 0.375
re118_k21 = 0.002
re118_K21 = 30.0
re119_K22 = 15.0
re119_V22 = 0.05
re120_kD = 2.8e-05
re121_kP = 0.02824
re121_K30 = 16.0
re122_unitime = 1.0
re123_V37 = 0.375
re123_K37 = 15.0
re124_k38 = 0.005
re124_K38 = 15.0
re125_V39 = 0.05
re125_K39 = 15.0
re126_k40 = 0.002
re126_K40 = 30.0
re127_K43 = 25.0
re127_k43 = 4e-05
re128_K44 = 5.0
re128_V44 = 0.002
re129_tau = 0.55
re129_tc = 3600.0
re129_uc = 3600.0
re129_delay = 5.0
re129_unity = 20.0
c1 = 1.0
c2 = 1.0
c3 = 1.0
default = 1.0


#initial values
s1 = 1500000.0
s10 = 0.0
s100 = 0.288
s101 = 1.486
s102 = 0.616
s103 = 119.384
s104 = 0.003
s105 = 0.0
s106 = 3.214
s107 = 0.0
s108 = 0.0
s11 = 0.0
s110 = 0.857
s111 = 3.23
s112 = 15.962
s113 = 5.577
s114 = 1.193
s115 = 0.09
s116 = 0.0
s117 = 0.0
s118 = 0.0
s119 = 0.0
s12 = 0.0
s13 = 0.0
s14 = 0.246
s15 = 150.0
s16 = 167.616
s17 = 0.345
s18 = 0.1
s19 = 6967.271
s2 = 2830000.0
s20 = 0.03
s21 = 0.0
s22 = 99.97
s23 = 0.0
s24 = 3.0
s25 = 999.754
s26 = 1.457
s27 = 1.723
s28 = 29.203
s3 = 117.2
s35 = 0.0
s37 = 10000.0
s38 = 0.0
s39 = 0.0
s4 = 3870.0
s42 = 50.0
s43 = 0.006
s44 = 0.0
s45 = 0.0415
s47 = 2.827
s48 = 347.52
s49 = 10.0
s5 = 0.0
s50 = 2.12
s51 = 34.98
s52 = 7635.36
s57 = 199987.0
s58 = 1.037
s6 = 0.0
s60 = 0.0089
s61 = 10.98
s62 = 0.106
s63 = 500000.0
s64 = 0.0
s65 = 0.643
s66 = 0.083
s7 = 0.0
s8 = 0.0
s9 = 0.0
s91 = 819.25
s92 = 18.0
s93 = 0.605
s94 = 72.0
s95 = 57.0
s96 = 299.706
s97 = 98.514
s98 = 157.162
s99 = 299.997


#assignment rules