Spinal cord stimulation (SCS) for pain is typically implemented in an open-loop manner using parameters that remain largely unchanged. To improve the overall efficacy and consistency of SCS, one closed-loop approach proposes to use evoked compound action potentials (ECAPs) recorded from the SCS lead(s) as a feedback control signal to guide parameter selection. The goal of this study was to use a computational modeling approach to investigate the source of these ECAP recordings and technical and physiological factors that affect their composition. We developed a computational model that coupled a finite element model of lower thoracic SCS with multicompartment models of sensory axons within the spinal cord. We used a reciprocity-based approach to calculate SCS-induced ECAPs recorded from the SCS lead. Our model ECAPs contained a triphasic, P1, N1, P2 morphology. The model P2-N1 amplitudes and conduction velocities agreed with previous experimental data from human subjects. Model results suggested that the ECAPs are dominated by the activation of axons with diameters 8.7-10.0 μm located in the dorsal aspect of the spinal cord. We also observed changes in the ECAP amplitude and shape due to the electrode location relative to the vertebrae and spinal cord. Our modeling results suggest that clinically effective SCS relies on the activation of numerous axons within a narrow fiber diameter range and that several factors affect the composition of the ECAP recordings. These results can improve how we interpret and implement these recordings in a potential closed-loop approach to SCS.
Read full abstract