Comparison of protein models minimized by the all-atom and united-atom models in the AMBER force field: correlation of RMS deviation with the crystallographic R factor and size.

We performed energy minimization of 25 protein structures, which vary significantly in their size, secondary structural content and crystallographic R factor, in the AMBER force field. We used an unconstrained path and the conjugate gradients algorithm. To determine the reliability of the united-atom approximation, we minimized all the proteins using both the all-atom and united-atom models. The RMS deviations of the minimized structures were plotted as a function of the crystallographic R factors of the initial structures. For the all-atom models, we found a strong linear relationship between the RMS deviations and the R factors (correlation coefficient of 0.78). The RMS deviations of protein structures minimized using united-atom models showed a wider range of distribution and had a correlation coefficient with the R factors of only 0.52. The RMS deviations decrease with an increase in the size of the protein, probably due to the decreased ratio of surface area to volume with increasing size of the protein. The surface atoms and residues showed higher RMS deviations than those in the interior of the protein. Even in these plots the united-atom models show a wide range of distribution of data points. From these results, we recommend the use of all-atom models for energy minimization of proteins in the AMBER force field.