Spiking neuron model of the basal ganglia (Humphries et al 2006)

 Download zip file 
Help downloading and running models
Accession:83559
A spiking neuron model of the basal ganglia (BG) circuit (striatum, STN, GP, SNr). Includes: parallel anatomical channels; tonic dopamine; dopamine receptors in striatum, STN, and GP; burst-firing in STN; GABAa, AMPA, and NMDA currents; effects of synaptic location. Model demonstrates selection and switching of input signals. Replicates experimental data on changes in slow-wave (<1 Hz) and gamma-band oscillations within BG nuclei following lesions and pharmacological manipulations.
Reference:
1 . Humphries MD, Stewart RD, Gurney KN (2006) A physiologically plausible model of action selection and oscillatory activity in the basal ganglia. J Neurosci 26:12921-42 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Basal ganglia;
Cell Type(s): Neostriatum medium spiny direct pathway GABA cell; Subthalamus nucleus projection neuron; Globus pallidus neuron; Abstract integrate-and-fire leaky neuron;
Channel(s):
Gap Junctions:
Receptor(s): Dopaminergic Receptor;
Gene(s):
Transmitter(s): Dopamine; Gaba; Glutamate;
Simulation Environment: MATLAB;
Model Concept(s): Oscillations; Parkinson's; Action Selection/Decision Making; Sleep; Rebound firing;
Implementer(s): Humphries, Mark D [m.d.humphries at shef.ac.uk];
Search NeuronDB for information about:  Neostriatum medium spiny direct pathway GABA cell; Dopaminergic Receptor; Dopamine; Gaba; Glutamate;
/
BG_spiking_model
LIFtools
Analysis
find_cur.m
find_Vm_cur.m
plot_fI.m
PSPtoPSC.m
                            
function I = find_cur(tm,R,theta,abs_ref,fire,varargin)

% FIND_CUR determines the current injection needed for a given firing rate for basic LIF neuron
%
%   find_cur(tm,R,theta,abs_ref,fire)
%       tm:         membrane time constant      (in seconds)
%       R:          resistance                  (in ohms)
%       theta:      threshold of firing         (in volts)          [value above resting potential (itself assumed to be 0)]
%       abs_ref:    absolute refractory period  (in seconds)
%       fire:       target firing rate          (in Hz)
%
%   find_cur(tm,R,theta,abs_ref,fire,'d') 
%       the 'd' switch displays the required current in "physiological" units to screen
%
%   Returns the value of the current (in amps) required to drive the LIF neuron to achieve the specified
%   firing rate. If the specified firing rate cannot be achieved with the given parameters, then
%   returns NaN and a warning given.
%
%   Mark Humphries. Last rev: 17/12/2004

I = -theta / (R * exp(-(1/tm * (1/fire - abs_ref))) - R);

if I < 0
    I = nan;
end

if nargin == 6
    if isnan(I)
        disp(['Cannot achieve firing rate of ' num2str(fire) 'Hz with the specified LIF neuron']);
    else
        disp(['A current of ~' num2str(I * 1e9) 'nA is required for an output of ' num2str(fire) 'Hz (for specified LIF neuron)']); 
    end
end

Loading data, please wait...