Computational modeling of peripheral nerve fibers is a key tool for designing improved neuromodulation therapies. The NEURON software is commonly used to create biophysical simulations of nerve fibers, often in the outdated HOC language. Whether written in HOC or Python, implementing fiber simulations involves a steep learning curve and requires a large amount of standard, boilerplate code that is typically written anew for each project. There is a need for a code package that standardizes and simplifies the creation of model fibers, the execution of simulations of electrical stimulation, and the analysis of the resulting data. We created PyFibers, a NEURON-Python package that provides tools for accomplishing all these tasks and supports the development of new fiber models and stimulation protocols. PyFibers includes 11 fiber models from prior publications under a shared framework, and we validated each model's implementation in PyFibers against the original results. Our open-source tool simplifies and standardizes computational modeling of peripheral nerve fiber responses to electrical stimulation, providing a platform for the development of novel therapies using electrical stimulation, block, and recording.

BACKGROUND: Peripheral nerve recordings can enhance the efficacy of neurostimulation therapies by providing a feedback signal to adjust stimulation settings for greater efficacy or reduced side effects. Computational models can accelerate the development of interfaces with high signal-to-noise ratio and selective recording. However, validation and tuning of model outputs against in vivo recordings remains computationally prohibitive due to the large number of fibers in a nerve. METHODS: We designed and implemented highly efficient modeling methods for simulating electrically evoked compound nerve action potential (CNAP) signals. The method simulated a subset of fiber diameters present in the nerve using NEURON, interpolated action potential templates across fiber diameters, and filtered the templates with a weighting function derived from fiber-specific conduction velocity and electromagnetic reciprocity outputs of a volume conductor model. We applied the methods to simulate CNAPs from rat cervical vagus nerve. RESULTS: Brute force simulation of a rat vagal CNAP with all 1,759 myelinated and 13,283 unmyelinated fibers in NEURON required 286 and 15,860 CPU hours, respectively, while filtering interpolated templates required 30 and 38 seconds on a desktop computer while maintaining accuracy. Modeled CNAP amplitude could vary by over two orders of magnitude depending on tissue conductivities and cuff opening within experimentally relevant ranges. Conduction distance and fiber diameter distribution also strongly influenced the modeled CNAP amplitude, shape, and latency. Modeled and in vivo signals had comparable shape, amplitude, and latency for myelinated fibers but not for unmyelinated fibers. CONCLUSIONS: Highly efficient methods of modeling neural recordings quantified the large impact that tissue properties, conduction distance, and nerve fiber parameters have on CNAPs. These methods expand the computational accessibility of neural recording models, enable efficient model tuning for validation, and facilitate the design of novel recording interfaces for neurostimulation feedback and understanding physiological systems.

Unmyelinated C-fibers constitute the vast majority of axons in peripheral nerves and play key roles in homeostasis and signaling pain. However, little is known about their ion channel expression, which controls their firing properties. Also, because of their small diameters (~ 1 μm), it has not been possible to characterize their membrane properties using voltage clamp. We developed a novel library of isoform-specific ion channel models to serve as the basis functions of our C-fiber models. We then developed a particle swarm optimization (PSO) framework that used the isoform-specific ion channel models to reverse engineer C-fiber membrane properties from measured autonomic and cutaneous C-fiber conduction responses. Our C-fiber models reproduced experimental conduction velocity, chronaxie, action potential duration, intracellular threshold, and paired pulse recovery cycle. The models also matched experimental activity-dependent slowing, a property not included in model optimization. We found that simple conduction responses, characterizing the action potential, were controlled by similar membrane properties in both the autonomic and cutaneous C-fiber models, but complicated conduction response, characterizing the afterpotenials, were controlled by differential membrane properties. The unmyelinated C-fiber models constitute important tools to study autonomic signaling, assess the mechanisms of pain, and design bioelectronic devices. Additionally, the novel reverse engineering approach can be applied to generate models of other neurons where voltage clamp data are not available.

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.

Electrical stimulation and block of peripheral nerves hold great promise for treatment of a range of disease and disorders, but promising results from preclinical studies often fail to translate to successful clinical therapies. Differences in neural anatomy across species require different electrodes and stimulation parameters to achieve equivalent nerve responses, and accounting for the consequences of these factors is difficult. We describe the implementation, validation, and application of a standardized, modular, and scalable computational modeling pipeline for biophysical simulations of electrical activation and block of nerve fibers within peripheral nerves. The ASCENT (Automated Simulations to Characterize Electrical Nerve Thresholds) pipeline provides a suite of built-in capabilities for user control over the entire workflow, including libraries for parts to assemble electrodes, electrical properties of biological materials, previously published fiber models, and common stimulation waveforms. We validated the accuracy of ASCENT calculations, verified usability in beta release, and provide several compelling examples of ASCENT-implemented models. ASCENT will enable the reproducibility of simulation data, and it will be used as a component of integrated simulations with other models (e.g., organ system models), to interpret experimental results, and to design experimental and clinical interventions for the advancement of peripheral nerve stimulation therapies.

The delivery of kilohertz frequency alternating current (KHFAC) generates rapid, controlled, and reversible conduction block in motor, sensory, and autonomic nerves, but causes transient activation of action potentials at the onset of the blocking current. We implemented a novel engineering optimization approach to design blocking waveforms that eliminated the onset response by moving voltage-gated Na+ channels (VGSCs) to closed-state inactivation (CSI) without first opening. We used computational models and particle swarm optimization (PSO) to design a charge-balanced 10 kHz biphasic current waveform that produced conduction block without onset firing in peripheral axons at specific locations and with specific diameters. The results indicate that it is possible to achieve onset-free KHFAC nerve block by causing CSI of VGSCs. Our novel approach for designing blocking waveforms and the resulting waveform may have utility in clinical applications of conduction block of peripheral nerve hyperactivity, for example in pain and spasticity.