begin name
bindschadler1
end name

begin reactions
InactivatedtoSCell1 {1.0}$EXTERNAL = {1.0}h1
InactivatedtoSCell2 {1.0}$EXTERNAL = {1.0}h2
JleakCell1 {1.0}$EXTERNAL = {1.0}c1
JleakCell2 {1.0}$EXTERNAL = {1.0}c2
JpumpCell1 {1.0}c1 = {1.0}$EXTERNAL
JpumpCell2 {1.0}c2 = {1.0}$EXTERNAL
JreceptorCell1 {1.0}$EXTERNAL = {1.0}c1
JreceptorCell2 {1.0}$EXTERNAL = {1.0}c2
OpentoInactivatedCell1 {1.0}h1 = {1.0}$EXTERNAL
OpentoInactivatedCell2 {1.0}h2 = {1.0}$EXTERNAL
diffusion {1.0}c2 = {1.0}c1
end reactions

begin rate equations
InactivatedtoSCell1 = compartment*Phi3c1*(1 - h1[t])
InactivatedtoSCell2 = compartment*Phi3c2*(1 - h2[t])
JleakCell1 = compartment*JleakCell1Jleak
JleakCell2 = compartment*JleakCell2Jleak
JpumpCell1 = (compartment*JpumpCell1Vp*c1[t]^2)/(JpumpCell1Kp^2 + c1[t]^2)
JpumpCell2 = (compartment*JpumpCell2Vp*c2[t]^2)/(JpumpCell2Kp^2 + c2[t]^2)
JreceptorCell1 = (compartment*JreceptorCell1kf*JreceptorCell1p^4*Phi1c1^4*h1[t]^4)/(JreceptorCell1p*Phi1c1 + Phiminus1c1)^4
JreceptorCell2 = (compartment*JreceptorCell2kf*JreceptorCell2p^4*Phi1c2^4*h2[t]^4)/(JreceptorCell2p*Phi1c2 + Phiminus1c2)^4
OpentoInactivatedCell1 = (compartment*OpentoInactivatedCell1p*Phi1c1*Phi2c1*h1[t])/(OpentoInactivatedCell1p*Phi1c1 + Phiminus1c1)
OpentoInactivatedCell2 = (compartment*OpentoInactivatedCell2p*Phi1c2*Phi2c2*h2[t])/(OpentoInactivatedCell2p*Phi1c2 + Phiminus1c2)
diffusion = compartment*diffusionD*(-c1[t] + c2[t])
end rate equations

begin 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
JreceptorCell1kf = 28.0
JreceptorCell1p = 0.2778
JpumpCell1Vp = 1.2
JpumpCell1Kp = 0.18
JleakCell1Jleak = 0.2
OpentoInactivatedCell1p = 0.2778
JreceptorCell2kf = 28.0
JreceptorCell2p = 0.2778
JpumpCell2Vp = 1.2
JpumpCell2Kp = 0.18
JleakCell2Jleak = 0.2
OpentoInactivatedCell2p = 0.2778
diffusionD = 0.01
compartment = 1.0
end parameters

begin initial conditions
c1[0] = c1i
c2[0] = c2i
h1[0] = h1i
h2[0] = h2i
end initial conditions

begin initial values
c1i = 0.3
c2i = 0.1
h1i = 0.8
h2i = 0.1
end initial values

begin assignment rules
Phi3c2 := k3/(R5 + c2[t])
Phi2c2 := (k2 + r4*c2[t])/(R3 + c2[t])
Phiminus1c2 := k1/(R3 + c2[t])
Phi1c2 := (r2*c2[t])/(R1 + c2[t])
Phi3c1 := k3/(R5 + c1[t])
Phi2c1 := (k2 + r4*c1[t])/(R3 + c1[t])
Phiminus1c1 := k1/(R3 + c1[t])
Phi1c1 := (r2*c1[t])/(R1 + c1[t])
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
Jpump_Cell1 http://identifiers.org/go/GO:0005388
Jpump_Cell2 http://identifiers.org/go/GO:0005388
Jreceptor_Cell1 http://identifiers.org/go/GO:0005220
Jreceptor_Cell1 http://identifiers.org/go/GO:0007204
Jreceptor_Cell2 http://identifiers.org/go/GO:0005220
Jreceptor_Cell2 http://identifiers.org/go/GO:0007204
end process annotations

begin species annotations
c1 http://identifiers.org/chebi/CHEBI:29108
c1 http://identifiers.org/kegg.compound/C00076
c2 http://identifiers.org/chebi/CHEBI:29108
c2 http://identifiers.org/kegg.compound/C00076
end species annotations

begin units
end units