#benary2



#reactions
reaction_1:
 $pool = {1.0}Dsha
  Compartment * k1 * piecewise(0, lt(time, 100000), exp(-ModelValue_27 * (time - 100000))) * Dshi
reaction_10:
 {1.0}APC__Axin__GSK3_b_catenin__HOS = {1.0}APC__Axin__GSK3 + {1.0}b_catenin__ub + {1.0}HOS
  Compartment * k10 * APC__Axin__GSK3_b_catenin__HOS
reaction_10b:
 {1.0}APC__Axin__GSK3_b_catenin__FWD1 = {1.0}APC__Axin__GSK3 + {1.0}b_catenin__ub + {1.0}FWD1
  Compartment * k10b * APC__Axin__GSK3_b_catenin__FWD1
reaction_11:
 {1.0}b_catenin__ub = $pool
  Compartment * k11 * b_catenin__ub
reaction_12:
 $pool = {1.0}b_catenin
  Compartment * Constant_flux__irreversible(v12)
reaction_13:
 {1.0}b_catenin = $pool
  Compartment * k13 * b_catenin
reaction_14:
 $pool = {1.0}Axin
  Compartment * Constant_flux__irreversible(v14)
reaction_15:
 {1.0}Axin = $pool
  Compartment * k15 * Axin
reaction_16:
 {1.0}TCF + {1.0}b_catenin = {1.0}b_catenin_TCF
  Compartment * (k16 * TCF * b_catenin - k16r * b_catenin_TCF)
reaction_17:
 {1.0}APC + {1.0}b_catenin = $pool
  Compartment * (k17 * APC * b_catenin - k17r * APC_b_catenin)
reaction_18:
 $pool = {1.0}HOSmRNA
  Compartment * Rate_Law_for_reaction_18(Khos, ki, vmax18, b_catenin_TCF)
reaction_19:
 {1.0}HOSmRNA = $pool
  Compartment * k19 * HOSmRNA
reaction_2:
 {1.0}Dsha = $pool
  Compartment * k2 * Dsha
reaction_20:
 $pool = {1.0}HOS
  Compartment * Rate_Law_for_reaction_20_25_29(k20, HOSmRNA)
reaction_21:
 {1.0}HOS = $pool
  Compartment * k21 * HOS
reaction_22:
 {1.0}HOS = {1.0}APC__Axin__GSK3_b_catenin__HOS
  Compartment * k22 * HOS * APC__Axin__GSK3_b_catenin_0
reaction_22b:
 {1.0}FWD1 = {1.0}APC__Axin__GSK3_b_catenin__FWD1
  Compartment * k22b * FWD1 * APC__Axin__GSK3_b_catenin_0
reaction_23:
 $pool = {1.0}CRD_BPmRNA
  Compartment * Rate_Law_for_reaction_23(vmax23, b_catenin_TCF, Kcrd_bp, vbasal)
reaction_24:
 {1.0}CRD_BPmRNA = $pool
  Compartment * k24 * CRD_BPmRNA
reaction_25:
 $pool = {1.0}CRD_BP
  Compartment * Rate_Law_for_reaction_20_25_29(k25, CRD_BPmRNA)
reaction_26:
 {1.0}CRD_BP = $pool
  Compartment * k26 * CRD_BP
reaction_27:
 {1.0}FWD1mRNA = $pool
  Compartment * Rate_Law_for_reaction_27(k27, Kfwd1, CRD_BP, ki2, FWD1mRNA)
reaction_28:
 $pool = {1.0}FWD1mRNA
  Compartment * Constant_flux__irreversible(v28)
reaction_29:
 $pool = {1.0}FWD1
  Compartment * Rate_Law_for_reaction_20_25_29(k29, FWD1mRNA)
reaction_3:
 {1.0}APC_Axin_GSK3 = {1.0}APC_Axin + {1.0}GSK3
  Compartment * Rate_Law_for_reaction_3(Dsha, APC_Axin_GSK3, k3)
reaction_30:
 {1.0}FWD1 = $pool
  Compartment * k30 * FWD1
reaction_31:
 $pool = {1.0}TCF
  Compartment * Constant_flux__irreversible(v31)
reaction_32:
 {1.0}TCF = $pool
  Compartment * k32 * TCF
reaction_4:
 {1.0}APC_Axin_GSK3 = {1.0}APC__Axin__GSK3
  Compartment * k4 * APC_Axin_GSK3
reaction_5:
 {1.0}APC__Axin__GSK3 = {1.0}APC_Axin_GSK3
  Compartment * k5 * APC__Axin__GSK3
reaction_6:
 {1.0}GSK3 + {1.0}APC_Axin = {1.0}APC_Axin_GSK3
  Compartment * (k6 * GSK3 * APC_Axin - k6r * APC_Axin_GSK3)
reaction_7:
 {1.0}APC + {1.0}Axin = {1.0}APC_Axin
  Compartment * (k7 * APC * Axin - k7r * APC_Axin)
reaction_8:
 {1.0}APC__Axin__GSK3 + {1.0}b_catenin = {1.0}APC__Axin__GSK3_b_catenin
  Compartment * (k8 * APC__Axin__GSK3 * b_catenin - k8r * APC__Axin__GSK3_b_catenin)
reaction_9:
 {1.0}APC__Axin__GSK3_b_catenin = $pool
  Compartment * k9 * APC__Axin__GSK3_b_catenin


#parameters
APCtotal = 100.0
Dshtotal = 100.0
GSKtotal = 50.0
Kcrd_bp = 585.0
Kfwd1 = 3e-06
Khos = 1.0
ModelValue_27 = lambda
k1 = 0.182
k10 = 206.0
k10b = 206.0
k11 = 0.417
k13 = 0.000257
k15 = 0.167
k16 = 1.0
k16r = 30.0
k17 = 1.0
k17r = 1200.0
k19 = 0.00048
k2 = 0.0182
k20 = 30.0
k21 = 0.038
k22 = 206.0
k22b = 206.0
k24 = 0.00246
k25 = 0.013
k26 = 7.35e-05
k27 = 0.003
k29 = 30.0
k3 = 0.05
k30 = 0.0058
k32 = 0.0028
k4 = 0.267
k5 = 0.133
k6 = 0.0909
k6r = 0.909
k7 = 1.0
k7r = 2000.0
k8 = 206.0
k8r = 123394.0
k9 = 206.0
ki = 1.667
ki2 = 1.667
lambda = 0.05
v12 = 0.423
v14 = 8.22e-05
v28 = 1e-06
v31 = 0.0228
vbasal = 9e-08
vmax18 = 1e-06
vmax23 = 4.22e-06
Compartment = 1.0


#initial values
APC = 97.89681
APC_Axin = 0.00096373
APC_Axin_GSK3 = 0.0048149
APC__Axin__GSK3 = 0.0096659
APC__Axin__GSK3_b_catenin = 0.0020218
APC__Axin__GSK3_b_catenin__FWD1 = 0.00013859
APC__Axin__GSK3_b_catenin__HOS = 0.0018833
Axin = 0.00049222
CRD_BP = 0.0065127
CRD_BPmRNA = 3.68215e-05
Dsha = 0
FWD1 = 0.006741
FWD1mRNA = 1.30327e-06
GSK3 = 49.96092
HOS = 0.091598
HOSmRNA = 0.00011602
TCF = 8.14286
b_catenin = 25.28956
b_catenin_TCF = 6.86431
b_catenin__ub = 0.9988


#assignment rules
!F APC_b_catenin = APCtotal - GSKtotal + GSK3 - APC_Axin - APC
!F Dshi = Dshtotal - Dsha
!F APC__Axin__GSK3_b_catenin_0 = GSKtotal - (APC_Axin_GSK3 + APC__Axin__GSK3 + GSK3 + APC__Axin__GSK3_b_catenin + APC__Axin__GSK3_b_catenin__HOS + APC__Axin__GSK3_b_catenin__FWD1)