Solvation parameters for predicting the structure of surface loops in proteins: transferability and entropic effects

Abstract

A new procedure for optimizing parameters of implicit solvation models introduced by us has been applied successfully first to cyclic peptides and more recently to three surface loops of ribonuclease A (Das and Meirovitch, Proteins 2001;43:303-314) using the simplified model E(tot) = E(FF)(epsilon = nr) + Sigma(i) sigma(i)A(i), where sigma(i) are atomic solvation parameters (ASPs) to be optimized, A(i) is the solvent accessible surface area of atom i, E(FF)(epsilon = nr) is the AMBER force-field energy of the loop-loop and loop-template interactions with a distance-dependent dielectric constant, epsilon = nr, where n is a parameter. The loop is free to move while the protein template is held fixed in its X-ray structure; an extensive conformational search for energy minimized loop structures is carried out with our local torsional deformation method. The optimal ASPs and n are those for which the structure with the lowest minimized energy [E(tot)(n,sigma(i))] becomes the experimental X-ray structure, or less strictly, the energy gap between these structures is within 2-3 kcal/mol. To check if a set of ASPs can be defined, which is transferable to a large number of loops, we optimize individual sets of ASPs (based on n = 2) for 12 surface loops from which an "averaged" best-fit set is defined. This set is then applied to the 12 loops and an independent "test" group of 8 loops leading in most cases to very small RMSD values; thus, this set can be useful for structure prediction of loops in homology modeling. For three loops we also calculate the free energy gaps to find that they are only slightly smaller than their energy counterparts, indicating that only larger n will enable reducing too large gaps. Because of its simplicity, this model allowed carrying out an extensive application of our methodology, providing thereby a large number of benchmark results for comparison with future calculations based on n > 2 as well as on more sophisticated solvation models with as yet unknown performance for loops.