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 sub-domain HP-36 wherewe find a lowest-energy configuration with a root-mean-square-deviation of less than 4 Angstroem to the experimentally determined structure.