Multi-Layer Free Energy Perturbation

Abstract: Free energy perturbation (FEP) is frequently used to evaluate the free energy change of a biological process, e.g. the drug binding free energy or the ligand solvation free energy. Due to the sampling inefficiency, FEP is often employed together with computationally expensive enhanced sampling methods. Here we show that this sampling inefficiency, which stems from not accounting for the environmental reorganization, is an inherent property of the single-ensemble ansatz of FEP, and hence simply prolonging the MD simulation can hardly alleviate the problem. Instead, we propose a new, multi-ensemble ansatz -- the multi-layer free energy perturbation (MLFEP), which allows environmental reorganization processes (relaxation) to occur automatically during the MD sampling. Our study paves the way toward a fast but accurate free energy calculation that can be employed in computer-aided drug design.