A kinetic model unifying presynaptic short-term facilitation and depression (Lee et al. 2009)

 Download zip file 
Help downloading and running models
"... Here, we propose a unified theory of synaptic short-term plasticity based on realistic yet tractable and testable model descriptions of the underlying intracellular biochemical processes. Analysis of the model equations leads to a closed-form solution of the resonance frequency, a function of several critical biophysical parameters, as the single key indicator of the propensity for synaptic facilitation or depression under repetitive stimuli. This integrative model is supported by a broad range of transient and frequency response experimental data including those from facilitating, depressing or mixed-mode synapses. ... the model provides the reasons behind the switching behavior between facilitation and depression observed in experiments. ..."
1 . Lee CC, Anton M, Poon CS, McRae GJ (2009) A kinetic model unifying presynaptic short-term facilitation and depression. J Comput Neurosci 26:459-73 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Synapse;
Brain Region(s)/Organism:
Cell Type(s): Hippocampus CA1 pyramidal GLU cell; Cerebellum Purkinje GABA cell; Auditory nerve;
Channel(s): I Calcium;
Gap Junctions:
Receptor(s): AMPA;
Transmitter(s): Glutamate;
Simulation Environment: MATLAB; Simulink;
Model Concept(s): Synaptic Plasticity; Facilitation; Depression; Calcium dynamics;
Implementer(s): Lee, Chuang-Chung [chchlee at alum.mit.edu]; Anton, Mihai [mianton at mit.edu]; Poon, Chi-Sang [cpoon at mit.edu]; McRae, Gregory [mcrae at mit.edu];
Search NeuronDB for information about:  Hippocampus CA1 pyramidal GLU cell; Cerebellum Purkinje GABA cell; AMPA; I Calcium; Glutamate;
function mixture_freq

% ========================================================================
% 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 Oct. '07.
% Latest modified in Oct. '08.
% ========================================================================

% ------------------------------------------------------------------------
% This function calculates the steady state frequency response of 
% synapses under mixed effects and plots both experimental & simulated
% results.  The Probability of release and Ratio of realeasable vehsicle
% resposes are also output.  
% Equations used: Eq. (8) - (11)

% Experimental Data source: Henry Markram et al. (1998), differential 
% signaling via the same axon of neocortical pyramidal neurons, PNAS, 
% 95:5323–5328.

% Output: Figure 5A and 5B.  Synaptic plasticity in rat pyramidal neurons 
% under mixed effects. A, The EPSP as a function of frequency, with maximum
% and bandwidth labeled. B, The frequency response of releasable vesicle 
% ratio and release probability by model.
% ------------------------------------------------------------------------

% -----------Plot the experimental frequency response first---------------
close all
freq_ex=[0.86 5.43 10.29 15.43 20.57 25.43 30.29 35.14 40.29 45.14 50.00...
 54.57 59.71 64.86 69.43];
EPSPss=[1.31 2.65 4.31 5.77 6.25 6.13 6.13 5.57 4.35 4.90 4.98 3.52 3.88...
 3.28 3.08];

plot(freq_ex, EPSPss,'ko','MarkerFaceColor','k','Markersize', 4);
hold on

xlabel('Frequency (Hz)','FontSize', 14);
ylabel ('EPSP_s_s (mV)','FontSize', 14);

% ---------Then Calculate and plot the model frequency response ----------
freq=linspace(1,100,40); % frequency range in the units of Hz
freqkHz=freq/1000; % frequency in the units of kHz
Ca0=0.075e2; % calcium dynamics
krecov0=0.75E-2; % recovery rate constants
Prelmax=1; % release probability

[sim_EPSP Rrelss Prelss] = calc_EPSP(Ca0,KCa,Krel,krecovmax,Prelmax,nHill,freqkHz);
[rresonance] = calc_rreson(Ca0,KCa,Krel,krecovmax_Prelmax,nHill);
rresonkHz=rresonance/1e3; % (kHz)
freq_start0_left=10/1e3;  freq_start0_right=60/1e3;
freq_limit_left= fzero(@(freqkHz) calc_band(Ca0,KCa,Krel,krecovmax...
freq_limit_right= fzero(@(freqkHz) calc_band(Ca0,KCa,Krel,krecovmax...

axis([0 100 0 8])

[sim_EPSP] = calc_EPSP(Ca0,KCa,Krel,krecovmax,Prelmax,nHill,rresonkHz);

% Add additional labels regarding resonance frequency and bandwidth
plot ([rresonance rresonance],[sim_EPSP 0],'b--');
plot ([100 rresonance],[sim_EPSP sim_EPSP],'b--');
plot ([freq_limit_left freq_limit_left]*1e3,[sim_EPSP/2^0.5 0],'r--');
plot ([freq_limit_right freq_limit_right]*1e3,[sim_EPSP/2^0.5 0],'r--');
plot ([freq_limit_left 0.1]*1e3,[sim_EPSP/2^0.5 sim_EPSP/2^0.5],'r--');
text(80,5,'V_m_a_x-3dB','FontSize',14); % text(2,5)
text(87,7,'V_m_a_x','FontSize',14); % text(2,7)

% -----Plot Prel and Rrel-----
axis([0 100 0 1])
ylabel ('R_r_e_l (r)','FontSize', 14);

axis([0 100 0 1])
xlabel('Frequency (Hz)','FontSize', 14);
ylabel ('P_r_e_l (r)','FontSize', 14);

% -----------------Functions for calculating rresonance------------------
function [rresonance] = calc_rreson(Ca0,KCa,Krel,krecovmax_Prelmax,nHill)
rresonance = (-Ca0./KCa + (nHill.*(Krel./KCa).^nHill.*...
krecovmax_Prelmax).^(1./(nHill+1)))*1e3; %resonance frequency

% -------------------Functions for calculating sim_EPSP-------------------
function [sim_EPSP Rrelss Prelss] = calc_EPSP(Ca0,KCa,Krel,krecovmax,Prelmax,nHill,freqkHz)
sim_EPSP=sim_EPSP*34; % normalize EPSPs wrt the lowest freq

% -------------------Functions for calculating bandwidth-------------------
function [EPSPzero] = calc_band(Ca0,KCa,Krel,krecovmax,Prelmax,nHill,freqkHz,rresonkHz)
[sim_EPSP] = calc_EPSP(Ca0,KCa,Krel,krecovmax,Prelmax,nHill,freqkHz);
[sim_EPSP_reson] = calc_EPSP(Ca0,KCa,Krel,krecovmax,Prelmax,nHill,rresonkHz);
EPSPzero = sim_EPSP_reson/2^0.5-sim_EPSP;