Computational investigation of circuit mechanisms underlying short-latency responses to cortical stimulation
- 1Duke University
Abstract
Transcranial magnetic stimulation (TMS) over the primary motor cortex (M1) elicits a series of high frequency volleys termed D- and I-waves measured epidurally in the corticospinal tract of awake humans. Further, intracortical microstimulation (ICMS) in M1 of non-human primates evokes D- and I-wave responses similar to those observed in TMS. The cortical circuits and mechanisms involved in the generation of D- and I-waves by stimulation of M1 remain unclear. Here, we implemented computational models of cortical columns with laminarly-organized biophysically-based neurons, following existing models published in the literature: (1) M1 - single compartment (SC), (2) M1 - multi-compartment (MC), (3) primary auditory cortex (A1) - MC, and (4) primary somatosensory cortex (S1) - MC. The network connectivity of each model represented wiring found in the respective cortical regions. The direct effects of stimulation-induced electric fields were modeled as activation of different proportions of pyramidal neurons (PNs) across layers, and dose response curves were constructed for layer 5 (L5) PNs. Both the M1 and A1 models reproduced D- and I-waves, with the magnitude of I-waves increasing with higher recruitment of layer 2/3 and layer 5 PNs. The S1-MC model evoked rhythmic firing activity but with timings mismatched to experimental I-waves. The models replicated the experimentally observed effects of pharmacological agents on I-waves, and virtual lesions of specific neural populations across models revealed plausible microcircuit explanations for the first and later I-waves. This comprehensive comparison of models across multiple cortical regions identified consistent mechanisms underlying the cortical response to TMS and contributes to the refinement of computational strategies for optimizing stimulation paradigms.
Related
- Predictability and controllability shape aversive learning and stress responses through independent computational mechanisms
- Mathematical model for visual cortical development based on molecular mechanisms underlying the generation and elimination of dendritic spines
- Dietary Modulation of Neurobiological Mechanisms Underlying Executive Function and Decision-Making: A Systematic Review
- Genetic mechanisms underlying the associations between brain structure and mental health
- Ecological Mechanisms Underlying the Effects of Forest–Medicinal Plant Agroforestry Systems on Soil Quality and the Understory Environment
BibTeX
@misc{kumaravelu2026computational,
title = {{Computational investigation of circuit mechanisms underlying short-latency responses to cortical stimulation}},
author = {K. Kumaravelu and G. J. Yu and A. S. Aberra and M. A. Sommer and A. V. Peterchev and W. M. Grill},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.11.750854},
url = {https://www.biorxiv.org/content/10.64898/2026.09.11.750854v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - Computational investigation of circuit mechanisms underlying short-latency responses to cortical stimulation
AU - K. Kumaravelu
AU - G. J. Yu
AU - A. S. Aberra
AU - M. A. Sommer
AU - A. V. Peterchev
AU - W. M. Grill
PY - 2026
DA - 2026/09/11
PB - bioRxiv
DO - 10.64898/2026.09.11.750854
UR - https://www.biorxiv.org/content/10.64898/2026.09.11.750854v1
AB - Transcranial magnetic stimulation (TMS) over the primary motor cortex (M1) elicits a series of high frequency volleys termed D- and I-waves measured epidurally in the corticospinal tract of awake humans. Further, intracortical microstimulation (ICMS) in M1 of non-human primates evokes D- and I-wave responses similar to those observed in TMS. The cortical circuits and mechanisms involved in the generation of D- and I-waves by stimulation of M1 remain unclear. Here, we implemented computational models of cortical columns with laminarly-organized biophysically-based neurons, following existing models published in the literature: (1) M1 - single compartment (SC), (2) M1 - multi-compartment (MC), (3) primary auditory cortex (A1) - MC, and (4) primary somatosensory cortex (S1) - MC. The network connectivity of each model represented wiring found in the respective cortical regions. The direct effects of stimulation-induced electric fields were modeled as activation of different proportions of pyramidal neurons (PNs) across layers, and dose response curves were constructed for layer 5 (L5) PNs. Both the M1 and A1 models reproduced D- and I-waves, with the magnitude of I-waves increasing with higher recruitment of layer 2/3 and layer 5 PNs. The S1-MC model evoked rhythmic firing activity but with timings mismatched to experimental I-waves. The models replicated the experimentally observed effects of pharmacological agents on I-waves, and virtual lesions of specific neural populations across models revealed plausible microcircuit explanations for the first and later I-waves. This comprehensive comparison of models across multiple cortical regions identified consistent mechanisms underlying the cortical response to TMS and contributes to the refinement of computational strategies for optimizing stimulation paradigms.
ER -