Abstract:Accelerated molecular dynamics (aMD) simulations greatly improve the efficiency of conventional molecular dynamics (cMD) for sampling biomolecular conformations, but they require proper reweighting for free energy calculation. In this work, we systematically compare the accuracy of different reweighting algorithms including the exponential average, Maclaurin series, and cumulant expansion on three model systems: alanine dipeptide, chignolin, and Trp-cage. Exponential average reweighting can recover the original free energy profiles easily only when the distribution of the boost potential is narrow (e.g., the range ≤20kBT) as found in dihedral-boost aMD simulation of alanine dipeptide. In dual-boost aMD simulations of the studied systems, exponential average generally leads to high energetic fluctuations, largely due to the fact that the Boltzmann reweighting factors are dominated by a very few high boost potential frames. In comparison, reweighting based on Maclaurin series expansion (equivalent to cumulant expansion on the first order) greatly suppresses the energetic noise but often gives incorrect energy minimum positions and significant errors at the energy barriers (∼2-3kBT). Finally, reweighting using cumulant expansion to the second order is able to recover the most accurate free energy profiles within statistical errors of ∼kBT, particularly when the distribution of the boost potential exhibits low anharmonicity (i.e., near-Gaussian distribution), and should be of wide applicability. A toolkit of Python scripts for aMD reweighting "PyReweighting" is distributed free of charge at http://mccammon.ucsd.edu/computing/amdReweighting/.

Probability Density Reweighting of High-Temperature Molecular Dynamics

On the Stochastic Dynamics of Molecular Conformation

Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules

RiteWeight: Randomized Iterative Trajectory Reweighting for Steady-State Distributions Without Discretization Error

Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation

Combining stochastic density functional theory with deep potential molecular dynamics to study warm dense matter

Adaptive Sampling Methods for Molecular Dynamics in the Era of Machine Learning

From data to noise to data for mixing physics across temperatures with generative artificial intelligence

A Comment on the Reweighting Method for Accelerated Molecular Dynamics Simulations

Warm dense matter simulation via electron temperature dependent deep potential molecular dynamics

An Integrate-over-temperature Approach for Enhanced Sampling.

Accurate Prediction for Protein-Peptide Binding Based on High-Temperature Molecular Dynamics Simulations

Equilibrium sampling by re-weighting non-equilibrium simulation trajectories

Coupling Accelerated Molecular Dynamics Methods with Thermodynamic Integration Simulations

A Discard-and-Restart MD algorithm for the sampling of realistic protein transition states and enhance structure-based drug discovery

Simulating rare events using a Weighted Ensemble-based string method

Enhanced sampling techniques in molecular dynamics simulations of biological systems

Lattice Models for Magnetic Fluids: Correlations Between Order Parameters