AutoDock & AutoDock Vina: Moleküler Docking'in Gücü — Adım Adım Kullanım Kılavuzu
Bu rehber, AutoDock ve AutoDock Vina yazılımları ile küçük moleküllerin (ligandların) proteinlere nasıl bağlandığını, bağlanma enerjilerini ve en olası konformasyonlarını (pose) nasıl hesaplayacağınızı adım adım açıklar. Bu açık kaynaklı yazılımlar, ilaç keşfi ve yapısal biyoloji araştırmalarında endüstri standardı haline gelmiştir.
💡 Bilmeniz Gereken: AutoDock Vina, orijinal AutoDock’tan çok daha hızlı ve çoğu durumda daha doğru sonuçlar verir. Yeni projelerinizde AutoDock Vina kullanmanız önerilir.
1. AutoDock ve AutoDock Vina Nedir?
AutoDock, Scripps Araştırma Enstitüsü tarafından geliştirilen, ligand-protein etkileşimlerini simüle eden ilk nesil açık kaynaklı docking yazılımıdır. AutoDock Vina, onun yerine geçen, daha hızlı, daha kolay kullanılabilir ve daha yüksek doğruluk sunan yeni nesil versiyondur.
Temel Özellikler:
Ligand esnekliği: Ligandın rotasyon ve konformasyon değişikliklerini otomatik olarak hesaplar.
Çoklu pose çıktısı: Birden fazla olası bağlanma modunu sıralar.
Açık kaynak ve ücretsiz: Ticari kullanım dahil herkes tarafından kullanılabilir.
Komut satırı tabanlı: Otomasyon ve yüksek verimli taramalar için idealdir.
GUI desteği: AutoDock Tools (ADT) veya PyRx gibi arayüzlerle kolaylaştırılabilir.
2. Adım Adım Kullanım Kılavuzu
Adım 1: Yazılımı İndirme ve Kurma
AutoDock Vina tamamen ücretsiz ve açık kaynaktır. Resmi GitHub sayfasından işletim sisteminize uygun versiyonu indirin:
Protein (.pdbqt): PDB formatındaki protein yapısını, AutoDock Tools (ADT) veya OpenBabel ile .pdbqt formatına dönüştürün. Su moleküllerini ve metal iyonlarını temizleyin, polar hidrojenleri ekleyin.
Ligand (.pdbqt): Ligandınızı çizip (Avogadro, ChemDraw), enerji minimizasyonu yaptıktan sonra .pdbqt formatına dönüştürün. Dönen bağları (rotatable bonds) tanımlayın.
Ayrıca, protein üzerindeki bağlanma bölgesini (search space) tanımlamanız gerekir. Bunun için:
Center_X, Center_Y, Center_Z: Aktif bölgenin merkez koordinatları (Å).
Size_X, Size_Y, Size_Z: Arama kutusunun boyutları (Å). Genellikle 20x20x20 Å yeterlidir.
Adım 3: Docking Hesaplamasını Çalıştırma
Terminal veya komut istemcisini açın ve aşağıdaki komutu çalıştırın:
Çıktı: En düşük bağlanma enerjisine (kcal/mol) sahip pose en iyi tahmin olarak kabul edilir.
Adım 4: Sonuçları Yorumlama ve Görselleştirme
Çıktı dosyası (output.pdbqt) birden fazla pose içerir. Bu dosyayı PyMOL, ChimeraX veya AutoDock Tools ile açarak:
En iyi 3-5 pose’u karşılaştırın.
Ligand ile protein arasındaki hidrojen bağlarını, hidrofobik etkileşimleri ve tuz köprülerini analiz edin.
Bağlanma ceperinin şekli ve ligandın uyumu hakkında yorum yapın.
Örnek Log Çıktısı:
MODEL 1
REMARK VINA RESULT: -8.2 0.000 0.000
MODEL 2
REMARK VINA RESULT: -7.5 0.000 0.000
Burada -8.2 kcal/mol, en düşük (en olası) bağlanma enerjisidir.
3. Pratik Uygulamalar ve Örnek Senaryolar
Senaryo 1: SARS-CoV-2 Proteaz İnhibitörü Taraması
1. Hedef protein: 6LU7 (PDB’den indirin, su ve inhibitörü temizleyin).
2. Ligand kütüphanesi: ZINC15 veya PubChem’den 100 adet küçük molekül seçin.
3. Her ligand için yukarıdaki adımları otomatikleştirin (bash/python script ile).
4. En düşük enerjiye sahip 5 ligandı seçip ChimeraX ile etkileşim analizi yapın.
Senaryo 2: Bağlanma Bölgesi Optimizasyonu
1. İlk docking’de bağlanma bölgesini geniş tutun (30x30x30 Å).
2. En iyi pose’un koordinatlarını alın ve yeni merkez olarak belirleyin.
3. Bölgeyi daraltın (15x15x15 Å) ve tekrar docking yapın. Daha hassas sonuçlar alırsınız.
4. Yaygın Sorunlar ve Çözümleri
“Command not found” hatası: Vina’nın kurulu olduğu klasörü PATH’e eklemeyi unutmuş olabilirsiniz. Linux/macOS’te export PATH=$PATH:/path/to/vina komutunu kullanın.
“Parse error” veya “File not found”: .pdbqt dosyalarının doğru formatlandığından ve yolda boşluk olmadığından emin olun.
Yavaş çalışma:--cpu parametresi ile işlemci sayısını artırabilirsiniz (örneğin, --cpu 4).
5. Bilgiyi Test Et
Soru 1: AutoDock Vina’nın orijinal AutoDock’a göre en büyük avantajı nedir?
Cevap: Çok daha hızlı ve çoğu durumda daha doğru sonuçlar verir.
Soru 2: Docking hesaplaması için gerekli olan iki temel girdi dosyası nedir?
Cevap:protein.pdbqt ve ligand.pdbqt
Soru 3: Çıktı dosyasında “REMARK VINA RESULT: -9.5 0.000 0.000” satırı ne anlama gelir?
Cevap: Bu pose’un bağlanma enerjisi -9.5 kcal/mol’dir. Düşük (daha negatif) değerler daha güçlü bağlanmayı gösterir.