Objective. Deep brain stimulation (DBS) is a growing treatment option for movement and psychiatric disorders. As DBS technology moves toward directional leads with increased numbers of smaller ...electrode contacts, trial-and-error methods of manual DBS programming are becoming too time-consuming for clinical feasibility. We propose an algorithm to automate DBS programming in near real-time for a wide range of DBS lead designs. Approach. Magnetic resonance imaging and diffusion tensor imaging are used to build finite element models that include anisotropic conductivity. The algorithm maximizes activation of target tissue and utilizes the Hessian matrix of the electric potential to approximate activation of neurons in all directions. We demonstrate our algorithm's ability in an example programming case that targets the subthalamic nucleus (STN) for the treatment of Parkinson's disease for three lead designs: the Medtronic 3389 (four cylindrical contacts), the direct STNAcute (two cylindrical contacts, six directional contacts), and the Medtronic-Sapiens lead (40 directional contacts). Main results. The optimization algorithm returns patient-specific contact configurations in near real-time-less than 10 s for even the most complex leads. When the lead was placed centrally in the target STN, the directional leads were able to activate over 50% of the region, whereas the Medtronic 3389 could activate only 40%. When the lead was placed 2 mm lateral to the target, the directional leads performed as well as they did in the central position, but the Medtronic 3389 activated only 2.9% of the STN. Significance. This DBS programming algorithm can be applied to cylindrical electrodes as well as novel directional leads that are too complex with modern technology to be manually programmed. This algorithm may reduce clinical programming time and encourage the use of directional leads, since they activate a larger volume of the target area than cylindrical electrodes in central and off-target lead placements.
Achieving deep brain stimulation (DBS) dose equivalence is challenging, especially with pulse width tuning and directional contacts. Further, the precise effects of pulse width tuning are unknown, ...and recent reports of the effects of pulse width tuning on neural selectivity are at odds with classic biophysical studies.
We created multicompartment neuron models for two axon diameters and used finite element modeling to determine extracellular influence from standard and segmented electrodes. We analyzed axon activation profiles and calculated volumes of tissue activated.
We find that long pulse widths focus the stimulation effect on small, nearby fibers, suppressing distant white matter tract activation (responsible for some DBS side effects) and improving battery utilization when equivalent activation is maintained for small axons. Directional leads enable similar benefits to a greater degree. Reexamining previous reports of short pulse stimulation reducing side effects, we explore a possible alternate explanation: non-dose equivalent stimulation may have resulted in reduced spread of neural activation. Finally, using internal capsule avoidance as an example in the context of subthalamic stimulation, we present a patient-specific model to show how long pulse widths could help increase the biophysical therapeutic window.
We find agreement with classic studies and predict that long pulse widths may focus the stimulation effect on small, nearby fibers and improve power consumption. While future pre-clinical and clinical work is necessary regarding pulse width tuning, it is clear that future studies must ensure dose equivalence, noting that energy- and charge-equivalent amplitudes do not result in equivalent spread of neural activation when changing pulse width.
•Long pulse widths are predicted to focus the stimulation effect on small, nearby axons, which may improve therapy.•Long pulse widths, up to several hundred microseconds, improve energetic efficiency.•Segmented electrodes may increase both therapeutic window and energy efficiency more effectively than pulse width tuning.•When tuning pulse width, energy- and charge-equivalent amplitudes generate differently sized volumes of tissue activated.•Future pulse width studies should select amplitudes via strength-duration metrics to maintain spread of neural activation.
We describe an electrical "running down" phenomenon and also a consistent spectral change (in the aperiodic component of the power spectrum) derived from chronic interictal electrocorticography ...(ECoG) after surgery in a patient with drug‐resistant epilepsy. These data were recorded using a closed‐loop neurostimulation system that was implanted after resection. The patient has been seizure‐free for 2.5 years since resection without requiring the neurostimulator to be turned on to stimulate. Concurrently, there was an exponential decrease in the number of epileptiform electrographic detections recorded by the device, particularly over the first 26 weeks, indicative of an electrical running down phenomenon as the brain adapted to an extended period of seizure freedom. We also find that the aperiodic exponent of the power spectrum gradually decreases over time. The aperiodic component of intracranial ECoG may represent a novel marker of epileptogenicity, independent of seizures.
In this study, we quantified the coverage of gray and white matter during intracranial electroencephalography in a cohort of epilepsy patients with surface and depth electrodes. We included 65 ...patients with strip electrodes (n = 12), strip and grid electrodes (n = 24), strip, grid, and depth electrodes (n = 7), or depth electrodes only (n = 22). Patient-specific imaging was used to generate probabilistic gray and white matter maps and atlas segmentations. Gray and white matter coverage was quantified using spherical volumes centered on electrode centroids, with radii ranging from 1 to 15 mm, along with detailed finite element models of local electric fields. Gray matter coverage was highly dependent on the chosen radius of influence (RoI). Using a 2.5 mm RoI, depth electrodes covered more gray matter than surface electrodes; however, surface electrodes covered more gray matter at RoI larger than 4 mm. White matter coverage and amygdala and hippocampal coverage was greatest for depth electrodes at all RoIs. This study provides the first probabilistic analysis to quantify coverage for different intracranial recording configurations. Depth electrodes offer increased coverage of gray matter over other recording strategies if the desired signals are local, while subdural grids and strips sample more gray matter if the desired signals are diffuse.
Accurate anatomical localization of intracranial electrodes is important for identifying the seizure foci in patients with epilepsy and for interpreting effects from cognitive studies employing ...intracranial electroencephalography. Localization is typically performed by coregistering postimplant computed tomography (CT) with preoperative magnetic resonance imaging (MRI). Electrodes are then detected in the CT, and the corresponding brain region is identified using the MRI. Many existing software packages for electrode localization chain together separate preexisting programs or rely on command line instructions to perform the various localization steps, making them difficult to install and operate for a typical user. Further, many packages provide solutions for some, but not all, of the steps needed for confident localization. We have developed software, Locate electrodes Graphical User Interface (LeGUI), that consists of a single interface to perform all steps needed to localize both surface and depth/penetrating intracranial electrodes, including coregistration of the CT to MRI, normalization of the MRI to the Montreal Neurological Institute template, automated electrode detection for multiple types of electrodes, electrode spacing correction and projection to the brain surface, electrode labeling, and anatomical targeting. The software is written in MATLAB, core image processing is performed using the Statistical Parametric Mapping toolbox, and standalone executable binaries are available for Windows, Mac, and Linux platforms. LeGUI was tested and validated on 51 datasets from two universities. The total user and computational time required to process a single dataset was approximately 1 h. Automatic electrode detection correctly identified 4362 of 4695 surface and depth electrodes with only 71 false positives. Anatomical targeting was verified by comparing electrode locations from LeGUI to locations that were assigned by an experienced neuroanatomist. LeGUI showed a 94% match with the 482 neuroanatomist-assigned locations. LeGUI combines all the features needed for fast and accurate anatomical localization of intracranial electrodes into a single interface, making it a valuable tool for intracranial electrophysiology research.
Directional deep brain stimulation (DBS) leads have recently been approved and used in patients, and growing evidence suggests that directional contacts can increase the therapeutic window by ...redirecting stimulation to the target region while avoiding side-effect-inducing regions. We outline the design, fabrication, and testing of a novel directional DBS lead, the µDBS, which utilizes microscale contacts to increase the spatial resolution of stimulation steering and improve the selectivity in targeting small diameter fibers. We outline the steps of fabrication of the µDBS, from an integrated circuit design to postprocessing and validation testing. We tested the onboard digital circuitry for programming fidelity, characterized impedance for a variety of electrode sizes, and demonstrated functionality in a saline bath. In a computational experiment, we determined that reduced electrode sizes focus the stimulation effect on small, nearby fibers. Smaller electrode sizes allow for a relative decrease in small-diameter axon thresholds compared to thresholds of large-diameter fibers, demonstrating a focusing of the stimulation effect within small, and possibly therapeutic, fibers. This principle of selectivity could be useful in further widening the window of therapy. The µDBS offers a unique, multiresolution design in which any combination of microscale contacts can be used together to function as electrodes of various shapes and sizes. Multiscale electrodes could be useful in selective neural targeting for established neurological targets and in exploring novel treatment targets for new neurological indications.
The degeneration of dopaminergic neurons in the substantia nigra pars compacta leads to parkinsonian motor symptoms via changes in electrophysiological activity throughout the basal ganglia. ...High-frequency deep brain stimulation (DBS) partially treats these symptoms, but the mechanisms are unclear. We hypothesize that motor symptoms of Parkinson's disease (PD) are associated with increased information transmission from basal ganglia output neurons to motor thalamus input neurons and that therapeutic DBS of the subthalamic nucleus (STN) treats these symptoms by reducing this extraneous information transmission. We tested these hypotheses in a unilateral, 6-hydroxydopamine-lesioned rodent model of hemiparkinsonism. Information transfer between basal ganglia output neurons and motor thalamus input neurons increased in both the orthodromic and antidromic directions with hemiparkinsonian (hPD) onset, and these changes were reversed by behaviorally therapeutic STN-DBS. Omnidirectional information increases in the parkinsonian state underscore the detrimental nature of that pathological information and suggest a loss of information channel independence. Therapeutic STN-DBS reduced that pathological information, suggesting an effective increase in the number of independent information channels. We interpret these data with a model in which pathological information and fewer information channels diminishes the scope of possible motor activities, driving parkinsonian symptoms. In this model, STN-DBS restores information-channel independence by eliminating or masking the parkinsonism-associated information, and thus enlarges the scope of possible motor activities, alleviating parkinsonian symptoms.
Objective. Computational models are a popular tool for predicting the effects of deep brain stimulation (DBS) on neural tissue. One commonly used model, the volume of tissue activated (VTA), is ...computed using multiple methodologies. We quantified differences in the VTAs generated by five methodologies: the traditional axon model method, the electric field norm, and three activating function based approaches-the activating function at each grid point in the tangential direction (AF-Tan) or in the maximally activating direction (AF-3D), and the maximum activating function along the entire length of a tangential fiber (AF-Max). Approach. We computed the VTA using each method across multiple stimulation settings. The resulting volumes were compared for similarity, and the methodologies were analyzed for their differences in behavior. Main results. Activation threshold values for both the electric field norm and the activating function varied with regards to electrode configuration, pulse width, and frequency. All methods produced highly similar volumes for monopolar stimulation. For bipolar electrode configurations, only the maximum activating function along the tangential axon method, AF-Max, produced similar volumes to those produced by the axon model method. Further analysis revealed that both of these methods are biased by their exclusive use of tangential fiber orientations. In contrast, the activating function in the maximally activating direction method, AF-3D, produces a VTA that is free of axon orientation and projection bias. Significance. Simulating tangentially oriented axons, the standard approach of computing the VTA, is too computationally expensive for widespread implementation and yields results biased by the assumption of tangential fiber orientation. In this work, we show that a computationally efficient method based on the activating function, AF-Max, reliably reproduces the VTAs generated by direct axon modeling. Further, we propose another method, AF-3D as a potentially superior model for representing generic neural tissue activation.
Objective. During deep brain stimulation (DBS), it is well understood that extracellular cathodic stimulation can cause activation of passing axons. Activation can be predicted from the second ...derivative of the electric potential along an axon, which depends on axonal orientation with respect to the stimulation source. We hypothesize that fiber orientation influences activation thresholds and that fiber orientations can be selectively targeted with DBS waveforms. Approach. We used bioelectric field and multicompartment NEURON models to explore preferential activation based on fiber orientation during monopolar or bipolar stimulation. Preferential fiber orientation was extracted from the principal eigenvectors and eigenvalues of the Hessian matrix of the electric potential. We tested cathodic, anodic, and charge-balanced pulses to target neurons based on fiber orientation in general and clinical scenarios. Main results. Axons passing the DBS lead have positive second derivatives around a cathode, whereas orthogonal axons have positive second derivatives around an anode, as indicated by the Hessian. Multicompartment NEURON models confirm that passing fibers are activated by cathodic stimulation, and orthogonal fibers are activated by anodic stimulation. Additionally, orthogonal axons have lower thresholds compared to passing axons. In a clinical scenario, fiber pathways associated with therapeutic benefit can be targeted with anodic stimulation at 50% lower stimulation amplitudes. Significance. Fiber orientations can be selectively targeted with simple changes to the stimulus waveform. Anodic stimulation preferentially activates orthogonal fibers, approaching or leaving the electrode, at lower thresholds for similar therapeutic benefit in DBS with decreased power consumption.