Spatial Averaging: Sampling Enhancement for Exploring Configurational Space of Atomic Clusters and Biomolecules

Spatial averaging Monte Carlo (SA-MC) is an efficient algorithm dedicated to the study of rare-event problems. At the heart of this method is the realization that from the equilibrium density a related, modified probability density can be constructed through a suitable transformation. This new density is more highly connected than the original density, which increases the probability for transitions between neighboring states, which in turn speeds up the sampling. The first successful investigations included the diffusion of small molecules in condensed phase environments and characterization of the metastable states of the migration of the CO ligand in myoglobin. In the present work, a general and robust implementation including rotational and torsional moves in the CHARMM molecular modeling software is introduced. Also, a procedure to estimate unbiased properties is proposed in order to compute thermodynamic observables. These procedures are suitable to study a range of topical systems including Lennard-Jones clusters of different sizes and the blocked alanine dipeptide (Ala)2 in implicit and explicit solvent. In all cases, SA-MC is found to outperform standard Metropolis simulations in sampling configurational space at little extra computational expense. The results for (Ala)2 in explicit solvent are in good agreement with previous umbrella sampling simulations.