Electric Field Model of Transcranial Electric Stimulation in Nonhuman Primates: Correspondence to Individual Motor Threshold

Bioelectronics for electrical stimulation: materials, devices and biomedical applications

Bioelectronics is a hot research topic, yet an important tool, as it facilitates the creation of advanced medical devices that interact with biological systems to effectively diagnose, monitor and treat a broad spectrum of health conditions. Electrical stimulation (ES) is a pivotal technique in bioelectronics, offering a precise, non-pharmacological means to modulate and control biological processes across molecular, cellular, tissue, and organ levels. This method holds the potential to restore or enhance physiological functions compromised by diseases or injuries by integrating sophisticated electrical signals, device interfaces, and designs tailored to specific biological mechanisms. This review explains the mechanisms by which ES influences cellular behaviors, introduces the essential stimulation principles, discusses the performance requirements for optimal ES systems, and highlights the representative applications. From this review, we can realize the potential of ES based bioelectronics in therapy, regenerative medicine and rehabilitation engineering technologies, ranging from tissue engineering to neurological technologies, and the modulation of cardiovascular and cognitive functions. This review underscores the versatility of ES in various biomedical contexts and emphasizes the need to adapt to complex biological and clinical landscapes it addresses.

A semi-automated pipeline for finite element modeling of electric field induced in nonhuman primates by transcranial magnetic stimulation

BACKGROUND

Transcranial magnetic stimulation (TMS) is used to treat a range of brain disorders by inducing an electric field (E-field) in the brain. However, the precise neural effects of TMS are not well understood. Nonhuman primates (NHPs) are used to model the impact of TMS on neural activity, but a systematic method of quantifying the induced E-field in the cortex of NHPs has not been developed.

NEW METHOD

The pipeline uses statistical parametric mapping (SPM) to automatically segment a structural MRI image of a rhesus macaque into five tissue compartments. Manual corrections are necessary around implants. The segmented tissues are tessellated into 3D meshes used in finite element method (FEM) software to compute the TMS induced E-field in the brain. The gray matter can be further segmented into cortical laminae using a volume preserving method for defining layers.

RESULTS

Models of three NHPs were generated with TMS coils placed over the precentral gyrus. Two coil configurations, active and sham, were simulated and compared. The results demonstrated a large difference in E-fields at the target. Additionally, the simulations were calculated using two different E-field solvers and were found to not significantly differ.

COMPARISON WITH EXISTING METHODS

Current methods segment NHP tissues manually or use automated methods for only the brain tissue. Existing methods also do not stratify the gray matter into layers.

CONCLUSION

The pipeline calculates the induced E-field in NHP models by TMS and can be used to plan implant surgeries and determine approximate E-field values around neuron recording sites.

Amplitude-determined seizure-threshold, electric field modeling, and electroconvulsive therapy antidepressant and cognitive outcomes

Electroconvulsive therapy (ECT) pulse amplitude, which dictates the induced electric field (E-field) magnitude in the brain, is presently fixed at 800 or 900 milliamperes (mA) without clinical or scientific rationale. We have previously demonstrated that increased E-field strength improves ECT's antidepressant effect but worsens cognitive outcomes. Amplitude-determined seizure titration may reduce the E-field variability relative to fixed amplitude ECT. In this investigation, we assessed the relationships among amplitude-determined seizure-threshold (STa), E-field magnitude, and clinical outcomes in older adults (age range 50 to 80 years) with depression. Subjects received brain imaging, depression assessment, and neuropsychological assessment pre-, mid-, and post-ECT. STa was determined during the first treatment with a Soterix Medical 4×1 High Definition ECT Multi-channel Stimulation Interface (Investigation Device Exemption: G200123). Subsequent treatments were completed with right unilateral electrode placement (RUL) and 800 mA. We calculated Ebrain defined as the 90th percentile of E-field magnitude in the whole brain for RUL electrode placement. Twenty-nine subjects were included in the final analyses. Ebrain per unit electrode current, Ebrain/I, was associated with STa. STa was associated with antidepressant outcomes at the mid-ECT assessment and bitemporal electrode placement switch. Ebrain/I was associated with changes in category fluency with a large effect size. The relationship between STa and Ebrain/I extends work from preclinical models and provides a validation step for ECT E-field modeling. ECT with individualized amplitude based on E-field modeling or STa has the potential to enhance neuroscience-based ECT parameter selection and improve clinical outcomes.

How electroconvulsive therapy works in the treatment of depression: is it the seizure, the electricity, or both?

We have known for nearly a century that triggering seizures can treat serious mental illness, but what we do not know is why. Electroconvulsive Therapy (ECT) works faster and better than conventional pharmacological interventions; however, those benefits come with a burden of side effects, most notably memory loss. Disentangling the mechanisms by which ECT exerts rapid therapeutic benefit from the mechanisms driving adverse effects could enable the development of the next generation of seizure therapies that lack the downside of ECT. The latest research suggests that this goal may be attainable because modifications of ECT technique have already yielded improvements in cognitive outcomes without sacrificing efficacy. These modifications involve changes in how the electricity is administered (both where in the brain, and how much), which in turn impacts the characteristics of the resulting seizure. What we do not completely understand is whether it is the changes in the applied electricity, or in the resulting seizure, or both, that are responsible for improved safety. Answering this question may be key to developing the next generation of seizure therapies that lack these adverse side effects, and ushering in novel interventions that are better, faster, and safer than ECT.

Cortical origin of theta error signals

A multi-scale approach elucidated the origin of the error-related-negativity (ERN), with its associated theta-rhythm, and the post-error-positivity (Pe) in macaque supplementary eye field (SEF). Using biophysical modeling, synaptic inputs to a subpopulation of layer-3 (L3) and layer-5 (L5) pyramidal cells (PCs) were optimized to reproduce error-related spiking modulation and inter-spike intervals. The intrinsic dynamics of dendrites in L5 but not L3 error PCs generate theta rhythmicity with random phases. Saccades synchronized the phases of the theta-rhythm, which was magnified on errors. Contributions from error PCs to the laminar current source density (CSD) observed in SEF were negligible and could not explain the observed association between error-related spiking modulation in L3 PCs and scalp-EEG. CSD from recorded laminar field potentials in SEF was comprised of multipolar components, with monopoles indicating strong electro-diffusion, dendritic/axonal electrotonic current leakage outside SEF, or violations of the model assumptions. Our results also demonstrate the involvement of secondary cortical regions, in addition to SEF, particularly for the later Pe component. The dipolar component from the observed CSD paralleled the ERN dynamics, while the quadrupolar component paralleled the Pe. These results provide the most advanced explanation to date of the cellular mechanisms generating the ERN.

Impact of modeled field of view in electroconvulsive therapy current flow simulations

Background

The field of view (FOV) considered in MRI-guided forward models of electroconvulsive therapy (ECT) are, as expected, limited to the MRI volume collected. Therefore, there is variation in model extent considered across simulation efforts. This study examines the impact of FOV on the induced electric field (E-field) due to two common electrode placements: right unilateral (RUL) and bilateral (BL).

Methods

A full-body dataset was obtained and processed for modeling relevant to ECT physics. Multiple extents were derived by truncating from the head down to four levels: upper head (whole-brain), full head, neck, and torso. All relevant stimulation and focality metrics were determined. The differences in the 99th percentile peak of stimulation strength in the brain between each extent to the full-body (reference) model were considered as the relative error (RE). We also determine the FOV beyond which the difference to a full-body model would be negligible.

Results

The 2D and 3D spatial plots revealed anticipated results in line with prior efforts. The RE for BL upper head was ~50% reducing to ~2% for the neck FOV. The RE for RUL upper head was ~5% reducing to subpercentage (0.28%) for the full-head FOV. As shown previously, BL was found to stimulate a larger brain volume-but restricted to the upper head and for amplitude up to ~480 mA. To some extent, RUL stimulated a larger volume. The RUL-induced volume was larger even when considering the neural activation threshold corresponding to brief pulse BL if ECT amplitude was >270 mA. This finding is explained by the BL-induced current loss through the inferior regions as more FOV is considered. Our result is a departure from prior efforts and raises questions about the focality metric as defined and/or inter-individual differences.

Conclusion

Our findings highlight that BL is impacted more than RUL with respect to FOV. It is imperative to collect full-head data at a minimum for any BL simulation and possibly more. Clinical practice resorts to using BL ECT when RUL is unsuccessful. However, the notion that BL is more efficacious on the premise of stimulating more brain volume needs to be revisited.

Responses of model cortical neurons to temporal interference stimulation and related transcranial alternating current stimulation modalities

Objective.Temporal interference stimulation (TIS) was proposed as a non-invasive, focal, and steerable deep brain stimulation method. However, the mechanisms underlying experimentally-observed suprathreshold TIS effects are unknown, and prior simulation studies had limitations in the representations of the TIS electric field (E-field) and cerebral neurons. We examined the E-field and neural response characteristics for TIS and related transcranial alternating current stimulation modalities.Approach.Using the uniform-field approximation, we simulated a range of stimulation parameters in biophysically realistic model cortical neurons, including different orientations, frequencies, amplitude ratios, amplitude modulation, and phase difference of the E-fields, and obtained thresholds for both activation and conduction block.Main results. For two E-fields with similar amplitudes (representative of E-field distributions at the target region), TIS generated an amplitude-modulated (AM) total E-field. Due to the phase difference of the individual E-fields, the total TIS E-field vector also exhibited rotation where the orientations of the two E-fields were not aligned (generally also at the target region). TIS activation thresholds (75-230 V m-1) were similar to those of high-frequency stimulation with or without modulation and/or rotation. For E-field dominated by the high-frequency carrier and with minimal amplitude modulation and/or rotation (typically outside the target region), TIS was less effective at activation and more effective at block. Unlike AM high-frequency stimulation, TIS generated conduction block with some orientations and amplitude ratios of individual E-fields at very high amplitudes of the total E-field (>1700 V m-1).Significance. The complex 3D properties of the TIS E-fields should be accounted for in computational and experimental studies. The mechanisms of suprathreshold cortical TIS appear to involve neural activity block and periodic activation or onset response, consistent with computational studies of peripheral axons. These phenomena occur at E-field strengths too high to be delivered tolerably through scalp electrodes and may inhibit endogenous activity in off-target regions, suggesting limited significance of suprathreshold TIS.

Detection of motor-evoked potentials below the noise floor: rethinking the motor stimulation threshold

Objective. Motor-evoked potentials (MEPs) are among the most prominent responses to brain stimulation, such as supra-threshold transcranial magnetic stimulation and electrical stimulation. Understanding of the neurophysiology and the determination of the lowest stimulation strength that evokes responses requires the detection of even smaller responses, e.g. from single motor units. However, available detection and quantization methods suffer from a large noise floor. This paper develops a detection method that extracts MEPs hidden below the noise floor. With this method, we aim to estimate excitatory activations of the corticospinal pathways well below the conventional detection level.Approach. The presented MEP detection method presents a self-learning matched-filter approach for improved robustness against noise. The filter is adaptively generated per subject through iterative learning. For responses that are reliably detected by conventional detection, the new approach is fully compatible with established peak-to-peak readings and provides the same results but extends the dynamic range below the conventional noise floor.Main results. In contrast to the conventional peak-to-peak measure, the proposed method increases the signal-to-noise ratio by more than a factor of 5. The first detectable responses appear to be substantially lower than the conventional threshold definition of 50µV median peak-to-peak amplitude.Significance. The proposed method shows that stimuli well below the conventional 50µV threshold definition can consistently and repeatably evoke muscular responses and thus activate excitable neuron populations in the brain. As a consequence, the input-output (IO) curve is extended at the lower end, and the noise cut-off is shifted. Importantly, the IO curve extends so far that the 50µV point turns out to be closer to the center of the logarithmic sigmoid curve rather than close to the first detectable responses. The underlying method is applicable to a wide range of evoked potentials and other biosignals, such as in electroencephalography.

Resolving the mesoscopic missing link: Biophysical modeling of EEG from cortical columns in primates

Event-related potentials (ERP) are among the most widely measured indices for studying human cognition. While their timing and magnitude provide valuable insights, their usefulness is limited by our understanding of their neural generators at the circuit level. Inverse source localization offers insights into such generators, but their solutions are not unique. To address this problem, scientists have assumed the source space generating such signals comprises a set of discrete equivalent current dipoles, representing the activity of small cortical regions. Based on this notion, theoretical studies have employed forward modeling of scalp potentials to understand how changes in circuit-level dynamics translate into macroscopic ERPs. However, experimental validation is lacking because it requires in vivo measurements of intracranial brain sources. Laminar local field potentials (LFP) offer a mechanism for estimating intracranial current sources. Yet, a theoretical link between LFPs and intracranial brain sources is missing. Here, we present a forward modeling approach for estimating mesoscopic intracranial brain sources from LFPs and predict their contribution to macroscopic ERPs. We evaluate the accuracy of this LFP-based representation of brain sources utilizing synthetic laminar neurophysiological measurements and then demonstrate the power of the approach in vivo to clarify the source of a representative cognitive ERP component. To that end, LFP was measured across the cortical layers of visual area V4 in macaque monkeys performing an attention demanding task. We show that area V4 generates dipoles through layer-specific transsynaptic currents that biophysically recapitulate the ERP component through the detailed forward modeling. The constraints imposed on EEG production by this method also revealed an important dissociation between computational and biophysical contributors. As such, this approach represents an important bridge between laminar microcircuitry, through the mesoscopic activity of cortical columns to the patterns of EEG we measure at the scalp.

TMS intensity and focality correlation with coil orientation at three non-motor regions

Objective. The aim of this study is to define the best coil orientations for transcranial magnetic stimulation (TMS) for three clinically relevant brain areas: pre-supplementary motor area (pre-SMA), inferior frontal gyrus (IFG), and posterior parietal cortex (PPC), by means of simulations in 12 realistic head models of the electric field (E-field). Methods. We computed the E-field generated by TMS in our three volumes of interest (VOI) that were delineated based on published atlases. We then analysed the maximum intensity and spatial focality for the normal and absolute components of the E-field considering different percentile thresholds. Lastly, we correlated these results with the different anatomical properties of our VOIs. Results. Overall, the spatial focality of the E-field for the three VOIs varied depending on the orientation of the coil. Further analysis showed that differences in individual brain anatomy were related to the amount of focality achieved. In general, a larger percentage of sulcus resulted in better spatial focality. Additionally, a higher normal E-field intensity was achieved when the coil axis was placed perpendicular to the predominant orientations of the gyri of each VOI. A positive correlation between spatial focality and E-field intensity was found for PPC and IFG but not for pre-SMA. Conclusions. For a rough approximation, better coil orientations can be based on the individual’s specific brain morphology at the VOI. Moreover, TMS computational models should be employed to obtain better coil orientations in non-motor regions of interest. Significance. Finding better coil orientations in non-motor regions is a challenge in TMS and seeks to reduce interindividual variability. Our individualized TMS simulation pipeline leads to fewer inter-individual variability in the focality, likely enhancing the efficacy of the stimulation and reducing the risk of stimulating adjacent, non-targeted areas.

Adeno-Associated Virus Capsid-Promoter Interactions in the Brain Translate from Rat to the Nonhuman Primate

Recently, we established an adeno-associated virus (AAV9) capsid-promoter interaction that directly determined cell-specific gene expression across two synthetic promoters, Cbh and CBA, in the rat striatum. These studies not only expand this capsid-promoter interaction to include another promoter in the rat striatum but also establish AAV capsid-promoter interactions in the nonhuman primate brain. When AAV serotype 9 (AAV9) vectors were injected into the rat striatum, the minimal synthetic promoter JetI drove green fluorescent protein (GFP) gene expression predominantly in oligodendrocytes. However, similar to our previous findings, the insertion of six alanines into VP1/VP2 of the AAV9 capsid (AAV9AU) significantly shifted JetI-driven GFP gene expression to neurons. In addition, previous retrograde tracing studies in the nonhuman primate brain also revealed the existence of a capsid-promoter interaction. When rAAV2-Retro vectors were infused into the frontal eye field (FEF) of rhesus macaques, local gene expression was prominent using either the hybrid chicken beta actin (CAG) or human synapsin (hSyn) promoters. However, only the CAG promoter, not the hSyn promoter, led to gene expression in the ipsilateral claustrum and contralateral FEF. Conversely, infusion of rAAV2-retro-hSyn vectors, but not rAAV2-retro-CAG, into the macaque superior colliculus led to differential and selective retrograde gene expression in cerebellotectal afferent cells. Clearly, this differential promoter/capsid expression profile could not be attributed to promoter inactivation from retrograde transport of the rAAV2-Retro vector. In summary, we document the potential for AAV capsid/promoter interactions to impact cell-specific gene expression across species, experimental manipulations, and engineered capsids, independent of capsid permissivity.

A Minimal Biophysical Model of Neocortical Pyramidal Cells: Implications for Frontal Cortex Microcircuitry and Field Potential Generation

Ca²⁺ spikes initiated in the distal trunk of layer 5 pyramidal cells (PCs) underlie nonlinear dynamic changes in the gain of cellular response, critical for top-down control of cortical processing. Detailed models with many compartments and dozens of ionic channels can account for this Ca²⁺ spike-dependent gain and associated critical frequency. However, current models do not account for all known Ca²⁺-dependent features. Previous attempts to include more features have required increasing complexity, limiting their interpretability and utility for studying large population dynamics. We overcome these limitations in a minimal two-compartment biophysical model. In our model, a basal-dendrites/somatic compartment included fast-inactivating Na⁺ and delayed-rectifier K⁺ conductances, while an apical-dendrites/trunk compartment included persistent Na⁺, hyperpolarization-activated cation (I _h ), slow-inactivating K⁺, muscarinic K⁺, and Ca²⁺ L-type. The model replicated the Ca²⁺ spike morphology and its critical frequency plus three other defining features of layer 5 PC synaptic integration: linear frequency-current relationships, back-propagation-activated Ca²⁺ spike firing, and a shift in the critical frequency by blocking I _h Simulating 1000 synchronized layer 5 PCs, we reproduced the current source density patterns evoked by Ca²⁺ spikes and describe resulting medial-frontal EEG on a male macaque monkey. We reproduced changes in the current source density when I _h was blocked. Thus, a two-compartment model with five crucial ionic currents in the apical dendrites reproduces all features of these neurons. We discuss the utility of this minimal model to study the microcircuitry of agranular areas of the frontal lobe involved in cognitive control and responsible for event-related potentials, such as the error-related negativity.SIGNIFICANCE STATEMENT A minimal model of layer 5 pyramidal cells replicates all known features crucial for distal synaptic integration in these neurons. By redistributing voltage-gated and returning transmembrane currents in the model, we establish a theoretical framework for the investigation of cortical microcircuit contribution to intracranial local field potentials and EEG. This tractable model will enable biophysical evaluation of multiscale electrophysiological signatures and computational investigation of cortical processing.

Using rAAV2-retro in rhesus macaques: Promise and caveats for circuit manipulation

BACKGROUND

Recent genetic technologies such as opto- and chemogenetics allow for the manipulation of brain circuits with unprecedented precision. Most studies employing these techniques have been undertaken in rodents, but a more human-homologous model for studying the brain is the nonhuman primate (NHP). Optimizing viral delivery of transgenes encoding actuator proteins could revolutionize the way we study neuronal circuits in NHPs. NEW METHOD: rAAV2-retro, a popular new capsid variant, produces robust retrograde labeling in rodents. Whether rAAV2-retro's highly efficient retrograde transport would translate to NHPs was unknown. Here, we characterized the anatomical distribution of labeling following injections of rAAV2-retro encoding opsins or DREADDs in the cortico-basal ganglia and oculomotor circuits of rhesus macaques.

RESULTS

rAAV2-retro injections in striatum, frontal eye field, and superior colliculus produced local labeling at injection sites and robust retrograde labeling in many afferent regions. In every case, however, a few brain regions with well-established projections to the injected structure lacked retrogradely labeled cells. We also observed robust terminal field labeling in downstream structures.

COMPARISON WITH EXISTING METHOD(S)

Patterns of labeling were similar to those obtained with traditional tract-tracers, except for some afferent labeling that was noticeably absent.

CONCLUSIONS

rAAV2-retro promises to be useful for circuit manipulation via retrograde transduction in NHPs, but caveats were revealed by our findings. Some afferently connected regions lacked retrogradely labeled cells, showed robust axon terminal labeling, or both. This highlights the importance of anatomically characterizing rAAV2-retro's expression in target circuits in NHPs before moving to manipulation studies.

Mixed Finite Element Method for Full-Wave Simulation of Bioelectromagnetism From DC to Microwave Frequencies

Bioelectromagnetism focuses on the study of electromagnetic fields in biological tissues from direct current (DC) to optical frequencies. It is challenging to develop an electromagnetics (EM) simulation method to cover this entire frequency band due to the electrically small/large scattering problem at extremely low/high frequencies. This paper focuses on the band from DC to microwave frequencies in bioelectromagnetism. Its main research objective is to develop a method that can overcome the low frequency breakdown problem at low frequencies (practically DC) and still stay stable at microwave frequencies. Based on the scattered field vector Helmholtz equation, the mixed finite element method (mixed FEM) is developed for the broadband electromagnetic field simulation in biological tissues. By imposing Gauss' law as the constraint condition, the mixed FEM overcomes the low frequency breakdown problem without resorting to the quasi-static approximation and remains effective and accurate at high frequencies. Extremely low frequency and high frequency numerical results are demonstrated to verify that the mixed FEM is a stable full-wave electromagnetic field simulation method for the full-bandwidth bioelectromagnetism.

Noninvasive Detection of Motor-Evoked Potentials in Response to Brain Stimulation Below the Noise Floor-How Weak Can a Stimulus Be and Still Stimulate

Motor-evoked potentials (MEP) are one of the most important responses to brain stimulation, such as supra-threshold transcranial magnetic stimulation (TMS) and electrical stimulation. The understanding of the neurophysiology and the determination of the lowest stimulation strength that evokes responses requires the detection of even smallest responses, e.g., from single motor units, but available detection and quantization methods are rather simple and suffer from a large noise floor. The paper introduces a more sophisticated matched-filter detection method that increases the detection sensitivity and shows that activation occurs well below the conventional detection level. In consequence, also conventional threshold definitions, e.g., as 50 μV median response amplitude, turn out to be substantially higher than the point at which first detectable responses occur. The presented method uses a matched-filter approach for improved sensitivity and generates the filter through iterative learning from the presented data. In contrast to conventional peak-to-peak measures, the presented method has a higher signal-to-noise ratio (≥14 dB). For responses that are reliably detected by conventional detection, the new approach is fully compatible and provides the same results but extends the dynamic range below the conventional noise floor. The underlying method is applicable to a wide range of well-timed biosignals and evoked potentials, such as in electroencephalography.

Statistical Model of Motor-Evoked Potentials

Motor-evoked potentials (MEPs) are widely used for biomarkers and dose individualization in transcranial stimulation. The large variability of MEPs requires sophisticated methods of analysis to extract information fast and correctly. Development and testing of such methods relies on the availability for realistic models of MEP generation, which are presently lacking. This paper presents a statistical model that can simulate long sequences of individualized MEP amplitude data with properties matching experimental observations. The MEP model includes three sources of trial-to-trial variability: excitability fluctuations, variability in the neural and muscular pathways, and physiological and measurement noise. It also generates virtual human subject data from statistics of population variability. All parameters are extracted as statistical distributions from experimental data from the literature. The model exhibits previously described features, such as stimulus-intensity-dependent MEP amplitude distributions, including bimodal ones. The model can generate long sequences of test data for individual subjects with specified parameters or for subjects from a virtual population. The presented MEP model is the most detailed to date and can be used for the development and implementation of dosing and biomarker estimation algorithms for transcranial stimulation.

Individual differences in TMS sensitivity influence the efficacy of tDCS in facilitating sensorimotor adaptation

BACKGROUND

Transcranial direct current stimulation (tDCS) can enhance cognitive function in healthy individuals, with promising applications as a therapeutic intervention. Despite this potential, variability in the efficacy of tDCS has been a considerable concern.

OBJECTIVE

/Hypothesis: Given that tDCS is always applied at a set intensity, we examined whether individual differences in sensitivity to brain stimulation might be one variable that modulates the efficacy of tDCS in a motor learning task.

METHODS

In the first part of the experiment, single-pulse transcranial magnetic stimulation (TMS) over primary motor cortex (M1) was used to determine each participant's resting motor threshold (rMT). This measure was used as a proxy of individual sensitivity to brain stimulation. In an experimental group of 28 participants, 2 mA tDCS was then applied during a motor learning task with the anodal electrode positioned over left M1. Another 14 participants received sham stimulation.

RESULTS

M1-Anodal tDCS facilitated learning relative to participants who received sham stimulation. Of primary interest was a within-group analysis of the experimental group, showing that the rate of learning was positively correlated with rMT: Participants who were more sensitive to brain stimulation as operationalized by our TMS proxy (low rMT), showed faster adaptation.

CONCLUSIONS

Methodologically, the results indicate that TMS sensitivity can predict tDCS efficacy in a behavioral task, providing insight into one source of variability that may contribute to replication problems with tDCS. Theoretically, the results provide further evidence of a role of sensorimotor cortex in adaptation, with the boost from tDCS observed during acquisition.

Immediate neurophysiological effects of transcranial electrical stimulation

Noninvasive brain stimulation techniques are used in experimental and clinical fields for their potential effects on brain network dynamics and behavior. Transcranial electrical stimulation (TES), including transcranial direct current stimulation (tDCS) and transcranial alternating current stimulation (tACS), has gained popularity because of its convenience and potential as a chronic therapy. However, a mechanistic understanding of TES has lagged behind its widespread adoption. Here, we review data and modelling on the immediate neurophysiological effects of TES in vitro as well as in vivo in both humans and other animals. While it remains unclear how typical TES protocols affect neural activity, we propose that validated models of current flow should inform study design and artifacts should be carefully excluded during signal recording and analysis. Potential indirect effects of TES (e.g., peripheral stimulation) should be investigated in more detail and further explored in experimental designs. We also consider how novel technologies may stimulate the next generation of TES experiments and devices, thus enhancing validity, specificity, and reproducibility.

A Computational Assessment of Target Engagement in the Treatment of Auditory Hallucinations with Transcranial Direct Current Stimulation

We use auditory verbal hallucinations (AVH) to illustrate the challenges in defining and assessing target engagement in the context of transcranial direct current stimulation (tDCS) for psychiatric disorders. We defined the target network as the cluster of regions of interest (ROIs) that are consistently implicated in AVH based on the conjunction of multimodal meta-analytic neuroimaging data. These were prescribed in the New York Head (a population derived model) and head models of four single individuals. We appraised two potential measures of target engagement, tDCS-induced peak electric field strength and tDCS-modulated volume defined as the percentage of the volume of the AVH network exposed to electric field magnitude stronger than the postulated threshold for neuronal excitability. We examined a left unilateral (LUL) montage targeting the prefrontal cortex (PFC) and temporoparietal junction (TPJ), a bilateral (BL) prefrontal montage, and a 2 × 1 montage targeting the left PFC and the TPJ bilaterally. Using computational modeling, we estimated the peak electric field strength and modulated volume induced by each montage for current amplitudes ranging 1-4 mA. We found that the LUL montage was inferior to both other montages in terms of peak electric field strength in right-sided AVH-ROIs. The BL montage was inferior to both other montages in terms of modulated volume of the left-sided AVH-ROIs. As the modulated volume is non-linear, its variability between montages reduced for current amplitudes above 3 mA. These findings illustrate how computational target engagement for tDCS can be tailored to specific networks and provide a principled approach for future study design.

Minimum Electric Field Exposure for Seizure Induction with Electroconvulsive Therapy and Magnetic Seizure Therapy

Lowering and individualizing the current amplitude in electroconvulsive therapy (ECT) has been proposed as a means to produce stimulation closer to the neural activation threshold and more focal seizure induction, which could potentially reduce cognitive side effects. However, the effect of current amplitude on the electric field (E-field) in the brain has not been previously linked to the current amplitude threshold for seizure induction. We coupled MRI-based E-field models with amplitude titrations of motor threshold (MT) and seizure threshold (ST) in four nonhuman primates (NHPs) to determine the strength, distribution, and focality of stimulation in the brain for four ECT electrode configurations (bilateral, bifrontal, right-unilateral, and frontomedial) and magnetic seizure therapy (MST) with cap coil on vertex. At the amplitude-titrated ST, the stimulated brain subvolume (23-63%) was significantly less than for conventional ECT with high, fixed current (94-99%). The focality of amplitude-titrated right-unilateral ECT (25%) was comparable to cap coil MST (23%), demonstrating that ECT with a low current amplitude and focal electrode placement can induce seizures with E-field as focal as MST, although these electrode and coil configurations affect differently specific brain regions. Individualizing the current amplitude reduced interindividual variation in the stimulation focality by 40-53% for ECT and 26% for MST, supporting amplitude individualization as a means of dosing especially for ECT. There was an overall significant correlation between the measured amplitude-titrated ST and the prediction of the E-field models, supporting a potential role of these models in dosing of ECT and MST. These findings may guide the development of seizure therapy dosing paradigms with improved risk/benefit ratio.

The influence of corticospinal activity on TMS-evoked activity and connectivity in healthy subjects: A TMS-EEG study

Combined transcranial magnetic stimulation (TMS) and electroencephalography (EEG) can be used to analyze cortical reactivity and connectivity. However, the effects of corticospinal and peripheral muscle activity on TMS-evoked potentials (TEPs) are not well understood. The aim of this paper is to evaluate the relationship between cortico-spinal activity, in the form of peripheral motor-evoked potentials (MEPs), and the TEPs from motor areas, along with the connectivity among activated brain areas. TMS was applied to left and right motor cortex (M1), separately, at motor threshold while multi-channel EEG responses were recorded in 17 healthy human subjects. Cortical excitability and source imaging analysis were performed for all trials at each stimulation location, as well as comparing trials resulting in MEPs to those without. Connectivity analysis was also performed comparing trials resulting in MEPs to those without. Cortical excitability results significantly differed between the MEP and no-MEP conditions for left M1 TMS at 60 ms (CP1, CP3, C1) and for right M1 TMS at 54 ms (CP6, C6). Connectivity analysis revealed higher outflow and inflow between M1 and somatosensory cortex bi-directionally for trials with MEPs than those without for both left M1 TMS (at 60, 100, 164 ms) and right M1 TMS (at 54, 100, and 164 ms). Both TEP amplitudes and connectivity measures related to motor and somatosensory areas ipsilateral to the stimulation were shown to correspond with peripheral MEP amplitudes. This suggests that cortico-spinal activation, along with the resulting somatosensory feedback, affects the cortical activity and dynamics within motor areas reflected in the TEPs. The findings suggest that TMS-EEG, along with adaptive connectivity estimators, can be used to evaluate the cortical dynamics associated with sensorimotor integration and proprioceptive manipulation along with the influence of peripheral muscle feedback.

Transcranial electric stimulation seen from within the brain

Computer models can make transcranial electric stimulation a better tool for research and therapy.

Novel Neuromodulation Techniques to Assess Interhemispheric Communication in Neural Injury and Neurodegenerative Diseases

Interhemispheric interaction has a major role in various neurobehavioral functions. Its disruption is a major contributor to the pathological changes in the setting of brain injury such as traumatic brain injury, peripheral nerve injury, and stroke, as well as neurodegenerative diseases. Because interhemispheric interaction has a crucial role in functional consequence in these neuropathological states, a review of noninvasive and state-of-the-art molecular based neuromodulation methods that focus on or have the potential to elucidate interhemispheric interaction have been performed. This yielded approximately 170 relevant articles on human subjects or animal models. There has been a recent surge of reports on noninvasive methods such as transcranial magnetic stimulation and transcranial direct current stimulation. Since these are noninvasive techniques with little to no side effects, their widespread use in clinical studies can be easily justified. The overview of novel neuromodulation methods and how they can be applied to study the role of interhemispheric communication in neural injury and neurodegenerative disease is provided. Additionally, the potential of each method in therapeutic use as well as investigating the pathophysiology of interhemispheric interaction in neurodegenerative diseases and brain injury is discussed. New technologies such as transcranial magnetic stimulation or transcranial direct current stimulation could have a great impact in understanding interhemispheric pathophysiology associated with acquired injury and neurodegenerative diseases, as well as designing improved rehabilitation therapies. Also, advances in molecular based neuromodulation techniques such as optogenetics and other chemical, thermal, and magnetic based methods provide new capabilities to stimulate or inhibit a specific brain location and a specific neuronal population.

Measurements and models of electric fields in the in vivo human brain during transcranial electric stimulation

Transcranial electric stimulation aims to stimulate the brain by applying weak electrical currents at the scalp. However, the magnitude and spatial distribution of electric fields in the human brain are unknown. We measured electric potentials intracranially in ten epilepsy patients and estimated electric fields across the entire brain by leveraging calibrated current-flow models. When stimulating at 2 mA, cortical electric fields reach 0.8 V/m, the lower limit of effectiveness in animal studies. When individual whole-head anatomy is considered, the predicted electric field magnitudes correlate with the recorded values in cortical (r = 0.86) and depth (r = 0.88) electrodes. Accurate models require adjustment of tissue conductivity values reported in the literature, but accuracy is not improved when incorporating white matter anisotropy or different skull compartments. This is the first study to validate and calibrate current-flow models with in vivo intracranial recordings in humans, providing a solid foundation to target stimulation and interpret clinical trials.

DOI:http://dx.doi.org/10.7554/eLife.18834.001

Sensorimotor Rhythm BCI with Simultaneous High Definition-Transcranial Direct Current Stimulation Alters Task Performance

BACKGROUND

Transcranial direct current stimulation (tDCS) has been used to alter the excitability of neurons within the cerebral cortex. Improvements in motor learning have been found in multiple studies when tDCS was applied to the motor cortex before or during task learning. The motor cortex is also active during the performance of motor imagination, a cognitive task during which a person imagines, but does not execute, a movement. Motor imagery can be used with noninvasive brain computer interfaces (BCIs) to control virtual objects in up to three dimensions, but to master control of such devices requires long training times.

OBJECTIVE

To evaluate the effect of high-definition tDCS on the performance and underlying electrophysiology of motor imagery based BCI.

METHODS

We utilize high-definition tDCS to investigate the effect of stimulation on motor imagery-based BCI performance across and within sessions over multiple training days.

RESULTS

We report a decreased time-to-hit with anodal stimulation both within and across sessions. We also found differing electrophysiological changes of the stimulated sensorimotor cortex during online BCI task performance for left vs. right trials. Cathodal stimulation led to a decrease in alpha and beta band power during task performance compared to sham stimulation for right hand imagination trials.

CONCLUSION

These results suggest that unilateral tDCS over the sensorimotor motor cortex differentially affects cortical areas based on task specific neural activation.

Animal models of transcranial direct current stimulation: Methods and mechanisms

The objective of this review is to summarize the contribution of animal research using direct current stimulation (DCS) to our understanding of the physiological effects of transcranial direct current stimulation (tDCS). We comprehensively address experimental methodology in animal studies, broadly classified as: (1) transcranial stimulation; (2) direct cortical stimulation in vivo and (3) in vitro models. In each case advantages and disadvantages for translational research are discussed including dose translation and the overarching "quasi-uniform" assumption, which underpins translational relevance in all animal models of tDCS. Terminology such as anode, cathode, inward current, outward current, current density, electric field, and uniform are defined. Though we put key animal experiments spanning decades in perspective, our goal is not simply an exhaustive cataloging of relevant animal studies, but rather to put them in context of ongoing efforts to improve tDCS. Cellular targets, including excitatory neuronal somas, dendrites, axons, interneurons, glial cells, and endothelial cells are considered. We emphasize neurons are always depolarized and hyperpolarized such that effects of DCS on neuronal excitability can only be evaluated within subcellular regions of the neuron. Findings from animal studies on the effects of DCS on plasticity (LTP/LTD) and network oscillations are reviewed extensively. Any endogenous phenomena dependent on membrane potential changes are, in theory, susceptible to modulation by DCS. The relevance of morphological changes (galvanotropy) to tDCS is also considered, as we suggest microscopic migration of axon terminals or dendritic spines may be relevant during tDCS. A majority of clinical studies using tDCS employ a simplistic dose strategy where excitability is singularly increased or decreased under the anode and cathode, respectively. We discuss how this strategy, itself based on classic animal studies, cannot account for the complexity of normal and pathological brain function, and how recent studies have already indicated more sophisticated approaches are necessary. One tDCS theory regarding "functional targeting" suggests the specificity of tDCS effects are possible by modulating ongoing function (plasticity). Use of animal models of disease are summarized including pain, movement disorders, stroke, and epilepsy.

Comparison of electric field strength and spatial distribution of electroconvulsive therapy and magnetic seizure therapy in a realistic human head model

BACKGROUND

This study examines the strength and spatial distribution of the electric field induced in the brain by electroconvulsive therapy (ECT) and magnetic seizure therapy (MST).

METHODS

The electric field induced by standard (bilateral, right unilateral, and bifrontal) and experimental (focal electrically administered seizure therapy and frontomedial) ECT electrode configurations as well as a circular MST coil configuration was simulated in an anatomically realistic finite element model of the human head. Maps of the electric field strength relative to an estimated neural activation threshold were used to evaluate the stimulation strength and focality in specific brain regions of interest for these ECT and MST paradigms and various stimulus current amplitudes.

RESULTS

The standard ECT configurations and current amplitude of 800-900mA produced the strongest overall stimulation with median of 1.8-2.9 times neural activation threshold and more than 94% of the brain volume stimulated at suprathreshold level. All standard ECT electrode placements exposed the hippocampi to suprathreshold electric field, although there were differences across modalities with bilateral and right unilateral producing respectively the strongest and weakest hippocampal stimulation. MST stimulation is up to 9 times weaker compared to conventional ECT, resulting in direct activation of only 21% of the brain. Reducing the stimulus current amplitude can make ECT as focal as MST.

CONCLUSIONS

The relative differences in electric field strength may be a contributing factor for the cognitive sparing observed with right unilateral compared to bilateral ECT, and MST compared to right unilateral ECT. These simulations could help understand the mechanisms of seizure therapies and develop interventions with superior risk/benefit ratio.

Individualized Low-Amplitude Seizure Therapy: Minimizing Current for Electroconvulsive Therapy and Magnetic Seizure Therapy

Electroconvulsive therapy (ECT) at conventional current amplitudes (800-900 mA) is highly effective but carries the risk of cognitive side effects. Lowering and individualizing the current amplitude may reduce side effects by virtue of a less intense and more focal electric field exposure in the brain, but this aspect of ECT dosing is largely unexplored. Magnetic seizure therapy (MST) induces a weaker and more focal electric field than ECT; however, the pulse amplitude is not individualized and the minimum amplitude required to induce a seizure is unknown. We titrated the amplitude of long stimulus trains (500 pulses) as a means of determining the minimum current amplitude required to induce a seizure with ECT (bilateral, right unilateral, bifrontal, and frontomedial electrode placements) and MST (round coil on vertex) in nonhuman primates. Furthermore, we investigated a novel method of predicting this amplitude-titrated seizure threshold (ST) by a non-convulsive measurement of motor threshold (MT) using single pulses delivered through the ECT electrodes or MST coil. Average STs were substantially lower than conventional pulse amplitudes (112-174 mA for ECT and 37.4% of maximum device amplitude for MST). ST was more variable in ECT than in MST. MT explained 63% of the ST variance and is hence the strongest known predictor of ST. These results indicate that seizures can be induced with less intense electric fields than conventional ECT that may be safer; efficacy and side effects should be evaluated in clinical studies. MT measurement could be a faster and safer alternative to empirical ST titration for ECT and MST.