Indexed on: 26 Oct '05Published on: 26 Oct '05Published in: Physical review. E, Statistical, nonlinear, and soft matter physics
It is shown in this paper that some simulations of protein folding in lattice models, which use an incorrect implementation of the Monte Carlo algorithm, do not converge towards thermal equilibrium. I developed a rigorous treatment for protein folding simulation on a lattice model relying on the introduction of a parameter standing for the rigidity of the conformations. Its properties are discussed and its role during the folding process is elucidated. The calculation of thermal properties of small chains living on a two-dimensional lattice is performed and a Bortz-Kalos-Lebowitz scheme is implemented in the presented method in order to study kinetics of chains at very low temperature. The coefficients of the Arrhenius law obtained with this algorithm are found to be in excellent agreement with the value of the main potential barrier of the system. Finally, a scenario of the mechanisms, including the rigidity parameters, that guide a protein towards its native structure, at medium temperature, is given.