Abstract:Free-energy perturbation calculation is frequently used to calculate free-energy differences because it is easy to implement and the computation is fast. However, the calculation is subject to large inaccuracies in some circumstances due to the insufficient sampling of the relevant tails of the energy-difference distributions. Here we expand this knowledge of insufficient sampling into a two-dimensional (2D) energy space using a model of harmonic oscillators. We show analytically the relation between the energies of the sampling system and those of the desired target energy spaces, which provide the basis to understand the difficulties in free-energy perturbation calculations. We clarify the reasons of the inaccurate calculation in the different harmonic cases that stem from the spatial separations of the reference and the target energy pairs located in the two-dimensional energy space. The potential-energy space introduced into this 2D energy-space model provides additional clues to improve the sampling efficiency. Based on this understanding, we propose two ways to calculate the free-energy differences using the two schemes of the distribution method. We show that the distribution method implemented in the appropriate energy space--the energy-difference space and the potential-energy space, respectively--can improve the calculation of free energies in different circumstances. This analysis implies that the sampling can be improved if it is directed toward the appropriate region in the potential-energy space, which is easily implemented in various types of free-energy calculations. To test this, we calculate the free-energy surface of alanine dipeptide in gas phase and in aqueous phase, respectively. We demonstrate that the free-energy surface calculation is improved when the biased sampling of the potential energy is integrated into the sampling scheme.

Bayesian Approach for Computing Free Energy on Perturbation Graphs with Cycles

Principled Approach for Computing Free Energy on Perturbation Graphs with Cycles

Energy-based Analysis of Biochemical Cycles using Bond Graphs

On oriented cycles in randomly perturbed digraphs

Construction of a Pathway Map on a Complicated Energy Landscape.

Free Energy Landscape for Biological Systems

A Method to Calculate the Number of Spanning Connected Unicyclic(bicyclic) Subgraphs in 2-Separable Networks

Dual-Topology Hamiltonian-Replica-Exchange Overlap Histogramming Method to Calculate Relative Free Energy Difference in Rough Energy Landscape

An Open Source Graph-Based Weighted Cycle Closure Method for Relative Binding Free Energy Calculations

A perturbation approach for refining Boolean models of cell cycle regulation

Convergent Message-Passing Algorithms for Inference over General Graphs with Convex Free Energies

An Improved Approach for Calculating Energy Landscape of Gene Networks from Moment Equations

Understanding Free-Energy Perturbation Calculations Through A Model of Harmonic Oscillators: Theory and Implications to Improve the Sampling Efficiency by Molecular Simulation

Full Conditional Free Energy Based Inference

Unstable Limit Cycles Estimation from Small Perturbations: A Model-Free Approach

Skewed Thermodynamic Geometry and Optimal Free Energy Estimation

Detecting Quantumness in the N-Cycle Exclusivity Graphs

Comparison among Classical, Probabilistic and Quantum Algorithms for Hamiltonian Cycle problem

Counting and correcting thermodynamically infeasible flux cycles in genome-scale metabolic networks

The energy landscape of the Kuramoto model in random geometric graphs in a circle

On the Convexity and Reliability of the Bethe Free Energy Approximation