A Dynamics-Informed Machine Learning Framework for Inhibitor Optimization Against Quickly Evolving Targets to Avoid Resistance.

Drug resistance often emerges through mutations that reshape protein conformational ensembles, allowing inhibitor potency to be altered without directly disrupting protein-ligand contacts. Consequently, the molecular mechanisms underlying these effects remain difficult to identify, hindering the design of inhibitors with durable potency against rapidly evolving drug targets. To address this challenge, we developed ROBUST, a physics-informed machine learning framework that integrates experimental measurements, molecular dynamics simulations, and interpretable statistical modeling to identify energetic determinants of inhibitor potency. Using HIV-1 protease as a stringent model system for drug resistance, we determined 157 new inhibition constants and integrated them with prior measurements to establish a comprehensive experimental dataset of 364 protease-inhibitor pairs. Analysis of the corresponding molecular dynamics simulations revealed a compact network of energetic interactions that quantitatively explains how resistance mutations and inhibitor modifications jointly reshape molecular recognition and inhibitor potency across the multidimensional resistance landscape. Models constructed from these energetic signatures predicted binding free energies with high accuracy (RMSE < 1 kcal mol