Moleküler dinamik

Moleküler dinamik (MD), atomların ve moleküllerin fiziksel hareketlerini incelemek için bir bilgisayar simülasyon yöntemidir. Atomların ve moleküllerin sabit bir süre boyunca etkileşime girmesine izin verilir ve bu da sistemin dinamik evrimi hakkında bilgi verir. En yaygın versiyonda, atomların ve moleküllerin yörüngeleri, parçacıklar ve bunların potansiyel enerjileri arasındaki kuvvetlerin çoğu zaman atomlararası potansiyeller veya moleküler mekanik kuvvet alanları kullanılarak hesaplandığı, etkileşen parçacıkların bir sistemi için Newton'un hareket denklemlerinin sayısal olarak çözülmesiyle belirlenir. Metot ilk olarak 1950'lerin sonunda teorik fizik alanında geliştirildi, ancak günümüzde çoğunlukla kimyasal fizik, malzeme bilimi ve biyomoleküllerin modellenmesinde uygulanmaktadır.
Moleküler sistemler tipik olarak çok sayıda parçacıktan oluştuğundan, bu tür karmaşık sistemlerin özelliklerini analitik olarak belirlemek mümkün değildir; MD simülasyonu, sayısal yöntemleri kullanarak bu sorunu çözmektedir. Bununla birlikte, uzun MD simülasyonları matematiksel olarak kötü koşullandırılmıştırlar (yüksek koşul sayısına sahiptir), yani sayısal entegrasyonda uygun algoritma ve parametrelerin seçimi ile minimize edilebilecek, ancak tamamen ortadan kalkamayan kümülatif hatalar oluştururlar.
Ergodik hipoteze uyan sistemler için, bir moleküler dinamik simülasyonunun evrimi, sistemin makroskopik termodinamik özelliklerini belirlemek için kullanılabilir: bir ergodik sistemin zaman ortalamaları, mikrokanonik topluluk ortalamalarına karşılık gelir. MD, aynı zamanda doğanın kuvvetlerini canlandırarak ve atomik bir ölçekte moleküler hareketi inceleme şansı sağlayarak "sayılarla istatistiksel mekanik" ve " Laplace'in Newton mekaniğinin vizyonu" olarak da adlandırılmıştır.
Tarih
Monte Carlo simülasyonlarının önceki başarılarından sonra, yöntem ilk kez 1950'lerin ortalarında Fermi, Pasta, Ulam ve Tsingou tarafından geliştirilmiştir. 1957'de Alder ve Wainwright, sert küreler arasında mükemmel elastik çarpışmaları simüle etmek için bir IBM 704 bilgisayarı kullandı. 1960 yılında Gibson ve ark. yapışkan bir yüzey kuvveti ile birlikte Born-Mayer tipi itici etkileşimi kullanarak katı bakırın radyasyon hasarını simüle etti. 1964 yılında, Rahman, bir Lennard-Jones potansiyeli kullanan, sıvı argonun dönüm noktası simülasyonlarını yayınladı. Kendi kendine yayılma katsayısı gibi sistem özelliklerinin hesaplanması deneysel verilerle deneysel verilerle iyi karşılaştırıldı.
Moleküler dinamikleri bilgisayarlarla simüle etmek mümkün olmadan önce bazı insanlar atom hareketinin fiziksel modeller (örneğin makroskopik küreler) oluşturarak modellenmesi için yoğun bir çalışma üstlendi. Amaç, bu küreleri bir sıvının yapısını taklit edecek şekilde düzenlemek ve davranışını incelemek için kullanmaktı. JD Bernal, 1962de şunları söyledi: ... Çok sayıda lastik top aldım ve bunları 2.75 ila 4 inç arasında değişen farklı uzunluklardaki çubuklarla bir araya getirdim. Bunu ilk etapta mümkün olduğunca rastgele yapmayı denedim, kendi ofisimde çalışıyor, her beş dakikada bir rahatsız ediliyor ve rahatsız edilmeden önce ne yaptığımı hatırlamıyordum.
Uygulama alanları
İlk olarak teorik fizikte kullanılan MD yöntemi, malzeme biliminde kısa süre sonra popülerlik kazanmıştır ve 1970'lerden beri biyokimya ve biyofizikte de yaygındır. MD, X ışını kristalografisi veya NMR spektroskopisinden elde edilen deneysel kısıtlamalara dayanarak proteinlerin ve diğer makromoleküllerin 3 boyutlu yapılarını geliştirmek/daha iyi hale getirmek için sıklıkla kullanılır. Fizikte, MD, ince film büyümesi ve iyon alt ekimi gibi doğrudan gözlemlenemeyen atomik seviye fenomenlerin dinamiklerini incelemek ve ayrıca henüz yaratılmayan veya henüz yaratılamayan nanoteknolojik cihazların fiziksel özelliklerini incelemek için kullanılır. Biyofizik ve yapısal biyolojide, bu yöntem proteinler ve nükleik asitler gibi makromoleküllerin hareketlerini incelemek için sıkça uygulanır. Burdan elde edilecek veri ligand yerleştirme gibi biyofiziksel deneylerin sonuçlarını yorumlamak ve diğer moleküller ile etkileşimleri modellemek için faydalı olabilir. Prensipte MD, polipeptit zincirinin rastgele katlanmasını simüle ederek protein yapısının ab initio tahmini için kullanılabilir.
Uygulama örnekleri
Moleküler dinamikler birçok bilim dalında kullanılmaktadır.
Basitleştirilmiş bir biyolojik katlama işleminin ilk MD simülasyonu 1975'te yayınlandı. Nature'da yayınlanan simülasyon, modern bilgisayarlı protein katlamanın alanının yolunu açtı.
Bir biyolojik sürecin ilk MD simülasyonu 1976'da yayınlandı. Nature'da yayınlanan simülasyonu, protein hareketini sadece yardımcı bir faktör olarak değil de bir temel işlev olarak gerekli görmenin yolunu açtı.
MD, enerjik nötron ve iyon ışınlanmasının katı ve katı yüzeyler üzerindeki etkileri gibi ısı sivri uç rejimindeki (heat spike regime) çarpışma basamaklarını (çarpışma şelalerini) ele almak için standart bir yöntemdir.
Gaucher Hastalığına neden olan en yaygın protein mutasyonu N370S'in moleküler temelini tahmin etmek için MD simülasyonları başarıyla uygulandı. Bunu takip eden daha sonraki bağımsız olarak yayınlanan yayında ise, bu kör tahminlerin aynı mutant üzerindeki deneysel çalışma ile yüksek bir korelasyon gösterdiği gösterilmiştir.
Yüzey yüklerinin ince su filmlerinin metal yüzeyler üzerindeki ayrışma basıncı üzerindeki etkisini araştırmak için MD simülasyonları kullanılmıştır.
MD simülasyonları transmisyon elektron mikroskobu görüntüsünün özelliklerini anlamak için çok kesitli görüntü simülasyonları ile birlikte kullanılır
Hamiltonian-ikilemi algoritmaları ile birleştirilmiş MD hesaplamaları, DNA çift iplikçiklerinin hidrofobik ve hidrofilik tek duvarlı karbon nanotüpler üzerine kapsülleme termodinamiklerini araştırmak için kullanılmıştır.
Aşağıdaki biyofiziksel örnekler, çok büyük boyutta bir sistem (bir virüsün tamamı) veya çok uzun simülasyon sürelerde (1.112 milisaniyeye varan süreler) simülasyonları üretmek için kayda değer çabaları göstermektedir:
Uydu tütün mozaik virüsünün tamamının MD simülasyonu (STMV) (2006, Boyut: 1 milyon atom, Simülasyon süresi: 50 ns, program: NAMD ): Bu virüs, Tütün Mozaik Virüsü (TMV) enfeksiyonunun belirtilerini kötüleştiren küçük, ikosahedral bitki virüsüdür. Viral düzeneğin mekanizmalarını araştırmak için moleküler dinamik simülasyonlar kullanılmıştır. STMV partikülünün tamamı, viral kapsidi (kaplama) oluşturan bir proteinin 60 özdeş kopyasından ve bir 1063 nükleotidli tek zincirli RNA genomundan oluşur. Bir anahtar bulgu, kapsidin içinde RNA olmadığında çok dengesiz olmasıdır. Simülasyonun tamamlanması 2006 model bir masaüstü bilgisayarda yaklaşık 35 yıl alacaktır. Bundan dolayı, görev birçok paralel işlemci ile aralarında sürekli iletişim oluşturarak yapıldı.
Villin Başlığının tüm atom detaylarında katlanması simülasyonları (2006, Boyut: 20.000 atom; Simülasyon süresi: 500 =s = 500.000 ns, Program: Folding@home ) Bu simülasyon, dünya genelinde 200.000 CPU'nun katıldığı kişisel bilgisayarlarda yürütülmüştür. Bu bilgisayarlarda, Viford Pande tarafından Stanford Üniversitesi'nde koordine edilen geniş çaplı bir dağıtılmış hesaplama çalışması olan Folding @ home programı kuruldu. Villin Başlık proteininin kinetik özellikleri sürekli gerçek zamanlı iletişimden yoksun CPU'lar tarafından birçok bağımsız kısa yörünge (trajectories) kulanılarak araştırılmıştır. Kullanılan bir yöntem de, belirli bir başlangıç konformasyonunun açılmasından önce katlanma olasılığını ölçen Pfold değer analiziydi. Pfold geçiş durumu yapıları ve katlama yolağı boyuncaki konformasyon sırası hakkında bilgi verir. Pfold hesaplamasındaki her bir yörünge nispeten kısa olabilir, ancak birçok bağımsız yörüngeye ihtiyaç duyulur.
Bu moleküler simülasyonlar, malzeme çıkarma mekanizmalarını, araç geometrisinin etkilerini, sıcaklık ve kesme hızı ve kesme kuvvetleri gibi işlem parametrelerini anlamak için kullanılmıştır.
Moleküler dinamik algoritmaları
Korunmuş Coulomb Potansiyelleri Örtük Solvent Modeli
Entegratörler
Simplektik entegratör
Verlet-Stoermer entegrasyonu
RungeKutta entegrasyonu
Beeman algoritması
Kısıtlama algoritmaları (kısıtlı sistemler için)
Kısa menzilli etkileşim algoritmaları
Hücre listeleri
Verlet listesi
Bağlı etkileşimler
Uzun menzilli etkileşim algoritmaları
Ewald toplamı
Parçacık ağı Ewald toplamı (PME)
Parçacık - partikül - partikül ağı (P3M )
Kaydırılmış kuvvet yöntemi
Paralelleştirme stratejileri
Etki alanı ayrıştırma yöntemi (Paralel hesaplama için sistem verilerinin dağıtımı)
Ab-initio moleküler dinamik
CarParrinello moleküler dinamiği
MD simülasyonları için özel donanım
Anton - MD simülasyonlarını yürütmek için tasarlanmış uzman, devasa paralel bir süper bilgisayar
MDGRAPE - Özellikle protein yapı tahmini için oluşturulmuş özel amaçlı bir moleküler dinamik simülasyonları sistemi
MD simülasyonları için bir donanım olarak grafik kartı
GPU'da moleküler modelleme