Individually optimized multi-channel tDCS for targeting somatosensory cortex

Comparing the performance of beamformer algorithms in estimating orientations of neural sources

The efficacy of transcranial electric stimulation (tES) to effectively modulate neuronal activity depends critically on the spatial orientation of the targeted neuronal population. Therefore, precise estimation of target orientation is of utmost importance. Different beamforming algorithms provide orientation estimates; however, a systematic analysis of their performance is still lacking. For fixed brain locations, EEG and MEG data from sources with randomized orientations were simulated. The orientation was then estimated (1) with an EEG and (2) with a combined EEG-MEG approach. Three commonly used beamformer algorithms were evaluated with respect to their abilities to estimate the correct orientation: Unit-Gain (UG), Unit-Noise-Gain (UNG), and Array-Gain (AG) beamformer. Performance depends on the signal-to-noise ratios for the modalities and on the chosen beamformer. Overall, the UNG and AG beamformers appear as the most reliable. With increasing noise, the UG estimate converges to a vector determined by the leadfield, thus leading to insufficient orientation estimates.

Individual and group-level optimization of electric field in deep brain region during multichannel transcranial electrical stimulation

Electrode montage optimization for transcranial electric stimulation (tES) is a challenging topic for targeting a specific brain region. Targeting the deep brain region is difficult due to tissue inhomogeneity, resulting in complex current flow. In this study, a simplified protocol for montage optimization is proposed for multichannel tES (mc-tES). The purpose of this study was to reduce the computational cost for mc-tES optimization and to evaluate the mc-tES for deep brain regions. Optimization was performed using a simplified protocol for montages under safety constraints with 20 anatomical head models. The optimization procedure is simplified using the surface EF of the deep brain target region, considering its small volume and non-concentric distribution of the electrodes. Our proposal demonstrated that the computational cost was reduced by >90%. A total of six-ten electrodes were necessary for robust EF in the target region. The optimization with surface EF is comparable to or marginally better than using conventional volumetric EF for deep brain tissues. An electrode montage with a mean injection current amplitude derived from individual analysis was demonstrated to be useful for targeting the deep region at the group level. The optimized montage and injection current were derived at the group level. Our proposal at individual and group levels showed great potential for clinical application.

Lattice layout and optimizer effect analysis for generating optimal transcranial electrical stimulation (tES) montages through the metaheuristic L1L1 method

Introduction

This study focuses on broadening the applicability of the metaheuristic L1-norm fitted and penalized (L1L1) optimization method in finding a current pattern for multichannel transcranial electrical stimulation (tES). The metaheuristic L1L1 optimization framework defines the tES montage via linear programming by maximizing or minimizing an objective function with respect to a pair of hyperparameters.

Methods

In this study, we explore the computational performance and reliability of different optimization packages, algorithms, and search methods in combination with the L1L1 method. The solvers from Matlab R2020b, MOSEK 9.0, Gurobi Optimizer, CVX's SeDuMi 1.3.5, and SDPT3 4.0 were employed to produce feasible results through different linear programming techniques, including Interior-Point (IP), Primal-Simplex (PS), and Dual-Simplex (DS) methods. To solve the metaheuristic optimization task of L1L1, we implement an exhaustive and recursive search along with a well-known heuristic direct search as a reference algorithm.

Results

Based on our results, and the given optimization task, Gurobi's IP was, overall, the preferable choice among Interior-Point while MOSEK's PS and DS packages were in the case of Simplex methods. These methods provided substantial computational time efficiency for solving the L1L1 method regardless of the applied search method.

Discussion

While the best-performing solvers show that the L1L1 method is suitable for maximizing either focality and intensity, a few of these solvers could not find a bipolar configuration. Part of the discrepancies between these methods can be explained by a different sensitivity with respect to parameter variation or the resolution of the lattice provided.

Normative tDCS over V5 and FEF reveals practice-induced modulation of extraretinal smooth pursuit mechanisms, but no specific stimulation effect

The neural networks subserving smooth pursuit eye movements (SPEM) provide an ideal model for investigating the interaction of sensory processing and motor control during ongoing movements. To better understand core plasticity aspects of sensorimotor processing for SPEM, normative sham, anodal or cathodal transcranial direct current stimulation (tDCS) was applied over visual area V5 and frontal eye fields (FEF) in sixty healthy participants. The identical within-subject paradigm was used to assess SPEM modulations by practice. While no specific tDCS effects were revealed, within- and between-session practice effects indicate plasticity of top-down extraretinal mechanisms that mainly affect SPEM in the absence of visual input and during SPEM initiation. To explore the potential of tDCS effects, individual electric field simulations were computed based on calibrated finite element head models and individual functional localization of V5 and FEF location (using functional MRI) and orientation (using combined EEG/MEG) was conducted. Simulations revealed only limited electric field target intensities induced by the applied normative tDCS montages but indicate the potential efficacy of personalized tDCS for the modulation of SPEM. In sum, results indicate the potential susceptibility of extraretinal SPEM control to targeted external neuromodulation (e.g., personalized tDCS) and intrinsic learning protocols.

Monopolar tDCS might affect brainstem reflexes: A computational and neurophysiological study

OBJECTIVE

To assess whether monopolar multi-electrode transcranial direct current stimulation (tDCS) montages might selectively affect deep brain structures through computational predictions and neurophysiological assessment.

METHODS

Electric field distribution in deep brain structures (i.e., thalamus and midbrain) were estimated through computational models simulating tDCS with two monopolar and two monopolar multi-electrode montages. Monopolar multi-electrode tDCS was then applied to healthy subject, and effects on pontine and medullary circuitries was evaluated studying changes in blink reflex (BR) and masseter inhibitory reflex (MIR).

RESULTS

Computational results suggest that tDCS with monopolar multi-electrode montages might induce electric field intensities in deep brain structure comparable to those in grey matter, while neurophysiological results disclosed that BR and MIR were selectively modulated by tDCS only when cathode was placed over the right deltoid.

CONCLUSIONS

Multi-electrode tDCS (anodes over motor cortices, cathode over right deltoid) could induce significant electric fields in the thalamus and midbrain, and selectively affect brainstem neural circuits.

SIGNIFICANCE

Multi-electrode tDCS (anodes over motor cortices, cathode over right deltoid) might be further explored to affect brainstem activity, also in the context of non-invasive deep brain stimulation.

Personalized alpha-tACS targeting left posterior parietal cortex modulates visuo-spatial attention and posterior evoked EEG activity

BACKGROUND

Covert visuo-spatial attention is marked by the anticipatory lateralization of neuronal alpha activity in the posterior parietal cortex. Previous applications of transcranial alternating current stimulation (tACS) at the alpha frequency, however, were inconclusive regarding the causal contribution of oscillatory activity during visuo-spatial attention.

OBJECTIVE

Attentional shifts of behavior and electroencephalography (EEG) after-effects were assessed in a cued visuo-spatial attention paradigm. We hypothesized that parietal alpha-tACS shifts attention relative to the ipsilateral visual hemifield. Furthermore, we assumed that modulations of behavior and neurophysiology are related to individual electric field simulations.

METHODS

We applied personalized tACS at alpha and gamma frequencies to elucidate the role of oscillatory neuronal activity for visuo-spatial attention. Personalized tACS montages were algorithmically optimized to target individual left and right parietal regions that were defined by an EEG localizer.

RESULTS

Behavioral performance in the left hemifield was specifically increased by alpha-tACS compared to gamma-tACS targeting the left parietal cortex. This hemisphere-specific effect was observed despite the symmetry of simulated electric fields. In addition, visual event-related potential (ERP) amplitudes showed a reduced lateralization over posterior sites induced by left alpha-tACS. Neuronal sources of this effect were localized in the left premotor cortex. Interestingly, accuracy modulations induced by left parietal alpha-tACS were directly related to electric field magnitudes in the left premotor cortex.

CONCLUSION

Overall, results corroborate the notion that alpha lateralization plays a causal role in covert visuo-spatial attention and indicate an increased susceptibility of parietal and premotor brain regions of the left dorsal attention network to subtle tACS-neuromodulation.

An Operational Approach for Optimizing Transcranial Direct Current Stimulation

Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation technique that has been utilized for treating brain disorders and improving cognitive function. In order to achieve targeted tDCS, many optimization methods of montages and electric currents have been proposed. However, these methods have some limitations. Most of them were proposed for single-objective optimization (focality or intensity) and have no constrain with the number of electrodes (Most devices only have less than 8 electrodes currently). In this study, we proposed an operational optimization approach for well-targeted tDCS, which aims to optimize for two objectives of electric field (EF) intensity and focality with constraints on the number of electrodes. Compared with traditional tDCS in our cohort (10 subjects), our method significantly improves the EF focality. When compared to commonly used 4×1 high-definition tDCS (HD-tDCS), our method can achieve higher EF intensity in the target region with less than 8 electrodes. Our method can balance the two objectives of EF and shorten optimization time, which is convenient for practical application.

Noninvasive Brain Stimulation Techniques for Treatment-Resistant Depression: Transcranial Magnetic Stimulation and Transcranial Direct Current Stimulation

Transcranial magnetic stimulation is a safe, effective, and well-tolerated intervention for depression; it is currently approved for treatment-resistant depression. This article summarizes the mechanism of action, evidence of clinical efficacy, and the clinical aspects of this intervention, including patient evaluation, stimulation parameters selection, and safety considerations. Transcranial direct current stimulation is another neuromodulation treatment for depression; although promising, the technique is not currently approved for clinical use in the United States. The final section outlines the open challenges and future directions of the field.

Review of individualized current flow modeling studies for transcranial electrical stimulation

There is substantial intersubject variability of behavioral and neurophysiological responses to transcranial electrical stimulation (tES), which represents one of the most important limitations of tES. Many tES protocols utilize a fixed experimental parameter set disregarding individual anatomical and physiological properties. This one-size-fits-all approach might be one reason for the observed interindividual response variability. Simulation of current flow applying head models based on available anatomical data can help to individualize stimulation parameters and contribute to the understanding of the causes of this response variability. Current flow modeling can be used to retrospectively investigate the characteristics of tES effectivity. Previous studies examined, for example, the impact of skull defects and lesions on the modulation of current flow and demonstrated effective stimulation intensities in different age groups. Furthermore, uncertainty analysis of electrical conductivities in current flow modeling indicated the most influential tissue compartments. Current flow modeling, when used in prospective study planning, can potentially guide stimulation configurations resulting in individually effective tES. Specifically, current flow modeling using individual or matched head models can be employed by clinicians and scientists to, for example, plan dosage in tES protocols for individuals or groups of participants. We review studies that show a relationship between the presence of behavioral/neurophysiological responses and features derived from individualized current flow models. We highlight the potential benefits of individualized current flow modeling.

Bioelectromagnetism in Human Brain Research: New Applications, New Questions

Bioelectromagnetism has contributed some of the most commonly used techniques to human neuroscience such as magnetoencephalography (MEG), electroencephalography (EEG), transcranial magnetic stimulation (TMS), and transcranial electric stimulation (TES). The considerable differences in their technical design and practical use give rise to the impression that these are quite different techniques altogether. Here, we review, discuss and illustrate the fundamental principle of Helmholtz reciprocity that provides a common ground for all four techniques. We show that, more than 150 years after its discovery by Helmholtz in 1853, reciprocity is important to appreciate the strengths and limitations of these four classical tools in neuroscience. We build this case by explaining the concept of Helmholtz reciprocity, presenting a methodological account of this principle for all four methods and, finally, by illustrating its application in practical clinical studies.

Can individually targeted and optimized multi-channel tDCS outperform standard bipolar tDCS in stimulating the primary somatosensory cortex?

BACKGROUND

Transcranial direct current stimulation (tDCS) has emerged as a non-invasive neuro-modulation technique. Most studies show that anodal tDCS increases cortical excitability, however, with variable outcomes. Previously, we have shown in computer simulations that our multi-channel tDCS (mc-tDCS) approach, the distributed constrained maximum intensity (D-CMI) method can potentially lead to better controlled tDCS results due to the improved directionality of the injected current at the target side for individually optimized D-CMI montages.

OBJECTIVE

In this study, we test the application of the D-CMI approach in an experimental study to stimulate the somatosensory P20/N20 target source in Brodmann area 3b and compare it with standard bipolar tDCS and sham conditions.

METHODS

We applied anodal D-CMI, the standard bipolar and D-CMI based Sham tDCS for 10 min to target the 20 ms post-stimulus somatosensory P20/N20 target brain source in Brodmann area 3b reconstructed using combined magnetoencephalography (MEG) and electroencephalography (EEG) source analysis in realistic head models with calibrated skull conductivity in a group-study with 13 subjects. Finger-stimulated somatosensory evoked fields (SEF) were recorded and the component at 20 ms post-stimulus (M20) was analyzed before and after the application of the three tDCS conditions in order to read out the stimulation effect on Brodmann area 3b.

RESULTS

Analysis of the finger stimulated SEF M20 peak before (baseline) and after tDCS shows a significant increase in source amplitude in Brodmann area 3b for D-CMI (6-16 min after tDCS), while no significant effects are found for standard bipolar (6-16 min after tDCS) and sham (6-16 min after tDCS) stimulation conditions. For the later time courses (16-26 and 27-37 min post-stimulation), we found a significant decrease in M20 peak source amplitude for standard bipolar and sham tDCS, while there was no effect for D-CMI.

CONCLUSION

Our results indicate that targeted and optimized, and thereby highly individualized, mc-tDCS can outperform standard bipolar stimulation and lead to better control over stimulation outcomes with, however, a considerable amount of additional work compared to standard bipolar tDCS.

L1-norm vs. L2-norm fitting in optimizing focal multi-channel tES stimulation: linear and semidefinite programming vs. weighted least squares

BACKGROUND AND OBJECTIVE

This study focuses on Multi-Channel Transcranial Electrical Stimulation, a non-invasive brain method for stimulating neuronal activity under the influence of low-intensity currents. We introduce a mathematical formulation for finding a current pattern that optimizes an L1-norm fit between a given focal target distribution and volumetric current density inside the brain. L1-norm is well-known to favor well-localized or sparse distributions compared to L2-norm (least-squares) fitted estimates.

METHODS

We present a linear programming approach that performs L1-norm fitting and penalization of the current pattern (L1L1) to control the number of non-zero currents. The optimizer filters a large set of candidate solutions using a two-stage metaheuristic search from a pre-filtered set of candidates.

RESULTS

The numerical simulation results obtained with both 8- and 20-channel electrode montages suggest that our hypothesis on the benefits of L1-norm data fitting is valid. Compared to an L1-norm regularized L2-norm fitting (L1L2) via semidefinite programming and weighted Tikhonov least-squares method (TLS), the L1L1 results were overall preferable for maximizing the focused current density at the target position, and the ratio between focused and nuisance current magnitudes.

CONCLUSIONS

We propose the metaheuristic L1L1 optimization approach as a potential technique to obtain a well-localized stimulus with a controllable magnitude at a given target position. L1L1 finds a current pattern with a steep contrast between the anodal and cathodal electrodes while suppressing the nuisance currents in the brain, hence, providing a potential alternative to modulate the effects of the stimulation, e.g., the sensation experienced by the subject.

Directionality of the injected current targeting the P20/N20 source determines the efficacy of 140 Hz transcranial alternating current stimulation (tACS)-induced aftereffects in the somatosensory cortex

Interindividual anatomical differences in the human cortex can lead to suboptimal current directions and may result in response variability of transcranial electrical stimulation methods. These differences in brain anatomy require individualized electrode stimulation montages to induce an optimal current density in the targeted area of each individual subject. We aimed to explore the possible modulatory effects of 140 Hz transcranial alternating current stimulation (tACS) on the somatosensory cortex using personalized multi-electrode stimulation montages. In two randomized experiments using either tactile finger or median nerve stimulation, we measured by evoked potentials the plasticity aftereffects and oscillatory power changes after 140 Hz tACS at 1.0 mA as compared to sham stimulation (n = 17, male = 9). We found a decrease in the power of oscillatory mu-rhythms during and immediately after tactile discrimination tasks, indicating an engagement of the somatosensory system during stimulus encoding. On a group level both the oscillatory power and the evoked potential amplitudes were not modulated by tACS neither after tactile finger stimulation nor after median nerve stimulation as compared to sham stimulation. On an individual level we could however demonstrate that lower angular difference (i.e., differences between the injected current vector in the target region and the source orientation vector) is associated with significantly higher changes in both P20/N20 and N30/P30 source activities. Our findings suggest that the higher the directionality of the injected current correlates to the dipole orientation the greater the tACS-induced aftereffects are.