Premium
Strike a balance: Optimization of backbone torsion parameters of AMBER polarizable force field for simulations of proteins and peptides
Author(s) -
Wang ZhiXiang,
Zhang Wei,
Wu Chun,
Lei Hongxing,
Cieplak Piotr,
Duan Yong
Publication year - 2006
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.20386
Subject(s) - force field (fiction) , polarizability , molecular dynamics , torsion (gastropod) , chemistry , computational chemistry , statistical physics , polarization (electrochemistry) , physics , molecular physics , chemical physics , atomic physics , molecule , quantum mechanics , medicine , surgery , organic chemistry
Based on the AMBER polarizable model (ff02), we have reoptimized the parameters related to the main‐chain (Φ, Ψ) torsion angles by fitting to the Boltzmann‐weighted average quantum mechanical (QM) energies of the important regions (i.e., β, P II , α R , and α L regions). Following the naming convention of the AMBER force field series, this release will be called ff02pol.rl The force field has been assessed both by energetic comparison against the QM data and by the replica exchange molecular dynamics simulations of short alanine peptides in water. For Ace‐Ala‐Nme, the simulated populations in the β, P II and α R regions were approximately 30, 43, and 26%, respectively. For Ace‐(Ala) 7 ‐Nme, the populations in these three regions were approximately 24, 49, and 26%. Both were in qualitative agreement with the NMR and CD experimental conclusions. In comparison with the previous force field, ff02pol.rl demonstrated good balance among these three important regions. The optimized torsion parameters, together with those in ff02, allow us to carry out simulations on proteins and peptides with the consideration of polarization. © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 781–790, 2006