Model-Based Analysis of Pathway Recruitment During Subthalamic Deep Brain Stimulation

BACKGROUND

Subthalamic deep brain stimulation (DBS) is an established clinical therapy, but an anatomically clear definition of the underlying neural target(s) of the stimulation remains elusive. Patient-specific models of DBS are commonly used tools in the search for stimulation targets, and recent iterations of those models are focused on characterizing the brain connections that are activated by DBS.

OBJECTIVE

The goal of this study was to quantify axonal pathway activation in the subthalamic region from DBS at different electrode locations and stimulation settings.

MATERIALS AND METHODS

We used an anatomically and electrically detailed computational model of subthalamic DBS to generate recruitment curves for eight different axonal pathways of interest, at three generalized DBS electrode locations in the subthalamic nucleus (STN) (ie, central STN, dorsal STN, posterior STN). These simulations were performed with three levels of DBS electrode localization uncertainty (ie, 0.5 mm, 1.0 mm, 1.5 mm).

RESULTS

The recruitment curves highlight the diversity of pathways that are theoretically activated with subthalamic DBS, in addition to the dependence of the stimulation location and parameter settings on the pathway activation estimates. The three generalized DBS locations exhibited distinct pathway recruitment curve profiles, suggesting that each stimulation location would have a different effect on network activity patterns. We also found that the use of anodic stimuli could help limit activation of the internal capsule relative to other pathways. However, incorporating realistic levels of DBS electrode localization uncertainty in the models substantially limits their predictive capabilities.

CONCLUSIONS

Subtle differences in stimulation location and/or parameter settings can impact the collection of pathways that are activated during subthalamic DBS.

Quantifying the Variability Associated with Postoperative Localization of Deep Brain Stimulation Electrodes

INTRODUCTION

Computational models of deep brain stimulation (DBS) have become common tools in clinical research studies that attempt to establish correlations between stimulation locations in the brain and behavioral outcome measures. However, the accuracy of any patient-specific DBS model depends heavily upon accurate localization of the DBS electrodes within the anatomy, which is typically defined via co-registration of clinical CT and MRI datasets. Several different approaches exist for this challenging registration problem, and each approach will result in a slightly different electrode localization. The goal of this study was to better understand how different processing steps (e.g., cost-function masking, brain extraction, intensity remapping) affect the estimate of the DBS electrode location in the brain.

METHODS

No "gold standard" exists for this kind of analysis, as the exact location of the electrode in the living human brain cannot be determined with existing clinical imaging approaches. However, we can estimate the uncertainty associated with the electrode position, which can be used to guide statistical analyses in DBS mapping studies. Therefore, we used high-quality clinical datasets from 10 subthalamic DBS subjects and co-registered their long-term postoperative CT with their preoperative surgical targeting MRI using 9 different approaches. The distances separating all of the electrode location estimates were calculated for each subject.

RESULTS

On average, electrodes were located within a median distance of 0.57 mm (0.49-0.74) of one another across the different registration approaches. However, when considering electrode location estimates from short-term postoperative CTs, the median distance increased to 2.01 mm (1.55-2.78).

CONCLUSIONS

The results of this study suggest that electrode location uncertainty needs to be factored into statistical analyses that attempt to define correlations between stimulation locations and clinical outcomes.

Comparison of methodologies for modeling directional deep brain stimulation electrodes

Deep brain stimulation (DBS) is an established clinical therapy, and directional DBS electrode designs are now commonly used in clinical practice. Directional DBS leads have the ability to increase the therapeutic window of stimulation, but they also increase the complexity of clinical programming. Therefore, computational models of DBS have become available in clinical software tools that are designed to assist in the identification of therapeutic settings. However, the details of how the DBS model is implemented can influence the predictions of the software. The goal of this study was to compare different methods for representing directional DBS electrodes within finite element volume conductor (VC) models. We evaluated 15 different DBS VC model variants and quantified how their differences influenced estimates on the spatial extent of axonal activation from DBS. Each DBS VC model included the same representation of the brain and head, but the details of the current source and electrode contact were different for each model variant. The more complex VC models explicitly represented the DBS electrode contacts, while the more simple VC models used boundary condition approximations. The more complex VC models required 2-3 times longer to mesh, build, and solve for the DBS voltage distribution than the more simple VC models. Differences in individual axonal activation thresholds across the VC model variants were substantial (-24% to +47%). However, when comparing total activation of an axon population, or estimates of an activation volume, the differences between model variants decreased (-7% to +8%). Nonetheless, the technical details of how the electrode contact and current source are represented in the DBS VC model can directly affect estimates of the voltage distribution and electric field in the brain tissue.

Deep brain stimulation of terminating axons

BACKGROUND

Deep brain stimulation (DBS) of the subthalamic region is an established treatment for the motor symptoms of Parkinson's disease. Several types of neural elements reside in the subthalamic region, including subthalamic nucleus (STN) neurons, fibers of passage, and terminating afferents. Recent studies suggest that direct activation of a specific population of subthalamic afferents, known as the hyperdirect pathway, may be responsible for some of the therapeutic effects of subthalamic DBS.

OBJECTIVE

The goal of this study was to quantify how axon termination affects neural excitability from DBS. We evaluated how adjusting different stimulation parameters influenced the relative excitability of terminating axons (TAs) compared to fibers of passage (FOPs).

METHODS

We used finite element electric field models of DBS, coupled to multi-compartment cable models of axons, to calculate activation thresholds for populations of TAs and FOPs. These generalized models were used to evaluate the response to anodic vs. cathodic stimulation, with short vs. long stimulus pulses.

RESULTS

Terminating axons generally exhibited lower thresholds than fibers of passage across all tested parameters. Short pulse widths accentuated the relative excitability of TAs over FOPs.

CONCLUSION(S)

Our computational results demonstrate a hyperexcitability of terminating axons to DBS that is robust to variation in the stimulation parameters, as well as the axon model parameters.