A test of lattice protein folding algorithms

We report a blind test of lattice-model-based search strategies for finding global minima of model protein chains. One of us (E.I.S.) selected 10 compact conformations of 48-mer chains on the three-dimensional cubic lattice and used their inverse folding algorithm to design HP (H, hydrophobic; P, polar) sequences that should fold to those 'target' structures. The sequences, but not the structures, were sent to the UCSF group (K.Y., K.M.F., P.D.T., H.S.C., and K.A.D.), who used two methods to attempt to find the globally optimal conformations: 'hydrophobic zippers' and a constraint-based hydrophobic core construction (CHCC) method. The CHCC method found global minima in all cases, and the hydrophobic zippers method found global minima in some cases, in minutes to hours on workstations. In 9 out of 10 sequences, the CHCC method found lower energy conformations than the 48-mers were designed to fold to. Thus the search strategies succeed for the HP model but the design strategy does not. For every sequence the global energy minimum was found to have multiple degeneracy with 10³ to 10⁶ conformations. We discuss the implications of these results for (i) searching conformational spaces of simple models of proteins and (ii) how these simple models relate to proteins.