Computational Modeling of Ultrasonic Subthalamic Nucleus Stimulation.

To explore the potential of ultrasonic modulation of plateau-potential generating subthalamic nucleus neurons (STN), by modeling their interaction with continuous and pulsed ultrasonic waves. A computational model for ultrasonic stimulation of the STN is created by combining the Otsuka-model with the bilayer sonophore model. The neuronal response to continuous and pulsed ultrasonic waves is computed in parallel for a range of frequencies, duty cycles, pulse repetition frequencies, and intensities. Ultrasonic intensity in continuous-wave stimulation determines the firing pattern of the STN. Three observed spiking modes in order of increasing intensity are low frequency spiking, high frequency spiking with significant spike-frequency and spike-amplitude adaptation, and a silenced mode. Continuous-wave stimulation has little capability to manipulate the saturated spiking rate in the high frequency spiking mode. In contrast, STN firing rates induced by pulsed ultrasound insonication will saturate to the pulse repetition frequency with short latencies, for sufficiently large intensity and repetition frequency. Computational results show that the activity of plateau-potential generating STN can be modulated by selection of the stimulus parameters. Low intensities result in repetitive firing, while higher intensities silence the STN. Pulsed ultrasonic stimulation results in a shorter saturation latency and is able to modulate spiking rates. Stimulation or suppresion of the STN is important in the treatment of Parkinson's disease, e.g., in deep brain stimulation. This explorative study on ultrasonic modulation of the STN, could be a step in the direction of minimally invasive alternatives to conventional deep brain stimulation.