Synchronicity of fast-spiking interneurons balances medium-spiny neurons (Damodaran et al. 2014)

 Download zip file 
Help downloading and running models
Accession:156260
This study investigates the role of feedforward and feedback inhibition in maintaining the balance between D1 and D2 MSNs of the striatum. The synchronized firing of FSIs are found to be critical in this mechanism and specifically the gap junction connections between FSIs.
Reference:
1 . Damodaran S, Evans RC, Blackwell KT (2014) Synchronized firing of fast-spiking interneurons is critical to maintain balanced firing between direct and indirect pathway neurons of the striatum. J Neurophysiol 111:836-48 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism:
Cell Type(s): Neostriatum medium spiny direct pathway GABA cell; Neostriatum medium spiny indirect pathway GABA cell; Neostriatum fast spiking interneuron;
Channel(s):
Gap Junctions: Gap junctions;
Receptor(s): NMDA; Gaba;
Gene(s):
Transmitter(s):
Simulation Environment: GENESIS;
Model Concept(s): Detailed Neuronal Models; Parkinson's;
Implementer(s): Blackwell, Avrama [avrama at gmu.edu]; Damodaran, Sriraman [dsriraman at gmail.com];
Search NeuronDB for information about:  Neostriatum medium spiny direct pathway GABA cell; Neostriatum medium spiny indirect pathway GABA cell; NMDA; Gaba;
/
striatalnetwork
MScell
channels
unusedChannels
AddCaSpines.g *
addchans.g *
addchans.g~ *
addinput.g *
addoutput.g *
addoutput.g~
addoutput_old.g *
AddSynapticChannels.g *
connectCaChannels.g *
DA_files.txt *
globals.g
globals.g~
globals_old.g
globalsCaComp.g *
include_channels.g *
MScell.g~
MScell.p *
MScell_D1.g *
MScell_D1.g~
MScell_D1_old.g
MScell_D2.g *
MScell_D2.g~
MScell_D2_old.g
MScell08.p *
MScell08CM3_D1.p *
MScell08CM3_D2.p *
MScell08CM3_old.p
MScellSpine.g *
MScellSpineCtx.g *
MScellSpineTh.g *
MScellSyn.g~ *
MScellSyn_D1.g *
MScellSyn_D2.g *
parametersA.g~
parametersA_D1.g *
parametersA_D2.g *
parametersB.g *
parametersC.g *
parametersD.g *
proto.g *
proto.g~
spines.g *
SynParams.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={99000}  //qfactor = 1.2   // 95%
        str gNaFmid_D1={3003}    //95%
        str gNaFdist_D1={1073}    //95%

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

        str gKAfprox={3214}   //qfactor=1.5 for inact -3214
        str gKAfmid={471}     //1/qfactor=1.5 for act!!!
        str gKAfdist={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