Neural pathways associated with reduced rigidity during pallidal deep brain stimulation for Parkinson's disease

Deep brain stimulation (DBS) of the internal segment of the globus pallidus (GPi) can markedly reduce muscle rigidity in people with Parkinson's disease (PD); however, the mechanisms mediating this effect are poorly understood. Computational modeling of DBS provides a method to estimate the relative contributions of neural pathway activations to changes in outcomes. In this study, we generated subject-specific biophysical models of GPi DBS (derived from individual 7-T MRI), including pallidal efferent, putamenal efferent, and internal capsule pathways, to investigate how activation of neural pathways contributed to changes in forearm rigidity in PD. Ten individuals (17 arms) were tested off medication under four conditions: off stimulation, on clinically optimized stimulation, and on stimulation specifically targeting the dorsal GPi or ventral GPi. Quantitative measures of forearm rigidity, with and without a contralateral activation maneuver, were obtained with a robotic manipulandum. Clinically optimized GPi DBS settings significantly reduced forearm rigidity (P < 0.001), which aligned with GPi efferent fiber activation. The model demonstrated that GPi efferent axons could be activated at any location along the GPi dorsal-ventral axis. These results provide evidence that rigidity reduction produced by GPi DBS is mediated by preferential activation of GPi efferents to the thalamus, likely leading to a reduction in excitability of the muscle stretch reflex via overdriving pallidofugal output.NEW & NOTEWORTHY Subject-specific computational models of pallidal deep brain stimulation, in conjunction with quantitative measures of forearm rigidity, were used to examine the neural pathways mediating stimulation-induced changes in rigidity in people with Parkinson's disease. The model uniquely included internal, efferent and adjacent pathways of the basal ganglia. The results demonstrate that reductions in rigidity evoked by deep brain stimulation were principally mediated by the activation of globus pallidus internus efferent pathways.

Effect of Anisotropic Brain Conductivity on Patient-Specific Volume of Tissue Activation in Deep Brain Stimulation for Parkinson Disease

OBJECTIVE

Deep brain stimulation (DBS) modeling can improve surgical targeting by quantifying the spatial extent of stimulation relative to subcortical structures of interest. A certain degree of model complexity is required to obtain accurate predictions, particularly complexity regarding electrical properties of the tissue around DBS electrodes. In this study, the effect of anisotropy on the volume of tissue activation (VTA) was evaluated in an individualized manner.

METHODS

Tissue activation models incorporating patient-specific tissue conductivity were built for 40 Parkinson disease patients who had received bilateral subthalamic nucleus (STN) DBS. To assess the impact of local changes in tissue anisotropy, one VTA was computed at each electrode contact using identical stimulation parameters. For comparison, VTAs were also computed assuming isotropic tissue conductivity. Stimulation location was considered by classifying the anisotropic VTAs relative to the STN. VTAs were characterized based on volume, spread in three directions, sphericity, and Dice coefficient.

RESULTS

Incorporating anisotropy generated significantly larger and less spherical VTAs overall. However, its effect on VTA size and shape was variable and more nuanced at the individual patient and implantation levels. Dorsal VTAs had significantly higher sphericity than ventral VTAs, suggesting more isotropic behavior. Contrastingly, lateral and posterior VTAs had significantly larger and smaller lateral-medial spreads, respectively. Volume and spread correlated negatively with sphericity.

CONCLUSION

The influence of anisotropy on VTA predictions is important to consider, and varies across patients and stimulation location.

SIGNIFICANCE

This study highlights the importance of considering individualized factors in DBS modeling to accurately characterize the VTA.

Lateral cerebellothalamic tract activation underlies DBS therapy for Essential Tremor

BACKGROUND

While deep brain stimulation (DBS) therapy can be effective at suppressing tremor in individuals with medication-refractory Essential Tremor, patient outcome variability remains a significant challenge across centers. Proximity of active electrodes to the cerebellothalamic tract (CTT) is likely important in suppressing tremor, but how tremor control and side effects relate to targeting parcellations within the CTT and other pathways in and around the ventral intermediate (VIM) nucleus of thalamus remain unclear.

METHODS

Using ultra-high field (7T) MRI, we developed high-dimensional, subject-specific pathway activation models for 23 directional DBS leads. Modeled pathway activations were compared with post-hoc analysis of clinician-optimized DBS settings, paresthesia thresholds, and dysarthria thresholds. Mixed-effect models were utilized to determine how the six parcellated regions of the CTT and how six other pathways in and around the VIM contributed to tremor suppression and induction of side effects.

RESULTS

The lateral portion of the CTT had the highest activation at clinical settings (p < 0.05) and a significant effect on tremor suppression (p < 0.001). Activation of the medial lemniscus and posterior-medial CTT was significantly associated with severity of paresthesias (p < 0.001). Activation of the anterior-medial CTT had a significant association with dysarthria (p < 0.05).

CONCLUSIONS

This study provides a detailed understanding of the fiber pathways responsible for therapy and side effects of DBS for Essential Tremor, and suggests a model-based programming approach will enable more selective activation of lateral fibers within the CTT.

Deep Brain Stimulation: Emerging Tools for Simulation, Data Analysis, and Visualization

Deep brain stimulation (DBS) is a well-established neurosurgical procedure for movement disorders that is also being explored for treatment-resistant psychiatric conditions. This review highlights important consideration for DBS simulation and data analysis. The literature on DBS has expanded considerably in recent years, and this article aims to identify important trends in the field. During DBS planning, surgery, and follow up sessions, several large data sets are created for each patient, and it becomes clear that any group analysis of such data is a big data analysis problem and has to be handled with care. The aim of this review is to provide an update and overview from a neuroengineering perspective of the current DBS techniques, technical aids, and emerging tools with the focus on patient-specific electric field (EF) simulations, group analysis, and visualization in the DBS domain. Examples are given from the state-of-the-art literature including our own research. This work reviews different analysis methods for EF simulations, tractography, deep brain anatomical templates, and group analysis. Our analysis highlights that group analysis in DBS is a complex multi-level problem and selected parameters will highly influence the result. DBS analysis can only provide clinically relevant information if the EF simulations, tractography results, and derived brain atlases are based on as much patient-specific data as possible. A trend in DBS research is creation of more advanced and intuitive visualization of the complex analysis results suitable for the clinical environment.

Modeling implanted metals in electrical stimulation applications

OBJECTIVE

Metal implants impact the dosimetry assessment in electrical stimulation techniques. Therefore, they need to be included in numerical models. While currents in the body are ionic, metals only allow electron transport. In fact, charge transfer between tissues and metals requires electric fields to drive electrochemical reactions at the interface. Thus, metal implants may act as insulators or as conductors depending on the scenario. The aim of this paper is to provide a theoretical argument that guides the choice of the correct representation of metal implants in electrical models while considering the electrochemical nature of the problem Approach: We built a simple model of a metal implant exposed to a homogeneous electric field of various magnitudes. The same geometry was solved using two different models: a purely electric one (with different conductivities for the implant), and an electrochemical one. As an example of application, we also modeled a transcranial electrical stimulation (tES) treatment in a realistic head model with a skull plate using a high and low conductivity value for the plate.

MAIN RESULTS

Metal implants generally act as electric insulators when exposed to electric fields up to around 100 V/m and they only resemble a perfect conductor for fields in the order of 1000 V/m and above. The results are independent of the implant's metal, but they depend on its geometry. tES modeling with implants incorrectly treated as conductors can lead to errors of 50% or more in the estimation of the induced fields Significance: Metal implants can be accurately represented by a simple electrical model of constant conductivity, but an incorrect model choice can lead to large errors in the dosimetry assessment. Our results can be used to guide the selection of the most appropriate model in each scenario.

Deep brain stimulation electrode modeling in rats

Deep Brain Stimulation (DBS) is an efficacious treatment option for an increasing range of brain disorders. To enhance our knowledge about the mechanisms of action of DBS and to probe novel targets, basic research in animal models with DBS is an essential research base. Beyond nonhuman primate, pig, and mouse models, the rat is a widely used animal model for probing DBS effects in basic research. Reconstructing DBS electrode placement after surgery is crucial to associate observed effects with modulating a specific target structure. Post-mortem histology is a commonly used method for reconstructing the electrode location. In humans, however, neuroimaging-based electrode localizations have become established. For this reason, we adapt the open-source software pipeline Lead-DBS for DBS electrode localizations from humans to the rat model. We validate our localization results by inter-rater concordance and a comparison with the conventional histological method. Finally, using the open-source software pipeline OSS-DBS, we demonstrate the subject-specific simulation of the VTA and the activation of axon models aligned to pathways representing neuronal fibers, also known as the pathway activation model. Both activation models yield a characterization of the impact of DBS on the target area. Our results suggest that the proposed neuroimaging-based method can precisely localize DBS electrode placements that are essentially rater-independent and yield results comparable to the histological gold standard. The advantages of neuroimaging-based electrode localizations are the possibility of acquiring them in vivo and combining electrode reconstructions with advanced imaging metrics, such as those obtained from diffusion or functional magnetic resonance imaging (MRI). This paper introduces a freely available open-source pipeline for DBS electrode reconstructions in rats. The presented initial validation results are promising.

Estimation of electric field impact in deep brain stimulation from axon diameter distribution in the human brain

Objective.Finite element method (FEM) simulations of the electric field magnitude (EF) are commonly used to estimate the affected tissue surrounding the active contact of deep brain stimulation (DBS) leads. Previous studies have found that DBS starts to noticeably activate axons at approximately 0.2 V mm^-1, corresponding to activation of 3.4μm axons in simulations of individual axon triggering. Most axons in the brain are considerably smaller however, and the effect of the electric field is thus expected to be stronger with increasing EF as more and more axons become activated. The objective of this study is to estimate the fraction of activated axons as a function of electric field magnitude.Approach. The EF thresholds required for axon stimulation of myelinated axon diameters between 1 and 5μm were obtained from a combined cable and Hodgkin-Huxley model in a FEM-simulated electric field from a Medtronic 3389 lead. These thresholds were compared with the average axon diameter distribution from literature from several structures in the human brain to obtain an estimate of the fraction of axons activated at EF levels between 0.1 and 1.8 V mm^-1.Main results. The effect of DBS is estimated to be 47·EF-8.8% starting at a threshold levelEF_t0 = 0.19 V mm^-1.Significance. The fraction of activated axons from DBS in a voxel is estimated to increase linearly with EF above the threshold level of 0.19 V mm^-1. This means linear regression between EF above 0.19 V mm^-1and clinical outcome is a suitable statistical method when doing improvement maps for DBS.

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.

Impact of brain shift on neural pathways in deep brain stimulation: a preliminary analysis via multi-physics finite element models

OBJECTIVE

The effectiveness of deep brain stimulation (DBS) depends on electrode placement accuracy, which can be compromised by brain shift during surgery. While there have been efforts in assessing the impact of electrode misplacement due to brain shift using preop- and postop- imaging data, such analysis using preop- and intraop- imaging data via biophysical modeling has not been conducted. This work presents a preliminary study that applies a multi-physics analysis framework using finite element biomechanical and bioelectric models to examine the impact of realistic intraoperative shift on neural pathways determined by tractography.

APPROACH

The study examined six patients who had undergone interventional magnetic resonance (iMR)-guided DBS surgery. The modeling framework utilized a biomechanical approach to update preoperative MR to reflect shift-induced anatomical changes. Using this anatomically deformed image and its undeformed counterpart, bioelectric effects from shifting electrode leads could be simulated and neural activation differences were approximated. Specifically, for each configuration, volume of tissue activation (VTA) was computed and subsequently used for tractography estimation. Total tract volume and overlapping volume with motor regions as well as connectivity profile were compared. In addition, volumetric overlap between different fiber bundles among configurations was computed and correlated to estimated shift.

MAIN RESULT

The study found deformation-induced differences in tract volume, motor region overlap, and connectivity behavior, suggesting the impact of shift. There is a strong correlation (R=-0.83) between shift from intended target and intended neural pathway recruitment, where at threshold of ~2.94 mm, intended recruitment completely degrades. The determined threshold is consistent with and provides quantitative support to prior observations and literature that deviations of 2-3 mm are detrimental.

SIGNIFICANCE

The findings support and advance prior studies and understanding to illustrate the need to account for shift in DBS and the potentiality of computational modeling for estimating influence of shift on neural activation.

Atlas-independent, N-of-1 tissue activation modeling to map optimal regions of subthalamic deep brain stimulation for Parkinson disease

•

Neuroanatomical variations among patients are obscured in atlas-based VTA modeling.

•

N-of-1 neuroanatomical and VTA modeling enables patient-level precision.

•

Mean optimal stimulation is dorsomedial to the STN, near its posterior half.

•

Individual VTAs deviate from optimal stimulation sites to varying degrees.

•

Optimal stimulation sites for rigidity, bradykinesia, and tremor partially overlap.

Background

Motor outcomes after subthalamic deep brain stimulation (STN DBS) for Parkinson disease (PD) vary considerably among patients and strongly depend on stimulation location. The objective of this retrospective study was to map the regions of optimal STN DBS for PD using an atlas-independent, fully individualized (N-of-1) tissue activation modeling approach and to assess the relationship between patient-level therapeutic volumes of tissue activation (VTAs) and motor improvement.

Methods

The stimulation-induced electric field for 40 PD patients treated with bilateral STN DBS was modeled using finite element analysis. Neurostimulation models were generated for each patient, incorporating their individual STN anatomy, DBS lead position and orientation, anisotropic tissue conductivity, and clinical stimulation settings. A voxel-based analysis of the VTAs was then used to map the optimal location of stimulation. The amount of stimulation in specific regions relative to the STN was measured and compared between STNs with more and less optimal stimulation, as determined by their motor improvement scores and VTA. The relationship between VTA location and motor outcome was then assessed using correlation analysis. Patient variability in terms of STN anatomy, active contact position, and VTA location were also evaluated. Results from the N-of-1 model were compared to those from a simplified VTA model.

Results

Tissue activation modeling mapped the optimal location of stimulation to regions medial, posterior, and dorsal to the STN centroid. These regions extended beyond the STN boundary towards the caudal zona incerta (cZI). The location of the VTA and active contact position differed significantly between STNs with more and less optimal stimulation in the dorsal-ventral and anterior-posterior directions. Therapeutic stimulation spread noticeably more in the dorsal and posterior directions, providing additional evidence for cZI as an important DBS target. There were significant linear relationships between the amount of dorsal and posterior stimulation, as measured by the VTA, and motor improvement. These relationships were more robust than those between active contact position and motor improvement. There was high variability in STN anatomy, active contact position, and VTA location among patients. Spherical VTA modeling was unable to reproduce these results and tended to overestimate the size of the VTA.

Conclusion

Accurate characterization of the spread of stimulation is needed to optimize STN DBS for PD. High variability in neuroanatomy, stimulation location, and motor improvement among patients highlights the need for individualized modeling techniques. The atlas-independent, N-of-1 tissue activation modeling approach presented in this study can be used to develop and evaluate stimulation strategies to improve clinical outcomes on an individual basis.

Orientation-selective and directional deep brain stimulation in swine assessed by functional MRI at 3T

Functional MRI (fMRI) has become an important tool for probing network-level effects of deep brain stimulation (DBS). Previous DBS-fMRI studies have shown that electrical stimulation of the ventrolateral (VL) thalamus can modulate sensorimotor cortices in a frequency and amplitude dependent manner. Here, we investigated, using a swine animal model, how the direction and orientation of the electric field, induced by VL-thalamus DBS, affects activity in the sensorimotor cortex. Adult swine underwent implantation of a novel 16-electrode (4 rows × 4 columns) directional DBS lead in the VL thalamus. A within-subject design was used to compare fMRI responses for (1) directional stimulation consisting of monopolar stimulation in four radial directions around the DBS lead, and (2) orientation-selective stimulation where an electric field dipole was rotated 0°−360° around a quadrangle of electrodes. Functional responses were quantified in the premotor, primary motor, and somatosensory cortices. High frequency electrical stimulation through leads implanted in the VL thalamus induced directional tuning in cortical response patterns to varying degrees depending on DBS lead position. Orientation-selective stimulation showed maximal functional response when the electric field was oriented approximately parallel to the DBS lead, which is consistent with known axonal orientations of the cortico-thalamocortical pathway. These results demonstrate that directional and orientation-selective stimulation paradigms in the VL thalamus can tune network-level modulation patterns in the sensorimotor cortex, which may have translational utility in improving functional outcomes of DBS therapy.

OSS-DBS: Open-source simulation platform for deep brain stimulation with a comprehensive automated modeling

In this study, we propose a new open-source simulation platform that comprises computer-aided design and computer-aided engineering tools for highly automated evaluation of electric field distribution and neural activation during Deep Brain Stimulation (DBS). It will be shown how a Volume Conductor Model (VCM) is constructed and examined using Python-controlled algorithms for generation, discretization and adaptive mesh refinement of the computational domain, as well as for incorporation of heterogeneous and anisotropic properties of the tissue and allocation of neuron models. The utilization of the platform is facilitated by a collection of predefined input setups and quick visualization routines. The accuracy of a VCM, created and optimized by the platform, was estimated by comparison with a commercial software. The results demonstrate no significant deviation between the models in the electric potential distribution. A qualitative estimation of different physics for the VCM shows an agreement with previous computational studies. The proposed computational platform is suitable for an accurate estimation of electric fields during DBS in scientific modeling studies. In future, we intend to acquire SDA and EMA approval. Successful incorporation of open-source software, controlled by in-house developed algorithms, provides a highly automated solution. The platform allows for optimization and uncertainty quantification (UQ) studies, while employment of the open-source software facilitates accessibility and reproducibility of simulations.

Development and testing of implanted carbon electrodes for electromagnetic field mapping during neuromodulation

PURPOSE

Deep brain stimulation electrodes composed of carbon fibers were tested as a means of administering and imaging magnetic resonance electrical impedance tomography (MREIT) currents. Artifacts and heating properties of custom carbon-fiber deep brain stimulation (DBS) electrodes were compared with those produced with standard DBS electrodes.

METHODS

Electrodes were constructed from multiple strands of 7-μm carbon-fiber stock. The insulated carbon electrodes were matched to DBS electrode diameter and contact areas. Images of DBS and carbon electrodes were collected with and without current flow and were compared in terms of artifact and thermal effects in phantoms or tissue samples in 7T imaging conditions. Effects on magnetic flux density and current density distributions were also assessed.

RESULTS

Carbon electrodes produced magnitude artifacts with smaller FWHM values compared to the magnitude artifacts around DBS electrodes in spin echo and gradient echo imaging protocols. DBS electrodes appeared 269% larger than actual size in gradient echo images, in sharp contrast to the negligible artifact observed in diameter-matched carbon electrodes. As expected, larger temperature changes were observed near DBS electrodes during extended RF excitations compared with carbon electrodes in the same phantom. Magnitudes and distribution of magnetic flux density and current density reconstructions were comparable for carbon and DBS electrodes.

CONCLUSION

Carbon electrodes may offer a safer, MR-compatible method for administering neuromodulation currents. Use of carbon-fiber electrodes should allow imaging of structures close to electrodes, potentially allowing better targeting, electrode position revision, and the facilitation of functional imaging near electrodes during neuromodulation.

A high-performance 4 nV (√Hz)-1 analog front-end architecture for artefact suppression in local field potential recordings during deep brain stimulation

OBJECTIVE

Recording of local field potentials (LFPs) during deep brain stimulation (DBS) is necessary to investigate the instantaneous brain response to stimulation, minimize time delays for closed-loop neurostimulation and maximise the available neural data. To our knowledge, existing recording systems lack the ability to provide artefact-free high-frequency (>100 Hz) LFP recordings during DBS in real time primarily because of the contamination of the neural signals of interest by the stimulation artefacts.

APPROACH

To solve this problem, we designed and developed a novel, low-noise and versatile analog front-end (AFE) that uses a high-order (8th) analog Chebyshev notch filter to suppress the artefacts originating from the stimulation frequency. After defining the system requirements for concurrent LFP recording and DBS artefact suppression, we assessed the performance of the realised AFE by conducting both in vitro and in vivo experiments using unipolar and bipolar DBS (monophasic pulses, amplitude ranging from 3 to 6 V peak-to-peak, frequency 140 Hz and pulse width 100 µs). A full performance comparison between the proposed AFE and an identical AFE, equipped with an 8th order analog Bessel notch filter, was also conducted.

MAIN RESULTS

A high-performance, 4 nV ([Formula: see text])^-1 AFE that is capable of recording nV-scale signals was designed in accordance with the imposed specifications. Under both in vitro and in vivo experimental conditions, the proposed AFE provided real-time, low-noise and artefact-free LFP recordings (in the frequency range 0.5-250 Hz) during stimulation. Its sensing and stimulation artefact suppression capabilities outperformed the capabilities of the AFE equipped with the Bessel notch filter.

SIGNIFICANCE

The designed AFE can precisely record LFP signals, in and without the presence of either unipolar or bipolar DBS, which renders it as a functional and practical AFE architecture to be utilised in a wide range of applications and environments. This work paves the way for the development of externalized research tools for closed-loop neuromodulation that use low- and higher-frequency LFPs as control signals.

Optimized programming algorithm for cylindrical and directional deep brain stimulation electrodes

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.

Memory Prosthesis: Is It Time for a Deep Neuromimetic Computing Approach?

Memory loss, one of the most dreaded afflictions of the human condition, presents considerable burden on the world's health care system and it is recognized as a major challenge in the elderly. There are only a few neuromodulation treatments for memory dysfunctions. Open loop deep brain stimulation is such a treatment for memory improvement, but with limited success and conflicting results. In recent years closed-loop neuroprosthesis systems able to simultaneously record signals during behavioral tasks and generate with the use of internal neural factors the precise timing of stimulation patterns are presented as attractive alternatives and show promise in memory enhancement and restoration. A few such strides have already been made in both animals and humans, but with limited insights into their mechanisms of action. Here, I discuss why a deep neuromimetic computing approach linking multiple levels of description, mimicking the dynamics of brain circuits, interfaced with recording and stimulating electrodes could enhance the performance of current memory prosthesis systems, shed light into the neurobiology of learning and memory and accelerate the progress of memory prosthesis research. I propose what the necessary components (nodes, structure, connectivity, learning rules, and physiological responses) of such a deep neuromimetic model should be and what type of data are required to train/test its performance, so it can be used as a true substitute of damaged brain areas capable of restoring/enhancing their missing memory formation capabilities. Considerations to neural circuit targeting, tissue interfacing, electrode placement/implantation, and multi-network interactions in complex cognition are also provided.

Energy Cost of Action Potential Generation and Propagation in Thalamocortical Relay Neurons During Deep Brain Stimulation

Thalamocortical (TC) relay neurons generate antidromic and orthodromic action potentials (APs) during thalamic deep brain stimulation (DBS). To maintain signaling, each AP requires Na⁺/K⁺ pump to expend adenosine triphosphate (ATP) to restore Na⁺ and K⁺ gradients. Our aim was to estimate the energy demand associated with AP generation and propagation within TC relay cells during DBS. We used a morphology-based computational model to simulate the APs at different locations. We determined AP energy cost by calculating the amount of ATP required to reverse Na⁺ influx during the spike and measured metabolic efficiency by using Na⁺/K⁺ charge overlap. The ATP cost for AP generation exhibited location dependence, which was determined by spike shape, spatial morphology, and heterogeneously distributed currents. The APs in the axonal initial segment (AIS) were energetically efficient, but backpropagation to the soma and forward propagation to the axon were inefficient. Due to large surface area, the soma and AIS dominated the overall ATP usage. The AP cost also depended on membrane potential, which controlled T-type Ca²⁺ conductance and degree of availability of Na⁺ and K⁺ channels. The excitatory/inhibitory synaptic inputs affected spike cost by increasing/reducing the excitability of local cells. There was a tradeoff between AP cost and firing rate at high firing frequencies. We explained a fundamental link between biophysics of ionic currents, spatial morphology of neural segments, and ATP cost per AP. The predictions should be considered when understanding the functional magnetic resonance imaging data of thalamic DBS.

Modeling of Brain Tissue Heating Caused by Direct Cortical Stimulation for Assessing the Risk of Thermal Damage

This paper aims to employ the numerical simulations to assess the risk of cellular damage during the application of a novel paradigm of electrical stimulation mapping (ESM) used in neurosurgery. The core principle of the paradigm is the use of short, high-intensity and high-frequency stimulation pulses. We developed a complex numerical model and performed coupled electro-thermal transient simulations. The model was optimized by incorporating ESM electrodes' resistance obtained during multiple intraoperative measurements and validated by comparing them with the results of temperature distribution measurement acquired by thermal imaging. The risk of heat-induced cellular damage was assessed by applying the Arrhenius equation integral on the computed time-dependent spatial distribution of temperature in the brain tissue. Our results suggest that the impact of the temperature increase during our novel ESM paradigm is thermally non-destructive. The presented simulation results match the previously published thermographic measurement and histopathological examination of the stimulated brain tissue and confirm the safety of the novel ESM.

Optimization-Based Contact Fault Alleviation in Deep Brain Stimulation Leads

Deep brain stimulation (DBS) is a neurosurgical treatment in, e.g., Parkinson's Disease. Electrical stimulation in DBS is delivered to a certain target through electrodes implanted into the brain. Recent developments aiming at better stimulation target coverage and lesser side effects have led to an increase in the number of contacts in a DBS lead as well as higher hardware complexity. This paper proposes an optimization-based approach to alleviation of the fault impact on the resulting therapeutical effect in field steering DBS. Faulty contacts could be an issue given recent trends of increasing number of contacts in DBS leads. Hence, a fault detection/alleviation scheme, such as the one proposed in this paper, is necessary ensure resilience in the chronic stimulation. Two alternatives are considered and compared with the stimulation prior to the fault: one using higher amplitudes on the remaining contacts and another with alleviating contacts in the neighborhood of the faulty one. Satisfactory compensation for a faulty contact can be achieved in both ways. However, to designate alleviating contacts, a model-based optimization procedure is necessary. Results suggest that stimulating with more contacts yields configurations that are more robust to contact faults, though with reduced selectivity.

Lead-DBS v2: Towards a comprehensive pipeline for deep brain stimulation imaging

Deep brain stimulation (DBS) is a highly efficacious treatment option for movement disorders and a growing number of other indications are investigated in clinical trials. To ensure optimal treatment outcome, exact electrode placement is required. Moreover, to analyze the relationship between electrode location and clinical results, a precise reconstruction of electrode placement is required, posing specific challenges to the field of neuroimaging. Since 2014 the open source toolbox Lead-DBS is available, which aims at facilitating this process. The tool has since become a popular platform for DBS imaging. With support of a broad community of researchers worldwide, methods have been continuously updated and complemented by new tools for tasks such as multispectral nonlinear registration, structural/functional connectivity analyses, brain shift correction, reconstruction of microelectrode recordings and orientation detection of segmented DBS leads. The rapid development and emergence of these methods in DBS data analysis require us to revisit and revise the pipelines introduced in the original methods publication. Here we demonstrate the updated DBS and connectome pipelines of Lead-DBS using a single patient example with state-of-the-art high-field imaging as well as a retrospective cohort of patients scanned in a typical clinical setting at 1.5T. Imaging data of the 3T example patient is co-registered using five algorithms and nonlinearly warped into template space using ten approaches for comparative purposes. After reconstruction of DBS electrodes (which is possible using three methods and a specific refinement tool), the volume of tissue activated is calculated for two DBS settings using four distinct models and various parameters. Finally, four whole-brain tractography algorithms are applied to the patient's preoperative diffusion MRI data and structural as well as functional connectivity between the stimulation volume and other brain areas are estimated using a total of eight approaches and datasets. In addition, we demonstrate impact of selected preprocessing strategies on the retrospective sample of 51 PD patients. We compare the amount of variance in clinical improvement that can be explained by the computer model depending on the preprocessing method of choice. This work represents a multi-institutional collaborative effort to develop a comprehensive, open source pipeline for DBS imaging and connectomics, which has already empowered several studies, and may facilitate a variety of future studies in the field.

Multi-objective particle swarm optimization for postoperative deep brain stimulation targeting of subthalamic nucleus pathways

OBJECTIVE

The effectiveness of deep brain stimulation (DBS) therapy strongly depends on precise surgical targeting of intracranial leads and on clinical optimization of stimulation settings. Recent advances in surgical targeting, multi-electrode designs, and multi-channel independent current-controlled stimulation are poised to enable finer control in modulating pathways within the brain. However, the large stimulation parameter space enabled by these technologies also poses significant challenges for efficiently identifying the most therapeutic DBS setting for a given patient. Here, we present a computational approach for programming directional DBS leads that is based on a non-convex optimization framework for neural pathway targeting.

APPROACH

The algorithm integrates patient-specific pre-operative 7 T MR imaging, post-operative CT scans, and multi-objective particle swarm optimization (MOPSO) methods using dominance based-criteria and incorporating multiple neural pathways simultaneously. The algorithm was evaluated on eight patient-specific models of subthalamic nucleus (STN) DBS to identify electrode configurations and stimulation amplitudes to optimally activate or avoid six clinically relevant pathways: motor territory of STN, non-motor territory of STN, internal capsule, superior cerebellar peduncle, thalamic fasciculus, and hyperdirect pathway.

MAIN RESULTS

Across the patient-specific models, single-electrode stimulation showed significant correlations across modeled pathways, particularly for motor and non-motor STN efferents. The MOPSO approach was able to identify multi-electrode configurations that achieved improved targeting of motor STN efferents and hyperdirect pathway afferents than that achieved by any single-electrode monopolar setting at equivalent power levels.

SIGNIFICANCE

These results suggest that pathway targeting with patient-specific model-based optimization algorithms can efficiently identify non-trivial electrode configurations for enhancing activation of clinically relevant pathways. However, the results also indicate that inter-pathway correlations can limit selectivity for certain pathways even with directional DBS leads.

Clinical deep brain stimulation strategies for orientation-selective pathway activation

OBJECTIVE

This study investigated stimulation strategies to increase the selectivity of activating axonal pathways within the brain based on their orientations relative to clinical deep brain stimulation (DBS) lead implants.

APPROACH

Previous work has shown how varying electrode shape and controlling the primary electric field direction through preclinical electrode arrays can produce orientation-selective axonal stimulation. Here, we significantly extend those results using computational models to evaluate the degree to which clinical DBS leads can direct stimulus-induced electric fields and generate orientation-selective activation of fiber pathways in the brain. Orientation-selective pulse paradigms were evaluated in conceptual models and in patient-specific models of subthalamic nucleus (STN)-DBS for treating Parkinson's disease.

MAIN RESULTS

Single-contact monopolar or two-contact bipolar stimulation through clinical DBS leads with cylindrical electrodes primarily activated axons orientated parallel to the lead. Conversely, multi-contact monopolar stimulation with a cathode-leading pulse waveform selectively activated axons perpendicular to the DBS lead. Clinical DBS leads with segmented rows of electrodes and a single current source provided additional angular resolution for activating axons oriented 0°, ±22.5°, ±45°, ±67.5°, or 90° relative to the lead shaft. Employing multiple independent current sources to deliver unequal amounts of current through these leads further increased the angular resolution of activation relative to the lead shaft. The patient-specific models indicated that multi-contact cathode configurations, which are rarely used in clinical practice, could increase activation of the hyperdirect pathway collaterals projecting into STN (a putative therapeutic target), while minimizing direct activation of the corticospinal tract of internal capsule, which can elicit sensorimotor side-effects when stimulated.

SIGNIFICANCE

When combined with patient-specific tissue anisotropy and patient-specific anatomical morphologies of neural pathways responsible for therapy and side effects, orientation-selective DBS approaches show potential to significantly improve clinical outcomes of DBS therapy for a range of existing and investigational clinical indications.

Electrical stimulation of the midbrain excites the auditory cortex asymmetrically

BACKGROUND

Auditory midbrain implant users cannot achieve open speech perception and have limited frequency resolution. It remains unclear whether the spread of excitation contributes to this issue and how much it can be compensated by current-focusing, which is an effective approach in cochlear implants.

OBJECTIVE

The present study examined the spread of excitation in the cortex elicited by electric midbrain stimulation. We further tested whether current-focusing via bipolar and tripolar stimulation is effective with electric midbrain stimulation and whether these modes hold any advantage over monopolar stimulation also in conditions when the stimulation electrodes are in direct contact with the target tissue.

METHODS

Using penetrating multielectrode arrays, we recorded cortical population responses to single pulse electric midbrain stimulation in 10 ketamine/xylazine anesthetized mice. We compared monopolar, bipolar, and tripolar stimulation configurations with regard to the spread of excitation and the characteristic frequency difference between the stimulation/recording electrodes.

RESULTS

The cortical responses were distributed asymmetrically around the characteristic frequency of the stimulated midbrain region with a strong activation in regions tuned up to one octave higher. We found no significant differences between monopolar, bipolar, and tripolar stimulation in threshold, evoked firing rate, or dynamic range.

CONCLUSION

The cortical responses to electric midbrain stimulation are biased towards higher tonotopic frequencies. Current-focusing is not effective in direct contact electrical stimulation. Electrode maps should account for the asymmetrical spread of excitation when fitting auditory midbrain implants by shifting the frequency-bands downward and stimulating as dorsally as possible.

Quantitative theory of deep brain stimulation of the subthalamic nucleus for the suppression of pathological rhythms in Parkinson's disease

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is modeled to explore the mechanisms of this effective, but poorly understood, treatment for motor symptoms of drug-refractory Parkinson's disease and dystonia. First, a neural field model of the corticothalamic-basal ganglia (CTBG) system is developed that reproduces key clinical features of Parkinson's disease, including its characteristic 4-8 Hz and 13-30 Hz electrophysiological signatures. Deep brain stimulation of the STN is then modeled and shown to suppress the pathological 13-30 Hz (beta) activity for physiologically realistic and optimized stimulus parameters. This supports the idea that suppression of abnormally coherent activity in the CTBG system is a major factor in DBS therapy for Parkinson's disease, by permitting normal dynamics to resume. At high stimulus intensities, nonlinear effects in the target population mediate wave-wave interactions between resonant beta activity and the stimulus pulse train, leading to complex spectral structure that shows remarkable similarity to that seen in steady-state evoked potential experiments.

Adaptive Estimation of the Neural Activation Extent in Computational Volume Conductor Models of Deep Brain Stimulation

OBJECTIVE

The aim of this study is to propose an adaptive scheme embedded into an open-source environment for the estimation of the neural activation extent during deep brain stimulation and to investigate the feasibility of approximating the neural activation extent by thresholds of the field solution.

METHODS

Open-source solutions for solving the field equation in volume conductor models of deep brain stimulation and computing the neural activation are embedded into a Python package to estimate the neural activation dependent on the dielectric tissue properties and axon parameters by employing a spatially adaptive scheme. Feasibility of the approximation of the neural activation extent by field thresholds is investigated to further reduce the computational expense.

RESULTS

The varying extents of neural activation for different patient-specific dielectric properties were estimated with the adaptive scheme. The results revealed the strong influence of the dielectric properties of the encapsulation layer in the acute and chronic phase after surgery. The computational time required to determine the neural activation extent in each studied model case was substantially reduced.

CONCLUSION

The neural activation extent is altered by patient-specific parameters. Threshold values of the electric potential and electric field norm facilitate a computationally efficient method to estimate the neural activation extent.

SIGNIFICANCE

The presented adaptive scheme is able to robustly determine neural activation extents and field threshold estimates for varying dielectric tissue properties and axon diameters while substantially reducing the computational expense.

Subthalamic deep brain stimulation sweet spots and hyperdirect cortical connectivity in Parkinson's disease

OBJECTIVES

Firstly, to identify subthalamic region stimulation clusters that predict maximum improvement in rigidity, bradykinesia and tremor, or emergence of side-effects; and secondly, to map-out the cortical fingerprint, mediated by the hyperdirect pathways which predict maximum efficacy.

METHODS

High angular resolution diffusion imaging in twenty patients with advanced Parkinson's disease was acquired prior to bilateral subthalamic nucleus deep brain stimulation. All contacts were screened one-year from surgery for efficacy and side-effects at different amplitudes. Voxel-based statistical analysis of volumes of tissue activated models was used to identify significant treatment clusters. Probabilistic tractography was employed to identify cortical connectivity patterns associated with treatment efficacy.

RESULTS

All patients responded well to treatment (46% mean improvement off medication UPDRS-III [p < 0.0001]) without significant adverse events. Cluster corresponding to maximum improvement in tremor was in the posterior, superior and lateral portion of the nucleus. Clusters corresponding to improvement in bradykinesia and rigidity were nearer the superior border in a further medial and posterior location. The rigidity cluster extended beyond the superior border to the area of the zona incerta and Forel-H₂ field. When the clusters where averaged, the coordinates of the area with maximum overall efficacy was X = -10(-9.5), Y = -13(-1) and Z = -7(-3) in MNI(AC-PC) space. Cortical connectivity to primary motor area was predictive of higher improvement in tremor; whilst that to supplementary motor area was predictive of improvement in bradykinesia and rigidity; and connectivity to prefrontal cortex was predictive of improvement in rigidity.

INTERPRETATION

These findings support the presence of overlapping stimulation sites within the subthalamic nucleus and its superior border, with different cortical connectivity patterns, associated with maximum improvement in tremor, rigidity and bradykinesia.

A comparative study of approaches to compute the field distribution of deep brain stimulation in the Hemiparkinson rat model

Computational modeling of the stimulating field distribution during Deep Brain Stimulation provides an opportunity to advance our knowledge of this neurosurgical therapy for Parkinson's disease. There exist several approaches to model the target region for Deep Brain Stimulation in Hemi-parkinson Rats with volume conductor models. We have described and compared the normalized mapping approach as well as the modeling with three-dimensional structures, which include curvilinear coordinates to assure an anatomically realistic conductivity tensor orientation.

An implementation for the simulation of cells on micro-post arrays

The mechanical interaction between cells and their underlying substrates is important in understanding the processes that take place at an interface between biological tissue and the surface of implants. There have been numerous studies that examine these interactions both by experimental and numerical modeling. The bio-chemo-mechanical model for cell contractility by Deshpande et al. [1] has numerous applications and advantages. This work shows a way to implement this model in COMSOL MULTIPHYSICS® so it can be easily modified or extended. This will allow us in a next step to couple the differential system with additional external stimuli.

Particle swarm optimization for programming deep brain stimulation arrays

OBJECTIVE

Deep brain stimulation (DBS) therapy relies on both precise neurosurgical targeting and systematic optimization of stimulation settings to achieve beneficial clinical outcomes. One recent advance to improve targeting is the development of DBS arrays (DBSAs) with electrodes segmented both along and around the DBS lead. However, increasing the number of independent electrodes creates the logistical challenge of optimizing stimulation parameters efficiently.

APPROACH

Solving such complex problems with multiple solutions and objectives is well known to occur in biology, in which complex collective behaviors emerge out of swarms of individual organisms engaged in learning through social interactions. Here, we developed a particle swarm optimization (PSO) algorithm to program DBSAs using a swarm of individual particles representing electrode configurations and stimulation amplitudes. Using a finite element model of motor thalamic DBS, we demonstrate how the PSO algorithm can efficiently optimize a multi-objective function that maximizes predictions of axonal activation in regions of interest (ROI, cerebellar-receiving area of motor thalamus), minimizes predictions of axonal activation in regions of avoidance (ROA, somatosensory thalamus), and minimizes power consumption.

MAIN RESULTS

The algorithm solved the multi-objective problem by producing a Pareto front. ROI and ROA activation predictions were consistent across swarms (<1% median discrepancy in axon activation). The algorithm was able to accommodate for (1) lead displacement (1 mm) with relatively small ROI (⩽9.2%) and ROA (⩽1%) activation changes, irrespective of shift direction; (2) reduction in maximum per-electrode current (by 50% and 80%) with ROI activation decreasing by 5.6% and 16%, respectively; and (3) disabling electrodes (n  =  3 and 12) with ROI activation reduction by 1.8% and 14%, respectively. Additionally, comparison between PSO predictions and multi-compartment axon model simulations showed discrepancies of  <1% between approaches.

SIGNIFICANCE

The PSO algorithm provides a computationally efficient way to program DBS systems especially those with higher electrode counts.

Investigation into Deep Brain Stimulation Lead Designs: A Patient-Specific Simulation Study

New deep brain stimulation (DBS) electrode designs offer operation in voltage and current mode and capability to steer the electric field (EF). The aim of the study was to compare the EF distributions of four DBS leads at equivalent amplitudes (3 V and 3.4 mA). Finite element method (FEM) simulations (n = 38) around cylindrical contacts (leads 3389, 6148) or equivalent contact configurations (leads 6180, SureStim1) were performed using homogeneous and patient-specific (heterogeneous) brain tissue models. Steering effects of 6180 and SureStim1 were compared with symmetric stimulation fields. To make relative comparisons between simulations, an EF isolevel of 0.2 V/mm was chosen based on neuron model simulations (n = 832) applied before EF visualization and comparisons. The simulations show that the EF distribution is largely influenced by the heterogeneity of the tissue, and the operating mode. Equivalent contact configurations result in similar EF distributions. In steering configurations, larger EF volumes were achieved in current mode using equivalent amplitudes. The methodology was demonstrated in a patient-specific simulation around the zona incerta and a "virtual" ventral intermediate nucleus target. In conclusion, lead design differences are enhanced when using patient-specific tissue models and current stimulation mode.

Model-Based Comparison of Deep Brain Stimulation Array Functionality with Varying Number of Radial Electrodes and Machine Learning Feature Sets

Deep brain stimulation (DBS) leads with radially distributed electrodes have potential to improve clinical outcomes through more selective targeting of pathways and networks within the brain. However, increasing the number of electrodes on clinical DBS leads by replacing conventional cylindrical shell electrodes with radially distributed electrodes raises practical design and stimulation programming challenges. We used computational modeling to investigate: (1) how the number of radial electrodes impact the ability to steer, shift, and sculpt a region of neural activation (RoA), and (2) which RoA features are best used in combination with machine learning classifiers to predict programming settings to target a particular area near the lead. Stimulation configurations were modeled using 27 lead designs with one to nine radially distributed electrodes. The computational modeling framework consisted of a three-dimensional finite element tissue conductance model in combination with a multi-compartment biophysical axon model. For each lead design, two-dimensional threshold-dependent RoAs were calculated from the computational modeling results. The models showed more radial electrodes enabled finer resolution RoA steering; however, stimulation amplitude, and therefore spatial extent of the RoA, was limited by charge injection and charge storage capacity constraints due to the small electrode surface area for leads with more than four radially distributed electrodes. RoA shifting resolution was improved by the addition of radial electrodes when using uniform multi-cathode stimulation, but non-uniform multi-cathode stimulation produced equivalent or better resolution shifting without increasing the number of radial electrodes. Robust machine learning classification of 15 monopolar stimulation configurations was achieved using as few as three geometric features describing a RoA. The results of this study indicate that, for a clinical-scale DBS lead, more than four radial electrodes minimally improved in the ability to steer, shift, and sculpt axonal activation around a DBS lead and a simple feature set consisting of the RoA center of mass and orientation enabled robust machine learning classification. These results provide important design constraints for future development of high-density DBS arrays.

Technical aspects of neurostimulation: Focus on equipment, electric field modeling, and stimulation protocols

Neuromodulation is a field of science, medicine, and bioengineering that encompasses implantable and non-implantable technologies for the purpose of improving quality of life and functioning of humans. Brain neuromodulation involves different neurostimulation techniques: transcranial magnetic stimulation (TMS), transcranial direct current stimulation (tDCS), vagus nerve stimulation (VNS), and deep brain stimulation (DBS), which are being used both to study their effects on cognitive brain functions and to treat neuropsychiatric disorders. The mechanisms of action of neurostimulation remain incompletely understood. Insight into the technical basis of neurostimulation might be a first step towards a more profound understanding of these mechanisms, which might lead to improved clinical outcome and therapeutic potential. This review provides an overview of the technical basis of neurostimulation focusing on the equipment, the present understanding of induced electric fields, and the stimulation protocols. The review is written from a technical perspective aimed at supporting the use of neurostimulation in clinical practice.

A latent force model for describing electric propagation in deep brain stimulation: a simulation study

Deep brain stimulation (DBS) is a neurosurgical method used to treat symptoms of movement disorders by implanting electrodes in deep brain areas. Often, the DBS modeling approaches found in the literature assume a quasi-static approximation, and discard any dynamic behavior. Nevertheless, in a real DBS system the stimulus corresponds to a wave that changes as a function of time. It is clear that DBS demands an approach that takes into account the time-varying behavior of the input stimulus. In this work, we present a novel latent force model for describing the dynamic electric propagation occurred during DBS. The performance of the proposed model was studied by simulations under different conditions. The results show that our approach is able to take into account the time variations of the source and the produced field. Moreover, by restricting our model it is possible to obtain solutions for electrostatic formulations, here experimental results were compared with the finite element method. Additionally, our approach allows a solution to the inverse problem, which is a valuable clinical application allowing the appropriate tuning of the DBS device by the expert physician.

Modeling sequence and quasi-uniform assumption in computational neurostimulation

Computational neurostimulation aims to develop mathematical constructs that link the application of neuromodulation with changes in behavior and cognition. This process is critical but daunting for technical challenges and scientific unknowns. The overarching goal of this review is to address how this complex task can be made tractable. We describe a framework of sequential modeling steps to achieve this: (1) current flow models, (2) cell polarization models, (3) network and information processing models, and (4) models of the neuroscientific correlates of behavior. Each step is explained with a specific emphasis on the assumptions underpinning underlying sequential implementation. We explain the further implementation of the quasi-uniform assumption to overcome technical limitations and unknowns. We specifically focus on examples in electrical stimulation, such as transcranial direct current stimulation. Our approach and conclusions are broadly applied to immediate and ongoing efforts to deploy computational neurostimulation.

Computational modeling of neurostimulation in brain diseases

Neurostimulation as a therapeutic tool has been developed and used for a range of different diseases such as Parkinson's disease, epilepsy, and migraine. However, it is not known why the efficacy of the stimulation varies dramatically across patients or why some patients suffer from severe side effects. This is largely due to the lack of mechanistic understanding of neurostimulation. Hence, theoretical computational approaches to address this issue are in demand. This chapter provides a review of mechanistic computational modeling of brain stimulation. In particular, we will focus on brain diseases, where mechanistic models (e.g., neural population models or detailed neuronal models) have been used to bridge the gap between cellular-level processes of affected neural circuits and the symptomatic expression of disease dynamics. We show how such models have been, and can be, used to investigate the effects of neurostimulation in the diseased brain. We argue that these models are crucial for the mechanistic understanding of the effect of stimulation, allowing for a rational design of stimulation protocols. Based on mechanistic models, we argue that the development of closed-loop stimulation is essential in order to avoid inference with healthy ongoing brain activity. Furthermore, patient-specific data, such as neuroanatomic information and connectivity profiles obtainable from neuroimaging, can be readily incorporated to address the clinical issue of variability in efficacy between subjects. We conclude that mechanistic computational models can and should play a key role in the rational design of effective, fully integrated, patient-specific therapeutic brain stimulation.

Subject-specific computational modeling of DBS in the PPTg area

Deep brain stimulation (DBS) in the pedunculopontine tegmental nucleus (PPTg) has been proposed to alleviate medically intractable gait difficulties associated with Parkinson's disease. Clinical trials have shown somewhat variable outcomes, stemming in part from surgical targeting variability, modulating fiber pathways implicated in side effects, and a general lack of mechanistic understanding of DBS in this brain region. Subject-specific computational models of DBS are a promising tool to investigate the underlying therapy and side effects. In this study, a parkinsonian rhesus macaque was implanted unilaterally with an 8-contact DBS lead in the PPTg region. Fiber tracts adjacent to PPTg, including the oculomotor nerve, central tegmental tract, and superior cerebellar peduncle, were reconstructed from a combination of pre-implant 7T MRI, post-implant CT, and post-mortem histology. These structures were populated with axon models and coupled with a finite element model simulating the voltage distribution in the surrounding neural tissue during stimulation. This study introduces two empirical approaches to evaluate model parameters. First, incremental monopolar cathodic stimulation (20 Hz, 90 μs pulse width) was evaluated for each electrode, during which a right eyelid flutter was observed at the proximal four contacts (−1.0 to −1.4 mA). These current amplitudes followed closely with model predicted activation of the oculomotor nerve when assuming an anisotropic conduction medium. Second, PET imaging was collected OFF-DBS and twice during DBS (two different contacts), which supported the model predicted activation of the central tegmental tract and superior cerebellar peduncle. Together, subject-specific models provide a framework to more precisely predict pathways modulated by DBS.

High-definition transcranial direct current stimulation induces both acute and persistent changes in broadband cortical synchronization: a simultaneous tDCS-EEG study

The goal of this study was to develop methods for simultaneously acquiring electrophysiological data during high-definition transcranial direct current stimulation (tDCS) using high-resolution electroencephalography (EEG). Previous studies have pointed to the after-effects of tDCS on both motor and cognitive performance, and there appears to be potential for using tDCS in a variety of clinical applications. However, little is known about the real-time effects of tDCS on rhythmic cortical activity in humans due to the technical challenges of simultaneously obtaining electrophysiological data during ongoing stimulation. Furthermore, the mechanisms of action of tDCS in humans are not well understood. We have conducted a simultaneous tDCS-EEG study in a group of healthy human subjects. Significant acute and persistent changes in spontaneous neural activity and event-related synchronization (ERS) were observed during and after the application of high-definition tDCS over the left sensorimotor cortex. Both anodal and cathodal stimulation resulted in acute global changes in broadband cortical activity which were significantly different than the changes observed in response to sham stimulation. For the group of eight subjects studied, broadband individual changes in spontaneous activity during stimulation were apparent both locally and globally. In addition, we found that high-definition tDCS of the left sensorimotor cortex can induce significant ipsilateral and contralateral changes in event-related desynchronization and ERS during motor imagination following the end of the stimulation period. Overall, our results demonstrate the feasibility of acquiring high-resolution EEG during high-definition tDCS and provide evidence that tDCS in humans directly modulates rhythmic cortical synchronization during and after its administration.

Coordinated reset stimulation in a large-scale model of the STN-GPe circuit

Synchronization of populations of neurons is a hallmark of several brain diseases. Coordinated reset (CR) stimulation is a model-based stimulation technique which specifically counteracts abnormal synchrony by desynchronization. Electrical CR stimulation, e.g., for the treatment of Parkinson's disease (PD), is administered via depth electrodes. In order to get a deeper understanding of this technique, we extended the top-down approach of previous studies and constructed a large-scale computational model of the respective brain areas. Furthermore, we took into account the spatial anatomical properties of the simulated brain structures and incorporated a detailed numerical representation of 2 · 10⁴ simulated neurons. We simulated the subthalamic nucleus (STN) and the globus pallidus externus (GPe). Connections within the STN were governed by spike-timing dependent plasticity (STDP). In this way, we modeled the physiological and pathological activity of the considered brain structures. In particular, we investigated how plasticity could be exploited and how the model could be shifted from strongly synchronized (pathological) activity to strongly desynchronized (healthy) activity of the neuronal populations via CR stimulation of the STN neurons. Furthermore, we investigated the impact of specific stimulation parameters especially the electrode position on the stimulation outcome. Our model provides a step forward toward a biophysically realistic model of the brain areas relevant to the emergence of pathological neuronal activity in PD. Furthermore, our model constitutes a test bench for the optimization of both stimulation parameters and novel electrode geometries for efficient CR stimulation.

Relationship between neural activation and electric field distribution during deep brain stimulation

Models and simulations are commonly used to study deep brain stimulation (DBS). Simulated stimulation fields are often defined and visualized by electric field isolevels or volumes of tissue activated (VTA). The aim of the present study was to evaluate the relationship between stimulation field strength as defined by the electric potential V, the electric field E, and the divergence of the electric field ∇(2) V, and neural activation. Axon cable models were developed and coupled to finite-element DBS models in three-dimensional (3-D). Field thresholds ( VT , ET, and ∇(2) VT ) were derived at the location of activation for various stimulation amplitudes (1 to 5 V), pulse widths (30 to 120 μs), and axon diameters (2.0 to 7.5 μm). Results showed that thresholds for VT and ∇(2) VT were highly dependent on the stimulation amplitude while ET were approximately independent of the amplitude for large axons. The activation field strength thresholds presented in this study may be used in future studies to approximate the VTA during model-based investigations of DBS without the need of computational axon models.