#bindschadler1



#reactions
Inactivated_to_S_Cell1:
 $pool = {1.0}h1
  compartment * Phi3_c1 * (1 - h1)
Inactivated_to_S_Cell2:
 $pool = {1.0}h2
  compartment * Phi3_c2 * (1 - h2)
Jleak_Cell1:
 $pool = {1.0}c1
  compartment * Jleak
Jleak_Cell2:
 $pool = {1.0}c2
  compartment * Jleak
Jpump_Cell1:
 {1.0}c1 = $pool
  compartment * (Vp * pow(c1, 2) / (pow(Kp, 2) + pow(c1, 2)))
Jpump_Cell2:
 {1.0}c2 = $pool
  compartment * Vp * pow(c2, 2) / (pow(Kp, 2) + pow(c2, 2))
Jreceptor_Cell1:
 $pool = {1.0}c1
  compartment * kf * pow(p * h1 * Phi1_c1 / (Phi1_c1 * p + Phi_minus1_c1), 4)
Jreceptor_Cell2:
 $pool = {1.0}c2
  compartment * kf * pow(p * h2 * Phi1_c2 / (Phi1_c2 * p + Phi_minus1_c2), 4)
Open_to_Inactivated_Cell1:
 {1.0}h1 = $pool
  compartment * Phi1_c1 * Phi2_c1 * h1 * p / (Phi1_c1 * p + Phi_minus1_c1)
Open_to_Inactivated_Cell2:
 {1.0}h2 = $pool
  compartment * Phi1_c2 * Phi2_c2 * h2 * p / (Phi1_c2 * p + Phi_minus1_c2)
diffusion:
 {1.0}c2 = {1.0}c1
  compartment * D * (c2 - c1)


#parameters
R1 = 6.0
R3 = 50.0
R5 = 1.6
k1 = 44.0
k2 = 26.5
k3 = 1.6
r2 = 100.0
r4 = 20.0
Jreceptor_Cell1_kf = 28.0
Jreceptor_Cell1_p = 0.2778
Jpump_Cell1_Vp = 1.2
Jpump_Cell1_Kp = 0.18
Jleak_Cell1_Jleak = 0.2
Open_to_Inactivated_Cell1_p = 0.2778
Jreceptor_Cell2_kf = 28.0
Jreceptor_Cell2_p = 0.2778
Jpump_Cell2_Vp = 1.2
Jpump_Cell2_Kp = 0.18
Jleak_Cell2_Jleak = 0.2
Open_to_Inactivated_Cell2_p = 0.2778
diffusion_D = 0.01
compartment = 1.0


#initial values
c1 = 0.3
c2 = 0.1
h1 = 0.8
h2 = 0.1


#assignment rules
!F Phi3_c2 = k3 / (R5 + c2)
!F Phi2_c2 = (k2 + r4 * c2) / (R3 + c2)
!F Phi_minus1_c2 = k1 / (R3 + c2)
!F Phi1_c2 = r2 * c2 / (R1 + c2)
!F Phi3_c1 = k3 / (R5 + c1)
!F Phi2_c1 = (k2 + r4 * c1) / (R3 + c1)
!F Phi_minus1_c1 = k1 / (R3 + c1)
!F Phi1_c1 = r2 * c1 / (R1 + c1)