%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Victor V. Matveev, Richard Bertram and Arthur Sherman (2011)
% Ca2+ Cooperativity of Exocytosis as a Measure of Ca2+ Channel Domain Overlap.
% Brain Research 1398: 126-138.
%
% This script produces data used in Fig. 4E
% Victor Matveev, January 8, 2011
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for M = 1 to 12 step 1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
cfile = "M" M "equidistConfig.par"
include cfile
file = "M" M "_ica" Imax "pA_B" BT "uM_" 1000*Dist "nm_" Dim1 "x" Dim2 ".dat"
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
verbose = 0
Dist0 = 0.03
Ca.bgr = 0.05
for Dist = Dist0 to 0.1 step 0.01
for conf = 1 to 2^M step 1
volume 0 L 0 L 0 L
L = 1
MID = L / 2
Dim1 = 70
Dim2 = 50
Imax = 0.05
BT = $6
grid Dim1 Dim1 Dim2
c0 = MID - Dist
c1 = MID + Dist
stretch x c0 c1
stretch y c0 c1
stretch z 0 0
stretch.factor = 1.05
ww = 0.004
Ca.D = 0.2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
T0 = 1; T1 = 1; T2 = 4
Run adaptive T0
Run adaptive T1
Run adaptive T2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
buffer B
B.total = BT
B.D = 0.05
B.KD = 1
B.kplus = 0.65
uptake = 0.004
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Felmy-Neher-Schenggenburger Ca-binding scheme (Neuron 37:801, 2003)
X(0) = 9.965526822168e-01;
X1(0) = 3.428235798835e-03;
X2(0) = 1.886953032458e-05;
X3(0) = 2.077213841252e-07;
X4(0) = 4.572912343658e-09
C := Ca[MID,MID,0]
dX/dt = - 5 C kon X + koff X1
dX1/dt = 5 kon C X - koff X1 - 4 C kon X1 + 2 koff b X2
dX2/dt = 4 kon C X1 - 2 koff b X2 - 3 C kon X2 + 3 koff b^2 X3
dX3/dt = 3 kon C X2 - 3 koff b^2 X3 - 2 C kon X3 + 4 koff b^3 X4
dX4/dt = 2 kon C X3 - 4 koff b^3 X4 - C kon X4 + 5 koff b^4 X5
X5 := 1 - X - X1 - X2 - X3 - X4
R := gamma X5
kon = 0.116; koff = 8.430; gamma = 6.96; b = 0.25
TM = T0 + T1 + T2
Tmax Rmax max R 0 TM
Cmax max C 0 TM
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if conf == 1
print file " "
endif
append file ind " " mult " " sumf " " Rmax " " Cmax " " Tmax
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
I = Imax pA