Abstract Nirmatrelvir, an oral antiviral targeting the 3CL protease of SARS-CoV-2, has been demonstrated to be clinically useful in reducing hospitalization or death due to COVID-19 1,2 . However, as SARS-CoV-2 has evolved to become resistant to other therapeutic modalities 3–9 , there is a concern that the same could occur for nirmatrelvir. Here, we have examined this possibility by in vitro passaging of SARS-CoV-2 in increasing concentrations of nirmatrelvir using two independent approaches, including one on a large scale in 480 wells. Indeed, highly resistant viruses emerged from both, and their sequences revealed a multitude of 3CL protease mutations. In the experiment done at a larger scale with many replicates, 53 independent viral lineages were selected with mutations observed at 23 different residues of the enzyme. Yet, several common mutational pathways to nirmatrelvir resistance were preferred, with a majority of the viruses descending from T21I, P252L, or T304I as precursor mutations. Construction and analysis of 13 recombinant SARS-CoV-2 clones, each containing a unique mutation or a combination of mutations showed that the above precursor mutations only mediated low-level resistance, whereas greater resistance required accumulation of additional mutations. E166V mutation conferred the strongest resistance (~100-fold), but this mutation resulted in a loss of viral replicative fitness that was restored by compensatory changes such as L50F and T21I. Structural explanations are discussed for some of the mutations that are proximal to the drug-binding site, as well as cross-resistance or lack thereof to ensitrelvir, another clinically important 3CL protease inhibitor. Our findings indicate that SARS-CoV-2 resistance to nirmatrelvir does readily arise via multiple pathways in vitro , and the specific mutations observed herein form a strong foundation from which to study the mechanism of resistance in detail and to inform the design of next generation protease inhibitors.
