Linking Profiles of Pathway Activation with Clinical Motor Improvements – a Retrospective Computational Study =============================================================================================================== * Konstantin Butenko * Ningfei Li * Clemens Neudorfer * Jan Roediger * Andreas Horn * Gregor R. Wenzel * Hazem Eldebakey * Andrea A. Kühn * Martin Reich * Jens Volkmann * Ursula van Rienen ## Abstract Deep brain stimulation (DBS) is an established therapy for patients with Parkinson’s disease. *In silico* computer models for DBS allow to pre-select a set of potentially optimal stimulation parameters. If efficacious, they could further carry insight into the mechanism of action of DBS and foster the development of more efficient stimulation approaches. In recent years, the focus has shifted towards DBS-induced firing in myelinated axons, deemed particularly relevant for the external modulation of neural activity. We use the concept of pathway activation modeling, which incorporates advanced volume conductor models and anatomically authentic fiber trajectories to estimate DBS-induced action potential initiation in anatomically plausible pathways that traverse in close proximity to targeted nuclei. We apply the method on a retrospective dataset with the aim of providing a model-based prediction of clinical improvement following DBS (as measured by the motor part of the Unified Parkinson’s Disease Rating Scale). Based on differences in outcome and activation rates for two DBS protocols in a training cohort, we compute a theoretical 100% improvement profile and enhance it by analyzing the importance of profile matching for individual pathways. Finally, we validate the performance of our profile-based predictive model in a test cohort. As a result, we demonstrate the clinical utility of pathway activation modeling in the context of motor symptom alleviation in Parkinson’s patients treated with DBS. ## Introduction Deep brain stimulation (DBS) is an effective treatment for various neurological and mental disorders, and it has become an established therapy for patients suffering from therapy-refractory Parkinson’s disease (PD). During DBS, short high-frequency pulses are delivered to subcortical brain structures via implanted electrodes. These electrodes usually have 4 or 8 contacts, each of which can be used as a current source. Modern DBS systems allow great flexibility in pulse modulation, including adjustment of width, amplitude, and frequency. Determining an optimal stimulation protocol in such a large parameter space is challenging, and to assist medical professionals in this procedure, *in silico* computational models for DBS could be of use. Besides, such models could provide insights into the action mechanism of the treatment, which in turn could drive the development of more efficient and effective stimulation paradigms. In the basal ganglia, DBS has been hypothesized to create an “informational lesion” partly facilitated by electrical stimulation of highly excitable myelinated axons [21] that form projections across DBS targets, such as the subthalamic nucleus (STN) and the globus pallidus internus (GPi). To quantify the effect, different concepts have been proposed, among which the volume of tissue activated [8] or its approximations [13] are most commonly employed. Multiple studies have investigated correlations of symptom alleviation with voxel metrics defined based on the volume of tissue activated and its interaction with structural and functional connectivity [28, 29, 38, 43]. More recently, the concept of pathway activation modeling was proposed [23] that comprises advanced volume conductor models and anatomically authentic fiber trajectories to estimate a DBS-induced action potential initiation along pathways residing in the vicinity of the stimulation targets. To date, studies on Tourette syndrome, obsessive-compulsive disorder, treatment-resistant depression, and PD have successfully employed pathway activation modeling to predict symptom alleviation in patients [20, 24, 30, 34]. However, to our best knowledge, the method has not been applied on a cohort level to investigate correlations between DBS-induced axonal activation and alleviation of aggregated motor symptoms in PD. In this retrospective computational study, we apply pathway activation modeling to identify the network correlates underlying the improvement of the Unified Parkinson’s Disease Rating Scale III score (motor examination, further referred to as UPRDS-III) and the Movement Disorder Society (MDS) sponsored revision of UPRDS-III in patients suffering from therapy-refractory PD. We base this prediction on pathway activation profiles defined by multiple pathways recruited by Subthalamic Nucleus Deep Brain Stimulation (STN-DBS). Correlating symptom alleviation with profiles instead of activation in individual pathways can be a more robust metric considering compensatory and adverse effects of fiber recruitment in the vicinity of the STN. Based on interprotocol scores, defined as the difference of two DBS-on UPDRS-III scores assessed in the training cohort for each patient, we construct a theoretical 100% UPDRS-III improvement profile and enhance it by analyzing the significance of activation levels in individual pathways. The performance of the resulting profile-based predictive model is then successfully tested in an independent cohort. ## Materials and Methods ### Patient Cohorts and Imaging Two cohorts from independent DBS centers, namely Charité – Universitätsmedizin Berlin and Würzburg University Hospital, were retrospectively analyzed for derivation and validation of a pathway activation-based predictive model. The cohorts were formed based on the following criteria: * availability of medical imaging data, such as preoperative T1-weighted magnetic resonance imaging (MRI) and postoperative computed tomography (CT), necessary for electrode localization and patient-specific modeling; * current-controlled stimulation mode, which allows to compensate for a voltage drop on the electrode-tissue interface, reducing the computational model complexity; * placement of the electrode inside the STN with a minimal penetration of the pallidus. The former restriction is imposed since only activation in pathways in the vicinity of the STN was investigated, and the latter criterion was applied to reduce effects of the pallidal lesioning. As a result, 15 patients for the training (Berlin) and 19 for the test (Würzburg) cohort were admitted. Although the first cohort contained fewer patients, it provided a total of 30 datapoints: two stimulation protocols were documented for each patient as a part of another study on an algorithm-guided DBS-programming [53]. Later, this aspect was used to derive the predictive model. In both cohorts, the patients constituted a representative sampling of a clinical PD-DBS cohort (see Table 1). All received octopolar DBS electrodes bilaterally to STN (Fig. 1), either using omnidirectional or directional DBS electrodes (Boston Scientific Vercise™, Marlborough, MA, USA). Stimulation was performed using a conventional current-controlled DBS signal: a rectangular pulse of 20–60 *µs* length, 79–185 *Hz* repetition rate, each followed by an extended charge balancing period at low amplitude. Motor performance in the training and the test cohorts was evaluated with either UPDRS-III or MDS-UPDRS-III, respectively. These scores are strongly (*r >* 0.95) correlated [42], and, for brevity, both will be referred to as UPRDS-III. They were taken at baseline (off medication, DBS-off) preoperatively in the test cohort, and at least 6 months after surgery (same day as DBS-on scores) in the training cohort. UPDRS-III scores under active DBS (off medication) were acquired at least after 6 months in both cohorts. View this table: [Table 1.](http://medrxiv.org/content/early/2021/10/11/2021.10.08.21264743/T1) Table 1. Cohort Demographics and Clinical Outcomes. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/11/2021.10.08.21264743/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/10/11/2021.10.08.21264743/F1) Figure 1. Reconstruction of DBS electrodes in the training (Berlin, A) and the test (Würzburg, B) cohorts. Electrodes are visualized in standard stereotactic Montreal Neurological Institute (MNI) space using Lead-DBS. Displayed brain structures are defined by the DISTAL atlas [15] and include subthalamic nucleus (orange), globus pallidus externus (blue) and globus pallidus internus (green). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/11/2021.10.08.21264743/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/10/11/2021.10.08.21264743/F2) Figure 2. Dataflow for computing pathway activation. Based on patient imaging and brain atlases, Lead-DBS (orange box) reconstructs the electrode and provides a description of tissueand water diffusion distributions in the brain. These data are used by OSS-DBS (green box) to create an accurate patient-specific volume conductor model. The model is then employed to compute electric potential distribution in space and time along axon models allocated on trajectories described by a pathway atlas. Finally, for the given distribution, the cable equation is solved to probe axonal activation, i.e., occurrence of an action potential in response to DBS. Imaging data was processed using Lead-DBS software [26, 27] ([lead-dbs.org](http://lead-dbs.org)). Postoperative volumes and preoperative weighted multimodal MRI scans were linearly co-registered using SPM12 [16] ([fil.ion.ucl.ac.uk/spm](https://fil.ion.ucl.ac.uk/spm/)). This was followed by a non-linear normalization step of co-registered patient scans to stereotactic Montreal Neurological Institute (MNI) space using Advanced Normalization Tools (ANTs, [stnava.github.io/ANTs/](http://stnava.github.io/ANTs/)) ’SyS’ algorithm [5]. We also accounted for potential non-linear displacements introduced by brain shift using an additional refinement of the co-registration step that focused on the subcortex. Electrodes were localized based on CT scans using PaCER [32], and results were visually evaluated and refined, if necessary. Bioelectrical effects of DBS were modeled in patient-specific (native) space. For that, patient-specific brain tissue segmentations (grey matter, white matter and cerebrospinal fluid) were obtained based on tissue probabilistic mapping [16] derived from the multispectral MNI template and T1-weighted MRI images. To account for tissue anisotropy, the mean diffusion tensor data of the human brain [56] were transformed into patient-specific space using the inverse deformation field derived from the diffeomorphic normalization procedure. The same procedure was applied to a basal ganglia pathway atlas [47] that describes fiber distribution and classification necessary for pathway activation modeling. Created under the guidance of expert neuroanatomists, it provides an accurate description of axonal trajectories affected by STN stimulation. ### Field and Pathway Activation Modeling Isotropic conductivity values for grey and white matter were evaluated in the frequency domain based on [17], but omitting *α*-dispersion as proposed in [57] and additionally upscaled. The conductivity *σ* of non-dispersive cerebrospinal fluid was fixed to 2.0 *S/m*. In the power spectrum of a conventional DBS signal (up to 1 *MHz*), the conductivities were monotonically increasing (*σ**grey*: 0.168–0.235 *S/m* and *σ**white*: 0.120–0.153 *S/m*). To account for this dispersion, the Fourier Finite Element Method [8] was applied to solve the quasistatic formulation of Maxwell’s equations that describes the spatial distribution of the electric potential *φ*(**r**): ![Formula][1] where *ω* = 2*πf* is the angular frequency of one of the harmonics that compose the DBS signal. Capacitive properties of brain tissue can be neglected in Eq. 1 due to its relatively low contribution after omitting *α*-dispersion. An octave band approximation method [6] was used to reduce the number of computations in the frequency domain. The solution in the time domain is then retrieved using an Inverse Fourier Transform. Anisotropy, especially prominent in white matter tracts [18], was modeled by expressing conductivity in terms of tensors defined according to the mean diffusion tensor data that were normalized voxel-wise following the volume conservation approach [22] and scaled by the isotropic conductivity of brain tissue. The electrode-tissue interface was neglected assuming its minor effect on current-controlled stimulation [8] with a charge density per phase below 0.03 *mC/cm*2. Nevertheless, the electrode’s encapsulation layer was accounted for by removing axons within a 0.1 *mm* vicinity, where neuron degeneration and glial scarring occur. The accuracy of the Finite Element Method computations was controlled based on the convergence of the electric field and the current, and elements with large deviations were refined. The obtained distribution the of extracellular potential in space and time was used to solve a double cable equation of a myelinated axon model described in [41]. The models were allocated on the trajectories delineated in [47] (for passing fibers, the closest point on the trajectory to active contacts was treated as the midpoint seed), and the pathway activation rate was computed as a fraction of axons in the pathway that elicited an action potential in response to the DBS pulse. In the present study, we modeled activation in the corticofugal and the hyperdirect pathways (HDP) originating in the primary and premotor cortex, as well as the supplementary motor area. To reduce the computational costs, the corticofugal pathway was uniformly downsampled from 5000 to 1250 streamlines. In addition, activation in the pallidosubthalamic pathways (sensorimotor portion) and the pallidothalamic pathways (ansa lenticularis and lenticular fasciculus) was investigated. Besides, we computed the extent of direct DBS recruitment of the passing cerebellothalamic tract associated with tremor suppression [10]. In total, 4000 axon models of 16 pathways per hemisphere were deployed for pathway activation modeling, with some axons being later removed due to their intersection with the electrode, the encapsulation layer or cerebrospinal fluid. In the present study, the fiber diameters and the number of nodes of Ranvier, which together defined the axonal length, were fixed for axons within one pathway. A fiber diameter of 3.0 *µm* was chosen for axons of the local STN-GPe and STN-GPi pathways with a length of 10 *mm* and 6 *mm*, respectively (35 nodes and 21 nodes of Ranvier). For other pathways, the fiber diameter was set to 5.7 *µm*, which is on a larger side for data reported in [40]. The same axon model was employed in [31], where it was noted that larger fiber diameters improve predictability of evoked potentials. In a preliminary analysis, we also tested a 12.0 *µm* fiber diameter as suggested by the authors, but this setup predicted a high activation in the corticofugal pathway, which is unlikely for clinically accepted protocols due to evoked motor contractions. For computational reasons, the length of the axons of passage and the HDP was limited to 20 *mm* (40 nodes of Ranvier), thus not covering the whole length of the corresponding projections. Nevertheless, this truncation is acceptable due to a minor effect of DBS on distant axonal compartments, and a preliminary analysis with longer axons yielded the same activation rates. All steps described in this subsection were carried out using the open-source simulation software OSS-DBS [7] that was developed for highly automated DBS modeling in heterogeneous anisotropic and dispersive volume conductor models. For the present study, the software was implemented as a computational backend, freely distributed along with Lead-DBS. The coupling allowed a seamless transition between state-of-the-art processing of medical imaging and accurate modeling of DBS-induced electric field with subsequent quantification of its effect on neural tissue using cable models. The OSS-DBS computations were encapsulated in Docker containers ([docker.com/](http://docker.com/)) deployed on Intel Xeon(R) Gold 6136 CPU @ 3.00 GHz x 48 cores with 376.6 GB of memory. ### Analysis of Pathway Activation Profiles As a first exploratory step, we calculated correlations between activation rates of individual pathways and the UPDRS-III improvement from baseline in the training cohort. We employed the Pearson correlation coefficient, assuming a linear interaction between the quantities. The Spearman rank correlation was comparable for all tests presented in this study. Secondly, we analyzed activation rate profiles **AR** (vector quantities composed of activation rates over all 16 pathways) and their correlation with alleviation of motor symptoms. Since the training cohort contained two DBS-on datapoints per patient, we decided to use these measurements to derive a vector of the interprotocol UPDRS-III improvement in the pathway activation space. The difference in UPDRS-III scores between two protocols allowed us to compute a distance to a theoretical pathway activation profile of 100% improvement, i.e., an optimal activation rate profile: ![Formula][2] where subscripts “+” and ” − “refer to the better and the worse performing stimulation protocols, respectively. Such a definition of the optimal pathway activation profile is preferable over derivations using baseline, where only endogenous activity is present that is not quantified in the model and assumed to be overwritten by the DBS-induced activation of fiber tracts. The difference between two profiles was quantified with the Canberra distance [36] ![Formula][3] where *w**i* is the weighting factor for the *i* th pathway, whose default value is 1.0. The Canberra distance was chosen based on a preliminary analysis within the training cohort where it showed the highest predictive ability in comparison with other metrics, such as Bray-Curtis dissimilarity and Euclidean distance. It is noteworthy that this normalized metric demonstrated the best performance, thus sparing us the argument of whether activation rates (percent of fibers present in the tract) or absolute activation (amount of fibers activated in each tract) should be employed when analyzing pathway activation results. The next problem was to determine which patients, i.e., pairs of datapoints, to use to derive **AR**optimal. At first glance, the patient with the best improvement between two DBS-on protocols would be a good candidate. However, the best improvement does not guarantee the shortest path to the optimal activation rate profile. Instead, we considered patients who had the highest interprotocol UPDRS-III improvement normalized by the mean of the difference vector (**AR**+ − **AR**−). To increase the robustness of the predictive model, **AR**optimal can be averaged among several patients. In the present study, we picked three patients with a prominently higher normalized UPDRS-III improvement. Additionally, we tested the predictive ability of **AR**optimal derived from all datapoints in the training cohort using a “leave-one-out” approach. Previous experimental and clinical studies showed that stimulation of specific pathways in the vicinity of the STN associated with symptom alleviation and occurrence of side-effects [9, 10, 48, 51]. Therefore, we presume that the importance of matching the optimal profile is not uniform among pathways. To test this hypothesis, we optimized the weighting factor *w**i* for each pathway using the training cohort, but excluding those recruited for **AR**optimal. The optimization was conducted according to the following procedure. First, a patient was excluded, and the rest were randomly shuffled and split. Next, for one set, the weighting factors were optimized to maximize the inverse correlation of the Canberra distance and the UPDRS-III improvement, both from baseline and between protocols. The second set was then used to test the predictive ability of the weighted metric. The procedure was repeated for all patients. Removal of a patient and random shuffling was employed to estimate convergence of the optimal weighting factors across the cohort, and multiple trials yielded no significant disparities introduced by the shuffling. Note that the limited number of datapoints overall necessitated use of the same patients (but in different combinations) in both sets. ## Results ### Pathway Activation and UPDRS-III Improvement Pathway activation rates computed for both cohorts are presented in Fig. 3. Analysis of activation in individual pathways and UPDRS-III improvement from baseline revealed only one pathway with a statistically significant positive correlation, namely, the HDP branch to the face-neck region of the primary motor cortex. However, it was not predictive when later applied to the test cohort. Additionally, we observed that increase in activation in the cerebellothalamic pathway correlated with interprotocol UPDRS improvement of tremor (see Suppl. 1). Interestingly, there was no correlation with the UPDRS-III improvement from baseline, possibly due to irrelevance of the metric activation for akinetic-rigid patients and contribution of the subthalamic projections to the tremor alleviation. The weak correlations for individual pathways might indicate an interplay of activation in different pathways. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/11/2021.10.08.21264743/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2021/10/11/2021.10.08.21264743/F3) Figure 3. Violin plots of pathway activation rates across both cohorts (two stimulation protocols for each patient in the training cohort) visualized with the correspnding pathways in Lead-DBS in MNI space. MC and PMC refer to the primary and premotor cortical regions, SMA – supplementary motor area; cf and hdp are the corticofugal and the hyperdirect pathways, respectively; Ansa – ansa lenticularis, Lent – lenticular fasciculus, CbTh – cerebellothalamic pathway; l, f, up - lower extremity, face-neck region, and upper extremity in the primary motor cortex. For clarity, only one direction of the reciprocal sensorimotor pallidosubthalamic projection is shown. Note that the cohorts’ datapoints provide comprehensive coverage of the pathway activation space. Fig. 4, A shows activation profiles of the best responders, i.e., patients with the highest UPRDS improvement from baseline, in both cohorts and the theoretical 100% improvement profile, derived from datapoints of three patients (from the training cohort) with the highest normalized UPDRS-III improvement. The Canberra distance of the datapoints from the rest of the training cohort to the optimal profile showed a significant inverse correlation with the correponding UPDRS-III improvements from baseline (see Fig. 5, A). In contrast, the optimal profile derived from all datapoints of the training cohort was not predictive (see Suppl. 2). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/11/2021.10.08.21264743/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2021/10/11/2021.10.08.21264743/F4) Figure 4. Theoretical 100% UPDRS-III improvement optimal profile. A: activation rates for the optimal profile and the best responders from each cohort. Note that the match between these datapoints and the optimal profile occurs for different pathways, suggesting that two distinct DBS mechanisms might be present. The mean values of optimized weighting factors below the bars can be interpreted as the importance of matching the optimal profile for a particular pathway. Red boxes highlight pathways with the highest effect on predictive ability. B: box plots of the optimized weighting factors computed with the ’leave-one-out’ approach. Note the low number of outliers. ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/11/2021.10.08.21264743/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2021/10/11/2021.10.08.21264743/F5) Figure 5. Correlation of UPDRS-III improvement in the training cohort with the distance to the optimal profile defined in the pathway activation space and based on three patients with the highest normalized interprotocol UPDRS-III improvement. Colored lines between datapoints depict interprotocol binary correlation (green – valid prediction, black – false.) ### Key Pathways of the Activation Profile Optimization results for the weighting factors *w**i* in the Canberra distance (see Eq. 3), conducted on datapoints of the training cohort, indicated that DBS-induced activation in specific pathways might play a prominent role in motor symptom alleviation (see Fig. 4, B). In particular, the optimization emphasized a moderate activation in both pallidothalamic pathways as well as two branches of the HDP to the primary motor cortex, while avoiding stimulation of the corticofugal pathway descending from the supplementary motor area. The weighted Canberra distances from activation profiles to **AR**optimal for patients who were successively left out showed a statistically significant correlation with UPDRS-III improvement from baseline. Notably, the optimized weights were comparable among patients and were averaged to draw a general conclusion on the importance of particular pathways, where deviations in activation rates from the optimal profile led to worse performance. The averagely weighted Canberra distance to **AR**optimal demonstrated a higher correlation with UPDRS-III improvement from baseline than the unweighted metric (Fig. 5, B), which is not, however, surprising, since the optimization was conducted on these datapoints. Therefore, the profile-based predictive model had to be validated on unseen patients from the hold-out test cohort. Critically, patients from the test cohort did not play a role in the model selection process and were entirely naive to the final metric derived from the training cohort. This metric, the weighted Canberra distance to **AR**optimal, defined using the training cohort, was predictive of UPDRS-III improvement from baseline in patients of the test cohort, while the unweighted distance was not (see Fig. 6.) Distribution of activation in the highly weighted pathways computed for the best responders is shown in Fig. 7. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/11/2021.10.08.21264743/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2021/10/11/2021.10.08.21264743/F6) Figure 6. Cross-predicting UPDRS-III improvement in the test cohort using the optimal activation rate profile derived from datapoints of the three patients in the training cohort. Note that only the weighted Canberra distance has a statistically significant correlation with the improvement. ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/11/2021.10.08.21264743/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2021/10/11/2021.10.08.21264743/F7) Figure 7. Pathway activation visualized in Lead-DBS. A: pathways which activation rates determine the UPDRS-III improvement (shown in MNI space.) Blue and cyan – the hyperdirect pathway, descending from face-neck and upper extremity of the primary motor cortex, red – the corticofugal pathway, descending from the supplementary motor area, green – lenticular fasciculus, orange – ansa lenticularis. B and C: pathway activation (shown on computational axon models in patient-specific space) for the best responders in the training and the test cohorts, respectively. Axons closer than 0.1 mm from the electrode (purple) are excluded due to glial scaring and neurodegeneration. Note the high recruitment of the HDP and the corticofugal pathway for the best responder in the test cohort. ## Discussion In this study, we conducted pathway activation modeling for two cohorts of Parkinson’s patients that underwent DBS surgery at two independent centers. Based on modeling results for a training cohort, we proposed a method to derive an optimal activation profile using differences in activation rates and UPDRS-III scores assessed for two protocols evaluated in each patient. Furthermore, we evaluated the importance of activation in particular pathways for the obtained profile. At last, we demonstrated that proximity to the optimal activation profile in the multidimensional pathway activation space has a statistically significant correlation with UPDRS-III improvement in a test cohort from an independent center. ### Effect of Activation in Individual Pathways Contrary to our prior expectations, the computed activation rates of individual axonal pathways demonstrated no clear correlation with clinical improvements. This result is especially remarkable for the hyperdirect pathway (HDP), associated with motor symptom alleviation in Parkinson’s disease (PD) [9, 39, 48]. Furthermore, even though activation in one branch of the HDP correlated with UPDRS-III improvement from baseline, the result was not reproducible in the test cohort. One possible explanation is that the correlation of improvement with HDP activation might be hindered by the nearby corticofugal pathway, whose stimulaton has been associated with tetanic motor contractions [51]. In addition, subthreshold stimulation of these tracts has been reported to exacerbate akinesia and bradykinesia [55], so patients may not necessarily show significant relief of motor symptoms despite high activation rates in the HDP. Admittedly, this was not the case for the best responder in the test cohort (see Fig. 3.), possibly due to its individual lower responsiveness to the stimulation. We also analyzed correlations of UPDRS-III improvement with activation in the HDP, the corticofugal, and the pallidothalamic pathways without differentiating between branches. However, this did not yield conclusive results. Although high activation in the corticofugal pathway induces motor contractions, such stimulation protocols are avoided in chronic DBS, so this side effect could not be quantified in the correlation tests. For the cerebellothalamic pathway, we did observe a correlation between activation rates and tremor suppression, which is in agreement with other studies [2, 4, 10]. However, the correlation was present only when comparing two DBS protocols, but not DBS against baseline (see Suppl. 1). This suggests involvement of other pathways in thesuppression of tremors [52]. Furthermore, our model did not assign a high weighting of the cerebellothalamic pathway for predicting optimal outcomes. Nevertheless, the correlation of interprotocol differences implies that more effective alleviation of tremor could be achieved by steering current to the cerebellothalamic pathway. ### Effect of Balanced Pathway Activation More conclusive results were obtained when considering motor improvement as a function of the whole activation profile. Defined by the weighted Canberra distance, discrepancy of test cohort profiles with the optimal activation rate profile demonstrated a statistically significant inverse correlation with UPRDS-III improvement from baseline. Importantly, the optimal profile was derived from datapoints of the training cohort, in particular, from differences of DBS-on intrapatient UPDRS-III scores and corresponding activation rates. When analyzing the optimal profile, it is notable that two branches of the HDP to the primary motor cortex and of the corticofugal pathway to the supplementary motor area were assigned high weighting factors (see Fig. 4, B). For the latter, the activation rate in the optimal profile was minimal, suggesting its detrimental effect. On the other hand, the moderate activation rates in the presumably beneficial HDP branches can be explained by the proximity to the corticofugal pathway (note the high activation rates of the HDP and the corticofugal pathway for the best responder in the test cohort.) The activation rates and the weighting factors for the optimal profile could be used to draw conclusions regarding possible mechanisms of action of STN-DBS. The most noteworthy is the presumed importance of the pallidothalamic projections: the ansa lenticularis and lenticular fasciculus, which have been considered a functional continuum [45, 46]. By comparing the weighted optimal profile to the best responders, we can hypothesize that these pathways as well as the HDP are beneficial for alleviation of motor symptoms in PD. The ansa lenticularis and lenticular fasciculus are the primary inhibitory outputs of the basal ganglia to the ventral anterior thalamic nucleus [37], which itself is projecting to motor-relevant cortical regions. In a PD affected network, the GPi exhibits bursting activity [50], which has been associated with generation of motor symptoms [35]. Alleviation could arise from a direct modulation of this pathological activity, e.g. via an “informational lesion” induced by the non-physiological pattern of DBS [21]. Furthermore, therapeutic effects of stimulation of the pallidothalamic pathways would also explain the comparable efficiency of GPi-DBS in treating motor symptoms of PD. In historical context of ablative surgery, lesioning of ansa lenticularis was shown to be beneficial for tremor alleviation [54]. Considering these pathways as potential targets has important clinical implications. For DBS patients who respond poorly to the treatment, such as those in whom the active electrode contacts are outside the dorsolateral “sweet spot” of the STN [11], clinicians could attempt to alleviate symptoms by thoroughly investigating stimulation responses for electrode contacts near the ansa lenticularis and lenticular fasciculus. Prospective studies could consider new implantation trajectories that allow both mechanisms with precise targeting of the hyperdirect and pallidothalamic pathways, preferably stimulated via different electrode contacts. This strategy would offer the possibility of inducing desynchronised DBS patterns. At the same time, one should be aware that the STN is a site of convergence of different functional circuits [1], which is the most probable reason for its efficiency as a DBS target. Apart from the obvious cases, such as the near-zero activation of the corticofugal branch to the supplementary motor area (Fig. 4), the weighted profile of optimal activation rates does not suggest which pathways are beneficial or detrimental, but rather indicates a theoretically optimal balance of DBS-induced activation. Such a metric is more feasible than separately defined optimal activation rates, especially since pathway-specific DBS is not yet possible with clinically approved electrodes. ### Limitations The most disputable result presented in this study is the high activation in the corticofugal pathway computed for the best responder in the test cohort. For clinical application, this activation profile would be undesirable due to possible side effects. At the same time, the activation of the corticofugal branch to the lower extremity of the primary motor cortex was nearly zero for this patient, as well as for others in both cohorts (see Fig. 3). Further studies are warranted to investigate pathway activation profiles associated with DBS-induced side-effects. Nevertheless, one should consider that the arguable results of correlation tests for individual pathways might originate from uncertainties introduced by the computational model. We suspect that the lack of data on pathway-specific axonal morphology is the primary source of the error in pathway activation modeling. Furthermore, additional inaccuracies arise from volume conductor modeling, processing of medical images and clinical evaluation. Investigating present results by means of uncertainty quantification, as was shown in [6, 49] in the context of DBS, could enhance our understanding of how results should be interpreted. In the optimal profile, noteworthy is the small weighting factor of the sensorimotor pallidosubthalamic and reverse projections. That might be attributed to the limitations of the pathway activation modeling when investigating this local circuit, which excitatory-inhibitory reciprocity is a probable source of synchronization in the basal ganglia [44]. These short projections follow nearly the same trajectory, hence nearly the same DBS-induced activation rates, and a neural activity of one pathway is directly modulated not only by the DBS-induced membrane polarization, but simultaneously via the activity of the reciprocal pathway, which is also recruited by DBS. Therefore, a simulation model that could more accurately represent the circuit would need to account for temporal dynamics and consider effects of extracellular stimulation on synaptic inputs [12]. Besides, due to the electrode implantation, these short pathways are relatively more exposed to glial scaring that potentially affects neurotransmission and extracellular ionic concentrations [33]. Apart from the limitations of the computational model mentioned above, the following study limitations deserve mentioning. First, not all pathways affected by STN-DBS were modeled. To decrease the computational effort, we deliberately excluded the anterior cingulate and the prefrontal cortex assuming that they play a minor role in the motor symptoms of Parkinson’s disease, although connectivity with the latter was associated with rigidity improvement [3]. Among the local projections not examined, pathways of the substantia nigra pars reticulata are of particular interest, as it plays a similar role in the basal ganglia as the GPi. Unfortunately, these projections were not included in the basal ganglia pathway atlas [47], employed here as an exclusive source of the structural connectivity. Developed by experienced neuroanatomists using advanced visualization techniques, this atlas is less prone to contain false-positive trajectories occurring due to a poor signal-to-noise ratio of diffusion imaging, commonly employed for individualized fiber tractography. However, the atlas is neither patient-specific nor able to account for disease-related changes. In the present study, individual anatomical variability was partially accounted for by translation of pathways to the patient-specific space using a subcortical normalization strategy [14]. It should be noted that the optimal profile was derived based on UPDRS-III improvement in the training cohort, while validated on the test cohort assessed with MDS-UPDRS-III score. However, the scores are strongly correlated (*r >* 0.95) [42] and were shown to be linearly dependent [19, 25], hence invariant in the Pearson correlation test. Furthermore, to generalize UPDRS-III improvement, we avoided differentiating between lateral and axial symptoms when analyzing activation rates, which were calculated separately for each hemisphere and then averaged. Differences in the pathway activation rates across both hemispheres are presented in Suppl. 3. Finally, the optimization procedure of the weighting factors, conducted in the training cohort, was developed rather intuitively, and future studies should consider more robust methods. It also remains unclear whether optimization should emphasize the interprotocol correlation of the improvement with the distance to the optimal profile or a general increase of correlation for all datapoints. For the given datapoints, we did not observe significantly higher predictability in the test cohort for either case. ## Supporting information Suppl. material [[supplements/264743_file02.pdf]](pending:yes) ## Data Availability The data that support the findings of this study are available from the first corresponding author. ## 1 Funding Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – SFB 1270/1 – 299150580 and TRR 295 – 424778381. A.H. was further supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Emmy Noether Stipend – 410169619 as well as Deutsches Zentrum für Luftund Raumfahrt (DynaSti grant within the EU Joint Programme Neurodegenerative Disease Research, JPND). A.H. is participant in the BIH-Charité Clinician Scientist Program funded by the Charité –Universitätsmedizin Berlin and the Berlin Institute of Health. ## 2 Competing interests J.V. have business relationships with Medtronic, Abbott, and Boston Scientific, which are makers of DBS devices, outside the submitted work. A.A.K. reports personal fees and non-financial support from Medtronic, personal fees from Boston Scientific, grants and personal fees from Abbott, outside the submitted work. M.R. reports grant support and honoraria for speaking from Medtronic and Boston Scientific, outside the submitted work. A.H. reports lecture fees for Medtronic and Boston Scientific and is a consultant for AlphaOmega. G.R.W received travel expenses and attendance fees by Boston Scientific and Ipsen Pharma. K.B., N.L., C.N., J.R., H.E. and U.v.R. have nothing to disclose. ## Footnotes * * konstantin.butenko{at}charite.de, ursula.van-rienen{at}uni-rostock.de * Received October 8, 2021. * Revision received October 8, 2021. * Accepted October 11, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1. E. A. Accolla, M. Herrojo Ruiz, A. Horn, G.-H. Schneider, T. Schmitz-Hübsch, B. Draganski, and A. A. Kühn. Brain networks modulated by subthalamic nucleus deep brain stimulation. Brain, 139(9):2503–2515, 07 2016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/aww182&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27412387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 2. 2. H. Akram, V. Dayal, P. Mahlknecht, D. Georgiev, J. Hyam, T. Foltynie, and et al. Connectivity derived thalamic segmentation in deep brain stimulation for tremor. NeuroImage: Clinical, 18:130–142, 2018. 3. 3. H. Akram, S. N. Sotiropoulos, S. Jbabdi, D. Georgiev, P. Mahlknecht, J. Hyam, T. Foltynie, P. Limousin, E. De Vita, M. Jahanshahi, M. Hariz, J. Ashburner, T. Behrens, and L. Zrinzo. Subthalamic deep brain stimulation sweet spots and hyperdirect cortical connectivity in Parkinson’s disease. NeuroImage, 158:332–345, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2017.07.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28711737&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 4. 4. B. Al-Fatly, S. Ewert, D. Kübler, D. Kroneberg, A. Horn, and A. A. Kühn. Connectivity profile of thalamic deep brain stimulation to effectively treat essential tremor. Brain, 142(10):3086–3098, 08 2019. 5. 5. B. Avants, C. Epstein, M. Grossman, and J. Gee. Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain. Medical Image Analysis, 12(1):26–41, 2008. Special Issue on The Third International Workshop on Biomedical Image Registration – WBIR 2006. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.media.2007.06.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17659998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254032800004&link_type=ISI) 6. 6. K. Butenko, C. Bahls, and U. v. Rienen. Evaluation of epistemic uncertainties for bipolar deep brain stimulation in rodent models. In 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), pages 2136–2140, 2019. 7. 7. K. Butenko, C. Bahls, M. Schröder, R. Köhling, and U. van Rienen. OSS-DBS: Open-source simulation platform for deep brain stimulation with a comprehensive automated modeling. PLoS Computational Biology, 16(7):1–18, 2020. 8. 8. C. Butson and C. C. McIntyre. Tissue and electrode capacitance reduce neural activation volumes during deep brain stimulation. Clinical Neurophysiology, 116(10):2490–2500, 2005. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2005.06.023&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16125463&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000232233300023&link_type=ISI) 9. 9. Y. Chen, S. Ge, Y. Li, N. Li, J. Wang, X. Wang, and et al. Role of the cortico-subthalamic hyperdirect pathway in deep brain stimulation for the treatment of Parkinson disease: A diffusion tensor imaging study. World Neurosurgery, 114:e1079–e1085, 2018. 10. 10. V. A. Coenen, B. Sajonz, T. Prokop, M. Reisert, T. Piroth, and et al. The dentato-rubro-thalamic tract as the potential common deep brain stimulation target for tremor of various origin: an observational case series. Acta Neurochirurgica, 162:1053–1066, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00701-020-04248-2&link_type=DOI) 11. 11. T. A. Dembek, J. Roediger, A. Horn, P. Reker, C. Oehrn, H. S. Dafsari, N. Li Kühn, G. R. Fink, V. Visser-Vandewalle, M. T. Barbe, and L. Timmermann. Probabilistic sweet spots predict motor outcome for deep brain stimulation in Parkinson disease. Annals of Neurology, 86(4):527–538, 2019. 12. 12. J. O. Dostrovsky, R. Levy, J. P. Wu, W. D. Hutchison, R. R. Tasker, and A. M. Lozano. Microstimulation-induced inhibition of neuronal firing in human globus pallidus. Journal of Neurophysiology, 84(1):570–574, 2000. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/jn.2000.84.1.570&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10899228&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000088215000051&link_type=ISI) 13. 13. G. Duffley, D. Anderson, J. Vorwerk, A. Dorval, and C. Butson. Evaluation of methodologies for computing the deep brain stimulation volume of tissue activated. Journal of Neural Engineering, 16(6):1–15, 2019. 14. 14. S. Ewert, A. Horn, F. Finkel, N. Li, A. A. Kühn, and T. M. Herrington. Optimization and comparative evaluation of nonlinear deformation algorithms for atlas-based segmentation of DBS target nuclei. NeuroImage, 184:586–598, 2019. 15. 15. S. Ewert, P. Plettig, N. Li, M. Chakravarty, D. Collins, T. Herrington, and et al. Toward defining deep brain stimulation targets in MNI space: A subcortical atlas based on multimodal MRI, histology and structural connectivity. NeuroImage, 170:271–282, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2017.05.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28536045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 16. 16. K. J. Friston, A. P. Holmes, K. J. Worsley, J.-P. Poline, C. D. Frith, and R. S. J. Frackowiak. Statistical parametric maps in functional imaging: A general linear approach. Human Brain Mapping, 2(4):189–210, 1994. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.460020402&link_type=DOI) 17. 17. S. Gabriel, R. Lau, and C. Gabriel. The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues. Physics in Medicine and Biology, 41(11):2271–2293, nov 1996. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/0031-9155/41/11/003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8938026&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VR53200003&link_type=ISI) 18. 18. L. Geddes and L. Baker. The specific resistance of biological material—A compendium of data for the biomedical engineer and physiologist. Medical and biological engineering, 5:271–293, 1967. 19. 19. C. G. Goetz, G. T. Stebbins, and B. C. Tilley. Calibration of unified Parkinson’s disease rating scale scores to movement disorder society-unified Parkinson’s disease rating scale scores. Movement Disorders, 27(10):1239–1242, 2012. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mds.25122&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22886777&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 20. 20. M. Goftari, J. Kim, E. Johnson, R. Patriat, T. Palnitkar, N. Harel, and et al. Pallidothalamic tract activation predicts suppression of stimulation-induced dyskinesias in a case study of Parkinson’s disease. Brain Stimulation, 13(6):1821–1823, 2020. 21. 21. W. M. Grill, A. Snyder, and S. Miocinovic. Deep brain stimulation creates an informational lesion of the stimulated nucleus. Neuroreport, 15:1137–1140, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00001756-200405190-00011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15129161&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225140400011&link_type=ISI) 22. 22. D. Güllmar, J. Haueisen, and J. Reichenbach. Influence of anisotropic electrical conductivity in white matter tissue on the EEG/MEG forward and inverse solution. A high-resolution whole head simulation study. NeuroImage, 51(1):145–163, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2010.02.014&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20156576&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 23. 23. K. Gunalan, A. Chaturvedi, B. Howell, Y. Duchin, S. Lempka, R. Patriat, G. Sapiro, N. Harel, and C. C. McIntyre. Creating and parameterizing patient-specific deep brain stimulation pathway-activation models using the hyperdirect pathway as an example. PLoS ONE, 12(4):1–19, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0179159&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28957441&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 24. 24. C. J. Hartmann, J. L. Lujan, A. Chaturvedi, W. K. Goodman, M. S. Okun, C. C. McIntyre, and I. U. Haq. Tractography activation patterns in dorsolateral prefrontal cortex suggest better clinical responses in OCD DBS. Frontiers in Neuroscience, 9:519, 2016. 25. 25. J. G. Hentz, S. H. Mehta, H. A. Shill, E. Driver-Dunckley, T. G. Beach, and C. H. Adler. Simplified conversion method for unified Parkinson’s disease rating scale motor examinations. Movement Disorders, 30(14):1967–1970, 2015. 26. 26. A. Horn and A. A. Kühn. Lead-DBS: A toolbox for deep brain stimulation electrode localizations and visualizations. NeuroImage, 107:127–135, 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2014.12.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25498389&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 27. 27. A. Horn, N. Li, T. Dembek, A. Kappel, C. Boulay, S. Ewert, and et al. Lead-DBS v2: Towards a comprehensive pipeline for deep brain stimulation imaging. NeuroImage, 184:293–316, 2019. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 28. 28. A. Horn, M. Reich, J. Vorwerk, N. Li, G. Wenzel, Q. Fang, T. Schmitz-Hübsch, R. Nickl, A. Kupsch, J. Volkmann, A. A. Kühn, and M. D. Fox. Connectivity predicts deep brain stimulation outcome in Parkinson disease. Annals of Neurology, 82(1):67–78, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.24974&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28586141&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 29. 29. A. Horn, G. Wenzel, F. Irmen, J. Huebl, N. Li, W.-J. Neumann, P. Krause, G. Bohner, M. Scheel, and A. A. Kühn. Deep brain stimulation induced normalization of the human functional connectome in Parkinson’s disease. Brain, 142(10):3129–3143, 08 2019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awz239&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31412106&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 30. 30. B. Howell, K. S. Choi, K. Gunalan, J. Rajendra, H. S. Mayberg, and C. C. McIntyre. Quantifying the axonal pathways directly stimulated in therapeutic subcallosal cingulate deep brain stimulation. Human Brain Mapping, 40(3):889–903, 2019. 31. 31. B. Howell, F. Isbaine, J. T. Willie, E. Opri, R. E. Gross, C. De Hemptinne, P. A. Starr, C. C. McIntyre, and S. Miocinovic. Image-based biophysical modeling predicts cortical potentials evoked with subthalamic deep brain stimulation. Brain Stimulation, 14(3):549–563, 2021. 32. 32. A. Husch, M. V. Petersen, P. Gemmar, J. Goncalves, and F. Hertel. PaCER - a fully automated method for electrode trajectory and contact reconstruction in deep brain stimulation. NeuroImage: Clinical, 17:80–89, 2018. 33. 33. M. Jakobs, A. Fomenko, A. M. Lozano, and K. L. Kiening. Cellular, molecular, and clinical mechanisms of action of deep brain stimulation—a systematic review on established indications and outlook on future developments. EMBO Molecular Medicine, 11(4):e9575, 2019. 34. 34. K. Johnson, G. Duffley, T. Foltynie, M. Hariz, L. Zrinzo, E. Joyce, and et al. Basal ganglia pathways associated with therapeutic pallidal deep brain stimulation for Tourette syndrome. Biological Psychiatry: Cognitive Neuroscience and Neuroimaging, 2020. 35. 35. J. Kim, Y. Kim, R. Nakajima, A. Shin, M. Jeong, A. H. Park, and et al. Inhibitory basal ganglia inputs induce excitatory motor signals in the thalamus. Neuron, 95(5):1181–1196.e8, 2017. 36. 36. G. Lance and W. T. Williams. Mixed-data classificatory programs i - agglomerative systems. Aust. Comput. J., 1:15–20, 1967. 37. 37. J. L. Lanciego, N. Luquin, and J. A. Obeso. Functional neuroanatomy of the basal ganglia. Cold Spring Harbor Perspectives in Medicine, 2(12), 2012. 38. 38. N. Li, J. Baldermann, A. Kibleur, S. Treu, H. Akram, G. Elias, and et al. A unified connectomic target for deep brain stimulation in obsessive-compulsive disorder. Nature Communications, 11, 2020. 39. 39. Q. Li, Z.-M. Qian, G. W. Arbuthnott, Y. Ke, and W.-H. Yung. Cortical Effects of Deep Brain Stimulation: Implications for Pathogenesis and Treatment of Parkinson Disease. JAMA Neurology, 71(1):100–103, 01 2014. 40. 40.Liewald, R. Miller, N. Logothetis, H.-J. Wagner, and A. Schüz. Distribution of axon diameters in cortical white matter: an electron-microscopic study on three human brains and a macaque. Biological Cybernetics, 108:541–557, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00422-014-0626-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25142940&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 41. 41. C. McIntyre, A. G. Richardson, and W. M. Grill. Modeling the excitability of mammalian nerve fibers: Influence of afterpotentials on the recovery cycle. Journal of Neurophysiology, 87(2):995–1006, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/jn.00353.2001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11826063&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174735100031&link_type=ISI) 42. 42. M. Merello, E. R. Gerschcovich, D. Ballesteros, and D. Cerquetti. Correlation between the movement disorders society unified Parkinson’s disease rating scale (MDS-UPDRS) and the unified Parkinson’s disease rating scale (UPDRS) during l-dopa acute challenge. Parkinsonism & Related Disorders, 17(9):705–707, 2011. 43. 43. E.H. Middlebrooks, I. Tuna, S. S. Grewal, L. Almeida, M. Heckman, E. Lesser, K. Foote, M. S. Okun, and V. Holanda. Segmentation of the globus pallidus internus using probabilistic diffusion tractography for deep brain stimulation targeting in Parkinson disease. American Journal of Neuroradiology, 39(6):1127–1134, 2018. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo5OiIzOS82LzExMjciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMC8xMS8yMDIxLjEwLjA4LjIxMjY0NzQzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 44. 44. A. Nambu, Y. Tachibana, and S. Chiken. Cause of parkinsonian symptoms: Firing rate, firing pattern or dynamic activity changes? Basal Ganglia, 5(1):1–6, 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.baga.2014.11.001&link_type=DOI) 45. 45. C. Neudorfer and M. Maarouf. Neuroanatomical background and functional considerations for stereotactic interventions in the H fields of Forel. Brain Structure and Function, 223:17–30, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00429-017-1570-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29152666&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 46. 46. M. Parent and A. Parent. The pallidofugal motor fiber system in primates Parkinsonism & Related Disorders, 10(4):203–211, 2004. 47. 47. M. V. Petersen, J. Mlakar, S. Haber, M. Parent, Y. Smith, P. Strick, M. Griswold, and C. C. McIntyre. Holographic reconstruction of axonal pathways in the human brain. Neuron, 104(6):1056–1064.e3, 2019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2019.09.030&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 48. 48. T. H. Sanders and D. Jaeger. Optogenetic stimulation of cortico-subthalamic projections is sufficient to ameliorate bradykinesia in 6-OHDA lesioned mice. Neurobiology of Disease, 95:225–237, 2016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nbd.2016.07.021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27452483&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 49. 49. C. Schmidt, P. Grant, M. Lowery, and U. van Rienen. Influence of uncertainties in the material properties of brain tissue on the probabilistic volume of tissue activated. IEEE Transactions on Biomedical Engineering, 60(5):1378–1387, 2013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi:10.1109/TBME.2012.2235835&link_type=DOI) 50. 50. Y. Tachibana, H. Iwamuro, H. Kita, M. Takada, and A. Nambu. Subthalamo-pallidal interactions underlying parkinsonian neuronal oscillations in the primate basal ganglia. European Journal of Neuroscience, 34(9):1470–1484, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1460-9568.2011.07865.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22034978&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 51. 51.Tommasi, P. Krack, V. Fraix, and et al. Pyramidal tract side effects induced by deep brain stimulation of the subthalamic nucleus. Journal of Neurology, Neurosurgery & Psychiatry, 79(7):813–819, 2008. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czo4OiI3OS83LzgxMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEwLzExLzIwMjEuMTAuMDguMjEyNjQ3NDMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 52. 52. T. Tsuboi, J. K. Wong, R. S. Eisinger, L. Okromelidze, M. R. Burns Ramirez-Zamora, L. Almeida, A. Wagle Shukla, K. D. Foote, M. S. Okun, S. S. Grewal, and E. H. Middlebrooks. Comparative connectivity correlates of dystonic and essential tremor deep brain stimulation. Brain, 144(6):1774–1786, 04 2021. 53. 53.Wenzel, J. Roediger, C. Brücke, A. Marcelino, E. Gülke, M. Pötter-Nerger, and et al. CLOVER-DBS: Algorithm-guided deep brain stimulation-programming based on external sensor feedback evaluated in a prospective, randomized, crossover, double-blind, two-center study. Journal of Parkinson’s Disease, 2021. 54. 54. H. Wycis and E. Spiegel. Ansotomy in paralysis agitans. Confinia neurologica, 12(4):245–246, 1952. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000105783&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=13010002&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 55. 55. W. Xu, S. Miocinovic, J. Zhang, K. B. Baker, C. C. McIntyre, and J. L. Vitek. Dissociation of motor symptoms during deep brain stimulation of the subthalamic nucleus in the region of the internal capsule. Experimental Neurology, 228(2):294–297, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.expneurol.2010.08.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20713049&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F11%2F2021.10.08.21264743.atom) 56. 56. S. Zhang and K. Arfanakis. Evaluation of standardized and study-specific diffusion tensor imaging templates of the adult human brain: Template characteristics, spatial normalization accuracy, and detection of small inter-group FA differences. NeuroImage, 172:40–50, 2018. 57. 57.Zimmermann and U. van Rienen. Ambiguity in the interpretation of the low-frequency dielectric properties of biological tissues. Bioelectrochemistry, 140:107773, 2021. [1]: /embed/graphic-4.gif [2]: /embed/graphic-5.gif [3]: /embed/graphic-6.gif