Molekulska mehanika koristi klasičnu mehaniku za modelovanje molekularnih sistema. Potencijalna energija svih sistema u molekulskoj mehanici se izračunava koristeći polja sila. Molekulska mehanika se može koristiti za izučavanje molekulskih sistema u opsegu veličina i kompleksnosti od malih do velikih bioloških sistema ili materijalnih skupina sa od više hiljada do više miliona atoma.
Sve-atomistički molekulsko mehanički metodi imaju sledeća svojstva:
Svaki atom se simulira kao jedna čestica
Svakoj čestici se dodeljuje prečnik (tipično Van der Valsov radijus), polarizabilnost, i konstantno neto naelektrisanje (generalno izvedeno iz kvantnih proračuna i/ili eksperimentalnih merenja)
Interakcije vezivanja se tretiraju kao opruge sa ravnotežnim rastojanjem jednakim eksperimentalnoj ili izračunatoj dužini veze
Varijante ovog pristupa su moguće. Na primer, mnoge simulacije su istorijski koristile ujedinjenu atomku reprezentaciju u kojoj je svaka terminalna metil grupa ili intermedijarna metilenska jedinica smatrana jednom česticom, i veliki proteinski sistemi su obično simulirani koristeći model perli koji pripisuje dve to četiri čestice svakoj aminokiselini.
Funkciona forma
Sledeća funkciona abstrakcija, koja se naziva funkcija potencijala ili polje sila u hemiji, daje potencijalnu energiju molekularnog sistema (E) u datoj konformaciji u obluku sume individualnih energijskih članova.
gde su komponente kovalentnih i nekovalentnih doprinosa date sledećim zbirom:
Precizna funkciona forma potencijalne funkcije, ili polja sila, zavisi od datog simulacionog programa. Generalno članovi veza i uglova se modeluju kao harmonijski potencijali centrirani oko ravnotežnih vrednosti veza i uglova na bazi eksperimenata ili teoretskih proračuna elektronske strukture izvedenih softverom koji vrši ab-initio tipove proračuna, kao što je Gaussian. Za preciznu reprodukciju vibracionog spektra, može se koristiti Morzov potencijal, uz znatan trošak izračunavanja. Diedralni ili torzioni članovi tipično imaju višestruke minimume i stoga se ne mogu modelovati kao harmonijski oscilatori, mada njihova specifična funkciona forma varira u zavisnosti od implementacije. Ova klasa članova može da sadrži neispravne diedralne članove, koji funkcionišu kao korekcioni faktori za devijacije izvan ravni (na primer, oni se mogu koristiti za održavanje planarnosti benzenovih prstena, ili za korekciju geometrije i hiralnosti tetraedralnih atoma u ujedinjeno-atomskoj reprezentaciji).
Članovi koji se ne odnose direktno na veze su daleko skuplji u pogledu potpunog proračuna, pošto je tipični atom vezan za samo nekoliko svojih suseda, a formira interakcije sa svakim drugim atomom molekula. Srećom van der Valsov član brzo opada. On se tipično izračunava koristeći 6–12 Lenard-Džounsov potencijal, što znači da privlačne sile opadaju sa rastojanjem kao r−6, a repulzivne kao r−12, gde r predstavlja rastojanje između atoma. Repulzivni deo r−12 međutim nije fizički, pošto se repulzija eksponencijalno povećava. Opisivanje van der Valsovih sila Lenard-Džounsovim 6–12 potencijalom unosi nepreciznosti, što postaje znatno na kratkim rastojanjima.[1] Generalno se koristi prekidni radijus da bi se ubrzao proračun, tako da atomski parovi čija su rastojanja veća od tog limita imaju energiju van der Valsove interakcije jednaku nuli.
Elektrostatičke članove je veoma teško precizno proračunati pošto oni ne opadaju brzo sa rastojanjem, i elektrostatičke interakcije dugog opsega su često važna osobina studiranog sistema (posebno za proteine). Osnovina funkcionalna forma je Kulonov potencijal, koji opada kao r−1. Različite metode se koriste za rešavanje ovog problema, najjednostavnija je prekidni radijus, slično van der Valsovim članovima. Međutim, time se uvodi oštra diskontinualnost između atoma unutar i izvan radijusa. Prekidačke ili skalirajuće funkcije koje moduliraju pojavnu elektrostatičku energiju su donekle precizniji metodi. Njima se množi izračunata energija glatko varirajućim skalnim faktorom u opsegu od 0 do 1 na spoljašnjim i unutrašnjim prekidnim radijusima. Drugi sofistikovaniji, ali računski intenzivniji metodi su Evaldove mreže čestica (PME) i multipolni algoritam.
Pored funkcionalnog oblika svakog energetskog člana, korisna funkcija energije mora da sadrži parametre za konstante sile, van der Valsove multiplikatore, i druge konstantne članove. Ti članovi, zajedno sa ravnotežnim vrednostima veza, uglova i diedrala, vrednostima parcijalnih naboja, atomskih masa i radijusa, i definicijama funkcija energije, se kolektivno nazivaju poljem sila. Parametrizacija se tipično vrši putem usklađivanja sa eksperimentalnim vrednostima i rezultatima teoretskih proračuna.
Svako polje sila je parameterizovano tako da je unutrašnji konzistentno, mada parameteri generalno nisu transferabilni iz jednog polja sila u drugo.
Oblasti primene
Molekulska mehanika se prvenstveno koristi u polju molekularne dinamike, gde se polje sila koristi za proračun sila koje deluju na svaku česticu, a podesni integrator za modelovanje dinamike čestica i predviđanje trajektorija. Ukoliko su simulacije dovoljno duge i ako je zadovoljena ergodička hipoteza, molekularno dinamičke trajektorije, mogu se koristiti za procenu hermodinamičkih parametara sistema ili za predviđanje kinetičkih svojstava, kao što su reakcione brzine i mehanizmi.
Molekulska mehanika se isto tako koristi za minimizaciju energije, pri čemu se polje sila koristi kao optimizacioni kriterijum. Ovaj metod koristi podesni algoritam (e.g. algoritam opadajućeg gradijenta) za nalaženje molekularne strukture u lokalnom energetskom minimumu. Minimumi odgovaraju stabilnim konformerima molekula (u izabranom polju sila) i molekularno kretanje se može modelovati u vidu vibracija oko i interkonverzija između tih stabilnih konformera. Stoga se često sreću metode lokalne minimizacije energije u kombinaciji sa globalnom optimizacijom energije, koje se koriste u nalaženju globalnog energijskog minimuma (i drugih stanja niske energije). Molekul provodi najveći deo vremena u tim niskim energetskim stanjima, koja stoga dominiraju molekularne osobine. Globalna optimizacija se može ostvariti koristeći algoritam simuliranog žarenja, Metropolisov algoritam i druge Monte Karlo metode, ili koristeći različite determinističke metode diskretne ili kontinuirane optimizacije. Dok polje sila predstavlja samo entalpijsku komponentu Gibsove slobodne energije (i samo ta komponenta se koristi tokom optimizacije energije), moguće je inkorporirati entropsku komponentu koristeći druge metode, kao što je analiza normalnog moda.
U molekulskoj mehanici postoji nekoliko načina da se definiše okruženje molekula i molekul od interesa. Sistem se može simulirati u vakuumu (tzv. simulacija u gasnoj fazi) bez okruženja, mada je to obično nepoželjno pošto uvodi artefakte u molekularnu geometriju, posebno kod naelektrisanih molekula. Površinska naelektrisanja koja bi normalno formirala interakcije sa molekulima rastvarača umesto toga formiraju interakcije jedno s drugim, proizvodeći molekularne konformacije koje se ne javljaju u bilo kom drugom okruženju. Najbolji način solvatacije sistema je da se stave eksplicitni molekuli vode u simulacionu kutiju sa studiranim molekulima i da se molekuli vode tretiraju kao čestice koje formiraju interakcije. Postoji mnoštvo modela vode sa različitim nivoima kompleksnosti, koji predstavljaju vodu kao jednostavnu tvrdu sferu (ujedinjeni atomski model), kao tri zasebne čestice sa fiksnim uglovima veza, ili čak sa četiri ili pet zasebnih interakcionih centara da bi se uzeli u obzir nespareni elektroni na atomu kiseonika. Što je model vode kompleksniji, to su simulacije računarski intenzivnije. Kompromisni metod je implicitna solvatacija, koji zamenjuje eksplicitno predstavljene molekule vode sa matematičkim izrazom koji reprodukuje prosečno ponašanje molekula vode (ili drugih solvenata kao što su lipidi). Ovaj method je koristan za sprečavanje artifakata koji se javljaju pri simulacijama u vakuumu i reprodukuje svojstva rastvarača, ali ne može da reprodukuje situacije u kojima pojedinačni molekuli vode imaju interesantne interakcije sa studiranim molekulima.
^Zgarbova M, et al. (2010). „Large-scale compensation of errors in pairwise-additive empirical force fields: comparison of AMBER intermolecular terms with rigorous DFT-SAPT calculations”. Phys. Chem. Chem. Phys. 12 (35): 10476—10493. Bibcode:2010PCCP...1210476Z. doi:10.1039/C002656E.
^Kuhn B, Kollman PA (2000). „Binding of a diverse set of ligands to avidin and streptavidin: an accurate quantitative prediction of their relative affinities by a combination of Molekulska mehanika and continuum solvent models”. Journal of Medicinal Chemistry. 43 (20): 3786—91. PMID11020294. doi:10.1021/jm000241h.
^Huo S, Massova I, Kollman PA (2002). „Computational alanine scanning of the 1:1 human growth hormone-receptor complex”. J Comput Chem. 23 (1): 15—27. PMID11913381. doi:10.1002/jcc.1153.
^Wang J, Kang X, Kuntz ID, Kollman PA (2005). „Hierarchical database screenings for HIV-1 reverse transcriptase using a pharmacophore model, rigid docking, solvation docking, and MM-PB/SA”. Journal of Medicinal Chemistry. 48 (7): 2432—44. PMID15801834. doi:10.1021/jm049606e.
^Kollman PA, Massova I, Reyes C, et al. (2000). „Calculating structures and free energies of complex molecules: combining molekulska mehanika and continuum models”. Acc Chem Res. 33 (12): 889—97. PMID11123888. doi:10.1021/ar000033j.
Mackerell AD (2004). „Empirical force fields for biological macromolecules: overview and issues”. J Comput Chem. 25 (13): 1584—604. PMID15264253. doi:10.1002/jcc.20082.