$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Bu protokol, ağ toksikolojisi, transkriptomik, makine öğrenimi ve moleküler kenetlenme gibi yöntemlerden yararlanarak asetaminofen kaynaklı karaciğer hasarının olası biyobelirteçlerini tanımlamak için hesaplamalı bir yöntem ortaya koymaktadır (Şekil 1). Protokol, biyoinformatik araçlara, transkriptomik veri setlerine ve moleküler kenetlenme yazılımlarına erişimi olan araştırmacılara yöneliktir.
Prosedür
Adım 1: Asetaminofen hedeflerinin tanımlanması
Asetaminofenin SMILES temsili (APAP) PubChem'den alın. APAP'ın potansiyel moleküler hedeflerini tahmin etmek için çevrimiçi platformları (ChEMBL, SwissTargetPrediction, STITCH, SEA) kullanın. Tahmin edilen hedefleri entegre edin ve deduplicate ederek 140 yüksek güvenliğe sahip APAP hedefinden oluşan bir liste oluşturun.
Adım 2: Hepatotoksisite hedeflerinin tanımlanması
GeneCards veritabanından hepatotoksisiteyle ilgili genleri alın. Listeyi derleyip deduplicate ederek hepatotoksisiteyle ilgili 657 genden oluşan yedeksiz bir set oluşturulur. APAP hedefleri ile hepatotoksisiteye bağlı genler arasında örtüşen genleri Venn diyagramı kullanarak tespit edin.
Adım 3: Transkriptomik Veri Ön İşlemi
GSE74000 veri setini GEO'dan indirin. DESeq2 kullanarak ham ifade verilerini ön işleyin: düşük ekspresyona sahip genleri çıkarın, boyut faktörleriyle normalize edin ve varyans stabilizasyon dönüşümü (VST) uygulayın. Limma ve DESeq2 kullanarak diferansiyel ifade analizi yapın; ayarlanmış p-değeri < 0.05 ve |log2FC| > 1.
Adım 4: Fonksiyonel Zenginleştirme Analizi
Üst üste binen genleri Gen Ontolojisi (GO), KEGG yolu, doku ifadesi ve hastalık korelasyon analizleri için STRING'e yükleyin. Fonksiyonel zenginleştirme sonuçlarını balon grafikleri ve ısı haritaları kullanarak görselleştirin.
Adım 5: Özellik Gen Seçimi İçin Makine Öğrenimi
Rastgele Orman sınıflandırıcısı (n_estimators=500, max_depth=10) uygulayarak örtüşen APAP ve hepatotoksisite genlerinden gelen özellik genlerini önceliklendirin. Modelin performansını Out-of-Bag (OOB) hata ve özellik önemi puanları kullanarak değerlendirin. Daha fazla analiz için en iyi 20 özellikli genden birini seçin.
Adım 6: Moleküler Kenetlenme
PubChem'den APAP yapısını (CID 1983) ve PDB'den protein reseptörlerini (ESR1: PDB ID 1SJ0, PNP: PDB ID 1V2H) alın. Ligand ve reseptör dosyalarını hazırlayın: PDB formatına dönüştürün, polar hidrojenler ekleyin, yükleri atayın ve PDBQT dosyaları olarak kaydedin. AutoDock Tools'ta proteinin aktif yerini kapsayan kenetlenme ızgarasını tanımlayın. AutoDock Vina ile moleküler kenetlenme gerçekleştirin, kapsamlılık 8'e ayarlanmış ve bağlanma eğilimleri ile etkileşimleri analiz edin. Bağlama konformasyonlarını ve anahtar etkileşimleri analiz etmek için PyMOL kullanarak kenetlenme sonuçlarını görselleştirin.
Adım 7: İstatistiksel Analiz
T-testleriyle istatistiksel anlamlılığı belirlemek ve Benjamini-Hochbergyöntemi 31 ile çoklu karşılaştırmalar için p-değerleri ayarlamak. Isı haritaları ve dağılma grafikleri kullanarak istatistiksel olarak anlamlı ilişkileri görselleştirin.
Malzemeler ve Yöntemler
Asetaminofen hedeflerinin tanımlanması
APAP'ın potansiyel moleküler hedeflerini belirlemek için önce PubChem veritabanından SMILES temsilini aldık. Daha sonra Kimyasal Avrupa Moleküler Biyoloji Laboratuvarı (ChemBL)32, İsviçre HedefTahmini 33, Kimyasallar ve Hedeflerin Etkileşimi Arama Aracı (STITCH)34 ve Benzerlik topluluk yaklaşımı (SEA)35 gibi çeşitli çevrimiçi platformları kullanarak olası hedefleri tahmin ettik. Bu araçlardan elde edilen sonuçları entegre ettikten sonra, APAP için yüksek güvenliğe sahip hedefler seti seçtik. Tablo 1, bu çalışmada kullanılan temel gen kategorilerini tanımlar, veri analizi ve biyolojik yorumlamadaki rollerini netleştirir. Bu terimlerin tutarlı kullanımı, sonuçlarımızın net bir şekilde iletmesini sağlar.
Hepatotoksisite hedeflerinin tanımlanması
Hepatotoksisiteyle ilişkili potansiyel genler GeneCards veritabanından alındı. Tüm tanımlanan genler derlendi, kopyalar çıkarıldı ve sonraki analizler için yedek olmayan bir liste oluşturuldu. GSE74000 veri seti 15 Mart 2024'te GEO deposundan indirildi. Ham ifade verileri limma paketi (quantil normalizasyon) kullanılarak işlenip normalleştirildi. Diferansiyel ifade analizi, ampirik Bayes küçülmesi ile lineer modelleme kullanılarak gerçekleştirilmiştir. Genler, ayarlanmış p-değeri < 0.05 (Benjamini-Hochberg FDR) ve |log₂FC| ile eşleşen genler. > 1 önemli kabul edildi. DEG'lerin görselleştirilmesi, ggplot2 ile oluşturulan volkan grafikleri ve ısı haritaları kullanılarak gerçekleştirildi.
Veri ön işleme
Ham sayım verileri DESeq2 kullanılarak önceden işlendi. Düşük ekspresyon genler, örneklerin en az %70'inde CPM >1 tespit eşiği kullanılarak çıkarıldı. Kütüphane boyutu normalizasyonu, denklem (1)'de tanımlandığı şekilde DESeq2 boyut faktörleri kullanılarak gerçekleştirildi:
(1)
Burada medyan oran boyut faktörü sj ile gösterilir. Ortalama-varyans ilişkilerini stabilize etmek için, (2) denkleminde varyans-stabilizator dönüşüm (VST) kullanılmıştır:
(2)
DEG'ler, Benjamini-Hochberg düzeltmesi ile Wald testi kullanılarak tanımlandı; denklemde (3) kullanılan genler anlamlı olarak değerlendirildi:
(3)
DESeq2'den türetilen DEG kümesi, denklem (4) olarak tanımlanmıştır:
(4)
Bu küme (X), Limma Trend (Y) ve Limma Voom (Z) ile birlikte uzlaşma stratejisinde kullanıldı.
PPI ağ inşası
Venn diyagramları, APAP ve hepatotoksisite hedefleri arasındaki ortak genleri belirlemek için kullanıldı. Üst üste düşen genler, daha sonra Etkileşimli Genler/Proteinlerin Geri Alınması Arama Aracı (STRING) veritabanına yüklenerek PPI ağları oluşturuldu.
Çok boyutlu fonksiyonel zenginleştirme analizi
İlk olarak STRING web sitesini kullanarak Gen Ontolojisi (GO), Kyoto Genler ve Genomlar Ansiklopedisi (KEGG), doku ifadesi ve örtüşen genlerin APAP ve hepatotoksisite için hastalıkla ilgili fonksiyonel analizlerini yaptık. Sonrasında, Sendo Academic Tools kullanarak hepatotoksisite için diferansiyel genlerin GO, KEGG ve Gen Seti Zenginleştirme Analizi (GSEA) (REACTOME) zenginleştirme analizlerini yaptık.
Rastgele orman analizi
APAP kaynaklı hepatotoksisite transkriptom verilerinde örtüşen APAP ve hepatotoksisite genlerinden en iyi 20 özellik genini belirlemek için makine öğrenimi uyguladık. Rastgele Orman sınıflandırıcı uygulandı; 500 ağaç (n_estimators=500), maksimum ağaç derinliği 10 (max_depth=10), bir düğümü bölmek için minimum 2 örnek gerekiyordu (min_samples_split=2) ve Gini safsızlık kriteri (kriter='gini')36. Modelin performansı, 0'a yakın bir değerin daha yüksek tahmin doğruluğunu gösterdiği Out-of-Bag (OOB) hatası kullanılarak değerlendirildi. Özellik önemi puanları, her genin katkısını sıralamak için Rastgele Orman analizine dayanarak hesaplandı ve görselleştirildi; denklem (5)'te olduğu gibi.
(5)
N toplam örnek sayısı, yi 1(·) gösterge fonksiyonu, 1 koşul doğruysa ve 0 aksi halde;
yalnızca III'ün eğitime dahil edilmediği ağaçlarla önerilen örnek iii etiketi olarak öngörüldü.
Karakterize edilmiş genlerin diferansiyel ifadesi
Transkriptomik verilerdeki özellik genlerinin ifade farklılıkları keman grafikleri kullanılarak görselleştirildi. İstatistiksel olarak anlamlı farklılıklar gösteren biyobelirteçler, APAP kaynaklı hepatotoksisitenin potansiyel yeni biyobelirteçleri olarak daha ileri inceleme için belirlendi.
Moleküler kenetleme
Küçük molekül bileşikleri (CID 1983) PubChem veritabanından alınmış, protein reseptörleri ESR1 ve PNP (PDB ID'leri 1SJ0 ve 1V2H) ise Protein Veri Bankası'ndan indirilmiştir. Ligand yapıları OpenBabel kullanılarak PDB formatına dönüştürüldü ve AutoDock Tools'ta polar hidrojenler ekleyerek, Gasteiger yükleri atama, dönebilir bağlar tanımlanarak ve PDBQT formatında tasarruf edilerek önceden işlendi. Protein reseptörleri, su molekülleri ve birlikte kristalleşmiş ligandlar çıkararak PyMOL kullanılarak hazırlanmış, ardından polar hidrojenler eklenmiş ve AutoDock Tools kullanılarak Kollman yüklerinin atanması yapılmış, PDBQT dosyaları olarak kaydedilmiştir.
Moleküler kenetlenmede, AutoDock yazılımı, protein37'nin aktif bölgesini kaplayan kenetlenme ızgarasını tanımlamak için kullanıldı. Grid kutusu, koordinatlarda (x = XX·XX, y = YY· YY,z = ZZ· ZZ) boyutları 40 × 40 × 40 şve ızgara aralığı 0.375 şile bağlantı cebinin tam kapsayışını garanti eder. AutoDock Vina, ligand-protein bağlanma modlarını ve bağlanma eğilimlerini hesaplamak için kullanıldı; tüketicilik parametresi 8'e ayarlandı ve her ligand için en üst dokuz bağlanma pozu oluşturuldu.
Kenetlenme protokolü doğrulaması, birlikte kristalleşmiş ligandın aktif siteye yeniden kenetlenmesiyle gerçekleştirildi ve RMSD değeri < 2.0 şelde edildi ve kenetlenme prosedürünün güvenilirliği doğrulandı. Kenetlenme sonuçları, PyMOL kullanılarak bağlanma konformasyonları ve hidrojen bağları dahil ana etkileşimleri analiz etmek için görselleştirildi.
Kenetlenme simülasyonları, Bileşik X'in protein Y'nin bağlanma bölgesine uyduğunu ve potansiyel hidrojen bağları ve hidrofobik temaslar oluşturduğunu ve tahmini bağlanma enerjisinin -8,5 kcal/mol olduğunu öngördü.
Sorun giderme ve olası değişiklikler
Önerilen iş akışının sağlamlığını ve tekrarlanabilirliğini artırmak için, birkaç sorun giderme değerlendirmesi ve olası değişiklikler not edilmelidir. Beklenmedik derecede düşük sayıda farklı şekilde ifade edilen gen (DEG) tespit edilirse, kullanıcıların normalizasyon prosedürünü doğrulaması, doğru grup etiketlemesini doğrulaması ve |log2 katlama değişimini ayarlamayı düşünmeleri tavsiye edilir. eşik olarak uygun yanlış keşif oranı (FDR) kontrolünü korurken. Buna karşılık, aşırı sayıda DEG'in elde edilmesi durumunda, daha katı FDR sınırları uygulamak veya düşük varyanslı genlerin filtrelenmesi diferansiyel ekspresyon analizinden önce spesifikliği artırabilir.
Toplu etkiler, ana bileşen analizi (PCA) gibi keşif analizlerinde kümeleme kalıplarını etkileyebilir. Örnekler ağırlıklı olarak biyolojik koşuldan ziyade toplu olarak kümelenirse, toplu düzeltme yöntemleri (örneğin, ComBat gibi ampirik Bayes yaklaşımları) uygulanmalı ve örnek meta verileri tutarlılık açısından dikkatlice yeniden değerlendirilmelidir.
Rastgele Orman tabanlı özellik seçimi için, yüksek Out-of-Bag (OOB) hata oranları veya kararsız özellik sıralamaları, model yapılandırmasının optimal olmadığını gösterebilir. Bu tür durumlarda, ağaç sayısını artırmak, mtry parametresini ayarlamak veya uzlaşma özellik seçimiyle tekrar edilen model çalışmaları model kararlılığını ve tahmin güvenilirliğini artırabilir. Ayrıca, model sağlamlığını daha iyi değerlendirmek için çapraz doğrulama stratejileri kullanılabilir.
Güvenilirliği artırmak için, kullanıcılar alternatif DEG eşikleri veya makine öğrenimi parametre ayarları kullanarak analizi isteğe bağlı olarak tekrarlayabilir ve tanımlanan özellik genlerinin tutarlılığını karşılaştırabilirler. Bu tür duyarlılık analizleri, temel bulguların belirli parametre seçimleriyle yönlendirilmemesini sağlar ve benzer transkriptomik veri setleri arasında iş akışının tekrarlanabilirliğini destekler.
İstatistiksel analizler
İstatistiksel anlamlılık t-testi kullanılarak belirlendi ve karşılaştırma için p-değerleri bildirildi. Gen ifade seviyeleri ile hepatotoksisiteye bağlı fenotipler arasındaki istatistiksel korelasyonlar, veri normalliğine bağlı olarak Pearson ve Spearman korelasyon katsayıları kullanılarak değerlendirildi. P-değerleri, Benjamini-Hochberg yöntemiyle birden fazla karşılaştırma için ayarlandı. Önemli ilişkiler, ısı haritaları ve saçılım grafikleriyle görselleştirildi ve transkriptomik ilişkilerin sağlam değerlendirilmesini sağladı.