Model-based deconstruction of cortical evoked potentials generated by subthalamic nucleus deep brain stimulation.

Parkinson's disease is associated with altered neural activity in the motor cortex. Chronic high-frequency deep brain stimulation (DBS) of the subthalamic nucleus (STN) is effective in suppressing parkinsonian motor symptoms and modulates cortical activity. However, the anatomical pathways responsible for STN DBS-mediated cortical modulation remain unclear. Cortical evoked potentials (cEP) generated by STN DBS reflect the response of cortex to subcortical stimulation, and the goal of this study was to determine the neural origin of STN DBS-generated cEP using a two-step approach. First, we recorded cEP over ipsilateral primary motor cortex during different frequencies of STN DBS in awake healthy and unilateral 6-OHDA-lesioned parkinsonian rats. Second, we used a detailed, biophysically based model of the thalamocortical network to deconstruct the neural origin of the recorded cEP. The in vivo cEP included short (R1)-, intermediate (R2)-, and long-latency (R3) responses. Model-based cortical responses to simulated STN DBS matched remarkably well the in vivo responses. The short-latency response was generated by antidromic activation of layer 5 pyramidal neurons, whereas recurrent activation of layer 5 pyramidal neurons via excitatory axon collaterals reproduced the intermediate-latency response. The long-latency response was generated by polysynaptic activation of layer 2/3 pyramidal neurons via the cortico-thalamic-cortical pathway. Antidromic activation of the hyperdirect pathway and subsequent intracortical and cortico-thalamo-cortical synaptic interactions were sufficient to generate cortical potential evoked by STN DBS, and orthodromic activation through basal ganglia-thalamus-cortex pathways was not required. These results demonstrate the utility of cEP to determine the neural elements activated by STN DBS that might modulate cortical activity and contribute to the suppression of parkinsonian symptoms. NEW & NOTEWORTHY Subthalamic nucleus (STN) deep brain stimulation (DBS) is increasingly used to treat Parkinson's disease (PD). Cortical potentials evoked by STN DBS in patients with PD exhibit consistent short-latency (1-3 ms), intermediate-latency (5-15 ms), and long-latency (18-25 ms) responses. The short-latency response occurs as a result of antidromic activation of the hyperdirect pathway comprising corticosubthalamic axons. However, the neural origins of intermediate- and long-latency responses remain elusive, and the dominant view is that these are produced through the orthodromic pathway (basal ganglia-thalamus-cortex). By combining in vivo electrophysiology with computational modeling, we demonstrate that antidromic activation of the cortico-thalamic-cortical pathway is sufficient to generate the intermediate- and long-latency cortical responses to STN DBS.

A biophysical model of the cortex-basal ganglia-thalamus network in the 6-OHDA lesioned rat model of Parkinson's disease.

Electrical stimulation of sub-cortical brain regions (the basal ganglia), known as deep brain stimulation (DBS), is an effective treatment for Parkinson's disease (PD). Chronic high frequency (HF) DBS in the subthalamic nucleus (STN) or globus pallidus interna (GPi) reduces motor symptoms including bradykinesia and tremor in patients with PD, but the therapeutic mechanisms of DBS are not fully understood. We developed a biophysical network model comprising of the closed loop cortical-basal ganglia-thalamus circuit representing the healthy and parkinsonian rat brain. The network properties of the model were validated by comparing responses evoked in basal ganglia (BG) nuclei by cortical (CTX) stimulation to published experimental results. A key emergent property of the model was generation of low-frequency network oscillations. Consistent with their putative pathological role, low-frequency oscillations in model BG neurons were exaggerated in the parkinsonian state compared to the healthy condition. We used the model to quantify the effectiveness of STN DBS at different frequencies in suppressing low-frequency oscillatory activity in GPi. Frequencies less than 40 Hz were ineffective, low-frequency oscillatory power decreased gradually for frequencies between 50 Hz and 130 Hz, and saturated at frequencies higher than 150 Hz. HF STN DBS suppressed pathological oscillations in GPe/GPi both by exciting and inhibiting the firing in GPe/GPi neurons, and the number of GPe/GPi neurons influenced was greater for HF stimulation than low-frequency stimulation. Similar to the frequency dependent suppression of pathological oscillations, STN DBS also normalized the abnormal GPi spiking activity evoked by CTX stimulation in a frequency dependent fashion with HF being the most effective. Therefore, therapeutic HF STN DBS effectively suppresses pathological activity by influencing the activity of a greater proportion of neurons in the output nucleus of the BG.

Neural mechanisms of the temporal response of cortical neurons to intracortical microstimulation.

BACKGROUND: Intracortical microstimulation (ICMS) is used to map neuronal circuitry in the brain and restore lost sensory function, including vision, hearing, and somatosensation. The temporal response of cortical neurons to single pulse ICMS is remarkably stereotyped and comprises short latency excitation followed by prolonged inhibition and, in some cases, rebound excitation. However, the neural origin of the different response components to ICMS are poorly understood, and the interactions between the three response components during trains of ICMS pulses remains unclear. OBJECTIVE: We used computational modeling to determine the mechanisms contributing to the temporal response to ICMS in model cortical neurons. METHODS: We implemented a biophysically based computational model of a cortical column comprising neurons with realistic morphology and synapses and quantified the temporal response of cortical neurons to different ICMS protocols. We characterized the temporal responses to single pulse ICMS across stimulation intensities and inhibitory (GABA-B/GABA-A) synaptic strengths. To probe interactions between response components, we quantified the response to paired pulse ICMS at different inter-pulse intervals and the response to short trains at different stimulation frequencies. Finally, we evaluated the performance of biomimetic ICMS trains in evoking sustained neural responses. RESULTS: Single pulse ICMS evoked short latency excitation followed by a period of inhibition, but model neurons did not exhibit post-inhibitory excitation. The strength of short latency excitation increased and the duration of inhibition increased with increased stimulation amplitude. Prolonged inhibition resulted from both after-hyperpolarization currents and GABA-B synaptic transmission. During the paired pulse protocol, the strength of short latency excitation evoked by a test pulse decreased marginally compared to those evoked by a single pulse for interpulse intervals (IPI) < 100 m s. Further, the duration of inhibition evoked by the test pulse was prolonged compared to single pulse for IPIs <50 m s and was not predicted by linear superposition of individual inhibitory responses. For IPIs>50 m s, the duration of inhibition evoked by the test pulse was comparable to those evoked by a single pulse. Short ICMS trains evoked repetitive excitatory responses against a background of inhibition. However, the strength of the repetitive excitatory response declined during ICMS at higher frequencies. Further, the duration of inhibition at the cessation of ICMS at higher frequencies was prolonged compared to the duration following a single pulse. Biomimetic pulse trains evoked comparable neural response between the onset and offset phases despite the presence of stimulation induced inhibition. CONCLUSIONS: The cortical column model replicated the short latency excitation and long-lasting inhibitory components of the stereotyped neural response documented in experimental studies of ICMS. Both cellular and synaptic mechanisms influenced the response components generated by ICMS. The non-linear interactions between response components resulted in dynamic ICMS-evoked neural activity and may play an important role in mediating the ICMS-induced precepts.

Stoney vs. Histed: Quantifying the spatial effects of intracortical microstimulation.

BACKGROUND: Intracortical microstimulation (ICMS) is used to map neural circuits and restore lost sensory modalities such as vision, hearing, and somatosensation. The spatial effects of ICMS remain controversial: Stoney and colleagues proposed that the volume of somatic activation increased with stimulation intensity, while Histed et al., suggested activation density, but not somatic activation volume, increases with stimulation intensity. OBJECTIVE: We used computational modeling to quantify the spatial effects of ICMS intensity and unify the apparently paradoxical findings of Histed and Stoney. METHODS: We implemented a biophysically-based computational model of a cortical column comprising neurons with realistic morphology and representative synapses. We quantified the spatial effects of single pulses and short trains of ICMS, including the volume of activated neurons and the density of activated neurons as a function of stimulation intensity. RESULTS: At all amplitudes, the dominant mode of somatic activation was by antidromic propagation to the soma following axonal activation, rather than via transsynaptic activation. There were no occurrences of direct activation of somata or dendrites. The volume over which antidromic action potentials were initiated grew with stimulation amplitude, while the volume of somatic activation increased marginally. However, the density of somatic activation within the activated volume increased with stimulation amplitude. CONCLUSIONS: The results resolve the apparent paradox between Stoney and Histed's results by demonstrating that the volume over which action potentials are initiated grows with ICMS amplitude, consistent with Stoney. However, the volume occupied by the activated somata remains approximately constant, while the density of activated neurons within that volume increase, consistent with Histed.

Characterizing the short-latency evoked response to intracortical microstimulation across a multi-electrode array.

Objective.Persons with tetraplegia can use brain-machine interfaces to make visually guided reaches with robotic arms. Without somatosensory feedback, these movements will likely be slow and imprecise, like those of persons who retain movement but have lost proprioception. Intracortical microstimulation (ICMS) has promise for providing artificial somatosensory feedback. ICMS that mimics naturally occurring neural activity, may allow afferent interfaces that are more informative and easier to learn than stimulation evoking unnaturalistic activity. To develop such biomimetic stimulation patterns, it is important to characterize the responses of neurons to ICMS.Approach.Using a Utah multi-electrode array, we recorded activity evoked by both single pulses and trains of ICMS at a wide range of amplitudes and frequencies in two rhesus macaques. As the electrical artifact caused by ICMS typically prevents recording for many milliseconds, we deployed a custom rapid-recovery amplifier with nonlinear gain to limit signal saturation on the stimulated electrode. Across all electrodes after stimulation, we removed the remaining slow return to baseline with acausal high-pass filtering of time-reversed recordings.Main results.After single pulses of stimulation, we recorded what was likely transsynaptically-evoked activity even on the stimulated electrode as early as ∼0.7 ms. This was immediately followed by suppressed neural activity lasting 10-150 ms. After trains, this long-lasting inhibition was replaced by increased firing rates for ∼100 ms. During long trains, the evoked response on the stimulated electrode decayed rapidly while the response was maintained on non-stimulated channels.Significance.The detailed description of the spatial and temporal response to ICMS can be used to better interpret results from experiments that probe circuit connectivity or function of cortical areas. These results can also contribute to the design of stimulation patterns to improve afferent interfaces for artificial sensory feedback.

A comprehensive model-based framework for optimal design of biomimetic patterns of electrical stimulation for prosthetic sensation.

OBJECTIVE: Touch and proprioception are essential to motor function as shown by the movement deficits that result from the loss of these senses, e.g. due to neuropathy of sensory nerves. To achieve a high-performance brain-controlled prosthetic arm/hand thus requires the restoration of somatosensation, perhaps through intracortical microstimulation (ICMS) of somatosensory cortex (S1). The challenge is to generate patterns of neuronal activation that evoke interpretable percepts. We present a framework to design optimal spatiotemporal patterns of ICMS (STIM) that evoke naturalistic patterns of neuronal activity and demonstrate performance superior to four previous approaches. APPROACH: We recorded multiunit activity from S1 during a center-out reach task (from proprioceptive neurons in Brodmann's area 2) and during application of skin indentations (from cutaneous neurons in Brodmann's area 1). We implemented a computational model of a cortical hypercolumn and used a genetic algorithm to design STIM that evoked patterns of model neuron activity that mimicked their experimentally-measured counterparts. Finally, from the ICMS patterns, the evoked neuronal activity, and the stimulus parameters that gave rise to it, we trained a recurrent neural network (RNN) to learn the mapping function between the physical stimulus and the biomimetic stimulation pattern, i.e. the sensory encoder to be integrated into a neuroprosthetic device. MAIN RESULTS: We identified ICMS patterns that evoked simulated responses that closely approximated the measured responses for neurons within 50 µm of the electrode tip. The RNN-based sensory encoder generalized well to untrained limb movements or skin indentations. STIM designed using the model-based optimization approach outperformed STIM designed using existing linear and nonlinear mappings. SIGNIFICANCE: The proposed framework produces an encoder that converts limb state or patterns of pressure exerted onto the prosthetic hand into STIM that evoke naturalistic patterns of neuronal activation.

Model-Based Evaluation of Closed-Loop Deep Brain Stimulation Controller to Adapt to Dynamic Changes in Reference Signal.

High-frequency deep brain stimulation (DBS) of the subthalamic nucleus (STN) is effective in suppressing the motor symptoms of Parkinson's disease (PD). Current clinically-deployed DBS technology operates in an open-loop fashion, i.e., fixed parameter high-frequency stimulation is delivered continuously, invariant to the needs or status of the patient. This poses two major challenges: (1) depletion of the stimulator battery due to the energy demands of continuous high-frequency stimulation, (2) high-frequency stimulation-induced side-effects. Closed-loop deep brain stimulation (CL DBS) may be effective in suppressing parkinsonian symptoms with stimulation parameters that require less energy and evoke fewer side effects than open loop DBS. However, the design of CL DBS comes with several challenges including the selection of an appropriate biomarker reflecting the symptoms of PD, setting a suitable reference signal, and implementing a controller to adapt to dynamic changes in the reference signal. Dynamic changes in beta oscillatory activity occur during the course of voluntary movement, and thus there may be a performance advantage to tracking such dynamic activity. We addressed these challenges by studying the performance of a closed-loop controller using a biophysically-based network model of the basal ganglia. The model-based evaluation consisted of two parts: (1) we implemented a Proportional-Integral (PI) controller to compute optimal DBS frequencies based on the magnitude of a dynamic reference signal, the oscillatory power in the beta band (13-35 Hz) recorded from model globus pallidus internus (GPi) neurons. (2) We coupled a linear auto-regressive model based mapping function with the Routh-Hurwitz stability analysis method to compute the parameters of the PI controller to track dynamic changes in the reference signal. The simulation results demonstrated successful tracking of both constant and dynamic beta oscillatory activity by the PI controller, and the PI controller followed dynamic changes in the reference signal, something that cannot be accomplished by constant open-loop DBS.