Development of a new physics‐based internal coordinate mechanics force field and its application to protein loop modeling

YA Arnautova, RA Abagyan… - … : Structure, Function, and …, 2011 - Wiley Online Library
YA Arnautova, RA Abagyan, M Totrov
Proteins: Structure, Function, and Bioinformatics, 2011Wiley Online Library
We report the development of internal coordinate mechanics force field (ICMFF), new force
field parameterized using a combination of experimental data for crystals of small molecules
and quantum mechanics calculations. The main features of ICMFF include:(a)
parameterization for the dielectric constant relevant to the condensed state (ε 5 2) instead of
vacuum,(b) an improved description of hydrogen-bond interactions using duplicate sets of
van der Waals parameters for heavy atom-hydrogen interactions, and (c) improved …
Abstract
We report the development of internal coordinate mechanics force field (ICMFF), new force field parameterized using a combination of experimental data for crystals of small molecules and quantum mechanics calculations. The main features of ICMFF include:(a) parameterization for the dielectric constant relevant to the condensed state (ε 5 2) instead of vacuum,(b) an improved description of hydrogen-bond interactions using duplicate sets of van der Waals parameters for heavy atom-hydrogen interactions, and (c) improved backbone covalent geometry and energetics achieved using novel backbone torsional potentials and inclusion of the bond angles at the C a atoms into the internal variable set. The performance of ICMFF was evaluated through loop modeling simulations for 4–13 residue loops. ICMFF was combined with a solvent-accessible surface area solvation model optimized using a large set of loop decoys. Conformational sampling was carried out using the biased probability Monte Carlo method. Average/median backbone root-mean-square deviations of the lowest energy conformations from the native structures were 0.25/0.21 Afor four residues loops, 0.84/0.46 Afor eight residue loops, and 1.16/0.73 Afor 12 residue loops. To our knowledge, these results are significantly better than or comparable with those reported to date for any loop modeling method that does not take crystal packing into account. Moreover, the accuracy of our method is on par with the best previously reported results obtained considering the crystal environment. We attribute this success to the high accuracy of the new ICM force field achieved by meticulous parameterization, to the optimized solvent model, and the efficiency of the search method.
Wiley Online Library