Axon trajectories in local circuits of the primary motor cortex in the macaque monkey (Macaca fuscata)

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.

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.

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.

Where does TMS Stimulate the Motor Cortex? Combining Electrophysiological Measurements and Realistic Field Estimates to Reveal the Affected Cortex Position

Much of our knowledge on the physiological mechanisms of transcranial magnetic stimulation (TMS) stems from studies which targeted the human motor cortex. However, it is still unclear which part of the motor cortex is predominantly affected by TMS. Considering that the motor cortex consists of functionally and histologically distinct subareas, this also renders the hypotheses on the physiological TMS effects uncertain. We use the finite element method (FEM) and magnetic resonance image-based individual head models to get realistic estimates of the electric field induced by TMS. The field changes in different subparts of the motor cortex are compared with electrophysiological threshold changes of 2 hand muscles when systematically varying the coil orientation in measurements. We demonstrate that TMS stimulates the region around the gyral crown and that the maximal electric field strength in this region is significantly related to the electrophysiological response. Our study is one of the most extensive comparisons between FEM-based field calculations and physiological TMS effects so far, being based on data for 2 hand muscles in 9 subjects. The results help to improve our understanding of the basic mechanisms of TMS. They also pave the way for a systematic exploration of realistic field estimates for dosage control in TMS.

How the brain tissue shapes the electric field induced by transcranial magnetic stimulation

In transcranial magnetic stimulation (TMS), knowledge of the distribution of the induced electric field is fundamental for a better understanding of the position and extent of the stimulated brain region. However, the different tissue types and the varying fibre orientation in the brain tissue result in an inhomogeneous and anisotropic conductivity distribution and distort the electric field in a non-trivial way. Here, the field induced by a figure-8 coil is characterized in detail using finite element calculations and a geometrically accurate model of an individual head combined with high-resolution diffusion-weighted imaging for conductivity mapping. It is demonstrated that the field strength is significantly enhanced when the currents run approximately perpendicular to the local gyral orientation. Importantly, the spatial distribution of this effect differs distinctly between gray matter (GM) and white matter (WM): While the field in GM is selectively enhanced at the gyral crowns and lips, high field strengths can still occur rather deep in WM. Taking the anisotropy of brain tissue into account tends to further boost this effect in WM, but not in GM. Spatial variations in the WM anisotropy affect the local field strength in a systematic way and result in localized increases of up to 40% (on average ~7% for coil orientations perpendicular to the underlying gyri). We suggest that these effects might create hot spots in WM that might contribute to the excitation of WM structures by TMS. However, our results also demonstrate the necessity of using realistic nerve models in the future to allow for more definitive conclusions.

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

OBJECTIVE

To determine which mechanisms lead to activation of neurons in the motor cortex during transcranial magnetic stimulation (TMS) with different current directions and pulse waveforms.

METHODS

The total electric field induced in a simplified model of a cortical sulcus by a figure-eight coil was calculated using the finite element method (FEM). This electric field was then used as the input to determine the response of compartmental models of several types of neurons.

RESULTS

The modeled neurons were stimulated at different sites: fiber bends for pyramidal tract neurons, axonal terminations for cortical interneurons and axon collaterals, and a combination of both for pyramidal association fibers. All neurons were more easily stimulated by a PA - directed electric field, except association fibers. Additionally, the second phase of a biphasic pulse was found to be more efficient than the first phase of either monophasic or biphasic pulses.

CONCLUSIONS

The stimulation threshold for different types of neurons depends on the pulse waveform and relative current direction. The reported results might account for the range of responses obtained in TMS of the motor cortex when using different stimulation parameters.

SIGNIFICANCE

Modeling studies combining electric field calculations and neuronal models may lead to a deeper understanding of the effect of the TMS-induced electric field on cortical tissue, and may be used to improve TMS coil and waveform design.

Modeling the Effects of Transcranial Magnetic Stimulation on Cortical Circuits

Transcranial magnetic stimulation (TMS) is commonly used to activate or inactivate specific cortical areas in a noninvasive manner. Because of technical constraints, the precise effects of TMS on cortical circuits are difficult to assess experimentally. Here, this issue is investigated by constructing a detailed model of a portion of the thalamocortical system and examining the effects of the simulated delivery of a TMS pulse. The model, which incorporates a large number of physiological and anatomical constraints, includes 33,000 spiking neurons arranged in a 3-layered motor cortex and over 5 million intra- and interlayer synaptic connections. The model was validated by reproducing several results from the experimental literature. These include the frequency, timing, dose response, and pharmacological modulation of epidurally recorded responses to TMS (the so-called I-waves), as well as paired-pulse response curves consistent with data from several experimental studies. The modeled responses to simulated TMS pulses in different experimental paradigms provide a detailed, self-consistent account of the neural and synaptic activities evoked by TMS within prototypical cortical circuits.

A voltage-gated potassium channel, Kv3.1b, is expressed by a subpopulation of large pyramidal neurons in layer 5 of the macaque monkey cortex

In the cerebral cortex, the voltage-gated potassium channel, Kv3.1b, a splicing variant of Kv3.1, has been associated with fast-firing interneurons. Here, we report strong expression of Kv3.1b-protein and mRNA in both Betz and Meynert pyramidal cells of the monkey, as shown by immunohistochemistry and in situ hybridization. Strong expression also occurs in large pyramidal neurons in layer 5 of several cortical areas. In addition, most of these Betz and layer 5 pyramids, and about 10% of the labeled Meynert cells weakly co-expressed the calcium binding protein parvalbumin. Electron microscopy shows that the expression of Kv3.1b is localized to the somal and proximal dendritic cytoplasmic membrane, as expected for a channel protein. These results suggest that some large pyramidal neurons may constitute a functional subpopulation, with a distinctive distribution of voltage-gated potassium channels capable of influencing their repetitive firing properties.