Subcallosal cingulate deep brain stimulation evokes two distinct cortical responses via differential white matter activation.

Subcallosal cingulate (SCC) deep brain stimulation (DBS) is an emerging therapy for refractory depression. Good clinical outcomes are associated with the activation of white matter adjacent to the SCC. This activation produces a signature cortical evoked potential (EP), but it is unclear which of the many pathways in the vicinity of SCC is responsible for driving this response. Individualized biophysical models were built to achieve selective engagement of two target bundles: either the forceps minor (FM) or cingulum bundle (CB). Unilateral 2 Hz stimulation was performed in seven patients with treatment-resistant depression who responded to SCC DBS, and EPs were recorded using 256-sensor scalp electroencephalography. Two distinct EPs were observed: a 120 ms symmetric response spanning both hemispheres and a 60 ms asymmetrical EP. Activation of FM correlated with the symmetrical EPs, while activation of CB was correlated with the asymmetrical EPs. These results support prior model predictions that these two pathways are predominantly activated by clinical SCC DBS and provide first evidence of a link between cortical EPs and selective fiber bundle activation.

Optimization of an anatomically and electrically detailed rodent subthalamic nucleus neuron model.

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is an effective treatment for Parkinson's disease, but its mechanisms of action remain unclear. Detailed multi-compartment computational models of STN neurons are often used to study how DBS electric fields modulate the neurons. However, currently available STN neuron models have some limitations in their biophysical realism. In turn, the goal of this study was to update a detailed rodent STN neuron model originally developed by Gillies & Willshaw [2006]. Our design requirements consisted of explicitly representing an axon connected to the neuron and updating the ion channel distributions based on the experimental literature to match established electrophysiological features of rodent STN neurons. We found that adding an axon to the STN neuron model substantially altered its firing characteristics. We then used a genetic algorithm to optimize the biophysical parameters of the model. The optimized model exhibited spontaneous firing, action potential shape, hyperpolarization response, and frequency-current curve that aligned well with experimental recordings from STN neurons. Subsequently, we evaluated the general compatibility of the updated biophysics by applying them to 26 different STN neuron morphologies derived from 3D anatomical reconstructions. The different morphologies affected the firing behavior of the model, but the updated biophysics were robustly capable of maintaining the desired electrophysiological features. The new STN neuron model developed in this work offers a valuable tool for studying STN neuron firing properties, and may find application in simulating STN local field potentials or analysis of the effects of STN DBS.

Coupled Activation of the Hyperdirect and Cerebellothalamic Pathways with Zona Incerta Deep Brain Stimulation.

BACKGROUND: Deep brain stimulation (DBS) of the subthalamic nucleus (STN) or ventral intermediate nucleus (VIM) are established targets for the treatment of Parkinson's disease (PD) or essential tremor (ET), respectively. However, DBS of the zona incerta (ZI) can be effective for both disorders. VIM DBS is assumed to achieve its therapeutic effect via activation of the cerebellothalamic (CBT) pathway, whereas the activation of the hyperdirect (HD) pathway likely plays a role in the mechanisms of STN DBS. Interestingly, HD pathway axons also emit collaterals to the ZI and red nucleus (RN) and the CBT pathway courses nearby to the ZI. OBJECTIVE: The aim was to examine the ability of ZI DBS to mutually activate the HD and CBT pathways in a detailed computational model of human DBS. METHODS: We extended a previous model of the human HD pathway to incorporate axon collaterals to the ZI and RN. The anatomical framework of the model system also included representations of the CBT pathway and internal capsule (IC) fibers of passage. We then performed detailed biophysical simulations to quantify DBS activation of the HD, CBT, and IC pathways with electrodes located in either the STN or ZI. RESULTS: STN DBS and ZI DBS both robustly activated the HD pathway. However, STN DBS was limited by IC activation at higher stimulus amplitudes. Alternatively, ZI DBS avoided IC activation while simultaneously activating the HD and CBT pathways. CONCLUSIONS: From both neuroanatomical and biophysical perspectives, ZI DBS represents an advantageous target for coupled activation of the HD and CBT pathways. © 2024 International Parkinson and Movement Disorder Society.

Does Vestibular Motion Perception Correlate with Axonal Pathways Stimulated by Subthalamic Deep Brain Stimulation in Parkinson's Disease?

Perception of our linear motion - heading - is critical for postural control, gait, and locomotion, and it is impaired in Parkinson's disease (PD). Deep brain stimulation (DBS) has variable effects on vestibular heading perception, depending on the location of the electrodes within the subthalamic nucleus (STN). Here, we aimed to find the anatomical correlates of heading perception in PD. Fourteen PD participants with bilateral STN DBS performed a two-alternative forced-choice discrimination task where a motion platform delivered translational forward movements with a heading angle varying between 0 and 30° to the left or to the right with respect to the straight-ahead direction. Using psychometric curves, we derived the heading discrimination threshold angle of each patient from the response data. We created patient-specific DBS models and calculated the percentages of stimulated axonal pathways that are anatomically adjacent to the STN and known to play a major role in vestibular information processing. We performed correlation analyses to investigate the extent of these white matter tracts' involvement in heading perception. Significant positive correlations were identified between improved heading discrimination for rightward heading and the percentage of activated streamlines of the contralateral hyperdirect, pallido-subthalamic, and subthalamo-pallidal pathways. The hyperdirect pathways are thought to provide top-down control over STN connections to the cerebellum. In addition, STN may also antidromically activate collaterals of hyperdirect pathway that projects to the precerebellar pontine nuclei. In select cases, there was strong activation of the cerebello-thalamic projections, but it was not consistently present in all participants. Large volumetric overlap between the volume of tissue activation and the STN in the left hemisphere positively impacted rightward heading perception. Altogether, the results suggest heavy involvement of basal ganglia cerebellar network in STN-induced modulation of vestibular heading perception in PD.

Climate control on terrestrial biospheric carbon turnover.

Terrestrial vegetation and soils hold three times more carbon than the atmosphere. Much debate concerns how anthropogenic activity will perturb these surface reservoirs, potentially exacerbating ongoing changes to the climate system. Uncertainties specifically persist in extrapolating point-source observations to ecosystem-scale budgets and fluxes, which require consideration of vertical and lateral processes on multiple temporal and spatial scales. To explore controls on organic carbon (OC) turnover at the river basin scale, we present radiocarbon (14C) ages on two groups of molecular tracers of plant-derived carbon-leaf-wax lipids and lignin phenols-from a globally distributed suite of rivers. We find significant negative relationships between the 14C age of these biomarkers and mean annual temperature and precipitation. Moreover, riverine biospheric-carbon ages scale proportionally with basin-wide soil carbon turnover times and soil 14C ages, implicating OC cycling within soils as a primary control on exported biomarker ages and revealing a broad distribution of soil OC reactivities. The ubiquitous occurrence of a long-lived soil OC pool suggests soil OC is globally vulnerable to perturbations by future temperature and precipitation increase. Scaling of riverine biospheric-carbon ages with soil OC turnover shows the former can constrain the sensitivity of carbon dynamics to environmental controls on broad spatial scales. Extracting this information from fluvially dominated sedimentary sequences may inform past variations in soil OC turnover in response to anthropogenic and/or climate perturbations. In turn, monitoring riverine OC composition may help detect future climate-change-induced perturbations of soil OC turnover and stocks.

Biophysical Principles and Computational Modeling of Deep Brain Stimulation.

BACKGROUND: Deep brain stimulation (DBS) has revolutionized the treatment of neurological disorders, yet the mechanisms of DBS are still under investigation. Computational models are important in silico tools for elucidating these underlying principles and potentially for personalizing DBS therapy to individual patients. The basic principles underlying neurostimulation computational models, however, are not well known in the clinical neuromodulation community. OBJECTIVE: In this study, we present a tutorial on the derivation of computational models of DBS and outline the biophysical contributions of electrodes, stimulation parameters, and tissue substrates to the effects of DBS. RESULTS: Given that many aspects of DBS are difficult to characterize experimentally, computational models have played an important role in understanding how material, size, shape, and contact segmentation influence device biocompatibility, energy efficiency, the spatial spread of the electric field, and the specificity of neural activation. Neural activation is dictated by stimulation parameters including frequency, current vs voltage control, amplitude, pulse width, polarity configurations, and waveform. These parameters also affect the potential for tissue damage, energy efficiency, the spatial spread of the electric field, and the specificity of neural activation. Activation of the neural substrate also is influenced by the encapsulation layer surrounding the electrode, the conductivity of the surrounding tissue, and the size and orientation of white matter fibers. These properties modulate the effects of the electric field and determine the ultimate therapeutic response. CONCLUSION: This article describes biophysical principles that are useful for understanding the mechanisms of neurostimulation.

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.

Automatic Segmentation of Parkinson Disease Therapeutic Targets Using Nonlinear Registration and Clinical MR Imaging: Comparison of Methodology, Presence of Disease, and Quality Control.

INTRODUCTION: Accurate and precise delineation of the globus pallidus pars interna (GPi) and subthalamic nucleus (STN) is critical for the clinical treatment and research of Parkinson's disease (PD). Automated segmentation is a developing technology which addresses limitations of visualizing deep nuclei on MR imaging and standardizing their definition in research applications. We sought to compare manual segmentation with three workflows for template-to-patient nonlinear registration providing atlas-based automatic segmentation of deep nuclei. METHODS: Bilateral GPi, STN, and red nucleus (RN) were segmented for 20 PD and 20 healthy control (HC) subjects using 3T MRIs acquired for clinical purposes. The automated workflows used were an option available in clinical practice and two common research protocols. Quality control (QC) was performed on registered templates via visual inspection of readily discernible brain structures. Manual segmentation using T1, proton density, and T2 sequences was used as "ground truth" data for comparison. Dice similarity coefficient (DSC) was used to assess agreement between segmented nuclei. Further analysis was done to compare the influences of disease state and QC classifications on DSC. RESULTS: Automated segmentation workflows (CIT-S, CRV-AB, and DIST-S) had the highest DSC for the RN and lowest for the STN. Manual segmentations outperformed automated segmentation for all workflows and nuclei; however, for 3/9 workflows (CIT-S STN, CRV-AB STN, and CRV-AB GPi) the differences were not statically significant. HC and PD only showed significant differences in 1/9 comparisons (DIST-S GPi). QC classification only demonstrated significantly higher DSC in 2/9 comparisons (CRV-AB RN and GPi). CONCLUSION: Manual segmentations generally performed better than automated segmentations. Disease state does not appear to have a significant effect on the quality of automated segmentations via nonlinear template-to-patient registration. Notably, visual inspection of template registration is a poor indicator of the accuracy of deep nuclei segmentation. As automatic segmentation methods continue to evolve, efficient and reliable QC methods will be necessary to support safe and effective integration into clinical workflows.

Subthalamic deep brain stimulation of an anatomically detailed model of the human hyperdirect pathway.

The motor hyperdirect pathway (HDP) is considered a key target in the treatment of Parkinson's disease with subthalamic deep brain stimulation (DBS). This hypothesis is partially derived from the association of HDP activation with evoked potentials (EPs) generated in the motor cortex and subthalamic nucleus (STN) after a DBS pulse. However, the biophysical details of how and when DBS-induced action potentials (APs) in HDP neurons reach their terminations in the cortex or STN remain unclear. Therefore, we used an anatomically detailed representation of the motor HDP, as well as the internal capsule (IC), in a model of human subthalamic DBS to explore AP activation and transmission in the HDP and IC. Our results show that small diameter HDP axons exhibited AP initiation in their subthalamic terminal arbor, which resulted in relatively long transmission latencies to cortex (∼3.5-8 ms). Alternatively, large diameter HDP axons were most likely to be directly activated in the capsular region, which resulted in short transmission times to the cortex (∼1-3 ms). However, those large diameter HDP antidromic APs would be indistinguishable from any other IC axons that were also activated by the stimulus. Conversely, DBS-induced APs in both small and large diameter HDP axons reached their synaptic boutons in the STN with similar timings, but both spanned a wide temporal range (∼0.5-5 ms). We also found that using anodic or bipolar stimulation helped to bias activation of the HDP over the IC. These computational results provide useful information for linking HDP activation with EP recordings in clinical experiments.NEW & NOTEWORTHY We used biophysical models to study pathway recruitment and conduction latencies of the hyperdirect pathway (HDP) in response to subthalamic deep brain stimulation (DBS). The model system allowed us to assess the influence of increased anatomical realism on pathway activity and the possibility of identifying HDP activity in evoked potentials (EPs) recorded in either the subthalamic nucleus (STN) or cortex. The model predicts that HDP activation is accentuated by complex axonal branching in the STN.

Temporal Patterns of Spontaneous Fixational Eye Movements: The Influence of Basal Ganglia.

BACKGROUND: Spontaneity is a unique feature of the nervous system. One of the fundamentally critical and recognized forms of spontaneous motor activity is witnessed in the visuomotor system. Microsaccades, the miniature spontaneous eye movements, are critical for the visual perception. We hypothesized that microsaccades follow specific temporal patterns that are modulated by the basal ganglia output. METHODS: We used high-resolution video-oculography to capture microsaccades in 48 subjects (31 healthy and 17 with Parkinson's disease) when subjects were asked to hold their gaze on a straight-ahead target projected on white background. We analyzed spontaneous discharge patterns of microsaccades. RESULTS: The first analysis considering coefficient of variation in intersaccadic interval distribution demonstrated that microsaccades in Parkinson's disease are more dispersed than the control group. The second analysis scrutinized microsaccades' temporal variability and revealed 3 distinct occurrence patterns: regular rhythmic, clustered, and randomly occurring following a Poisson-like process. The regular pattern was relatively more common in Parkinson's disease. Subthalamic DBS modulated this temporal pattern. The amount of change in the temporal variability depended on the DBS-induced volume of tissue activation and its overlap with the subthalamic nucleus. The third analysis determined the autocorrelations of microsaccades within 2-second time windows. We found that Parkinson's disease altered local temporal organization in microsaccade generation, and DBS had a modulatory effect. CONCLUSION: The microsaccades occur in 3 temporal patterns. The basal ganglia are one of the modulators of the microsaccade spontaneity.

Levodopa Versus Dopamine Agonist after Subthalamic Stimulation in Parkinson's Disease.

BACKGROUND: No clinical trials have been specifically designed to compare medical treatments after surgery in Parkinson's disease (PD). OBJECTIVE: Study's objective was to compare the efficacy and safety of levodopa versus dopamine agonist monotherapy after deep brain stimulation (DBS) in PD. METHODS: Thirty-five surgical candidates were randomly assigned to receive postoperative monotherapy with either levodopa or dopamine agonist in a randomized, single-blind study. All patients were reevaluated in short- (3 months), mid- (6 months), and long-term (2.5 years) follow-up after surgery. The primary outcome measure was the change in the Non-Motor Symptoms Scale (NMSS) 3 months after surgery. Secondary outcome measures were the percentage of patients maintaining monotherapy, change in motor symptoms, and specific non-motor symptoms (NMS). Analysis was performed primarily in the intention-to-treat population. RESULTS: Randomization did not significantly affect the primary outcome (difference in NMSS between treatment groups was 4.88 [95% confidence interval: -11.78-21.53, P = 0.566]). In short- and mid-term follow-up, monotherapy was safe and feasible in more than half of patients (60% in short- and 51.5% in mid-term follow-up), but it was more often possible for patients on levodopa. The ability to maintain dopamine agonist monotherapy was related to optimal contact location. In the long term, levodopa monotherapy was feasible only in a minority of patients (34.2%), whereas dopamine agonist monotherapy was not tolerated due to worsening of motor conditions or occurrence of impulse control disorders. CONCLUSIONS: This trial provides evidence for simplifying pharmacological treatment after functional neurosurgery for PD. The reduction in dopamine receptor agonists should be attempted while monitoring for occurrence of NMSs, such as apathy and sleep disturbances. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Effects of subthalamic deep brain stimulation on fixational eye movements in Parkinson's disease.

Miniature yoked eye movements, fixational saccades, are critical to counteract visual fading. Fixational saccades are followed by a return saccades forming squarewaves. Present in healthy states, squarewaves, if too many or too big, affect visual stability. Parkinson's disease (PD), where visual deficits are not uncommon, is associated with the squarewaves that are excessive in number or size. Our working hypothesis is that the basal ganglia are at the epicenter of the abnormal fixational saccades and squarewaves in PD; the effects are manifested through their connections to the superior colliculus (affecting saccade frequency and amplitude) and the cerebellum (affecting velocity and amplitude). We predict that the subthalamic deep brain stimulation (DBS) variably affects the amplitude, frequency, and velocity of fixational saccade and that the effect depends on the electrode's proximity or the volume of activated tissue in the subthalamic nucleus' connections with the superior colliculus or the cerebellum. We found that DBS modulated saccade amplitude, frequency, and velocity in 11 PD patients. Although all three parameters were affected, the extent of the effects varied amongst subjects. The modulation was dependent upon the location and size of the electrically activated volume of the subthalamic region.

Feasibility of Interferential and Pulsed Transcranial Electrical Stimulation for Neuromodulation at the Human Scale.

OBJECTIVES: Transcranial electrical stimulation (tES) is a promising tool for modulating neural activity, but tES has poor penetrability and spatiotemporal resolution compared to invasive techniques like deep brain stimulation (DBS). Interferential strategies for alternating-current stimulation (IF-tACS) and pulsed/intersectional strategies for transcranial direct-current stimulation (IS-tDCS) address some of the limitations of tES, but the comparative advantages and disadvantages of these new techniques is not well understood. This study's objective was to evaluate the suprathreshold and subthreshold membrane dynamics of neurons in response to IF-tACS and IS-tDCS. MATERIALS AND METHODS: We analyzed the biophysics of IF-tACS and IS-tDCS using a bioelectric field model of tES. Neural responses were quantified for suprathreshold generation of action potentials in axons and for subthreshold modulation of membrane dynamics in spiking pyramidal neurons. RESULTS: IF-tACS and IS-tDCS could not directly activate axons at or below 10 mA, but within this current range, these fields were able to modulate, albeit indirectly, spiking activity in the neuron model. IF-tACS facilitated phase synchronization similar to tACS, and IS-tDCS enhanced and suppressed spiking activity similar to tDCS; however, in either case, the modulatory effects of these fields were less potent than their standard counterparts at a matched field intensity. Moreover, neither IF-tACS nor IS-tDCS improved the spatial selectivity of their parent strategies. CONCLUSIONS: Enhancing the spatiotemporal precision and penetrability of tES with interferential and intersectional strategies is possible at the human scale. However, IF-tACS or IS-tDCS will likely require spatial multiplexing with multiple simultaneous sources to counteract their reduced potency, compared to standard techniques, to maintain stimulation currents at tolerable levels.

StimVision v2: Examples and Applications in Subthalamic Deep Brain Stimulation for Parkinson's Disease.

OBJECTIVE: Subthalamic deep brain stimulation (DBS) is an established therapy for Parkinson's disease. Connectomic DBS modeling is a burgeoning subfield of research aimed at characterizing the axonal connections activated by DBS. This article describes our approach and methods for evolving the StimVision software platform to meet the technical demands of connectomic DBS modeling in the subthalamic region. MATERIALS AND METHODS: StimVision v2 was developed with Visualization Toolkit (VTK) libraries and integrates four major components: 1) medical image visualization, 2) axonal pathway visualization, 3) electrode positioning, and 4) stimulation calculation. RESULTS: StimVision v2 implemented two key technological advances for connectomic DBS analyses in the subthalamic region. First was the application of anatomical axonal pathway models to patient-specific DBS models. Second was the application of a novel driving-force method to estimate the response of those axonal pathways to DBS. Example simulations with directional DBS electrodes and clinically defined therapeutic DBS settings are presented to demonstrate the general outputs of StimVision v2 models. CONCLUSIONS: StimVision v2 provides the opportunity to evaluate patient-specific axonal pathway activation from subthalamic DBS using anatomically detailed pathway models and electrically detailed electric field distributions with interactive adjustment of the DBS electrode position and stimulation parameter settings.

Quantifying the axonal pathways directly stimulated in therapeutic subcallosal cingulate deep brain stimulation.

Deep brain stimulation (DBS) of the subcallosal cingulate (SCC) is an emerging experimental therapy for treatment-resistant depression. New developments in SCC DBS surgical targeting are focused on identifying specific axonal pathways for stimulation that are estimated from patient-specific computational models. This connectomic-based biophysical modeling strategy has proven successful in improving the clinical response to SCC DBS therapy, but the DBS models used to date have been relatively simplistic, limiting the precision of the pathway activation estimates. Therefore, we used the most detailed patient-specific foundation for DBS modeling currently available (i.e., field-cable modeling) to evaluate SCC DBS in our most recent cohort of six subjects, all of which were responders to the therapy. We quantified activation of four major pathways in the SCC region: forceps minor (FM), cingulum bundle (CB), uncinate fasciculus (UF), and subcortical connections between the frontal pole and the thalamus or ventral striatum (FP). We then used the percentage of activated axons in each pathway as regressors in a linear model to predict the time it took patients to reach a stable response, or TSR. Our analysis suggests that stimulation of the left and right CBs, as well as FM are the most likely therapeutic targets for SCC DBS. In addition, the right CB alone predicted 84% of the variation in the TSR, and the correlation was positive, suggesting that activation of the right CB beyond a critical percentage may actually protract the recovery process.

A Driving-Force Predictor for Estimating Pathway Activation in Patient-Specific Models of Deep Brain Stimulation.

OBJECTIVE: Detailed biophysical modeling of deep brain stimulation (DBS) provides a theoretical approach to quantify the cellular response to the applied electric field. However, the most accurate models for performing such analyses, patient-specific field-cable (FC) pathway-activation models (PAMs), are so technically demanding to implement that their use in clinical research is greatly limited. Predictive algorithms can simplify PAM calculations, but they generally fail to reproduce the output of FC models when evaluated over a wide range of clinically relevant stimulation parameters. Therefore, we set out to develop a novel driving-force (DF) predictive algorithm (DF-Howell), customized to the study of DBS, which can better match FC results. METHODS: We developed the DF-Howell algorithm and compared its predictions to FC PAM results, as well as to the DF-Peterson algorithm, which is currently the most accurate and generalizable DF-based method. Comparison of the various methods was quantified within the context of subthalamic DBS using activation thresholds of axons representing the internal capsule, hyperdirect pathway, and cerebellothalamic tract for various combinations of fiber diameters, stimulus pulse widths, and electrode configurations. RESULTS: The DF-Howell predictor estimated activation of the three axonal pathways with less than a 6.2% mean error with respect to the FC PAM for all 21 cases tested. In 15 of the 21 cases, DF-Howell outperformed DF-Peterson in estimating pathway activation, reducing mean-errors up to 22.5%. CONCLUSIONS: DF-Howell represents an accurate predictor for estimating axonal pathway activation in patient-specific DBS models, but errors still exist relative to FC PAM calculations. Nonetheless, the tractability of DF algorithms helps to reduce the technical barriers for performing accurate biophysical modeling in clinical DBS research studies.

Quantifying axonal responses in patient-specific models of subthalamic deep brain stimulation.

Medical imaging has played a major role in defining the general anatomical targets for deep brain stimulation (DBS) therapies. However, specifics on the underlying brain circuitry that is directly modulated by DBS electric fields remain relatively undefined. Detailed biophysical modeling of DBS provides an approach to quantify the theoretical responses to stimulation at the cellular level, and has established a key role for axonal activation in the therapeutic mechanisms of DBS. Estimates of DBS-induced axonal activation can then be coupled with advances in defining the structural connectome of the human brain to provide insight into the modulated brain circuitry and possible correlations with clinical outcomes. These pathway-activation models (PAMs) represent powerful tools for DBS research, but the theoretical predictions are highly dependent upon the underlying assumptions of the particular modeling strategy used to create the PAM. In general, three types of PAMs are used to estimate activation: 1) field-cable (FC) models, 2) driving force (DF) models, and 3) volume of tissue activated (VTA) models. FC models represent the "gold standard" for analysis but at the cost of extreme technical demands and computational resources. Consequently, DF and VTA PAMs, derived from simplified FC models, are typically used in clinical research studies, but the relative accuracy of these implementations is unknown. Therefore, we performed a head-to-head comparison of the different PAMs, specifically evaluating DBS of three different axonal pathways in the subthalamic region. The DF PAM was markedly more accurate than the VTA PAMs, but none of these simplified models were able to match the results of the patient-specific FC PAM across all pathways and combinations of stimulus parameters. These results highlight the limitations of using simplified predictors to estimate axonal stimulation and emphasize the need for novel algorithms that are both biophysically realistic and computationally simple.

Biophysical basis of subthalamic local field potentials recorded from deep brain stimulation electrodes.

Clinical deep brain stimulation (DBS) technology is evolving to enable chronic recording of local field potentials (LFPs) that represent electrophysiological biomarkers of the underlying disease state. However, little is known about the biophysical basis of LFPs, or how the patient's unique brain anatomy and electrode placement impact the recordings. Therefore, we developed a patient-specific computational framework to analyze LFP recordings within a clinical DBS context. We selected a subject with Parkinson's disease implanted with a Medtronic Activa PC+S DBS system and reconstructed their subthalamic nucleus (STN) and DBS electrode location using medical imaging data. The patient-specific STN volume was populated with 235,280 multicompartment STN neuron models, providing a neuron density consistent with histological measurements. Each neuron received time-varying synaptic inputs and generated transmembrane currents that gave rise to the LFP signal recorded at DBS electrode contacts residing in a finite element volume conductor model. We then used the model to study the role of synchronous beta-band inputs to the STN neurons on the recorded power spectrum. Three bipolar pairs of simultaneous clinical LFP recordings were used in combination with an optimization algorithm to customize the neural activity parameters in the model to the patient. The optimized model predicted a 2.4-mm radius of beta-synchronous neurons located in the dorsolateral STN. These theoretical results enable biophysical dissection of the LFP signal at the cellular level with direct comparison to the clinical recordings, and the model system provides a scientific platform to help guide the design of DBS technology focused on the use of subthalamic beta activity in closed-loop algorithms. NEW & NOTEWORTHY The analysis of deep brain stimulation of local field potential (LFP) data is rapidly expanding from scientific curiosity to the basis for clinical biomarkers capable of improving the therapeutic efficacy of stimulation. With this growing clinical importance comes a growing need to understand the underlying electrophysiological fundamentals of the signals and the factors contributing to their modulation. Our model reconstructs the clinical LFP from first principles and highlights the importance of patient-specific factors in dictating the signals recorded.

StimVision Software: Examples and Applications in Subcallosal Cingulate Deep Brain Stimulation for Depression.

OBJECTIVE: Create a software tool to facilitate tractography-based deep brain stimulation (DBS) electrode targeting within the patient-specific stereotactic coordinate system used in the operating room. APPROACH: StimVision was developed with Visualization Toolkit libraries and integrates four major components: 1) medical image visualization, 2) tractography visualization, 3) DBS electrode positioning, and 4) DBS activation volume calculation with tractography intersection. RESULTS: Initial applications of StimVision are focused on the study of subcallosal cingulate (SCC) DBS for the treatment of depression. Retrospective modeling results on SCC DBS have suggested that direct stimulation of a specific collection of tractographic pathways are necessary for therapeutic benefit; thereby creating a tractography-based DBS surgical targeting hypotheses. StimVision is the tool we created to facilitate prospective clinical evaluation of that hypothesis. SIGNIFICANCE: Retrospective tractography-based analyses are common in DBS research; however, intraoperative software tools for interactive selection of a tractography-based DBS target are not readily available. StimVision provides an academic research tool to assist clinical implementation of new DBS targeting strategies and postoperative evaluation of targeting outcome.