Mechanisms of magnetic stimulation of central nervous system neurons (Pashut et al. 2011)

 Download zip file   Auto-launch 
Help downloading and running models
Transcranial magnetic stimulation (TMS) is a widely applied tool for probing cognitive function in humans and is one of the best tools for clinical treatments and interfering with cognitive tasks. Surprisingly, while TMS has been commercially available for decades, the cellular mechanisms underlying magnetic stimulation remain unclear. Here we investigate these mechanisms using compartmental modeling. We generated a numerical scheme allowing simulation of the physiological response to magnetic stimulation of neurons with arbitrary morphologies and active properties. Computational experiments using this scheme suggested that TMS affects neurons in the central nervous system (CNS) primarily by somatic stimulation.
1 . Pashut T, Wolfus S, Friedman A, Lavidor M, Bar-Gad I, Yeshurun Y, Korngreen A (2011) Mechanisms of magnetic stimulation of central nervous system neurons. PLoS Comput Biol 7:e1002022 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Neocortex L5/6 pyramidal GLU cell; Squid axon;
Gap Junctions:
Simulation Environment: NEURON; MATLAB;
Model Concept(s): Action Potential Initiation; Magnetic stimulation;
Implementer(s): Korngreen, Alon [alon.korngreen at]; Pashut, Tamar [tamar.pashut at];
Search NeuronDB for information about:  Neocortex L5/6 pyramidal GLU cell;
 create  dend[2]
 objref sec_list
  access dend
  proc celldef() {
  access dend[0] {distance(0, 0)}

proc basic_shape() {
   connect dend[0](0), dend[1](0)
  dend[0] {
  	pt3dadd(0, 0, 0, 1) 
  	pt3dadd(-20000, 0, 0, 1)
  	pt3dadd(-40000, 0, 0, 1)
  	pt3dadd(-60000, 0, 0, 1)
  	pt3dadd(-80000,0 , 0, 1)
  dend[1] {
  	pt3dadd(0, 0, 0, 1) 
  	pt3dadd(20000, 0, 0, 1)
  	pt3dadd(40000, 0, 0, 1)
  	pt3dadd(60000, 0, 0, 1)
  	pt3dadd(79999, 0, 0, 1)
  forall { L =arc3d(n3d()-1)   diam = 100  }  
  forall { nseg = 100  }
  sec_list = new SectionList()
  for i=0, 1 dend[i] sec_list.append()

proc biophys() {
  forall {
    Ra = 35.4
    cm = 1
    insert hh
    gnabar_hh = 0.12
    gkbar_hh = 0.036
    gl_hh = 0.0003
    el_hh = -54.387
    insert xtra