Podsumowując, przedstawiamy protokół umożliwiający monitorowanie neuromodulacji FUS in vitro z wykorzystaniem neuronów hodowanych z HiPSC. Ogólna platforma systemowa służąca do stymulacji neuronów pochodzących z HiPSC oraz rejestrowania odpowiadających im odpowiedzi elektrycznych do analizy została przedstawiona na Rysunku 1. Niniejsze badanie koncentruje się na stymulacji neuronów za pomocą FUS i rejestrowaniu odpowiedzi elektrycznych w systemie MEA, co pokazano na Rysunku 2. Komponenty peryferyjne systemów FUS i MEA oraz ich połączenia zostały zilustrowane na Rysunku 3.
Charakteryzację ogniska wykonuje się przed eksperymentami neuronalnymi, aby upewnić się, że dno dołka jest w pełni pokryte ogniskiem FUS. W celu oceny systemu FUS należy przeprowadzić wizualizację plamki ogniskowej na arkuszach termochromowych, jak pokazano na Ryc. 4. Po charakteryzacji plamki ogniskowej należy przeprowadzić etapy postprocessingu, w tym filtrowanie, progowanie i obliczanie częstotliwości wyładowań, które zostały podsumowane na Ryc. 5 oraz Ryc. 6. Kroki te są niezbędne do odzyskania użytecznych sygnałów poprzez odfiltrowanie szumów z otoczenia, a tym samym do uzyskania wglądu w zmiany aktywności neuronalnej wywołane przez FUS. Wykresy rastrowe na Ryc. 6A-B przedstawiają wykryte kolce w każdym kanale. Ponieważ całe dno dołka znajduje się w punkcie ogniskowym przetwornika FUS, oczekuje się, że FUS zmieni częstotliwość wyładowań we wszystkich elektrodach. Zmiana ta jest przedstawiona na wykresie częstotliwości wyładowań na Ryc. 6C, który pokazuje, że wybrane parametry stymulacji doprowadziły do zwiększenia częstotliwości wyładowań neuronalnych. Konkretnie, częstotliwość wyładowań przed FUS (tzn. bazowa) wynosiła 140 Hz ± 116,7 Hz, podczas gdy częstotliwość wyładowań po FUS przy zastosowaniu ciągłej fali FUS wynosiła 786 Hz ± 419,4 Hz. Dodatkowo Ryc. 6C pokazuje, w jaki sposób zmiana parametrów FUS (np. zastosowanie falowanej fali FUS zamiast ciągłej) może zmienić wielkość zmiany częstotliwości wyładowań, a także czas potrzebny na powrót neuronów do stanu bazowego. Niskointensywny ultradźwięk skupiony (LIFU) nie powoduje znaczącego ogrzewania kultur, szczególnie w porównaniu z wysokointensywnym ultradźwiękiem skupionym, którego celem jest uzyskanie zmiany termicznej. Brak klinicznie istotnej zmiany temperatury jest potwierdzony obliczeniami teoretycznymi i symulacjami (Ryc. uzupełniająca 2). Nawet w ekstremalnych przypadkach parametrów eksperymentalnych FUS wymienionych w Tabeli 1, można było zaobserwować jedynie minimalny wzrost temperatury o około 0,04 °C.
Zastosowanie wykresu częstości wyładowań pozwala na ilościowe określenie efektów neuromodulacyjnych FUS i może służyć do rozróżnienia odpowiedzi pobudzających od hamujących. Istotną zaletą wielodołkowej płytki MEA jest możliwość jej wielokrotnego wykorzystania do badania różnych stanów neuronalnych oraz parametrów stymulacji w sposób wysokoprzepustowy.

Rysunek 1: Przegląd platformy in vitro do neuromodulacji neuronów w studience za pomocą skupionego ultradźwięku (FUS) oraz pomiaru ich aktywności neuronalnej z wykorzystaniem matrycy mikrolektrod. Każda elektroda (linie czerwona, zielona i niebieska) rejestruje sygnały z populacji neuronów w obrębie jednej studienki. Zaimplementowano potok przetwarzania w celu przekształcenia surowych rejestracji elektrycznych neuronów w detekcję wzorców wyładowań neuronalnych. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Neuromodulacja FUS z użyciem wielodołkowej macierzy mikroelektrod (MEA). (A) Schemat układu do neuromodulacji FUS z wielodołkową macierzą MEA. Fale akustyczne generowane przez przetwornik FUS rozprzestrzeniają się przez stożek FUS wypełniony odgazowaną wodą i są sprzęgane za pomocą żelu ultradźwiękowego. Parafilm jest przymocowany do dołka za pomocą gumki recepturki, aby zapobiec kontaminacji. Płytka MEA przesyła zapisy elektryczne z neuronów do systemu MEA. (B) Zdjęcie przetwornika FUS na płytce wielodołkowej umieszczonej w systemie MEA. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Konfiguracja platformy in vitro. (A) Widok z przodu konfiguracji platformy in vitro. Wyjście mocy przetwornika (TPO; po lewej) służy do programowania parametrów FUS. System MEA (po prawej) rejestruje aktywność elektryczną neuronów w płytce wielodołkowej, które są neuromodulowane przez przetwornik FUS. (B) Widok z tyłu konfiguracji platformy in vitro z połączeniami z sieci dopasowującej (1) do TPO oraz (2) do przetwornika. (3) Połączenie z systemu MEA do TPO synchronizuje akwizycję danych. (4) Połączenie z systemu MEA do komputera w celu przesyłania danych. (5) Połączenie zasilania systemu MEA. (6) Połączenie zasilania systemu FUS. (7) Przycisk sonikacji. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Charakterystyka przetwornika FUS. (A) Mapa ciśnienia ogniska przy użyciu parametrów FUS szczegółowo opisanych w Tabeli 1, zmierzona za pomocą systemu AMPLITUDE15. (B) Stan przed i po sonikacji arkusza termochromowego umieszczonego na dnie dołka przy użyciu układu eksperymentalnego przedstawionego na Rycynie 3. Arkusz termochromowy zmienia kolor w odpowiedzi na zmiany temperatury, co zapewnia wizualną walidację skutecznej stymulacji w miejscu położenia neuronów. Maksymalna średnia intensywność impulsu w szczycie przestrzennym (ISPPA) wynosząca 30 W/cm2 oraz ciągła sonikacja przez 3 min zostały dobrane tak, aby drastycznie zmienić temperaturę lokalną w celu lepszej wizualizacji takiego punktu ogniskowania. Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.

Rysunek 5: Potok przetwarzania. Krok 1: Surowe zapisy elektryczne są rejestrowane z N = 16 kanałów. Kolejne kroki przedstawiają proces z wykorzystaniem kanału 16 (obramowanego na czerwono). Krok 2: Dla każdego kanału stosowany jest filtr pasmowoprzepustowy Butterwortha (pasmo od 5 Hz do 3 kHz), a następnie filtr Gaussa (σ = 3). Próg jest ustawiany na pięciokrotność odchylenia standardowego sygnału w oknie 2 s wycentrowanym na początku sonikacji. Krok 3: Sygnały powyżej lub poniżej progu są klasyfikowane jako kolce. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 6: Wykresy rastrowe i częstotliwości wyładowań. (A) Wykres rastrowy wykrytych impulsów w każdym kanale w funkcji czasu sonikacji. Czas stymulacji FUS zaznaczono czerwoną linią. (B) Wykres rastrowy neuronów przy różnych ustawieniach FUS wraz z ciągłym FUS dla porównania. (C) Częstotliwość wyładowań obliczono przy użyciu przesuwnego okna o szerokości 50 ms. Średnie częstotliwości wyładowań przed i po neuromodulacjach FUS wynosiły odpowiednio 140 Hz i 786 Hz. Przy pulsacyjnym FUS średnie częstotliwości wyładowań wynosiły 230 Hz i 540 Hz. Zaobserwowano, że ten zestaw zmiennych stymulacji FUS indukuje krótszą aktywację i mniejszą zmianę częstotliwości. Proces obliczania częstotliwości wyładowań szczegółowo opisano na Rysunku uzupełniającym 1. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
| Parametr | Wartość |
| Maksymalna moc/kanał | 1,200 W |
| Przeczywisty | 0,749 W/kanał. |
| ISPPA | 10,79 W/cm2 |
| ISPTA | 0,05 W/cm2 |
| Długość serii | 0,100 ms |
| Częstotliwość | 250,00 kHz |
| Ogniskowanie | 39,800 mm |
| Okres | 20 000 ms |
| minutnik | 60,000 s |
Tabela 1: Parametry skupionej ultradźwiękowej stymulacji (FUS) ustawione na urządzeniu TPO dla badania przedstawionego na Rysunku 4.
Rysunek uzupełniający 1: Przetwarzanie wykresu rastrowego na częstotliwość wyładowań. Krok 1: Policz kolce we wszystkich kanałach, aby uzyskać liczbę zdarzeń w określonym przesuwnym oknie czasowym. Uwaga: W tym przypadku, dla lepszej ilustracji, wybrano większe okno przesuwne (ustawione na 0.1 s). Krok 2: Przelicz liczbę kolców na długość okna na liczbę kolców na sekundę (np. tutaj pomnożono liczbę zdarzeń przez 10, aby przeliczyć je na herce [Hz], a następnie podzielono przez 1 000, aby uzyskać wartość w kilohertcach [kHz]). Krok 3: Krzywa częstotliwości wyładowań uzyskana w wyniku tych działań. Otwarty zestaw narzędzi wraz z próbkami danych jest dostępny na GitHubie (https://github.com/Rxliang/FUSNeuromod). Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający 2: Profil temperatury z wyników symulacji K-wave dla LIFU16. Na podstawie mapy intensywności akustycznej przedstawionej na Rysunku 4, wynik symulacji K-wave sugeruje maksymalny wzrost temperatury o 0,04 °C w obszarze centralnym strefy ogniskowania (promień: 2 mm), przy zastosowaniu ekstremalnych parametrów eksperymentalnego FUS wymienionych w Tabeli 1. Kliknij tutaj, aby pobrać ten plik.