Calyx of Held, short term plasticity (Yang Z et al. 2009)

 Download zip file 
Help downloading and running models
Accession:118554
This model investigates mechanisms contributing to short term plasticity at the calyx of Held, a giant glutamatergic synapse in the mammalian brainstem auditory system. It is a stochastic version of the model described in: Hennig, M., Postlethwaite, M., Forsythe, I.D. and Graham, B.P. (2007). A biophysical model of short-term plasticity at the calyx of Held. Neurocomputing, 70:1626-1629. This version introduces stochastic vesicle recycling and release. It has been used to investigate the information transmission properties of this synapse, as detailed in: Yang, Z., Hennig, M., Postlethwaite, M., Forsythe, I.D. and Graham, B.P. (2008). Wide-band information transmission at the calyx of Held. Neural Computation, 21(4):991-1018.
References:
1 . Yang Z, Hennig MH, Postlethwaite M, Forsythe ID, Graham BP (2009) Wide-band information transmission at the calyx of Held. Neural Comput 21:991-1017 [PubMed]
2 . Hennig MH, Postlethwaite M, Forsythe ID, Graham BP (2007) A biophysical model of short-term plasticity at the calyx of Held, Neurocomputing 70(12):1626-1629
Model Information (Click on a link to find other models with that property)
Model Type: Synapse;
Brain Region(s)/Organism: Auditory brainstem;
Cell Type(s): Medial Nucleus of the Trapezoid Body (MNTB) neuron;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s): Glutamate;
Simulation Environment: MATLAB;
Model Concept(s): Short-term Synaptic Plasticity; Vestibular;
Implementer(s): Hennig, Matthias H [mhhennig at gmail.com];
Search NeuronDB for information about:  Glutamate;
function [spikes, stimes, isi] = inhreg(t, dt, f)
% function [spikes, stimes isi] = inhreg(t, dt, f)
% Inhomogenous regular distributed ISIs
% t - time vector
% dt - time step
% f - instantaneous rate vector (per timebase)
% Basic assumptions are:
% (1) constant rate (frequency) over a time step
% (2) only a single arrival possible in a time step
% (so time step should be small relative to the rate of change in
% frequency and arrival rate)
% BPG 14-1-08

spikes=zeros(1,length(t));    
spikes(1)=0;
tp=0;   % index of previous spike
for i=2:length(t)
    cisi=1/f(i);    % current ISI (seconds)
    if (i-tp)*dt >= cisi
        spikes(i)=1;
        tp=i;       % index of new spike
    end;
end;
stimes=t(spikes==1);
isi=stimes(2:length(stimes))-stimes(1:length(stimes)-1);

Loading data, please wait...