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.

Role of the volume conductor on simulations of local field potential recordings from deep brain stimulation electrodes.

OBJECTIVE: Local field potential (LFP) recordings from deep brain stimulation (DBS) electrodes are commonly used in research analyses, and are beginning to be used in clinical practice. Computational models of DBS LFPs provide tools for investigating the biophysics and neural synchronization that underlie LFP signals. However, technical standards for DBS LFP model parameterization remain to be established. Therefore, the goal of this study was to evaluate the role of the volume conductor (VC) model complexity on simulated LFP signals in the subthalamic nucleus (STN). APPROACH: We created a detailed human head VC model that explicitly represented the inhomogeneity and anisotropy associated with 12 different tissue structures. This VC model represented our "gold standard" for technical detail and electrical realism. We then incrementally decreased the complexity of the VC model and quantified the impact on the simulated LFP recordings. Identical STN neural source activity was used when comparing the different VC model variants. Results Ignoring tissue anisotropy reduced the simulated LFP amplitude by ~12%, while eliminating soft tissue heterogeneity had a negligible effect on the recordings. Simplification of the VC model to consist of a single homogenous isotropic tissue medium with a conductivity of 0.215 S/m contributed an additional ~3% to the error. SIGNIFICANCE: Highly detailed VC models do generate different results than simplified VC models. However, with errors in the range of ~15%, the use of a well-parameterized simple VC model is likely to be acceptable in most contexts for DBS LFP modeling.

Mitigating Mismatch Compression in Differential Local Field Potentials.

Deep brain stimulation (DBS) devices capable of measuring differential local field potentials ( ∂ LFP) enable neural recordings alongside clinical therapy. Efforts to identify oscillatory correlates of various brain disorders, or disease readouts, are growing but must proceed carefully to ensure readouts are not distorted by brain environment. In this report we identified, characterized, and mitigated a major source of distortion in ∂ LFP that we introduce as mismatch compression (MC). Using in vivo, in silico, and in vitro models of MC, we showed that impedance mismatches in the two recording electrodes can yield incomplete rejection of stimulation artifact and subsequent gain compression that distorts oscillatory power. We then developed and validated an opensource mitigation pipeline that mitigates the distortions arising from MC. This work enables more reliable oscillatory readouts for adaptive DBS applications.

Time-frequency signatures evoked by single-pulse deep brain stimulation to the subcallosal cingulate.

Precision targeting of specific white matter bundles that traverse the subcallosal cingulate (SCC) has been linked to efficacy of deep brain stimulation (DBS) for treatment resistant depression (TRD). Methods to confirm optimal target engagement in this heterogenous region are now critical to establish an objective treatment protocol. As yet unexamined are the time-frequency features of the SCC evoked potential (SCC-EP), including spectral power and phase-clustering. We examined these spectral features-evoked power and phase clustering-in a sample of TRD patients (n = 8) with implanted SCC stimulators. Electroencephalogram (EEG) was recorded during wakeful rest. Location of electrical stimulation in the SCC target region was the experimental manipulation. EEG was analyzed at the surface level with an average reference for a cluster of frontal sensors and at a time window identified by prior study (50-150 ms). Morlet wavelets generated indices of evoked power and inter-trial phase clustering. Enhanced phase clustering at theta frequency (4-7 Hz) was observed in every subject and was significantly correlated with SCC-EP magnitude, but only during left SCC stimulation. Stimulation to dorsal SCC evinced stronger phase clustering than ventral SCC. There was a weak correlation between phase clustering and white matter density. An increase in evoked delta power (2-4 Hz) was also coincident with SCC-EP, but was less consistent across participants. DBS evoked time-frequency features index mm-scale changes to the location of stimulation in the SCC target region and correlate with structural characteristics implicated in treatment optimization. Results also imply a shared generative mechanism (inter-trial phase clustering) between evoked potentials evinced by electrical stimulation and evoked potentials evinced by auditory/visual stimuli and behavioral tasks. Understanding how current injection impacts downstream cortical activity is essential to building new technologies that adapt treatment parameters to individual differences in neurophysiology.

Isolating two sources of variability of subcortical stimulation to quantify fluctuations of corticospinal tract excitability.

OBJECTIVE: Investigate the variability previously found with cortical stimulation and handheld transcranial magnetic stimulation (TMS) coils, criticized for its high potential of coil position fluctuations, bypassing the cortex using deep brain electrical stimulation (DBS) of the corticospinal tract with fixed electrodes where both latent variations of the coil position of TMS are eliminated and cortical excitation fluctuations should be absent. METHODS: Ten input-output curves were recorded from five anesthetized cats with implanted DBS electrodes targeting the corticospinal tract. Goodness of fit of regressions with a conventional single variability source as well as a dual variability source model was quantified using a Schwarz Bayesian Information approach to avoid overfitting. RESULTS: Motor evoked potentials (MEPs) through DBS of the corticospinal tract revealed short-term fluctuations in excitability of the targeted neuron pathway reflecting endogenous input-side variability at similar magnitude as TMS despite bypassing cortical networks. CONCLUSION: Input-side variability, i.e., variability resulting in changing MEP amplitudes as if the stimulation strength was modulated, also emerges in electrical stimulation at a similar degree and is not primarily a result of varying stimulation, such as minor coil movements in TMS. More importantly, this variability component is present, although the cortex is bypassed. Thus, it may be of spinal origin, which can include cortical input from spinal projections. Further, the nonlinearity of the compound variability entails complex heteroscedastic non-Gaussian distributions and typically does not allow simple linear averages in statistical analysis of MEPs. As the average is dominated by outliers, it risks bias. With appropriate regression, the net effects of excitatory and inhibitory inputs to the targeted neuron pathways become noninvasively observable and quantifiable. SIGNIFICANCE: The neural responses evoked by artificial stimulation in the cerebral cortex are variable. For example, MEPs in response to repeated presentations of the same stimulus can vary from no response to saturation across trials. Several sources of such variability have been suggested, and most of them may be technical in nature, but localization is missing.

Comparison of methodologies for modeling directional deep brain stimulation electrodes.

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

Image-based biophysical modeling predicts cortical potentials evoked with subthalamic deep brain stimulation.

BACKGROUND: Subthalamic deep brain stimulation (DBS) is an effective surgical treatment for Parkinson's disease and continues to advance technologically with an enormous parameter space. As such, in-silico DBS modeling systems have become common tools for research and development, but their underlying methods have yet to be standardized and validated. OBJECTIVE: Evaluate the accuracy of patient-specific estimates of neural pathway activations in the subthalamic region against intracranial, cortical evoked potential (EP) recordings. METHODS: Pathway activations were modeled in eleven patients using the latest advances in connectomic modeling of subthalamic DBS, focusing on the hyperdirect pathway (HDP) and corticospinal/bulbar tract (CSBT) for their relevance in human research studies. Correlations between pathway activations and respective EP amplitudes were quantified. RESULTS: Good model performance required accurate lead localization and image fusions, as well as appropriate selection of fiber diameter in the biophysical model. While optimal model parameters varied across patients, good performance could be achieved using a global set of parameters that explained 60% and 73% of electrophysiologic activations of CSBT and HDP, respectively. Moreover, restricted models fit to only EP amplitudes of eight standard (monopolar and bipolar) electrode configurations were able to extrapolate variation in EP amplitudes across other directional electrode configurations and stimulation parameters, with no significant reduction in model performance across the cohort. CONCLUSIONS: Our findings demonstrate that connectomic models of DBS with sufficient anatomical and electrical details can predict recruitment dynamics of white matter. These results will help to define connectomic modeling standards for preoperative surgical targeting and postoperative patient programming applications.

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.

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.

Field potential 1/f activity in the subcallosal cingulate region as a candidate signal for monitoring deep brain stimulation for treatment-resistant depression.

Subcallosal cingulate cortex deep brain stimulation (SCC-DBS) is an experimental therapy for treatment-resistant depression (TRD). Refinement and optimization of SCC-DBS will benefit from increased study of SCC electrophysiology in context of ongoing high-frequency SCC-DBS therapy. The study objective was a 7-mo observation of frequency-domain 1/f slope in off-stimulation local field potentials (SCC-LFPs) alongside standardized measurements of depression severity in 4 patients undergoing SCC-DBS. SCC was implanted bilaterally with a combined neurostimulation-LFP recording system. Following a 1-mo off-stimulation postoperative phase with multiple daily recordings, patients received bilateral SCC-DBS therapy (130 Hz, 90 µs) and weekly resting-state SCC-LFP recordings over a 6-mo treatment phase. 1/f slopes for each time point were estimated via linear regression of log-transformed Welch periodograms. General linear mixed-effects models were constructed to estimate pretreatment sources of 1/f slope variance, and 95% bootstrap confidence intervals were constructed to estimate treatment phase 1/f slope association with treatment response (50% decrease in preimplantation symptom severity). Results show the time of recording was a prominent source of pretreatment 1/f slope variance bilaterally, with increased 1/f slope magnitude observed during night hours (2300-0659). Increase in right 1/f slope was observed in the setting of treatment response, with bootstrap analysis supporting this observation in 3 of 4 subjects. We conclude that 1/f slope can be measured longitudinally in a combined SCC-DBS/LFP recording system and likely conforms to known 1/f circadian variability. The preliminary evidence of 1/f slope increase during treatment response suggests a potential utility as a candidate biomarker for ongoing development of adaptive TRD-neuromodulation strategies.NEW & NOTEWORTHY In four patients with treatment-resistant depression undergoing therapeutic deep brain stimulation (DBS), we present the first longitudinal observations of local field potentials (LFP) from the subcallosal cingulate region outside the postoperative period. Specifically, our results demonstrate that frequency-domain 1/f activity is measurable in a combined DBS-LFP recording system and that right hemisphere recordings appear sensitive to mood state, thus suggesting a potential readout suitable for consideration in ongoing efforts to develop adaptive DBS delivery systems.

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 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.

Test-retest reliability of a stimulation-locked evoked response to deep brain stimulation in subcallosal cingulate for treatment resistant depression.

Deep brain stimulation (DBS) to the subcallosal cingulate cortex (SCC) is an emerging therapy for treatment resistant depression. Precision targeting of specific white matter fibers is now central to the model of SCC DBS treatment efficacy. A method to confirm SCC DBS target engagement is needed to reduce procedural variance across treatment providers and to optimize DBS parameters for individual patients. We examined the reliability of a novel cortical evoked response that is time-locked to a 2 Hz DBS pulse and shows the propagation of signal from the DBS target. The evoked response was detected in four individuals as a stereotyped series of components within 150 ms of a 6 V DBS pulse, each showing coherent topography on the head surface. Test-retest reliability across four repeated measures over 14 months met or exceeded standards for valid test construction in three of four patients. Several observations in this pilot sample demonstrate the prospective utility of this method to confirm surgical target engagement and instruct parameter selection. The topography of an orbital frontal component on the head surface showed specificity for patterns of forceps minor activation, which may provide a means to confirm DBS location with respect to key white matter structures. A divergent cortical response to unilateral stimulation of left (vs. right) hemisphere underscores the need for feedback acuity on the level of a single electrode, despite bilateral presentation of therapeutic stimulation. Results demonstrate viability of this method to explore patient-specific cortical responsivity to DBS for brain-circuit pathologies.

Action potential initiation, propagation, and cortical invasion in the hyperdirect pathway during subthalamic deep brain stimulation.

BACKGROUND: High frequency (∼130 Hz) deep brain stimulation (DBS) of the subthalamic region is an established clinical therapy for the treatment of late stage Parkinson's disease (PD). Direct modulation of the hyperdirect pathway, defined as cortical layer V pyramidal neurons that send an axon collateral to the subthalamic nucleus (STN), has emerged as a possible component of the therapeutic mechanisms. However, numerous questions remain to be addressed on the basic biophysics of hyperdirect pathway stimulation. OBJECTIVE: Quantify action potential (AP) initiation, propagation, and cortical invasion in hyperdirect neurons during subthalamic stimulation. METHODS: We developed an anatomically and electrically detailed computational model of hyperdirect neuron stimulation with explicit representation of the stimulating electric field, axonal response, AP propagation, and synaptic transmission. RESULTS: We found robust AP propagation throughout the complex axonal arbor of the hyperdirect neuron. Even at therapeutic DBS frequencies, stimulation induced APs could reach all of the intracortical axon terminals with ∼100% fidelity. The functional result of this high frequency axonal driving of the thousands of synaptic connections made by each directly stimulated hyperdirect neuron is a profound synaptic suppression that would effectively disconnect the neuron from the cortical circuitry. CONCLUSIONS: The synaptic suppression hypothesis integrates the fundamental biophysics of electrical stimulation, axonal transmission, and synaptic physiology to explain a generic mechanism of DBS.

Characterization of the stimulus waveforms generated by implantable pulse generators for deep brain stimulation.

OBJECTIVE: To determine the circuit elements required to theoretically describe the stimulus waveforms generated by an implantable pulse generator (IPG) during clinical deep brain stimulation (DBS). METHODS: We experimentally interrogated the Medtronic Activa PC DBS IPG and defined an equivalent circuit model that accurately captured the output of the IPG. We then compared the detailed circuit model of the clinical stimulus waveforms to simplified representations commonly used in computational models of DBS. We quantified the errors associated with these simplifications using theoretical activation thresholds of myelinated axons in response to DBS. RESULTS: We found that the detailed IPG model generated substantial differences in activation thresholds compared to simplified models. These differences were largest for bipolar stimulation with long pulse widths. Average errors were ∼3 to 24% for voltage-controlled stimulation and ∼2 to 11% for current-controlled stimulation. CONCLUSIONS: Our results demonstrate the importance of including basic circuit elements (e.g. blocking capacitors, lead wire resistance, electrode capacitance) in model analysis of DBS. SIGNIFICANCE: Computational models of DBS are now commonly used in academic research, industrial technology development, and in the selection of clinical stimulation parameters. Incorporating a realistic representation of the IPG output is necessary to improve the accuracy and utility of these clinical and scientific tools.

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.

Stimulation Efficiency With Decaying Exponential Waveforms in a Wirelessly Powered Switched-Capacitor Discharge Stimulation System.

The purpose of this study was to test the feasibility of using a switched-capacitor discharge stimulation (SCDS) system for electrical stimulation, and, subsequently, determine the overall energy saved compared to a conventional stimulator. We have constructed a computational model by pairing an image-based volume conductor model of the cat head with cable models of corticospinal tract (CST) axons and quantified the theoretical stimulation efficiency of rectangular and decaying exponential waveforms, produced by conventional and SCDS systems, respectively. Subsequently, the model predictions were tested in vivo by activating axons in the posterior internal capsule and recording evoked electromyography (EMG) in the contralateral upper arm muscles. Compared to rectangular waveforms, decaying exponential waveforms with time constants >500 µs were predicted to require 2%-4% less stimulus energy to activate directly models of CST axons and 0.4%-2% less stimulus energy to evoke EMG activity in vivo. Using the calculated wireless input energy of the stimulation system and the measured stimulus energies required to evoke EMG activity, we predict that an SCDS implantable pulse generator (IPG) will require 40% less input energy than a conventional IPG to activate target neural elements. A wireless SCDS IPG that is more energy efficient than a conventional IPG will reduce the size of an implant, require that less wireless energy be transmitted through the skin, and extend the lifetime of the battery in the external power transmitter.

Creating and parameterizing patient-specific deep brain stimulation pathway-activation models using the hyperdirect pathway as an example.

BACKGROUND: Deep brain stimulation (DBS) is an established clinical therapy and computational models have played an important role in advancing the technology. Patient-specific DBS models are now common tools in both academic and industrial research, as well as clinical software systems. However, the exact methodology for creating patient-specific DBS models can vary substantially and important technical details are often missing from published reports. OBJECTIVE: Provide a detailed description of the assembly workflow and parameterization of a patient-specific DBS pathway-activation model (PAM) and predict the response of the hyperdirect pathway to clinical stimulation. METHODS: Integration of multiple software tools (e.g. COMSOL, MATLAB, FSL, NEURON, Python) enables the creation and visualization of a DBS PAM. An example DBS PAM was developed using 7T magnetic resonance imaging data from a single unilaterally implanted patient with Parkinson's disease (PD). This detailed description implements our best computational practices and most elaborate parameterization steps, as defined from over a decade of technical evolution. RESULTS: Pathway recruitment curves and strength-duration relationships highlight the non-linear response of axons to changes in the DBS parameter settings. CONCLUSION: Parameterization of patient-specific DBS models can be highly detailed and constrained, thereby providing confidence in the simulation predictions, but at the expense of time demanding technical implementation steps. DBS PAMs represent new tools for investigating possible correlations between brain pathway activation patterns and clinical symptom modulation.

Role of Soft-Tissue Heterogeneity in Computational Models of Deep Brain Stimulation.

BACKGROUND: Bioelectric field models of deep brain stimulation (DBS) are commonly utilized in research and industrial applications. However, the wide range of different representations used for the human head in these models may be responsible for substantial variance in the stimulation predictions. OBJECTIVE: Determine the relative error of ignoring cerebral vasculature and soft-tissue heterogeneity outside of the brain in computational models of DBS. METHODS: We used a detailed atlas of the human head, coupled to magnetic resonance imaging data, to construct a range of subthalamic DBS volume conductor models. We incrementally simplified the most detailed base model and quantified changes in the stimulation thresholds for direct activation of corticofugal axons. RESULTS: Ignoring cerebral vasculature altered predictions of stimulation thresholds by <10%, whereas ignoring soft-tissue heterogeneity outside of the brain altered predictions between -44 % and 174%. CONCLUSIONS: Heterogeneity in the soft tissues of the head, if unaccounted for, introduces a degree of uncertainty in predicting electrical stimulation of neural elements that is not negligible and thereby warrants consideration in future modeling studies.

Analyzing the tradeoff between electrical complexity and accuracy in patient-specific computational models of deep brain stimulation.

OBJECTIVE: Deep brain stimulation (DBS) is an adjunctive therapy that is effective in treating movement disorders and shows promise for treating psychiatric disorders. Computational models of DBS have begun to be utilized as tools to optimize the therapy. Despite advancements in the anatomical accuracy of these models, there is still uncertainty as to what level of electrical complexity is adequate for modeling the electric field in the brain and the subsequent neural response to the stimulation. APPROACH: We used magnetic resonance images to create an image-based computational model of subthalamic DBS. The complexity of the volume conductor model was increased by incrementally including heterogeneity, anisotropy, and dielectric dispersion in the electrical properties of the brain. We quantified changes in the load of the electrode, the electric potential distribution, and stimulation thresholds of descending corticofugal (DCF) axon models. MAIN RESULTS: Incorporation of heterogeneity altered the electric potentials and subsequent stimulation thresholds, but to a lesser degree than incorporation of anisotropy. Additionally, the results were sensitive to the choice of method for defining anisotropy, with stimulation thresholds of DCF axons changing by as much as 190%. Typical approaches for defining anisotropy underestimate the expected load of the stimulation electrode, which led to underestimation of the extent of stimulation. More accurate predictions of the electrode load were achieved with alternative approaches for defining anisotropy. The effects of dielectric dispersion were small compared to the effects of heterogeneity and anisotropy. SIGNIFICANCE: The results of this study help delineate the level of detail that is required to accurately model electric fields generated by DBS electrodes.