Abstract
Quantifying the effects of point mutations is of utmost interest for pharmaceutical and biotechnological applications. Reliable computational methods range from statistical and AI-based to physics-based approaches, with the optimal balance between accurate and fast predictions remaining a challenge. Free energy perturbation (FEP) simulations, a powerful physics-based approach available for decades, constitutes nowadays a method of common application in protein mutational studies. We present QresFEP-2, a novel hybrid-topology FEP protocol benchmarked on a comprehensive protein stability dataset of 10 protein systems, encompassing almost 600 mutations. QresFEP-2 combines excellent accuracy with the highest computational efficiency among available FEP protocols, and its robustness is further validated through comprehensive domain-wide mutagenesis, assessing the thermodynamic stability of over 400 mutations generated by a systematic mutation scan of the 56-residue B1 domain of streptococcal protein G (Gβ1). We also demonstrate the applicability domain of QresFEP-2 on evaluating site-directed mutagenesis effects on protein-ligand binding, tested on a GPCR, as well as on protein-protein interactions examined on the barnase/barstar complex. QresFEP-2 emerges as an open-source, physics-based alternative for advancing protein engineering, drug design, and elucidating the impact of mutations on human health.