Abstract:Computing free energy differences between metastable states characterized by non-overlapping Boltzmann distributions is often a computationally intensive endeavour, usually requiring chains of intermediate states to connect these metastable states. Targeted free energy perturbation (TFEP) can significantly lower the computational cost of FEP calculations by choosing a set of invertible maps used to directly transform the distributions of interest, achieving the necessary statistically significant overlaps without sampling any intermediate states. Probabilistic generative models (PGMs) based on normalising-flow architectures can make it much easier via machine learning to train invertible maps needed for TFEP. However, the accuracy and applicability of approaches based on empirically learned maps depend crucially on the choice of reweighting method adopted to estimate the free energy differences. In this work, we assess the accuracy, rate of convergence, and data efficiency of different free energy estimators, including exponential averaging, BAR, and MBAR, in reweighting PGMs trained by maximum likelihood on limited amounts of molecular dynamics data sampled only from end-states of interest. We carry out the comparisons on a set of simple but representative case studies, including conformational ensembles of alanine dipeptide and ibuprofen. Our results indicate that BAR and MBAR are both data efficient and robust, even in the presence of significant model overfitting in the generation of invertible maps. This analysis can serve as a stepping stone for the deployment of efficient and quantitatively accurate ML-based FE calculation methods in complex systems.

Bayesian Multistate Bennett Acceptance Ratio Methods

Bayesian Multistate Bennett Acceptance Ratio Methods

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

Understanding the sources of error in MBAR through asymptotic analysis

Fast Solver for Large Scale Multistate Bennett Acceptance Ratio Equations.

Bayesian Parameter Estimation for Biochemical Reaction Networks Using Region-Based Adaptive Parallel Tempering

DeepBAR: A Fast and Exact Method for Binding Free Energy Computation.

Complex-valued Bayesian parameter estimation via Markov chain Monte Carlo

Principled Approach for Computing Free Energy on Perturbation Graphs with Cycles

Bayesian mixtures of common factor analyzers: Model, variational inference, and applications

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

Multimode process data modeling: A Dirichlet process mixture model based Bayesian robust factor analyzer approach

A Bayesian approach for parameter estimation and prediction using a computationally intensive model

An easy-to-use empirical likelihood ABC method

Bayesian strategies for uncertainty quantification of the thermodynamic properties of materials

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

Adaptive, delayed-acceptance MCMC for targets with expensive likelihoods

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

An adaptive scheduling scheme for calculating Bayes factors with thermodynamic integration using Simpson’s rule

Enhancing MM/P(G)BSA Methods: Integration of Formulaic Entropy for Improved Binding Free Energy Calculations and Virtual Screening Efficiency

Multi-Statistic Approximate Bayesian Computation with Multi-Armed Bandits