Striatal GABAergic microcircuit, dopamine-modulated cell assemblies (Humphries et al. 2009)

 Download zip file 
Help downloading and running models
To begin identifying potential dynamically-defined computational elements within the striatum, we constructed a new three-dimensional model of the striatal microcircuit's connectivity, and instantiated this with our dopamine-modulated neuron models of the MSNs and FSIs. A new model of gap junctions between the FSIs was introduced and tuned to experimental data. We introduced a novel multiple spike-train analysis method, and apply this to the outputs of the model to find groups of synchronised neurons at multiple time-scales. We found that, with realistic in vivo background input, small assemblies of synchronised MSNs spontaneously appeared, consistent with experimental observations, and that the number of assemblies and the time-scale of synchronisation was strongly dependent on the simulated concentration of dopamine. We also showed that feed-forward inhibition from the FSIs counter-intuitively increases the firing rate of the MSNs.
1 . Humphries MD, Wood R, Gurney K (2009) Dopamine-modulated dynamic cell assemblies generated by the GABAergic striatal microcircuit. Neural Netw 22:1174-88 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Neostriatum fast spiking interneuron;
Gap Junctions: Gap junctions;
Receptor(s): D1; D2; GabaA; AMPA; NMDA; Dopaminergic Receptor;
Transmitter(s): Dopamine; Gaba; Glutamate;
Simulation Environment: MATLAB;
Model Concept(s): Activity Patterns; Temporal Pattern Generation; Synchronization; Spatio-temporal Activity Patterns; Parkinson's; Action Selection/Decision Making; Connectivity matrix;
Implementer(s): Humphries, Mark D [m.d.humphries at]; Wood, Ric [ric.wood at];
Search NeuronDB for information about:  D1; D2; GabaA; AMPA; NMDA; Dopaminergic Receptor; Dopamine; Gaba; Glutamate;
function out = Experiment_RandomInput(DA, fname)
if nargin == 0
    DA = 0;
    fname = {'RandomInput'};

% set the model parameters
SIMPARAMS = StriatumNetworkParameters;

% name for the log file
SIMPARAMS.sim.logfname = [char(fname) '.log'];

% -------------------------------------------------------------------------
% set the DA level
SIMPARAMS.physiology.DA = DA; 

% set all the GABA the weights to 0 = ones(length(,1) .* 6.1; = ones(length(,1) .* 6.1; = ones(length(,1) .* 4.36; = ones(length(,1) .* (4.36 * 5); = ones(length(,1) .* (4.36 * 5); = ones(length(,1).* (150/5); 

% set the input parameters
SIMPARAMS.sim.tfinal = 10000; % length of simulation in msec
SIMPARAMS.input.CTX.r_MSSEG = ones(,1) .* 1.9; % Hz
SIMPARAMS.input.CTX.N_MSSEG = int32(ones(,1) .* 250);
SIMPARAMS.input.CTX.r_FSSEG = ones(,1) .* 1.9; % Hz
SIMPARAMS.input.CTX.N_FSSEG = int32(ones(,1) .* 250);
SIMPARAMS.input.CTX.alpha_FSSEG = ones(,1) .* 0.0;

save([char(fname) '_SIMPARAMS'],'SIMPARAMS');
% -------------------------------------------------------------------------
% Run the simulation
out = RunSimulation(SIMPARAMS);

% -------------------------------------------------------------------------
% Save the results to disc
save(char(fname), 'out');

% -------------------------------------------------------------------------
% plot the results
close all
figure(1); clf; plot(out.STms(:,2), out.STms(:,1), '.'); title('MSN raster plot')
figure(2); clf; plot(out.STfs(:,2), out.STfs(:,1), '.'); title('FSI raster plot')
figure(3); clf; plot(out.RecordChan_MS(:,1:25))
figure(4); clf; plot(out.RecordChan_MS(:,26:end))