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;
/*===========================================
BUILD THE AXON ACCORDING TO SPECIFIED INPUTS
=============================================*/

proc buildAxon() {
	fiberDiameter=$1
	xElec=$2
	yElec=$3  
	zElec=$4
	rho=$5			
	k=(rho)/(4*PI)

	xyzPos=new Matrix(totalSections+1, 3)
	xyzPos.zero()
	xyPos()
	setPosMatrix()

	forall disconnect()
	createaxon()
	
	forall nseg = spatial_factor_sections

	distVec = new Vector(totalSections+1,0)
	xtrares = new Vector(totalSections+1,0)
	computeResistance()
} 

buildAxon(deffiberDiameter,defxElec,defyElec,defzElec,defrho)

access node[ (nodeSections-1)/2 ]

Loading data, please wait...