  An optimization technique is used to determine the pairwise interactions
between amino acids in globular proteins. A numerical strategy is applied to a
set of proteins for maximizing the native fold stability with respect to
alternative structures obtained by gapless threading. The extracted parameters
are shown to be very reliable for identifying the native states of proteins
(unrelated to those in the training set) among thousands of conformations. The
only poor performers are proteins with heme groups and/or poor compactness
whose complexity cannot be captured by standard pairwise energy functionals.
