Re: 전두엽 세타파(4–7Hz)의 증가 .. 만성통증치료 효과.. 2024 네이처!! VNS는?

[문형철]

1. nature  
2. scientific reports  
3. articles  
4. article

A novel theta-controlled vibrotactile brain–computer interface to treat chronic pain: a pilot study

Download PDF

• Article
• Open access
• Published: 10 February 2024

A novel theta-controlled vibrotactile brain–computer interface to treat chronic pain: a pilot study

• Phillip Demarest, 
• Nabi Rustamov, 
• James Swift, 
• Tao Xie, 
• Markus Adamek, 
• Hohyun Cho, 
• Elizabeth Wilson, 
• Zhuangyu Han, 
• Alexander Belsten, 
• Nicholas Luczak, 
• Peter Brunner, 
• Simon Haroutounian & 
• Eric C. Leuthardt 

Scientific Reports volume 14, Article number: 3433 (2024) Cite this article

• 6329 Accesses

• 17 Citations

• 9 Altmetric

• Metricsdetails

Abstract

Limitations in chronic pain therapies necessitate novel interventions that are effective, accessible, and safe. Brain–computer interfaces (BCIs) provide a promising modality for targeting neuropathology underlying chronic pain by converting recorded neural activity into perceivable outputs. Recent evidence suggests that increased frontal theta power (4–7 Hz) reflects pain relief from chronic and acute pain. Further studies have suggested that vibrotactile stimulation decreases pain intensity in experimental and clinical models. This longitudinal, non-randomized, open-label pilot study's objective was to reinforce frontal theta activity in six patients with chronic upper extremity pain using a novel vibrotactile neurofeedback BCI system. Patients increased their BCI performance, reflecting thought-driven control of neurofeedback, and showed a significant decrease in pain severity (1.29 ± 0.25 MAD, p = 0.03, q = 0.05) and pain interference (1.79 ± 1.10 MAD p = 0.03, q = 0.05) scores without any adverse events. Pain relief significantly correlated with frontal theta modulation. These findings highlight the potential of BCI-mediated cortico-sensory coupling of frontal theta with vibrotactile stimulation for alleviating chronic pain.

초록

만성 통증 치료의 한계는

효과적이고 접근 가능하며 안전한 새로운 중재법을 필요로 합니다.

뇌-컴퓨터 인터페이스(BCI)는

기록된 신경 활동을 인지 가능한 출력으로 변환함으로써

만성 통증의 근본적인 신경병리를 표적으로 삼는 유망한 방식을 제공합니다.

최근 증거는

전두엽 세타파(4–7Hz)의 증가가

만성 및 급성 통증의 완화를 반영함을 시사합니다.

추가 연구들은 진동 촉각 자극이 실험적 및 임상 모델에서 통증 강도를 감소시킨다고 제안했습니다. 본 종단적, 비무작위, 공개 라벨 파일럿 연구의 목적은 새로운 진동 촉각 신경 피드백 BCI 시스템을 사용하여 만성 상지 통증 환자 6명에서 전두엽 세타 활동을 강화하는 것이었다. 환자들은 신경 피드백에 대한 사고 기반 제어를 반영하여 BCI 수행 능력을 향상시켰으며, 통증 심각도(1.29±0.25 MAD, p=0.03, q=0.05) 및 통증 간섭(1.79±1.10 MAD, p=0.03, q=0.05) 점수가 유의하게 감소했으며, 어떠한 부작용도 발생하지 않았다.

통증 완화는

전두엽 세타 조절과 유의한 상관관계를 보였다.

이러한 결과는

진동 촉각 자극을 통한 BCI 매개 전두엽 세타의 피질-감각 결합이

만성 통증 완화에 잠재력을 지닌다는 점을 강조한다.

Similar content being viewed by others

First-in-human prediction of chronic pain state using intracranial neural biomarkers

Article 22 May 2023

The modulation of neural insular activity by a brain computer interface differentially affects pain discrimination

Article Open access07 May 2021

An externally validated resting-state brain connectivity signature of pain-related learning

Article Open access17 July 2024

Introduction

Currently, available approaches for treating chronic pain have limited therapeutic benefits and are often associated with side effects. Only 30–40% of patients with chronic pain achieve meaningful pain relief with pharmacotherapy1,2. Many of these patients receive opioids, which are associated with increased risks of dependence and substance use disorders3,4. Invasive approaches for treating chronic pain, such as deep-brain stimulation and spinal cord stimulation, are effective in some patients but are expensive and can lead to complications5,6. Given the current limitations in treating chronic pain, there is a critical need to develop novel, safe, non-invasive, and affordable therapies.

서론

현재 이용 가능한 만성 통증 치료법은

치료 효과가 제한적이며 종종 부작용과 연관됩니다.

만성 통증 환자의 30~40%만이

약물 치료로 의미 있는 통증 완화를 달성합니다1,2.

이들 중 다수는 의존성 및 물질 사용 장애 위험 증가와 연관된 오피오이드를 투여받습니다3,4.

심부 뇌 자극 및 척수 자극과 같은 만성 통증 치료를 위한 침습적 접근법은

일부 환자에게 효과적이지만

비용이 많이 들고 합병증을 유발할 수 있다5,6.

만성 통증 치료의 현재 한계를 고려할 때,

새롭고 안전하며 비침습적이고 경제적인 치료법을 개발할 필요가 절실하다.

Chronic pain is associated with pathological changes in neural circuitry, including maladaptive neuroplasticity, cortical reorganization, and changes to descending pain modulation pathways7,8,9,10. In humans, brain–computer interface (BCI) systems have emerged as a promising tool for remodeling neural circuits11. This premise of BCI relies on Hebbian learning—i.e., activity-dependent synaptic plasticity where concurrent neuronal activations strengthen neural connections. BCIs can decode cortical physiology in real-time and provide temporally precise outputs to enable the necessary conditions for Hebbian learning12. This physiologic and temporal precision enables powerful conditions for enhancing neuroplastic changes. Specifically, BCIs enable robust cortico-sensory coupling, where a desired brain physiology is amplified with functionally relevant sensory feedback that is provided in real-time and with high temporal precision13. This BCI-mediated method has enabled novel therapies for what have been considered intractable neurologic conditions, such as chronic stroke-induced hemiparesis14,15,16. Early efforts in using BCI as a pain intervention have been attempted with mixed results and rarely report relationships between physiology and outcome17,18,19,20. This variable effect is partly due to the brain signals used and the feedback provided. Specifically, past BCIs have commonly used cortical sensorimotor rhythms associated with motor imagery and visual feedback17,18,20,21,22. Neither of these inputs and outputs are directly relevant to pain perception. Thus, in the setting of chronic pain, it will be essential to couple cortical signals associated with pain relief with pain-alleviating sensory input, such as vibrotactile stimulation23,24,25, to best enable neural remodeling to occur. Further, the variation in brain regions previously reported in pain processing justifies targeted reinforcement of physiology consistently observed in multiple contexts26. The current work builds on the previous demonstrations that relief of both chronic pain and experimental acute pain is associated with increased frontal theta (θ) rhythms (4–7 Hz)27,28.

만성 통증은

부적응적 신경가소성,

피질 재조직화,

하행성 통증 조절 경로의 변화 등

신경 회로의 병리적 변화와 연관되어 있다7,8,9,10.

maladaptive neuroplasticity,

cortical reorganization, and

changes to descending pain modulation pathways

인간에서 뇌-컴퓨터 인터페이스(BCI) 시스템은

신경 회로 재구성을 위한 유망한 도구로 부상했습니다11.

BCI의 이러한 전제는

헤비안 학습,

즉 동시 신경 세포 활성화가 신경 연결을 강화하는 활동 의존적 시냅스 가소성에 기반합니다.

BCI는

실시간으로 피질 생리학을 해독하고 시간적으로 정밀한 출력을 제공하여

헤비안 학습에 필요한 조건을 가능하게 합니다12.

이러한 생리학적·시간적 정밀도는

신경가소성 변화를 증진시키는 강력한 조건을 제공합니다.

특히 BCI는

강력한 대뇌피질-감각 결합(cortico-sensory coupling)을 가능하게 하는데,

이는 기능적으로 관련 있는 감각 피드백을 실시간으로 높은 시간적 정밀도로 제공함으로써

원하는 뇌 생리학을 증폭시킵니다13.

이러한 BCI 매개 방식은

만성 뇌졸중으로 인한 편마비14,15,16과 같이 난치성 신경학적 질환으로

여겨져 온 상태에 대한 새로운 치료법을 가능하게 했습니다.

통증 개입 수단으로 BCI를 활용한 초기 시도들은 혼합된 결과를 보였으며,

생리학적 상태와 치료 결과 간의 연관성을 보고한 사례는 드물다17,18,19,20.

이러한 효과의 변동성은

부분적으로 사용된 뇌 신호와 제공된 피드백에 기인한다.

구체적으로,

기존 BCI들은 주로 운동 상상과 시각적 피드백과 연관된 피질 감각운동 리듬을 활용해왔다17,18,20,21,22.

이러한 입력과 출력 모두 통증 지각과 직접적으로 관련이 없습니다.

따라서

만성 통증 환경에서는 통증 완화와 관련된 피질 신호를

진동 촉각 자극23,24,25과 같은 통증 완화 감각 입력과 결합하여

신경 재구성이 최적으로 이루어지도록 하는 것이 필수적입니다.

이 무작위 이중맹검 위약 대조 교차 연구에서, 우리는 말초 신경 차단을 통해 절단 후 발생하는 유령통 및 절단부 통증을 일시적으로 제거할 수 있는지 조사하였다. 지속적인 절단 후 통증을 가진 절단 환자를 포함시켜, 교차 설계로 리도카인 2%에 아드레날린을 첨가한 신경 차단액 또는 생리식염수를 투여받도록 무작위 배정하였다. 기저선 및 리도카인 또는 생리식염수 주사 후 120분까지 고정된 시간점에 자발적 환부통 및 절단부 통증과 유발 반응을 평가하였다. 주요 결과는 기저선과 리도카인 또는 생리식염수 주사 30분 후의 최악 통증 강도(환부통 또는 절단부 통증) 간 절대 변화량의 차이였다. 12명의 절단 환자가 무작위 배정되었으며 9명의 환자가 분석에 포함되었다. 리도카인 주사와 식염수 주사 간 최악의 통증 강도 중앙값의 절대적 변화는 −2.0(사분위수 범위, −4.0~0.0)이었다(n = 9, P = 0.12). 9명의 환자 중 9명이 리도카인 주사 후 어느 정도 통증 완화를 보고한 반면, 식염수 주사 후에는 9명 중 2명만이 통증 완화를 보고했다(P = 0.02). 리도카인 주사 후 유령통 강도는 식염수 주사에 비해 유의하게 감소했으나(P = 0.04), 두 중재 간 절단부 통증 강도에는 유의한 차이가 없었다(P = 0.17). 9명의 절단 환자 모두에서 리도카인 주사 후 유발 반응이 소실되었다.

따라서
본 연구 결과는
말초 신경계의 구심성 입력이 절단 후 통증에 중요한 역할을 함을 시사한다.

요약
신경병성 통증에서 구심성 감각 입력의 역할을 두 그룹에서 조사하였다. 말초 신경 차단은 모든 환자에서 자발적 및 유발 통증을 소멸시켰다. 전신 리도카인은 신경 손상보다 다발성 신경병증으로 인한 통증에 더 효과적이었다. 중추 감작은 자율적인 자발적 통증 생성 기전으로 작용하지 않는다. 말초 신경 손상 후 통증 유지에는 말초 감각 입력이 중요하다.

말초 신경 손상 후 중추 감작은 척수 후각에서 이소성 신경 세포 활동을 유발할 수 있으며, 이는 잠재적인 자율적 통증 생성 기전을 시사한다. 본 연구는 말초 신경 차단 및 전신 리도카인 투여와 함께 상세한 체성감각 평가를 통해 말초 신경병성 통증 유지에 일차 구심성 입력의 기여도를 확인하였다. 신경병성 통증 환자 14명(말초 신경 손상으로 인한 일측성 발 통증 환자 7명, 말초 다발신경병증으로 인한 양측성 발 통증 환자 7명)을 대상으로 정량적 감각 검사를 통해 체성감각 기능을 종합적으로 평가하였다. 이후 환자들에게 초음파 유도 말초 신경 차단술과 리도카인 정맥 주입을 무작위 순서로 시행하였다. 각 세션에서 이러한 중재가 자발적 통증 강도 및 유발된 냉각, 온열, 바늘자극, 붓자극 반응에 미치는 효과를 평가하였다. 모든 환자는 기준선에서 감각 장애를 보였다. 말초 신경 차단은 모든 환자에서 10분(중앙값) 이내에 동측 통증을 완전히 소멸시켰으며, 리도카인 혈장 농도는 약물의 전신 효과를 설명하기에는 너무 낮았다. 정맥 내 리도카인 주입은 자발적 통증을 45.5%(±31.7%) 감소시켰으며, 기계적 및 열적 과민증을 보인 대부분의 환자에서 이러한 증상을 완화시켰다. 그러나 유발성 과민성의 개선은 자발적 통증 강도에 대한 약물 효과와 관련이 없었다. 본 연구는 개별 체성감각 표현형 및 중추 감작 징후와 무관하게, 말초 신경 손상 및 말초 다발성 신경병증에서 일차 감각신경 입력은 신경병성 통증 유지에 결정적임을 입증하였다.

또한

통증 처리 과정에서 이전에 보고된 뇌 영역의 변이는

다양한 맥락에서 일관되게 관찰되는 생리학적 반응을

표적 강화할 필요성을 정당화합니다26.

본 연구는

만성 통증과 실험적 급성 통증 완화가

모두 전두엽 세타(θ) 리듬(4–7Hz) 증가와 연관된다는

기존 연구 결과27,28을 기반으로 한다.

Increases in frontal θ have been previously reported during cognitive tasks involving attention and concentration (such as working memory tasks, mental arithmetic, and meditation)29,30,31,32. Some of these tasks have previously been shown to have implications for pain management, such as the positive clinical outcomes associated with meditation or the phenomenon of distraction-induced analgesia32,33,34,35,36,37. Recent reviews of studies that have elucidated the neural correlates of acute and chronic pain have revealed divergent evidence. Studies that reported changes in θ observed variable results, with some reporting increases, decreases, or no changes. The locations of reported changes in θ were also highly variable, with most studies only reporting changes in global activity38,39. Further, these studies are primarily associative and are unable to decouple biomarkers of pathophysiology (i.e., active contributors of pathology) from biomarkers of homeostatic responses (i.e., compensatory mechanisms that reduce pathology). This limitation of previous studies prevents inferences about how a particular EEG activity contributes to symptoms of pain and only provides evidence of co-occurrence. As such, approaches that reinforce pain-related biomarkers and subsequently observe their effect on pain-related symptoms are crucial for interpreting the relationship between EEG features and pain symptoms.

전두엽 세타 증가 현상은

주의력과 집중력을 요하는 인지 과제(작업 기억 과제, 정신 산술, 명상 등) 수행 중에도 보고된 바 있다29,30,31,32.

이러한 과제 중 일부는

명상과 연관된 긍정적 임상 결과나 주의 분산에 의한 진통 현상32,33,34,35,36,37과 같이

통증 관리에 시사점을 제공하는 것으로 이전에 입증된 바 있다.

급성 및 만성 통증의 신경적 상관관계를 규명한 연구에 대한 최근의 검토는 상이한 증거를 드러냈다. θ

변화에 대해 보고한 연구들은 증가, 감소 또는 변화 없음 등 다양한 결과를 보였다. 보고된 θ파 변화의 위치 역시 매우 다양했으며, 대부분의 연구는 전역적 활동 변화만을 보고했다38,39. 또한, 이러한 연구들은 주로 연관성 연구이며 병리생리학적 생체표지자(즉, 병리의 능동적 기여 요소)와 항상성 반응 생체표지자(즉, 병리를 감소시키는 보상 기전)를 분리해낼 수 없다. 이러한 기존 연구의 한계로 인해 특정 EEG 활동이 통증 증상에 어떻게 기여하는지에 대한 추론은 불가능하며, 단지 동시 발생의 증거만 제공할 뿐이다. 따라서 통증 관련 생체표지자를 강화하고 그 후 통증 관련 증상에 미치는 영향을 관찰하는 접근법은 EEG 특징과 통증 증상 간의 관계를 해석하는 데 매우 중요하다.

To this end, we created a novel BCI system that is controlled by frontal θ rhythms and uses vibrotactile stimulation of the affected area as a means of neurofeedback in patients with medically refractory chronic upper extremity pain. In order to justify placebo-controlled clinical trials that leverage this novel electrophysiological target for neuromodulation in the future, the present study aimed to test the feasibility of BCI-mediated reinforcement of frontal θ rhythms, and whether this physiology can be linked to symptom relief in chronic pain patients. Furthermore, using secondary analysis, we aimed to identify mechanistic evidence related to BCI therapy that warrants future elucidation and study design refinement. We hypothesized that a longitudinal intervention with such a BCI-controlled vibrotactile system would be feasible, safe, and would reduce pain severity and pain interference as measured by the Brief Pain Inventory (BPI). Therefore, the objective of this pilot trial was to test the feasibility and determine the initial efficacy of treating chronic upper extremity pain with a BCI-controlled vibrotactile system and identify potential neurophysiological correlates of the observed clinical outcomes. Taken together, the findings in this work provide early compelling evidence that BCI-mediated cortico-sensory coupling may provide a novel approach to alleviating chronic pain in the upper extremity and justify conducting randomized, controlled clinical trials in the future.

이를 위해 본 연구는 의학적 치료에 반응하지 않는 만성 상지 통증 환자를 대상으로 전두엽 θ 리듬에 의해 제어되며, 신경피드백 수단으로 통증 부위의 진동 촉각 자극을 활용하는 새로운 BCI 시스템을 개발하였다. 향후 이 새로운 신경조절 전기생리학적 표적을 활용한 위약 대조 임상시험의 타당성을 입증하기 위해, 본 연구는 BCI 매개 전두엽 θ 리듬 강화의 실행 가능성과 이 생리학적 변화가 만성 통증 환자의 증상 완화와 연관될 수 있는지 검증하는 것을 목표로 하였다. 또한 부차적 분석을 통해 향후 규명과 연구 설계 개선이 필요한 BCI 치료 관련 기전적 증거를 규명하고자 하였다. 우리는 이러한 BCI 제어 진동 촉각 시스템을 이용한 장기적 개입이 실행 가능하고 안전하며, 간이 통증 평가표(BPI)로 측정된 통증 심각도와 통증 간섭을 감소시킬 것이라는 가설을 세웠다. 따라서 본 파일럿 시험의 목적은 BCI 제어 진동 촉각 시스템으로 만성 상지 통증을 치료하는 실행 가능성과 초기 효능을 검증하고, 관찰된 임상 결과의 잠재적 신경생리학적 상관관계를 규명하는 것이었다. 종합하면, 본 연구 결과는 BCI 매개 대뇌-감각 결합이 상지 만성 통증 완화를 위한 새로운 접근법을 제공할 수 있다는 초기 설득력 있는 증거를 제시하며, 향후 무작위 대조 임상시험 수행의 타당성을 입증한다.

Results

Patient characteristics

This study included six patients diagnosed with upper extremity chronic pain who met the inclusion criteria (3 females, three males; age: 59.5 ± 5.5 [median ± median absolute deviation (MAD)]; range 18–74 years). The median pre-BCI intervention Visual Analog Scale (VAS) pain rating was 47.5 (± 18.5 MAD). Patient demographics, characteristics, and comorbidities are shown in Table 1. Participant Hospital Anxiety and Depression Scale (HADS) score was 5.5 (± 2.0 MAD) for depression and 5.00 (± 2.50 MAD) for anxiety. Two of the six patients had depression and anxiety scores above the normal 0–7 range. The total score for Pain Catastrophizing Scale (PCS) was 19.5 (± 13.0 MAD). Table 2 shows the HADS and PCS results for all patients. Four of the six patients had hypersensitivity in their affected area. One had heat and pinprick hypersensitivity, two had pinprick hypersensitivity, and one had brush and heat hypersensitivity. The location of the pain area and type of hypersensitivity for each patient are shown in Supp. Fig. 1. Three of the patient’s affected areas were bilateral, with the average reported pain intensity being the same on both sides. Hypersensitivity, Quantitative Sensory Testing (QST), and conditioned pain modulation (CPM) tests were conducted on the affected side with greater pain intensity on the day of the baseline visit. The results of QST and CPM are reported in Supp. Tables 1 and 2, respectively. The baseline median scores on the BPI, the primary outcome measure, were 5.25 (± 1.63 MAD) for the Pain Severity Score (PSS), and 4.86 (± 2.29 MAD) for Pain Interference Score (PIS). The median Neuropathic Pain Symptom Inventory (NPSI) total score before BCI intervention was 17.9 (± 9.60 MAD).

Table 1 Patient demographics and medical history.

Full size table

Table 2 Patient HADS and pain catastrophizing scale scores.

Full size table

Intervention retention, adherence, safety, and comfort

Patients underwent three to five one-half to 1-h BCI training sessions per week over 5–6 weeks. The intervention duration was specified to be at least 5 weeks and no more than 6 weeks. The BCI system consisted of a dry electrode electroencephalogram (EEG) headset, a laptop, a custom-made hand vibration stimulation array (HVSA) for functionally relevant sensory feedback, and a video screen for visual feedback (Fig. 1a,b). Figure 1c shows all electrode locations used during recordings, and the integration of the EEG system and the HVSA, where vibrotactile stimulation is delivered to the affected area as a function of patient-generated frontal θ power.

Figure 1

Schematic and description of BCI system. (a) System overview schematic of EEG-driven BCI system consisting of EEG signals as input and dual visual and vibrotactile sensory feedback. Electrophysiological activity is recorded using a 24-channel dry EEG headset connected to a PC laptop via Bluetooth. The PC uses BCI2000 for real-time signal processing (i.e., extrapolating frontal midline θ power changes) and command execution (i.e., initiating visual and tactile feedback in response to changes in frontal θ power). Visual feedback was delivered on an external monitor, and tactile feedback was delivered using the custom HVSA. (b) System set up with BCI system with vibrotactile feedback system enclosing the affected hand. (c) Overview of electrode configuration (Top-left) and vibrotactile neurofeedback. Increased F3 θ modulation (teal) during relevant phases of BCI therapy leads to vibrotactile feedback of the affected area (green). No modulation of θ or decreased modulation of θ leads to no vibrotactile feedback (red).

Full size image

The overall design of the intervention is shown in Fig. 2a. Patients meeting inclusion and exclusion criteria attended an initial baseline visit to determine their pain symptom characteristics, record baseline EEG, and perform an initial assessment of task-mediated frontal θ modulation (Supplemental Methods)30,38. VAS pain scores were obtained before and after BCI training (main intervention) in each of these sessions. BPI, consisting of pain severity score and pain interference score (PSS and PIS, respectively), and NPSI metrics were obtained once per week via survey after all of each week’s training sessions. A specific breakdown of the structure of BCI training sessions is shown in Fig. 2b.

Figure 2

Overview of BCI therapy. (a) Overall treatment timeline. Solid black shapes indicate the main therapy sections. White shapes indicate specific breakdowns of each main section. (b) Overview of a single BCI trial. The teal box shows the numerical stimulus used in early BCI sessions to induce θ modulation. “Start” and “Stop” cues are indicated by red arrows. Cursor movement is indicated by the orange dotted arrow. Trial states progress from left to right. θ modulation direction leads to movement of the cursor during neurofeedback. The neurofeedback period was the only period with a variable duration, with 20 s being the maximum duration of the neurofeedback. If the cursor reaches the green target at the top of the screen, several pulses of vibration are delivered to the affected area over the duration of the “Stop” cue.

Full size image

A total of eight patients meeting the inclusion/exclusion criteria attended an initial baseline visit. One patient could not meet the retention criteria and withdrew from the study after the baseline visit (recruitment rate 87.5%, 7/8). One other patient could not adhere to the treatment regimen due to scheduling conflicts and withdrew from the study after the fourth week (retention 85.7%, 6/7). The regimens followed by each of the six patients who completed the pilot study are visualized in Fig. 3a. Five out of the six patients completed 6 weeks of BCI intervention. One patient completed 5 weeks. The retention and recruitment rates of this study are comparable to the rates found in larger-scale randomized controlled trials (median recruitment rate 72%, median retention rate 88%)40 and comparable studies using longitudinal therapeutic BCIs (recruitment range 68.4–80.5%, retention range 55.6–76.9%)15,21,41.

Figure 3

Intervention feasibility. (a) Each patient’s intervention time. Black bars represent the number of weeks each patient underwent BCI training. Gray bars represent the mean number of BCI training days per week. Error bars represent the range of days per week completed. All patients underwent at least 3 days of BCI training per week for 5 weeks. Four of the patients completed 3 days of BCI training every week. One patient completed 3 days of training per week for 5 weeks, then completed 2 days of training on the sixth week (range 2–3, where only two sessions were completed on the sixth week). One patient completed 4–5 days of training per week for 6 weeks (range 4–5). (b) Median BCI performance across all participants across 6 weeks of BCI therapy. The shaded region depicts the median absolute deviation. (c) Violin plots reflecting session duration across all patients each week. The dotted line represents the median session time across all patients each week. There was no significant difference in session duration between weeks (n = 19, 20, 19, 19, 19, 16; Kruskal–Wallis test, df = 5, χ2 = 8.28, p = 0.14). (d) Exemplar power spectral density plots comparing average θ range power during rest and neurofeedback phases of BCI training for patient 1. Shaded regions represent standard error. θ modulation increases over BCI training sessions. The red rectangle indicates frequencies used for BCI real-time signal processing.

Full size image

The median BCI performance across all 6 weeks of treatment is shown in Fig. 3b. Patients underwent a median BCI training session duration of 43 min, and completed a median of 1,192 BCI trials over the intervention period. There were no significant changes in weekly session duration throughout the pilot (n = 19, 20, 19, 19, 19, 16; Kruskal–Wallis test, df = 5, χ2 = 8.28, p = 0.14; Fig. 3c). Each patient’s daily session duration is visualized in Supp. Fig. 2. The total number of trials each patient underwent during each BCI training day and across the entirety of the intervention is shown in Supp. Table 3. Use of the BCI system did not lead to any apparent complications or discomfort for any of the patients. No patients reported any issues or additional symptoms throughout the course of the intervention.

Feasibility of frontal theta modulation

To confirm‎ the feasibility of our intervention in creating observable differences in task-mediated frontal (F3) θ power, patients performed a serial subtraction mental arithmetic task during their baseline visit, and differences in power spectra were analyzed. The results from this task are shown in Supp. Fig. 2. These data demonstrate that the serial subtraction task can create qualitative differences in θ power in patients; however, the effect sizes, specificity, and directionality vary. Figure 3d shows an exemplar progression of power spectral density in θ range for a single patient during BCI training for neurofeedback and rest periods at the F3 electrode. There is a clear qualitative separation of power in θ range between the neurofeedback period and baseline (rest) period of BCI training. Furthermore, the magnitude of the θ modulation during neurofeedback increases with each exemplar BCI training session. The comparison of differences in θ power between the pre-intervention baseline and BCI intervention is shown in Supp. Fig. 2. The normalized θ power change over each BCI session neurofeedback period across all patients is shown in Supp. Fig. 2. In sum, these results reflect the feasibility of our pre-screening task and BCI intervention in producing and detecting modulation in frontal (F3) θ power.

Clinical outcomes and BCI performance

At the group level, for the primary endpoint, there was a significant decrease in PSS and PIS (N = 6, Wilcoxon signed-rank, median decrease 1.29 ± 0.25 MAD, p = 0.03, q = 0.05 and 1.79 ± 1.10 MAD p = 0.03, q = 0.05, respectively, Fig. 4a and b, Supp. Table 4) after BCI intervention. For the secondary endpoint, the decrease in NPSI total score was not statistically significant (N = 6, Wilcoxon signed-rank, median decrease 5.25 ± 5 MAD p = 0.09, q = 0.09). Changes in NPSI and BPI throughout BCI intervention for each patient are visualized in Supp. Fig. 4. These results indicate improvements in pain severity and the extent to which pain interfered with day-to-day functions after participation in the BCI intervention. Changes in primary and secondary outcome (BPI, NPSI) for each patient before and after BCI intervention is shown in Supp. Table 4.

Figure 4

Vibrotactile BCI therapy decreases chronic pain symptoms. Individual patient outcomes: (a) PSS (N = 6, Wilcoxon signed-rank, median decrease 1.29 ± 0.25 MAD p = 0.03, q = 0.05) and (b) PIS (N = 6, Wilcoxon signed-rank, median decrease 1.79 ± 1.10 MAD p = 0.03, q = 0.05) significantly decreases after the BCI intervention duration. (c) Median BCI performance and VAS rating for all participants across the BCI therapy intervention duration. Across sessions, BCI performance improves while VAS pain rating decreases. The shaded region depicts the median absolute deviation (MAD). (d) The median baseline θ power over the course of BCI therapy at electrode location F3 across participants with either significant pain reduction or consistent BCI performance over the ROC calculated optimal threshold (n = 4). Topography plots show median θ power across all electrodes during the first week (upper) and fourth week (lower, where baseline F3 θ was highest). Data were z-scored to baseline θ power, with shaded regions representing MAD.

Full size image

There were observed differences in BCI performance trends between patients, and thus, it was expected that changes in pain ratings would be different between patients. It is therefore important to acknowledge the individual changes in pain ratings over time and the overall relationship between these trends and frontal θ modulation. Median pain ratings decreased as median BCI performance accuracy increased across all patients over the BCI intervention (Fig. 4c). The overall trends in BCI accuracy, BCI bit rate, and VAS pain ratings for each patient are shown in Supp. Fig. 3. VAS was used to sample pain ratings during independent BCI training sessions, providing a transient metric of pain that we relate to electrophysiology and BCI performance. This enabled additional analysis to be conducted to relate electrophysiology, BCI performance metrics, and transient measures of pain. On an individual level, three patients showed a significant decrease in pain over time (n = 18, 15, 27, Spearman’s ρ = − 0.60, − 0.75, − 0.74; p = 0.01, 4.43e−3, 1.37e−4, q = 0.02, 0.01, 8.20e−4). Three of the patients significantly increased BCI performance over the intervention period (n = 18, 15, 27, Spearman’s ρ = 0.49, 0.58, 0.43; p = 0.04, 0.03, 0.03, q = 0.07, 0.07, 0.07), and one patient began the intervention with a near-perfect performance which decreased mildly over time. These four patients all showed a decrease in pain ratings before and after treatment (Supp. Fig. 3). Notably, while four of the six participants were classified as either performers (i.e. patients who had a significant increase in BCI performance over time) or VAS pain responders (i.e. patients who had a significant decrease in VAS ratings over time), all six patients had an improvement in primary outcome measures. Taken together, these results suggest a relationship between BCI performance, theta modulation, and pain ratings.

Changes in power spectral density across BCI therapy runs

The study was designed to assess whether 6 weeks of BCI training led to significant resting-state frequency-band power changes. Patients classified as VAS pain responders and performers were those with significant pain decrease over time or with a significant increase in BCI performance over time, respectively (Supp. Fig. 3). For these patients, resting θ power at F3 increased over the course of therapy (Fig. 4d). Week four has the greatest median resting θ power across VAS pain responders and performers, which was notably the same week with the greatest BCI performance across all patients (Fig. 4c and d). There was no significant change in the resting power of any frequency band when comparing measures across BCI training sessions for all patients. Topography plots showing qualitative changes in resting power across all electrodes are shown in Supplementary Fig. 5.

Specificity of theta modulation during BCI training

In order to identify the specificity of θ modulation and investigate whether frequency-specific modulation was different between high BCI control and the absence of BCI control, we compared the progression of frequency-specific power during BCI training trials. Figure 5 shows the median progression of normalized θ power across BCI training trials. The normalized θ power, averaged across all BCI sessions, was compared between high BCI performance and the absence of BCI control. These results exemplify the qualitative difference in θ power modulation specificity between sessions where patients did not yet acquire BCI control with sessions where subjects achieved BCI control. There were significant differences between the two power traces only during neurofeedback for F3 θ power, specifically 0.5 s to 1 s and 2.5 s to 4 s relative to neurofeedback onset. For α frequency band at F3, there was a significant difference in power from − 4 to − 3 s and 3 to 3.5 s. For β frequency band at F3, there was a significant difference from 2.5 to 3 s (for θ, α, and β, 500 ms window, n = 60, Wilcoxon rank sum test, p < 0.0005, Bonferroni corrected). There were no significant differences between power traces for δ and γ frequency bands. Traces for all power bands, including topography plots at each time point, are provided as videos in supplementary materials (Supp. Movie 1, Supp. Movie 2, Supp. Movie 3, Supp. Movie 4, Supp. Movie 5). Qualitatively, θ modulation was specific to the frontal region during high BCI performance, indicating the specificity of neuromodulation. Notably, there was a consistent decrease in α power during the neurofeedback period, irrespective of performance, which may indicate patient engagement in mental tasks (Supp. Movie 3)42.

Figure 5

Differential θ modulation before and after BCI control. Top and Center: Patient normalized θ power during sessions with performance above optimal threshold (solid line) and first-day sessions below optimal performance threshold (dashed line). The optimal threshold was determined using ROC analysis to estimate the performance cutoff for BCI sessions with the greatest likelihood of significant θ modulation during neurofeedback. Black dashed lines indicate the onset of the start cue and neurofeedback spans. Orange shading indicates p values at durations without standard error overlap. Light orange shows duration without overlap between traces (500 ms window, n = 60, p = 0.006). Dark orange indicates durations of statistically significant differences between time traces (500 ms window, n = 60, Wilcoxon rank sum test, p < 0.0005, Bonferroni Corrected). The line shows averaged power values across each within-subject average, with the shaded region showing average standard error. Bottom: Topography of the normalized θ power during baseline and peak normalized θ during neurofeedback before and after BCI control. Theta power increases are localized to frontal channels during high BCI performance.

Full size image

BCI performance as a proxy for theta modulation

To confirm‎ that BCI performance was a proxy for θ modulation, and thus further validating the feasibility of this intervention, the distribution of BCI training sessions was plotted across all patients by performance with and without a significant increase of each frequency band at F3 during neurofeedback. Receiver operating characteristic (ROC) analysis revealed an area under the curve (AUC) of 0.73 for δ, 0.90 for θ, 0.84 for α, 0.76 for β, and 0.78 for γ (Fig. 6a, Supp. Fig. 6). BCI performance was the greatest predictor of power increase for θ frequency band at the BCI classification electrode F3. Also of note, the distributions of training session performance tended to congregate at either lower or higher performance values, with higher BCI performance corresponding to significant θ power increase during neurofeedback. This trend was not observed for any other frequency band. Sessions with significant F3 θ increase during neurofeedback also accounted for the greatest number of sessions with significant power increases across all frequency bands. There were 9 sessions with a significant δ increase, 38 with a significant θ increase, 23 with a significant α increase, 9 with a significant β increase, and 17 with a significant γ increase (Supp. Fig. 6). We found that θ modulation strength was positively correlated with BCI performance (n = 106, Spearman’s ρ = 0.67, p < 0.001; Fig. 6b).

Figure 6

The relationship between θ modulation, BCI performance, and pain relief. (a) Top: Distribution of all BCI sessions binned by performance. Blue shaded bars represent sessions where participants significantly increased their θ power during the neurofeedback period relative to the pre-intervention baseline (Rank-Biserial Spearman’s Correlation, cyan: p < 0.05 and ρ > 0). The dashed line shows the optimal threshold calculated via ROC analysis. Bottom: ROC curve of the patient performance distribution. AUC of 0.90 and threshold of 0.83 or 83% for optimally classifying sessions with significant θ modulation. (b,c) θ modulation correlates with BCI performance as well as pain relief. (b) Relationship between BCI performance and θ modulation (n = 106, Spearman’s ρ = 0.67, p < 0.001) across all subjects and sessions. (c) Relationship between pain relief (percent pre-intervention baseline) and θ modulation (n = 106, Spearman’s ρ = 0.32, p = 0.001, Bonferroni Corrected alpha = 0.005) across all subjects and sessions.

Full size image

Correlation between pain reduction and theta power change at F3

After confirmation that the BCI intervention was able to specifically serve as a proxy for θ modulation, the relationship between pain relief and θ modulation was then examined. There was a direct positive correlation between pain relief and θ modulation increase during BCI training (n = 106, Spearman’s ρ = 0.32, p = 0.001, Bonferroni Corrected alpha = 0.005, Fig. 6c). This relationship was not significant for any other power band (Supp. Fig. 6). We tested this correlation for both pain ratings acquired before and after each BCI training session (n = 106, Spearman’s ρ = 0.29 and 0.32, p = 0.0009 and 0.001, respectively, Bonferroni Corrected alpha = 0.005). We found a significant correlation between pain relief and F3 θ modulation strength in both cases, which alludes to both direct pain relief as a function of each day’s performance, as well as an overall pain reduction as a function of increasing BCI performance throughout the duration of the therapy.

Discussion

In this 6-week pilot study of patients with chronic upper extremity pain, a vibrotactile BCI therapy was used to increase frontal theta power and measure patient-reported changes in pain symptoms. Feasibility goals of patient retention and treatment protocol adherence were met. The results further indicate that this intervention reliably modulated frontal theta power and demonstrated feasibility in reducing pain. Specifically, this study met its proposed endpoint of a significant decrease in both pain interference score and pain severity score after BCI intervention in this group of patients. This result may allude to the generalizability of this treatment, as this patient group had heterogeneous demographics, initial symptoms, and trends in VAS pain ratings; however, future placebo-controlled clinical studies will be required to conclude the efficacy and limitations of this intervention. Throughout the vibrotactile BCI therapy, the magnitude of frontal θ modulation was significantly and positively correlated with pain relief as measured by VAS across BCI sessions. This was unique to frontal θ rhythms (F3 location) and was not present in any other frequency-specific power bands, indicating the specificity of θ modulation in pain relief. These results support our initial hypothesis that a non-invasive vibrotactile BCI device can feasibly target and reinforce brain activity associated with pain relief. These results warrant future work to validate this approach clinically and to better elucidate the mechanism that underlies how cortico-sensory coupling impacts chronic pain patient populations.

The experience of pain and pain relief has been shown to involve a complex neuronal network, often referred to as the pain neuromatrix. This model suggests contributions of bottom-up and top-down neural regulation of pain experience, including features beyond pain sensation such as emotional and salient components26,43. In neuromodulation studies using transcranial magnetic stimulation (TMS) or transcranial direct-current stimulation (tDCS), regions associated with top-down modulatory mechanisms of pain have been targeted, yielding successful reductions in pain symptoms44. TMS, applied to the dorsolateral prefrontal cortex in multiple experimental and chronic pain studies, led to pain reduction45,46,47,48,49,50,51. Similarly, various studies using tDCS of the prefrontal cortex have previously been shown to reduce pain. It is hypothesized that these techniques target the activity of pain perception and modulation areas, including the cingulate cortex, insula, amygdala, and thalamus44. These regions have previously been shown to play a major role in the top-down modulation of pain and are perturbed during chronic pain7,52,53,54,55. More recently, the electrophysiology of pain relief in acute and chronic settings has demonstrated that EEG frontal θ power increases are associated with the remittance of pain27,28. Source localization of these brain signals converges with previously identified anatomic sites from TMS, tDCS, and functional imaging studies27,28,44,46,47,48,49,50,51,56,57,58. These data motivated the electrode (F3) and frequency band θ configuration of this BCI system. The positive pilot results of this study support that the reinforcement of θ activity in this region may enhance the top-down capacity of critical sites such as the dorsolateral prefrontal and cingulate cortex to dampen noxious and aversive sensory perceptions.

Using a BCI to enhance the cortical physiology of pain relief as a therapeutic strategy is unique from prior approaches. Most studies using non-invasive neurofeedback for chronic pain have targeted pathological profiles of the sensory perception of pain, specifically by attempting to decrease activity associated with nociceptive processing17,18,22. A review by Patel et al. describes targeted neurofeedback systems, showing that the majority of existing work using neuromodulation for chronic pain has reinforced or inhibited EEG power at a variety of frequencies mainly over somatosensory areas (C3, Cr, Cz, T3, T4, P3, and P4 electrodes), yielding variable success in relieving symptoms of chronic pain18. The strategy for frequency band modulation was typically to increase α power (8–15 Hz), decrease β power (18–22 Hz), and occasionally decrease θ power. Only one study reinforced activity at frontal electrodes, but this was α at FP159. These studies approached the problem of chronic pain by targeting pathological activity patterns associated with the encoding of pain over somatosensory areas. While a logical approach, the pathological profiles of pain are highly variable and can differ depending on the type of chronic pain10,39,60. Our novel approach targeted reinforcement of activity patterns associated with pain relief (i.e., left-frontal θ power increase), which was shown to be present in both relief from experimentally induced tonic pain and relief from chronic pain27,28. Given the substantial similarities of pain relief cortical physiology between these two studies with very different pain experiences, this physiology may be less variable across patients experiencing pain and enable a more broadly applicable BCI-mediated neurofeedback strategy.

Beyond the choice of cortical location and frequency band, the type and character of feedback are also critical. Prior BCI systems applied to chronic pain have used visual feedback that changed with brain signals associated with somatomotor intentions. For a system to be effective in neurorehabilitation of pain circuits, visual feedback is not likely to be optimal. To optimally induce remodeling via Hebbian learning, it is important to temporally couple a desired cortical physiology with functionally relevant sensory feedback to the affected anatomic location12,13. For example, in the setting of chronic stroke, it was critical to couple brain signals associated with motor intentions with proprioceptive kinematic feedback provided by a robotic exoskeleton to achieve a functional improvement14,15,61. In the setting of chronic pain, we aimed to couple cortical signals associated with pain relief (e.g., increased frontal θ power)27,28 with pain alleviating vibrotactile sensory stimulation in the distribution of pain25,62,63, to enable neural remodeling to occur. To this end, in addition to providing visual cues of the patient’s performance in controlling their frontal θ power, patients also received sensory-relevant vibrotactile feedback with stimulus intensity proportional to the magnitude of θ modulation. Previous evidence has shown that tactile feedback provides sensory-relevant feedback and can improve BCI control64,65. Similarly, the use of vibrotactile feedback, which has previously been shown to lead to transient pain relief62,63,66, can act as a potent reinforcer to the central physiology of pain relief. From a Hebbian standpoint, the use of a BCI is especially salient because transient pain relief co-occurs precisely in time with the θ changes in the frontal lobe. While an important consideration for pain, this general configuration of BCI-mediated cortical-sensory coupling has yet to be optimized. It is unclear what are the optimal vibrotactile parameters in terms of intensity, timescale, and vibration frequency to best support central neural remodeling to lead to long-term pain relief.

While the central results of this pilot study are clinical improvements with the use of the system, there are also intriguing neurophysiology findings that support the notion of neural remodeling occurring. Specifically, baseline θ power (in the absence of task) changed over time, primarily in the left hemisphere. This change occurred both in frontal and central regions (EEG electrodes F3, C3, P3). There is prior evidence that resting-state dynamics of θ rhythms have been considered a biomarker of BCI-mediated changes in the brain67. In the setting of BCI-induced motor rehabilitation in chronic stroke, the degree of change in resting-state θ rhythms was closely correlated with the degree of motor improvement. Similarly, those patients with significant improvement in their pain demonstrated notable changes in their fronto-central baseline power. While too early to make definitive conclusions, these θ power changes merit future exploration as a potential EEG biomarker for assessing the impact of BCI-mediated pain therapy.

While the results demonstrate the potential of a novel BCI-mediated neuromodulation therapy for chronic pain, there are several limitations that merit consideration. Primarily, this clinical pilot included only six patients making more generalized claims limited. This study also lacked a control group to differentiate between the effect of θ-driven vibrotactile stimulation and placebo effects, particularly considering the setting of in-person visits to our research facility. In the future, it will be important to increase the number of patients to power the study adequately and include an additional pseudo-neurofeedback control group which will receive the visual and vibrotactile feedback, in comparable dosages, without any correlation to brain activity modulation. This will better define the role of frontal θ modulation in reducing pain, and whether coupling the vibrotactile and visual neurofeedback is effective in mitigating pain symptoms. Further, this control condition will be able to decouple the effect of frontal θ reinforcement and vibrotactile feedback over a longitudinal intervention period on observed therapeutic outcomes. A larger study will also be necessary to define electrophysiological biomarkers such as general changes in power, connectivity, and phase-amplitude coupling, which have previously been shown to be associated with pathology in chronic pain. The dosage (i.e., how long patients should be treated) and the durability of the treatment (i.e., how long the effect lasts) will need to be defined in future studies. Future iterations of this BCI will also involve the refinement of parameters, including feature discrimination and feedback optimization, to maximize therapeutic effect. As an initial pilot study, our approach to feature selection and feedback parameterization was simple. In the future, machine learning approaches can be used in several ways to increase BCI-literacy and therapeutic output68. Variance in spectral properties between patient EEG implies feature weights can be modified in order to best tailor BCIs to individual patients. Specifically, classification training can be used to select the best spectral targets for reinforcement in order to maximize therapeutic benefit. Given the importance of feedback in BCIs, machine learning methods can also be used to converge on vibrotactile stimulation feedback parameters to maximize reinforcement of frontal θ power and therapeutic outcome.

This pilot study aimed to assess the feasibility of a novel frontal θ driven BCI mediated vibrotactile therapy for chronic pain. We demonstrated the effectiveness of frontal θ reinforcement using our BCI, and that reductions in patient pain severity and pain interference scores can be achieved through use of our system.

Methods

Study participants

The study was approved by the IRB of Washington University School of Medicine, and all experimental procedures were conformed to the standards set by the latest revision of the Declaration of Helsinki. All participants signed an informed consent form prior to participation. Inclusion criteria were (1) age ≥ 18; (2) upper extremity pain located distally to the elbow (e.g. hand, wrist, forearm); (3) continuation of painful symptoms despite at least 2 previous different interventions, including topical or oral NSAID treatment, local corticosteroid injection, or surgery; (4) pain duration ≥ 4 months; (5) daily moderate to severe pain with an average intensity of ≥ 4 on 0–10 Numerical Rating Scale (NRS); (6) ability to attend training sessions for 3–5 h a week. Exclusion criteria were (1) inability or unwillingness to provide informed consent; (2) significant psychiatric or neurological condition (e.g., dementia, epilepsy, stroke) known to cause changes in EEG; (3) inability to wear a head cap due to non-healing wounds or other scalp disorders; (4) Inability to have vibrating elements put on the skin due to skin ulceration, non-healing wounds, or other skin conditions; (5) Complete lack of sensation in the testable limb; (6) Pain located proximally to the elbow; (7) The use of caffeinated products or nicotine three hours prior to each study session.

Feasibility assessment

For the scope of this pilot study, intervention feasibility was assessed on various criteria. Patient adherence and retention were reflected in participation and treatment frequency. Patients were instructed to attend BCI training sessions for a minimum of 3 days a week over the intervention duration. Patients were encouraged to participate in each session for 30–45 min but were accommodated for scheduling conflicts. During treatment sessions, patients were instructed to allocate as much attention as possible to the BCI training task (main intervention), to turn off communication devices, and to talk only during the allotted break periods. Those unable to adhere to this regimen during the intervention were given the option of restarting the treatment (i.e., restarting the protocol from the first week). The inability to adhere to the treatment schedule resulted in withdrawal from the study. Data from patients who did not complete 5 weeks of training were not analyzed. Successful patient retention was defined as completing at least 5 weeks of intervention for at least 3 days a week. Weekly session duration was assessed to determine patient adherence and retention. Intervention safety, risk, and comfort were determined by patient-reported feedback after each BCI training session. BCI patients were instructed to verbalize any discomfort they experienced and whether they felt any part of the intervention was unreasonable.

The success of this pilot’s objectives and hypotheses are predicated on the intervention’s design in successfully inducing modulation in patient θ power during BCI training. Two main steps were used to assess the feasibility of the intervention in producing and detecting modulation in left frontal θ. First, a pre-screening frontal θ modulation task was implemented. This task aimed to assess each participant’s ability to modulate θ activity and whether this modulation could be detected via analysis. Second, the qualitative changes in θ power during consecutive days of BCI training were assessed and compared to BCI performance scores and pain reduction. To this end, we eval‎uated whether BCI performance could serve as a proxy for θ modulation strength and whether the degree of θ modulation was related to changes in pain ratings.

Baseline metrics

Patients completed self-report questionnaires, underwent cognitive testing, baseline brain activity recordings, and screening for task-specific θ modulation prior to any BCI-intervention. Specifically, participants rated their pain intensity on VAS, completed BPI, NPSI, and the PCS. In addition, HADS was used to determine patient anxiety and depressive symptoms. QST was performed to assess patients’ somatosensory profile69,70,71. The EEG recordings took place in a quiet study room with minimal outside disturbance. Participants were comfortably seated in a stationary chair and were asked to limit their movements during brain activity recordings. The details of the methods used to obtain these metrics are described in the Supplemental Methods.

System overview

The BCI system consisted of a DSI-24 EEG recording system (Wearable Sensing, San Diego, CA, USA) connected to a Windows PC Laptop via Bluetooth, enabling easy and comfortable measurement of high-fidelity EEG signals in a laboratory environment, with minimal setup time. The HVSA provided real-time vibrotactile neurofeedback, and a monitor provided visual feedback in the form of a vertical cursor task. BCI2000, a general-purpose, open-source modular software, was used to manage data acquisition, real-time signal processing, application execution, event logging, and system synchronization72. This device consisted of a custom-designed adjustable frame able to fit most hand sizes, with two (upper and lower) grids of 24 small vibrating motors contacting both sides of the hand. This device executed vibration commands via an Arduino Microcontroller connected to the PC.

Real-time signal processing

A real-time spatial and spectral filter was applied to acquired EEG data to extract frontal-midline θ power. A 1 Hz high pass and 50 Hz low pass filter were applied in real-time to avoid the influence of artifacts from line noise or channel drift on BCI classification accuracy. Signals were common average re-referenced, then spatially filtered such that F3, Fz, and F4 electrodes were selected as feature input channels. θ power was calculated in real-time by an autoregressive (AR) spectral estimation using a 500 ms sliding window. The AR estimation was selected as it provides a computationally efficient means for real-time spectral estimation and has been shown to have high spectral resolution even for shorter data segments73. 4–6 Hz power was defined as the as the driving spectral feature. 4–6 Hz was chosen over 4–7 Hz to minimize the effect of α band variance on the spectral filter. Previous evidence has suggested task-specific and resting-state variance of α power both within and between populations of subjects74. Further, the proximity of commonly defined α frequencies to our desired spectral feature of θ could have potentially influenced the specificity of our spectral estimations due to spectral leakage75. To avoid these confounds, we chose not to include the frequencies closest to either θ or α as part of our real-time feature extraction. Further elaboration on the selection of 4–6 Hz can be found in Supplemental Results. This real-time processing method enabled computationally efficient processing in order to minimize the latency between real-time processing of electrophysiology and execution of neurofeedback commands.

Classifier output calibration and normalization

To define baseline θ power for real-time classification during BCI training for each patient, a zero-mean distribution of the subject’s resting θ power at the F3 channel was generated. This was done during a set of ten “adaptation” BCI trials, where the BCI2000 adaptation function was enabled. This function optimized normalization values for calculating relative changes in θ power by determining a zero-mean value for θ power during rest, and an arbitrary gain value which aided in magnifying output commands. This gain value was also optimized such that patient eyeblink artifacts, which lead to transient increases in δ-θ power, did not significantly influence BCI performance metrics (primarily by regulating the rate of trial completion, such that consecutive eyeblinks alone could not lead to trial success). After ten of these adaptation trials, a set zero mean value was uniquely used for the remainder of the BCI therapy for each patient. The output of the classifier was thus a θ power value normalized to the baseline zero-mean value, multiplied by the determined gain value:

The output value was used to drive neurofeedback, with positive values indicating increases in θ power at F3 and negative values indicating decreases. The overall magnitude of the output value was reflected in the velocity of the visual cursor and the intensity of the vibrotactile feedback, as described below.

BCI training intervention

Patients underwent at least three one-hour BCI training sessions per week for at least 5 weeks, but no more than 6 weeks (intervention duration). This duration and weekly session number were selected based on intervention durations of previous longitudinal BCI or neurofeedback therapies18,76. VAS pain scores were obtained from patients upon arrival to the study room and immediately after each day’s BCI training session. Baseline recordings were obtained via EEG during the pre-BCI training period each day (5 min), during the BCI training session (up to 45 min), and during each day’s post-BCI training session (5 min). The study design is shown in Fig. 2a. During baseline recordings (pre and post BCI training), patients were asked not to perform specific mental tasks and to remain calm and relaxed while remaining awake. A fixation cross was presented on the display, and patients were asked to fix their gaze on the cross. Patients were instructed to blink at a natural rate, ensuring that they did not blink inconsistently between parts of the BCI trials. After the pre-BCI training recording, patients were asked to place their affected hand into the HVSA.

BCI training consisted of three 15-min blocks. Each trial consisted of a ten-second baseline period, a three-second start cue period, a 20-s limit neurofeedback period, a three-second stop cue, and a seven-second post stop-cue pre-baseline buffer period (Fig. 2b). The subjects were presented with a yellow on-screen cursor which remained stationary at the center of the screen during non-neurofeedback periods. At the top of the screen, a green rectangle indicated the target and remained on top of the screen during all periods. The green rectangle served as the target and would always appear at the top of the screen in the same location. During the start cue, the word “Start” was presented in the center of the cursor, and a random prime number between 3 and 19 appeared in the center of the top rectangle in cyan. This prime number remained in the rectangle for the duration of the neurofeedback period, as an option for patients to use serial subtraction as the driving task for frontal θ modulation if needed. Additional methods regarding the role of the prime numbers are provided in Supplemental Methods. After three seconds, the word “Start” would disappear, and the neurofeedback period would begin. During the neurofeedback period, the center cursor would either move up if the subject increased their left frontal-midline θ power above baseline or move down if their left frontal-midline θ power decreased below baseline average. The speed of cursor movement corresponded to the magnitude of increase or decrease (see description of intervention). During neurofeedback, the HVSA would stimulate the participant’s affected area when there was an increase in θ power. The intensity of the vibration corresponded to the magnitude of θ power increase. The vibration was pulsed at a random frequency selected from 4, 5, 6, or 7 Hz. If the participant was able to reach the target by moving the cursor to the top rectangle within the allocated 20-s period, the stop cue would initiate, and a burst of vibration at maximum intensity was delivered for the duration of the three-second stop period. If the participant could not reach the target within the allocated time limit, the stop cue would initiate with no vibration delivered. The stop cue consisted of the cursor returning to the center of the screen, and the word “Stop” was displayed in the center of the cursor. The seven-second buffer period allowed for any brain activity changes resulting from the stop period and the neurofeedback period to wash out such that the activity recorded during the baseline period was unaffected. Figure 2b shows the design of a single trial within a training session. Additional information describing the methodology and characteristics of the continuous BCI used in this study is shown in Supp. Table 5. At the end of each week on Friday, patients received an email containing a survey. The survey was used to obtain self-reported BPI and NPSI metrics each week outside of the experimental environment.

Data collection and outcome measures

Electrophysiological activity recordings throughout the therapy duration were carried out using dry-surface EEG. Primary and secondary outcome measures were self-reported by patients.

EEG recording

EEG was recorded using a 24 wireless dry electrode headset in an international 10–20 system (DSI 24, Wearable Sensing, San Diego, CA, USA). EEG was referenced to the Pz electrode and sampled at 300 Hz with a ground electrode placed on the earlobe. Electrode impedance was maintained below 10 kΩ.

Primary and secondary outcome measures

The primary outcome measure was the BPI PIS sub-score. Secondary outcome measures were the NPSI total score (50 points total, consisting of sub-scores of 10 points each for burning pain, pressing pain, paroxysmal pain, evoked pain, and paraesthesia/dysaesthesia, where higher scores reflect more severe symptoms of a specific quality) and BPI PSS sub-score, which were recorded at the end of each week for each patient77,78. These outcomes, especially BPI, are considered robust and reliable outcome measures for trials in patients with chronic pain, where higher scores indicate worse or more severe symptoms related to the quality and intensity of pain (PSS, 0–10 with 0 = no pain and 10 = worst imaginable pain), and the degree that the pain symptoms interfere with everyday life (PIS, 0–10 with 0 = no interference and 10 = interferes completely)79,80. These metrics have previously been shown to have high test–retest reliability and have been adapted to a variety of diseases and types of pain77,81,82,83. For reference, previous work has reported that 1-point changes in PIS reflect a typical 0.5 standard deviation change, and is a benchmark for minimal clinical importance, and improvements of 2 points or greater reflect moderate or meaningful changes. For pain severity, recorded as PSS, 10% point changes were considered minimally important, while reductions of 30% or more were considered to be clinically important83. These metrics were acquired via patient response on a survey emailed to each patient at the end of each week. Survey outreach was handled by the REDCap database/survey management software. VAS scores (additional secondary outcome, where 0 = ‘no pain’ and 100 = ‘worst pain imaginable’) were acquired from patients before and after each BCI session and used for secondary analysis to relate changes in pain ratings, BCI performance, and electrophysiological features84.

Offline data analysis

Processing of electrophysiological data and statistical analysis was conducted in MATLAB environment (Mathworks, Nattick, MA, USA).

EEG processing

EEG data were preprocessed by first applying a 1 Hz high pass and 40 Hz low pass Butterworth filter to the signal. Next, high amplitude signal artifacts were removed by thresholding above each subject’s manually determined eyeblink artifact threshold. This was done for each data file qualitatively by manually identifying the largest magnitude eyeblink artifact and then setting the threshold above this value. The threshold was set above the magnitude of the eyeblink artifact because independent component analysis (ICA) was used to remove eyeblink artifacts, which minimized the amount of electrophysiological activity removed from the data. Signals containing samples above this threshold were deleted and then smoothed using a five-point moving average interpolation. The signal was then common average re-referenced, and ICA was performed to remove eyeblinks, eye saccades, lateral eye movements, muscle-related, and cardiac-related artifacts from the data85.

Band-limited amplitude time series was extracted by applying a Butterworth with frequency cutoffs as follows: 1–3 Hz for delta (δ), 4–7 Hz for θ, 7–13 Hz for alpha (α), 13–30 Hz for beta (β), and 30–40 Hz for gamma (γ). Power envelopes for each signal were calculated by squaring the absolute value of the Hilbert transform of the time series. The mean power envelope within each relevant timeframe was used for overall power calculations.

For electrophysiological comparisons within and between subjects across multiple days, power values were normalized to baseline recordings by z-scoring all power envelope time series to power envelope from EEG data recorded before any BCI intervention.

For analysis of continuous, relative power changes within single BCI sessions, a time–frequency analysis was conducted by time locking pre-processed signals acquired from each session across each subject to trial start and neurofeedback onset. Time–frequency was determined by the mean Fast Fourier Transform of all trials with Hamming window tapering85. Time–frequency signals were binned into power bands as described above, then normalized to the minimum and maximum power values within each session, exemplifying the spatial–temporal modulation of different frequency-specific activities within each BCI training session.

For analysis of spectral properties between conditions, specifically, to visualize qualitative differences in power spectral density plots, Welch’s power spectral density estimate was used with a Hamming window equal to ten percent of the epoch duration86. This analysis was applied to demonstrate qualitative differences in θ power modulation during the pre-screening mental arithmetic task and BCI training.

Statistical analysis

Statistical tests were conducted using non-parametric methods. None of the implemented statistics were used for concluding clinical or mechanistic outcomes. Rather, the tests employed for this study aimed at quantitatively supporting the feasibility and potential of the device for larger, placebo-controlled clinical trials and for identifying potentially relevant mechanisms evoked by BCI intervention that justify further investigation. For correlation analysis conducted between pain intensity and treatment time, BCI performance (measured as BCI accuracy: # trials goals reached/# total trials) and treatment time, BCI bit rate and treatment time, θ power modulation and BCI performance, and θ power modulation and pain relief, Spearman’s correlation was used. Statistical power analysis showed that 87 sessions were ample to detect a significant correlation strength of 0.3 (Moderate Correlation) with a power of 95%. The present study includes data from 106 BCI sessions. BCI throughput (bit rate) was calculated using Fitts’s metrics in bits per second as an alternative means of measuring BCI performance in addition to accuracy87. Throughput (bits per second) was calculated as the cursor task’s index of difficulty which was equal to the log2(1 + 238) over the time to reach the target for each trial. 238 was the total distance the cursor traveled in pixels to reach the goal. For trials where the goal was not reached, the time to reach the target used for the calculation was 20 s, as this was the maximum duration participants had to reach the goal. Patients with a significant increase in BCI performance accuracy over time were classified as performers. Patients with significant decreases in pain over time were classified as VAS pain responders. To determine the significance and directionality of power modulation from baseline, we used a non-parametric rank-biserial Spearman’s correlation88. This method enabled correlation between continuous ranked variables (i.e., EEG data from two conditions) and a dichotomy (i.e., a generated data vector of − 1 s and 1 s corresponding to the lengths of either EEG data conditions). We chose this method over a Wilcoxon rank-sum test as it reflected the directionality and magnitude of significantly different frequency-specific power. To compare group-level changes in the outcome measures, a two-sided Wilcoxon signed-rank test was used. For statistical comparisons of power differences at single channels between conditions, a two-sided Wilcoxon rank-sum test was used. To account for multiple comparisons for these statistical tests, a Benjamini and Hochberg procedure was used to calculate the false discovery rate (q-value), as this test has shown to be more conservative and better fit for fewer comparisons89. Due to the small study size, the median was used as the central tendency value for comparisons across all subjects, and the median absolute deviation (MAD) was used for variance. The mean and standard error were used for comparisons of electrophysiological data consisting of multiple trials. For detecting whether there were significant differences between weekly intervention participation time, a Kruskal–Wallis test was used90. Repeated measurements were taken across all participants for both clinical and electrophysiological variables.

ROC analysis was performed to assess whether BCI performance was an accurate proxy for θ modulation. For each session, we eval‎uated whether significant θ modulation was achieved by performing a rank-biserial Spearman’s correlation between pre-BCI intervention baseline θ power and each session’s neurofeedback period θ power at F3. Baseline epochs and neurofeedback epochs were the same lengths. We defined sessions with successful θ power increases as those with a p < 0.05 and Spearman’s ρ > 0. We then performed ROC analysis on the distribution of all sessions, including sessions with or without significant increases in θ power at F3. An optimal performance threshold for identifying a BCI performance value most likely to reflect θ power increase at F3 was determined by maximizing the true positive probability (i.e., session performance truly reflects significant θ power increase) and minimizing false positive probability (i.e., session performance reflects significant θ increase without actual increase). Sessions with high BCI performance were defined as sessions with performance above the optimal threshold calculated with ROC analysis (83%). Sessions with no BCI control (naïve sessions) were defined as each patient’s first BCI training session with performance below the optimal threshold.

For comparisons of normalized band-power differences between traces generated by these groups, we compared data within 500 ms bins using a two-sided Wilcoxon rank-sum test. Permutation tests (1000 iterations) were performed for all correlations, a biserial correlation analysis. We next compared the relationship between magnitude θ increase (Spearman’s ρ) with BCI performance and pain relief using Spearman’s correlation. Significance levels were Bonferroni corrected for these tests.

Data availability

The authors would like to indicate that our data is available upon reasonable request by contacting the corresponding author or the first author, Phillip Demarest (d.phillip@wustl.edu).

References

1. Hylands-White, N., Duarte, R. V. & Raphael, J. H. An overview of treatment approaches for chronic pain management. Rheumatol. Int. 37, 29–42 (2017).

   Article CAS PubMed Google Scholar 

2. Finnerup, N. B. et al. Neuropathic pain clinical trials: Factors associated with decreases in estimated drug efficacy. Pain 159, 2339–2346 (2018).

   Article CAS PubMed PubMed Central Google Scholar 

3. Vowles, K. E. et al. Rates of opioid misuse, abuse, and addiction in chronic pain: A systematic review and data synthesis. Pain 156, 569–576 (2015).

   Article PubMed Google Scholar 

4. Ballantyne, J. C. Opioids for the treatment of chronic pain: Mistakes made, lessons learned, and future directions. Anesth. Analg. 125, 1769–1778 (2017).

   Article PubMed Google Scholar 

5. Frizon, L. A. et al. Deep brain stimulation for pain in the modern era: A systematic review. Neurosurgery 86, 191–202 (2020).

   Article PubMed Google Scholar 

6. D’Souza, R. S. et al. Adverse events associated with 10-kHz dorsal column spinal cord stimulation: A 5-year analysis of the manufacturer and user facility device experience (MAUDE) database. Clin. J. Pain 38, 320–327 (2022).

   Article PubMed Google Scholar 

7. Bliss, T. V. P., Collingridge, G. L., Kaang, B.-K. & Zhuo, M. Synaptic plasticity in the anterior cingulate cortex in acute and chronic pain. Nat. Rev. Neurosci. 17, 485–496 (2016).

   Article CAS PubMed Google Scholar 

8. Čeko, M. et al. Partial recovery of abnormal insula and dorsolateral prefrontal connectivity to cognitive networks in chronic low back pain after treatment. Hum. Brain Mapp. 36, 2075–2092 (2015).

   Article PubMed PubMed Central Google Scholar