Short-range conformational energies, secondary structure propensities, and recognition of correct sequence-structure matches

Abstract

A statistical analysis of known structures is made for an assessment of the utility of short-range energy considerations. For each type of amino acid, the potentials governing (1) the torsions and bond angle changes of virtual C alpha-C alpha bonds and (2) the coupling between torsion and bond angle changes are derived. These contribute approximately -2 RT per residue to the stability of native proteins, approximately half of which is due to coupling effects. The torsional potentials for the alpha-helical states of different residues are verified to be strongly correlated with the free-energy change measurements made upon single-site mutations at solvent-exposed regions. Likewise, a satisfactory correlation is shown between the beta-sheet potentials of different amino acids and the scales from free-energy measurements, despite the role of tertiary context in stabilizing beta-sheets. Furthermore, there is excellent agreement between our residue-specific potentials for alpha-helical state and other thermodynamic based scales. Threading experiments performed by using an inverse folding protocol show that 50 of 62 test structures correctly recognize their native sequence on the basis of short-range potentials. The performance is improved to 55, upon simultaneous consideration of short-range potentials and the nonbonded interaction potentials between sequentially distant residues. Interactions between near residues along the primary structure, i.e., the local or short-range interactions, are known to be insufficient, alone, for understanding the tertiary structural preferences of proteins alone. Yet, knowledge of short-range conformational potentials permits rationalizing the secondary structure propensities and aids in the discrimination between correct and incorrect tertiary folds.