Abstract
The use of antiviral drugs can promote the appearance of mutations in the target protein that increase the resistance of the virus to the treatment. This is also the case of nirmatrelvir, a covalent inhibitor of the 3CL protease, or main protease, of SARS-CoV-2. In this work we show how the by-residue decomposition of noncovalent interactions established between the drug and the enzyme, in combination with an analysis of naturally occurring mutations, can be used to detect potential mutations in the 3CL protease conferring resistance to nirmatrelvir. We also investigate the consequences of these mutations on the reaction mechanism to form the covalent enzyme-inhibitor complex using QM/MM methods. In particular, we show that the E166V variant of the protease displays smaller binding affinity to nirmatrelvir and larger activation free energy for the formation of the covalent complex, both factors contributing to the observed resistance to the treatment with this drug. The conclusions derived from our work can be used to anticipate the consequences of the introduction of nirmatrelvir in the fitness landscape of the virus and to design new inhibitors adapted to some of the possible resistance mechanisms.