Geny związane z egzemą odbytu, geny docelowe SDG oraz wspólne cele
W bazie danych Genecards przesiewowo zbadano łącznie 958 potencjalnych genów kandydackich, a w bazie OMIM 634, wykluczając duplikaty. Aby uzyskać kompleksowe zrozumienie genów związanych z egzemą odbytu, połączono wyniki z wielu baz danych, co pozwoliło uzyskać łącznie 958 odrębnych genów. W konsekwencji precyzyjnie opracowano sieć oddziaływań białko-białko (PPI) specyficzną dla egzemy odbytu. SDG składa się z pięciu tradycyjnych chińskich leków, mianowicie: indigo naturalis, golden cypress, gipsu kalcynowanego, kalaminy oraz Chinese Gall15,16. Głównym składnikiem gipsu kalcynowanego jest bezwodny siarczan wapnia (CaSO4), natomiast głównym składnikiem kalaminy jest węglan cynku (ZnCO3). Indigo naturalis, golden cypress i Chinese Gall mają złożony skład. Z bazy danych TCMSP wynika, że leki te zawierają 92 składniki związkowe, co pozwoliło uzyskać łącznie 867 wiarygodnych celów lekowych (Tabela 1).
Poprzez nałożenie na siebie obu zbiorów danych genów docelowych zidentyfikowano łącznie 149 często współwystępujących genów docelowych (Rysunek 2A), a następnie skonstruowano istotną sieć oddziaływań białko-białko (PPI) białek docelowych (Rysunek 2B). Wykorzystując metodę przesiewową opartą na medianie dla stopnia, bliskości i pośredniości, wyłoniono 59 kluczowych celów jako potencjalne cele lekowe w leczeniu egzemy odbytu. Mediany wyników stopnia, bliskości i pośredniości dla kluczowych celów wynosiły odpowiednio 49, 40,31947 oraz 0,522. 10 genów z najwyższym wynikiem stopnia obejmowało AKT1, TNF, TP53, EGFR, STAT3, SRC, JUN, CASP3, HRAS i PTGS2 (Tabela 2). Geny te wykazują wysoką istotność w kontekście egzemy odbytu.
Ścieżki i sieci obejmujące wspólne cele
W celu analizy 59 kluczowych celów wykorzystano metody wzbogacania KEGG i GO, co ujawniło 218 powiązanych ścieżek oraz ponad 3000 powiązanych procesów biologicznych. Analiza wykazała ścieżki silnie korelujące z białkami w SDG i egzemy odbytu, w tym zakażenie wirusem Cherry simplex virus 1, shigellozę, ścieżkę sygnałową TNF, oporność na inhibitory tyrozynowej kinazy EGFR, zakażenie ludzkim cytomegalowirusem oraz ścieżkę sygnałową receptora limfocytów T (Rysunek 3A). Ścieżki te dotyczą głównie genów takich jak AKT1, TNF, TP53, STAT3, SRC, EGFR i CASP3. Rysunek 3B przedstawia szczegółowy opis genów docelowych i ścieżek. Analizę GO przeprowadzono w odniesieniu do procesów biologicznych (BP), składu komórkowego (CC) i funkcji molekularnej (MF) (Rysunek 4A). Wyniki sugerują, że niniejsze badanie koncentruje się głównie na wspólnych celach dla SDG i egzemy odbytu w procesach biologicznych, z nielicznymi odniesieniami do CC i MF. Szczególnie istotne funkcje biologiczne obejmowały fosforylację peptydylo-tyrozyny, modyfikację peptydylo-tyrozyny, regulację adhezji komórka-komórka, pozytywną regulację adhezji komórkowej, aktywację limfocytów T oraz regulację adhezji komórka-komórka leukocytów (Rysunek 4B-D).
Przewidywanie wiązania aktywnych składników SDG z celami w egzemie odbytu
Na podstawie median wartości stopnia, bliskości i pośredniości wyselekcjonowano 59 kluczowych celów, w tym AKT1, TNF, TP53, EGFR, STAT3, SRC, JUN, CASP3, HRAS oraz PTGS2. Dalsza analiza bazy danych GEO wykazała nadekspresję PPARG, EGFR i TNF, podczas gdy PTPRC, MMP9, MAPK14 i CASP3 wykazywały niedoekspresję w grupie doświadczalnej (zatoka atopowe) (Rysunek 5). Poprzez analizę wzbogacenia wspólnych szlaków genowych ustalono, że geny te uczestniczyły głównie w kaskadzie sygnalizacyjnej TNF oraz szlaku sygnalizacyjnym MAPK. W szlaku sygnalizacyjnym TNF stwierdzono nadekspresję TNF, natomiast niedoekspresję MMP9, MAPK14 i CASP3. W szlaku sygnalizacyjnym MAPK nadekspresję wykazywały EGFR i TNF, natomiast niedoekspresję MAPK14 i CASP3 (Rysunek 6). Na podstawie tych wyników TNF, MAPK14 i CASP3 uznano za potencjalne cele w terapii SDG.
W celu walidacji potencjalnych celów w aktywnych składnikach SDG wykorzystano analizę dokowania, aby sprawdzić dokładność dopasowania między strukturą aktywnego składnika a potencjalnymi białkami docelowymi. Białka te są zaangażowane w różnorodne powiązania funkcjonalne i stanowią węzły o wysokim stopniu w sieci, co sugeruje, że odgrywają one kluczową rolę w odpowiedzi SDG na przeciwodbytniczy egzema. Ujemna wartość energii wiązania w dokowaniu wskazuje na zdolność SDG do dokowania z celami chorobowymi in vivo, przy czym bardziej ujemna wartość oznacza łatwiejsze dokowanie. W niniejszym badaniu uzyskano pomyślne dokowanie molekularne głównych aktywnych składników z kluczowym celem, a energia wiązania była ujemna, z wartościami mniejszymi niż -1 kcal/mol. Indygo i berberubina wykazują dobrą aktywność wiązania, z energią wiązania mniejszą niż -5 kcal/mol (Tabela 3, Rycina 7). Podsumowując, wyniki te stanowią dalszy dowód na to, że białka odpowiadające tym loci genowym mogą pełnić rolę celów SDG w egzemie przeciwodbytniczej.

Rycina 1: Schemat analizy farmakologii sieciowej. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; TCMSP, Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform; GEO, Gene Expression Omnibus. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Diagram Venna i sieć PPI wspólnych celów. (A) Diagram Venna przedstawiający część wspólną celów leku i celów choroby. (B) Sieć PPI wspólnych celów opracowana za pomocą bazy STRING. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Analiza wzbogacenia szlaków KEGG. (A) Analiza wzbogacenia szlaków KEGG. 10 najważniejszych szlaków KEGG zostało uszeregowanych według rosnących wartości P. (B) Połączenie między szlakiem a celem: szlak (żółty), cele (czerwone). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Analiza wzbogacenia GO. (A) Wyniki GO dla trzech ontologii. (B) Wykres pęcherzykowy dla procesów biologicznych (BP). (C) Wykres pęcherzykowy dla komponentów komórkowych (CC). (D) Wykres pęcherzykowy dla funkcji molekularnych (MF). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 5: Wyniki przewidywania potencjalnych celów. (A) Mapa ciepła ekspresji genów centralnych (hub genes) w bazie danych GEO; grupa A to grupa badawcza (atopowe zapalenie skóry), a grupa B to grupa kontrolna (brak atopowego zapalenia skóry); (B) węzły sieci PPI reprezentują białka, a krawędzie reprezentują relacje między nimi. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 6: Szlak sygnalizacyjny. (A) Szlak sygnalizacyjny MAPK. (B) Szlak sygnalizacyjny TNF. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 7: Dokowanie molekularne genów kluczowych i składników. Kolor magenta reprezentuje kluczowe składniki SDG, a kolor niebieski reprezentuje reszty genów kluczowych. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Tradycyjne leki chińskie | Składniki aktywne |
| Indigo naturalis | 9α,13α-dihydroksyizopropylidenyloisatyna, a, bisindigotyna, indykant, isatan B, isatyna, a, izoorientyna, izoskoparyna, izowiteksyna, (+)-izolarycyrezynol, 10h-indolo[3,2-b]chinolina, izoindygo, saponaryna, indygo, tryptantryna, 6-(3-oksoindolin-2-ylideno)indolo[2,1-b]chinazolin-12-on |
| Indirubina, beta-sitosterol, laricyrezynol, nonakosan, isowiteksyna, dotriakontanol |
| Cyprys złoty | berberyna, koptyzyna, daurycyna (8CI), jawanicyna, (±)-lyoniresinol, kihadalakton A, kwas obakunowy, obakunon, felawina, Phellavin_qt, felodendryna, delta 7-stigmastenol, felopteryna, wanillolosyd, koniferyna, dehydrotanszynon II A, delta7-dehydrosoforamina, amurensyna, Amurensin_qt, dhydronilotycyna, hispidol B, kihadalakton B, kihadanina A, nilotycyna, nomilina, rutakarpina, skimmianina, chelerytryna, stigmasterol, worenina, ferulan kampesterolu, kawidyna, candletoksyna A, hericenon H, hispidon, syrygina, beta-sitosterol, magnograndiolid, (2S,3S)-3,5,7-trihydroksy-2-(4-hydroksyfenylo)chroman-4-on, palmidyna A, magnofloryna, menisporfina, palmatyna, fumaryna, izokorypalmina, kwercetyna, sitoglusyd, friedelina |
| STOCK1N-14407, jatrorryzyna, menisperyna, phellamurina_qt, (S)-kanadyna, kolumbamina, poriferast-5-en-3beta-ol, magnofloryna, berberrubina, phellodendryna, limonina, hyperyna, kampesterol, SMR000232320, kantyna-6-on, 4-[(1R,3aS,4R,6aS)-4-(4-hydroksy-3,5-dimetoksifenylo)-1,3,3a,4,6,6a-heksahydrofuro[4,3-c]furan-1-ylo]-2,6-dimetoksifenol, dihydronilotycyna, melianon, phellochina, thalifendyna, wanillolosyd, aurapten |
| Gips wypalany | bezwodny siarczan wapnia (CaSO4) |
| Kalamina | węglan cynku (ZnCO3) |
| Gal chiński | digallat |
Tabela 1: Składniki aktywne w SDG.
| Gen | Stopień | Centralność pośrednictwa | Centralność bliskości |
| AKT1 | 204 | 1669.1692 | 0.765625 |
| TNF | 202 | 1988.4543 | 0.761658 |
| TP53 | 190 | 1590.9288 | 0.73134327 |
| EGFR | 174 | 686.3063 | 0.7033493 |
| STAT3 | 168 | 673.03723 | 0.6869159 |
| SRC | 162 | 568.1574 | 0.69014084 |
| JUN | 162 | 435.33737 | 0.6805556 |
| CASP3 | 156 | 483.45276 | 0.67431194 |
| HRAS | 148 | 515.28815 | 0.65625 |
| PTGS2 | 134 | 761.34094 | 0.6447368 |
Tabela 2: Charakterystyka 10 najważniejszych genów hubowych.
| Cel (PDB ID) | Powinowactwo (kcal/mol) |
| Indygo | Berberrubina | Digallat |
| TNF (1A8M) | -5.96 | -5.19 | -2.22 |
| MAPK14 (1A9U) | -5.51 | -5.41 | -1.93 |
| CASP3 (1CP3) | -5.77 | -4.98 | -1.06 |
Tabela 3: Energia wiązania w dokowaniu molekularnym składników i genów kluczowych.