Driving Structural Transitions in Molecular Simulations Using the Nonequilibrium Candidate Monte Carlo

Area of Science:

• Computational Chemistry
• Molecular Modeling
• Biophysics

Background:

• Hybrid simulation methods, such as combining molecular dynamics (MD) and Monte Carlo (MC), are gaining traction for enhancing sampling efficiency in molecular simulations.
• Nonequilibrium candidate protocols, where MC moves are gradually applied and interleaved with equilibration, have shown promising results for improving simulation speed.

Purpose of the Study:

• To systematically investigate the efficiency of the nonequilibrium candidate Monte Carlo (NCMC) method for sampling structural transitions in macromolecules.
• To evaluate NCMC's performance on atomistic simulations of peptide systems, focusing on rotameric distributions in both vacuum and explicit solvent environments.

Main Methods:

• Atomistic molecular simulations were performed on two peptide systems.
• The nonequilibrium candidate Monte Carlo (NCMC) procedure was applied to sample rotameric distributions.
• Simulations were conducted in both vacuum and explicit solvent to assess environmental effects.

Main Results:

• NCMC demonstrated the potential for order-of-magnitude gains in sampling efficiency for certain scenarios involving slow degrees of freedom.
• The study identified specific challenges associated with applying NCMC in explicit solvent simulations.
• The sensitivity of NCMC efficiency to its free parameters was explored.

Conclusions:

• Nonequilibrium candidate Monte Carlo is a promising technique for enhancing molecular simulation sampling efficiency, particularly for systems with well-separated free-energy minima.
• Practical application in explicit solvent requires careful consideration of parameter choices and potential challenges.
• Further research is needed to optimize NCMC for complex macromolecular systems and solvent environments.