Easily searched protein folding potentials

In order to calculate the tertiary structure of a protein from its amino acid sequence, the thermodynamic approach requires a potential function of sequence and conformation that has its global minimum at the native conformation for many different proteins. Here we study the behavior of such functions for the simplest model system that still has the essential features of the protein folding problem, namely two-dimensional square lattice chain configurations involving two residue types. First we demonstrate a method for accurately recovering the given contact potential from only a knowledge of which sequences fold to which structures and what the non-native structures are. Second, we show how to derive from the same information more general potential functions having much better positive correlations between potential function value and conformational deviation from the native. These functions consequently permit faster and more reliable searches for the native conformation, given the native sequence. Furthermore, the method for finding such potentials is easily applied to more realistic protein models.