Striatal NN model of MSNs and FSIs investigated effects of dopamine depletion (Damodaran et al 2015)

 Download zip file 
Help downloading and running models
Accession:169984
This study investigates the mechanisms that are affected in the striatal network after dopamine depletion and identifies potential therapeutic targets to restore normal activity.
Reference:
1 . Damodaran S, Cressman JR, Jedrzejewski-Szmek Z, Blackwell KT (2015) Desynchronization of Fast-Spiking Interneurons Reduces beta-Band Oscillations and Imbalance in Firing in the Dopamine-Depleted Striatum. J Neurosci 35:1149-59 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Neuron or other electrically excitable cell; Axon; Dendrite;
Brain Region(s)/Organism:
Cell Type(s): Neostriatum spiny direct pathway neuron; Neostriatum spiny indirect pathway neuron; Neostriatum fast spiking interneuron;
Channel(s): I Sodium; I Potassium; Kir;
Gap Junctions: Gap junctions;
Receptor(s): D1; D2; GabaA; Glutamate;
Gene(s):
Transmitter(s): Gaba; Glutamate;
Simulation Environment: GENESIS;
Model Concept(s): Synchronization; Detailed Neuronal Models; Parkinson's;
Implementer(s): Damodaran, Sriraman [dsriraman at gmail.com];
Search NeuronDB for information about:  Neostriatum spiny direct pathway neuron; Neostriatum spiny indirect pathway neuron; D1; D2; GabaA; Glutamate; I Sodium; I Potassium; Kir; Gaba; Glutamate;
/
DamodaranEtAl2015
MScell
channels
unusedChannels
AddCaSpines.g *
addchans.g *
addinput.g *
addoutput.g *
addoutput_old.g *
AddSynapticChannels.g *
connectCaChannels.g *
DA_files.txt *
globals.g
globalsCaComp.g *
include_channels.g *
MScell.p *
MScell_D1.g *
MScell_D1_old.g *
MScell_D2.g *
MScell08.p *
MScell08CM3_D1.p *
MScell08CM3_D2.p *
MScellSpine.g *
MScellSpineCtx.g *
MScellSpineTh.g *
MScellSyn_D1.g *
MScellSyn_D2.g *
parametersA_D1.g *
parametersA_D2.g *
parametersB.g *
parametersC.g *
parametersD.g *
proto.g *
spines.g *
SynParams_D1.g
SynParams_D2.g
SynParamsCtx.g *
SynParamsTh.g *
                            
//genesis

/***************************		MS Model, Version 5.10	**********************
**************************** 	      	globals.g 			**********************
Tom Sheehan tsheeha2@gmu.edu	thsheeha@vt.edu	703-538-8361
*******************************************************************************
	The capatilized parameters defined below are global and visable to all files
	Capatilized parameters should be treated as constants
******************************************************************************/
/* AB: Comments indicate the qfactor used with the various channels
*  conductance of CaN and CaL12 sum the inactivating and non-inactivating channel
*  conductances from the previous model
*/
        float ELEAK = -0.070
        float PI = 3.1415926
        float RA = 1.0;
        float RM = 8.69565217;
        float CM = 0.01;
        float EREST_ACT = -0.085
 	float TEMPERATURE = 35
	str  CA_BUFF_1 = "Ca_difshell_1"     // L and T type channels
	str  CA_BUFF_2 = "Ca_difshell_2"     // coupled to the other channels
	str  CA_BUFF_3 = "Ca_difshell_3"     // all calcium channels
	
	int CaDyeFlag = 0    // flags of calcium dye. "0" means NO calcium dyes.
                     // flag =2 : Fluo-4
                     // flag =3 : Fluo-5F
	int shellMode = 1     // we  have two shell-modes:
                     //  mode = 0 : detailed multi-shell model, using "difshell" object
                     //  mode = 1 : simple calcium pool adopted from  Sabatini's work(Sabatini, 2001, 2004)

//parameters determined by hand tuning to match spike width, AHP shape &amp, fI curve
//spike width with these globals plus spines = 0.88 ms
        str gNaFprox_D1={85500}  //qfactor = 1.2   // 95%
        str gNaFmid_D1={2594}    //95%
        str gNaFdist_D1={927}    //95%

        str gNaFprox_D2={99000}  //qfactor = 1.2 //110% 
        str gNaFmid_D2={3003}    //110%
        str gNaFdist_D2={1073}    //110%

       	str gKAfprox_D1={3214}   //qfactor=1.5 for inact -3214
        str gKAfmid_D1={571}     //1/qfactor=1.5 for act!!!   471  571- pretty good 591-better 571
	str gKAfdist_D1={554}    // 314 414 614 554
        str gKAfdist2_D1={314}    // 314

	str gKAfprox_D2={3214}   //qfactor=1.5 for inact 3214
        str gKAfmid_D2={471}     //1/qfactor=1.5 for act!!!  471  371
        str gKAfdist_D2={414}    //314 234 314
	str gKAfdist2_D2={314}    //314 234 314

        str gKAsprox_D1={277}    //qfactor=2	 
        str gKAsdist_D1={22.9}

        str gKAsprox_D2={304.7}    //qfactor=2  //110%	 
        str gKAsdist_D2={25.19}   //110%

        str gKIR_D1=5.25          //qfactor = 0.5 //125%
        str gKDR_D1={7.25}       //qfactor = 0.5  

        str gKIR_D2=4.2          //qfactor = 0.5
        str gKDR_D2={7.25}       //qfactor = 0.5  

	float gCaL13_D1 = 1.0625e-7  //qfactor=2
        float gCaL13_D2 = 0.796875e-7  //qfactor=2  //75%
	float gCaT  =  0.5875e-8
	float gCaR  =  6.5e-7
//	float gCaQ  =  1.5e-7
	float gCaN_D1 =   0.8e-7       //qfactor=2 //20%
	float gCaL12_D1 = 1.675e-7    //qfactor=2  //200%
        float gCaN_D2 =   2.5e-7       //qfactor=2
	float gCaL12_D2 = 0.8375e-7    //qfactor=2



Loading data, please wait...