Specific inhibition of dendritic plateau potential in striatal projection neurons (Du et al 2017)

 Download zip file 
Help downloading and running models
We explored dendritic plateau potentials in a biophysically detailed SPN model. We coupled the dendritic plateaus to different types of inhibitions (dendritic fast and slow inhibitions, perisomatic inhibition from FS interneurons , etc.) We found the inhibition provides precise control over the plateau potential, and thus the spiking output of SPNs.
1 . Du K, Wu YW, Lindroos R, Liu Y, Rózsa B, Katona G, Ding JB, Kotaleski JH (2017) Cell-type-specific inhibition of the dendritic plateau potential in striatal spiny projection neurons. Proc Natl Acad Sci U S A 114:E7612-E7621 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Dendrite; Channel/Receptor;
Brain Region(s)/Organism: Striatum;
Cell Type(s): Neostriatum spiny neuron;
Channel(s): I A; I A, slow; Kir; I Calcium; I K,Ca; I L high threshold; I Na,p; I Na,t; I Q; I R; I K;
Gap Junctions:
Receptor(s): AMPA; GabaA; NMDA;
Gene(s): Cav3.2 CACNA1H;
Transmitter(s): Glutamate; Gaba;
Simulation Environment: GENESIS;
Model Concept(s):
Implementer(s): Du, Kai [kai.du at ki.se];
Search NeuronDB for information about:  GabaA; AMPA; NMDA; I Na,p; I Na,t; I L high threshold; I A; I K; I K,Ca; I Calcium; I A, slow; I R; I Q; Kir; Gaba; Glutamate;

include MScell/globals_1		  		// Defines & initializes cell specific parameters
include clamps.g

/***************************		MS Model, Version 5.17	***************
****************************** Begin includes *******************************/
int usingHsolve = 1
int usingNewKaF = 0
int usingNaP    = 1
int addCa2Spine = 1
int useAMPA = 1
int use_F_factor = 1
int usingElectro_Leak = 0      // to simulate eletro-leak in the soma 
int usingClamp = 0
int usingCaT = 1
int usingCaR = 1
int usingKIR = 1
int usingKAs = 1
int usingKDR = 1
int showMSoutput = 1

str 	neuronname = "/cell"

//str compt_list2={findCompt 90e-6 110e-6 5e-6 {neuronname} } 

/* F is the factor for spine compensation
*** you can ignore this part*******************/
if ({{use_F_factor}==1})
float F = 1.38  // 1.38
float F = 1
float F2 = 1     // RA in the distal dendrite = 4, Surmeier's paper 

include MScell/MScellSpine_test                // access make_MS_cell this function is only called from MSsim.g
str pfile="MScell/swc5_SPN.g"

include xcell		        // access functions "record_channel"
				// These two functions are only called from MSsim.g

/****************************** End includes *********************************/
float speedA = 20
	setclock 0 {speedA*1e-6}   // Simulation time step (Second)       
	setclock 1 {speedA*1e-6}   // clock for graphic output

/******************* Compensate for spines***********************************/
make_MS_cell {neuronname} {pfile} 30.0e-6 160.0e-6 {F} {F2} 	// MS_cell.g 
//echo The current cell is {neuronname}

/*****************************Visualize the neuron**********************************/
str list = {findCompt {neuronname} 1e-6 200e-6 0.5e-6 }
increaseXcellDia {neuronname} {list} 2
  create xform /cellform1 [620,50,1000,1000]
    create xdraw /cellform1/draw [0,0,100%,100%] \
                             -wx 3e-3  \
                             -wy 3e-3  \
                             -transform ortho3d \
                             -bg black
    //setfield {neuronname}/soma dia 8e-6
    setfield /cellform1/draw xmin -1.2e-4 xmax 1.2e-4 ymin -1.2e-4 ymax 1.2e-4 \
        zmin -1e-3 zmax 1e-3 

     setfield /cellform1/draw transform z

    xshow /cellform1
    echo creating xcell
    create xcell /cellform1/draw/cell
    setfield /cellform1/draw/cell colmin -0.09 colmax -0.01 \
                                  path {neuronname}/##[TYPE=compartment]   field Vm \
                                  script "echo widget clicked on = <w> value = <v>" \
                                  diarange -50

     xcolorscale hot

    str above = "cell" 
    str parent  =  "/cellform1/draw"
  //make_colorbar {parent} {above}

str compt_1 = "1409"     // distal: 1409; proximal: 1146
add_spines {neuronname} {compt_1} spine 15         // located 90-100 um from the soma


// insert different types of inhibitions at the locations of interest 
// "MSN_GABA_channel" : dendritic fast inhibition 
//"NPY_NGF_GABA_channel": dendritic slow inhibition 
//"FS_GABAA"            : perisomatic inhibitions from FS interneurons
// Note: everytime you can only test one type of inhibition
addGABAchannel   {neuronname} "1409 1306 1245 1574 1437 soma"  "MSN_GABA_channel" "MSN_GABAA" 1500e-12   2
//addGABAchannel   {neuronname} "1409 1306 1245 1574 1437 1351"  "NPY_NGF_GABA_channel" "NPY_NGF_GABAA" 1500e-12   2
//addGABAchannel   {neuronname} "soma"  "GABA_channel" "FS_GABAA" 1500e-12   2

if ({usingHsolve == 1})
pushe {neuronname}
  create hsolve solve
  setfield  solve  \
           chanmode 1 \
           path     {neuronname}/##[][TYPE=compartment]

  call ^ SETUP

 setmethod 11


/****************************** End MSsim.g **********************************/

// // /********************************************************
// // ************* begin to draw results**********************
// // *********************************************************/

/*****function make_xcells(xcell, label, title1,title2,title3,title4,tmin,tmax)********/
str xcell ="xcell"
str label = "Dendritic Plateau at "@{compt_1}
str title1 = "Membrane_Voltage"
str xcellPath = "output/"@{xcell}
int chanmode = 0  
float tmin = 0.2
float tmax = 0.5

make_VerticalXcells {xcellPath} {label} {title1} {tmin} {tmax} 
/*********function record_channel(cellpath,compt,channel,xcell,title, color,ymin,ymax)********/
/*********function record_voltage (cellpath, compt, xcell,lable,color,ymin,ymax,chanmode)****/
/******************** record voltage ****************************************************/
float ymin_volt = -0.09
float ymax_volt = -0.0
record_voltage   {neuronname}   soma                 {xcellPath}  {title1} "soma"        "blue"    {ymin_volt} {ymax_volt} {chanmode}    
record_voltage {neuronname}     {compt_1}         {xcellPath}  {title1}    "dendrite"     "red"     {ymin_volt} {ymax_volt} {chanmode}   // to record voltage of the compartment            
//record_voltage {neuronname}    1437               {xcellPath} {title1}     "distal tip"   "black"  {ymin_volt} {ymax_volt} {chanmode}
//record_voltage {neuronname}    746                {xcellPath} {title1}     "another branch"   "pink" {ymin_volt} {ymax_volt} {chanmode}



str subpath
int n, n0,n_max,num_spines,i,j,jj,num_gaba
float dt,t1,t2,isi,inject,dur,delay,dt_2

num_spines = 15 
isi = 0.001

// Simulation protocols 

// we now test inhibitions at different locations. "1409" is compartment name for the "on-site" inhibition
// if you want to test FS-inhibition, you need to change those compartments in the arglist to "soma"  

foreach subpath ({arglist "1409 "})

setfield /{neuronname}/soma inject 100e-12
step 0.3 -time
activate_spines {neuronname} {compt_1} "spine" 15 {isi} 1 
step 0.005 -time                                                // delay between the plateau potential and the inhibition
activateGABAchannels {neuronname} {subpath} 2 0.0
step 0.3 -time

//step 1