TY - JOUR
T1 - Performance of four different force fields for simulations of dipeptide conformations
T2 - GlyGly, GlyGly-, GlyGly·Cl-, GlyGly·Na+and GlyGly·(H2O)2
AU - Dong, Chen
AU - Yong-Zhi, Li
AU - Zhi-Chao, Wei
AU - Bo, Liu
PY - 2014/6
Y1 - 2014/6
N2 - Monte-Carlo conformational searches with four kinds of force fields (AMBER94, MM3*, MMFFs, and OPLS-2005) were performed on glycylglycine (GlyGly), deprotonated glycylglycine (GlyGly-), glycylglycine chloride anion complex (GlyGly·Cl-), glycylglycine sodium cation complex (GlyGly·Na+) and glycylglycine dihydrate [GlyGly·(H2O)2]. Combined with Hartree-Fock (HF) and second-order Møller-Plesset (MP2) optimizations, conformations within an energy of 20 kJ mol-1 were predicted. After MP2 calculations, the geometries and relative energies of the predicted structures were the same regardless of the force field used. Therefore, the performance of different force fields reflects mainly the conformational search process. For GlyGly, there was practically no difference among the four force fields. Due to the complex hydrogen bonding network when involving water, the total number of resulting conformers for GlyGly·(H2O)2 increased drastically. Moreover, the MMFFs force field fared best in finding the global minimum compared to the remaining three force fields. In describing hydrogen bonded and inter-molecular complexes, we recommend application of the MMFFs and AMBER94 force fields. Furthermore, the MMFFs and OPLS-2005 force fields have a good description of electrostatic interactions. This work will contribute to helping the reader make an optimal choice of force field, taking into account the latter's strengths and limitations.
AB - Monte-Carlo conformational searches with four kinds of force fields (AMBER94, MM3*, MMFFs, and OPLS-2005) were performed on glycylglycine (GlyGly), deprotonated glycylglycine (GlyGly-), glycylglycine chloride anion complex (GlyGly·Cl-), glycylglycine sodium cation complex (GlyGly·Na+) and glycylglycine dihydrate [GlyGly·(H2O)2]. Combined with Hartree-Fock (HF) and second-order Møller-Plesset (MP2) optimizations, conformations within an energy of 20 kJ mol-1 were predicted. After MP2 calculations, the geometries and relative energies of the predicted structures were the same regardless of the force field used. Therefore, the performance of different force fields reflects mainly the conformational search process. For GlyGly, there was practically no difference among the four force fields. Due to the complex hydrogen bonding network when involving water, the total number of resulting conformers for GlyGly·(H2O)2 increased drastically. Moreover, the MMFFs force field fared best in finding the global minimum compared to the remaining three force fields. In describing hydrogen bonded and inter-molecular complexes, we recommend application of the MMFFs and AMBER94 force fields. Furthermore, the MMFFs and OPLS-2005 force fields have a good description of electrostatic interactions. This work will contribute to helping the reader make an optimal choice of force field, taking into account the latter's strengths and limitations.
KW - Electrostatic interaction
KW - Force field
KW - Hydrogen bonding interaction
KW - The hybrid MCMM/QM approach
UR - http://www.scopus.com/inward/record.url?scp=84901576197&partnerID=8YFLogxK
U2 - 10.1007/s00894-014-2279-4
DO - 10.1007/s00894-014-2279-4
M3 - Article
C2 - 24863532
AN - SCOPUS:84901576197
SN - 1610-2940
VL - 20
JO - Journal of Molecular Modeling
JF - Journal of Molecular Modeling
IS - 6
M1 - 2279
ER -