Determining which mechanisms lead to activation in the motor cortex: a modeling study of transcranial magnetic stimulation using realistic stimulus waveforms and sulcal geometry

Axon morphology and intrinsic cellular properties determine repetitive transcranial magnetic stimulation threshold for plasticity

Introduction

Repetitive transcranial magnetic stimulation (rTMS) is a widely used therapeutic tool in neurology and psychiatry, but its cellular and molecular mechanisms are not fully understood. Standardizing stimulus parameters, specifically electric field strength, is crucial in experimental and clinical settings. It enables meaningful comparisons across studies and facilitates the translation of findings into clinical practice. However, the impact of biophysical properties inherent to the stimulated neurons and networks on the outcome of rTMS protocols remains not well understood. Consequently, achieving standardization of biological effects across different brain regions and subjects poses a significant challenge.

Methods

This study compared the effects of 10 Hz repetitive magnetic stimulation (rMS) in entorhino-hippocampal tissue cultures from mice and rats, providing insights into the impact of the same stimulation protocol on similar neuronal networks under standardized conditions.

Results

We observed the previously described plastic changes in excitatory and inhibitory synaptic strength of CA1 pyramidal neurons in both mouse and rat tissue cultures, but a higher stimulation intensity was required for the induction of rMS-induced synaptic plasticity in rat tissue cultures. Through systematic comparison of neuronal structural and functional properties and computational modeling, we found that morphological parameters of CA1 pyramidal neurons alone are insufficient to explain the observed differences between the groups. Although morphologies of mouse and rat CA1 neurons showed no significant differences, simulations confirmed that axon morphologies significantly influence individual cell activation thresholds. Notably, differences in intrinsic cellular properties were sufficient to account for the 10% higher intensity required for the induction of synaptic plasticity in the rat tissue cultures.

Conclusion

These findings demonstrate the critical importance of axon morphology and intrinsic cellular properties in predicting the plasticity effects of rTMS, carrying valuable implications for the development of computer models aimed at predicting and standardizing the biological effects of rTMS.

Probing intrahemispheric interactions with a novel dual-site TMS setup

OBJECTIVE

Using dual-site transcranial magnetic stimulation (dsTMS), the effective connectivity between the primary motor cortex (M1) and adjacent brain areas such as the dorsal premotor cortex (PMd) can be investigated. However, stimulating two brain regions in close proximity (e.g., ±2.3 cm for intrahemispheric PMd-M1) is subject to considerable spatial restrictions that potentially can be overcome by combining two standard figure-of-eight coils in a novel dsTMS setup.

METHODS

After a technical evaluation of its induced electric fields, the dsTMS setup was tested in vivo (n = 23) by applying a short-interval intracortical inhibition (SICI) protocol. Additionally, the intrahemispheric PMd-M1 interaction was probed. E-field modelling was performed using SimNIBS.

RESULTS

The technical evaluation yielded no major alterations of the induced electric fields due to coil overlap. In vivo, the setup reliably elicited SICI. Investigating intrahemispheric PMd-M1 interactions was feasible (inter-stimulus interval 6 ms), resulting in modulation of M1 output.

CONCLUSIONS

The presented dsTMS setup provides a novel way to stimulate two adjacent brain regions with fewer technical and spatial limitations than previous attempts.

SIGNIFICANCE

This dsTMS setup enables more accurate and repeatable targeting of brain regions in close proximity and can facilitate innovation in the field of effective connectivity.

Novel approaches to motoneuron disease/ALS treatment using non-invasive brain and spinal stimulation: IFCN handbook chapter

Non-invasive brain stimulation techniques have been exploited in motor neuron disease (MND) with multifold objectives: to support the diagnosis, to get insights in the pathophysiology of these disorders and, more recently, to slow down disease progression. In this review, we consider how neuromodulation can now be employed to treat MND, with specific attention to amyotrophic lateral sclerosis (ALS), the most common form with upper motoneuron (UMN) involvement, taking into account electrophysiological abnormalities revealed by human and animal studies that can be targeted by neuromodulation techniques. This review article encompasses repetitive transcranial magnetic stimulation methods (including low-frequency, high-frequency, and pattern stimulation paradigms), transcranial direct current stimulation as well as experimental findings with the newer approach of trans-spinal direct current stimulation. We also survey and discuss the trials that have been performed, and future perspectives.

Brain stimulation-on-a-chip: a neuromodulation platform for brain slices

Electrical stimulation of ex vivo brain tissue slices has been a method used to understand mechanisms imparted by transcranial direct current stimulation (tDCS), but there are significant direct current electric field (dcEF) dosage and electrochemical by-product concerns in conventional experimental setups that may impact translational findings. Therefore, we developed an on-chip platform with fluidic, electrochemical, and magnetically-induced spatial control. Fluidically, the chamber geometrically confines precise dcEF delivery to the enclosed brain slice and allows for tissue recovery in order to monitor post-stimulation effects. Electrochemically, conducting hydrogel electrodes mitigate stimulation-induced faradaic reactions typical of commonly-used metal electrodes. Magnetically, we applied ferromagnetic substrates beneath the tissue and used an external permanent magnet to enable in situ rotational control in relation to the dcEF. By combining the microfluidic chamber with live-cell calcium imaging and electrophysiological recordings, we showcased the potential to study the acute and lasting effects of dcEFs with the potential of providing multi-session stimulation. This on-chip bioelectronic platform presents a modernized yet simple solution to electrically stimulate explanted tissue by offering more environmental control to users, which unlocks new opportunities to conduct thorough brain stimulation mechanistic investigations.

Axon morphology and intrinsic cellular properties determine repetitive transcranial magnetic stimulation threshold for plasticity

Repetitive transcranial magnetic stimulation (rTMS) is a widely used therapeutic tool in neurology and psychiatry, but its cellular and molecular mechanisms are not fully understood. Standardizing stimulus parameters, specifically electric field strength and direction, is crucial in experimental and clinical settings. It enables meaningful comparisons across studies and facilitating the translation of findings into clinical practice. However, the impact of biophysical properties inherent to the stimulated neurons and networks on the outcome of rTMS protocols remains not well understood. Consequently, achieving standardization of biological effects across different brain regions and subjects poses a significant challenge. This study compared the effects of 10 Hz repetitive magnetic stimulation (rMS) in entorhino-hippocampal tissue cultures from mice and rats, providing insights into the impact of the same stimulation protocol on similar neuronal networks under standardized conditions. We observed the previously described plastic changes in excitatory and inhibitory synaptic strength of CA1 pyramidal neurons in both mouse and rat tissue cultures, but a higher stimulation intensity was required for the induction of rMS-induced synaptic plasticity in rat tissue cultures. Through systematic comparison of neuronal structural and functional properties and computational modeling, we found that morphological parameters of CA1 pyramidal neurons alone are insufficient to explain the observed differences between the groups. However, axon morphologies of individual cells played a significant role in determining activation thresholds. Notably, differences in intrinsic cellular properties were sufficient to account for the 10 % higher intensity required for the induction of synaptic plasticity in the rat tissue cultures. These findings demonstrate the critical importance of axon morphology and intrinsic cellular properties in predicting the plasticity effects of rTMS, carrying valuable implications for the development of computer models aimed at predicting and standardizing the biological effects of rTMS.

Revisiting the Rotational Field TMS Method for Neurostimulation

Transcranial magnetic stimulation (TMS) is a non-invasive technique that has shown high efficacy in the treatment of major depressive disorder (MDD) and is increasingly utilized for various neuropsychiatric disorders. However, conventional TMS is limited to activating only a small fraction of neurons that have components parallel to the induced electric field. This likely contributes to the significant variability observed in clinical outcomes. A novel method termed rotational field TMS (rfTMS or TMS 360°) enables the activation of a greater number of neurons by reducing the sensitivity to orientation. Recruitment of a larger number of neurons offers the potential to enhance efficacy and reduce variability in the treatment of clinical indications for which neuronal recruitment and organization may play a significant role, such as MDD and stroke. The potential of the method remains to be validated in clinical trials. Here, we revisit and describe in detail the rfTMS method, its principles, mode of operation, effects on the brain, and potential benefits for clinical TMS.

Transcranial Magnetic Stimulation and its Imaging Features in Patients With Depression, Post-traumatic Stress Disorder, and Traumatic Brain Injury

Transcranial magnetic stimulation (TMS) is a type of noninvasive neurostimulation used increasingly often in clinical medicine. While most studies to date have focused on TMS's ability to treat major depressive disorder, it has shown promise in several other conditions including post-traumatic stress disorder (PTSD) and traumatic brain injury (TBI). As different treatment protocols are often used across studies, the ability to predict patient outcomes and evaluate immediate and long-term changes using imaging becomes increasingly important. Several imaging features, such as thickness, connectedness, and baseline activity of a variety of cortical and subcortical areas, have been found to be correlated with a greater response to TMS therapy. Intrastimulation imaging can reveal in real time how TMS applied to superficial areas activates or inhibits activity in deeper brain regions. Functional imaging performed weeks to months after treatment can offer an understanding of how long-term effects on brain activity relate to clinical improvement. Further work should be done to expand our knowledge of imaging features relevant to TMS therapy and how they vary across patients with different neurological and psychiatric conditions.

Neuron matters: neuromodulation with electromagnetic stimulation must consider neurons as dynamic identities

Neuromodulation with electromagnetic stimulation is widely used for the control of abnormal neural activity, and has been proven to be a valuable alternative to pharmacological tools for the treatment of many neurological diseases. Tremendous efforts have been focused on the design of the stimulation apparatus (i.e., electrodes and magnetic coils) that delivers the electric current to the neural tissue, and the optimization of the stimulation parameters. Less attention has been given to the complicated, dynamic properties of the neurons, and their context-dependent impact on the stimulation effects. This review focuses on the neuronal factors that influence the outcomes of electromagnetic stimulation in neuromodulation. Evidence from multiple levels (tissue, cellular, and single ion channel) are reviewed. Properties of the neural elements and their dynamic changes play a significant role in the outcome of electromagnetic stimulation. This angle of understanding yields a comprehensive perspective of neural activity during electrical neuromodulation, and provides insights in the design and development of novel stimulation technology.

Inter-individual variability in current direction for common tDCS montages

The direction of applied electric current relative to the cortical surface is a key determinant of transcranial direct current stimulation (tDCS) effects. Inter-individual differences in anatomy affect the consistency of current direction at a cortical target. However, the degree of this variability remains undetermined. Using current flow modelling (CFM), we quantified the inter-individual variability in tDCS current direction at a cortical target (left primary motor cortex, M1). Three montages targeting M1 using circular electrodes were compared: PA-tDCS directed current perpendicular to the central sulcus in a posterior-anterior direction relative to M1, ML-tDCS directed current parallel to the central sulcus in a medio-lateral direction, and conventional-tDCS applied electrodes over M1 and the contralateral forehead. In 50 healthy brain scans from the Human Connectome Project, we extracted current direction and intensity from the grey matter surface in the sulcal bank (M1BANK) and gyral crown (M1CROWN), and neighbouring primary somatosensory cortex (S1BANK and S1CROWN). Results confirmed substantial inter-individual variability in current direction (50%-150%) across all montages. Radial inward current produced by PA-tDCS was predominantly located in M1BANK, whereas for conventional-tDCS it was clustered in M1CROWN. The difference in radial inward current in functionally distinct subregions of M1 raises the testable hypothesis that PA-tDCS and conventional-tDCS modulate cortical excitability through different mechanisms. We show that electrode locations can be used to closely approximate current direction in M1 and precentral gyrus, providing a landmark-based method for tDCS application to address the hypothesis without the need for MRI. By contrast, ML-tDCS current was more tangentially orientated, which is associated with weaker somatic polarisation. Substantial inter-individual variability in current direction likely contributes to variable neuromodulation effects reported for these protocols, emphasising the need for individualised electrode montages, including the control of current direction.

Transcranial magnetic stimulation of the brain: What is stimulated? - A consensus and critical position paper

Transcranial (electro)magnetic stimulation (TMS) is currently the method of choice to non-invasively induce neural activity in the human brain. A single transcranial stimulus induces a time-varying electric field in the brain that may evoke action potentials in cortical neurons. The spatial relationship between the locally induced electric field and the stimulated neurons determines axonal depolarization. The induced electric field is influenced by the conductive properties of the tissue compartments and is strongest in the superficial parts of the targeted cortical gyri and underlying white matter. TMS likely targets axons of both excitatory and inhibitory neurons. The propensity of individual axons to fire an action potential in response to TMS depends on their geometry, myelination and spatial relation to the imposed electric field and the physiological state of the neuron. The latter is determined by its transsynaptic dendritic and somatic inputs, intrinsic membrane potential and firing rate. Modeling work suggests that the primary target of TMS is axonal terminals in the crown top and lip regions of cortical gyri. The induced electric field may additionally excite bends of myelinated axons in the juxtacortical white matter below the gyral crown. Neuronal excitation spreads ortho- and antidromically along the stimulated axons and causes secondary excitation of connected neuronal populations within local intracortical microcircuits in the target area. Axonal and transsynaptic spread of excitation also occurs along cortico-cortical and cortico-subcortical connections, impacting on neuronal activity in the targeted network. Both local and remote neural excitation depend critically on the functional state of the stimulated target area and network. TMS also causes substantial direct co-stimulation of the peripheral nervous system. Peripheral co-excitation propagates centrally in auditory and somatosensory networks, but also produces brain responses in other networks subserving multisensory integration, orienting or arousal. The complexity of the response to TMS warrants cautious interpretation of its physiological and behavioural consequences, and a deeper understanding of the mechanistic underpinnings of TMS will be critical for advancing it as a scientific and therapeutic tool.

Cellular mechanisms underlying state-dependent neural inhibition with magnetic stimulation

Novel stimulation protocols for neuromodulation with magnetic fields are explored in clinical and laboratory settings. Recent evidence suggests that the activation state of the nervous system plays a significant role in the outcome of magnetic stimulation, but the underlying cellular and molecular mechanisms of state-dependency have not been completely investigated. We recently reported that high frequency magnetic stimulation could inhibit neural activity when the neuron was in a low active state. In this paper, we investigate state-dependent neural modulation by applying a magnetic field to single neurons, using the novel micro-coil technology. High frequency magnetic stimulation suppressed single neuron activity in a state-dependent manner. It inhibited neurons in slow-firing states, but spared neurons from fast-firing states, when the same magnetic stimuli were applied. Using a multi-compartment NEURON model, we found that dynamics of voltage-dependent sodium and potassium channels were significantly altered by the magnetic stimulation in the slow-firing neurons, but not in the fast-firing neurons. Variability in neural activity should be monitored and explored to optimize the outcome of magnetic stimulation in basic laboratory research and clinical practice. If selective stimulation can be programmed to match the appropriate neural state, prosthetic implants and brain-machine interfaces can be designed based on these concepts to achieve optimal results.

Numerical analysis of transcranial magnetic stimulation application in patients with orofacial pain

In this paper, we monitored the accuracy of non-navigated application of repetitive Transcranial Magnetic Stimulation (rTMS) in 10 patients suffering from orofacial pain by using functional magnetic resonance (fMRI), computer modeling and numerical simulation. Through a unique process, each fMRI scan was used to define a Region of Interest (ROI) where the source of the orofacial pain was located, which was to be stimulated using rTMS. For each patient, MRI scans with a spatial resolution of 0.7 mm were converted into an anatomically accurate head model. The head model including the ROI was then co-registered with a model of the stimulation coil in an electromagnetic field numerical simulator. The accuracy of rTMS application was evaluated based on the calculations of electric field intensity distribution in the ROI. The research has yielded unique insight into ROIs (with average volume 904mm3) in patients with orofacial pain and has also extended further possibilities of human head MRI image semi-automatic segmentation. According to the calculations performed, the average ROI volume that was stimulated by an electric field with an intensity of over 80 V/m was only 4.4%, with the maximum ROI volume being 20.5%. Furthermore, a numerical study of the impact of coil rotation and translation was performed. It demonstrated a) the optimal placement of the stimulation coil can significantly increase the volume of the stimulated ROI up to 60% and b) patients with orofacial pain would need precise coil positioning with a navigation error lower than 10 mm. Due to an acceptable proccessing time of up to 6 hours, described numerical simulation opens up new options for precise rTMS treatment planning. This planning platform together withpatient-specific navigated rTMS, could lead to significant increase of treatment outcomes in patients suffering from orofacial pain.

Multi-scale modeling toolbox for single neuron and subcellular activity under Transcranial Magnetic Stimulation

BACKGROUND

Transcranial Magnetic Stimulation (TMS) is a widely used non-invasive brain stimulation method. However, its mechanism of action and the neural response to TMS are still poorly understood. Multi-scale modeling can complement experimental research to study the subcellular neural effects of TMS. At the macroscopic level, sophisticated numerical models exist to estimate the induced electric fields. However, multi-scale computational modeling approaches to predict TMS cellular and subcellular responses, crucial to understanding TMS plasticity inducing protocols, are not available so far.

OBJECTIVE

We develop an open-source multi-scale toolbox Neuron Modeling for TMS (NeMo-TMS) to address this problem.

METHODS

NeMo-TMS generates accurate neuron models from morphological reconstructions, couples them to the external electric fields induced by TMS, and simulates the cellular and subcellular responses of single-pulse and repetitive TMS.

RESULTS

We provide examples showing some of the capabilities of the toolbox.

CONCLUSION

NeMo-TMS toolbox allows researchers a previously not available level of detail and precision in realistically modeling the physical and physiological effects of TMS.

Degree of improving TMS focality through a geometrically stable solution of an inverse TMS problem

The Transcranial Magnetic Stimulation (TMS) inverse problem (TMS-IP) investigated in this study aims to focus the TMS induced electric field close to a specified target point defined on the gray matter interface in the M1_HAND area while otherwise minimizing it. The goal of the study is to numerically evaluate the degree of improvement of the TMS-IP solutions relative to the well-known sulcus-aligned mapping (a projection approach with the 90^∘ local sulcal angle). In total, 1536 individual TMS-IP solutions have been analyzed for multiple target points and multiple subjects using the boundary element fast multipole method (BEM-FMM) as the forward solver. Our results show that the optimal TMS inverse-problem solutions improve the focality - reduce the size of the field "hot spot" and its deviation from the target - by approximately 21-33% on average for all considered subjects, all observation points, two distinct coil types, two segmentation types, two intracortical observation surfaces under study, and three tested values of the field threshold. The inverse-problem solutions with the maximized focality simultaneously improve the TMS mapping resolution (differentiation between neighbor targets separated by approximately 10 mm) although this improvement is quite modest. Coil position/orientation and conductivity uncertainties have been included into consideration as the corresponding de-focalization factors. The present results will change when the levels of uncertainties change. Our results also indicate that the accuracy of the head segmentation critically influences the expected TMS-IP performance.

Multiscale Computational Model Reveals Nerve Response in a Mouse Model for Temporal Interference Brain Stimulation

There has been a growing interest in the non-invasive stimulation of specific brain tissues, while reducing unintended stimulation in surrounding regions, for the medical treatment of brain disorders. Traditional methods for non-invasive brain stimulation, such as transcranial direct current stimulation (tDCS) or transcranial magnetic stimulation (TMS), can stimulate brain regions, but they also simultaneously stimulate the brain and non-brain regions that lie between the target and the stimulation site of the source. Temporal interference (TI) stimulation has been suggested to selectively stimulate brain regions by superposing two alternating currents with slightly different frequencies injected through electrodes attached to the scalp. Previous studies have reported promising results for TI applied to the motor area in mice, but the mechanisms are yet to be clarified. As computational techniques can help reveal different aspects of TI, in this study, we computationally investigated TI stimulation using a multiscale model that computes the generated interference current pattern effects in a neural cortical model of a mouse head. The results indicated that the threshold increased with the carrier frequency and that the beat frequency did not influence the threshold. It was also found that the intensity ratio between the alternating currents changed the location of the responding nerve, which is in agreement with previous experiments. Moreover, particular characteristics of the envelope were investigated to predict the stimulation region intuitively. It was found that regions with high modulation depth (| maximum| − | minimum| values of the envelope) and low minimum envelope (near zero) corresponded with the activation region obtained via neural computation.

The Myelin Content of the Human Precentral Hand Knob Reflects Interindividual Differences in Manual Motor Control at the Physiological and Behavioral Level

The primary motor cortex hand area (M1_HAND) and adjacent dorsal premotor cortex (PMd) form the so-called motor hand knob in the precentral gyrus. M1_HAND and PMd are critical for dexterous hand use and are densely interconnected via corticocortical axons, lacking a sharp demarcating border. In 24 young right-handed volunteers, we performed multimodal mapping to delineate the relationship between structure and function in the right motor hand knob. Quantitative structural magnetic resonance imaging (MRI) at 3 tesla yielded regional R1 maps as a proxy of cortical myelin content. Participants also underwent functional MRI (fMRI). We mapped task-related activation and temporal precision, while they performed a visuomotor synchronization task requiring visually cued abduction movements with the left index or little finger. We also performed sulcus-aligned transcranial magnetic stimulation of the motor hand knob to localize the optimal site (hotspot) for evoking a motor evoked potential (MEP) in two intrinsic hand muscles. Individual motor hotspot locations varied along the rostrocaudal axis. The more rostral the motor hotspot location in the precentral crown, the longer were corticomotor MEP latencies. “Hotspot rostrality” was associated with the regional myelin content in the precentral hand knob. Cortical myelin content also correlated positively with task-related activation of the precentral crown and temporal precision during the visuomotor synchronization task. Together, our results suggest a link among cortical myelination, the spatial cortical representation, and temporal precision of finger movements. We hypothesize that the myelination of cortical axons facilitates neuronal integration in PMd and M1_HAND and, hereby, promotes the precise timing of movements.

SIGNIFICANCE STATEMENT Here we used magnetic resonance imaging and transcranial magnetic stimulation of the precentral motor hand knob to test for a link among cortical myelin content, functional corticomotor representations, and manual motor control. A higher myelin content of the precentral motor hand knob was associated with more rostral corticomotor presentations, with stronger task-related activation and a higher precision of movement timing during a visuomotor synchronization task. We propose that a high precentral myelin content enables fast and precise neuronal integration in M1 (primary motor cortex) and dorsal premotor cortex, resulting in higher temporal precision during dexterous hand use. Our results identify the degree of myelination as an important structural feature of the neocortex that is tightly linked to the function and behavior supported by the cortical area.

The impact of individual electrical fields and anatomical factors on the neurophysiological outcomes of tDCS: A TMS-MEP and MRI study

BACKGROUND

Transcranial direct current stimulation (tDCS), a neuromodulatory non-invasive brain stimulation technique, has shown promising results in basic and clinical studies. The known interindividual variability of the effects, however, limits the efficacy of the technique. Recently we reported neurophysiological effects of tDCS applied over the primary motor cortex at the group level, based on data from twenty-nine participants who received 15min of either sham, 0.5, 1.0, 1.5 or 2.0 mA anodal, or cathodal tDCS. The neurophysiological effects were evaluated via changes in: 1) transcranial magnetic stimulation (TMS)-induced motor evoked potentials (MEP), and 2) cerebral blood flow (CBF) measured by functional magnetic resonance imaging (MRI) via arterial spin labeling (ASL). At the group level, dose-dependent effects of the intervention were obtained, which however displayed interindividual variability.

METHOD

In the present study, we investigated the cause of the observed inter-individual variability. To this end, for each participant, a MRI-based realistic head model was designed to 1) calculate anatomical factors and 2) simulate the tDCS- and TMS-induced electrical fields (EFs). We first investigated at the regional level which individual anatomical factors explained the simulated EFs (magnitude and normal component). Then, we explored which specific anatomical and/or EF factors predicted the neurophysiological outcomes of tDCS.

RESULTS

The results highlight a significant negative correlation between regional electrode-to-cortex distance (rECD) as well as regional CSF (rCSF) thickness, and the individual EF characteristics. In addition, while both rCSF thickness and rECD anticorrelated with tDCS-induced physiological changes, EFs positively correlated with the effects.

CONCLUSION

These results provide novel insights into the dependency of the neuromodulatory effects of tDCS on individual physical factors.

Direct comparison of efficacy of the motor cortical plasticity induction and the interindividual variability between TBS and QPS

BACKGROUND

Theta burst stimulation (TBS) and quadripulse stimulation (QPS) are known to induce synaptic plasticity in humans. There have been no head-to-head comparisons of the efficacy and variability between TBS and QPS.

OBJECTIVE

To compare the efficacy and interindividual variability between the original TBS and QPS protocols. We hypothesized that QPS would be more effective and less variable than TBS.

METHODS

Forty-six healthy subjects participated in this study. Thirty subjects participated in the main comparison experiment, and the other sixteen subjects participated in the experiment to obtain natural variation in motor-evoked potentials. The facilitatory effects were compared between intermittent TBS (iTBS) and QPS5, and the inhibitory effects were compared between continuous TBS (cTBS) and QPS50. The motor-evoked potential amplitudes elicited by transcranial magnetic stimulation over the primary motor cortex were measured before the intervention and every 5 min after the intervention for 1 h. To investigate the interindividual variability, the responder/nonresponder/opposite-responder rates were also analyzed.

RESULTS

The facilitatory effects of QPS5 were greater than those of iTBS, and the inhibitory effects of QPS50 were much stronger than those of cTBS. The responder rate of QPS was significantly higher than that of TBS. QPS had a smaller number of opposite responders than TBS.

CONCLUSION

QPS is more effective and stable for synaptic plasticity induction than TBS.

Optimal stimulation site for rTMS to improve motor function: Anatomical hand knob vs. hand motor hotspot

Repetitive transcranial magnetic stimulation (rTMS) is used to modulate neuronal excitability of the human brain. Distant effects on contralateral corticomotor excitability can be exerted by interhemispheric modulation by low-frequency rTMS on ipsilateral hemisphere. To modulate corticospinal excitability, accurate determination of the stimulation site is important to maximize the effects of rTMS. In the present study, we investigated the difference in the distant effect of 1 Hz rTMS with respect to inducing functional improvement in the non-dominant hand by inhibiting the dominant hemisphere depending on cortical target areas. Ten healthy right-handed volunteers without any neurological disorders were enrolled. The anatomical hand knob (HK) identified from individual magnetic resonance imaging and the transcranial magnetic stimulation (TMS) induced hand motor hotspot (hMHS) by recording motor evoked potentials (MEPs) in the contralateral first dorsal interosseous muscle were determined. All participants underwent three conditions of 1 Hz rTMS on left hemisphere intervention; rTMS application over the HK, rTMS application over the hMHS, and sham-rTMS. Before and after each intervention, all participants undergone motor function assessments with their left hand. The cortical mapping showed that the hMHS was located anteriorly and laterally compared to the HK. Motor function tests showed the most significant improvements after the hMHS stimulation. When we compared the distant effects of target site on corticospinal excitability and motor behavior, delivering 1 Hz rTMS to the hMHS was more effective than delivering it to the HK for improving corticomotor excitability, motor skill, and dexterity. These results suggest that TMS-induced hMHS is an optimal target area to induce distant effect of low-frequency rTMS in motor function.

Variance in cortical depth across the brain surface: Implications for transcranial stimulation of the brain

The distance between the surface of the scalp and the surface of the grey matter of the brain is a key factor in determining the effective dose of non-invasive brain stimulation for an individual person. The highly folded nature of the cortical surface means that the depth of a particular brain area is likely to vary between individuals. The question addressed here is: what is the variability of this measure of cortical depth? Ninety-four anatomical MRI images were taken from the OASIS database. For each image, the minimum distance from each point in the grey matter to the scalp surface was determined. Transforming these estimates into standard space meant that the coefficient of variation could be determined across the sample. The results indicated that depth variability is high across the cortical surface, even when taking sulcal depth into account. This was true even for the primary visual and motor areas, which are often used in setting TMS dosage. The correlation of the depth of these areas and the depth of other brain areas was low. The results suggest that dose setting of TMS based on visual or evoked potentials may offer poor reliability, and that individual brain images should be used when targeting non-primary brain areas.

A software toolkit for TMS electric-field modeling with boundary element fast multipole method: an efficient MATLAB implementation

Objective

To present and disseminate our transcranial magnetic stimulation (TMS) modeling software toolkit, including several new algorithmic developments, and to apply this software to realistic TMS modeling scenarios given a high-resolution model of the human head including cortical geometry and an accurate coil model.

Approach

The recently developed charge-based boundary element fast multipole method (BEM-FMM) is employed as an alternative to the 1st order finite element method (FEM) most commonly used today. The BEM-FMM approach provides high accuracy and unconstrained numerical field resolution close to and across cortical interfaces. Here, the previously proposed BEM-FMM algorithm has been improved in several novel ways.

Main results

The improvements resulted in a threefold increase in computational speed while maintaining the same solution accuracy. The computational code based on the MATLAB® platform is made available to all interested researchers, along with a coil model repository and examples to create custom coils, head model repository, and supporting documentation. The presented software toolkit may be useful for post-hoc analyses of navigated TMS data using high-resolution subject-specific head models as well as accurate and fast modeling for the purposes of TMS coil/hardware development.

Significance

TMS is currently the only non-invasive neurostimulation modality that enables painless and safe supra-threshold stimulation by employing electromagnetic induction to efficiently penetrate the skull. Accurate, fast, and high resolution modeling of the electric fields may significantly improve individualized targeting and dosing of TMS and therefore enhance the efficiency of existing clinical protocols as well as help establish new application domains.

Obtaining accurate and calibrated coil models for transcranial magnetic stimulation using magnetic field measurements

Currently, simulations of the induced currents in the brain produced by transcranial magnetic stimulation (TMS) are used to elucidate the regions reached by stimuli. However, models commonly found in the literature are too general and neglect imperfections in the windings. Aiming to predict the stimulation sites in patients requires precise modeling of the electric field (E-field), and a proper calibration to adequate to the empirical data of the particular coil employed. Furthermore, most fabricators do not provide precise information about the coil geometries, and even using X-ray images may lead to subjective interpretations. We measured the three components of the vector magnetic field induced by a TMS figure-8 coil with spatial resolutions of up to 1 mm. Starting from a computerized tomography-based coil model, we applied a multivariate optimization algorithm to automatically modify the original model and obtain one that optimally fits the measurements. Differences between models were assessed in a human brain mesh using the finite-elements method showing up to 6% variations in the E-field magnitude. Our calibrated model could increase the precision of the estimated E-field induced in the brain during TMS, enhance the accuracy of delivered stimulation during functional brain mapping, and improve dosimetry for repetitive TMS. Graphical Abstract Geometrical model of TMS coil based on TAC images is optimally deformed to match magnetic field measurements. The calibrated model's induced electric field in the brain differs from the original.

Rotational field TMS: Comparison with conventional TMS based on motor evoked potentials and thresholds in the hand and leg motor cortices

BACKGROUND

Transcranial magnetic stimulation (TMS) is a rapidly expanding technology utilized in research and neuropsychiatric treatments. Yet, conventional TMS configurations affect primarily neurons that are aligned parallel to the induced electric field by a fixed coil, making the activation orientation-specific. A novel method termed rotational field TMS (rfTMS), where two orthogonal coils are operated with a 90° phase shift, produces rotation of the electric field vector over almost a complete cycle, and may stimulate larger portion of the neuronal population within a given brain area.

OBJECTIVE

To compare the physiological effects of rfTMS and conventional unidirectional TMS (udTMS) in the motor cortex.

METHODS

Hand and leg resting motor thresholds (rMT), and motor evoked potential (MEP) amplitudes and latencies (at 120% of rMT), were measured using a dual-coil array based on the H7-coil, in 8 healthy volunteers following stimulation at different orientations of either udTMS or rfTMS.

RESULTS

For both target areas rfTMS produced significantly lower rMTs and much higher MEPs than those induced by udTMS, for comparable induced electric field amplitude. Both hand and leg rMTs were orientation-dependent.

CONCLUSIONS

rfTMS induces stronger physiologic effects in targeted brain regions at significantly lower intensities. Importantly, given the activation of a much larger population of neurons within a certain brain area, repeated application of rfTMS may induce different neuroplastic effects in neural networks, opening novel research and clinical opportunities.

Patient Semi-specific Computational Modeling of Electromagnetic Stimulation Applied to Neuroprotective Treatments in Acute Ischemic Stroke

Neuroprotective effects of pulsed electromagnetic fields (PEMFs) have been demonstrated both in vivo and in vitro. Moreover, preliminary clinical studies have been conducted and suggested PEMFs as a possible alternative therapy to treat acute ischemic stroke. In this work, we show that it's possible to build-up a patient semi-specific head model, where the 3D reconstruction of the ischemic lesion of the patient under treatment is inserted in the head of the human body model "Duke" (v.1.0, Zurich MedTech AG). The semi-specific model will be used in the randomized, placebo-controlled, double-blind study currently ongoing. Three patients were modelled and simulated, and results showed that each ischemic lesion experiences a magnetic flux density field comparable to the one for which biological effects have been attested. Such a kind of dosimetric analysis reveals a reliable tool to assess the correlation between levels of exposure and the beneficial effect. Thus, once the on-going double blind study is complete it will prove if PEMFs treatment triggers a clinical effect, and we will then be able to characterize a dose-response curve with the methodology arranged in this study.

The effect of coil placement and orientation on the assessment of focal excitability in motor mapping with navigated transcranial magnetic stimulation

BACKGROUND

Navigated transcranial magnetic stimulation (nTMS) is used for mapping muscle representations in the primary motor cortex. We used sulcus-aligned mapping and electric field (E-field) modeling to investigate the excitability of the motor hand area for further understanding the methodological limitations of nTMS.

NEW METHOD

We studied 10 healthy volunteers to locate the cortical target eliciting the largest responses (the hotspot) in the first dorsal interosseous (FDI) muscle. Six additional targets were placed along the central sulcus at 5-mm distances. Resting motor thresholds (rMTs) and optimal coil orientations were determined at all targets, and a conventional motor mapping was conducted. The cortical E-fields, induced by stimulating the targets with rMT intensities and optimal coil orientations, were modeled in a realistic head geometry to estimate the activated cortical sites.

RESULTS

The rMTs increased with increasing distance from the hotspot (p < 0.001). The greatest motor-evoked potential (MEP) amplitudes occurred with the coil perpendicular to the sulcus and with the coil pointing towards the hotspot or the center of gravity of the motor map. The E-field strengths at the hotspot (99±26 V/m) remained above previously estimated thresholds for activation.

COMPARISON WITH EXISTING METHODS

Depending on the target location, optimal coil orientations may deviate significantly from the conventional perpendicular-to-sulcus angle, which is often assumed optimal. These orientations seem to maintain the E-field stable in the hand knob, regardless of the sulcal shape near the stimulated target.

CONCLUSIONS

The coil orientation is crucial for the accuracy of motor mapping, and the apparent motor map may extend due to remote hotspot activation.

Cervical trans-spinal direct current stimulation: a modelling-experimental approach

Background

Trans-spinal direct current stimulation (tsDCS) is a non-invasive technique with promising neuromodulatory effects on spinal cord (SC) circuitry. Computational studies are essential to guide effective tsDCS protocols for specific clinical applications. This study aims to combine modelling and experimental studies to determine the electrode montage that maximizes electric field (E-field) delivery during cervical tsDCS.

Methods

Current and E-field distributions in the cervical SC were predicted for four electrode montages in a human realistic model using computational methods. A double-blind crossover and randomized exploratory study was conducted using the montage that maximized E-field delivery. tsDCS was applied for 15 min in 10 healthy subjects (anodal, cathodal, sham, with polarity assigned to the cervical electrode), with a current intensity of 2.5 mA, resulting in a total current charge density delivery of 90 mC/cm². Upper limb motor (transcranial magnetic stimulation) and sensory evoked potentials (MEP, SEP), M-waves, H-reflex and F-wave responses were analysed. Central and peripheral conduction times were determined using MEP. Repeated measures ANOVA and Friedman test were used for statistical analysis (significance level α = 0.05).

Results

All montages presented higher current density and E-field magnitudes in the cervical SC region between the electrodes. However, electrodes at C3 and T3 spinous processes (C3-T3) originated the highest E-field magnitude (0.50 V/m). Using C3-T3 montage we observed significant changes in N9 SEP latency (p = 0.006), but significance did not persist in pairwise comparisons (sham-anodal: p = 0.022; sham-cathodal: p = 0.619; anodal-cathodal: p = 0.018; α = 0.017, Bonferroni corrected). MEP latency and central motor conduction time (CMCT) modified significantly on stimulation (p = 0.007 and p = 0.015, respectively). In addition, pairwise comparisons confirmed significant differences between sham and cathodal conditions after Bonferroni correction for MEP latency (sham-anodal: p = 0.868; sham-cathodal: p = 0.011; anodal-cathodal: p = 0.023) and CMCT (sham-anodal: p = 0.929; sham-cathodal: p = 0.010; anodal-cathodal: p = 0.034).

Conclusions

Computational models predicted higher E-field delivery in the cervical SC for the C3-T3 montage. Polarity-dependent effects in motor responses were reported using this montage consistent with spinal motor modulation. tsDCS experimental protocol designs should be guided by modelling studies to improve effectiveness.

Neuromodulation of Axon Terminals

Understanding which cellular compartments are influenced during neuromodulation underpins any rational effort to explain and optimize outcomes. Axon terminals have long been speculated to be sensitive to polarization, but experimentally informed models for CNS stimulation are lacking. We conducted simultaneous intracellular recording from the neuron soma and axon terminal (blebs) during extracellular stimulation with weak sustained (DC) uniform electric fields in mouse cortical slices. Use of weak direct current stimulation (DCS) allowed isolation and quantification of changes in axon terminal biophysics, relevant to both suprathreshold (e.g., deep brain stimulation, spinal cord stimulation, and transcranial magnetic stimulation) and subthreshold (e.g., transcranial DCS and transcranial alternating current stimulation) neuromodulation approaches. Axon terminals polarized with sensitivity (mV of membrane polarization per V/m electric field) 4 times than somas. Even weak polarization (<2 mV) of axon terminals significantly changes action potential dynamics (including amplitude, duration, conduction velocity) in response to an intracellular pulse. Regarding a cellular theory of neuromodulation, we explain how suprathreshold CNS stimulation activates the action potential at terminals while subthreshold approaches modulate synaptic efficacy through axon terminal polarization. We demonstrate that by virtue of axon polarization and resulting changes in action potential dynamics, neuromodulation can influence analog-digital information processing.

Simulation of transcranial magnetic stimulation in head model with morphologically-realistic cortical neurons

BACKGROUND

Transcranial magnetic stimulation (TMS) enables non-invasive modulation of brain activity with both clinical and research applications, but fundamental questions remain about the neural types and elements TMS activates and how stimulation parameters affect the neural response.

OBJECTIVE

To develop a multi-scale computational model to quantify the effect of TMS parameters on the direct response of individual neurons.

METHODS

We integrated morphologically-realistic neuronal models with TMS-induced electric fields computed in a finite element model of a human head to quantify the cortical response to TMS with several combinations of pulse waveforms and current directions.

RESULTS

TMS activated with lowest intensity intracortical axonal terminations in the superficial gyral crown and lip regions. Layer 5 pyramidal cells had the lowest thresholds, but layer 2/3 pyramidal cells and inhibitory basket cells were also activated at most intensities. Direct activation of layers 1 and 6 was unlikely. Neural activation was largely driven by the field magnitude, rather than the field component normal to the cortical surface. Varying the induced current direction caused a waveform-dependent shift in the activation site and provided a potential mechanism for experimentally observed differences in thresholds and latencies of muscle responses.

CONCLUSIONS

This biophysically-based simulation provides a novel method to elucidate mechanisms and inform parameter selection of TMS and other cortical stimulation modalities. It also serves as a foundation for more detailed network models of the response to TMS, which may include endogenous activity, synaptic connectivity, inputs from intrinsic and extrinsic axonal projections, and corticofugal axons in white matter.

The Quasi-uniform assumption for Spinal Cord Stimulation translational research

BACKGROUND

Quasi-uniform assumption is a general theory that postulates local electric field predicts neuronal activation. Computational current flow model of spinal cord stimulation (SCS) of humans and animal models inform how the quasi-uniform assumption can support scaling neuromodulation dose between humans and translational animal.

NEW METHOD

Here we developed finite element models of cat and rat SCS, and brain slice, alongside SCS models. Boundary conditions related to species specific electrode dimensions applied, and electric fields per unit current (mA) predicted.

RESULTS

Clinically and across animal, electric fields change abruptly over small distance compared to the neuronal morphology, such that each neuron is exposed to multiple electric fields. Per unit current, electric fields generally decrease with body mass, but not necessarily and proportionally across tissues. Peak electric field in dorsal column rat and cat were ∼17x and ∼1x of clinical values, for scaled electrodes and equal current. Within the spinal cord, the electric field for rat, cat, and human decreased to 50% of peak value caudo-rostrally (C5-C6) at 0.48 mm, 3.2 mm, and 8 mm, and mediolaterally at 0.14 mm, 2.3 mm, and 3.1 mm. Because these space constants are different, electric field across species cannot be matched without selecting a region of interest (ROI).

COMPARISON WITH EXISTING METHOD

This is the first computational model to support scaling neuromodulation dose between humans and translational animal.

CONCLUSIONS

Inter-species reproduction of the electric field profile across the entire surface of neuron populations is intractable. Approximating quasi-uniform electric field in a ROI is a rational step to translational scaling.

Modified cable equation incorporating transverse polarization of neuronal membranes for accurate coupling of electric fields

OBJECTIVE

We present a theory and computational methods to incorporate transverse polarization of neuronal membranes into the cable equation to account for the secondary electric field generated by the membrane in response to transverse electric fields. The effect of transverse polarization on nonlinear neuronal activation thresholds is quantified and discussed in the context of previous studies using linear membrane models.

APPROACH

The response of neuronal membranes to applied electric fields is derived under two time scales and a unified solution of transverse polarization is given for spherical and cylindrical cell geometries. The solution is incorporated into the cable equation re-derived using an asymptotic model that separates the longitudinal and transverse dimensions. Two numerical methods are proposed to implement the modified cable equation. Several common neural stimulation scenarios are tested using two nonlinear membrane models to compare thresholds of the conventional and modified cable equations.

MAIN RESULTS

The implementations of the modified cable equation incorporating transverse polarization are validated against previous results in the literature. The test cases show that transverse polarization has limited effect on activation thresholds. The transverse field only affects thresholds of unmyelinated axons for short pulses and in low-gradient field distributions, whereas myelinated axons are mostly unaffected.

SIGNIFICANCE

The modified cable equation captures the membrane's behavior on different time scales and models more accurately the coupling between electric fields and neurons. It addresses the limitations of the conventional cable equation and allows sound theoretical interpretations. The implementation provides simple methods that are compatible with current simulation approaches to study the effect of transverse polarization on nonlinear membranes. The minimal influence by transverse polarization on axonal activation thresholds for the nonlinear membrane models indicates that predictions of stronger effects in linear membrane models with a fixed activation threshold are inaccurate. Thus, the conventional cable equation works well for most neuroengineering applications, and the presented modeling approach is well suited to address the exceptions.

Coupling Magnetically Induced Electric Fields to Neurons: Longitudinal and Transverse Activation

We present a theory and computational models to couple the electric field induced by magnetic stimulation to neuronal membranes. Based on the characteristics of magnetically induced electric fields and the modified cable equation that we developed previously, quasipotentials are derived as a simple and accurate approximation for coupling of the electric fields to neurons. The conventional and modified cable equations are used to simulate magnetic stimulation of long peripheral nerves by circular and figure-8 coils. Activation thresholds are obtained over a range of lateral and vertical coil positions for two nonlinear membrane models representing unmyelinated and myelinated straight axons and also for undulating myelinated axons. For unmyelinated straight axons, the thresholds obtained with the modified cable equation are significantly lower due to transverse polarization, and the spatial distributions of thresholds as a function of coil position differ significantly from predictions by the activating function. However, the activation thresholds of unmyelinated axons obtained with either cable equation are very high and beyond the output capabilities of conventional magnetic stimulators. For myelinated axons, threshold values are similar for both cable equations and within the range of magnetic stimulators. Whereas the transverse field contributes negligibly to the activation thresholds of myelinated fibers, axonal undulation can significantly increase or decrease thresholds depending on coil position. The analysis provides a rigorous theoretical foundation and implementation methods for the use of the cable equation to model neuronal response to magnetically induced electric fields. Experimentally observed stimulation with the electric fields perpendicular to the nerve trunk cannot be explained by transverse polarization and is likely due to nerve fiber undulation and other geometrical inhomogeneities.

Rewiring the Addicted Brain: Circuits-Based Treatment for Addiction

The advent of the noninvasive brain stimulation (NIBS) technique has paved the way for neural circuit-based treatments for addiction. Recently, evidence from both preclinical and clinical studies has evaluated the use of transcranial magnetic stimulation (TMS) as a safe and cost-effective therapeutic tool for substance use disorders (SUDs). Indeed, repetitive TMS impacts on neural activity inducing short- and long-term effects involving neuroplasticity mechanisms locally within the target area of stimulation and the network level throughout the brain. Here, we provide an integrated view of evidence highlighting the mechanisms of TMS-induced effects on modulating the maladaptive brain circuitry of addiction. We then review the preclinical and clinical findings suggesting rTMS as an effective interventional tool for the treatment of SUDs.

Neuronal tuning: Selective targeting of neuronal populations via manipulation of pulse width and directionality

INTRODUCTION

Motor evoked potentials (MEP) in response to anteroposterior transcranial (AP) magnetic stimulation (TMS) are sensitive to the TMS pulse shape. We are now able to isolate distinct pulse properties, such as pulse width and directionality and evaluate them individually. Different pulse shapes induce different effects, likely by stimulating different populations of neurons. This implies that not all neurons respond in the same manner to stimulation, possibly, because individual segments of neurons differ in their membrane properties.

OBJECTIVES

To investigate the effect of different pulse widths and directionalities of TMS on MEP latencies, motor thresholds and plastic aftereffects of rTMS.

METHODS

Using a controllable pulse stimulator TMS (cTMS), we stimulated fifteen subjects with quasi-unidirectional TMS pulses of different pulse durations (40 μs, 80 μs and 120 μs) and determined thresholds and MEP AP latencies. We then compared the effects of 80 μs quasi-unidirectional pulses to those of 80 μs pulses with different pulse directionality characteristics (0.6 and 1.0 M ratios). We applied 900 pulses of the selected pulse shapes at 1 Hz.

RESULTS

The aftereffects of 1 Hz rTMS depended on pulse shape and duration. 40 and 80 μs wide unidirectional pulses induced inhibition, 120 μs wide pulses caused excitation. Bidirectional pulses induced inhibition during the stimulation but had facilitatory aftereffects. Narrower pulse shapes caused longer latencies and higher resting motor thresholds (RMT) as compared to wider pulse shapes.

CONCLUSIONS

We can tune the aftereffects of rTMS by manipulating pulse width and directionality; this may be due to the different membrane properties of the various neuronal segments such as dendrites.

SIGNIFICANCE

To date, rTMS frequency has been the main determinant of the plastic aftereffects. However, we showed that pulse width also plays a major role, probably by recruiting novel neuronal targets.

Fluctuation of primary motor cortex excitability during cataplexy in narcolepsy

Objective

Cataplexy is a complicated and dynamic process in narcolepsy type 1 (NT1) patients. This study aimed to clarify the distinct stages during a cataplectic attack and identify the changes of the primary motor cortex (PMC) excitability during these stages.

Methods

Thirty-five patients with NT1 and 29 healthy controls were recruited to this study. Cataplectic stages were distinguished from a cataplectic attack by video-polysomnogram monitoring. Transcranial magnetic stimulation motor-evoked potential (TMS-MEP) was performed to measure the excitability of PMC during quiet wakefulness, laughter without cataplexy, and each cataplectic stage.

Results

Based on the video and electromyogram observations, a typical cataplectic attack (CA) process is divided into four stages: triggering (CA1), resisting (CA2), atonic (CA3), and recovering stages (CA4). Compared with healthy controls, NT1 patients showed significantly decreased intracortical facilitation during quiet wakefulness. During the laughter stage, both patients and controls showed increased MEP amplitude compared with quiet wakefulness. The MEP amplitude significantly increased even higher in CA1 and 2, and then dramatically decreased in CA3 accompanied with prolonged MEP latency compared with the laughter stage and quiet wakefulness. The MEP amplitude and latency gradually recovered during CA4.

Interpretation

This study identifies four stages during cataplectic attack and reveals the existence of a resisting stage that might change the process of cataplexy. The fluctuation of MEP amplitude and MEP latency shows a potential participation of PMC and motor control pathway during cataplexy, and the increased MEP amplitude during CA1 and 2 strongly implies a compensatory mechanism in motor control that may resist or avoid cataplectic attack.

Biophysically realistic neuron models for simulation of cortical stimulation

OBJECTIVE

We implemented computational models of human and rat cortical neurons for simulating the neural response to cortical stimulation with electromagnetic fields.

APPROACH

We adapted model neurons from the library of Blue Brain models to reflect biophysical and geometric properties of both adult rat and human cortical neurons and coupled the model neurons to exogenous electric fields (E-fields). The models included 3D reconstructed axonal and dendritic arbors, experimentally-validated electrophysiological behaviors, and multiple, morphological variants within cell types. Using these models, we characterized the single-cell responses to intracortical microstimulation (ICMS) and uniform E-field with dc as well as pulsed currents.

MAIN RESULTS

The strength-duration and current-distance characteristics of the model neurons to ICMS agreed with published experimental results, as did the subthreshold polarization of cell bodies and axon terminals by uniform dc E-fields. For all forms of stimulation, the lowest threshold elements were terminals of the axon collaterals, and the dependence of threshold and polarization on spatial and temporal stimulation parameters was strongly affected by morphological features of the axonal arbor, including myelination, diameter, and branching.

SIGNIFICANCE

These results provide key insights into the mechanisms of cortical stimulation. The presented models can be used to study various cortical stimulation modalities while incorporating detailed spatial and temporal features of the applied E-field.

The Effects of Waveform and Current Direction on the Efficacy and Test-Retest Reliability of Transcranial Magnetic Stimulation

The pulse waveform and current direction of transcranial magnetic stimulation (TMS) influence its interactions with the neural substrate; however, their role in the efficacy and reliability of single- and paired-pulse TMS measures is not fully understood. We investigated how pulse waveform and current direction affect the efficacy and test-retest reliability of navigated, single- and paired-pulse TMS measures. 23 healthy adults (aged 18-35 years) completed two identical TMS sessions, assessing resting motor threshold (RMT), motor-evoked potentials (MEPs), cortical silent period (cSP), short- and long-interval intra-cortical inhibition (SICI and LICI), and intracortical facilitation (ICF) using either monophasic posterior-anterior (monoPA; n = 9), monophasic anterior-posterior (monoAP; n = 7), or biphasic (biAP-PA; n = 7) pulses. Averages of each TMS measure were compared across the three groups and intraclass correlation coefficients were calculated to assess test-retest reliability. RMT was the lowest and cSP was the longest with biAP-PA pulses, whereas MEP latency was the shortest with monoPA pulses. SICI and LICI had the largest effect with monoPA pulses, whereas only monoAP and biAP-PA pulses resulted in significant ICF. MEP amplitude was more reliable with either monoPA or monoAP than with biAP-PA pulses. LICI was the most reliable with monoAP pulses, whereas ICF was the most reliable with biAP-PA pulses. Waveform/current direction influenced RMT, MEP latency, cSP, SICI, LICI, and ICF, as well as the reliability of MEP amplitude, LICI, and ICF. These results show the importance of considering TMS pulse parameters for optimizing the efficacy and reliability of TMS neurophysiologic measures.

Micro-Coil Design Influences the Spatial Extent of Responses to Intracortical Magnetic Stimulation

OBJECTIVE

Electrical stimulation via cortically implanted electrodes has been proposed to treat a wide range of neurological disorders. Effectiveness has been limited, however, in part due to the inability of conventional electrodes to activate specific types of neurons while avoiding other types. Recent demonstrations that magnetic stimulation from a micro-coil can selectively activate pyramidal neurons (PNs) while avoiding passing axons suggest the possibility that such an approach can overcome some this limitation and here we use computer simulations to explore how the micro-coil design influences the selectivity with which neurons are activated.

METHODS

A computational model was developed to compare the selectivity of magnetic stimulation induced by rectangular-, V-, and W-shaped coil designs. The more promising designs (V- and W-shapes) were fabricated for use in electrophysiological experiments including in vitro patch-clamp recording and calcium imaging (GCaMP6f) of mouse brain slices.

RESULTS

Both V- and W-shaped coils reliably activated layer 5 (L5) PNs but V-coils were more effective while W-coils were more selective. Activation thresholds with double-loop coils were approximately one-half those of single-loop coils. Calcium imaging revealed that both V- and W-coils better confine activation than electrodes.

CONCLUSION

Individual design features can influence both the strength as well as the selectivity of micro-coils and can be accurately predicted by computer simulations.

SIGNIFICANCE

Our results show that how coil design influences the response of cortical neurons to stimulation and are an important step toward the development of next-generation cortical prostheses.

Adults who stutter lack the specialised pre-speech facilitation found in non-stutterers

OBJECTIVES

Persistent developmental stuttering is a speech fluency disorder defined by its symptoms, where the underlying neurophysiological causes remain uncertain. This study examined the underlying neurophysiological mechanisms of the speech planning process, using facilitation in the motor cortex during speech preparation as an analogue.

METHODS

transcranial magnetic stimulation (TMS) pulses induced motor evoked potentials (MEPs), which were recorded from the tongue. Eighteen adults who stutter (AWS) and 17 adults who do not stutter (ANS) completed three experiments, which involved reading a German prefix+verb utterance from a screen. Each experiment involved 120 trials with three distinct levels of speech production: immediate speech, delayed speech without pacing and delayed speech with predefined pacing. TMS was applied shortly before speech onset. Trial MEPs were normalised to average non-speech MEPs. MEP amplitude, MEP facilitation ratio (amplitude: pre-speech offset) and group difference were the outcomes of interest analysed by multiple regression, as well as speech reaction time analysed by correlation.

RESULTS

MEP values were 11·1%-23·4% lower in AWS than ANS (by standardised Beta), across all three experiments. MEP facilitation ratio slopes were also 4·9%-18·3% flatter in AWS than ANS across all three experiments. Reaction times for AWS were only significantly slower than for ANS in immediate speech and predefined pacing experiments. No stuttering was detected during the trials. The group difference in immediate speech was 100% and 101% greater than the other two experiments respectively.

DISCUSSION

While performance of both ANS and AWS worsens under disturbed speech conditions, greater disturbance conditions affected controls worse than AWS. Future research and therapy in stuttering should focus on non-disturbed speech.

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.

Neuromodulation of lower limb motor responses with transcutaneous lumbar spinal cord direct current stimulation

OBJECTIVE

Trans-spinal direct current stimulation (tsDCS) is a promising technique to modulate spinal circuits. Combining clinical with modelling studies can improve effectiveness of tsDCS protocols. The aim of this study is to measure the effects of lumbar tsDCS on motor spinal responses and observe if these are consistent with the electric field (E-field) predicted from a computational model.

METHODS

The exploratory study design was double-blind crossover and randomized. tsDCS was delivered for 15 min (anodal, cathodal, sham) at L2 vertebra level (2.5 mA, 90 C/cm²) in 14 healthy subjects. F-wave, H-reflex, cortical silent period, motor evoked potential and sympathetic skin response were analyzed. Statistical methods were applied with Bonferroni correction for multiple comparisons, a p < 0.05 was set as significant. A human volume conductor model was obtained from available databases. E-field distributions in the spinal grey matter (GM) and white matter (WM) were calculated.

RESULTS

No tsDCS effects were observed. E-field magnitude predicted in the lumbosacral spinal GM and WM was <0.15 V/m, insufficient to ensure neuromodulation, which is consistent with the absence of effects.

CONCLUSION

The tsDCS protocol applied did not change motor response to delivered stimulus, thus we observed no effect on motor spinal circuits.

SIGNIFICANCE

Future tsDCS protocols should be supported by computational models.