Learning intrinsic excitability in Medium Spiny Neurons (Scheler 2014)

 Download zip file 
Help downloading and running models
Accession:155131
"We present an unsupervised, local activation-dependent learning rule for intrinsic plasticity (IP) which affects the composition of ion channel conductances for single neurons in a use-dependent way. We use a single-compartment conductance-based model for medium spiny striatal neurons in order to show the effects of parameterization of individual ion channels on the neuronal membrane potential-curent relationship (activation function). We show that parameter changes within the physiological ranges are sufficient to create an ensemble of neurons with significantly different activation functions. ... "
Reference:
1 . Scheler G (2014) Learning intrinsic excitability in medium spiny neurons F1000Research 2:88 [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: Striatum;
Cell Type(s): Neostriatum medium spiny direct pathway GABA cell; Neostriatum medium spiny indirect pathway GABA cell;
Channel(s): I A; I K; I h; I K,Ca; I Calcium; I A, slow; I Cl, leak; I Ca,p;
Gap Junctions:
Receptor(s): GabaA; AMPA; NMDA;
Gene(s): Kv4.2 KCND2; Kv1.1 KCNA1; Kv1.2 KCNA2; Kv4.3 KCND3; Kv1.4 KCNA4; Kv1.3 KCNA3; Kv1.5 KCNA5; Kv3.3 KCNC3; Cav3.2 CACNA1H; Cav3.1 CACNA1G; Cav3.3 CACNA1I; Cav1.3 CACNA1D; Cav1.1 CACNA1S; Cav1.2 CACNA1C; KCa2.1 KCNN1; Kv2.1 KCNB1; Kv3.1 KCNC1; HCN Cnga1; Cav2.1 CACNA1A; Cav2.2 CACNA1B; KCa2.2 KCNN2; Kv1.9 Kv7.1 KCNQ1; IRK; NR2A GRIN2A; NR2B GRIN2B; Kv3.4 KCNC4; Kv4.1 KCND1;
Transmitter(s): Gaba; Glutamate; Ions;
Simulation Environment: MATLAB;
Model Concept(s): Intrinsic plasticity;
Implementer(s): Schumann, Johann [johann.schumann at gmail.com];
Search NeuronDB for information about:  Neostriatum medium spiny direct pathway GABA cell; Neostriatum medium spiny indirect pathway GABA cell; GabaA; AMPA; NMDA; I A; I K; I h; I K,Ca; I Calcium; I A, slow; I Cl, leak; I Ca,p; Gaba; Glutamate; Ions;
%
%	FN = filename to print (w/o extension)
%	sim: = sim structure 
%	nn_inputs = NN inputs
%	LW = line width for plotting
%	nn_pars: not used
%	off	temporal offset
%
%	$Revision:$
%
function plot_sr_conductIG(FN, sim, nn_inputs, LW, nn_pars, off)

[m_isi, s_isi] =calc_isi_plot(sim, off, 0);

m_isi
s_isi

N =sim.N_nn;

colors=['r','g','b'];
%colors=[{'k-'},{'k--'},{'k:'}];
%colors=['kk-','k--','kk:'];
colors=['k- ';'k-.';'k- '];
%colors(1,:)
%colors(2,:)
%colors(3,:)
LWs=[1,1,2]

figure;
i=1;
for nn=1,
	subplot(3,1,i);
	vm=reshape(sim.instrument.allvm(1,nn,1:sim.T_upd),sim.T_upd,1);
	plot(vm(off:end), 'k-', 'Linewidth',LWs(1));
	hold on;
	vm=reshape(sim.instrument.allvm(1,nn+1,1:sim.T_upd),sim.T_upd,1);
	plot(vm(off:end), 'k-.' , 'Linewidth',LWs(2));
	vm=reshape(sim.instrument.allvm(1,nn+2,1:sim.T_upd),sim.T_upd,1);
	plot(vm(off:end), 'k-', 'Linewidth',LWs(3));
	axis([0,sim.T_upd-off,-90,40]);
	set(gca,'Visible','off');
	subplot(3,1,i+2);
	plot(-nn_inputs(nn,off:end), 'Linewidth',LW);
	hold on;

	vm=reshape(sim.instrument.allvm(1,nn,1:sim.T_upd),sim.T_upd,1);
	vm(find(vm == 0)) = 1;
	in_g = nn_inputs(nn,off:end)./vm(off:end)';
min(in_g), max(in_g)

%	in_gs = (in_g + 0.88)* 8/(0.88+0.11) + -2.5;
%	plot(in_gs(1:end-5), 'Linewidth',LW);
%
%	plot([-1,-1],[8*0.88-3.5,8*(0.88+0.2)-3.5],'r','Linewidth',3);
%	text(-85,3,'.2mS/cm^2','FontSize',[16]);
%
	axis([-1,sim.T_upd-off,-3,8]);
	set(gca,'Visible','off');
	i=i+2;
	end;

subplot(3,1,3);
hold on;

plot([0,sim.T_upd-off],[-2.5,-2.5],'Linewidth',LW);
for i=0:50:sim.T_upd-off,
	plot([i,i],[-2.6,-2.1],'Linewidth',2);
	end;
plot([sim.T_upd-50-off, sim.T_upd],[-2.9,-2.9],'Linewidth',2);
text(sim.T_upd-50-off,-5.3,'50ms','FontSize',[16]);

%text1=sprintf('mu(%s) = [%.2f %.2f %.2f]', sim.exp, ...
%		nn_pars(1), nn_pars(2), nn_pars(3));
%	
%text(0,-5.3,text1,'FontSize',[16]);

subplot(3,1,1);
hold on;
%plot([0,0],[-20,30],'Linewidth',3);
%text(-100,10,'50mV','FontSize',[16]);
%plot([-15,-15],[-20,30],'r', 'Linewidth',3);
%text(-10,10,'50mV','FontSize',[16]);
%text(5,45,'50mV','FontSize',[16]);
text(-50,45,'50mV','FontSize',[16]);
plot([0,0],[-10,40],'r', 'Linewidth',3);

subplot(3,1,3);
hold on;
%plot([0,0],[4,6],'Linewidth',3);
%text(5,5,'2\muA/cm^2','FontSize',[16]);
plot([-1,-1],[-2,-0],'r','Linewidth',3);
%text(5,-1,'2\muA/cm^2','FontSize',[16]);
text(-80,-1,'2\muA/cm^2','FontSize',[16]);

%--------------------------------------------
%axis([1,sim.T_upd-off,0,sim.N_nn+1]);
%set(gca,'Visible','off');
xx=reshape(sim.instrument.allconduct,3,sim.T_upd);

subplot(3,1,2);
for i=1:3,
	fi_conduct=filter(ones(1,5)/5,1,xx(i,:));
%size(fi_conduct)
%	semilogy(abs(fi_conduct(off:end)),colors(i));
	plot(abs(fi_conduct(off:end-10)),colors(i),'LineWidth',LWs(i));
	hold on;

%	xxx=abs(fi_conduct(off:end));
%	mean(xxx(find(xxx <0.04)))
%	plot(sim.T_upd-off,mean(xxx(find(xxx <0.04))),'r*');
end;

%axis([1,sim.T_upd-off,1e-5,1e2]);
axis([1,sim.T_upd-off,0,0.04]);
%ylabel('Conductance [mS/cm^2]');
ylabel('mS/cm^2','FontSize',[16]);
set(gca,'FontSize',[16]);
%set(gca,'Visible','off');
grid;


subplot(3,1,3);
ti=1:sim.T_upd-off;
y_off_spikes = 2;
y_off_spikes = 5;
y_off_spikes = 6;
for sel=1:3,
        sp=find(sim.instrument.allvm(1,sel,off:end)> sim.activity_thr);
        if (length(sp) > 0),
                spp=zeros(1,length(sp))+y_off_spikes+0.6*sel;
                plot(ti(sp),spp,'r.','MarkerSize',8);
        end;
end;

%--------------------------------------------

fn_eps =sprintf('%s.eps', FN);
print('-deps', fn_eps);
fn_jpg =sprintf('%s.jpg', FN);
print('-djpeg', fn_jpg);

Loading data, please wait...