Premium
Monte Carlo simulations of biomolecules: The MC module in CHARMM
Author(s) -
Hu Jie,
Ma Ao,
Dinner Aaron R.
Publication year - 2005
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.20327
Subject(s) - monte carlo method , statistical physics , biomolecule , computer science , materials science , mathematics , physics , nanotechnology , statistics
We describe the implementation of a general and flexible Monte Carlo (MC) module for the program CHARMM, which is used widely for modeling biomolecular systems with empirical energy functions. Construction and use of an almost arbitrary move set with only a few commands is made possible by providing several predefined types of moves that can be combined. Sampling can be enhanced by noncanonical acceptance criteria, automatic optimization of step sizes, and energy minimization. A systematic procedure for improving MC move sets is introduced and applied to simulations of two peptides. The resulting move sets allow MC to sample the configuration spaces of these systems much more rapidly than Langevin dynamics. The rate of convergence of the difference in free energy between ethane and methanol in explicit solvent is also examined, and comparable performances are observed for MC and the Nosé–Hoover algorithm. Its ease of use combined with its sampling efficiency make the MC module in CHARMM an attractive alternative for exploring the behavior of biomolecular systems. © 2005 Wiley Periodicals, Inc. J Comput Chem 27: 203–216, 2006