Protein Language Model-Based Fitness Estimates Facilitate Resistance Mutation Identification

Drug resistance is a major challenge in cancer therapy. Cancer cells with pre-existing or acquired mutations that confer resistance to a given drug treatment outgrow the susceptible cell population and cause cancer recurrence after an initial successful treatment response. Knowledge about resistance mutations before they occur in the clinic could prevent unnecessary patient treatment with ineffective drugs, in clinical trials as well as clinical practice, or potentially speed up the development of follow-up compounds. Here, we focused on on-target amino acid mutations that confer resistance to an inhibitor compound with a known binding mode. We evaluated whether a combination of physics-based free energy perturbation (FEP) affinity estimates and protein language model-based protein fitness estimates could improve the in silico identification of resistance mutations. Validation was done with data from deep mutational scanning (DMS) experiments that tested for resistance to single amino acid mutations. A public data set testing ERK2 resistance against the inhibitor SCH772984 and an internal data set testing resistance of an EGFR_exon20 mutant against a Bayer small molecule inhibitor were used. Our results show that protein fitness estimates can facilitate the identification of resistance mutations by filtering mutations with a low estimated fitness. Even though FEP has flagged such mutations as affinity-decreasing and thus potentially resistant, they were not resistant according to the DMS experiment and therefore correctly filtered out. This indicates that protein language model-based protein fitness estimates could be a computationally efficient method to filter mutations without having to model the negative impact of mutations on native function or protein stability, which is error-prone and computationally expensive.