Abstract:Two improved versions of the pruned-enriched-Rosenbluth method (PERM) are proposed and tested on simple models of lattice heteropolymers. Both are found to outperform not only the previous version of PERM, but also all other stochastic algorithms which have been employed on this problem, except for the core directed chain growth method (CG) of Beutler & Dill. In nearly all test cases they are faster in finding low-energy states, and in many cases they found new lowest energy states missed in previous papers. The CG method is superior to our method in some cases, but less efficient in others. On the other hand, the CG method uses heavily heuristics based on presumptions about the hydrophobic core and does not give thermodynamic properties, while the present method is a fully blind general purpose algorithm giving correct Boltzmann-Gibbs weights, and can be applied in principle to any stochastic sampling problem.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop more effective algorithms to simulate and analyze the growth process of lattice heteropolymers (especially the protein - folding model) at low temperatures. Specifically, the researchers proposed two improved versions of the Pruned - Enriched - Rosenbluth Method (PERM) to more efficiently find low - energy states and provide a complete thermodynamic description.
### Core problems of the paper
1. **Finding low - energy states**: The protein - folding problem is a core problem in computational biology, that is, predicting the lowest - energy state (or native state) of a protein from its amino - acid sequence. Due to the progress of DNA analysis techniques, the number of known amino - acid sequences has increased dramatically, but understanding their three - dimensional structures and functions still lags behind, mainly because of the complexity of the folding problem.
2. **Improving algorithm efficiency**: Traditional Monte Carlo methods and other stochastic algorithms have limitations when dealing with complex energy landscapes, especially under low - temperature conditions. Therefore, new algorithms need to be developed to more efficiently explore these low - energy states and ensure the generality and accuracy of the algorithms.
### Research background
To simplify the description of proteins, lattice models were used in the study, replacing each amino acid with a simple point particle on a regular lattice. Although this simplification reduces complexity, it is still an example of combinatorial optimization and helps to understand complex energy landscapes.
### Models and methods
The models used in the paper mainly include:
- **HP model**: Only two types of monomers are considered, H (hydrophobic) and P (hydrophilic). Hydrophobic monomers tend to avoid water and achieve this by attracting each other.
- **Energy parameters**: Usually, \( \epsilon_{HH}=- 1\), \( \epsilon_{HP}=\epsilon_{PP} = 0\) are used.
### Improved PERM algorithms
The researchers proposed two new PERM variants:
- **nPERMss** (New PERM simple sampling): Simple random uniform sampling is used when selecting branches.
- **nPERMis** (New PERM importance sampling): Importance sampling is used when selecting branches, taking into account the number of free neighbors and energy of different branches.
### Results and advantages
- **Efficiency improvement**: The new algorithms found low - energy states faster than previous PERM versions and other stochastic algorithms in almost all test cases.
- **Discovering new low - energy states**: In some cases, the new algorithms discovered lower - energy states that had not been reported before.
- **Generality and thermodynamic description**: The new algorithms are general and blind, can give the correct Boltzmann - Gibbs weights, and are suitable for any random sampling problem.
### Summary
The main objective of this paper is to improve the efficiency of finding low - energy states of lattice heteropolymers at low temperatures by improving the PERM algorithm and to provide a complete thermodynamic description. This not only helps to understand the protein - folding problem, but also provides effective methods for other similar combinatorial optimization problems.