Encoding and retrieval in a model of the hippocampal CA1 microcircuit (Cutsuridis et al. 2009)

 Download zip file 
Help downloading and running models
Accession:123815
This NEURON code implements a small network model (100 pyramidal cells and 4 types of inhibitory interneuron) of storage and recall of patterns in the CA1 region of the mammalian hippocampus. Patterns of PC activity are stored either by a predefined weight matrix generated by Hebbian learning, or by STDP at CA3 Schaffer collateral AMPA synapses.
Reference:
1 . Cutsuridis V, Cobb S, Graham BP (2010) Encoding and retrieval in a model of the hippocampal CA1 microcircuit. Hippocampus 20:423-46 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Hippocampus;
Cell Type(s): Hippocampus CA1 pyramidal GLU cell; Hippocampus CA1 basket cell;
Channel(s):
Gap Junctions:
Receptor(s): GabaA; AMPA; NMDA;
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Pattern Recognition; Activity Patterns; Temporal Pattern Generation; Learning; STDP; Connectivity matrix; Storage/recall;
Implementer(s): Graham, Bruce [B.Graham at cs.stir.ac.uk]; Cutsuridis, Vassilis [vcutsuridis at gmail.com];
Search NeuronDB for information about:  Hippocampus CA1 pyramidal GLU cell; GabaA; AMPA; NMDA;
% Plot voltage traces from example PC and INs
% Summary diagram for Hippocampus paper
% BPG 23-1-09

FSTEM = 'HAM_P5R1_';   % data file stem

dt = 0.025;

figure;
ms=8;
lw=1;
nr = 5;

STIME = 600;
ETIME = 1100;
VMIN = -90;
VMAX = 50;

subplot(nr,1,1);
FV = [FSTEM 'pvsoma.dat'];   % voltage file
v = load(FV);  % load spike times
t = (0:length(v)-1)*dt;       % extract times
hold on;
plot(t, v, 'k-');   % voltage trace
title('(a) Pattern pyramidal cell');
ylabel('V (mV)');
axis([STIME ETIME VMIN+10 VMAX-30]);

subplot(nr,1,2);
FV = [FSTEM 'AAC.dat'];   % voltage file
v = load(FV);  % load spike times
t = (0:length(v)-1)*dt;       % extract times
hold on;
plot(t, v, 'k-');   % voltage trace
title('(b) Axo-axonic cell');
ylabel('V (mV)');
axis([STIME ETIME VMIN VMAX]);

subplot(nr,1,3);
FV = [FSTEM 'BC.dat'];   % voltage file
v = load(FV);  % load spike times
t = (0:length(v)-1)*dt;       % extract times
hold on;
plot(t, v, 'k-');   % voltage trace
title('(c) Basket cell');
ylabel('V (mV)');
axis([STIME ETIME VMIN VMAX]);

subplot(nr,1,4);
FV = [FSTEM 'BSC.dat'];   % voltage file
v = load(FV);  % load spike times
t = (0:length(v)-1)*dt;       % extract times
hold on;
plot(t, v, 'k-');   % voltage trace
title('(d) Bistratified cell');
ylabel('V (mV)');
axis([STIME ETIME VMIN VMAX]);

subplot(nr,1,5);
FV = [FSTEM 'OLM.dat'];   % voltage file
v = load(FV);  % load spike times
t = (0:length(v)-1)*dt;       % extract times
hold on;
plot(t, v, 'k-');   % voltage trace
title('(e) OLM cell');
ylabel('V (mV)');
xlabel('Time (msecs)');
axis([STIME ETIME VMIN VMAX]);

print('-dpng', ['Images/' FSTEM 'v']);

    

Loading data, please wait...