A stimulation-orthogonal neural biomarker of therapeutic response in deep brain stimulation

Technische Samenvatting: Een stimulatie-orthogonale neurale biomarker voor therapeutische respons in diepe hersenstimulatie

Probleemstelling
Closed-loop neuromodulatie, met name bij diepe hersenstimulatie (DBS) voor de ziekte van Parkinson, kampt met een fundamenteel "verstrengelingsprobleem". In huidige systemen wordt de neurale readout die voor feedback wordt gebruikt, opgenomen vanuit hetzelfde circuit dat door de stimulatie wordt verstoord. Bijgevolg zijn neurale kenmerken die veranderen tijdens stimulatie mogelijk een reflectie van de loutere aanwezigheid van de interventie (circuit-engagement) in plaats van echte therapeutische verbetering. Standaard biomarkers, zoals bèta-bandvermogen in de subthalamische nucleus (STN), zijn zeer gevoelig voor de aanwezigheid van stimulatie, maar slagen er niet in om therapeutische van niet-therapeutische staten te onderscheiden. Zo kan laagfrequente stimulatie het circuit weliswaar activeren zonder klinisch voordeel te bieden, terwijl standaard decoders dit engagement nog kunnen interpreteren als een therapeutisch signaal. De kernuitdaging bestaat eruit om de dimensies van de neurale representatie die de aanwezigheid van stimulatie identificeren, te scheiden van de dimensies die het therapeutische effect van de stimulatie meten, binnen een enkele representatie, zonder daarbij te vertrouwen op gesuperviseerde klinische uitkomstlabels tijdens het scheidingsproces.

Methodologie
De auteurs stellen een oplossing voor met behulp van Generative Causal Explanations (GCE), een deep learning-framework dat hoogdimensionale neurale kenmerken partitioneert in een gesuperviseerde latente component en orthogonale hulpcomponenten.

• Experimenteel ontwerp: De studie maakte gebruik van subthalamische lokale veldpotentialen (LFP's) van negen patiënten met gevorderde Parkinson (ontwikkelingscohort) en één prospectief uitgesloten patiënt (P11). Opnames werden verkregen bij vier stimulatiefrequenties (125, 110, 85, 55 Hz) en twee amplitudes (klinisch geprogrammeerde volledige amplitude en halve amplitude).
• Trainingsstrategie (De "Quarantaine"-aanpak): In plaats van de gesuperviseerde dimensie te trainen om de klinische uitkomst te voorspellen (wat geconfundeerd wordt door de aanwezigheid van stimulatie), trainden de auteurs het GCE-model om de aanwezigheid van stimulatie (ON vs. OFF) te detecteren met behulp van uitsluitend half-amplitude opnames. Bij een halve amplitude wordt het circuit weliswaar betrouwbaar geactiveerd, maar zijn de klinische voordelen minimaal en inconsistent.
  • Doel: Het classifier-gealigneerde latente dimensie (Dim 0) dwingen om het "stimulatiesignatuur" te absorberen.
  • Restrictie: De resterende hulpdimensies (Dims 1–5) worden beperkt om de neurale structuur te behouden die orthogonaal staat aan de taak van stimulatie-detectie.
• Kenmerkenruimte: De input bestond uit 70 kenmerken per segment van 10 seconden, inclus': spectrale vermogens (alfa-, bèta-, gamma-banden), burst-dynamiek, fase-amplitude koppeling, interhemisferische partiële gerichte coherentie (PDC) en niet-lineaire dynamische maten.
• Evaluatie: De getrainde encoder werd bevroren en toegepast zonder hertraining op:
  1. Volledige-amplitude opnames over alle frequenties in het ontwikkelingscohort.
  2. De prospectief uitgesloten patiënt (P11).
  3. Een onafhankelijke cohort van drie patiënten die via video-kinematica werden beoordeeld (ander protocol, andere uitkomstmodaliteit).
• Validatie: De stimulatie-gealigneerde dimensie werd getest op haar vermogen om ON/OFF-statussen te classificeren (generaliseerbaarheid) en op haar gebrek aan correlatie met motorische verbetering (negatieve controle). De hulp-subruimte werd getest op haar vermogen om gegradeerde motorische verbetering te volgen en tussen effectieve en ineffectieve stimulatie-instellingen te onderscheiden.

Belangrijkste Resultaten

1. Succesvolle ontstrengeling: De stimulatie-gealigneerde dimensie (Dim 0) generaliseerde succesvol als een marker voor de aanwezigheid van stimulatie over verschillende patiënten en amplitudes (Leave-One-Subject-Out AUC = 0,919 bij halve amplitude; 0,942 bij volledige amplitude). Cruciaal was dat Dim 0 geen associatie vertoonde met motorische verbetering (totale score, rigiditeit, bradykinesie of tremor), wat bevestigt dat het enkel de interventiesignatuur vastlegde.
2. Therapeutisch signaal in de orthogonale ruimte: De hulp-subruimte (met name de gecombineerde gesigneerde afstand van Dims 1 en 2) volgde de gegradeerde motorische verbetering bij volledige amplitude (Spearman $\rho = 0,50$, $p=0,002$).
   • Symptoomspecificiteit: Het signaal was het sterkst voor rigiditeit ($\rho = 0,479$) en vertoonde een zwakkere, readout-afhankelijke associatie met bradykinesie. Het vertoonde geen robuuste relatie met tremor.
   • Frequentieafhankelijkheid: De hulp-metriek identificeerde laagfrequente (55 Hz) stimulatie als een grote "anti-therapeutische" verplaatsing. Ondanks het feit dat 55 Hz het circuit sterk activeert (grote Euclidische verplaatsing in de hulp-ruimte), bewoog het alle patiënten naar de minst therapeutische regio van de latente ruimte, een onderscheid dat door standaard bèta- of gammapower niet werd opgepikt.
3. Generalisatie:
   • Prospectieve patiënt: In de uitgesloten patiënt (P11) correleerde de hulp-metriek significant met de verbetering van rigiditeit ($\rho = 0,731$, $p=0,040$) zonder hertraining.
   • Onafhankelijke kinematische cohort: Het model generaliseerde naar een aparte cohort die werd beoordeeld via markerloze video-kinematica. De hulp-metriek correleerde met de consistentie van vingerklopjes en de snelheid van handbewegingen, wat de therapeutische as valideert over verschillende uitkomstmodaliteiten heen.
4. Kenmerktoeschrijving: De stimulatie-gealigneerde as werd gedreven door interhemisferische PDC en hoge bèta-power. In tegenstelling hiertoe werden de therapeutische hulp-dimensies gedreven door lage-bèta burst-dynamiek, fase-amplitude koppeling en alfa-band vermogen, wat hen linkt aan de parkinsoniaanse motorische pathofysiologie in plaats van aan stimulatie-artefacten.
5. Superieuriteit ten opzichte van baselines: Onder hetzelfde lineaire mixed-effects framework faalden ongecontroleerde Principal Component Analysis (PCA) en standaard spectrale band-vermogensmaten in het terugwinnen van de therapeutische structuur of het onderscheiden van de anti-therapeutische 55 Hz conditie.

Betekenis en Claims
Het artikel beweert een nieuw architecturaal sjabloon te bieden voor het ontdekken van biomarkers in closed-loop settings waar de opname- en stimulatieplaatsen samenvallen. Door de stimulatiesignatuur te "quarantineren" via gesuperviseerd leren op een conditie met weinig voordeel, legt de methode een therapeutische structuur bloot in de orthogonale subruimte zonder dat daar tijdens de scheidingsfase uitkomstlabels voor nodig zijn.

• Ontkoppeling van Engagement en Voordeel: De aanpak laat zien dat circuit-engagement (gedetecteerd door Dim 0) en therapeutisch voordeel (gecodeerd in de hulp-ruimte) verschillende neurale fenomenen zijn die wiskundig gescheiden kunnen worden.
• Voorbij Bèta-power: De resultaten suggereren dat therapeutische effectiviteit niet uitsluitend wordt gedefinieerd door bèta-suppressie. De hulp-metriek ving nuances op, zoals de anti-therapeutische aard van laagfrequente stimulatie, die standaard spectrale samenvattingen misten.
• Generaliseerbaarheid: De bevroren representatie generaliseerde over verschillende stimulatie-amplitudes, frequenties, patiënten en zelfs naar een andere uitkomstmodaliteit (kinematica versus klinische beoordeling), wat suggereert dat de geleerde geometrie fundamentele therapeutische principes vastlegt in plaats van dataset-specifieke artefacten.
• Toekomstige toepassing: De auteurs stellen dat deze hulp-as als een continu feedbacksignaal kan dienen voor adaptieve DBS, waardoor controllers de volledige stimulatieparameterruimte (frequentie en amplitude) kunnen optimaliseren in plaats van alleen amplitude op basis van bèta-suppressie. Dit kan systemen voorkomen die anti-therapeutische instellingen versterken die slechts het circuit activeren.

De auteurs blijven bescheiden over de reikwijdte en merken op dat de bevindingen gebaseerd zijn op single-center data met kleine validatiecohorten en dat de longitudinale stabiliteit nog niet is getest. Ze stellen echter dat de methode een robuust bewijs van concept biedt voor het scheiden van de aanwezigheid van een interventie van het effect van die interventie in neurale opnames.