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/.

Reweighting from the mixture distribution as a better way to describe the Multistate Bennett Acceptance Ratio

Bayesian Multistate Bennett Acceptance Ratio Methods

Understanding the sources of error in MBAR through asymptotic analysis

Optimal Reduction of Multivariate Dirac Mixture Densities

An automatic adaptive method to combine summary statistics in approximate Bayesian computation

Fast Solver for Large Scale Multistate Bennett Acceptance Ratio Equations.

A Parzen-based distance between probability measures as an alternative of summary statistics in Approximate Bayesian Computation

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

Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation

An easy-to-use empirical likelihood ABC method

Assessing the accuracy and efficiency of free energy differences obtained from reweighted flow-based probabilistic generative models

Reweighting in nonequilibrium simulations

A minimalistic and general weighted averaging method for inconsistent data

Using Autodiff to Estimate Posterior Moments, Marginals and Samples

A Random Weighting Approach for Posterior Distributions

Adaptive, delayed-acceptance MCMC for targets with expensive likelihoods

Improved Nested Sampling and Surrogate‐Enabled Comparison with Other Marginal Likelihood Estimators

Ensemble Bayesian model averaging using Markov Chain Monte Carlo sampling