Influence of side chain conformations on local conformational features of amino acids and implication for force field development.

Statistical analysis of coil regions in protein structures has been used to obtain the local backbone phi, psi preferences of amino acids, which agree well with the NMR experiments of unfolded peptides and proteins. We analyzed the conformational features of amino acid residues in a restricted coil library of 4220 high-resolution protein crystal structures. In addition to Gly, Ala, and Pro, the phi, psi distribution (Ramachandran plot) of each amino acid is analyzed with respect to three side chain conformers: g+ (chi(1) approximately -60 degrees), g- (chi(1) approximately 60 degrees), and t (chi(1) approximately 180 degrees). The statistical study indicates that the effect of side chain conformations on phi, psi distributions is even greater than the effect of amino acid types. On the basis of the chi(1), phi, psi conformational preferences, the amino acids in addition to Gly, Pro, and Ala can be divided into five types: (1) ordinary amino acids, (2) Ser, (3) Asp and Asn, (4) Val and Ile, and (5) Thr, each with distinguished chi(1) rotamers. The alpha-helix, beta-sheet, and type-I beta-turn preferences of the different rotamers of various amino acid types can be captured by their intrinsic phi, psi preferences from our coil library. Molecular dynamics simulations of dipeptide Ac-X-NHMe and tetrapeptide Ac-A-X-A-NHMe models give nearly the same side chain rotamer distributions. However, for many amino acids, both OPLS-AA/L and AMBER-FF03 force fields give very different chi(1) rotamer distributions from the coil library. This may partially explain why dipeptide models sometimes cannot reproduce those of protein structures well. The current coil library analysis may be valuable in improving the force field for protein simulations.