ModelDB is moving. Check out our new site at The corresponding page is

Gap junction coupled network of striatal fast spiking interneurons (Hjorth et al. 2009)

 Download zip file 
Help downloading and running models
Gap junctions between striatal FS neurons has very weak ability to synchronise spiking. Input uncorrelated between neighbouring neurons is shunted, while correlated input is not.
1 . Hjorth J, Blackwell KT, Kotaleski JH (2009) Gap junctions between striatal fast-spiking interneurons regulate spiking activity and synchronization as a function of cortical activity. J Neurosci 29:5276-86 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Neuron or other electrically excitable cell; Synapse; Channel/Receptor; Dendrite;
Brain Region(s)/Organism: Basal ganglia;
Cell Type(s): Neostriatum fast spiking interneuron;
Channel(s): I A; I_K,Na;
Gap Junctions: Gap junctions;
Transmitter(s): Gaba; Glutamate;
Simulation Environment: GENESIS; MATLAB;
Model Concept(s): Activity Patterns; Ion Channel Kinetics; Synchronization; Detailed Neuronal Models;
Implementer(s): Hjorth, Johannes [hjorth at];
Search NeuronDB for information about:  I A; I_K,Na; Gaba; Glutamate;
% This version of the file only calculates the IF plots
% (It is a reduced version of the full script which uses current
% injection instead of synaptic inputs)
% To generate Figure 1A run the following MATLAB scripts:
% runIFscanTenFS.m
% readIFscanTenFS.m
% makeIFplots.m

clear all, format compact

% Matlab helper scripts are located here

% Genesis model is located here

nReps = 1 % 10 
nPoints = 20;

maxTime = 10;
numCells = 30;

channelMask = {'A_channel'} % Channel conductances we want to vary
cellVar = 0.5 % 0.2 % 0.1;
lenVar = 0.5 % 0.2 %0.1;

curIFrange = linspace(0,0.1e-9,nPoints);

dataFileCur = 'TenInhomoFS-curInject-IFscan'

dataDir = 'UTDATA/SAVED/TenFSIFscanHighRES/';

for rep=1:nReps

    disp('Pausing for 1 second, press Ctrl+C to abort')
    %% Generate input! First, seed...
    randSeed = floor(sum(clock)*1e5);

    % Generera FS morphologin
    makeFSMorph(numCells, cellVar, channelMask, lenVar)
    % Overwrite first cell morphology with the standard cell
    makeFSMorph(1, 0, channelMask, 0)

    maxCurTime = maxTime;
    % No gap junctions
    for j=1:length(curIFrange)
      for k=1:numCells
        curStart{k} = 0;
        curEnd{k} = maxCurTime + 1;
        curAmp{k} = curIFrange(j);
        curLoc{k} = sprintf('/fs[%d]/soma', k-1);
      % Write parameters to file
      writeCurrentInputInfo(curStart, curEnd, curAmp, curLoc)

      % This one ignores all synaptic input
      system('genesis ../genesisScripts/simFsMultiInhomogeneCurrentInjection');
      saveFileData = [dataDir dataFileCur ...
                      '-ID' num2str(randSeed) ...
                      '-cur-' num2str(curIFrange(j)) ...
      saveFileInfo = strrep(saveFileData,'.data','.info');

      system(['cp UTDATA/' dataFileCur '.data ' saveFileData]);
      system(['cp INDATA/parameters.txt ' saveFileInfo]);
      system(['cat INDATA/currentInputInfo.txt >> ' saveFileInfo]);

%% Add code to generate figures here


Loading data, please wait...