Premium
Hierarchical energy‐based approach to protein‐structure prediction: Blind‐test evaluation with CASP3 targets
Author(s) -
Lee Jooyoung,
Liwo Adam,
Ripoll Daniel R.,
Pillardy Jaroslaw,
Saunders Jeffrey A.,
Gibson Kenneth D.,
Scheraga Harold A.
Publication year - 2000
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/(sici)1097-461x(2000)77:1<90::aid-qua10>3.0.co;2-l
Subject(s) - globular protein , force field (fiction) , protein secondary structure , protein structure prediction , statistical potential , protein design , protein structure , simulated annealing , threading (protein sequence) , monte carlo method , energy landscape , chemistry , physics , crystallography , algorithm , statistical physics , mathematics , thermodynamics , statistics , quantum mechanics , nuclear magnetic resonance
A hierarchical approach based exclusively on finding the global minimum of an appropriate potential energy function, without the aid of secondary structure prediction, multiple‐sequence alignment, or threading, is proposed. The procedure starts from an extensive search of the conformational space of a protein, using our recently developed united‐residue off‐lattice UNRES force field and the conformational space annealing (CSA) method. The structures obtained in the search are clustered into families and ranked according to their UNRES energy. Structures within a preassigned energy cutoff are gradually converted into an all‐atom representation, followed by a limited conformational search at the all‐atom level, using the electrostatically driven Monte Carlo (EDMC) method and the ECEPP/3 force field including hydration. The approach was tested (in the CASP3 experiment) in blind predictions on seven targets, five of which were globular proteins with sizes ranging from 89 to 140 amino acid residues. Comparison of the computed lowest‐energy structures, with the experimental structures, made available after the predictions were submitted, shows that large fragments (∼60 residues, representing 45–80% of the proteins) of those five globular proteins were predicted with the root mean square deviations (RMSDs) ranging from 4 to 7 Å for the C α atoms, with correct secondary structure and topology. These results constitute an important step toward the prediction of protein structure based solely on global optimization of a potential energy function for a given amino acid sequence. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 90–117, 2000