Specific frequency THz photons resonate Kv channel and decrease the AP firing rate in cortical neurons through molecular dynamics simulation. (a) Absorbance spectra of voltage-gated potassium/sodium ion channels and the bulk water. (b) The dynamic attributes of the Kv1.2 filter structure in pre- and post-exposure to HFTS. Purple balls represent the K+, blue balls represent the Cl-. (c) The alterations in potassium/sodium ion conductance consequent to the influence of HFTS. (d) Changes of the firing rate of APs of cortical neuron models before and after HFTS. (e) The FWHM of an AP pre- and post HFTS. (f) Changes in FWHM and firing frequency with or without HFTS. HFTS, high frequency terahertz stimulation. AP, action potential. FWHM, Full Wide of Half Maximum.

HFTS enhances Kv currents of pyramidal neurons in SNI mouse in vitro. (a) Anatomical location of ACC region in mice and a recorded PYR neuron (biocytin-labeled, green). (b) Representative Nav currents without (orange) or with HFTS (red) under the given step voltage protocol. (c) The application of HFTS does not change the activation and inactivation curves of Nav currents. (d-e) The corresponding slopes of the activation and inactivation curves (d) and the comparison of the half-activation and inactivation voltages (e). (f) Representative Kv currents evoked by a series of step voltages (inset) without (orange) or with HFTS (red). (g) I-V plots constructed from the values of traces shown in (f) display a stronger K+ current with HFTS. (SNI vs. SNI + HFTS: F(1, 66) = 67.97, P < 0.0001, nSNI = 7, nSNI+HFTS = 7; Two-way ANOVA followed by post hoc comparison using the Šídák’s multiple comparisons test (Supplement Table 1). (h) The application of HFTS does not change the activation curves of the Kv currents. (i-j) The corresponding slopes of the activation curves (i) and the half-activation voltages (j). SNI vs. SNI + HFTS: t = 5.872, P = 0.0011, n = 7, unpaired t-test. **, P < 0.01, ****, P < 0.0001.

HFTS reduces the APs firing rate of pyramidal neurons in SNI mice in vitro. (a-c) Representative traces (upper panels) and line-charts (lower panels) showing the changes of evoked spikes of pyramidal neurons in different groups Sham vs. SNI: F(1, 40) = 124.2, P < 0.001, nsham= 6, nSNI=6; SNI vs. SNI + HFTS: F(1, 40) = 23.13, P < 0.0001, nSNI=6, nSNI+HFTS=6; SNI vs. SNI + BLS: F(1, 40) = 0.1401, P = 0.7101, nSNI = 6, nSNI+BLS = 6. Two-way ANOVA followed by post hoc comparison using the Šídák’s multiple comparisons test (Supplement Table 2-4). (d) Superimposed traces showing the single AP evoked by threshold current stimulation in different groups. (e) Histograms showing the statistical comparison of rheobase in each group. (Sham vs. SNI: P < 0.001, nsham = 6, nSNI = 19; SNI vs. SNI + HFTS: P < 0.01, nSNI = 19, nSNI+HFTS = 6, one-way ANOVA followed by post hoc comparison using the Tukey’s multiple comparisons test). (f) The RMP in each group (Sham vs. SNI: q = 4.5, P < 0.05, nsham = 6, nSNI = 19; SNI vs. SNI + HFTS: P < 0.05, nSNI = 19, nSNI+HFTS= 6, one-way ANOVA followed by post hoc comparison using the Tukey’s multiple comparisons test). (g-i) HFTS has no significant effect on the threshold, amplitude and half-width of APs in pyramidal neurons. *, P < 0.05, **, P < 0.01, ***, P < 0.001, ****, P < 0.0001. BLS, blue light stimulation.

HFTS decreases the mean firing rate of pyramidal neurons in the ACC in both sham and SNI awake mice. (a) Schematic showing the timeline and the stimulating pattern of HFTS on an awake mouse, as well as the single-unit recording of the ACC using an in vivo multi-channel recording technique. (b) Example recording signals of ACC neurons before and after HFTS application in sham and SNI groups, respectively. (c) ACC neurons are classified as pyramidal (PYR) cells and interneurons (INT) using k-means cluster-separation algorithm based on their electrophysiological properties. (d) Histograms of the inter-spike intervals (ISI) from the spikes of a PYR and an INT in pre- and post-HFTS recording period. Insets at the top right corner show the waveforms of the detected single unit. (e) Pie charts summarize the changes in firing rate of PYR and INT in sham and SNI groups. Pre vs. post HFTS, Wilcoxon rank-sum test. (f) The mean firing rate of all recorded neurons in sham and SNI groups before and after HFTS. Sham (P < 0.0001, n = 73, paired t-test,), SNI (P < 0.0001, n = 130, paired t-test). (g) The mean firing rate of PYR neurons in sham and SNI groups before and after HFTS. Sham group (P < 0.0001, n = 93, paired t-test), SNI group (P < 0.0001, n = 123, paired t-test). (h) The mean firing rate of INT neurons in sham and SNI groups before and after HFTS. Sham group (P = 0.065, n = 15, paired t-test), SNI group (P = 0.0115, n = 7, paired t-test). *, P < 0.05, ****, P < 0.0001.

HFTS alleviates neuropathic pain of SNI mice through pain behavior tests. (a) Schematic of the establishment of NP model, the application of HFTS in ACC region and the following behavior tests including Von Frey test and Catwalk analysis. (b) HFTS increases the paw withdrawal mechanical thresholds (PWMTs) compared to the SNI model (F(18, 140) = 12.65. P < 0.0001. Sham vs. SNI: P < 0.0001; SNI vs. SNI + HFTS: P < 0.0001; n = 6 in each group. Two-way ANOVA repeated measures followed by post hoc comparison using the Šídák’s multiple comparisons test). (c-d) The print view of a mouse (c) and the step sequence of a sham mouse who passing through the glass pane (d). The red line represents the maximum intensity of each foot, the color box represents the mean intensity of the corresponding print during walking. (e) The feet view of the left front (LF), left hind (LH), right front (RF) and right hind (RH) in the groups of sham, SNI, SNI + HFTS and SNI + BLS, respectively. (f) HFTS increases the LH stand time of SNI mice (sham vs. SNI: P < 0.0001; SNI vs. SNI + HFTS: P < 0.01). (g) HFTS increases the LH stand index of SNI mice (sham vs. SNI: P < 0.001; SNI vs. SNI + HFTS: P < 0.05). (h) HFTS increases the LH max contact area of SNI mice (sham vs. SNI: P < 0.0001; SNI vs. SNI + HFTS: P < 0.05). (i) HFTS increases the LH mean print area of SNI mice (sham vs. SNI: P < 0.0001; SNI vs. SNI + HFTS: P < 0.05; SNI vs. SNI + BLS: P < 0.05). (j) HFTS increases the LH mean intensity of SNI mice (sham vs. SNI: P < 0.05; SNI vs. SNI + HFTS: P < 0.05). (k) HFTS have no significant for the pain behavior parameter of the duty cycle. *, P < 0.05, **, P < 0.01, ***, P < 0.001, ****, P < 0.0001. One-way ANOVA (f-k) followed by post hoc comparison using the Tukey’s multiple comparisons test. nSham = 38, nSNI = 35, nSNI+HFTS= 34, nSNI+BLS = 9.

Schematic diagram shows the mechanism of HFTS in alleviating neuropathic pain. Left panel shows the group with BLS and the right panel shows the group with HFTS.