Abstract
OBJECTIVES: To develop the first high-resolution, multi-scale model of cervical non-invasive vagus nerve stimulation (nVNS) and to predict vagus fiber type activation, given clinically relevant rheobase thresholds. METHODS: An MRI-derived Finite Element Method (FEM) model was developed to accurately simulate key macroscopic (e.g., skin, soft tissue, muscle) and mesoscopic (cervical enlargement, vertebral arch and foramen, cerebral spinal fluid [CSF], nerve sheath) tissue components to predict extracellular potential, electric field (E-Field), and activating function along the vagus nerve. Microscopic scale biophysical models of axons were developed to compare axons of varying size (Aα-, Aβ- and Aδ-, B-, and C-fibers). Rheobase threshold estimates were based on a step function waveform. RESULTS: Macro-scale accuracy was found to determine E-Field magnitudes around the vagus nerve, while meso-scale precision determined E-field changes (activating function). Mesoscopic anatomical details that capture vagus nerve passage through a changing tissue environment (e.g., bone to soft tissue) profoundly enhanced predicted axon sensitivity while encapsulation in homogenous tissue (e.g., nerve sheath) dulled axon sensitivity to nVNS. CONCLUSIONS: These findings indicate that realistic and precise modeling at both macroscopic and mesoscopic scales are needed for quantitative predictions of vagus nerve activation. Based on this approach, we predict conventional cervical nVNS protocols can activate A- and B- but not C-fibers. Our state-of-the-art implementation across scales is equally valuable for models of spinal cord stimulation, cortex/deep brain stimulation, and other peripheral/cranial nerve models.