Monte Carlo Update for Chain Molecules: Biased Gaussian Steps in Torsional Space

Abstract

We develop a new elementary move for simulations of polymer chains in torsion angle space. The method is flexible and easy to implement. Tentative updates are drawn from a (conformation-dependent) Gaussian distribution that favors approximately local deformations of the chain. The degree of bias is controlled by a parameter b. The method is tested on a reduced model protein with 54 amino acids and the Ramachandran torsion angles as its only degrees of freedom, for different b. Without excessive fine tuning, we find that the effective step size can be increased by a factor of three compared to the unbiased b=0 case. The method may be useful for kinetic studies, too.

14 pages, 4 figures

Cited by in corpus (11)

• Oligomerization of amyloid Abeta peptides using hydrogen bonds and hydrophobicity forces
• Folding thermodynamics of peptides
• Thermodynamics of alpha- and beta-structure formation in proteins
• Changing the mechanical unfolding pathway of FnIII10 by tuning the pulling strength
• Advanced multicanonical Monte Carlo methods for efficient simulations of nucleation processes of polymers
• Heat Bath Efficiency with Metropolis-Type Updating
• Mechanical resistance in unstructured proteins
• Two-state folding over a weak free-energy barrier
• Rugged Metropolis Sampling with Simultaneous Updating of Two Dynamical Variables
• Markov modeling of peptide folding in the presence of protein crowders
• Study of a model for the folding of a small protein