Computational Lead Optimization on BACE1: Relative Binding Free Energy Perturbation as the Terminal Refinement Layer

Abstract

Structure-based drug discovery is known to apply computational methods in a tiered hierarchy, with each layer narrowing the candidate set and refining the binding picture before committing to the next, more expensive step. We present a four-tiered computational benchmarking study evaluating five engines against a panel of 36 compounds targeting β-secretase 1 (BACE1), a validated Alzheimer’s disease target with extensive co-crystal ground truth. This study evaluates Flexible Docking and Boltz2 Cofolding as the primary tier, followed by Ensemble Docking, and then Protein-Ligand MD with MM/PBSA and MM/GBSA post-processing. This is then concluded with Relative Binding Free Energy Perturbation (RevFEP) as the terminal refinement layer. Each method was benchmarked against the experimental binding free energies derived from the co-crystal structures spanning −7.85 to −11.35 kcal/mol. Our findings revealed that Flexible Docking reproduced the co-crystal binding mode for 35 of 36 ligands (97.2% within 2.0 Å RMSD) but did not rank potency at this resolution. Boltz2 CoFolding provided an orthogonal structural cross check with a receptor backbone RMSD of 0.293 Å against the experimental co-crystal structure. Ensemble Docking identified the optimal receptor conformation for downstream FEP setup. MD with MM/GBSA decomposition identified van der Waals complementarity as the primary potency driver (Pearson r = +0.855, R2 = 0.732 on a 10-compound subset). RevFEP delivered the highest affinity correlation of any method (Pearson r = +0.662, R2 = 0.438, Spearman ρ = +0.624, mean absolute error 1.02 kcal/mol across all 36 ligands), resolving potency differences within a narrow 3.5 kcal/mol congeneric window that no other engine could discriminate. We characterize what each engine contributes independently and where RevFEP delivers signals no other engine achieves.