Investigation of different targets in deep brain stimulation for Parkinson`s (Pirini et al. 2009)

 Download zip file 
Help downloading and running models
"We investigated by a computational model of the basal ganglia the different network effects of deep brain stimulation (DBS) for Parkinson’s disease (PD) in different target sites in the subthalamic nucleus (STN), the globus pallidus pars interna (GPi), and the globus pallidus pars externa (GPe). A cellular-based model of the basal ganglia system (BGS), based on the model proposed by Rubin and Terman (J Comput Neurosci 16:211–235, 2004), was developed. ... Our results suggest that DBS in the STN could functionally restore the TC relay activity, while DBS in the GPe and in the GPi could functionally over-activate and inhibit it, respectively. Our results are consistent with the experimental and the clinical evidences on the network effects of DBS."
1 . Pirini M, Rocchi L, Sensi M, Chiari L (2009) A computational modelling approach to investigate different targets in deep brain stimulation for Parkinson's disease. J Comput Neurosci 26:91-107 [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: Neocortex; Thalamus; Basal ganglia; Subthalamic Nucleus;
Cell Type(s): Thalamus geniculate nucleus/lateral principal neuron; Subthalamus nucleus projection neuron; Globus pallidus neuron;
Channel(s): I Na,t; I T low threshold; I K; I Calcium;
Gap Junctions:
Simulation Environment: MATLAB;
Model Concept(s): Activity Patterns; Parkinson's; Deep brain stimulation;
Implementer(s): Pirini, Marco [marco.pirini at];
Search NeuronDB for information about:  Thalamus geniculate nucleus/lateral principal neuron; I Na,t; I T low threshold; I K; I Calcium;
Author: Marco Pirini
July 7th, 2009
(Please see Readme.doc for a nicer formatted version of this text file)


-	the model has been created in a Matlab-Simulink environment;
-       the model has been divided in 3, functionally non re-entrant,
        modules: the STNGPE one, the GPi one and the Thalamic one
        (possibility to re-use some simulation results from the
        upstream modules without the need to re-simulate the entire
        model ? reduced computation time & more control).
-	The flow to simulate is:

1 launch launch_sysGPESTN1.m

o       upload IC, parameters, additional currents, DBS currents if
        needed, etc.
o       launch the STNGPE module (GPESTN1.mdl - in this .mdl you will
        see the architecture used by Rubin and Terman as arised from a
        private communication with them) o save results in a .mat file
        for the STNGPE simulation;


1 launch launch_sysGPESTN1_GPEDBS.m to simulate DBS on GPe target if
        needed. The steps performed by this script are the same as by
        launch_sysGPESTN1.m. Remember: when asked, choose the
2 launch launch_sysGPI_mod.m
o       upload IC, parameters, additional currents, DBS currents if
        needed, etc
o	launch the  GPI module (GPI1.mdl)
o	save results in a .mat file for the GPI simulation;

3 launch launch_sysTAL_mod.m

o       upload IC, parameters, additional currents, the cortical
        input, etc
o	launch the  TAL module (TAL8.mdl)
o       analyze the results (thalamic spike sorting and
o	save results in a .mat file for the TAL simulation;


o       once simulated all the possible (of interest) sets of
        parameters and conditions for the SNTGPe and the GPi modules
        (different conditions (norm, park, DBS on the 3 targets),
        different values of the ggpe-->gpi conductance, different
        values of the Istriatum-->Gpi, different DBS frequencies...)

it`s possible to perform in a single step the simulation of the final
thalamic module for all of the set of parameters (and for several
realisations of the cortical input sequence to thalamus). You could do
this by launching CN_launch_CN_analyze_TAL2.m that , in turn, recalls

-	Other files

o       data_STN.m, data_GPe.m, data_GPi.m, data_TAL.m,
        data_synapses.m: sets of parameters for populations and
o	in_val.mat: some sets of Initial conditions;
o	sp_rev_thresh.m: thalamic spikes identification
o	spike_contr.m: thalamic spike classification
o       CN_calculate_sm.m: scripts that calculate a Poisson
        realization of the cortical input sequence to thalamus.