Circuits and Mechanisms for TMS-Induced Corticospinal Waves: Connecting Sensitivity Analysis to the Network Graph

ABSTRACT Transcranial magnetic stimulation (TMS) is a non-invasive, FDA-cleared treatment for neuropsychiatric disorders with broad potential for new applications, but the neural circuits that are engaged during TMS are still poorly understood. Recordings of neural activity from the corticospinal tract provide a direct readout of the response of motor cortex to TMS, and therefore a new opportunity to model neural circuit dynamics. The study goal was to use epidural recordings from the cervical spine of human subjects to develop a computational model of a motor cortical macrocolumn through which the mechanisms underlying the response to TMS, including direct and indirect waves, could be investigated. An in-depth sensitivity analysis was conducted to identify important pathways, and machine learning was used to identify common circuit features among these pathways. Sensitivity analysis identified neuron types that preferentially contributed to single corticospinal waves. Single wave preference could be predicted using the average connection probability of all possible paths between the activated neuron type and L5 pyramidal tract neurons (PTNs). For these activations, the total conduction delay of the shortest path to L5 PTNs determined the latency of the corticospinal wave. Finally, there were multiple neuron type activations that could preferentially modulate a particular corticospinal wave. The results support the hypothesis that different pathways of circuit activation contribute to different corticospinal waves with participation of both excitatory and inhibitory neurons. Moreover, activation of both afferents to the motor cortex as well as specific neuron types within the motor cortex initiated different I-waves, and the results were interpreted to propose the cortical origins of afferents that may give rise to certain I-waves. The methodology provides a workflow for performing computationally tractable sensitivity analyses on complex models and relating the results to the network structure to both identify and understand mechanisms underlying the response to acute stimulation. AUTHOR SUMMARY Understanding circuit mechanisms underlying the response to transcranial magnetic stimulation remains a significant challenge for translational and clinical research. Computational models can reconstruct network activity in response to stimulation, but basic sensitivity analyses are insufficient to identify the fundamental circuit properties that underly an evoked response. We developed a data-driven neuronal network model of motor cortex, constrained with human recordings, that reproduced the corticospinal response to magnetic stimulation. The model supported several hypotheses, e.g., the importance of stimulating incoming fibers as well as neurons within the cortical column and the relevance of both excitatory and inhibitory neurons. Following a sensitivity analysis, we conducted a secondary structural analysis that linked the results of the sensitivity analysis to the network using machine learning. The structural analysis pointed to anatomical mechanisms that contributed to specific peaks in the response. Generally, given the anatomy and circuit of a neural region, identifying strongly connected paths in the network and the conduction delays of these paths can screen for important contributors to response peaks. This work supports and expands on hypotheses explaining the response to transcranial magnetic stimulation and adds a novel method for identifying generalizable neural circuit mechanisms.