CA1 pyramidal neuron: dendritic Ca2+ inhibition (Muellner et al. 2015)

 Download zip file 
Help downloading and running models
Accession:206244
In our experimental study, we combined paired patch-clamp recordings and two-photon Ca2+ imaging to quantify inhibition exerted by individual GABAergic contacts on hippocampal pyramidal cell dendrites. We observed that Ca2+ transients from back-propagating action potentials were significantly reduced during simultaneous activation of individual nearby GABAergic synapses. To simulate dendritic Ca2+ inhibition by individual GABAergic synapses, we employed a multi-compartmental CA1 pyramidal cell model with detailed morphology, voltage-gated channel distributions, and calcium dynamics, based with modifications on the model of Poirazi et al., 2003, modelDB accession # 20212.
Reference:
1 . Müllner FE, Wierenga CJ, Bonhoeffer T (2015) Precision of Inhibition: Dendritic Inhibition by Individual GABAergic Synapses on Hippocampal Pyramidal Cells Is Confined in Space and Time. Neuron 87:576-89 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism: Hippocampus;
Cell Type(s): Hippocampus CA1 pyramidal GLU cell;
Channel(s): I Calcium; I Sodium; I Potassium; I h;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s): Gaba;
Simulation Environment: NEURON;
Model Concept(s): Action Potentials; Dendritic Action Potentials; Active Dendrites; Calcium dynamics;
Implementer(s): Muellner, Fiona E [fiona.muellner at gmail.com];
Search NeuronDB for information about:  Hippocampus CA1 pyramidal GLU cell; I h; I Sodium; I Calcium; I Potassium; Gaba;
% Plot spatiotemporal profile of Ca2+-inhibition for a selected simulation
% ---
% Fiona Müllner, MPI Neurobiology
% Email: fiona.muellner at gmail.com

clearvars -except rootfolder
%  Choose file:
cd(sprintf('%s\\CA1_multi\\experiment\\regular',rootfolder))
load('simulationRes_2.mat')
% Choose conductance:
cond1 = 4;

figure
distv1 = [distv(1:n/2)',distv(1:n/2)'];
posv1 = posv;
pos1 = min(3,length(posv1));
fbp=find(branchid(1:n/2-1)>1,1,'first')+1;
if isempty(fbp)
    chan=n/2+1:n;
else
    chan=n/2+1:n/2+fbp-1;
end
dist = distv1(chan);
[~,s]=sort(dist);

Inh=1-(squeeze(M(pos1,cond1,1:length(timev),chan(s)))./squeeze(repmat(M(pos1,1,1,chan(s)),[1,1,nj,1])))';
if size(Inh,1)==1
    Inh=Inh';
end
newdist=floor(min(dist(s)-posv1(pos1))):ceil(max(dist(s)-posv1(pos1)));
Int=zeros(length(newdist),size(Inh,2));
for i = 1:size(Inh,2)
    for j = 1:length(newdist);
        f = find(dist(s)-posv1(pos1)>newdist(j),1,'first');
        if f == 1
            Int(j,i)=Inh(f,i);
        elseif isempty(f)
            Int(j,i)=Inh(end,i);
        elseif f>1 && dist(s(f))-posv1(pos1)-newdist(j)<=newdist(j)-dist(s(f-1))+posv1(pos1)
            Int(j,i)=Inh(f,i);
        elseif f>1 && dist(s(f))-posv1(pos1)-newdist(j)>newdist(j)-dist(s(f-1))+posv1(pos1)
            Int(j,i)=Inh(f-1,i);
        end
    end
end
iptsetpref('ImshowAxesVisible','on');
if size(Inh,2)==1
    imshow(Int','XData',newdist/20,'InitialMagnification','fit');
else
    imshow(Int','XData',newdist/20,'YData',timev-stimstart,'InitialMagnification','fit');
end
set(gca,'clim',[-1 1])
colormap(jet);
xlabel('Distance [20 um]')
ylabel('Timing [ms]')
colorbar
set(gca,'clim',[-0.3 0.5])

Loading data, please wait...