Run-time Estimates for Protein Folding Simulation--in the H-P Model

Peer reviewe

Run-time Estimates for Protein Folding Simulation in the H-P Model

The hydrophobic-hydrophilic (H-P) model for protein folding was introduced by Dill et al. [6]. A problem instance consists of a sequence of amino acids, each labeled as either hydrophobic (H) or hydrophilic (P). The sequence must be placed on a 2D or 3D grid without overlapping, so that adjacent amino acids in the sequence remain adjacent in the grid. The goal is to minimize the energy, which in the simplest variation corresponds to maximizing the number of adjacent hydrophobic pairs. Although the model is extremely simple, it captures the main features of the problem. The protein folding problem in the H-P model is NP-hard in both 2D and 3D. Recently, Fu and Wang [9] proved an exp(O(n ln n) algorithm for d-dimensional protein folding simulation in the HP-model. Our preliminary results on stochastic search applied to protein folding utilize complete move sets proposed by Lesh et al. [15] and Blazewicz et al. [3]. We obtain that after (n/# ) Markov chain transitions, the probability to be in a minimum energy conformation is at least 1 where n is the length of the instance, # is the maximum value of the minimum escape height from local minima of the underlying energy landscape, and c is a (small) constant