Optimized Wang-Landau sampling of lattice polymers: Ground state search and folding thermodynamics of HP model proteins
arXiv:1207.3974 · doi:10.1063/1.4742969
Abstract
Coarse-grained (lattice-) models have a long tradition in aiding efforts to decipher the physical or biological complexity of proteins. Despite the simplicity of these models, however, numerical simulations are often computationally very demanding and the quest for efficient algorithms is as old as the models themselves. Expanding on our previous work [T. Wüst and D. P. Landau, Phys. Rev. Lett. 102, 178101 (2009)], we present a complete picture of a Monte Carlo method based on Wang-Landau sampling in combination with efficient trial moves (pull, bond-rebridging and pivot moves) which is particularly suited to the study of models such as the hydrophobic-polar (HP) lattice model of protein folding. With this generic and fully blind Monte Carlo procedure, all currently known putative ground states for the most difficult benchmark HP sequences could be found. For most sequences we could also determine the entire energy density of states and, together with suitably designed structural observables, explore the thermodynamics and intricate folding behavior in the virtually inaccessible low-temperature regime. We analyze the differences between random and protein-like heteropolymers for sequence lengths up to 500 residues. Our approach is powerful both in terms of robustness and speed, yet flexible and simple enough for the study of many related problems in protein folding.
15 pages, 7 figures
References in corpus (19)
- An efficient, multiple range random walk algorithm to calculate the density of states
- Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram
- A Portable High-Quality Random Number Generator for Lattice Field Theory Simulations
- Fast Algorithm to Calculate Density of States
- Understanding and Improving the Wang-Landau Algorithm
- Wang-Landau algorithm for continuous models and joint density of states
- Multicanonical Chain Growth Algorithm
- Growth Algorithms for Lattice Heteropolymers at Low Temperatures
- Thermodynamics of Lattice Heteropolymers
- Versatile approach to access the low temperature thermodynamics of lattice polymers and proteins
- Multi-Self-Overlap Ensemble for protein folding: ground state search and thermodynamics
- A New Monte Carlo Algorithm for Protein Folding
- Simulation of Lattice Polymers with Multi-Self-Overlap Ensemble
- A Growth-based Optimization Algorithm for Lattice Heteropolymers
- Optimal Modification Factor and Convergence of the Wang-Landau Algorithm
- Long Range Moves for High Density Polymer Simulations
- Substrate Specificity of Peptide Adsorption: A Model Study
- Monte Carlo simulations of the HP model (the "Ising model" of protein folding)
- Sampling the two-dimensional density of states g(E,M) of a giant magnetic molecule using the Wang-Landau method
Cited by in corpus (13)
- Scalable replica-exchange framework for Wang-Landau sampling
- Sequence determines degree of knottedness in a coarse-grained protein model
- Generic folding and transition hierarchies for surface adsorption of hydrophobic-polar lattice model proteins
- Multicanonical MCMC for Sampling Rare Events
- Control of accuracy in the Wang-Landau algorithm
- On the non-convergence of the Wang-Landau algorithms with multiple random walkers
- Convergence of Stochastic Approximation Monte Carlo and modified Wang-Landau algorithms: Tests for the Ising model
- Predicting Non-Equilibrium Folding Behavior of Polymer Chains using the Steepest-Entropy-Ascent Quantum Thermodynamic Framework
- Random matrix ensembles: Wang-Landau algorithm for spectral densities
- Predicting Polymer Brush Behavior in Solvents using the Steepest-Entropy-Ascent Quantum Thermodynamic Framework
- Spectral properties of unimodular lattice triangulations
- Conformational transitions in random heteropolymer models
- A Parallel Trajectory Swapping Wang - Landau Study Of The HP Protein Model