Stochastic protein folding simulation in the d-dimensional HP-model

  • Authors:
  • K. Steinhöfel;A. Skaliotis;A. A. Albrecht

  • Affiliations:
  • Department of Computer Science, King's College London, Strand, London, UK;Department of Computer Science, King's College London, Strand, London, UK;School of Computer Science, University of Hertfordshire, Hatfield, Herts, UK

  • Venue:
  • BIRD'07 Proceedings of the 1st international conference on Bioinformatics research and development
  • Year:
  • 2007

Quantified Score

Hi-index 0.01

Visualization

Abstract

We present results from two- and three-dimensional protein folding simulations in the HP-model on selected benchmark problems. The importance of the HP-model for investigating general complexity issues of protein folding has been recently demonstrated by Fu & Wang (LNCS 3142:630-644, 2004) in proving an exp(O(n1-1/d ċ ln n)) time bound for d-dimensional protein folding simulation of sequences of length n. The time bound is close to the approximation of real folding times of exp(λ ċ n2/3 ± Χ ċ n1/2/2)ns by Finkelstein & Badretdinov (FOLD DES 2:115-121, 1997), where λ and Χ are constants close to unity. We utilise a stochastic local search procedure that is based on logarithmic simulated annealing. We obtain that after (m/δ)a.D Markov chain transitions the probability to be in a minimum energy conformation is at least 1 - δ, where m ≤ b(d) ċ n is the maximum neighbourhood size for a small integer b(d), a is a small constant, and D is the maximum value of the minimum escape height from local minima of the underlying energy landscape. We note that the time bound is sequence-specific, and we conjecture D n1 -1/d as a worst case upper bound. We analyse D n1-1/d experimentally on selected HP-model benchmark problems.