Constraint-based hydrophobic core construction for protein structure prediction in the face-centered-cubic lattice

We present an algorithm for exact protein structure prediction in the FCC-HP-model. This model is a lattice protein model on the face-centered-cubic lattice that models the main force of protein folding, namely the hydrophobic force. The structure prediction for this model can be based on the construction of hydrophobic cores. The main focus of the paper is on an algorithm for constructing maximally and submaximally compact hydrophobic cores of a given size. This algorithm treats core construction as a constraint satisfaction problem (CSP), and the paper describes its constraint model. The algorithm employs symmetry excluding constraint-based search and relies heavily on good upper bounds on the number of contacts. Here, we use and strengthen upper bounds presented earlier. The resulting structure prediction algorithm (including previous work) handles sequences of sizes in the range of real proteins fast, i.e. we predict a first structure often within a few minutes. The algorithm is the first exact one for the FCC, besides full enumeration which is impracticable for chain lengths greater than about 15. We tested the algorithm successfully up to sequence length of 160, which is far beyond the capabilities even of previous heuristic approaches.