Tissue Temperature Increases by a 10 kHz Spinal Cord Stimulation System: Phantom and Bioheat Model.

OBJECTIVE: A recently introduced Spinal Cord Stimulation (SCS) system operates at 10 kHz, faster than conventional SCS systems, resulting in significantly more power delivered to tissues. Using a SCS heat phantom and bioheat multi-physics model, we characterized tissue temperature increases by this 10 kHz system. We also evaluated its Implanted Pulse Generator (IPG) output compliance and the role of impedance in temperature increases. MATERIALS AND METHODS: The 10 kHz SCS system output was characterized under resistive loads (1-10 KΩ). Separately, fiber optic temperature probes quantified temperature increases (ΔTs) around the SCS lead in specially developed heat phantoms. The role of stimulation Level (1-7; ideal pulse peak-to-peak of 1-7mA) was considered, specifically in the context of stimulation current Root Mean Square (RMS). Data from the heat phantom were verified with the SCS heat-transfer models. A custom high-bandwidth stimulator provided 10 kHz pulses and sinusoidal stimulation for control experiments. RESULTS: The 10 kHz SCS system delivers 10 kHz biphasic pulses (30-20-30 µs). Voltage compliance was 15.6V. Even below voltage compliance, IPG bandwidth attenuated pulse waveform, limiting applied RMS. Temperature increased supralinearly with stimulation Level in a manner predicted by applied RMS. ΔT increases with Level and impedance until stimulator compliance was reached. Therefore, IPG bandwidth and compliance dampen peak heating. Nonetheless, temperature increases predicted by bioheat multi-physic models (ΔT = 0.64°C and 1.42°C respectively at Level 4 and 7 at the cervical segment; ΔT = 0.68°C and 1.72°C respectively at Level 4 and 7 at the thoracic spinal cord)-within ranges previously reported to effect neurophysiology. CONCLUSIONS: Heating of spinal tissues by this 10 kHz SCS system theoretically increases quickly with stimulation level and load impedance, while dampened by IPG pulse bandwidth and voltage compliance limitations. If validated in vivo as a mechanism of kHz SCS, bioheat models informed by IPG limitations allow prediction and optimization of temperature changes.

Electrophysiology equipment for reliable study of kHz electrical stimulation.

Characterizing the cellular targets of kHz (1-10 kHz) electrical stimulation remains a pressing topic in neuromodulation because expanding interest in clinical application of kHz stimulation has surpassed mechanistic understanding. The presumed cellular targets of brain stimulation do not respond to kHz frequencies according to conventional electrophysiology theory. Specifically, the low-pass characteristics of cell membranes are predicted to render kHz stimulation inert, especially given the use of limited-duty-cycle biphasic pulses. Precisely because kHz frequencies are considered supra-physiological, conventional instruments designed for neurophysiological studies such as stimulators, amplifiers and recording microelectrodes do not operate reliably at these high rates. Moreover, for pulsed waveforms, the signal frequency content is well above the pulse repetition rate. Thus, the very tools used to characterize the effects of kHz electrical stimulation may themselves be confounding factors. We illustrate custom equipment design that supports reliable electrophysiological recording during kHz-rate stimulation. Given the increased importance of kHz stimulation in clinical domains and compelling possibilities that mechanisms of actions may reflect yet undiscovered neurophysiological phenomena, attention to suitable performance of electrophysiological equipment is pivotal.

Quasi-static pipeline in electroconvulsive therapy computational modeling.

BACKGROUND: Computational models of current flow during Electroconvulsive Therapy (ECT) rely on the quasi-static assumption, yet tissue impedance during ECT may be frequency specific and change adaptively to local electric field intensity. OBJECTIVES: We systematically consider the application of the quasi-static pipeline to ECT under conditions where 1) static impedance is measured before ECT and 2) during ECT when dynamic impedance is measured. We propose an update to ECT modeling accounting for frequency-dependent impedance. METHODS: The frequency content on an ECT device output is analyzed. The ECT electrode-body impedance under low-current conditions is measured with an impedance analyzer. A framework for ECT modeling under quasi-static conditions based on a single device-specific frequency (e.g., 1 kHz) is proposed. RESULTS: Impedance using ECT electrodes under low-current is frequency dependent and subject specific, and can be approximated at >100 Hz with a subject-specific lumped parameter circuit model but at <100 Hz increased non-linearly. The ECT device uses a 2 µA 800 Hz test signal and reports a static impedance that approximate 1 kHz impedance. Combined with prior evidence suggesting that conductivity does not vary significantly across ECT output frequencies at high-currents (800-900 mA), we update the adaptive pipeline for ECT modeling centered at 1 kHz frequency. Based on individual MRI and adaptive skin properties, models match static impedance (at 2 µA) and dynamic impedance (at 900 mA) of four ECT subjects. CONCLUSIONS: By considering ECT modeling at a single representative frequency, ECT adaptive and non-adaptive modeling can be rationalized under a quasi-static pipeline.

Novel Evoked Synaptic Activity Potentials (ESAPs) Elicited by Spinal Cord Stimulation.

Spinal cord stimulation (SCS) evokes fast epidural evoked compound action potential (ECAP) that represent activity of dorsal column axons, but not necessarily a spinal circuit response. Using a multimodal approach, we identified and characterized a delayed and slower potential evoked by SCS that reflects synaptic activity within the spinal cord. Anesthetized female Sprague Dawley rats were implanted with an epidural SCS lead, epidural motor cortex stimulation electrodes, an epidural spinal cord recording lead, an intraspinal penetrating recording electrode array, and intramuscular electromyography (EMG) electrodes in the hindlimb and trunk. We stimulated the motor cortex or the epidural spinal cord and recorded epidural, intraspinal, and EMG responses. SCS pulses produced characteristic propagating ECAPs (composed of P1, N1, and P2 waves with latencies <2 ms) and an additional wave ("S1") starting after the N2. We verified the S1-wave was not a stimulation artifact and was not a reflection of hindlimb/trunk EMG. The S1-wave has a distinct stimulation-intensity dose response and spatial profile compared with ECAPs. 6-Cyano-7-nitroquinoxaline-2,3-dione (CNQX; a selective competitive antagonist of AMPA receptors (AMPARs)] significantly diminished the S1-wave, but not ECAPs. Furthermore, cortical stimulation, which did not evoke ECAPs, produced epidurally detectable and CNQX-sensitive responses at the same spinal sites, confirming epidural recording of an evoked synaptic response. Finally, applying 50-Hz SCS resulted in dampening of S1-wave but not ECAPs. Therefore, we hypothesize that the S1-wave is synaptic in origin, and we term the S1-wave type responses: evoked synaptic activity potentials (ESAPs). The identification and characterization of epidurally recorded ESAPs from the dorsal horn may elucidate SCS mechanisms.

Effects of direct current stimulation on synaptic plasticity in a single neuron.

BACKGROUND: Transcranial direct current stimulation (DCS) has lasting effects that may be explained by a boost in synaptic long-term potentiation (LTP). We hypothesized that this boost is the result of a modulation of somatic spiking in the postsynaptic neuron, as opposed to indirect network effects. To test this directly we record somatic spiking in a postsynaptic neuron during LTP induction with concurrent DCS. METHODS: We performed rodent in-vitro patch-clamp recordings at the soma of individual CA1 pyramidal neurons. LTP was induced with theta-burst stimulation (TBS) applied concurrently with DCS. To test the causal role of somatic polarization, we manipulated polarization via current injections. We also used a computational multi-compartment neuron model that captures the effect of electric fields on membrane polarization and activity-dependent synaptic plasticity. RESULTS: TBS-induced LTP was enhanced when paired with anodal DCS as well as depolarizing current injections. In both cases, somatic spiking during the TBS was increased, suggesting that evoked somatic activity is the primary factor affecting LTP modulation. However, the boost of LTP with DCS was less than expected given the increase in spiking activity alone. In some cells, we also observed DCS-induced spiking, suggesting DCS also modulates LTP via induced network activity. The computational model reproduces these results and suggests that they are driven by both direct changes in postsynaptic spiking and indirect changes due to network activity. CONCLUSION: DCS enhances synaptic plasticity by increasing postsynaptic somatic spiking, but we also find that an increase in network activity may boost but also limit this enhancement.

Adaptive current-flow models of ECT: Explaining individual static impedance, dynamic impedance, and brain current density.

BACKGROUND: Improvements in electroconvulsive therapy (ECT) outcomes have followed refinement in device electrical output and electrode montage. The physical properties of the ECT stimulus, together with those of the patient's head, determine the impedances measured by the device and govern current delivery to the brain and ECT outcomes. OBJECTIVE: However, the precise relations among physical properties of the stimulus, patient head anatomy, and patient-specific impedance to the passage of current are long-standing questions in ECT research and practice. To this end, we develop a computational framework based on diverse clinical data sets. METHODS: We developed anatomical MRI-derived models of transcranial electrical stimulation (tES) that included changes in tissue conductivity due to local electrical current flow. These "adaptive" models simulate ECT both during therapeutic stimulation using high current (∼1 A) and when dynamic impedance is measured, as well as prior to stimulation when low current (∼1 mA) is used to measure static impedance. We modeled two scalp layers: a superficial scalp layer with adaptive conductivity that increases with electric field up to a subject-specific maximum (σSS¯), and a deep scalp layer with a subject-specific fixed conductivity (σDS). RESULTS: We demonstrated that variation in these scalp parameters may explain clinical data on subject-specific static impedance and dynamic impedance, their imperfect correlation across subjects, their relationships to seizure threshold, and the role of head anatomy. Adaptive tES models demonstrated that current flow changes local tissue conductivity which in turn shapes current delivery to the brain in a manner not accounted for in fixed tissue conductivity models. CONCLUSIONS: Our predictions that variation in individual skin properties, rather than other aspects of anatomy, largely govern the relationship between static impedance, dynamic impedance, and ECT current delivery to the brain, themselves depend on assumptions about tissue properties. Broadly, our novel modeling pipeline opens the door to explore how adaptive-scalp conductivity may impact transcutaneous electrical stimulation (tES).