Synchronicity of fast-spiking interneurons balances medium-spiny neurons (Damodaran et al. 2014)

 Download zip file 
Help downloading and running models
This study investigates the role of feedforward and feedback inhibition in maintaining the balance between D1 and D2 MSNs of the striatum. The synchronized firing of FSIs are found to be critical in this mechanism and specifically the gap junction connections between FSIs.
1 . Damodaran S, Evans RC, Blackwell KT (2014) Synchronized firing of fast-spiking interneurons is critical to maintain balanced firing between direct and indirect pathway neurons of the striatum. J Neurophysiol 111:836-48 [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): Neostriatum medium spiny direct pathway GABA cell; Neostriatum medium spiny indirect pathway GABA cell; Neostriatum fast spiking interneuron;
Gap Junctions: Gap junctions;
Receptor(s): NMDA; Gaba;
Simulation Environment: GENESIS;
Model Concept(s): Detailed Neuronal Models; Parkinson's;
Implementer(s): Blackwell, Avrama [avrama at]; Damodaran, Sriraman [dsriraman at];
Search NeuronDB for information about:  Neostriatum medium spiny direct pathway GABA cell; Neostriatum medium spiny indirect pathway GABA cell; NMDA; Gaba;
Inputwithcorrelation.asv *
Inputwithcorrelation.m *
InputwithCorrelation2.m *
makeDaughterInput.m *
makeDaughterInsignal_no_upstate.m *
makeTrainInput.m *
makeTrainInsignal.m *
poissonMaxTime.m *
writeInput.asv *
writeInput.m *
% Generates input for each cell with Nampa(nSyn)columns and rows proportional to freq*maxTime
% nSyn in the number of synapses
% corr_syn is the correlation between the synapses

function noise = makeDaughterInput(corr_syn, nSyn, freq, maxTime)

	nShare = nSyn - sqrt(corr_syn)*(nSyn-1);
	pShare = 1/nShare;
	motherSpikes = poissonMaxTime(freq*nShare, maxTime);
	len = length(motherSpikes); % # rows of final array

	v = (rand(len, nSyn) < pShare).*repmat(motherSpikes,1,nSyn); %0 indicates don't assign spikes and 1 assigns
	v(find(v == 0)) = inf;
	v = sort(v,1); % sorts each column in ascending order (2nd argument is ascending or descending)
	vlen = 1+max(mod(find(v < inf) - 1, len));
	v(find(v == inf)) = 0;
	noise = v(1:vlen,:);