Multiple mechanisms of short term plasticity at the calyx of Held (Hennig et al. 2008)

 Download zip file 
Help downloading and running models
This is a new model of the short-term dynamics of glutamatergic synaptic transmission, which incorporates multiple mechanisms acting at differing sites and across a range of different time scales (ms to tens of seconds). In the paper, we show that this model can accurately reproduce the experimentally measured time-course of short term depression across different stimulus frequencies at the calyx of Held. The model demonstrates how multiple forms of activity-dependent modulation of release probability and vesicle pool depletion interact, and shows how stimulus-history-dependent recovery from synaptic depression can arise from dynamics on multiple time scales.
1 . Hennig MH, Postlethwaite M, Forsythe ID, Graham BP (2008) Interactions between multiple sources of short-term plasticity during evoked and spontaneous activity at the rat calyx of Held. J Physiol 586:3129-46 [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): Medial Nucleus of the Trapezoid Body (MNTB) neuron;
Channel(s): I Calcium;
Gap Junctions:
Receptor(s): AMPA; mGluR;
Transmitter(s): Glutamate;
Simulation Environment: MATLAB;
Model Concept(s): Short-term Synaptic Plasticity;
Implementer(s): Hennig, Matthias H [mhhennig at];
Search NeuronDB for information about:  AMPA; mGluR; I Calcium; Glutamate;
% Interactions between multiple sources of short term plasticity
% during evoked and spontaneous activity at the rat calyx of Held
% J Physiol, 2008
% Matthias H. Hennig, Michael Postlethwaite, Ian D. Forsythe, Bruce
% P. Graham
% MHH:; BPG:
% This function contains the set of ODEs of the synapse model:
% y(1) is the 1st inactivated state (i1)
% y(2) is the calcium transient amplitude (c1)
% y(3) represents effects of  mGluR/AMPAR activation (b)
% y(4) is the 2nd inactivated state (i2)
% y(5) is the vesicle pool occupancy (n)
% y(6) is the rate of activity-dependent vesicle recruitment (ke)

function dy=synOde(t,y)

global gcairel gcairel2 gfrel gglrel gkd gke kek ke

dy = zeros(6,1);

dy(1) = -y(1)/gcairel+y(4)/gcairel2;
dy(2) = -1/gfrel*(y(2) - 1+y(1)+y(3)+y(4) );
dy(3) = -y(3)/gglrel;
dy(4) = -y(4)/gcairel2;
dy(5) = (ke*y(6)+1/gkd)*(1-y(5));
dy(6) = -y(6)/kek;