$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Reprezentatywne wyniki zostały uzyskane zgodnie z tym protokołem w celu zbadania powiązań między białkami mitochondrialnymi (Tabela 2) a ośmioma kategoriami chorób sercowo-naczyniowych (Tabela 3). W tych kategoriach znaleźliśmy 363 567 publikacji opublikowanych od 2012 r. do października 2022 r. (362 878 skategoryzowanych według metadanych MeSH, 6 923 skategoryzowanych według imputacji etykiet). Wszystkie publikacje miały tytuły, 276 524 miały abstrakty, a 51 065 miało dostępne pełne teksty. Ogólnie rzecz biorąc, 584 z 1687 badanych białek mitochondrialnych zostało zidentyfikowanych w publikacjach, podczas gdy 3284 z 8026 badanych funkcjonalnie powiązanych białek zostało zidentyfikowanych. W sumie zidentyfikowano 14 unikalnych białek ze znaczącymi wynikami we wszystkich kategoriach choroby, z progiem z-score wynoszącym 3,0 (Ryc. 5). Analiza szlaku Reactome tych białek ujawniła 12 szlaków istotnych dla wszystkich chorób (Ryc. 6). Wszystkie białka, szlaki, choroby i wyniki zostały zintegrowane w grafie wiedzy (Tabela 4). Ten wykres wiedzy został wykorzystany do przewidzenia 12 688 nowych powiązań białkowych z chorobami i przefiltrowany z wynikiem prawdopodobieństwa 0,90, aby uzyskać 1 583 prognozy o wysokim poziomie pewności. Wyróżniony przykład dwóch powiązań białko-choroba jest pokazany na Rysunek 7, zilustrowany w kontekście innych istotnych jednostek biologicznych funkcjonalnie związanych z białkami. Metryki oceny modelu przedstawiono w tabeli 5.

Rysunek 1: Dynamiczny widok przepływu pracy. Na tym rysunku przedstawiono cztery główne kroki w tym przepływie pracy. Po pierwsze, odpowiednie białka są selekcjonowane na podstawie terminów GO dostarczonych przez użytkownika (np. składniki komórkowe), a kategorie chorób są przygotowywane na podstawie dostarczonych przez użytkownika identyfikatorów choroby MeSH. Po drugie, powiązania między białkami a chorobami są obliczane na etapie eksploracji tekstu. Publikacje w określonym zakresie dat są pobierane i indeksowane. Publikacje badające chorobę są identyfikowane (za pomocą etykiet MeSH i opcjonalnie za pomocą etykiet imputowanych), a ich pełne teksty są pobierane i indeksowane. Nazwy białek są wyszukiwane w publikacjach i wykorzystywane do obliczania wyników asocjacji białko-choroba. Następnie, po eksploracji tekstu, wyniki te pomagają zidentyfikować najważniejsze powiązania białek i szlaków. Na koniec konstruowany jest wykres wiedzy obejmujący te białka, choroby i ich relacje w bazie wiedzy biomedycznej. Nowe powiązania białko-choroba są przewidywane na podstawie skonstruowanego wykresu wiedzy. W tych krokach wykorzystywane są najnowsze dostępne dane z baz wiedzy biomedycznej i PubMed. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Techniczna architektura przepływu pracy. Szczegóły techniczne tego przepływu pracy przedstawiono na poniższym rysunku. Użytkownik podaje numery drzew MeSH kategorii chorób i terminów GO. Dokumenty tekstowe są pobierane z PubMed, dokumenty związane z chorobą są identyfikowane na podstawie dostarczonych etykiet MeSH, a dokumenty bez etykiet MeSH wskazujących temat otrzymują przypisane etykiety kategorii. Białka związane z podanymi terminami GO są pozyskiwane. Ten zestaw białek został rozszerzony o białka, które są funkcjonalnie powiązane poprzez interakcje białko-białko, wspólne szlaki biologiczne i zależność od czynników transkrypcyjnych. Białka te są wyszukiwane w dokumentach istotnych dla choroby i oceniane za pomocą CaseOLAP. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Przykład przetworzonego dokumentu. Poniżej przedstawiono przykład przetworzonego, indeksowanego dokumentu tekstowego. W odpowiedniej kolejności odpowiednie pola wskazują nazwę indeksu (_index, _type), identyfikator PubMed (_id, pmid), podsekcje dokumentu (tytuł, streszczenie, full_text, wprowadzenie, metody, wyniki, dyskusja) oraz inne metadane (rok, MeSH, lokalizacja, czasopismo). Tylko w celu wyświetlenia podsekcje dokumentu są obcinane za pomocą wielokropka. Pole MeSH zawiera tematy dokumentów, które czasami mogą być dostarczane przez nasz krok imputacji etykiety. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Schemat grafu wiedzy i zasoby biomedyczne. Ten rysunek przedstawia schemat grafu wiedzy. Każdy węzeł i krawędź reprezentują odpowiednio węzeł lub typ krawędzi. Granice między chorobami sercowo-naczyniowymi (CVD) a białkami są ważone za pomocą wyników CaseOLAP. Krawędzie interakcji białko-białko (PPI) są ważone przez wyniki ufności STRING. Krawędzie zależności od czynnika transkrypcyjnego (TFD) pochodzące z GRNdb/GTEx, krawędzie drzewa chorób pochodzące z MeSH i krawędzie szlaku pochodzące z reaktomu są nieważone. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Najważniejsze powiązania białkowo-chorobowe. Rysunek ten przedstawia białka mitochondrialne istotne dla każdej kategorii choroby. Transformację Z-score zastosowano do wyników CaseOLAP w każdej kategorii w celu zidentyfikowania istotnych białek przy użyciu progu 3,0. (Góra) Liczba białek mitochondrialnych istotnych dla każdej choroby: Te wykresy skrzypcowe przedstawiają rozkład z-score dla białek w każdej kategorii choroby. Całkowita liczba białek istotnych dla każdej kategorii choroby jest przedstawiona powyżej każdego wykresu skrzypiec. W sumie 14 unikalnych białek zostało zidentyfikowanych jako istotne we wszystkich chorobach, a niektóre białka były istotne dla wielu chorób. (Na dole) Białka z najwyższą oceną: Mapa cieplna wyświetla 10 najlepszych białek, które uzyskały najwyższe średnie wyniki z we wszystkich chorobach. Wartości ślepej próby nie reprezentują uzyskanego wyniku między białkiem a chorobą. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Najważniejsze powiązania między szlakiem a chorobą. Rysunek ten ilustruje najważniejsze szlaki biologiczne związane z badanymi kategoriami chorób, określone za pomocą analizy szlaków reaktomowych. Wszystkie analizy ścieżek zostały przefiltrowane za pomocą p < 0,05. Wartości mapy cieplnej reprezentują średni wskaźnik z wszystkich białek w szlaku. (Góra) Szlaki zachowane we wszystkich chorobach: Ogólnie rzecz biorąc, zidentyfikowano 14 białek mających znaczenie dla wszystkich kategorii chorób i ujawniono 12 konserwatywnych szlaków we wszystkich kategoriach chorób. Dendrogram został skonstruowany w oparciu o hierarchiczną strukturę szlaku, aby połączyć szlaki o podobnych funkcjach biologicznych. Wysokość dendrogramu reprezentuje względną głębokość w hierarchii szlaku; Szerokie funkcje biologiczne mają dłuższe kończyny, a bardziej specyficzne ścieżki mają krótsze kończyny. (Na dole) Szlaki odrębne dla danej kategorii choroby: Analizę szlaków przeprowadzono przy użyciu białek, które osiągnęły znaczący wskaźnik Z w każdej chorobie. Trzy najważniejsze szlaki z najniższymi wartościami p związanymi z każdą chorobą są pokazane i oznaczone gwiazdkami. Szlaki mogą znajdować się w pierwszej trójce w wielu chorobach. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 7: Zastosowanie głębokiego uczenia do uzupełniania grafów wiedzy. Przykład zastosowania uczenia głębokiego do wykresu wiedzy specyficznego dla choroby przedstawiono na poniższym rysunku. Przewiduje się ukryte związki między białkami a chorobą, które są oznaczone kolorem niebieskim. Wyświetlane są obliczone prawdopodobieństwa dla obu przewidywań, z wartościami z zakresu od 0,0 do 1,0, przy czym 1,0 oznacza silną prognozę. Uwzględniono kilka białek o znanych interakcjach, reprezentujących interakcje białko-białko, zależność od czynnika transkrypcyjnego i wspólne szlaki biologiczne. W celu wizualizacji pokazano podgraf kilku węzłów związanych z wyróżnionym przykładem. Klucz: IHD = choroba niedokrwienna serca; R-HSA-1430728 = metabolizm; O14949 = podjednostka 8 kompleksu cytochromu b-c1; P17568 = dehydrogenaza NADH (ubichinon) 1 podjednostka beta podkompleksu 7; Q9NYF8 Czynnik transkrypcyjny 1 związany z Bcl-2, wynik: 7,24 x 10−7; P49821 = dehydrogenaza NADH (ubichinon) flawoproteina 1, mitochondrialna, wynik: 1,06 x 10−5; P31930 = podjednostka 1 kompleksu cytochromu b-c1, mitochondrialna, wynik: 4,98 x 10−5; P99999 = cytochrom c, wynik: 0,399. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Tabela 1: Etapy przepływu pracy i ograniczania szybkości. W poniższej tabeli przedstawiono przybliżone szacunki czasu obliczeniowego dla każdego etapu przepływu pracy. Opcje dołączania składników potoku spowodują zmianę całkowitego czasu wykonywania wymaganego do ukończenia analizy. Całkowity szacowany czas różni się w zależności od dostępnych zasobów obliczeniowych, w tym specyfikacji sprzętu i ustawień oprogramowania. Z grubsza szacuje się, że wykonanie protokołu na naszym serwerze obliczeniowym zajęło 36 godzin aktywnego czasu działania, z sześcioma rdzeniami, 32 GB pamięci RAM i 2 TB pamięci masowej, ale na innych urządzeniach może to być szybsze lub wolniejsze. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Automatyczne składanie białek składnika komórkowego. Poniższa tabela przedstawia liczbę białek związanych z danym komponentem komórkowym (tj. termin GO), białek funkcjonalnie z nimi związanych poprzez interakcje białko-białko (PPI), wspólne szlaki (PW) i zależność od czynnika transkrypcyjnego (TFD). Liczba białek ogółem to liczba białek ze wszystkich poprzednich kategorii łącznie. Wszystkie funkcjonalnie powiązane białka uzyskano przy użyciu domyślnych parametrów CaseOLAP LIFT. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 3: Statystyki imputacji etykiet MeSH. W tej tabeli przedstawiono kategorie chorób, numery drzew MeSH używane jako termin nadrzędny dla wszystkich chorób uwzględnionych w kategorii, liczbę artykułów PubMed znalezionych w każdej kategorii w latach 2012-2022 oraz liczbę dodatkowych artykułów uwzględnionych na podstawie etapu imputacji etykiety. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 4: Statystyki budowy grafu wiedzy. W tej tabeli opisano statystyki dotyczące rozmiaru konstruowanego grafu wiedzy, w tym różnych węzłów i typów krawędzi. Wyniki CaseOLAP reprezentują związek między białkiem a kategorią chorób sercowo-naczyniowych (CVD). Kliknij tutaj, aby pobrać tę tabelę.
Tabela 5: Statystyki przewidywania wykresów wiedzy i walidacje. W tej tabeli przedstawiono metryki oceny dla przewidywania powiązań grafu wiedzy dotyczących nowych/ukrytych powiązań białko-choroba. Krawędzie grafu wiedzy zostały podzielone na zestawy danych treningowych i testowych 70/30, a łączność grafów krawędzi została zachowana w obu zestawach danych. Dokładność wskazuje proporcję przewidywań poprawnie sklasyfikowanych, podczas gdy dokładność zrównoważona koryguje nierównowagę klas. Swoistość wskazuje na odsetek negatywnych przewidywań prawidłowo sklasyfikowanych. Precyzja wskazuje proporcję prawidłowych pozytywnych przewidywań spośród wszystkich pozytywnych przewidywań, podczas gdy przywołanie wskazuje proporcję poprawnych pozytywnych przewidywań spośród wszystkich pozytywnych krawędzi (tj. asocjacji białko-choroba zidentyfikowanych za pomocą eksploracji tekstu). Wynik F1 jest średnią harmoniczną precyzji i kompletności. Obszar pod krzywą charakterystyki operacyjnej odbiornika (AUROC) opisuje, jak dobrze model rozróżnia przewidywania pozytywne i negatywne, przy czym 1,0 oznacza doskonały klasyfikator. Obszar pod krzywą precyzja-kompletność (AUPRC) mierzy kompromis między precyzją a kompletnością przy różnych progach prawdopodobieństwa, przy czym wyższe wartości wskazują na lepszą wydajność. Kliknij tutaj, aby pobrać tę tabelę.