/
depression_temp.m
88 lines (75 loc) · 2.92 KB
/
depression_temp.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
function depression_temp
% ========================================================================
% This file is part of the Supplemental Codes of the manuscript
% entitled "A Kinetic Model Unifying Presynaptic Short-Term Facilitation
% and Depression" accepted by Journal of Compuational Neuroscience.
% (Manucript No. #JCNS583R2).
% Authors: Chuang-Chung J. Lee, Mihai Anton, Chi-Sang Poon, Gregory McRae
%
% Created by Chuang-Chung J. Lee
% Created in July '07.
% Latest modified in Oct. '08.
% ========================================================================
% ------------------------------------------------------------------------
% This function calculates the transient response of
% synapses under depression and plots both experimental & simulated
% results. The Probability of release and Ratio of realeasable vehsicle
% resposes are also output.
% Equations used: Eq. (1) - (7)
% Experimental Data source: Henrique von Gersdorff et al. (1997) J.
% Neurosci., 17(21):8137–8146
% Output: Figure 3C and 3D. Responses of depressing synapses. C, The
% transient EPSC caused by stimuli at 10 Hz in the rat calyx of Held
% synapse. D, The corresponding transient releasable vesicle ratio and
% release probability by model.
% ------------------------------------------------------------------------
close all
% ---Retrieve and plot the experimental and simulation data of EPSC-------
cd Output
load sim_depression;
load RPrel_depression;
sim_data=RPrel_depression;
cd ..
data1= [1,-8.04;2,-4.83;3,-3.66;4,-3.32;5,-2.68;6,-3.1;7,-2.37;8,-2.12;
9,-1.98;10,-2.37;11,-2.46;12,-2.46;13,-2.15;14,-2.35;15,-1.9;16,-2.29;
17,-1.84;18,-2.15;19,-2.35;20,-1.87;21,-2.15;22,-1.73;23,-2.04;24,-2.37;
25,-2.18;26,-2.1;27,-2.04;28,-1.84;29,-2.04;30,-1.82];
sim_data1=sim_depression;
% Get the local maximum of I(EPSC) for 10Hz
local_min1=zeros(30,2);
local_min1(:,1)=linspace(1,30,30)';
grid_n=length(sim_data1(:,1));
f=10;
T=1000/f;
delay_t=0;
for i=1:30
for n=1:grid_n
if (sim_data1(n,1)<=delay_t+(i-1)*T) && (sim_data1(n+1,1)>=delay_t+(i-1)*T)
lower_marker=n;
end
if (sim_data1(n,1)<=delay_t+i*T) && (sim_data1(n+1,1)>=delay_t+i*T)
upper_marker=n;
end
end
local_min1(i,2)=min(sim_data1(lower_marker:upper_marker,2));
end
figure
plot(data1(:,1)*1e2,data1(:,2),'o','MarkerEdgeColor','m','MarkerFaceColor','m','MarkerSize',5);
hold on
plot(sim_data1(:,1),sim_data1(:,2)*1e6,'b-','LineWidth',1.5);
xlabel('Time (ms)','FontSize', 14);
ylabel ('EPSC (nA)','FontSize', 14)
axis([0 3100 -10 1])
legend('experiment','model','Location','southeast');
legend('boxoff')
% ------------------Plot Prel and Rrel---------------------
figure
subplot(2,1,1)
plot(sim_data(:,1),sim_data(:,2),'LineWidth',1.3)
axis([0 3100 0 1])
ylabel ('R_r_e_l (t)','FontSize', 14);
subplot(2,1,2)
plot(sim_data(:,1),sim_data(:,3),'r','LineWidth',1.3)
axis([0 3100 0 1])
xlabel('Time (ms)','FontSize', 14);
ylabel ('P_r_e_l (t)','FontSize', 14);