ASCENT (Automated Simulations to Characterize Electrical Nerve Thresholds): A pipeline for sample-specific computational modeling of electrical stimulation of peripheral nerves.

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.

Validated computational models predict vagus nerve stimulation thresholds in preclinical animals and humans.

Objective.We demonstrated how automated simulations to characterize electrical nerve thresholds, a recently published open-source software for modeling stimulation of peripheral nerves, can be applied to simulate accurately nerve responses to electrical stimulation.Approach.We simulated vagus nerve stimulation (VNS) for humans, pigs, and rats. We informed our models using histology from sample-specific or representative nerves, device design features (i.e. cuff, waveform), published material and tissue conductivities, and realistic fiber models.Main results.Despite large differences in nerve size, cuff geometry, and stimulation waveform, the models predicted accurate activation thresholds across species and myelinated fiber types. However, our C fiber model thresholds overestimated thresholds across pulse widths, suggesting that improved models of unmyelinated nerve fibers are needed. Our models of human VNS yielded accurate thresholds to activate laryngeal motor fibers and captured the inter-individual variability for both acute and chronic implants. For B fibers, our small-diameter fiber model underestimated threshold and saturation for pulse widths >0.25 ms. Our models of pig VNS consistently captured the range ofin vivothresholds across all measured nerve and physiological responses (i.e. heart rate, Aδ/B fibers, Aγfibers, electromyography, and Aαfibers). In rats, our smallest diameter myelinated fibers accurately predicted fast fiber thresholds across short and intermediate pulse widths; slow unmyelinated fiber thresholds overestimated thresholds across shorter pulse widths, but there was overlap for pulse widths >0.3 ms.Significance.We elevated standards for models of peripheral nerve stimulation in populations of models across species, which enabled us to model accurately nerve responses, demonstrate that individual-specific differences in nerve morphology produce variability in neural and physiological responses, and predict mechanisms of VNS therapeutic and side effects.

Measuring and modeling the effects of vagus nerve stimulation on heart rate and laryngeal muscles.

BACKGROUND: Reduced heart rate (HR) during vagus nerve stimulation (VNS) is associated with therapy for heart failure, but stimulation frequency and amplitude are limited by patient tolerance. An understanding of physiological responses to parameter adjustments would allow differential control of therapeutic and side effects. To investigate selective modulation of the physiological responses to VNS, we quantified the effects and interactions of parameter selection on two physiological outcomes: one related to therapy (reduced HR) and one related to side effects (laryngeal muscle EMG). METHODS: We applied a broad range of stimulation parameters (mean pulse rates (MPR), intra-burst frequencies, and amplitudes) to the vagus nerve of anesthetized mice. We leveraged the in vivo recordings to parameterize and validate computational models of HR and laryngeal muscle activity across amplitudes and temporal patterns of VNS. We constructed a finite element model of excitation of fibers within the mouse cervical vagus nerve. RESULTS: HR decreased with increased amplitude, increased MPR, and decreased intra-burst frequency. EMG increased with increased MPR. Preferential HR effects over laryngeal EMG effects required combined adjustments of amplitude and MPR. The model of HR responses highlighted contributions of ganglionic filtering to VNS-evoked changes in HR at high stimulation frequencies. Overlap in activation thresholds between small and large modeled fibers was consistent with the overlap in dynamic ranges of related physiological measures (HR and EMG). CONCLUSION: The present study provides insights into physiological responses to VNS required for informed parameter adjustment to modulate selectively therapeutic effects and side effects.

Asynchronous axonal firing patterns evoked via continuous subthreshold kilohertz stimulation.

Objective.Transcutaneous electrical stimulation of peripheral nerves is a common technique to assist or rehabilitate impaired muscle activation. However, conventional stimulation paradigms activate nerve fibers synchronously with action potentials time-locked with stimulation pulses. Such synchronous activation limits fine control of muscle force due to synchronized force twitches. Accordingly, we developed a subthreshold high-frequency stimulation waveform with the goal of activating axons asynchronously.Approach.We evaluated our waveform experimentally and through model simulations. During the experiment, we delivered continuous subthreshold pulses at frequencies of 16.67, 12.5, or 10 kHz transcutaneously to the median and ulnar nerves. We obtained high-density electromyographic (EMG) signals and fingertip forces to quantify the axonal activation patterns. We used a conventional 30 Hz stimulation waveform and the associated voluntary muscle activation for comparison. We modeled stimulation of biophysically realistic myelinated mammalian axons using a simplified volume conductor model to solve for extracellular electric potentials. We compared the firing properties under kHz and conventional 30 Hz stimulation.Main results.EMG activity evoked by kHz stimulation showed high entropy values similar to voluntary EMG activity, indicating asynchronous axon firing activity. In contrast, we observed low entropy values in EMG evoked by conventional 30 Hz stimulation. The muscle forces evoked by kHz stimulation also showed more stable force profiles across repeated trials compared with 30 Hz stimulation. Our simulation results provide direct evidence of asynchronous firing patterns across a population of axons in response to kHz frequency stimulation, while 30 Hz stimulation elicited synchronized time-locked responses across the population.Significance.We demonstrate that the continuous subthreshold high-frequency stimulation waveform can elicit asynchronous axon firing patterns, which can lead to finer control of muscle forces.

Fibers in smaller fascicles have lower activation thresholds with cuff electrodes due to thinner perineurium and smaller cross-sectional area.

Objective. In nerve stimulation therapies, fibers in larger fascicles generally have higher activation thresholds, but the mechanisms are not well understood. We implemented and analyzed computational models to uncover the effects of morphological parameters on activation thresholds.Approach. We implemented finite element models of human vagus nerve stimulation to quantify the effects of morphological parameters on thresholds in realistic nerves. We also implemented simplified models to isolate effects of perineurium thickness, endoneurium diameter, fiber diameter, and fascicle location on current density, potential distributions (Ve), and activation thresholds across cuff geometries and stimulation waveforms. UsingVefrom each finite element model, we simulated activation thresholds in biophysical cable models of mammalian axons.Main results. Perineurium thickness increases with fascicle diameter, and both thicker perineurium and larger endoneurial diameter contributed to higher activation thresholds via lower peak and broader longitudinal potentials. Thicker perineurium caused less current to enter the fascicle transversely, decreasing peakVe. Thicker perineurium also inhibited current from leaving the fascicle, causing more constant longitudinal current density, broadeningVe. With increasing endoneurial diameter, intrafascicular volume increased faster than surface area, thereby decreasing intrafascicular current density and peakVe. Additionally, larger fascicles have greater cross-sectional area, thereby facilitating longitudinal intrafascicular current flow and broadeningVe. A large neighboring fascicle could increase activation thresholds, and for a given fascicle, fiber diameter had the greatest effect on thresholds, followed by fascicle diameter, and lastly, fascicle location within the epineurium. The circumneural cuff elicited robust activation across the nerve, whereas a bipolar transverse cuff with small contacts delivering a pseudo-monophasic waveform enabled more selective activation across fiber diameters and locations.Significance. Our computational studies provide mechanistic understanding of neural responses across relevant morphological parameters of peripheral nerves, thereby informing rational design of effective therapies.

Spatially selective stimulation of the pig vagus nerve to modulate target effect versus side effect.

Electrical stimulation of the cervical vagus nerve using implanted electrodes (VNS) is FDA-approved for the treatment of drug-resistant epilepsy, treatment-resistant depression, and most recently, chronic ischemic stroke rehabilitation. However, VNS is critically limited by the unwanted stimulation of nearby neck muscles-a result of non-specific stimulation activating motor nerve fibers within the vagus. Prior studies suggested that precise placement of small epineural electrodes can modify VNS therapeutic effects, such as cardiac responses. However, it remains unclear if placement can alter the balance between intended effect and limiting side effect. We used an FDA investigational device exemption approved six-contact epineural cuff to deliver VNS in pigs and quantified how epineural electrode location impacts on- and off-target VNS activation. Detailed post-mortem histology was conducted to understand how the underlying neuroanatomy impacts observed functional responses. Here we report the discovery and characterization of clear neuroanatomy-dependent differences in threshold and saturation for responses related to both effect (change in heart rate) and side effect (neck muscle contractions). The histological and electrophysiological data were used to develop and validate subject-specific computation models of VNS, creating a well-grounded quantitative framework to optimize electrode location-specific activation of nerve fibers governing intended effect versus unwanted side effect.

Quantified Morphology of the Cervical and Subdiaphragmatic Vagus Nerves of Human, Pig, and Rat.

It is necessary to understand the morphology of the vagus nerve (VN) to design and deliver effective and selective vagus nerve stimulation (VNS) because nerve morphology influences fiber responses to electrical stimulation. Specifically, nerve diameter (and thus, electrode-fiber distance), fascicle diameter, fascicular organization, and perineurium thickness all significantly affect the responses of nerve fibers to electrical signals delivered through a cuff electrode. We quantified the morphology of cervical and subdiaphragmatic VNs in humans, pigs, and rats: effective nerve diameter, number of fascicles, effective fascicle diameters, proportions of endoneurial, perineurial, and epineurial tissues, and perineurium thickness. The human and pig VNs were comparable sizes (∼2 mm cervically; ∼1.6 mm subdiaphragmatically), while the rat nerves were ten times smaller. The pig nerves had ten times more fascicles-and the fascicles were smaller-than in human nerves (47 vs. 7 fascicles cervically; 38 vs. 5 fascicles subdiaphragmatically). Comparing the cervical to the subdiaphragmatic VNs, the nerves and fascicles were larger at the cervical level for all species and there were more fascicles for pigs. Human morphology generally exhibited greater variability across samples than pigs and rats. A prior study of human somatic nerves indicated that the ratio of perineurium thickness to fascicle diameter was approximately constant across fascicle diameters. However, our data found thicker human and pig VN perineurium than those prior data: the VNs had thicker perineurium for larger fascicles and thicker perineurium normalized by fascicle diameter for smaller fascicles. Understanding these differences in VN morphology between preclinical models and the clinical target, as well as the variability across individuals of a species, is essential for designing suitable cuff electrodes and stimulation parameters and for informing translation of preclinical results to clinical application to advance the therapeutic efficacy of VNS.

Empirically Based Guidelines for Selecting Vagus Nerve Stimulation Parameters in Epilepsy and Heart Failure.

Vagus nerve stimulation (VNS) is a promising therapy to treat patients with epilepsy and heart failure. Outcomes of preclinical studies and clinical trials indicate that the selection of stimulation parameters has a direct impact on therapeutic efficacy and patient tolerability, suggesting that both the efficacy and tolerability of VNS could potentially be improved with a change in stimulation parameters. In this review, the success of translating stimulation parameters for epilepsy and heart failure from preclinical studies in animal models to human use in the clinic is evaluated on the basis of patient outcomes and stimulation-induced side effects. Data suggest that patients receiving VNS for epilepsy may experience improved seizure reduction by increasing the frequency and/or duty cycle of stimulation as well as incorporating closed-loop systems to deliver stimulation closer to seizure onset. Further, data suggest that VNS for heart failure is limited by the inability to activate the nerve fibers mediating therapeutic benefit without co-activation of side effect-inducing fibers. This may explain why pivotal trials of VNS for heart failure failed to meet primary efficacy outcomes despite promising preclinical outcomes in animal models. Improved characterization of the relationship between the stimulation parameter space and recruitment of the underlying fiber populations will likely expand the use of VNS to treat a variety of diseases and also improve upon current understanding of the mechanisms of action underlying VNS.