Free energy calculations in molecular modeling: from classical methods to machine learning

Free energy calculations have evolved from specialized theoretical tools to essential components of modern molecular design, particularly in pharmaceutical research where alchemical binding free energy methods now guide lead optimization decisions. Despite significant methodological advances over the past two decades, including statistically optimal estimators, enhanced sampling techniques, multi-dimensional nonequilibrium protocols, and improved force fields, current methods achieve reliable accuracy only for favorable systems involving rigid proteins and uncharged ligands. Challenges include inadequate sampling of protein flexibility, systematic biases from classical force fields for charged species, and convergence failures for systems with slow conformational transitions, limiting broader application to challenging but industrially relevant targets.

This review provides a comprehensive analysis of free energy calculation methodologies, from foundational approaches (thermodynamic integration, Bennett acceptance ratio) to advanced enhanced sampling strategies (replica exchange, metadynamics) and cutting-edge developments including machine learning integration. Current accuracy benchmarks are assessed through systematic analysis of community challenges and industrial validation studies, identifying specific failure modes and their underlying physical origins. This review emphasizes practical method selection criteria and realistic accuracy expectations while also examining how machine learning approaches-including neural network potentials, automated collective variable discovery, and active learning protocols-address limitations in sampling efficiency and force field accuracy.