Submyelin Potassium accumulation in Subthalamic neuron (STN) axons (Bellinger et al. 2008)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:121253
"To better understand the direct effects of DBS (Deep brain stimulation) on central neurons, a computational model of a myelinated axon has been constructed which includes the effects of K+ accumulation within the peri-axonal space. Using best estimates of anatomic and electrogenic model parameters for in vivo STN axons, the model predicts a functional block along the axon due to K+ accumulation in the submyelin space. ... These results suggest that therapeutic DBS of the STN likely results in a functional block for many STN axons, although a subset of STN axons may also be activated at the stimulating frequency. "
Reference:
1 . Bellinger SC, Miyazawa G, Steinmetz PN (2008) Submyelin potassium accumulation may functionally block subsets of local axons during deep brain stimulation: a modeling study. J Neural Eng 5:263-74 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Axon;
Brain Region(s)/Organism:
Cell Type(s): Subthalamus nucleus projection neuron;
Channel(s): I Na,p; I K; I Sodium; I_Ks; Na/K pump;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Axonal Action Potentials; Action Potentials; Deep brain stimulation; Sodium pump; Depolarization block;
Implementer(s): Bellinger, Steven [Steve.Bellinger at asu.edu];
Search NeuronDB for information about:  I Na,p; I K; I Sodium; I_Ks; Na/K pump;
objref xtraClamp, tstim, tpstim
xtraClamp=new trainIClamp()
create electrode

	DELAY=0
	PW=0.1
	DURATION=10000
	FREQ=130
	AMPLITUDE=-10000 //nA

proc stimulus() {
	electrode {
	  xtraClamp.loc(0.5)
	  xtraClamp.del=$1
          xtraClamp.PW=$2
          xtraClamp.train=$3
          xtraClamp.freq=$4
          xtraClamp.amp=$5
          }

	//Record stimulation times
	tstim=new Vector(DURATION/(10^3/FREQ))
	for(ej=0; ej<tstim.size(); ej=ej+1){
		tstim.x[ej]=ej*1/FREQ*10^3
		}
	
	//Record stimulation times + psuedo stimulation times
	tpstim=new Vector(tstop/(10^3/FREQ))
	for(ej=0; ej<tpstim.size(); ej=ej+1){
		tpstim.x[ej]=ej*1/FREQ*10^3
		}
	}

stimulus(DELAY,PW,DURATION,FREQ,AMPLITUDE)

Loading data, please wait...