CA1 pyramidal cells, basket cells, ripples (Malerba et al 2016)

 Download zip file 
Help downloading and running models
Accession:188977
Model of CA1 pyramidal layer Ripple activity, triggered when receiving current input (to represent CA3 sharp-waves). Cells are Adaptive-Exponential Integrate and Fire neurons, receiving independent OU noise.
Reference:
1 . Malerba P, Krishnan GP, Fellous JM, Bazhenov M (2016) Hippocampal CA1 Ripples as Inhibitory Transients. PLoS Comput Biol 12:e1004880 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism:
Cell Type(s): Hippocampus CA1 pyramidal GLU cell; Hippocampus CA1 PV+ fast-firing interneuron;
Channel(s):
Gap Junctions:
Receptor(s): GabaA; AMPA;
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Oscillations; Sleep; Brain Rhythms;
Implementer(s):
Search NeuronDB for information about:  Hippocampus CA1 pyramidal GLU cell; GabaA; AMPA;
%% summary plots - samples of default case
close all
clear
s = 'example_1.mat';
load(s);
rpick = 2; % which ripple to plot
% inputs, rastergram, lfp and filteredLFP
tout = 0:0.001:T;
[ripples,spcount,recruit,filtLFP] = CountRipples(T,lfp,tsp_E,tsp_I,NE,NI,inpseq);
rpt = ripples.time(rpick); %[=]s
disp(rpt)
rpl = ripples.length(rpick)/1000; %[=]s

figure(1)
plot(tout,veg.E,'-k',tout,veg.I,'-r')
title('Example voltage')
ylabel('mV')
xlabel('time [=] s');
xlim([xl xr]);

xl = rpt-0.02;
xr = rpt+rpl+0.02;
k1 = find(tout>=xl,1);
k2 = find(tout>=xr,1);
X = k1:k2;
figure(2)
subplot(311)
plot(tout(X),inp.Itrace(1,X),'-r',tout(X),inp.Etrace(1,X),'-k');
xlim([xl xr]);
title('current input from CA3');
subplot(312)
plot(tsp_E.times,tsp_E.celln+NI,'+k',tsp_I.times,tsp_I.celln,'+r');
xlim([xl xr]);
subplot(313)
tk = tout(X);
l1 = lfp(X);
l2 = filtLFP(X);
%plot(tk,[l1; l2]')
% set(gca,'xlimmode','manual');
xlim([xl xr]);
plotyy(tout(X),lfp(X),tout(X),filtLFP(X))
title('LFP');
xlabel('time [=] s');

tx = 0:0.0005:T;
figure(3)
subplot(311)
hist(tsp_E.times,tx);
xlim([3.5 3.6]);
subplot(312)
hist(tsp_I.times,tx);
xlim([3.5 3.6]);
subplot(313)
plot(tout,filtLFP);
xlim([3.5 3.6]);