#model = 'combined'
#model = 'kpr'
model = 'two_ligand_kpr'
#for model in ['mode_1', 'mode_2', 'combined', 'kpr']:
for model in [model]:
print(model)
num_traj = 20000
num_steps = 1000
init_bound = 0.0
# model specification
#params = DEFAULT_PARAMS
params = Params(k_on=1.0, k_p=20, k_f=1.0)
params.c1 = 10.0
params.c2 = 4.0
params.k_off_1 = 2.0
params.k_off_2 = 1.0
# Steady state receptor occupancies
D1 = (params.k_off_1 * params.k_off_2 +
params.c2 * params.k_off_1 * params.k_on +
params.c1 * params.k_off_2 * params.k_on)
P0 = params.k_off_1 * params.k_off_2 / D1
P1 = params.c1 * params.k_off_1 * params.k_off_2 * params.k_on / (
(params.k_f + params.k_off_1) * D1)
P2 = params.c1 * params.k_f * params.k_off_2 * params.k_on / (
(params.k_f + params.k_off_1) * D1)
P3 = params.c2 * params.k_off_1 * params.k_off_2 * params.k_on / (
(params.k_f + params.k_off_2) * D1)
P4 = 1 - P0 - P1 - P2 - P3
