Multi-scale model of axonal and dendritic polarization by transcranial direct current stimulation in realistic head geometry

1.BackgroundTranscranial direct current stimulation (tDCS) is a non-invasive brain stimulation modality that can alter cortical excitability. However, it remains unclear how the subcellular elements of different neuron types are polarized by specific electric field (E-field) distributions. ObjectiveTo quantify neuronal polarization generated by tDCS in a multi-scale computational model. MethodsWe embedded layer-specific, morphologically-realistic cortical neuron models in a finite element model of the E-field in a human head and simulated steady-state polarization generated by conventional primary-motor-cortex-supraorbital (M1-SO) and 4x1 high-definition (HD) tDCS. We quantified somatic, axonal, and dendritic polarization of excitatory pyramidal cells in layers 2/3, 5, and 6, as well as inhibitory interneurons in layers 1 and 4 of the hand knob. ResultsAxonal and dendritic terminals were polarized more than the soma in all neurons, with peak axonal and dendritic polarization of 0.92 mV and 0.21 mV, respectively, compared to peak somatic polarization of 0.07 mV for 1.8 mA M1-SO stimulation. Both montages generated regions of depolarization and hyperpolarization beneath the M1 anode; M1-SO produced slightly stronger, more diffuse polarization peaking in the central sulcus, while 4x1 HD produced higher peak polarization in the gyral crown. Simulating polarization by uniform local E-field approximated the spatial distribution of tDCS polarization but produced large errors in some regions. ConclusionsPolarization of pre- and postsynaptic compartments of excitatory and inhibitory cortical neurons may play a significant role in tDCS neuromodulation. These effects cannot be predicted from the E-field distribution alone but rather require calculation of the neuronal response.