Nirmatrelvir, which targets the SARS-CoV-2 main protease (Mpro), is the first-in-line drug for prevention and treatment of severe COVID-19, and additional Mpro inhibitors are in development. However, the risk of resistance development threatens the future efficacy of such direct-acting antivirals. To gain knowledge on viral correlates of resistance to Mpro inhibitors, we selected resistant SARS-CoV-2 under treatment with the nirmatrelvir-related protease inhibitor boceprevir. SARS-CoV-2 selected during five escape experiments in VeroE6 cells showed cross-resistance to nirmatrelvir with up to 7.3-fold increased half-maximal effective concentration compared to original SARS-CoV-2, determined in concentration–response experiments. Sequence analysis revealed that escape viruses harbored Mpro substitutions L50F and A173V. For reverse genetic studies, these substitutions were introduced into a cell-culture-infectious SARS-CoV-2 clone. Infectivity titration and analysis of genetic stability of cell-culture-derived engineered SARS-CoV-2 mutants showed that L50F rescued the fitness cost conferred by A173V. In the concentration–response experiments, A173V was the main driver of resistance to boceprevir and nirmatrelvir. Structural analysis of Mpro suggested that A173V can cause resistance by making boceprevir and nirmatrelvir binding less favorable. This study contributes to a comprehensive overview of the resistance profile of the first-in-line COVID-19 treatment nirmatrelvir and can thus inform population monitoring and contribute to pandemic preparedness.

Various SARS-CoV-2 remdesivir resistance-associated substitutions (RAS) have been reported, but a comprehensive comparison of their resistance levels is lacking. We identified novel RAS and performed head-to-head comparisons with known RAS in Vero E6 cells. A remdesivir escape polyclonal virus exhibited a 3.6-fold increase in remdesivir EC₅₀ and mutations throughout the genome, including substitutions in nsp12 (E796D) and nsp14 (A255S). However, in reverse-genetics infectious assays, viruses harboring both these substitutions exhibited only a slight decrease in remdesivir susceptibility (1.3-fold increase in EC₅₀). The nsp12-E796D substitution did not impair viral fitness (Vero E6 cells or Syrian hamsters) and was reported in a remdesivir-treated COVID-19 patient. In replication assays, a subgenomic replicon containing nsp12-E796D+nsp14-A255S led to a 16.1-fold increase in replication under remdesivir treatment. A comparison with known RAS showed that S759A, located in the active site of nsp12, conferred the highest remdesivir resistance (106.1-fold increase in replication). Nsp12-RAS V166A/L, V792I, E796D or C799F, all adjacent to the active site, caused intermediate resistance (2.0- to 11.5-fold), whereas N198S, D484Y, or E802D, located farther from the active site, showed no resistance (≤2.0-fold). In conclusion, our classification system, correlating replication under remdesivir treatment with RAS location in nsp12, shows that most nsp12-RAS cause moderate resistance.