Protein Engineering, Design and Selection vol. 17 no. 1 pp. 29-36, 2004
© 2004 Oxford University Press
Prediction of the structures of proteins with the UNRES force field, including dynamic formation and breaking of disulfide bonds
aw Odziej1,2
1Baker Laboratory of Chemistry, Cornell University, Ithaca, NY 14853-1301, USA,
2Faculty of Chemistry, University of Gda
sk, ul. Sobieskiego 18, 80-952 Gdask and
3Academic Computer Center in Gdask TASK, ul. Narutowicza 11/12, 80-952 Gdansk, Poland
4 To whom correspondence should be addressed. e-mail: has5{at}cornell.edu
The presence of disulfide bonds is essential for maintaining the structure and function of many proteins. The disulfide bonds are usually formed dynamically during folding. This process is not accounted for in present algorithms for protein-structure prediction, which either deduce the possible positions of disulfide bonds only after the structure is formed or assume fixed disulfide bonds during the course of simulated folding. In this work, the conformational space annealing (CSA) method and the UNRES united-residue force field were extended to treat dynamic formation of disulfide bonds. A harmonic potential is imposed on the distance between disulfide-bonded cysteine side-chain centroids to describe the energetics of bond distortion and an energy gain of 5.5 kcal/mol is added for disulfide-bond formation. Formation, breaking and rearrangement of disulfide bonds are included in the CSA search by introducing appropriate operations; the search can also be carried out with a fixed disulfide-bond arrangement. The algorithm was applied to four proteins: 1EI0 (
), 1NKL (), 1L1I (ß-helix) and 1ED0 ( + ß). For 1EI0, a low-energy structure with correct fold was obtained both in the runs without and with disulfide bonds; however, it was obtained as the lowest in energy only with the native disulfide-bond arrangement. For the other proteins studied, structures with the correct fold were obtained as the lowest (1NKL and 1L1I) or low-energy structures (1ED0) only in runs with disulfide bonds, although the final disulfide-bond arrangement was non-native. The results demonstrate that, by including the possibility of formation of disulfide bonds, the predictive power of the UNRES force field is enhanced, even though the disulfide-bond potential introduced here rarely produces disulfide bonds in native positions. To the best of our knowledge, this is the first algorithm for energy-based prediction of the structure of disulfide-bonded proteins without any assumption as to the positions of native disulfides or human intervention. Directions for improving the potentials and the search method are suggested.
Accepted October 16, 2003 Edited by Valerie Daggett
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  N. Brouwers, K. Sleegers, S. Engelborghs, S. Maurer-Stroh, I. Gijselinck, J. van der Zee, B. A. Pickut, M. Van den Broeck, M. Mattheijssens, K. Peeters, et al. Genetic variability in progranulin contributes to risk for clinically diagnosed Alzheimer disease Neurology, August 26, 2008; 71(9): 656 - 664. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
K. Sleegers, N. Brouwers, S. Maurer-Stroh, M. A. van Es, P. V. Damme, P.W.J. van Vught, J. van der Zee, S. Serneels, T. D. Pooter, M. Van den Broeck, et al. Progranulin genetic variability contributes to amyotrophic lateral sclerosis Neurology, July 22, 2008; 71(4): 253 - 259. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. Oldziej, C. Czaplewski, A. Liwo, M. Chinchio, M. Nanias, J. A. Vila, M. Khalili, Y. A. Arnautova, A. Jagielska, M. Makowski, et al. Physics-based protein-structure prediction using a hierarchical protocol based on the UNRES force field: Assessment in two blind tests PNAS, May 24, 2005; 102(21): 7547 - 7552. [Abstract] [Full Text] [PDF]
|
|