Name: Celowanie w szlaki włókien nerwowych w celu terapii głębokiej stymulacji mózgu przy użyciu interaktywnych, specyficznych dla pacjenta modeli (Wideo) | JoVE Journal
Uploaded: 2018-08-12T12:30:00.000+00:00
Duration: 14 min 14 s
Description: Poznaj interaktywne modelowanie głębokiej stymulacji mózgu przy użyciu danych MRI i oprogramowania SCIRun. Ta metoda umożliwia symulację w czasie rzeczywistym wpływu stymulacji na szlaki nerwowe u pacjentów z chorobą Parkinsona.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

Ten projekt był wspierany przez National Institute of Health Grants UH3, NS095554. Wsparcie techniczne zostało zapewnione przez Centrum Integratywnych Obliczeń Biomedycznych w Naukowym Instytucie Informatyki i Obrazowania i było możliwe częściowo dzięki oprogramowaniu opracowanemu przez NIH P41-GM103545, Center of Integrative Biomedical Computing.
Podziękowania kierujemy do Lexie Floor i Nathana Galli z Scientific Computing and Imaging Institute za produkcję i edycję filmu z przesłaniem, a także do Theresy Lins za pomoc w przygotowaniu manuskryptu.

1. Przetwarzanie obrazu
<ol>
	<li>Segmentacja rezonansu magnetycznego T1
	<ol>
		<li>Pobierz i zainstaluj FreeSurfer17. Zakłada się, że uzyskano rezonans magnetyczny T1 i jest on w formacie DICOM lub NIFTI.</li>
		<li>Wpisz następujące polecenie, aby utworzyć katalog pacjentów i dodać ich wolumin T1 w FreeSurfer: recon-all -s nazwa_pacjenta -i /Full/Path/To/nii</li>
		<li>Wpisz następujące polecenie, aby uruchomić automatyczną segmentację FreeSurfer: recon-all -s nazwa_pacjenta -all</li>
		<li>Użyj mri_convert, aby przekonwertować "aseg.auto.mgz" z zastrzeżonego formatu .mgz FreeSurfer.</li>
	</ol>
	</li>
	<li>Obrazowanie dyfuzyjne ważone (DWI)
	<ol>
		<li>Pobierz i zainstaluj Slicer18. Zakłada się, że akwizycja DWI została przeprowadzona i jest dostępna jako stos plików DICOM. 
		UWAGA: Dobrym narzędziem wiersza poleceń do wykonania tej rekonstrukcji jest FMRIB Software Library (FSL)19. FSL został użyty w tym protokole do wykonania korekcji zniekształceń echo-płaskich, ruchowych i prądów wirowych przed rekonstrukcją. Dane DWI użyte w tej analizie to seria plików DICOM dla 41 różnych kierunków gradientu uwrażliwiających na dyfuzję.</li>
		<li>Użyj modułu DWIConverter w narzędziu Slicer, aby zrekonstruować 41 oddzielnych nabytków w jednym woluminie. Określ katalog danych wejściowych Dicom, w którym znajdują się dane DWI. Zaznacz przycisk radiowy DicomToNrrd i podaj nazwę wyjściowego woluminu DWI.</li>
		<li>Utwórz maskę tensorową za pomocą modułu Dyfuzja ważona maskowaniem objętości, aby wyeliminować szum tła z oszacowania DTI. Pozostaw parametr threshold na wartość domyślną 0,5.</li>
		<li>Przekonwertuj ten pojedynczy wolumin DWI na wolumin tensora dyfuzji (DTI) za pomocą oszacowania DWI do DTI, używając maski utworzonej w poprzednim kroku jako maska tensora dyfuzji. Ustaw parametry estymacji na Najmniejsze kwadraty i upewnij się, że opcja Przesuń ujemne wartości własne jest zaznaczona.</li>
		<li>Zapisz nowo utworzony wolumin DTI w formacie NRRD.</li>
	</ol>
	</li>
</ol>
2. Generowanie modelu elementów skończonych
<ol>
	<li>SCIRun
	<ol>
		<li>Pobierz i zainstaluj SCIRun w wersji 5 (http://www.sci.utah.edu/software/scirun.html). SCIRun to środowisko rozwiązywania problemów do modelowania, symulacji i wizualizacji problemów naukowych. Jest to obliczeniowe środowisko pracy używane do tworzenia, rozwiązywania i interakcji z modelem opracowanym w tym projekcie. 
		UWAGA: Sieci SCIRun są budowane poprzez łączenie sekwencji poszczególnych modułów, które wykonują określone zadania. Przekazywanie danych przez sieć złożoną z prostych modułów pozwala na bardziej złożone przetwarzanie i symulację.</li>
	</ol>
	</li>
	<li>Budowa geometrii elektrody 
	UWAGA: Elektroda modelowana w tym projekcie to Medtronic 3387 DBS lead20. Jest to cylindryczny przewód z czterema stykami o wysokości 1,5 mm, średnicy 1,27 mm i odstępach 1,5 mm od siebie. Materiał nieprzewodzący oddziela cztery styki.
	<ol>
		<li>Utwórz siatkę o zamkniętej powierzchni dla każdego komponentu w wyprowadzeniu DBS pasującą do podanych wymiarów. Utwórz cały potencjalny klient, używając dwóch geometrii pierwotnych, cylindrów i sfer. 
		UWAGA: Geometrię wiodącego można tworzyć w różnych programach do modelowania 3D.</li>
		<li>Utwórz siatkę powierzchniową z warstwy enkapsulacyjnej o grubości 0,5 mm otaczającej cały wyprowadzenie.</li>
	</ol>
	</li>
	<li>Tworzenie całej siatki mózgu
	<ol>
		<li>Załaduj powierzchnię mózgu utworzoną w części 1, aby służyła jako zewnętrzna granica MES.</li>
		<li>Utwórz dwie koncentryczne powierzchnie pudełkowe otaczające przewód DBS, aby kontrolować gęstość siatki wokół elektrody. 
		UWAGA: Powierzchnia przewodu DBS i skrzynki powinny być wyrównane wzdłuż dodatniej osi z dolną końcówką wału elektrody znajdującą się w początku. Jest to ważne dla późniejszej rotacji i translacji elektrody.</li>
		<li>Utwórz chmurę punktów z pojedynczym punktem znajdującym się w każdym odrębnym regionie siatki. Obszary siatki to: segmenty trzonu, cztery kontakty, pudełko i powierzchnia mózgu. Wartości dołączone do tych punktów reprezentują maksymalne ograniczenia objętości dla każdego elementu czworościennego w siatce dla tego regionu. Ograniczenia te są dostosowywane w celu kontrolowania gęstości siatki w elektrodzie i wokół niej.</li>
		<li>Użyj modułu "InterfaceWithTetGen", aby stworzyć całą siatkę mózgu. Ustaw porty wejściowe w następujący sposób: 
		Port 1: powierzchnia mózgu (od 1.1) 
		Port 2: chmura punktów atrybutu regionalnego, ograniczenia objętości (od 2.3.3) 
		Port 3: none 
		Port 4: powierzchnie modelu elektrody (od 2.2)</li>
	</ol>
	</li>
	<li>Interaktywne umieszczanie elektrod 
	UWAGA: Celem tego kroku jest umożliwienie łatwego przemieszczania elektrody w dowolne miejsce w mózgu i automatyczna aktualizacja systemu o nową siatkę i symulację.
	<ol>
		<li>Użyj dwóch kaskadowych modułów CreateGeometricTransform, aby obrócić elektrodę w płaszczyźnie X i Y. Kliknij Obróć i ustaw pierwszy moduł na 1.00 na osi Y, a drugi na 1.00 na osi X. Kąt obrotu ustawia się za pomocą dolnego suwaka, Obróć Theta (stopnie).</li>
		<li>Połącz moduł TransformMeshWithTransform dla każdego przekształcenia i obiektu, który musi zostać przekształcony.</li>
		<li>Utwórz moduł GenerateSinglePointProbeFromField i połącz jego port geometrii z modułem ViewScene. Przesuń sferę w oknie ViewScene, przytrzymując Shift i klikając lewym przyciskiem myszy. Moduł sondy punktowej aktualizuje położenie kuli po zwolnieniu. 
		UWAGA: To położenie służy do określenia transformacji translacyjnej zastosowanej do geometrii elektrody.</li>
		<li>Użyj modułu TransformMeshWithTransform, aby przetłumaczyć każdą część geometrii elektrody z lokalizacją sondy punktowej.</li>
		<li>Zastosuj dwie transformacje obrotu i translację z widżetu sondy punktowej do każdego obiektu z wyjątkiem powierzchni mózgu. 
		UWAGA: Wszystkie trzy z tych transformacji powinny być wykonane na każdym obiekcie, zanim zostanie on przesłany do modułu InterfaceWithTetGen, aby algorytm tworzenia siatki wiedział, jak zasiać geometrię elektrody wewnątrz mózgu w nowym miejscu i orientacji. Za każdym razem, gdy elektroda zostanie przesunięta w mózgu, zostanie zbudowana nowa siatka.</li>
	</ol>
	</li>
</ol>
3. Obliczanie pola bioelektrycznego
<ol>
	<li>Ustawianie parametrów przewodności
	<ol>
		<li>Załaduj wolumin DTI utworzony w kroku 1 i ustaw dane na wyjściu siatki czworościennej z InterfaceWithTetGen utworzonego w kroku 2 przy użyciu MapFieldDataFromSourceToDestination z opcją linear('weighted'). 
		UWAGA: Na tym etapie można zaimplementować dowolne oszacowanie przewodności anizotropowej.</li>
		<li>Ustaw przewodność obszaru wału na 1e-6, a obszar kontaktu na 1e6 za pomocą CalculateFieldData. Po ustawieniu wszystkich wartości przewodności przelej pole do BuildFEMatrix.</li>
	</ol>
	</li>
	<li>Ustawianie źródeł prądu/napięcia i warunków brzegowych
	<ol>
		<li>Przelej powierzchnię mózgu używaną w procesie tworzenia siatki z kroku 2 do SetFieldDataToConstantValue i ustaw wartość na 0. Spowoduje to pochłanianie napięcia na całej powierzchni mózgu. Następnie użyj polecenia InsertVoltageSource, aby umieścić wartości napięcia na wyjściu siatki czworościennej z InterfaceWithTetGen.</li>
		<li>W przypadku stymulacji monopolarnej należy utworzyć pojedynczy punkt w środku jednej z lokalizacji kontaktu, obliczając wartość FieldCenter z modułu ReportFieldInfo na jednej z powierzchni kontaktowych. Ustaw dane pola w tym punkcie na -1 wolt. Następnie użyj InsertVoltageSource, aby umieścić źródło punktowe w siatce czworościennej. 
		UWAGA: Istnieją teraz dwie kopie siatki czworościennej. Jeden z wartościami napięcia ustawionymi na powierzchni mózgu, a drugi ze źródłem punktowym wewnątrz kontaktu.</li>
		<li>Połącz dane wyjściowe dwóch modułów InsertVoltageSource i przelej dane wyjściowe do ApplyFEMVoltageSource wraz z danymi wyjściowymi BuildFEMatrix.</li>
	</ol>
	</li>
	<li>Rozwiązywanie problemu z przechodzeniem do przodu 
	UWAGA: Dane wyjściowe ApplyFEMVoltageSource zawierają macierz sztywności na porcie 1 i macierz po prawej stronie na porcie 2. Są to dwa składniki potrzebne do rozwiązania układu równań.
	<ol>
		<li>Przelej dwa wyjścia ApplyFEMVoltageSource do SolveLinearSystem. Wybierz metodę rozwiązywania gradientu sprzężonego i Jacobiego jako kondycjoner wstępny. Wyjściem tego modułu jest macierz rozwiązań, wartość napięcia w każdym węźle w siatce czworościennej.</li>
		<li>Użyj SetFieldData, aby umieścić macierz rozwiązania napięcia na siatce czworościennej w celu uzyskania reprezentacji 3D rozkładu napięcia.</li>
	</ol>
	</li>
	<li>Wizualizacja Isosurface
	<ol>
		<li>Potok wyniku SetFieldData do ExtractSimpleIsosurface. Otwórz interfejs użytkownika na karcie Ilość i ustaw liczbę równomiernie rozmieszczonych wartości ISO na 10.</li>
		<li>Użyj modułu ShowField i połącz go z modułem ViewScene, aby wyświetlić Isosurfaces jako przezroczyste z mapą kolorów tęczy przeskalowaną do min i max rozwiązania napięcia.</li>
	</ol>
	</li>
</ol>

<ol><li>Benabid, A. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Chronic+electrical+stimulation+of+the+ventralis+intermedius+nucleus+of+the+thalamus+as+a+treatment+of+movement+disorders%5BTitle%5D)+AND+%22Journal+of+Neurosurgery%22%5BJournal%5D)">Chronic electrical stimulation of the ventralis intermedius nucleus of the thalamus as a treatment of movement disorders.</a> Journal of Neurosurgery. 84 (2), 203-214 (1996).</li>
<li>Limousin, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Effect+of+parkinsonian+signs+and+symptoms+of+bilateral+subthalamic+nucleus+stimulation%5BTitle%5D)+AND+%22The+Lancet%22%5BJournal%5D)">Effect of parkinsonian signs and symptoms of bilateral subthalamic nucleus stimulation.</a> The Lancet. 345 (8942), 91-95 (1995).</li>
<li>Schiff, N. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Behavioural+improvements+with+thalamic+stimulation+after+severe+traumatic+brain+injury%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Behavioural improvements with thalamic stimulation after severe traumatic brain injury.</a> Nature. 448 (7153), 600-603 (2007).</li>
<li>Vandewalle, V., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Stereotactic+treatment+of+Gilles+de+la+Tourette+syndrome+by+high+frequency+stimulation+of+thalamus%5BTitle%5D)+AND+%22Lancet%22%5BJournal%5D)">Stereotactic treatment of Gilles de la Tourette syndrome by high frequency stimulation of thalamus.</a> Lancet. 353 (9154), 724(1999).</li>
<li>Mayberg, H. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Deep+brain+stimulation+for+treatment-resistant+depression%5BTitle%5D)+AND+%22Neuron%22%5BJournal%5D)">Deep brain stimulation for treatment-resistant depression.</a> Neuron. 45 (5), 651-660 (2005).</li>
<li>Hashimoto, T., Elder, C. M., Okun, M. S., Patrick, S. K., Vitek, J. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Stimulation+of+the+subthalamic+nucleus+changes+the+firing+pattern+of+pallidal+neurons%5BTitle%5D)+AND+%22Journal+of+Neuroscience%22%5BJournal%5D)">Stimulation of the subthalamic nucleus changes the firing pattern of pallidal neurons.</a> Journal of Neuroscience. 23 (5), 1916-1923 (2003).</li>
<li>York, M. K., Wilde, E. A., Simpson, R., Jankovic, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Relationship+between+Neuropsychological+Outcome+and+DBS+Surgical+Trajectory+and+Electrode+Location%5BTitle%5D)+AND+%22J.+Neurol.+Sci%22%5BJournal%5D)">Relationship between Neuropsychological Outcome and DBS Surgical Trajectory and Electrode Location.</a> J. Neurol. Sci. 287 (1-2), 159-171 (2009).</li>
<li>Machado, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Deep+brain+stimulation+for+Parkinson%27s+disease%3A+surgical+technique+and+perioperative+management%5BTitle%5D)+AND+%22Movement+Disorders%22%5BJournal%5D)">Deep brain stimulation for Parkinson's disease: surgical technique and perioperative management.</a> Movement Disorders. 21, S247-S258 (2006).</li>
<li>Volkmann, J., Moro, E., Pahwa, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Basic+algorithms+for+the+programming+of+deep+brain+stimulation+in+Parkinson%27s+disease%5BTitle%5D)+AND+%22Movement+Disorders%22%5BJournal%5D)">Basic algorithms for the programming of deep brain stimulation in Parkinson's disease.</a> Movement Disorders. 21 (S14), S284-S289 (2006).</li>
<li>Contarino, M. F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Directional+steering%3A+A+novel+approach+to+deep+brain+stimulation%5BTitle%5D)+AND+%22Neurology%22%5BJournal%5D)">Directional steering: A novel approach to deep brain stimulation.</a> Neurology. 83 (13), 1163-1169 (2014).</li>
<li>Pollo, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Directional+deep+brain+stimulation%3A+An+intraoperative+double-blind+pilot+study%5BTitle%5D)+AND+%22Brain%22%5BJournal%5D)">Directional deep brain stimulation: An intraoperative double-blind pilot study.</a> Brain. 137 (7), 2015-2026 (2014).</li>
<li>Willsie, A. C., Dorval, A. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fabrication+and+initial+testing+of+the+%C2%B5DBS%3A+a+novel+deep+brain+stimulation+electrode+with+thousands+of+individually+controllable+contacts%5BTitle%5D)+AND+%22Biomedical+Microdevices%22%5BJournal%5D)">Fabrication and initial testing of the µDBS: a novel deep brain stimulation electrode with thousands of individually controllable contacts.</a> Biomedical Microdevices. 17 (56), 9961(2015).</li>
<li>Oliveria, S. F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Safety+and+efficacy+of+dual-lead+thalamic+deep+brain+stimulation+for+patients+with+treatment-refractory+multiple+sclerosis+tremor%3A+a+single-centre%2C+randomised%2C+single-blind%2C+pilot+trial%5BTitle%5D)+AND+%22The+Lancet+Neurology%22%5BJournal%5D)">Safety and efficacy of dual-lead thalamic deep brain stimulation for patients with treatment-refractory multiple sclerosis tremor: a single-centre, randomised, single-blind, pilot trial.</a> The Lancet Neurology. 16 (9), 691-700 (2017).</li>
<li>Butson, C. R., Cooper, S. E., Henderson, J. M., Wolgamuth, B., Mcintyre, C. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Probabilistic+Analysis+of+Activation+Volumes+Generated+During+Deep+Brain+Stimulation%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">Probabilistic Analysis of Activation Volumes Generated During Deep Brain Stimulation.</a> Neuroimage. 54 (3), 2096-2104 (2011).</li>
<li>Basser, P. J., Mattiello, J., LeBihan, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MR+diffusion+tensor+spectroscopy+and+imaging%5BTitle%5D)+AND+%22Biophys.+J%22%5BJournal%5D)">MR diffusion tensor spectroscopy and imaging.</a> Biophys. J. 66 (1), 259-267 (1994).</li>
<li>Gullmar, D., Haueisen, J., Reichenbach, J. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Influence+of+anisotropic+electrical+conductivity+in+white+matter+tissue+on+the+EEG%2FMEG+forward+and+inverse+solution.+A+high-resolution+whole+head+simulation+study%5BTitle%5D)+AND+%22NeuroImage%22%5BJournal%5D)">Influence of anisotropic electrical conductivity in white matter tissue on the EEG/MEG forward and inverse solution. A high-resolution whole head simulation study.</a> NeuroImage. 51 (1), 145-163 (2010).</li>
<li>Fischl, B. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/22248573">FreeSurfer.</a> Neuroimage. 62 (2), 774-781 (2012).</li>
<li>Fedorov, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((3D+Slicer+as+an+Image+Computing+Platform+for+the+Quantitative+Imaging+Network%5BTitle%5D)+AND+%22Magn.+Reson.+Imaging%22%5BJournal%5D)">3D Slicer as an Image Computing Platform for the Quantitative Imaging Network.</a> Magn. Reson. Imaging. 30 (9), 1323-1341 (2012).</li>
<li>Jenkinson, M., Beckmann, C. F., Behrens, T. E., Woolrich, M. W., Smith, S. M. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/21979382">FSL.</a> Neuroimage. 62 (2), 782-790 (2012).</li>
<li><a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=%22Medtronic+DBS+3387%2F3389+Lead+Kit+Manual%22">Medtronic DBS 3387/3389 Lead Kit Manual</a>. , (2018).</li>
<li>Baker, J. L., et al. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/27582298">Robust modulation of arousal regulation, performance and frontostriatal activity through central thalamic deep brain stimulation in healthy non-human primates.</a> Journal of Neurophysiology. 116 (5), 2383-2404 (2016).</li>
<li>Fregni, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Treatment+of+major+depression+with+transcranial+direct+current+stimulation%5BTitle%5D)+AND+%22Bipolar+Disorders%22%5BJournal%5D)">Treatment of major depression with transcranial direct current stimulation.</a> Bipolar Disorders. 8 (2), 203-204 (2006).</li>
<li>Hodaie, M., Wennberg, R. A., Dostrovsky, J. O., Lozano, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Chronic+anterior+thalamus+stimulation+for+intractable+epilepsy%5BTitle%5D)+AND+%22Epilepsia%22%5BJournal%5D)">Chronic anterior thalamus stimulation for intractable epilepsy.</a> Epilepsia. 43 (6), 603-608 (2002).</li>
<li>Rosenow, F., Lüders, H. <a target="_blank" href="https://academic.oup.com/brain/article/124/9/1683/303186">Presurgical evaluation of epilepsy.</a> Brain. 124 (9), 1683-1700 (2001).</li>
</ol>

Dr Christopher R. Butson pełnił funkcję konsultanta w firmach NeuroPace, Advanced Bionics, Boston Scientific, Intelect Medical, St. Jude Medical i Functional Neuromodulation.

W protokole tym zademonstrowano techniki mające na celu zmniejszenie obciążenia czasowego i kosztów obliczeniowych związanych z tworzeniem modeli specyficznych dla pacjenta dla DBS w stopniu, który umożliwia uzyskanie informacji zwrotnej o wynikach symulacji w czasie zbliżonym do rzeczywistego. Szybkie sprzężenie zwrotne umożliwia eksplorację dużej przestrzeni parametrów w celu lepszego zrozumienia, w jaki sposób zmiany parametrów wpływają na wyniki modelu. Parametry te obejmują lokalizację przewodu, wybór aktywnych styków oraz amplitudę, szerokość impulsu i częstotliwość kształtu fali stymulacji na tych stykach. Kluczowe cechy proponowanego narzędzia to: 1) prosty interfejs użytkownika do dostosowywania parametrów modelu z wizualizacją w czasie zbliżonym do rzeczywistego, w jaki sposób te parametry wpływają na symulację, oraz 2) automatyzacja tworzenia modelu na podstawie niewielkiego zestawu danych wejściowych: powierzchni mózgu, tensorów przewodnictwa tkanki mózgowej i powierzchniowej reprezentacji geometrii elektrody. Ta automatyzacja przyspiesza tworzenie modeli specyficznych dla pacjenta dla wielu osób, które mają różną geometrię mózgu i przewodnictwo tkanek, a także ocenę skutków wprowadzenia różnych wzorów elektrod do istniejących modeli. Etapy wstępnego przetwarzania obrazu opisane w tym protokole nie były w pełni zautomatyzowane i mogą trwać do jednego dnia. Jednak po zakończeniu dane wygenerowane na tych etapach są uważane za statyczne, co oznacza, że dane te nie są modyfikowane podczas symulacji. Automatyzacja tworzenia modeli wynika ze zdolności systemu do zastosowania tych danych do MES specyficznego dla pacjenta bez wysiłku ręcznego. Sieć SCIRun do wykonywania wszystkich kroków generowania modelu, symulacji i wizualizacji musi być zbudowana tylko raz. W związku z tym należy ponownie wykonać tylko etapy wstępnego przetwarzania obrazu, aby wygenerować model specyficzny dla pacjenta przychodzącego.
Wzrost wydajności generowania wyników z potoku modelowania wynika z integracji generowania siatki, obliczeń pola bioelektrycznego i wizualizacji rozwiązania w jednym środowisku oprogramowania. Istniejące techniki modelowania, takie jak adaptacyjne generowanie siatki, zostały wykorzystane do uzyskania większej gęstości siatki wokół elektrody i mniejszej gęstości dalej od elektrody, co skraca czas budowy i rozwiązania MES. Oprogramowanie SCIRun umożliwia również automatyzację generowania siatek i obliczeń pola bioelektrycznego. Zdefiniowany przez użytkownika ruch elektrody za pomocą interaktywnych widżetów uruchamia budowę nowej siatki ze zaktualizowanym położeniem elektrody. Obejmuje to modyfikację warunków brzegowych i wartości przewodności dla nowego położenia elektrody.
Geometria elektrody jest traktowana jako swobodnie poruszający się obiekt wewnątrz objętości mózgu, zanim jej położenie zostanie zintegrowane z MES. Kluczową implikacją tego podejścia dla konstrukcji siatki jest to, że do modelu można łatwo wprowadzić wiele elektrod. Na przykład druga kopia geometrii elektrody może być umieszczona w odległości kilku milimetrów i obie zostaną uwzględnione w MES. W ostatnich badaniach dwie elektrody zostały wszczepione w bliskiej odległości w celu leczenia drżenia stwardnienia rozsianego13 i były wykorzystywane w eksperymentach na naczelnych innych niż ludzie w celu zbadania skutecznych celów stymulacji21. Zaletą stosowania wielu elektrod jest zapewnienie lepszej kontroli pola elektrycznego generowanego w tkance na większym obszarze. Stymulacja za pomocą aktywnych styków na obu elektrodach może skierować prąd w kierunku obszaru docelowego i z dala od obszarów, które prowadziłyby do negatywnych skutków ubocznych. Dokładniejsza kontrola stymulacji na większym obszarze jest również przydatna do badania różnych miejsc stymulacji, gdy dokładna lokalizacja zamierzonego celu jest nieznana, jak ma to miejsce w przypadku wielu pojawiających się zastosowań terapeutycznych DBS. Jednak określenie parametrów w celu uzyskania stymulacji terapeutycznej jest trudniejsze niż w przypadku pojedynczej elektrody ze względu na zwiększenie i tak już dużej przestrzeni parametrów.
Przewidujemy, że to interaktywne narzędzie do modelowania może przynieść korzyści podczas planowania przedoperacyjnego implantacji DBS. Informacje zwrotne na temat zakresu stymulacji w tkance mózgowej mogą pozwolić chirurgom na modyfikację położenia elektrod w ich planie chirurgicznym w celu zapewnienia stymulacji terapeutycznej w regionie docelowym. Podczas gdy terapia DBS była głównym czynnikiem motywującym do opracowania tego narzędzia, techniki przedstawione w tym artykule można zastosować do dowolnego modelu pola bioelektrycznego MES o różnych paradygmatach stymulacji lub rejestracji. Terapie stymulacyjne, takie jak przezczaszkowa stymulacja prądem stałym w przypadku depresji22 lub stosowanie elektrod głębinowych w leczeniu padaczki23 , mają te same wyzwania co DBS w określaniu najlepszej lokalizacji stymulacji w celu osiągnięcia wyników terapeutycznych. Elektrokortykografia, technika rejestracji z układem elektrod na powierzchni mózgu w celu identyfikacji regionów początku napadu, stoi przed wyzwaniem określenia, gdzie umieścić elektrody, aby rejestrować z docelowych obszarów mózgu24. Wszystkie te zastosowania są zależne od położenia elektrody, jednocześnie mając do czynienia z niepewnością dotyczącą przepływu prądu przez tkankę mózgową. Techniki przedstawione w tym artykule zmniejszają obciążenie związane z generowaniem i używaniem modeli obliczeniowych, zapewniając jednocześnie istotne informacje zwrotne badaczom i klinicystom korzystającym z tych urządzeń, którzy nie są ekspertami w dziedzinie modelowania.

Głęboka stymulacja mózgu (DBS) jest uznaną metodą leczenia zaburzeń ruchowych, takich jak drżenie samoistne1 i choroba Parkinsona2. Terapia ta jest również badana jako potencjalna metoda leczenia rosnącej liczby zaburzeń, w tym urazowego uszkodzenia mózgu3, zespołu Tourette'a4 i depresji5. Systemy DBS wymagają chirurgicznego wszczepienia elektrody w celu dostarczenia stymulacji w zlokalizowanym obszarze mózgu w celu modulowania trwającej aktywności neuronalnej6. Zarówno lokalizacja elektrod, jak i parametry stymulacji mają wpływ na modulację obwodów neuronalnych, które zapewniają korzyści terapeutyczne. Niewielkie różnice w lokalizacji mogą wpływać na okno terapeutyczne, co może zwiększyć prawdopodobieństwo wystąpienia działań niepożądanych przed osiągnięciem korzyści terapeutycznych7,8,9. W praktyce często trudno jest przewidzieć wpływ, jaki stymulacja będzie miała na aktywność neuronalną; W związku z tym to okno korzyści terapeutycznych jest identyfikowane indywidualnie dla każdego pacjenta, ponieważ urządzenie stymulujące jest programowane przez lekarza8,9. Proces ten staje się coraz bardziej złożony w miarę pojawiania się nowych generacji urządzeń DBS. Na przykład, nowe projekty leadów są wprowadzane z większą liczbą kontaktów10,11,12, a w niektórych przypadkach wiele przewodów jest wszczepianych w bliskiej odległości od siebie13. W związku z tym istnieje potrzeba badania i przewidywania skutków DBS w dużej i rosnącej przestrzeni parametrów.
Modelowanie i analiza obliczeniowa mogą być używane do przewidywania fizjologicznych i klinicznych skutków DBS w zależności od pacjenta. Modele te wykorzystują modelowanie metodą elementów skończonych (MES) do budowy trójwymiarowych reprezentacji tkanki mózgowej i cech biofizycznych wszczepionej elektrody. Modele pola bioelektrycznego MES były z powodzeniem wykorzystywane do przewidywania efektów DBS14, ale do tej pory były one czasochłonne i kosztowne obliczeniowo. Istnieje zapotrzebowanie na nowatorskie techniki modelowania, aby nadążyć za rosnącą złożonością urządzeń DBS. Te specyficzne dla pacjenta modele powinny zapewniać wizualną informację zwrotną w czasie zbliżonym do rzeczywistego na temat skutków DBS w miarę zmiany lokalizacji odprowadzania lub parametrów stymulacji. Użytkownik w ciągu kilku sekund uzyskałby informację zwrotną na temat lokalizacji leada i ustawienia stymulacji, co umożliwiłoby ciągłe doskonalenie umieszczania leadów w ciągu kilku minut. Specyficzność pacjenta uzyskuje się poprzez uwzględnienie anatomii pacjenta, kształtu i wielkości jego mózgu podczas budowy MES i zastosowania biofizycznych właściwości mózgu, takich jak anizotropowe przewodnictwo tkankowe. Przewodność anizotropowa opisuje, w jaki sposób prąd będzie rozprzestrzeniał się w różnych obszarach mózgu i może być nieinwazyjnie mierzony dla całego mózgu, podobnie jak typowy obraz rezonansu magnetycznego (MRI).
Metody modelowania DBS, które nie wykorzystują informacji specyficznych dla pacjenta, mogą dostarczyć szybkich, ale mniej dokładnych prognoz efektów stymulacji, ze względu na uogólnione geometrie i wartości przewodnictwa dla tkanki mózgowej. W tym podejściu stosuje się pojedynczy MES dla wszystkich pacjentów, a przewidywaną aktywność neuronalną można obliczyć z wyprzedzeniem. Modele specyficzne dla pacjenta nie mogą być uogólnione i wstępnie obliczone, ponieważ dla każdej osoby tworzony jest nowy MES. Zbudowanie tych modeli wymaga więcej wysiłku, ale może być dokładniejsze. Kilka czynników ogranicza szybkość, z jaką te modele mogą być budowane i używane: 1) Modyfikacja parametrów na początku procesu budowy modelu, takich jak położenie elektrody, wymaga ręcznego wysiłku w celu zaktualizowania wszystkich kolejnych kroków; oraz 2) Etapy procesu modelowania nie są łatwo ze sobą zintegrowane, co wymaga przekazywania danych między wieloma pakietami oprogramowania. Często chcemy ocenić wiele różnych sytuacji, takich jak położenie elektrod, parametry stymulacji lub konstrukcje elektrod. Aby zapewnić znaczącą informację zwrotną na temat wpływu tych zmian na efekt terapeutyczny, jaki otrzyma pacjent, wyniki te powinny być zarówno dokładne, jak i szybko wygenerowane.
Naszym celem jest przedstawienie nowych technik budowania modeli specyficznych dla pacjenta, które wykorzystują szybkość uzyskaną w modelach uogólnionych i automatyzują wiele kroków lejka, aby stworzyć interaktywne środowisko modelowania, które dostarcza wizualne informacje zwrotne w czasie zbliżonym do rzeczywistego na temat efektów DBS. Interaktywna symulacja pozwala użytkownikowi testować przewidywania i szybko uzyskiwać wyniki bez skupiania się na szczegółach budowy modelu. Jest to korzystne, gdy istnieje duża przestrzeń parametrów do zbadania, a sposób, w jaki te parametry wpływają na symulację, jest niepewny. Opiszemy etapy procesu przetwarzania w celu wygenerowania interaktywnych, specyficznych dla pacjenta modeli MES na podstawie akwizycji obrazowania metodą rezonansu magnetycznego (MRI). Wykorzystanie narzędzi i technik opisanych w tym artykule zmniejszy koszt tworzenia modeli pola bioelektrycznego MES i zapewni sposób na udostępnienie tych modeli badaczom i klinicystom, którzy nie są ekspertami w dziedzinie modelowania.
Ten protokół opisuje, jak zbudować specyficzny dla pacjenta model elementów skończonych na podstawie uzyskanych objętości MRI, a następnie symulować pole elektryczne indukowane przez elektrodę DBS. Głównymi etapami tworzenia tych modeli są: 1) zbudowanie modelu elementów skończonych (MES) reprezentującego mózg pacjenta i wszczepioną elektrodę, 2) dodanie do MES właściwości biofizycznych mózgu i parametrów stymulacji urządzenia DBS oraz 3) rozwiązanie napięcia generowanego przez elektrodę w modelu. Do zbudowania specyficznego dla pacjenta modelu symulacji DBS potrzebne są dwie metody obrazowania. Rezonans magnetyczny T1 służy do konstruowania segmentacji powierzchni mózgu, komór i określonych jąder. Obrazowanie zależne dyfuzją (DWI), miara dyfuzyjności wody, służy do szacowania tensorów dyfuzji w tkance mózgowej15. Tensory dyfuzji są przekształcane w tensory przewodności, które określają ilościowo niejednorodne, anizotropowe właściwości biofizyczne tkanki na podstawie woksela po wokselu16. Rozkład napięcia w mózgu indukowany przez elektrodę oblicza się poprzez rozwiązanie równania Poissona, które poprzez zastosowanie MES upraszcza się do liniowego układu równań Ax = b, gdzie A jest macierzą sztywności, która reprezentuje przewodność i geometrię siatki, x jest rozwiązaniem napięciowym w każdym węźle siatki, a b jest modyfikowany w oparciu o warunki brzegowe i źródła prądu.

<table><tbody><tr><td>FreeSurfer</td><td>Athinoula A. Martinos Centrum Obrazowania</td><td></td> Biomedycznego <td><a href="https://surfer.nmr.mgh.harvard.edu/">https://surfer.nmr.mgh.harvard.edu/</a></td></tr><tr><td>3D Slicer</td><td>BWH, Uniwersytet</td><td></td> <td>Harvarda<a href="https://www.slicer.org/">https://www.slicer.org/</a></td></tr><tr><td>SCIRun&nbsp;</td><td>Centrum Integracyjnych Obliczeń</td><td></td> Biomedycznych Uniwersytetu Utah<td><a href="http://www.sci.utah.edu/cibc-software/scirun.html">http://www.sci.utah.edu/cibc-software/scirun.html</a></td></tr></tbody></table>

celowanie w szlaki włókien nerwowych w celu terapii głębokiej stymulacji mózgu przy użyciu interaktywnych, specyficznych dla pacjenta modeli

Targeting Neuronal Fiber Tracts for Deep Brain Stimulation Therapy Using Interactive, Patient-Specific Models

Głęboka stymulacja mózgu (DBS), która polega na wprowadzeniu elektrody w celu dostarczenia stymulacji do zlokalizowanego obszaru mózgu, jest uznaną terapią zaburzeń ruchowych i jest stosowana w coraz większej liczbie zaburzeń. Modelowanie obliczeniowe zostało z powodzeniem wykorzystane do przewidywania skutków klinicznych DBS; jednak istnieje zapotrzebowanie na nowatorskie techniki modelowania, aby nadążyć za rosnącą złożonością urządzeń DBS. Modele te muszą również szybko i dokładnie generować prognozy. Celem tego projektu jest opracowanie procesu przetwarzania obrazu w celu włączenia obrazowania strukturalnego rezonansu magnetycznego (MRI) i obrazowania zależnego dyfuzją (DWI) do interaktywnego, specyficznego dla pacjenta modelu symulującego efekty DBS. Wirtualny przewód DBS można umieścić w modelu pacjenta, wraz z aktywnymi kontaktami i ustawieniami stymulacji, w których zmiany pozycji lub orientacji przewodu generują nową siatkę elementów skończonych i rozwiązanie problemu pola bioelektrycznego w czasie zbliżonym do rzeczywistego, w przedziale czasowym wynoszącym około 10 sekund. System ten umożliwia również symulację wielu przewodów w bliskiej odległości, aby umożliwić sterowanie prądem za pomocą różnych anod i katod na różnych przewodach. Techniki przedstawione w tym artykule zmniejszają obciążenie związane z generowaniem i wykorzystywaniem modeli obliczeniowych, zapewniając jednocześnie znaczące informacje zwrotne na temat wpływu położenia elektrody, konstrukcji elektrod i konfiguracji stymulacji dla badaczy lub klinicystów, którzy mogą nie być ekspertami w dziedzinie modelowania.

Na koniec niniejszego protokołu dostarczone zostaną wszystkie niezbędne dane wejściowe do stworzenia modelu specyficznego dla pacjenta: powierzchnia mózgu, geometria elektrod oraz tensory przewodnictwa. Powinna zostać również stworzona sieć SCIRun, która integruje wszystkie dane wejściowe w celu zbudowania modelu elementów kończonych i symulacji indukowanego pola bioelektrycznego. Model symulacyjny umożliwia przemieszczanie elektrody w obrębie mózgu oraz modyfikację parametrów, takich jak aktywne kontakty lub amplituda stymulacji.
Rycina 1 przedstawia możliwość symulacji i porównywania przewidywań aktywacji włókien dla wielu pozycji elektrod. Generowanie traktów włókien z obrazowania dyfuzyjnego nie zostało przedstawione w niniejszym protokole, jednak analiza ta może być przeprowadzona przy użyciu dowolnej deterministycznej metody traktografii. Pozycja 2 przewiduje silniejszą aktywację docelowego wiązki włókien przy niższych amplitudach stymulacji we wszystkich kontaktach w porównaniu z pozycją 1. Analiza ta jest wykorzystywana na etapie planowania operacji DBS w celu wyznaczenia trajektorii elektrody, która skutecznie zastymuluje docelową drogę włókien.
Nową cechą tego systemu jest możliwość szybkiej zmiany geometrii elektrody w modelu oraz możliwość symulowania wielu elektrod w bliskim sąsiedztwie. Obie te metody są stosowane w celu uzyskania większej kontroli nad kształtem i kierunkiem pola elektrycznego wokół elektrody. Rysunek 2 przedstawia porównanie przewidywań aktywacji włókien pomiędzy osiowo-symetryczną a kierunkową elektrodą DBS podczas badania lokalizacji elektrod w celu stymulacji przebiegów włókien wzgórza środkowego i unikania obszarów sąsiednich, takich jak jądra czuciowe. Przy tej samej trajektorii elektrody widać, że elektroda kierunkowa pozwala na skierowanie aktywacji włókien w stronę obszaru docelowego, unikając jednoczesnej stymulacji niepożądanych włókien. Analiza ta może być wykorzystana do ustalenia, kiedy konieczna byłaby zmiana urządzenia wszczepianego pacjentowi w ramach terapii DBS, aby uniknąć stymulacji obszarów sąsiednich, która wywoływałaby niepożądane skutki uboczne. Podejście to jest przydatne przy celowaniu w jądro podwzgórzowe w celu leczenia objawów choroby Parkinsona, przy jednoczesnym unikaniu pobliskiej torebki wewnętrznej.
<img alt="Wykres aktywacji włókien nerwowych w funkcji amplitudy stymulacji; pozycje elektrod; badania z zakresu neuronauki." class="xfigimg" src="/files/ftp_upload/57292/57292fig1.jpg" data-altupdatedat="1788250015000" data-alt="Rycina 1"> 
Rysunek 1Porównanie położenia elektrod i przewidywań aktywacji włókien docelowych. Przewidywanie całkowitej aktywacji wiązki włókien dla amplitud stymulacji (od -0,5 V do -5,0 V) i wszystkich czterech kontaktów dla dwóch różnych pozycji elektrody. <a href="https://www.jove.com/files/ftp_upload/57292/57292fig1large.jpg" target="_blank">Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.</a>
<img alt="Schemat porównawczy aktywacji elektrod do głębokiej stymulacji mózgu z aktywacją włókien przy różnych napięciach." class="xfigimg" src="/files/ftp_upload/57292/57292fig2.jpg" data-altupdatedat="1788250015000" data-alt="Rysunek 2"> 
Rycina 2Porównanie przewidywań aktywacji pęczków włókien brzuszno-przyśrodzeniowych wzgórza dla osiowo-symetrycznych oraz kierunkowych elektrod DBS. (A) Wizualizacja cylindrycznej elektrody DBS Medtronic 3387 oraz kierunkowej elektrody DBS Sapiens wraz z obliczeniami funkcji aktywacji w docelowym pęczku włókien wzgórza środkowego dla pojedynczej amplitudy stymulacji monopolarnej. (B)(i) Przekrój poprzeczny elektrody DBS i pęczków włókien. (ii) Dwuwymiarowa projekcja pola przekroju poprzecznego z identyfikacją elektrody DBS, włókien docelowych oraz włókien, których należy unikać. (C) Rozprzestrzenianie się aktywacji w obrębie zarówno włókien docelowych, jak i włókien unikowych wraz ze wzrostem amplitudy stymulacji dla osiowo-symetrycznej oraz kierunkowej elektrody DBS. Włókna aktywowane przy danej amplitudzie stymulacji zaznaczono na czerwono, natomiast włókna nieaktywowane – na niebiesko. (D) Kompilacja wszystkich symulacji przedstawionych w panelu C, ukazująca przewidywaną amplitudę progu stymulacji we wszystkich obszarach pęczków włókien. <a href="https://www.jove.com/files/ftp_upload/57292/57292fig2large.jpg">Prosimy o kliknięcie tutaj, aby wyświetlić powiększoną wersję tej ryciny.</a>

Celem tego projektu jest opracowanie interaktywnego, specyficznego dla pacjenta procesu modelowania, aby symulować efekty głębokiej stymulacji mózgu w czasie zbliżonym do rzeczywistego i dostarczyć istotnych informacji zwrotnych na temat tego, jak te urządzenia wpływają na aktywność neuronalną w mózgu.

Celowanie w szlaki włókien nerwowych w celu terapii głębokiej stymulacji mózgu przy użyciu interaktywnych, specyficznych dla pacjenta modeli

Deep brain stimulation (DBS), which involves insertion of an electrode to deliver stimulation to a localized brain region, is an established therapy for ...

Ogólnym celem tego projektu jest opracowanie interaktywnego procesu modelowania specyficznego dla pacjenta i stymulowanie wpływu głębokiej stymulacji mózgu na szlaki włókien nerwowych w czasie zbliżonym do rzeczywistego. Głęboka stymulacja mózgu to interwencja terapeutyczna, która została zatwierdzona przez FDA i jest stosowana od ponad dekady w zaburzeniach ruchowych, takich jak choroba Parkinsona, i została zbadana pod kątem wielu innych schorzeń neurologicznych i psychiatrycznych. Modelowanie obliczeniowe głębokiej stymulacji mózgu jest narzędziem, które jest bardzo przydatne, a jego wykorzystanie jest coraz częściej wykorzystywane w celu lepszego zrozumienia i przewidywania efektów stymulacji u poszczególnych osób.Jednym z uporczywych problemów związanych z modelowaniem obliczeniowym jest ilość czasu potrzebnego na konstruowanie i uruchamianie eksperymentów z modelami obliczeniowymi. Celem tego protokołu jest wyjaśnienie sposobów, w jakie modelowanie obliczeniowe może być wykonywane w czasie zbliżonym do rzeczywistego i w znacznie bardziej interaktywny sposób, a osiąga się to dzięki kilku innowacjom. Główną z nich jest możliwość interaktywnego tworzenia siatki całych głów poszczególnych obiektów.Protokół przedstawiony w tym filmie opisuje, jak zbudować specyficzny dla pacjenta model elementów skończonych na podstawie uzyskanych objętości MRI, a następnie symulować pole elektryczne indukowane przez wszczepioną elektrodę do głębokiej stymulacji mózgu. Główne kroki w tworzeniu tych modeli to, po pierwsze, zbudowanie modelu elementów skończonych reprezentującego mózg pacjenta i wszczepioną elektrodę, po drugie, dodanie właściwości biofizycznych mózgu i parametrów stymulacji urządzenia do głębokiej stymulacji mózgu do siatki elementów skończonych, a po trzecie, rozwiązanie napięcia generowanego przez elektrodę w modelu. Elektroda modelowana w tym projekcie to przewód DBS firmy Medtronic 3387.Jednak przy tym podejściu można użyć dowolnej geometrii elektrody. 3387 to cylindryczny przewód z czterema stykami oddzielonymi materiałem nieprzewodzącym. Reprezentacje powierzchni dla każdego z tych komponentów są tworzone dla procesu tworzenia siatki.Zamknięte powierzchnie są tworzone dla czterech styków, pięciu segmentów wału, warstwy enkapsulacji i dwóch otaczających skrzynek, które służą do kontrolowania gęstości siatki wokół elektrody. Ważne jest, aby stworzyć siatkę powierzchniową styków o dużej gęstości, aby zachować cylindryczną geometrię i dokładnie rozwiązać problem pola bioelektrycznego. Rozkład napięcia w mózgu indukowany przez elektrodę jest obliczany przez rozwiązanie równania Poissona, które dzięki zastosowaniu modelowania elementów skończonych upraszcza się do liniowego układu równań, Ax równa się B, gdzie A jest macierzą sztywności, która reprezentuje przewodność i geometrię siatki, x jest rozwiązaniem napięciowym w każdym węźle siatki, a B jest modyfikowany w oparciu o warunki brzegowe i źródła prądu.Warstwa hermetyzacji, segmenty wału i cztery styki również muszą mieć przypisane wartości przewodności w oparciu o ich właściwości materiałowe pokazane w tabeli na ekranie. Do obliczenia pola bioelektrycznego stosuje się stałe warunki brzegowe. Powierzchnie punktowe są umieszczone w środku każdego z czterech styków DBS, a stymulacja monopolarna lub bipolarna może być stosowana za pomocą dowolnej kombinacji.Zewnętrzna powierzchnia mózgu jest ustawiona na zero woltów, aby symulować odległą elektrodę powrotną. Warunki te reprezentują znane napięcie w modelu elementów skończonych. Znamy napięcie na każdym z czterech styków oraz na elektrodzie powrotnej, która jest wszczepialnym generatorem impulsów.Liniowy układ równań jest rozwiązywany przy użyciu metody sprzężonego gradientu w celu obliczenia rozwiązania napięcia w modelu dla zestawu warunków brzegowych. SCIRun to zaawansowane środowisko rozwiązywania problemów do modelowania, symulacji i wizualizacji kilku klas problemów naukowych. Wizualna składnia programowania pozwala użytkownikom łączyć ze sobą modele oprogramowania w celu przetwarzania i wizualizacji danych naukowych.W przypadku tego protokołu SCIRun służy do tworzenia siatki elementów skończonych mózgu i głębokiej stymulacji mózgu oraz rozwiązywania problemu pola bioelektrycznego w interaktywnym środowisku. Jednym z kluczowych aspektów SCIRun jest możliwość interaktywnej wizualizacji danych w dowolnym momencie potoku, co zademonstrujemy. Pierwszym krokiem w budowie modelu elementów skończonych jest określenie anatomii strukturalnej mózgu.Aby wykonać ten krok, zdecydowaliśmy się użyć FreeSurfer, pakietu oprogramowania typu open source do przetwarzania i analizowania objętości MRI dla ludzkich mózgów. Aby rozpocząć, rezonans magnetyczny pacjenta T1 jest przekazywany do FreeSurfer jako stos DICOM za pomocą następującego polecenia. Następnie następujące polecenie służy do rozpoczęcia automatycznego przetwarzania i segmentacji.FreeSurfer tworzy nowy obraz objętości tkanki mózgowej z usuniętą czaszką. Progowanie intensywności i izosurfacing są wykonywane na tej nowej objętości w celu uzyskania powierzchniowej reprezentacji mózgu, która jest następnie wykorzystywana do zdefiniowania zewnętrznej granicy modelu elementów skończonych. Rurociąg FreeSurfer segmentuje również istotę białą, istotę szarą, jądra podkorowe i komory.Każdy region mózgu ma przypisaną unikalną etykietę liczby całkowitej. Kolejnym etapem przetwarzania wstępnego jest oszacowanie tensorów dyfuzji na podstawie danych obrazowania ważonych dyfuzją i jednoczesne zarejestrowanie tej objętości za pomocą TI MRI. DWI to pomiar dyfuzyjności wody przez tkankę mózgową w wielu kierunkach.Na tym etapie zdecydowaliśmy się skorzystać z 3DSlicer, bezpłatnej platformy oprogramowania open source do przetwarzania obrazów medycznych i wizualizacji 3D danych obrazowych. Najpierw przekonwertuj serię DICOM na pojedynczy wolumin DWI za pomocą modułu konwertera DWI w Slicerze. W celu wyeliminowania tła i szumów z rekonstrukcji tensora należy utworzyć maskę z modułem maskowania objętości ważonej dyfuzją.Pozwoli to na obliczanie tensorów tylko w tkance mózgowej. Następnym krokiem jest użycie modułu DWI do DTI Estimation do wygenerowania tensorów dyfuzji. Ustaw algorytm estymacji tensora na najmniejsze kwadraty i zaznacz opcję Przesunięcie ujemnych wartości własnych.Spowoduje to usunięcie tensorów o ujemnych wartościach własnych związanych z szumem lub błędem akwizycji. Geometrie powierzchni 3D dla modelu elektrody można zbudować przy użyciu dowolnego oprogramowania do modelowania 3D, które nie zostanie zademonstrowane w tym filmie. W tym miejscu utworzyliśmy siatkę o zamkniętej powierzchni dla każdego ze styków elektrody, segmentów wału, warstwy enkapsulacyjnej i dwóch otaczających pól, które służą do kontrolowania gęstości siatki.Powierzchnia mózgu utworzona w części pierwszej i geometria elektrody są łączone w celu stworzenia modelu mózgu metodą elementów skończonych z wszczepioną elektrodą za pomocą narzędzia do tworzenia siatki TetGen. Przeprowadź powierzchnię mózgu i geometrię elektrody do modułu InterfaceWithTetGen. Zastosuj ograniczenia objętości do generowania elementów czworościennych w każdym regionie modelu, jak pokazano w tabeli na ekranie.Ta technika tworzy siatkę o adaptacyjnym rozmiarze i większej gęstości oczek wokół granicy faz tkanki elektrody, gdzie gradient pola elektrycznego jest wysoki. Rezultatem jest czworościenna siatka elementów skończonych modelu elektrody osadzona w tkance mózgowej. SCIRun umożliwia przesuwanie obiektów za pomocą modułów transformacji geometrycznych i interaktywnych widżetów w oknie przeglądarki.Aby zmienić kąt geometrii elektrody, użyj modułu CreateGeometricTransform w połączeniu z modułem TransformMesh. Ustaw oś obrotu X na jedną i ustaw stopień obrotu dla theta, aby kontrolować kąt. Podłącz moduł GenerateSinglePointProbe bezpośrednio do ViewScene, aby utworzyć widżet sfery, a także do innego modułu TransformMesh, aby kontrolować położenie końcówki elektrody.W scenie widoku kliknij z wciśniętym Shift i przeciągnij widżet, aby przesunąć położenie końcówki elektrody. Transformacje te są stosowane do geometrii elektrody przed przekazaniem jej do TetGen, więc zaktualizowana lokalizacja i orientacja elektrody jest wbudowana w model elementów skończonych. Kolejnym krokiem jest przypisanie wartości przewodności do modelu elementów skończonych.Załaduj dane DTI, które zostały przetworzone w 3DSlicer i ustaw wartości przewodności zarówno dla modelu elektrody, jak i tkanki mózgowej na czworościennej siatce utworzonej w poprzedniej sekcji. Po zastosowaniu wartości przewodności do siatki, skieruj pole do BuildFEMatrix. Jest to pierwszy krok w tworzeniu układu równań do obliczania pola bioelektrycznego.W przypadku stymulacji sterowanej napięciem ustawimy teraz znane warunki brzegowe napięcia na siatce elementów skończonych. Ustaw powierzchnię mózgu używaną jako zewnętrzna granica naszej siatki na zero woltów, aby zasymulować odległą elektrodę powrotną. InsertVoltageSource jest następnie używany do zastosowania tego warunku brzegowego do siatki.W przypadku stymulacji monopolarnej przypisz wartość napięcia do punktu pośrodku dolnego styku. Ponownie użyj InsertVoltageSource, aby zastosować ten warunek brzegowy do siatki. Połącz dane wyjściowe BuildFEMatrix i InsertVoltageSource w SolveLinearSystem.Użyj metody rozwiązywania sprzężonego gradientu i wykonaj moduł, aby obliczyć roztwór napięcia. Umieść rozwiązanie układu równań z powrotem w węzłach siatki elementów skończonych za pomocą modułu SetFieldData. Roztwór napięcia można oglądać bezpośrednio na węzłach siatki, ale nie zapewnia znaczącej reprezentacji jego dystrybucji w tkance mózgowej.Lepszą wizualizacją jest utworzenie powierzchni izonapięciowych podobnych do warstwic na mapie, która dostarcza również informacji o orientacji i natężeniu pola elektrycznego. Przekaż rozwiązanie napięcia siatki elementów skończonych do ExtractSimpleIsosurface i ustaw ilość na 10. Spowoduje to utworzenie 10 równomiernie rozmieszczonych powierzchni izonapięciowych wokół styku.Dzięki tej technice możemy zwizualizować efekt dodania do modelu tensorów przewodnictwa anizotropowego oraz to, jak ten parametr wpływa na rozkład napięcia w tkance mózgowej. Integracja konstrukcji elementów skończonych, obliczeń pola bioelektrycznego i wizualizacji wyników w jednym środowisku umożliwia uzyskanie informacji zwrotnej w czasie zbliżonym do rzeczywistego, w jaki sposób zmiany parametrów, takich jak położenie elektrody i amplituda stymulacji, wpływają na to, jak pole bioelektryczne oddziałuje z otaczającą tkanką. Po zmodyfikowaniu położenia elektrody i ponownym obliczeniu pola bioelektrycznego, roztwór jest mapowany na docelowe wiązki włókien.Funkcja aktywacji jest obliczana w celu przewidywania aktywacji wiązki włókien dla danej amplitudy stymulacji. Zbieramy prognozy aktywacji wiązki włókien we wszystkich kontaktach i wielokrotnych amplitudach stymulacji i szybko określamy, która konfiguracja stymulacji byłaby najbardziej skuteczna w aktywacji docelowych włókien dla danej pozycji elektrody. Następnie porównujemy te profile aktywacji dla wielu pozycji elektrod na etapach planowania operacji DBS, aby zdefiniować trajektorię elektrody, która będzie skutecznie stymulować docelową ścieżkę włókien.Tutaj widzimy, że pozycja druga przewiduje bardziej solidną aktywację wiązki włókien docelowych, ma niższe amplitudy stymulacji we wszystkich kontaktach w porównaniu z pozycją pierwszą. Takie podejście można wykorzystać do modelowania znacznie bardziej złożonych scenariuszy. Nowatorską cechą tego systemu jest możliwość szybkiej zmiany geometrii elektrody w modelu oraz możliwość symulacji wielu przewodów w bliskiej odległości od siebie.Tutaj badamy lokalizacje prowadzące, aby stymulować szlaki włókien wzgórzowych biorących udział w regulacji pobudzenia jako potencjalny cel dla pacjentów z urazowym uszkodzeniem mózgu, unikając pobliskich obszarów, takich jak jądra czuciowe. W przypadku statycznego położenia elektrody przeprowadziliśmy symulację skutków zmiany na kierunkowy model elektrody, aby uzyskać bardziej skoncentrowaną kontrolę nad miejscami aktywacji włókien i uniknąć stymulacji innych wiązek włókien. Patrząc na pole przekroju poprzecznego wokół elektrody, widzimy, że elektroda kierunkowa jest w stanie skutecznie stymulować docelowe włókna przy minimalnym rozprzestrzenianiu się aktywacji na włókna, których chcemy uniknąć.Takie podejście byłoby również przydatne w przypadku celowania w jądro podwzgórzowe w celu leczenia objawów choroby Parkinsona, unikając jednocześnie pobliskiej torebki wewnętrznej. Tutaj symulujemy rozmieszczenie trzech niezależnych przewodów DBS i przewidujemy aktywację podobnych ścieżek centralnych włókien wzgórzowych do planowania chirurgicznego w eksperymentach na naczelnych innych niż ludzie. Zaletą stosowania wielu elektrod jest zapewnienie kontroli nad stymulacją na większym obszarze.Zarówno elektrody kierunkowe, jak i wieloodprowadzeniowe zwiększają trudność określenia parametrów stymulacji elektrod. Jednak wykorzystanie podejścia opisanego w tym filmie umożliwia szybszą i bardziej ukierunkowaną eksplorację złożonej przestrzeni parametrów w celu zidentyfikowania ustawień stymulacji terapeutycznej. Procedura modelowania specyficznego dla pacjenta, przedstawiająca etapy od obrazowania pacjenta do interaktywnych modeli, którą pokazaliśmy dzisiaj, pomoże wypełnić lukę między inżynierami a klinicystami.Opracowanie nowych urządzeń DBS z więcej niż czterema stykami, w tym wprowadzenie wielu elektrod, zwiększa złożoność i tak już trudnego problemu. Możliwość szybkiej wizualizacji efektów różnych konstrukcji elektrod w konfiguracjach stymulacji może dostarczyć cennych informacji zwrotnych, aby lepiej zrozumieć, jak te urządzenia działają w mózgu.

targeting-neuronal-fiber-tracts-for-deep-brain-stimulation-therapy

Research

JoVE Journal

Bioengineering

Bioinżynieria

Mechanizm leżący u podstaw terapeutycznych efektów operacji Głębokiej Stymulacji Mózgu (DBS) wymaga zbadania. Metody przedstawione w tym manuskrypcie opisują eksperymentalne podejście do badania zdarzeń komórkowych wywoływanych przez DBS poprzez analizę profilu ekspresji genów kandydujących, które mogą ułatwić neurogenezę po operacji DBS.

Analiza zmian ekspresji genów w hipokampie szczura po głębokiej stymulacji mózgu przedniego jądra wzgórza

Tutaj przedstawiamy procedurę powtarzalnej przezczaszkowej stymulacji magnetycznej pod kontrolą MRI do grzbietowo-przyśrodkowej kory przedczołowej jako eksperymentalnego leczenia dużego zaburzenia depresyjnego.

dmPFC-rTMS pod kontrolą MRI jako metoda leczenia opornego na leczenie dużego zaburzenia depresyjnego

Protokół pokazuje prototyp domowej platformy do zbierania danych multimodalnych, która wspiera badania optymalizujące adaptacyjną głęboką stymulację mózgu (aDBS) dla osób z neurologicznymi zaburzeniami ruchu. Przedstawiamy również kluczowe wnioski z wdrażania platformy przez ponad rok w domu osoby z chorobą Parkinsona.

Przeniesienie kliniki do domu: multimodalny ekosystem gromadzenia danych w domu wspierający adaptacyjną głęboką stymulację mózgu

Ten protokół opisuje metodę umieszczania elektrod sterowanych neuronami do ogniskowej, przezczaszkowej stymulacji prądem stałym (tDCS) podawanej podczas funkcjonalnego rezonansu magnetycznego (fMRI).

Neuronavigated Zogniskowana przezczaszkowa stymulacja prądem stałym podawana podczas funkcjonalnego rezonansu magnetycznego

Nowatorski protokół spersonalizowanej przezczaszkowej stymulacji magnetycznej w połączeniu z mapowaniem sieci korowych za pomocą elektroencefalografii (TMS–EEG), oparty na wielu modalnościach obrazowania metodą rezonansu magnetycznego (MRI).

Neuroobrazowanie TMS-EEG do mapowania sieci korowej w czasie rzeczywistym

Ten artykuł opisuje, jak lokalizować cele specyficzne dla funkcji dla powtarzających się przezczaszkowych interwencji lub zabiegów, gdy sprzęt nawigacyjny jest niedostępny.

Lokalizacja celów specyficznych dla funkcji przezczaszkowej stymulacji magnetycznej w przypadku braku sprzętu nawigacyjnego

Tutaj opisujemy standaryzowany protokół mapowania motorycznego z wykorzystaniem nTMS połączony z rekonstrukcją drogi korowo-rdzeniowej (CST) opartą na obrazowaniu tensorowym dyfuzji (DTI). Protokół jest powtarzalny, klinicznie wykonalny i łatwo integrowalny z rutynowym procesem klinicznym, zapewniając solidne i wartościowe ramy do oceny szlaków ruchowych, badań nad neuroplastycznością oraz planowania rehabilitacji.

Standaryzowany protokół do funkcjonalnego mapowania motorycznego z wykorzystaniem nawigowanej przezczaszkowej stymulacji magnetycznej

Protokół ten opisuje sposób pracy wykonywanej przez robotycznie sterowaną terapię magnetyczną stymulacji przezczaszkowej w leczeniu depresji opornej na leczenie. Kliniczne wdrożenie z wykorzystaniem tego protokołu jest wykonalne, a wyniki otwartej etykiety potwierdzają skuteczność tego podejścia w realnych warunkach.

Robotycznie dostarczana spersonalizowana terapia stymulacją magnetyczną przezczaszki w leczeniu depresji opornej na leczenie

Tu opisujemy protokół dotyczący ex vivo nagrań z wysoko gęstej matrycy mikroelektrody z ostrych kroi móżdżku u modelu chomika z dystonią, który otrzymał in vivo ciągłą stymulację głęboką mózgu w zakresie 11 dni.

Badanie mechanizmów działania stymulacji głębokiej mózgu obszaru ciemieniowego za pomocą wysokoczułowych mikroelektrodowych tablic

W wielu sytuacjach biologicznych i klinicznych korzystne jest badanie procesów komórkowych w ich naturalnym mikrośrodowisku. W niniejszym materiale opisujemy montaż i zastosowanie niedrogiego mikroskopu światłowodowego, który umożliwia obrazowanie w czasie rzeczywistym w hodowlach komórkowych, badaniach na zwierzętach oraz w badaniach klinicznych z udziałem pacjentów.

Wysokorozdzielcza mikroendoskopia światłowodowa do obrazowania komórkowego in situ

Poznaj interaktywne modelowanie głębokiej stymulacji mózgu przy użyciu danych MRI i oprogramowania SCIRun. Ta metoda umożliwia symulację w czasie rzeczywistym wpływu stymulacji na szlaki nerwowe u pacjentów z chorobą Parkinsona.

Celowanie w szlaki włókien nerwowych w celu terapii głębokiej stymulacji mózgu przy użyciu interaktywnych, specyficznych dla pacjenta modeli (Wideo) | JoVE Journal

Celowanie w szlaki włókien nerwowych w celu terapii głębokiej stymulacji mózgu przy użyciu interaktywnych, specyficznych dla pacjenta modeli

Przeglądaj więcej filmów