$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Podejście obliczeniowe staje się kluczowe dla postępu w zrozumieniu funkcji serca w zdrowiu i chorobie. Najnowocześniejsze modele elektrofizjologii i elektromechaniki całego serca są obecnie wykorzystywane do badania szerokiego zakresu zjawisk, takich jak prawidłowa propagacja komór, arytmia, defibrylacja, sprzężenie elektromechaniczne i resynchronizacja serca1. Jednak, aby podejście obliczeniowe mogło mieć bezpośrednie zastosowanie w środowisku klinicznym, konieczne jest, aby modele były specyficzne dla pacjenta, tj. modele muszą być oparte na określonej architekturze i właściwościach elektrofizjologicznych lub elektromechanicznych chorego serca pacjenta. Symulacja z wykorzystaniem takich modeli pomoże lekarzom w podejmowaniu wysoce spersonalizowanych decyzji dotyczących interwencji elektrofizjologicznych, a także profilaktyki, co znacznie poprawi opiekę zdrowotną nad pacjentami2-4 .
Tworzenie realistycznych modeli serca wymaga zrozumienia geometrii i struktury włókien serca pacjenta. Orientacje włókien określają kierunki propagacji elektrycznej i rozkładu odkształceń w sercu, a zatem ich pozyskiwanie jest niezbędne do modelowania serca5, 6. Dzięki niedawnym postępom w obrazowaniu medycznym możliwe jest obecnie uzyskanie geometrii serca pacjenta, w tym przebudowy strukturalnej, takiej jak zawał, in vivo z wysoką rozdzielczością przy użyciu technologii rezonansu magnetycznego (MRI) i tomografii komputerowej (CT). Nie ma jednak praktycznej metody pozyskiwania struktury włókien serca pacjenta in vivo. Tensor dyfuzji (DT) MRI7, 8, jedyna technika uzyskiwania orientacji włókien w nienaruszonym sercu, nie jest powszechnie dostępna in vivo ze względu na pewne ograniczenia 9. Krótki opis wcześniejszych prób przełożenia DTMRI na warunki kliniczne można znaleźć w innym miejscu 2. Chociaż metodologie, takie jak przypisywanie orientacji włókien w oparciu o reguły, stanowią alternatywę dla DTMRI, metodologie te mają pewne poważne ograniczenia 2, 10. W związku z tym trudności w uzyskaniu struktury włókien serca in vivo utrudniają obecnie zastosowanie elektrofizjologicznych i elektromechanicznych symulacji serca w warunkach klinicznych. Celem tych badań było bezpośrednie odniesienie się do tej potrzeby.
Postawiliśmy hipotezę, że orientacje włókien komorowych serca można dokładnie przewidzieć, biorąc pod uwagę geometrię serca i atlas, gdzie atlas to serce, którego geometria i orientacje włókien są dostępne. W związku z tym wykorzystaliśmy najnowocześniejsze techniki do opracowania metodologii szacowania orientacji włókien sercowych in vivo i przetestowaliśmy hipotezę w normalnych i niewydolnych komorach psów2. Główną ideą naszej metodologii szacowania włókien jest wykorzystanie podobieństw w orientacji włókien w odniesieniu do geometrii między różnymi sercami w celu przybliżenia struktury włókien (docelowego) serca, dla którego dostępne są tylko informacje o geometrii. Sercem naszej metodologii szacowania jest rejestracja geometrii atlasu z geometrią docelową za pomocą dyfeomorficznego mapowania metrycznego dużych deformacji (LDDMM) 11 oraz morfing orientacji włókien atlasu przy użyciu zachowania głównych składowych (PPD) 2, 12. Właściwości dyfeomorficzne LDDMM gwarantują, że atlas nie "zagina się" podczas deformacji, zachowując w ten sposób integralność struktur anatomicznych. Rysunek 1 ilustruje proces przetwarzania w ramach naszej metodologii. W sekcji §1 tekstu protokołu opisano różne elementy rurociągu, pokazując, w jaki sposób można przeprowadzić oszacowanie dla przykładowego pacjenta. Numery wewnątrz niektórych bloków na rysunku 1 odnoszą się do odpowiednich podsekcji w sekcji §1 tekstu protokołu.
Oceniliśmy wydajność proponowanej metodologii, kwantyfikując błąd oszacowania i mierząc wpływ tego błędu na symulacje elektrofizjologii serca, symulując obliczeniowo lokalne mapy aktywacji elektrycznej, jak również pseudo-elektrokardiogramy (pseudo-EKG). Ze względu na niedostępność ludzkich serc, ocenę działania przeprowadzono przy użyciu serc psów dostępnych z poprzednich badań 13-15 . Błąd oszacowania obliczono za pomocą kątów nachylenia 16, zgodnie z tradycją histologii, w której pomiary kątowe wykonuje się na odcinkach tkanki, które są cięte równolegle do powierzchni nasierdzia. Ponieważ kąt między kierunkiem włókna a płaszczyzną styczną nasierdzia jest na ogół mały 17, 18, utrata informacji przy opisywaniu kierunku włókna wyłącznie za pomocą jego kąta nachylenia jest nieznaczna. Na potrzeby symulacji obliczeniowych zbudowano modele oparte na obrazach, jak opisano wcześniej 19, 20, a tkanka serca w modelach została przedstawiona w oparciu o ustalone techniki matematyczne i dane eksperymentalne 21-25 . Rytm zatokowy symulowano przez replikację aktywacji pochodzącej z sieci Purkinjego26 i częstoskurcz komorowy za pomocą protokołu stymulacji S1-S2 27. Pseudo-EKG obliczono 28 i porównano przy użyciu wskaźnika średniego odchylenia bezwzględnego (MAD) 29.