Ten protokół opisuje metodę umieszczania elektrod sterowanych neuronami do ogniskowej, przezczaszkowej stymulacji prądem stałym (tDCS) podawanej podczas funkcjonalnego rezonansu magnetycznego (fMRI).
Artykuł metodologiczny
* These authors contributed equally
Ten protokół opisuje metodę umieszczania elektrod sterowanych neuronami do ogniskowej, przezczaszkowej stymulacji prądem stałym (tDCS) podawanej podczas funkcjonalnego rezonansu magnetycznego (fMRI).
Przezczaszkowa stymulacja prądem stałym (tDCS) to nieinwazyjna technika stymulacji mózgu, która pozwala na modulację pobudliwości i plastyczności ludzkiego mózgu. Zogniskowane konfiguracje tDCS wykorzystują określone układy elektrod, aby ograniczyć przepływ prądu do opisanych obszarów mózgu. Jednak skuteczność zogniskowanego tDCS może być zagrożona przez błędy w ustawieniu elektrod na skórze głowy, co skutkuje znacznym zmniejszeniem dawki prądu docierającej do docelowych obszarów mózgu dla tDCS. Rozmieszczenie elektrod pod kierunkiem neuronawigacji w oparciu o anatomię głowy i mózgu danej osoby uzyskaną z danych obrazowania strukturalnego rezonansu magnetycznego (MRI) może być odpowiednie do poprawy dokładności pozycjonowania.
Ten protokół opisuje metodę umieszczania elektrod sterowanych neuronami dla zogniskowanej konfiguracji tDCS, która jest odpowiednia do jednoczesnego podawania podczas funkcjonalnego MRI (fMRI). Określamy również ilościowo dokładność umieszczenia elektrod i badamy dryf elektrod w równoległym eksperymencie tDCS-fMRI. Krytyczne kroki obejmują optymalizację pozycji elektrod w oparciu o obecne modelowanie, które uwzględnia anatomię głowy i mózgu danej osoby, wdrożenie neuronażowanego umieszczania elektrod na skórze głowy oraz podawanie zoptymalizowanego i ogniskowego tDCS podczas fMRI.
Regionalna precyzja umieszczenia elektrod jest określana ilościowo za pomocą normy euklidesowej (norma L2) w celu określenia odchyleń rzeczywistych od zamierzonych pozycji elektrod podczas równoległego badania tDCS-fMRI. Wszelkie potencjalne przemieszczenia elektrod (dryft) podczas eksperymentu są badane poprzez porównanie rzeczywistych pozycji elektrod przed i po akwizycji fMRI. Ponadto bezpośrednio porównujemy dokładność umieszczania neuronażowanych tDCS z dokładnością osiąganą dzięki podejściu do celowania opartemu na skórze głowy (system elektroencefalografii 10-20 (EEG)). Analizy te wykazują lepszą dokładność umieszczania elektrod w neuronawigacji w porównaniu z umieszczaniem elektrod na skórze głowy i znikomym dryfem elektrody w ciągu 20-minutowego okresu skanowania.
Przezczaszkowa stymulacja prądem stałym (tDCS) to nieinwazyjna technika stymulacji mózgu, która pozwala na modyfikację funkcji poznawczych i fizjologicznych mózgu w kontekstach eksperymentalnych i klinicznych1,2,3. Ostre podanie tDCS może powodować przejściowe zmiany w pobudliwości neuronów, a skutki uboczne utrzymują się od kilku minut do godzin po stymulacji4,5. Przyłożony prąd nie indukuje potencjałów czynnościowych, ale raczej przejściowo przesuwa spoczynkowy potencjał błonowy neuronu w kierunku depolaryzacji lub hiperpolaryzacji, co skutkuje zwiększoną lub zmniejszoną pobudliwością neuronów na poziomie makroskopowym przy użyciu standardowych protokołów4,5,6. Ponadto, jeśli chodzi o efekty plastyczności synaptycznej tDCS, badania na zwierzętach i ludziach wykazały, że tDCS indukuje długotrwałe procesy podobne do wzmocnienia i depresji (LTP i LTD)4,5.
W układzie motorycznym, modulacja motorycznych potencjałów wywołanych (MEP) pozwala na bezpośrednią ocenę neurofizjologicznych efektów tDCS na lokalną pobudliwość korową7. Jednak to podejście nie może określić ilościowo neuronalnego wpływu tDCS na funkcje poznawcze wyższego rzędu wspierane przez funkcjonalne sieci mózgowe na dużą skalę8. Wpływ na sieci mózgowe można zbadać, łącząc tDCS z nowoczesnymi technikami obrazowania funkcjonalnego9,10. Wśród nich funkcjonalne obrazowanie metodą rezonansu magnetycznego (fMRI) stało się najczęściej stosowanym podejściem, ponieważ zapewnia doskonałą rozdzielczość przestrzenną i wystarczającą czasowo, aby ujawnić mechanizmy neuronalne, za pomocą których tDCS wpływa na lokalną aktywność mózgu w miejscu stymulacji i wielkoskalowych sieci neuronowych11,12,13,14.
Do tej pory, połączone badania fMRI-tDCS głównie wykorzystywały tak zwane konwencjonalne konfiguracje tDCS, które wykorzystują stosunkowo duże gumowe elektrody o wymiarach od 25 do 35 cm2 (5 x 5 cm2 i 5 x 7 cm2) umieszczone w kieszeniach z gąbki nasączonej solą fizjologiczną15,16. Konfiguracje te rzutują prąd między dwiema elektrodami, które są zwykle przymocowane do (a) docelowego obszaru mózgu dla tDCS i (b) elektrody powrotnej nad niedocelowymi obszarami mózgu lub obszarami pozaczaszkowymi (np. Barkiem). Powoduje to rozległy przepływ prądu przez mózg, wpływając na regiony inne niż region docelowy, komplikując w ten sposób założenia przyczynowe i interpretacje dotyczące neuronalnego pochodzenia efektów tDCS17.
Bardziej precyzyjne celowanie przestrzenne można osiągnąć za pomocą ogniskowego tDCS18. Konfiguracje te wykorzystują układy mniejszych elektrod rozmieszczonych w bliskiej odległości od siebie lub za pomocą katody w kształcie pierścienia umieszczonej wokół centralnej anody w celu ograniczenia przepływu prądu do obszaru docelowego18,19. Symulacje komputerowe przepływu prądu elektrycznego sugerują, że zogniskowany tDCS może skutkować wyższą precyzją przestrzenną przepływu prądu do obszaru docelowego niż konwencjonalne montaże20. Co więcej, badania behawioralne wykazały regionalną i specyficzną dla zadania modulację behawioralną przy użyciu konfiguracji ogniskowych19,21,22. Jednak tylko w kilku badaniach wykorzystano zogniskowany tDCS podczas fMRI. Badania te były w stanie ustalić wykonalność tego podejścia i dostarczyły pierwszych dowodów na modulację neuronalną specyficzną dla regionu19,23.
Jednakże, ze względu na regionalnie precyzyjne dostarczanie prądu, konfiguracje tDCS mogą być bardziej wrażliwe na błędy pozycjonowania elektrod na skórze głowy niż konwencjonalne montaże. Na przykład Seo i wsp. wykazali, że błędy pozycjonowania wynoszące 5 mm w ogniskowej korze ruchowej zmniejszyły szczytową polaryzację somatyczną w gałce ręki nawet o 87%24. Co więcej, niedawne badanie modelowania obliczeniowego wykazało, że przemieszczenie elektrody z pozycji zamierzonych dla ogniskowej w porównaniu z konwencjonalnymi konfiguracjami spowodowało znaczne zmniejszenie dawki prądu w obszarach docelowych dla tDCS, w zakresie od 26% do 43%25. W związku z tym stwierdzono, że przyszłe badania powinny rutynowo obejmować odpowiednie metody poprawy pozycjonowania elektrod i weryfikacji pozycjonowania elektrod przed i po fMRI5.
W niniejszym badaniu opisujemy metodę pozycjonowania elektrod sterowanych neuronami dla nowatorskiego, kompatybilnego z fMRI ogniskowego ustawienia 3 x 1 tDCS (tj. trzy pojedyncze katody, które są ułożone w okręgu wokół pojedynczej centralnej anody), która jest obecnie wykorzystywana we wspólnym konsorcjum badawczym finansowanym przez Niemiecką Fundację Naukową (DFG Research Unit 5429, https://www.memoslap.de). Konsorcjum bada behawioralny i neuronalny wpływ zogniskowanego tDCS na uczenie się i pamięć oraz predyktory reakcji na stymulację w czterech domenach funkcjonalnych (tj. funkcjach wzrokowo-przestrzennych, językowych, motorycznych i wykonawczych). Strukturalne dane MRI T1 i T2 uczestników badania są pozyskiwane podczas skanowania wyjściowego. Dane te są wykorzystywane do zindywidualizowanych symulacji przepływu prądu26 w celu określenia pozycji elektrod na skórze głowy, które maksymalizują przepływ prądu do obszaru docelowego u poszczególnych uczestników badania. Jako przykład, protokół ten będzie opisywał kierowane neuronami celowanie w indywidualnie określone pozycje elektrod wyśrodkowane nad prawą grzbietowo-boczną korą przedczołową (rDLPFC) u jednego uczestnika.
Sekcja reprezentatywnych wyników jest oparta na danych obrazowania strukturalnego uzyskanych przed i po równoczesnym tDCS-fMRI w trzech podprojektach jednostki badawczej. Badania te dotyczyły prawej kory potyliczno-skroniowej (rOTC), lewej kory skroniowo-ciemieniowej (lTPC) i rDLPFC. Dane uzyskano w Zakładzie Neurologii Uniwersytetu Medycznego w Greifswaldzie. Korzystając z tych danych, chcieliśmy osiągnąć dwa główne cele: (1) określić ilościowo precyzję przestrzenną rozmieszczenia elektrod sterowanych przez neurony poprzez porównanie "zamierzonych" i empirycznie określonych "rzeczywistych" pozycji elektrod 25, oraz (2) zbadanie stopnia przemieszczenia elektrody w trakcie sesji fMRI (tj. dryf elektrody). Czynniki te mają kluczowe znaczenie dla poprawy dokładności i wiarygodności efektów tDCS w równoległych badaniach tDCS-fMRI27. Ponadto dokładność celowania tDCS sterowanego neuronami jest porównywana z podejściem opartym na skórze głowy przy użyciu danych z poprzedniego badania tDCS-fMRI naszej grupy25.
Wszystkie procedury eksperymentalne przedstawione w tym protokole zostały sprawdzone i zatwierdzone przez komisję etyczną Uniwersytetu Medycznego w Greifswaldzie. Wszyscy uczestnicy wyrazili świadomą zgodę przed włączeniem do badania i wyrazili zgodę na anonimową publikację swoich danych.
1. Badanie przeciwwskazań i uwagi ogólne
2. Podstawowy skan MRI i zindywidualizowane modelowanie prądu
3. Neuronawigacja
4. tDCS-fMRI
Uwzględniono dane od 43 zdrowych młodych uczestników (20 mężczyzn/23 kobiety, wiek 24,74 ± 5,50 lat). Uczestnicy ukończyli do czterech sesji fMRI. Przed każdą sesją fMRI przeprowadzono neuronavigacyjne rozmieszczenie elektrod. W analizach danych ujęto łącznie 338 zestawów danych reprezentujących pozycje centralnych anod przed i po badaniu fMRI.
Aby wyznaczyć zamierzone pozycje elektrod, przeprowadzono indywidualne modelowanie prądu z wykorzystaniem danych strukturalnego MRI pozyskanych podczas oddzielnej sesji MRI (przed sesjami fMRI). Modelowanie prądu miało na celu optymalizację pozycji elektrod na skórze głowy dla układów 3 x 1, aby zmaksymalizować przepływ prądu do rDLPFC (N = 19 uczestników/N = 137 punktów danych), rOTC (N = 12/104) lub lTPJ (N = 12/97). Rzeczywiste pozycje środkowych anod w układach 3 x 1 zostały wyznaczone empirycznie na podstawie danych strukturalnego MRI (obrazy PETRA) pozyskanych przed i po każdej sesji fMRI (Uwaga: współrzędne celowania podano w przestrzeni badanego. W związku z tym skany PETRA poszczególnych uczestników zostały zarejestrowane wspólnie z ich własnym obrazem T1-zależnym uzyskanym podczas skanu bazowego, na którym oparto symulacje pola). Ponieważ identyfikacja pozycji elektrod w 3D jest podatna na błędy, do wyznaczenia sześciu punktów krawędziowych środkowych anod wykorzystano zmodyfikowany dwuwymiarowy widok „pancake” zrekonstruowanej skóry głowy z elektrodami (patrz Rysunek 537,38). Następnie obliczono środek masy, aby uzyskać współrzędne centrum anod dla każdej sesji fMRI i punktu czasowego (tj. przed lub po tDCS). W kolejnym kroku zastosowano odwrotną transformację 2D>3D, aby uzyskać współrzędne elektrod w przestrzeni 3D.
Precyzję przestrzenną rozmieszczenia elektrod z wykorzystaniem neuronawigacji zbadano poprzez porównanie pozycji „zamierzonych” i „rzeczywistych” elektrod25, a wyniki przedstawiono jako odległość L2 Norm (norma euklidesowa)39. L2 Norm to standardowa norma wektorowa obliczana jako pierwiastek kwadratowy z sumy kwadratów składowych wektora. Różnicę między współrzędnymi zamierzonymi a rzeczywistymi obliczono jako różnicę wszystkich składowych wektora, czyli każdego wymiaru przestrzeni 3D (x/y/z). Przesunięcie elektrod w trakcie eksperymentu (tj. dryft) zbadano poprzez porównanie rzeczywistych pozycji elektrod określonych przed i po fMRI.
Dodatkowo porównaliśmy dokładność celowania przy neuronavigacyjnym umiejscawianiu elektrod z danymi z poprzedniego badania naszej grupy, w którym zastosowano podejście celowania oparte na skórze głowy (system EEG 10-2025). Badanie to obejmowało dane od 60 uczestników (31 mężczyzn/29 kobiet, w wieku 22,17 ± 2,33 lat), dwa obszary docelowe (N = 30/30, kora ruchowa/dolny zakręt czołowy) oraz dwie sesje fMRI, co łącznie dało 240 punktów danych. Dane i kod wykorzystane we wszystkich analizach są publicznie dostępne pod adresem https://github.com/LawsOfForm/VerFlu_Simulation_Niemann_2023.
Do analizy (a) odchylenia rzeczywistych pozycji elektrod od zamierzonych lokalizacji oraz (b) dryfu w trakcie sesji fMRI wykorzystano liniowe modele mieszane z następującą formułą, przyjmując normę L² jako zmienną zależną:
yij = β0 + POZYCJAx1ij + REGIONx2ij + POZYCJA * REGIONx3ij + µ0j + εij
yij to odległość w normie L2 między zamierzonym a rzeczywistym położeniem elektrody w POZYCJI i u badanego j. β0 to stały efekt przecięcia, X1 to wartość dla POZYCJI (1: zamierzona, 2: rzeczywista przed fMRI; 3: po fMRI) w pozycji i u badanego j. X2 to wartość dla REGIONU zainteresowania (1: lTPC, 2: rDLPFC, 3: rOTC) w regionie i u badanego j. X3 to wartość dla interakcji POZYCJI i REGIONU dla pozycji i regionu i oraz badanego j, u0 to losowy efekt przecięcia dla badanego j (µ0j ~ N(0, τ002)); εij to residua dla punktu czasowego i u badanego j (εij ~ N(0,σ2)).
Porównanie dokładności celowania dla rozmieszczenia elektrod z nawigacją neuronawigacyjną i metodą opartą na skórze głowy przeprowadzono za pomocą testu t w oparciu o różnice normy L2 między pozycjami zamierzonymi a rzeczywistymi w poszczególnych grupach.
Dokładność rozmieszczenia z neuronavigacją oraz przesunięcia w trakcie eksperymentu
Model referencyjny (zamierzone rozmieszczenie anod centralnych; obszar: lTPC) wykazał istotne odchylenia pozycji elektrod od pozycji zamierzonych zarówno podczas oceny przed-fMRI, jak i po-fMRI (8,71 mm, 95%-CI: [6,57-10,84], p < .001; 8,51 mm, 95%-CI: [6,35-10,66], p < .001) Tabela uzupełniająca S1. Analiza post-hoc dla szacowanych średnich marginalnych (EMMs) została obliczona dla wszystkich głównych i prostych efektów będących przedmiotem zainteresowania (patrz Tabela uzupełniająca S2, Tabela uzupełniająca S3, Rycina 6, Rycina 7, Rycina uzupełniająca S3, oraz Rycina uzupełniająca S4). Przewidywany średni błąd rozmieszczenia dla wszystkich sesji przed-fMRI wynosił 8,54 mm 95%-CI: [7,66-9,42] i różnił się istotnie od pozycji zamierzonych (t (296) = -8,54, p < 0,0001). Sytuacja ta miała miejsce również w przypadku współrzędnych po-fMRI (8,61 mm, 95%-CI: [7,73-9,50]; t (296) = -10,97, p < 0,0001, patrz Rycina 6 i Tabela uzupełniająca S2). Ponieważ błędy rozmieszczenia mogą różnić się w zależności od obszaru25,40,41, błędne rozmieszczenie obliczono również pomiędzy obszarami docelowymi dla każdej pozycji. Jednak nie stwierdzono różnic specyficznych dla obszarów rOTC, lTPC i rDLPFC (Tabela uzupełniająca S3, Rycina 7 oraz Rycina uzupełniająca S4).
Przesunięcie anod centralnych (tzn. dryft), ocenione poprzez porównanie rzeczywistych pozycji elektrod przed i po badaniu fMRI, było pomijalne (-0.08 ± 0.39 mm; t (291) =-0.20, p = 0.98). Szczegóły znajdują się w Rysunku 6 oraz w Tabeli uzupełniającej S2.
Porównanie błędów rozmieszczenia elektrod pomiędzy podejściem z wykorzystaniem neuronawigacji a podejściem opartym na pomiarach skalpu
Aby zbadać istotność rozmieszczenia z wykorzystaniem neuronawigacji, różnice L² Norm pomiędzy pozycjami zamierzonymi a rzeczywistymi w oparciu o rozmieszczenie elektrod pod kontrolą neuronawigacji porównano retrospektywnie z danymi uzyskanymi w poprzednim badaniu, w którym zastosowano podejście celownicze oparte na pomiarach skalpu25. W tym badaniu różnica L² Norm pomiędzy zamierzoną a rzeczywistą pozycją anody wynosiła 12,33 ± 5,67 mm. Bezpośrednie porównanie podejścia z neuronawigacją (8,83 ± 4,68) i podejścia opartego na pomiarach skalpu wykazało istotny spadek dystansu L² Norm o 3,50 mm na korzyść rozmieszczenia elektrod z wykorzystaniem neuronawigacji (t (397) =4,91, p < 0,001).

Rysunek 1: Uproszczony opis sekwencji kroków przedstawiający metodę indywidualnego modelowania prądu. Metoda ta ma na celu wyznaczenie zoptymalizowanych pozycji elektrod na skórze głowy nad obszarem docelowym (prawy DLPFC) przy użyciu ogniskowej konfiguracji tDCS 3 x 1. (A) Anatomiczny obszar docelowy zaznaczony na szaro. (B) Zoptymalizowana pozycja elektrod przedstawiona na zrekonstruowanej głowie uczestnika (anoda = czarny, katoda 1 = czerwony, katoda 2 = zielony, katoda 3 = niebieski. (C) Symulowane pole elektryczne na podstawie indywidualnego modelowania prądu dla tego montażu (w V/m). (D) Współrzędne docelowe dla wszystkich zamierzonych pozycji elektrod wykorzystywanych podczas celowania neuronavigacyjnego (wymiary: x/y/z). Skróty: DLPFC = prawy grzbietowo-boczna kora przedczołowa; tDCS = przezczaszkowa stymulacja prądem stałym. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 2: Sprzęt do neuronawigacji. 1) Tracker badanego, 2) Śrubokręt, 3) Pręt sześciokątny, 4) Gogle, 5) Wskaźnik. UWAGA: Oprogramowanie do neuronawigacji komunikuje się z czujnikiem pozycji (kamerą śledzącą), aby uzyskać lokalizację trackera badanego (zamontowanego na goglach). Wykorzystuje ono macierz rejestracji (uzyskaną poprzez identyfikację punktów orientacyjnych na głowie badanego za pomocą wskaźnika), aby dopasować pozycję głowy badanego ze świata rzeczywistego do przestrzeni obrazu uczestnika (na podstawie obrazu T1-zależnego uczestnika). Przegląd całej konfiguracji przedstawiono na Rysunku dodatkowym S1. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 3: Walidacja obszaru docelowego. Przedstawia procedurę walidacji obszarów docelowych stymulacji za pomocą neuronawigacji. 1) Tracker badanego zamontowany po lewej stronie gogli. 2) Wskaźnik służący do lokalizacji zamierzonej pozycji anody (rDLPFC). 3) Zielony celownik i czerwona kropka (reprezentująca końcówkę wskaźnika) na ekranie komputera. 3a) Numer wskaźnika jest zielony, jeśli wskaźnik odchyla się od celownika o mniej niż 3 mm; 3b) numer wskaźnika jest pomarańczowy, jeśli odchylenie wskaźnika wynosi od 3 do 5 mm; oraz 3c) jest czerwony, jeśli wskaźnik odchyla się o więcej niż 5 mm. Skrót: rDLPFC = prawą grzbietowo-boczna kora przedczołowa. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 4: Wybrane komponenty wielokanałowego stymulatora prądu stałego (DC) oraz *komponenty niestandardowe. (A) 1) Adapter wpinany do gniazda stymulatora DC, 2) kabel LAN obudowy zewnętrznej, 3) *obudowa zewnętrzna (pudełko filtra) umieszczona wewnątrz falowodu skanera podczas tDCS-fMRI. Zacisk jest przymocowany do falowodu w celu połączenia pudełka filtra z pływającą masą ekranowanej komory MRI, aby zredukować potencjalne artefakty. 4) kabel obudowy wewnętrznej, 5) obudowa wewnętrzna z gniazdami dla kabli elektrod, 6) *adapter bezpiecznika łączący trzy oddzielne kanały w jeden dla centralnej anody, 7) *zabezpieczone kable elektrod oraz elektroda centralnej anody (promień 10 mm) wykonana z elastycznego i przewodzącego elektrycznie silikonu, 8) pojedyncza katoda i kabel, 9) przykład elektrody włożonej do wydrukowanego w technologii 3D pomocnika wypełniającego, 10) przykład elektrody umieszczonej w dystansie 3D, 11) taśma miernicza, 12) czepek EEG bez plastikowych wstawek do mocowania elektrod na skórze głowy. W celu zapoznania się z pełną konfiguracją i standardowymi komponentami niewykorzystanymi w tym protokole, czytelników odsyła się do strony internetowej producenta (link w Tabeli Materiałów). (B) Model 3D pomocnika wypełniającego elektrodę. (C) Model 3D dystansu. Skróty: tDCS-fMRI = transcranial direct current stimulation-functional magnetic resonance imaging (przezczaszkowa stymulacja prądem stałym – funkcjonalny rezonans magnetyczny); EEG = elektroencefalografia. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 5: Widok płaski (pancake view) skóry głowy i elektrod na podstawie obrazu PETRA. Kolory reprezentują różne rozmiary spłaszczonej maski skóry głowy (czerwony = 2 mm; zielony = 4 mm; niebieski = 6 mm nadmuchanej skóry głowy dla wizualizacji głębokości głowy). Orientacje przekształconego obrazu są oznaczone białymi literami. Białe kropki reprezentują ręcznie wybrane punkty krawędzi anody centralnej umieszczonej nad prawą grzbieniowo-boczną korą przedczołową. Środek masy obliczono jako współrzędną środka punktów krawędzi przed ponownym przekształceniem w przestrzeń 3D. Skróty: PETRA = Pointwise Encoding Time reduction with Radial Acquisition; L = lewy; R = prawy; P = tylny; A = przedni. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 6: Szacowane średnie marginalne odległości w normie L2 (norma euklidesowa) między zamierzonymi a rzeczywistymi pozycjami elektrody anodowej we wszystkich sesjach i obszarach (rOTC, lTPC, rDLPFC). Oryginalne punkty danych zaznaczono kolorem turkusowym. 95% przedziały ufności EMMs zaznaczono kolorem szarym. Skróty: EMMs = szacowane średnie marginalne; rOTC = prawostronna kora potyliczno-skroniowa; lTPC = lewostronna kora skroniowo-ciemieniowa; rDLPFC = prawostronna grzbietowo-boczna kora przedczołowa; fMRI = funkcjonalny rezonans magnetyczny. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 7: Szacowane średnie marginalne (EMM) odległości w normie L2 (norma euklidesowa) między zamierzonymi a rzeczywistymi pozycjami elektrody anodowej we wszystkich sesjach pomiędzy poszczególnymi obszarami (rOTC, lTPC, rDLPFC) dla każdej pozycji (zamierzonej, przed/po fMRI).Oryginalne punkty danych zaznaczono kolorem turkusowym. 95% przedziały ufności EMM zaznaczono kolorem szarym. Skróty: EMMs = szacowane średnie marginalne; rOTC = prawa kora potyliczno-skroniowa; lTPC = lewa kora skroniowo-ciemieniowa; rDLPFC = prawa grzbietowo-boczna kora przedczołowa; fMRI = funkcjonalny rezonans magnetyczny. Aby wyświetlić powiększoną wersję tego rysunku, kliknij tutaj.
Rycina dodatkowa S1: Przegląd systemu neuronawigacji. 1) Komputer sterujący systemem neuronawigacji z zainstalowanym oprogramowaniem do neuronawigacji; 2) Kamera Polaris; 3) Pedał sterujący; 4) Gogle z trackerem badanego (szczegóły patrz Rycina 2); 5) Wskaźnik. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca S2: Zrekonstruowana skóra z punktami orientacyjnymi zaznaczonymi na żółto. 1) Nasion, 2) lewa nozdrze; 3) prawa nozdrze; 4) lewy dołek przeduszny; 5) prawy dołek przeduszny (patrz krok protokołu 3.1.7). Skróty: LPA = lewy dołek przeduszny; RPA = prawy dołek przeduszny. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający S3: Szacowane średnie marginalne (EMMs) dystansu normy L² (normy euklidesowej) między zamierzonymi a rzeczywistymi pozycjami elektrody anodowej we wszystkich sesjach pomiędzy pozycjami (zamierzoną, przed/po fMRI) dla każdego obszaru (rOTC, lTPC, rDLPFC). Oryginalne punkty danych zaznaczono kolorem turkusowym. 95% przedziały ufności EMMs zaznaczono kolorem szarym. Skróty: EMMs = szacowane średnie marginalne; rOTC = prawa kora potyliczno-skroniowa; lTPC = lewa kora skroniowo-ciemieniowa; rDLPFC = prawa grzbietowo-boczna kora przedczołowa; fMRI = funkcjonalny rezonans magnetyczny. Kliknij tutaj, aby pobrać ten plik.
Rycina dodatkowa S4: Szacowane średnie marginalne (EMM) odległości normy L² (normy euklidesowej) między zamierzonymi a rzeczywistymi pozycjami elektrody anodowej we wszystkich sesjach i pozycjach (zamierzonej, przed/po fMRI) dla każdego obszaru zainteresowania (rOTC, lTPC, rDLPFC). Oryginalne punkty danych zaznaczono kolorem turkusowym. 95% przedziały ufności EMM zaznaczono kolorem szarym. Prawa kora potyliczno-skroniowa (rOTC), lewa kora skroniowo-ciemieniowa (lTPC), prawa grzbietowo-boczna kora przedczołowa (rDLPFC). Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca S1: Wyniki analizy liniowego modelu mieszanego dla odległości wg normy L² (normy euklidesowej) między współrzędnymi zamierzonymi a rzeczywistymi. Przedstawiono współczynniki regresji liniowych modeli mieszanych (modele z losową przeciętnością) oraz dwustronne wartości p. Skróty: CI = przedział ufności; rOTC = prawa kora potyliczno-skroniowa; lTPC = lewa kora skroniowo-ciemieniowa; rDLPFC = prawa grzbietowo-boczna kora przedczołowa; df = stopnie swobody; σ2 = połączona wariancja resztowa; τ = połączona wariancja wyjaśniona przez badanych; ICC = współczynnik korelacji wewnątrzklasowej; N = liczba badanych. Aby pobrać ten plik, kliknij tutaj.
Tabela uzupełniająca S2: Post-hoc EMMs dla analizy modelu mieszanego liniowego dla różnic w normie L² (norma euklidesowa). Kontrasty główne i proste dla głównego efektu pozycji (zamierzona, przed/po fMRI, patrz również Rysunek uzupełniający S2) oraz efekty proste pozycji dla każdego regionu (lTPC, rDLPFC, rOTC, patrz również Rysunek 6). W przypadku efektów prostych przedstawiono wyłącznie efekty istotne. Skróty: EMMs = Estimated Marginal Means (szacowane średnie marginalne); SE = Standard Error (błąd standardowy); df = stopnie swobody; rOTC = prawa kora potyliczno-skroniowa; lTPC = lewa kora skroniowo-ciemieniowa; rDLPFC = prawa grzbietowo-boczna kora przedczołowa. Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca S3: Analiza post-hoc EMMs dla liniowego modelu mieszanego różnicy normy L² (normy euklidesowej). Kontrasty główne i proste dla głównego efektu regionu (lTPC, rDLPFC, rOTC, patrz również Rysunek uzupełniający S3) oraz proste efekty regionu dla każdej pozycji (zamierzonej, przed/po fMRI, patrz również Rysunek 7). Skróty: EMMs = szacowane średnie marginalne; SE = błąd standardowy; df = stopnie swobody; rOTC = prawa kora potyliczno-skroniowa; lTPC = lewa kora skroniowo-ciemieniowa; rDLPFC = prawa grzbietowo-boczna kora przedczołowa. Kliknij tutaj, aby pobrać ten plik.
Kroki krytyczne, potencjalne modyfikacje i rozwiązywanie problemów z metodą
Dokładne pozycjonowanie elektrod jest kluczowym czynnikiem technicznym w eksperymentach tDCS, a odchylenia od zamierzonych pozycji skóry głowy lub dryf elektrod mogą wpływać na przepływ prądu do zamierzonych docelowych obszarów mózgu42,43. Jest to szczególnie istotne w przypadku zogniskowanego tDCS, ponieważ regionalna specyfika podawanego prądu sprawia, że układy te są szczególnie podatne na skutki błędów pozycjonowania 5,25. W obecnym protokole opisaliśmy metodę neuronażowanego, zogniskowanego tDCS podawanego podczas fMRI, która ma na celu zwiększenie dokładności umieszczania indywidualnie określonych pozycji skóry głowy. Kluczowe kroki obejmują (1) wykorzystanie modelowania przepływu prądu elektrycznego w celu optymalizacji pozycji skóry głowy w konfiguracji 3 x 1 dla każdego uczestnika badania, (2) identyfikację zoptymalizowanych pozycji skóry głowy za pomocą neuronawigacji, (3) pozycjonowanie i zabezpieczanie elektrod na skórze głowy przy użyciu standardowych procedur oraz (4) badanie dokładności umieszczania i dryfu elektrod przed i po równoczesnym tDCS-fMRI.
Pierwszy krok w naszym protokole opisuje optymalizację pozycji elektrod dla indywidualnego uczestnika badania naszej trwającej Jednostki Badawczej (RU) za pomocą zindywidualizowanego modelowania prądu (kod do odbudowania optymalizacji można znaleźć tutaj https://github.com/memoslap/SimNIBS_Optimization). Podejście to jest wysoce specyficzne dla projektów RU i obejmowało dwa etapy przygotowawcze, które nie zostały tutaj opisane: (1) Zdefiniowanie anatomicznych regionów docelowych dla ośmiu empirycznych projektów RU (www.memoslap.de) w oparciu o badania fMRI i/lub tDCS, w których wykorzystano te same/podobne zadania motoryczne lub poznawcze; (2) wyrównywanie średniej dawki prądu w różnych projektach poprzez regulację odstępów między anodą środkową a otaczającymi katodami. Wiązało się to z badaniem modelowania komputerowego i danymi obrazowania strukturalnego dostępnymi dla naszej grupy w celu określenia specyficznych dla projektu odstępów między elektrodami. Te specyficzne dla projektu montaże zostały następnie wykorzystane do optymalizacji pozycji skóry głowy dla poszczególnych uczestników projektów RU, aby osiągnąć maksymalne natężenie prądu w odpowiednich regionach docelowych (w tym protokole: po prawej stronie DLPFC). Warto zauważyć, że istnieje wiele możliwych sposobów określenia pożądanych pozycji skóry głowy, a zależą one od celów badania (tj. kontrola regionalna lub dawkowa, celowanie anatomiczne lub funkcjonalne).
Podobnie, istnieje wiele różnych dostępnych na rynku urządzeń do neuronawigacji i stymulatorów prądu stałego, które można wykorzystać do umieszczania neuronów i przeprowadzania badania tDCS-fMRI. W związku z tym modyfikacje procedur opisanych w niniejszym protokole będą wymagane, jeśli używany jest sprzęt różnych producentów lub tDCS-fMRI jest przeprowadzany przy użyciu innego typu skanera. Niemniej jednak uważamy, że protokół można dostosować w elastyczny sposób, aby spełnić wymagania specyficzne dla badania. Na przykład w przypadku konfiguracji 3 x 1 zalecamy użycie wydrukowanej w 3D przekładki, aby zapewnić osiągnięcie i utrzymanie żądanej odległości między elektrodami podczas umieszczania i skanowania. Standaryzacja ilości pasty przewodzącej może być osiągnięta za pomocą pomocy do napełniania elektrod, jak opisano w niniejszym protokole. Aby umożliwić personalizację, modele 3D używane do drukowania zostały udostępnione publicznie (patrz krok 4.1.3.2 protokołu).
Złożoność eksperymentów tDCS fMRI w obrębie skanera jest kolejnym kluczowym czynnikiem, który wymaga od uczestników przejścia do elektrod przymocowanych do głowy. Badania te niosą ze sobą znaczne ryzyko późniejszego przemieszczenia elektrody, nawet przy precyzyjnym ustawieniu5. Dlatego zalecamy rutynowe wdrażanie odpowiednich metod w celu zminimalizowania przemieszczenia i dryfu elektrody (np. stosowanie standardowych nadmuchiwanych poduszek kompatybilnych z MRI, aby uniknąć przemieszczania się elektrod wewnątrz skanera MRI, czepków EEG lub wyprofilowanych pasków utrzymujących elektrody na miejscu oraz innych środków wyszczególnionych w powyższym protokole).
Wreszcie, w porównaniu z behawioralnymi badaniami tDCS, podawanie tDCS podczas fMRI pozwala na weryfikację dokładnych pozycji elektrod na skórze głowy (na podstawie obrazów 3D strukturalnego rezonansu magnetycznego) oraz ilościowe określenie potencjalnych odchyleń od zamierzonych pozycji. Jednak identyfikacja pozycji elektrod w 3D jest podatna na błędy. Dlatego sugerujemy transformację 2D skóry głowy i elektrod, aby ułatwić weryfikację rzeczywistych pozycji elektrod. Potrzebne są jednak przyszłe badania, aby potwierdzić to nowatorskie podejście i porównać wyniki z ekstrakcją danych z obrazów 3D.
Znaczenie metody w stosunku do istniejących innych metod
Metody celowania oparte na skórze głowy (np. system 10-20 EEG44) były stosowane w większości dotychczasowych badań nad tDCS. Wiąże się to z ręczną lub automatyczną identyfikacją anatomicznych punktów orientacyjnych i dodatkowymi pomiarami (np. obwodem głowy, obliczaniem przecięć między punktami orientacyjnymi) w celu określenia zamierzonych pozycji elektrod na skórze głowy. Chociaż podejście to do pewnego stopnia uwzględnia rozmiar głowy poszczególnych uczestników, zaniedbuje inne cechy morfologii mózgu i czaszki, co skutkuje utratą dokładności 5,41,45. Zwiększona precyzja, którą można osiągnąć dzięki neuronawigacji pod kontrolą MRI, jest zilustrowana w sekcji dotyczącej reprezentatywnych wyników niniejszego protokołu. W tym miejscu określiliśmy ilościowo precyzję kierowania sterowanego neuronami, korzystając z danych z trwającego wieloośrodkowego badania tDCS-fMRI i retrospektywnie porównaliśmy dokładność umieszczania celowania sterowanego neuronami z podejściem opartym na skórze głowy. Ponadto zbadaliśmy stopień przemieszczenia elektrod podczas fMRI. Główne wyniki empirycznej części tego protokołu były takie, że (1) celowanie sterowane neuronami jest lepsze niż umieszczanie elektrod na skórze głowy, z około 40% wyższą precyzją w oparciu o normęL2; oraz (2) wystąpił minimalny dryf elektrod podczas sesji fMRI, co podkreśla, że wdrożone kroki techniczne i proceduralne w celu utrzymania elektrod na miejscu (tj. nasadka, przekładka, nadmuchiwane poduszki) zakończyły się sukcesem.
Ograniczenia i przyszłe zastosowania metody
Istnieje kilka czynników ograniczających stosowanie podejścia opisanego w tym protokole. Na przykład optymalizacja zamierzonych pozycji elektrod poprzez modelowanie obliczeniowe, celowanie za pomocą neuronów i weryfikację elektrod wymaga dostępu do danych strukturalnych MRI, specjalistycznej wiedzy technicznej i kosztownego sprzętu, który może być niedostępny. W związku z tym podejście to ogranicza się do wyspecjalizowanych ośrodków badawczych wykonujących badania, które wymagają wysokiej precyzji przestrzennej dostarczania prądu (np. badania eksperymentalne w populacjach neurotypowych badające związki przyczynowo-skutkowe między mózgiem a zachowaniem lub badania modelowania komputerowego badające zależności między dawką a reakcją). W tym kontekście warto wspomnieć, że trafność obecnego modelowania opartego na MRI może być zmniejszona u pacjentów z operacjami głowy lub zmianami postępującymi, co należy wziąć pod uwagę przy stosowaniu metod opisanych w tym protokole do populacji klinicznych.
Ponadto optymalizacja pozycji elektrod poprzez modelowanie obliczeniowe i umieszczanie neuronów w tym protokole ma na celu maksymalizację precyzji anatomicznej. Jak wspomniano powyżej, możliwe są inne sposoby zindywidualizowanego celowania (w tym podejścia inwazyjne ograniczone do określonych populacji klinicznych46), a precyzja anatomiczna niekoniecznie oznacza funkcjonalne znaczenie stymulowanego obszaru. To ostatnie wymagałoby zadania lokalizacyjnego (na poziomie indywidualnego uczestnika lub grupy), uzyskanego przed tDCS-fMRI, w celu zidentyfikowania funkcjonalnie istotnych regionów mózgu do kierowania neuronami. Ponadto przyszłe badania muszą określić funkcjonalne znaczenie przemieszczenia elektrod dla funkcji behawioralnych i neuronalnych, w tym bezpośrednie porównanie efektów między montażami ogniskowymi i konwencjonalnymi.
Należy zauważyć, że symulacje komputerowe opierają się na anatomicznej precyzji (pół)automatycznych segmentacji tkanek uzyskanych z MRI i szacują pole elektryczne na podstawie założeń dotyczących przewodności elektrycznej różnych klas tkanek 5,47. W związku z tym, aby zwiększyć dokładność modeli obliczeniowych, konieczne jest zweryfikowanie tych założeń i ocena dokładności segmentacji. W związku z tym nowe metody badania przewodnictwa przepływu prądu indukowanego przez tDCS w tkance mózgowej podczas jednoczesnych pomiarów tDCS-MRI obejmują obrazowanie gęstości prądu rezonansu magnetycznego (MRCDI) i tomografię impedancji elektrycznej rezonansu magnetycznego (MREIT)5,48,49. Podejścia te mają potencjał w zakresie walidacji symulacji pola elektrycznego tDCS i optymalizacji zindywidualizowanych obliczeń dawki prądu. Połączenie takich zoptymalizowanych podejść do modelowania obliczeniowego z tDCS sterowanym neuronami prawdopodobnie poprawi precyzję i skuteczność tDCS w kontekstach eksperymentalnych i klinicznych w najbliższej przyszłości.
Firma MAN jest członkiem naukowych komitetów doradczych Neuroelectrics i Précis. AH jest częściowo zatrudniony przez neuroConn GmbH. Pozostali autorzy nie mają do zadeklarowania konfliktu interesów.
To badanie zostało sfinansowane przez Niemiecką Fundację Badawczą (granty projektowe: FL 379/26-1; ME 3161/3-1; CRC INST 276/741-2 i 292/155-1, Jednostka badawcza 5429/1 (467143400), FL 379/34-1, FL 379/35-1, FL 379/37-1, Fl 379/22-1, Fl 379/26-1, ME 3161/5-1, ME 3161/6-1, AN 1103/5-1, TH 1330/6-1, TH 1330/7-1). AT był wspierany przez Fundację Lundbeck (grant R313-2019-622). Dziękujemy Sophie Dabelstein i Kirze Hering za pomoc w ekstrakcji danych.
| Nazwa | Firma | Numer katalogowy | Komentarze |
|---|---|---|---|
| System neuronawigacji Brainsight | Brainsight; Rogue Research Inc., Montré al, Kanada | ||
| Wysokotemperaturowa drukarka 3D CR-5 Pro | CREALITY, Shenzhen, Chiny | ||
| DC-STIMULATOR MC | NeuroConn GmbH, Ilmenau, Niemcy | https://www.neurocaregroup.com/technology/dc-stimulator-mc | |
| EMLA Cream 5% | Aspen, Dublin, Irlandia | ||
| MAGNETOM Vida 3T, syngo_MR_XA50 oprogramowanie | Siemens Healthineers AG, Forchheim, Niemcy | ||
| Kamera | PolarisPolaris Vicra; Northern Digital Inc., Waterloo, Kanada | ||
| Przewodząca pasta EEG Ten20 | Weaver and Company, Aurora, Stany Zjednoczone | ||
| Filament do drukarki 3D | TPU SUNLU International, Hong-Kong, Chiny | ||
| Example of alternatives | |||
| Ingenia 3.0T (MR-scanner) | Phillips, Amsterdam, Holandia | ||
| Localite TMS Navigator (sprzęt do neuronawigacji) | Localite, Bonn, Niemcy | ||
| Neural Navigator (sprzęt do neuronawigacji) | Soterix, New Jersey, USA | ||
| Filament do drukarki 3D PEBA | Kimya, Nantes, Francja | ||
| Filament do drukarki 3D PLA | Filamentworld, Neu-Ulm, Deutschland | ||
| StarStim (stymulator) | Neuroelektryka, Barcelona, Hiszpania |