Optimized parallel tempering simulations of proteins.

We apply a recently developed adaptive algorithm that systematically improves the efficiency of parallel tempering or replica exchange methods in the numerical simulation of small proteins. Feedback iterations allow us to identify an optimal set of temperatures/replicas which are found to concentrate at the bottlenecks of the simulations. A measure of convergence for the equilibration of the parallel tempering algorithm is discussed. We test our algorithm by simulating the 36-residue villin headpiece subdomain HP-36 where we find a lowest-energy configuration with a root-mean-square deviation of less than 4 A to the experimentally determined structure.

[1]  H. Scheraga,et al.  Intermolecular potentials from crystal data. 6. Determination of empirical potentials for O-H...O = C hydrogen bonds from packing configurations , 1984 .

[2]  H. Scheraga,et al.  Accessible surface areas as a measure of the thermodynamic parameters of hydration of peptides. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Chris Sander,et al.  The double cubic lattice method: Efficient approaches to numerical integration of surface area and volume and to dot surface contouring of molecular assemblies , 1995, J. Comput. Chem..

[4]  P S Kim,et al.  A thermostable 35-residue subdomain within villin headpiece. , 1996, Journal of molecular biology.

[5]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[6]  U. Hansmann.  Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[7]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[8]  U H Hansmann,et al.  New Monte Carlo algorithms for protein folding. , 1999, Current opinion in structural biology.

[9]  Ulrich H. E. Hansmann,et al.  SMMP) A modern package for simulation of proteins , 2001 .

[10]  David A. Kofke,et al.  ARTICLES On the acceptance probability of replica-exchange Monte Carlo trials , 2002 .

[11]  M. Shirts,et al.  Native-like mean structure in the unfolded ensemble of small proteins. , 2002, Journal of molecular biology.

[12]  Ulrich H E Hansmann,et al.  Global optimization by energy landscape paving. , 2002, Physical review letters.

[13]  Karl F Freed,et al.  All‐atom fast protein folding simulations: The villin headpiece , 2002, Proteins.

[14]  Ulrich H. E. Hansmann,et al.  Protein folding in silico: an overview , 2003, Comput. Sci. Eng..

[15]  Ulrich H E Hansmann,et al.  Parallel tempering simulations of HP‐36 , 2003, Proteins.

[16]  M. Troyer,et al.  Performance limitations of flat-histogram methods. , 2003, Physical review letters.

[17]  Harold A Scheraga,et al.  Folding of the villin headpiece subdomain from random structures. Analysis of the charge distribution as a function of pH. , 2004, Journal of molecular biology.

[18]  D. Huse,et al.  Optimizing the ensemble for equilibration in broad-histogram Monte Carlo simulations. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  C. Predescu,et al.  The incomplete beta function law for parallel tempering sampling of classical canonical systems. , 2003, The Journal of chemical physics.

[20]  Simulations of a small protein in a specifically designed generalized ensemble. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  D. Kofke,et al.  Selection of temperature intervals for parallel-tempering simulations. , 2005, The Journal of chemical physics.

[22]  J. D. de Pablo,et al.  Optimal allocation of replicas in parallel tempering simulations. , 2005, The Journal of chemical physics.

[23]  C. Predescu,et al.  On the efficiency of exchange in parallel tempering monte carlo simulations. , 2004, The journal of physical chemistry. B.

[24]  Ulrich H E Hansmann,et al.  Efficient sampling of protein structures by model hopping. , 2005, Physical review letters.