Modeling the volume of tissue activated in deep brain stimulation and its clinical influence: a review

Deep brain stimulation (DBS) is a neuromodulatory therapy that has been FDA approved for the treatment of various disorders, including but not limited to, movement disorders (e.g., Parkinson's disease and essential tremor), epilepsy, and obsessive-compulsive disorder. Computational methods for estimating the volume of tissue activated (VTA), coupled with brain imaging techniques, form the basis of models that are being generated from retrospective clinical studies for predicting DBS patient outcomes. For instance, VTA models are used to generate target-and network-based probabilistic stimulation maps that play a crucial role in predicting DBS treatment outcomes. This review defines the methods for calculation of tissue activation (or modulation) including ones that use heuristic and clinically derived estimates and more computationally involved ones that rely on finite-element methods and biophysical axon models. We define model parameters and provide a comparison of commercial, open-source, and academic simulation platforms available for integrated neuroimaging and neural activation prediction. In addition, we review clinical studies that use these modeling methods as a function of disease. By describing the tissue-activation modeling methods and highlighting their application in clinical studies, we provide the neural engineering and clinical neuromodulation communities with perspectives that may influence the adoption of modeling methods for future DBS studies.

Parkinson's Disease Tremor Differentially Responds to Levodopa and Subthalamic Stimulation

Background

Tremor in Parkinson's disease (PD) has an inconsistent response to levodopa and subthalamic deep brain stimulation (STN-DBS).

Objectives

To identify predictive factors of PD tremor responsiveness to levodopa and STN-DBS.

Material and Methods

PD patients with upper limb tremor who underwent STN-DBS were included. The levodopa responsiveness of tremor (overall, postural, and rest sub-components), was assessed using the relevant Unified Parkinson's Disease Rating Scale-III items performed during the preoperative assessment. Post-surgical outcomes were similarly assessed ON and OFF stimulation. A score for the rest/postural tremor ratio was used to determine the influence of rest and postural tremor severity on STN-DBS outcome. Factors predictive of tremor responsiveness were determined using multiple linear regression modeling. Volume of tissue activated measurement coupled to voxel-based analysis was performed to identify anatomical clusters associated with motor symptoms improvement.

Results

One hundred and sixty five patients were included in this study. Male gender was negatively correlated with tremor responsiveness to levodopa, whereas the ratio of rest/postural tremor was positively correlated with both levodopa responsiveness and STN-DBS tremor outcome. Clusters corresponding to improvement of tremor were in the subthalamic nucleus, the zona incerta and the thalamus, whereas clusters corresponding to improvement for akinesia and rigidity were located within the subthalamic nucleus.

Conclusion

More severe postural tremor and less severe rest tremor were associated with both poorer levodopa and STN-DBS response. The different locations of clusters associated with best correction of tremor and other parkinsonian features suggest that STN-DBS effect on PD symptoms is underpinned by the modulation of different networks.

Experimental Verification for Numerical Simulation of Thalamic Stimulation-Evoked Calcium-Sensitive Fluorescence and Electrophysiology with Self-Assembled Multifunctional Optrode

Owing to its capacity to eliminate a long-standing methodological limitation, fiber photometry can assist research gaining novel insight into neural systems. Fiber photometry can reveal artifact-free neural activity under deep brain stimulation (DBS). Although evoking neural potential with DBS is an effective method for mediating neural activity and neural function, the relationship between DBS-evoked neural Ca²⁺ change and DBS-evoked neural electrophysiology remains unknown. Therefore, in this study, a self-assembled optrode was demonstrated as a DBS stimulator and an optical biosensor capable of concurrently recording Ca²⁺ fluorescence and electrophysiological signals. Before the in vivo experiment, the volume of tissue activated (VTA) was estimated, and the simulated Ca²⁺ signals were presented using Monte Carlo (MC) simulation to approach the realistic in vivo environment. When VTA and the simulated Ca²⁺ signals were combined, the distribution of simulated Ca²⁺ fluorescence signals matched the VTA region. In addition, the in vivo experiment revealed a correlation between the local field potential (LFP) and the Ca²⁺ fluorescence signal in the evoked region, revealing the relationship between electrophysiology and the performance of neural Ca²⁺ concentration behavior. Concurrent with the VTA volume, simulated Ca²⁺ intensity, and the in vivo experiment, these data suggested that the behavior of neural electrophysiology was consistent with the phenomenon of Ca²⁺ influx to neurons.

Patient-specific structural connectivity informs outcomes of responsive neurostimulation for temporal lobe epilepsy

OBJECTIVE

Responsive neurostimulation is an effective therapy for patients with refractory mesial temporal lobe epilepsy. However, clinical outcomes are variable, few patients become seizure-free, and the optimal stimulation location is currently undefined. The aim of this study was to quantify responsive neurostimulation in the mesial temporal lobe, identify stimulation-dependent networks associated with seizure reduction, and determine if stimulation location or stimulation-dependent networks inform outcomes.

METHODS

We modeled patient-specific volumes of tissue activated and created probabilistic stimulation maps of local regions of stimulation across a retrospective cohort of 22 patients with mesial temporal lobe epilepsy. We then mapped the network stimulation effects by seeding tractography from the volume of tissue activated with both patient-specific and normative diffusion-weighted imaging. We identified networks associated with seizure reduction across patients using the patient-specific tractography maps and then predicted seizure reduction across the cohort.

RESULTS

Patient-specific stimulation-dependent connectivity was correlated with responsive neurostimulation effectiveness after cross-validation (p = .03); however, normative connectivity derived from healthy subjects was not (p = .44). Increased connectivity from the volume of tissue activated to the medial prefrontal cortex, cingulate cortex, and precuneus was associated with greater seizure reduction.

SIGNIFICANCE

Overall, our results suggest that the therapeutic effect of responsive neurostimulation may be mediated by specific networks connected to the volume of tissue activated. In addition, patient-specific tractography was required to identify structural networks correlated with outcomes. It is therefore likely that altered connectivity in patients with epilepsy may be associated with the therapeutic effect and that utilizing patient-specific imaging could be important for future studies. The structural networks identified here may be utilized to target stimulation in the mesial temporal lobe and to improve seizure reduction for patients treated with responsive neurostimulation.

Correlation between Electrode Location and Anxiety Depression of Subthalamic Nucleus Deep Brain Stimulation in Parkinson's Disease

Objectives: our group explored the correlation between postoperative coordinates of the electrode contacts, VTA, and anxiety and depression symptoms in Parkinson’s disease (PD) patients after subthalamic nucleus deep brain stimulation (STN-DBS). Methods: STN-DBS was conducted on PD patients (n = 57) for six months with follow-up. Clinical outcomes were explored using the unified Parkinson’s disease rating scale Part III (UPDRS-III), the Hamilton Anxiety Rating Scale (HAM-A), and the Hamilton Depression Rating Scale (HAM-D) before and after surgery. At the Montreal Neurological Institute (MNI), the location of active contacts and the volume of tissue activated (VTA) were calculated. Results: patient evaluations took place preoperatively and follow-ups took place at 1 month, 3 months, and 6 months. The average patient improvement rates for HAM-A and HAM-D scores at the 6-month follow-up were 41.7% [interquartile range (IQR) 34.9%] and 37.5% (IQR 33.4%), respectively (both p < 0.001). In medication-off, there were negative correlations between the HAM-A improvement rate and the Z-axis coordinate of the active contact (left side: r = −0.308, p = 0.020; right side: r = −0.390, p = 0.003), and negative correlations between the HAM-D improvement rate and the Z-axis coordinate of the active contact (left side: r = −0.345, p = 0.009; right side: r = −0.521, p = 0.001). There were positive correlations between the HAM-A and HAM-D scores improvement rate at 6 months after surgery and bilateral VTA in the right STN limbic subregion (HAM-A: r = 0.314, p = 0.018; HAM-D: r = 0.321, p = 0.015). Conclusion: bilateral STN-DBS can improve anxiety and depression symptoms in PD patients. The closer the stimulation to the ventral limbic region of the STN, the more significant the improvement in anxiety and depression symptoms of PD patients.

Subthalamic Stimulation for Camptocormia in Parkinson's Disease: Association of Volume of Tissue Activated and Structural Connectivity with Clinical Effectiveness

BACKGROUND

Deep brain stimulation of the subthalamic nucleus (STN-DBS) has been reported to be effective for camptocormia in Parkinson's disease (PD). However, the association between clinical effectiveness and the stimulated volumes or structural connectivity remains unexplored.

OBJECTIVE

To investigate the effectiveness of STN-DBS for treating camptocormia in PD and its association with volumes of tissue activated (VTA) and structural connectivity.

METHODS

We reviewed video recordings of patients who had undergone STN-DBS. The total and upper camptocormia (TCC and UCC) angles were measured to quantify changes in camptocormia. The Movement Disorders Society Unified Parkinson's Disease Rating Scale III (MDS-UPDRS III) was used to assess motor symptoms. Pre- and postoperative brain images were collected for modeling volume of VTA and structural connectivity using Lead-DBS software.

RESULTS

Participants included 36 patients with PD (8 with TCC-camptocormia and 2 with UCC-camptocormia) treated with bilateral STN-DBS. After surgery, patients showed a significant improvement in postural alignment at follow-up (mean follow-up duration: 6.0±2.2 months). In the entire sample, higher structural connectivity to the right supplementary motor area (SMA) and right lateral premotor cortex along the dorsal plane (PMd) was associated with larger postsurgical improvements in axial signs and TCC angles after stimulation was turned on. In patients diagnosed with camptocormia, larger improvement in camptocormia angles after STN-DBS was associated with a larger VTA overlap with STN (R = 0.75, p = 0.032).

CONCLUSION

This study suggests that both VTA overlap with STN and structural connectivity to cortical motor regions are associated with the effectiveness of STN-DBS for managing camptocormia in PD.

Long-Term Outcome and Neuroimaging of Deep Brain Stimulation in Holmes Tremor: A Case Series

BACKGROUND

Different deep brain stimulation (DBS) targets have been suggested as treatment for patients with pharmacologically refractory Holmes tremor (HT). We report the clinical and quality of life (QoL) long-term (up to nine years) outcome in four patients with HT treated with DBS (in thalamic ventral intermediate nucleus-VIM or in dentato-rubro-thalamic tract-DRTT).

MATERIALS AND METHODS

The patients underwent routine clinical evaluations before and after DBS (typically annually). Tremor severity and activities of daily living (ADL) were quantified by the Fahn-Tolosa-Marin Tremor-Rating-Scale (FTMTRS). QoL was assessed using the RAND SF-36-item Health Survey (RAND SF-36). In addition, we computed, in all four patients, the VTA based on the best stimulation settings using heuristic approaches included in the open source toolbox LEAD-DBS.

RESULTS

In all patients, tremor and ADL improved significantly at one-year post-DBS follow-up (34-61% improvement in FTMTRS total score compared to baseline). In three out of four patients, the improvement of tremor was sustained no longer than two to three years and only in one patient was sustained up to nine years. In this patient, the largest intersection between VTA and DBS target has been observed. Scores for ADL deteriorated over the course of time, reaching worse levels compared to baseline already during the three-year post-DBS follow-up, in three out of four patients. Physical and mental health component scores of RAND SF-36 had very different outcome between patients and follow-ups and were not associated with tremor-related outcomes.

CONCLUSIONS

The benefits of DBS in HT might not be always long lasting. Although QoL slightly improved, this change seemed to be independent of the motor outcome following DBS. The estimation of DBS target and VTA proximity could be a useful tool for DBS clinicians in order to facilitate the DBS programming process and optimize DBS treatment.

Selecting the Most Effective DBS Contact in Essential Tremor Patients Based on Individual Tractography

Postoperative choice of the most effective deep brain stimulation (DBS) contact in patients with essential tremor (ET) so far relies on lengthy clinical testing. Previous studies showed that the postoperative effectiveness of DBS contacts depends on the distance to the dentatorubrothalamic tract (DRTT). Here, we investigated whether the most effective DBS contact could be determined from calculating stimulation overlap with the individual DRTT. Seven ET patients with bilateral thalamic deep brain stimulation were included retrospectively. Tremor control was assessed for each contact during test stimulation with 2mA. Individual DRTTs were identified from diffusion tensor imaging and contacts were ranked by their stimulation overlap with the respective DRTT in relation to their clinical effectiveness. A linear mixed-effects model was calculated to determine the influence of the DRTT overlap on tremor control. In 92.9% of investigated DBS leads, the contact with the best clinical effect was the contact with the highest or second-highest DRTT-overlap. At the group level, the DRTT-overlap explained 26.7% of the variance in the clinical outcomes (p < 0.001). Our data suggest that the overlap with the DRTT based on individual tractography may serve as a marker to determine the most effective DBS contact in ET patients and reduce burdensome clinical testing in the future.

Cross-frequency coupling between gamma oscillations and deep brain stimulation frequency in Parkinson's disease

The disruption of pathologically enhanced beta oscillations is considered one of the key mechanisms mediating the clinical effects of deep brain stimulation on motor symptoms in Parkinson's disease. However, a specific modulation of other distinct physiological or pathological oscillatory activities could also play an important role in symptom control and motor function recovery during deep brain stimulation. Finely tuned gamma oscillations have been suggested to be prokinetic in nature, facilitating the preferential processing of physiological neural activity. In this study, we postulate that clinically effective high-frequency stimulation of the subthalamic nucleus imposes cross-frequency interactions with gamma oscillations in a cortico-subcortical network of interconnected regions and normalizes the balance between beta and gamma oscillations. To this end we acquired resting state high-density (256 channels) EEG from 31 patients with Parkinson's disease who underwent deep brain stimulation to compare spectral power and power-to-power cross-frequency coupling using a beamformer algorithm for coherent sources. To show that modulations exclusively relate to stimulation frequencies that alleviate motor symptoms, two clinically ineffective frequencies were tested as control conditions. We observed a robust reduction of beta and increase of gamma power, attested in the regions of a cortical (motor cortex, supplementary motor area, premotor cortex) and subcortical network (subthalamic nucleus and cerebellum). Additionally, we found a clear cross-frequency coupling of narrowband gamma frequencies to the stimulation frequency in all of these nodes, which negatively correlated with motor impairment. No such dynamics were revealed within the control posterior parietal cortex region. Furthermore, deep brain stimulation at clinically ineffective frequencies did not alter the source power spectra or cross-frequency coupling in any region. These findings demonstrate that clinically effective deep brain stimulation of the subthalamic nucleus differentially modifies different oscillatory activities in a widespread network of cortical and subcortical regions. Particularly the cross-frequency interactions between finely tuned gamma oscillations and the stimulation frequency may suggest an entrainment mechanism that could promote dynamic neural processing underlying motor symptom alleviation.

In silico Accuracy and Energy Efficiency of Two Steering Paradigms in Directional Deep Brain Stimulation

Background: In Deep Brain Stimulation (DBS), stimulation field steering is used to achieve stimulation spatial specificity, which is critical to obtain clinical benefits and avoid side effects. Multiple Independent Current Control (MICC) and Interleaving/Multi Stim Set (Interleaving/MSS) are two stimulation field steering paradigms in commercially available DBS systems. This work investigates the stimulation field steering accuracy and energy efficiency of these two paradigms in directional DBS.

Methods: Volumes of Tissue Activated (VTAs) were generated in silico using pulse widths of 60 μs and five pulse amplitude fractionalizations intended to steer the VTAs radially in 12° steps. For each fractionalization, VTAs were generated with nine pre-defined target radii. Stimulation field steering accuracy was assessed based on the VTAs rotation angle. Energy efficiency was inferred from current draw from battery values, which were calculated based on the pulse amplitudes needed to generate and steer the VTAs, as well as electrode impedance measurements of clinically implanted directional leads.

Results: For radial steering, MICC needed a single VTA. In contrast, Interleaving/MSS required the generation of two VTAs, whose union and intersection created an Interleaving/MSS VTA and an Intersection VTA, respectively. MICC VTAs were 6.8 (−3.2–11.8)% larger than Interleaving/MSS VTAs. The Intersection VTAs accounted for 26.2 (16.0–32.8)% of Interleaving/MSS VTAs and were exposed to a higher stimulation frequency. For all VTA radius-fractionalization combinations, steering accuracy was 7.0 (4.5–10.5)° for MICC and 24.0 (9.0–25.3)° for Interleaving/MSS. Pulse amplitudes were 16.1 (9.2–28.6)% lower for MICC than for Interleaving/MSS, leading to a 45.9 (18.8–72.6)% lower current draw from battery for MICC.

Conclusions: The results of this work show that in silico, MICC achieves a significantly better stimulation field steering accuracy and has a significantly higher energy efficiency than Interleaving/MSS. Although direct evidence still needs to be generated to translate the results of this work to clinical practice, clinical outcomes may profit from the better stimulation field steering accuracy of MICC and longevity of DBS systems may profit from its higher energy efficiency.

Selection of deep brain stimulation contacts using volume of tissue activated software following subthalamic nucleus stimulation

OBJECTIVE

High-frequency deep brain stimulation (DBS) of the subthalamic nucleus (STN) is effective in the treatment of motor symptoms of Parkinson's disease. Using a patient-specific lead and volume of tissue activated (VTA) software, it is possible to visualize contact positions in the context of the patient's own anatomy. In this study, the authors' aim was to demonstrate that VTA software can be used in clinical practice to help determine the clinical effectiveness of stimulation in patients with Parkinson's disease undergoing DBS of the STN.

METHODS

Brain images of 26 patients undergoing STN DBS were analyzed using VTA software. Preoperative clinical and neuropsychological data were collected. Contacts were chosen by two experts in DBS blinded to the clinical data. A therapeutic window of amplitude was determined. These results were compared with the parameter settings for each patient. Data were obtained at 3 months and 1 year postsurgery.

RESULTS

In 90.4% (95% CI 82%-98%) of the patients, the contacts identified by the VTA software were concordant with the clinically effective contacts or with an effective contact in contact-by-contact testing. The therapeutic window of amplitude selected virtually included 81.3% of the clinical amplitudes.

CONCLUSIONS

VTA software appears to present significant concordance with clinical data for selecting contacts and stimulation parameters that could help in postoperative follow-up and programming.

Modeling of Electric Fields in Individual Imaging Atlas for Capsular Threshold Prediction of Deep Brain Stimulation in Parkinson's Disease: A Pilot Study

Background: Modeling of deep brain stimulation electric fields and anatomy-based software might improve post-operative management of patients with Parkinson's disease (PD) who have benefitted from subthalamic nucleus deep brain stimulation (STN-DBS).

Objective: We compared clinical and software-guided determination of the thresholds for current diffusion to the pyramidal tract, the most frequent limiting side effect in post-operative management of STN-DBS PD patients.

Methods: We assessed monopolar reviews in 16 consecutive STN-DBS PD patients and retrospectively compared clinical capsular thresholds, which had been assessed according to standard clinical practice, to those predicted by volume of tissue activated (VTA) model software. All the modeling steps were performed blinded from patients' clinical evaluations.

Results: At the group level, we found a significant correlation (p = 0.0001) when performing statistical analysis on the z-scored capsular thresholds, but with a low regression coefficient (r = 0.2445). When considering intra-patient analysis, we found significant correlations (p < 0.05) between capsular threshold as modeled with the software and capsular threshold as determined clinically in five patients (31.2%).

Conclusions: In this pilot study, the VTA model software was of limited assistance in identifying capsular thresholds for the whole cohort due to a large inter-patient variability. Clinical testing remains the gold standard in selecting stimulation parameters for STN-DBS in PD.

Non-motor outcomes depend on location of neurostimulation in Parkinson's disease

Deep brain stimulation of the subthalamic nucleus is an effective and established therapy for patients with advanced Parkinson's disease improving quality of life, motor symptoms and non-motor symptoms. However, there is a considerable degree of interindividual variability for these outcomes, likely due to variability in electrode placement and stimulation settings. Here, we present probabilistic mapping data from a prospective, open-label, multicentre, international study to investigate the influence of the location of subthalamic nucleus deep brain stimulation on non-motor symptoms in patients with Parkinson's disease. A total of 91 Parkinson's disease patients undergoing bilateral deep brain stimulation of the subthalamic nucleus were included, and we investigated NMSScale, NMSQuestionnaire, Scales for Outcomes in Parkinson's disease-motor examination, -activities of daily living, and -motor complications, and Parkinson's disease Questionnaire-8 preoperatively and at 6-month follow-up after surgery. Leads were localized in standard space using the Lead-DBS toolbox and individual volumes of tissue activated were calculated based on clinical stimulation settings. Probabilistic stimulation maps and non-parametric permutation statistics were applied to identify voxels with significant above or below average improvement for each scale and analysed using the DISTAL atlas. All outcomes improved significantly at follow-up. Significant spatial distribution patterns of neurostimulation were observed for NMSScale total score and its mood/apathy and attention/memory domains. For both domains, voxels associated with below average improvement were mainly located dorsal to the subthalamic nucleus. In contrast, above average improvement for mood/apathy was observed in the ventral border region of the subthalamic nucleus and in its sensorimotor subregion and for attention/memory in the associative subregion. A trend was observed for NMSScale sleep domain showing voxels with above average improvement located ventral to the subthalamic nucleus. Our study provides evidence that the interindividual variability of mood/apathy, attention/memory, and sleep outcomes after subthalamic nucleus deep brain stimulation depends on the location of neurostimulation. This study highlights the importance of holistic assessments of motor and non-motor aspects of Parkinson's disease to tailor surgical targeting and stimulation parameter settings to patients' personal profiles.

Interleaved deep brain stimulation for dyskinesia management in Parkinson's disease

BACKGROUND

In patients with Parkinson's disease, stimulation above the subthalamic nucleus (STN) may engage the pallidofugal fibers and directly suppress dyskinesia.

OBJECTIVES

The objective of this study was to evaluate the effect of interleaving stimulation through a dorsal deep brain stimulation contact above the STN in a cohort of PD patients and to define the volume of tissue activated with antidyskinesia effects.

METHODS

We analyzed the Core Assessment Program for Surgical Interventional Therapies dyskinesia scale, Unified Parkinson's Disease Rating Scale parts III and IV, and other endpoints in 20 patients with interleaving stimulation for management of dyskinesia. Individual models of volume of tissue activated and heat maps were used to identify stimulation sites with antidyskinesia effects.

RESULTS

The Core Assessment Program for Surgical Interventional Therapies dyskinesia score in the on medication phase improved 70.9 ± 20.6% from baseline with noninterleaved settings (P < 0.003). With interleaved settings, dyskinesia improved 82.0 ± 27.3% from baseline (P < 0.001) and 61.6 ± 39.3% from the noninterleaved phase (P = 0.006). The heat map showed a concentration of volume of tissue activated dorsally to the STN during the interleaved setting with an antidyskinesia effect.

CONCLUSION

Interleaved deep brain stimulation using the dorsal contacts can directly suppress dyskinesia, probably because of the involvement of the pallidofugal tract, allowing more conservative medication reduction. © 2019 International Parkinson and Movement Disorder Society.