Powerful simulated‐annealing algorithm locates global minimum of protein‐folding potentials from multiple starting conformations

Protein‐folding potentials, designed with the explicit goal that the global energy minimum correspond to crystallographically observed conformations of protein molecules, may offer great promise toward calculating native protein structures. Achieving this promise, however, depends on finding an effective means of dealing with the multiple‐minimum problem inherent in such potentials. In this study, a protein‐folding‐potential test system has been developed that exhibits the properties of general protein‐folding potentials yet has a unique well‐defined global energy minimum corresponding to the crystallographically determined conformation of the test molecule. A simulated‐annealing algorithm is developed that locates the global minimum of this potential in four of eight test runs from random starting conformations. Exploration of the energy‐conformation surface of the potential indicates that it contains the numerous local minima typical of protein‐folding potentials and that the global minimum is not easily located by conventional minimization procedures. When the annealing algorithm is applied to a previously developed actual folding potential to analyze the conformation of avian pancreatic polypeptide, a new conformer is located that is lower in energy than any conformer located in previous studies using a variety of minimization techniques.

[1]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[2]  C. M. Reeves,et al.  Function minimization by conjugate gradients , 1964, Comput. J..

[3]  J. Deisenhofer,et al.  Crystallographic refinement of the structure of bovine pancreatic trypsin inhibitor at l.5 Å resolution , 1975 .

[4]  T. Blundell,et al.  Conformational flexibility in a small globular hormone: X‐ray analysis of avian pancreatic polypeptide at 0.98‐Å resolution , 1983, Biopolymers.

[5]  M. Karplus,et al.  An analysis of incorrectly folded protein models. Implications for structure predictions. , 1984, Journal of molecular biology.

[6]  A Wlodawer,et al.  Comparison of two highly refined structures of bovine pancreatic trypsin inhibitor. , 1987, Journal of molecular biology.

[7]  G. Crippen,et al.  Determination of an empirical energy function for protein conformational analysis by energy embedding , 1987 .

[8]  N Go,et al.  Normal modes of vibration in bovine pancreatic trypsin inhibitor and its mechanical property , 1987, Proteins.

[9]  L. Haines The application of the annealing algorithm to the construction of exact optimal designs for linear-regression models , 1987 .

[10]  R. Bruccoleri,et al.  Criteria that discriminate between native proteins and incorrectly folded models , 1988, Proteins.

[11]  Gordon M. Crippen,et al.  Distance Geometry and Molecular Conformation , 1988 .

[12]  R. H. J. M. Otten,et al.  The Annealing Algorithm , 1989 .

[13]  Gordon M. Crippen,et al.  Global energy minimization by rotational energy embedding , 1990, J. Chem. Inf. Comput. Sci..

[14]  G M Crippen,et al.  A 1.8 Å resolution potential function for protein folding , 1990, Biopolymers.

[15]  Dimensional oscillation. A fast variation of energy embedding gives good results with the AMBER potential energy function. , 2009, International journal of peptide and protein research.

[16]  Akbar Nayeem,et al.  A comparative study of the simulated‐annealing and Monte Carlo‐with‐minimization approaches to the minimum‐energy structures of polypeptides: [Met]‐enkephalin , 1991 .

[17]  D. Romero,et al.  Applications of simulated annealing to the multiple-minima problem in small peptides. , 1991, Journal of biomolecular structure & dynamics.

[18]  Stephen R. Wilson,et al.  Applications of simulated annealing to the conformational analysis of flexible molecules , 1991 .

[19]  P. Seetharamulu,et al.  A potential function for protein folding , 1991 .