Next Article in Journal
The Disruption of a Nuclear Export Signal in the C-Terminus of the Herpes Simplex Virus 1 Determinant of Pathogenicity UL24 Protein Leads to a Syncytial Plaque Phenotype
Next Article in Special Issue
Exploring Conformational Landscapes and Cryptic Binding Pockets in Distinct Functional States of the SARS-CoV-2 Omicron BA.1 and BA.2 Trimers: Mutation-Induced Modulation of Protein Dynamics and Network-Guided Prediction of Variant-Specific Allosteric Binding Sites
Previous Article in Journal
Transcriptional Stochasticity as a Key Aspect of HIV-1 Latency
Previous Article in Special Issue
Characterizing Risk Factors for Hospitalization and Clinical Characteristics in a Cohort of COVID-19 Patients Enrolled in the GENCOV Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 15
Issue 9
10.3390/v15091970
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Substitutions in SARS-CoV-2 Mpro Selected by Protease Inhibitor Boceprevir Confer Resistance to Nirmatrelvir

by
Karen Anbro Gammeltoft
1,2,†,
Yuyong Zhou
1,2,†,
Line Abildgaard Ryberg
1,2,
Long V. Pham
1,2,
Alekxander Binderup
1,2,
Carlos Rene Duarte Hernandez
1,2,
Anna Offersgaard
1,2,
Ulrik Fahnøe
1,2,
Günther Herbert Johannes Peters
3,
Santseharay Ramirez
1,2,
Jens Bukh
1,2 and
Judith Margarete Gottwein
1,2,*
1
Copenhagen Hepatitis C Program (CO-HEP), Department of Infectious Diseases, Copenhagen University Hospital-Hvidovre, Kettegård Alle 30, 2650 Hvidovre, Denmark
2
Copenhagen Hepatitis C Program (CO-HEP), Department of Immunology and Microbiology, Faculty of Health and Medical Sciences, University of Copenhagen, Blegdamsvej 3B, 2200 Copenhagen, Denmark
3
Department of Chemistry, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2023, 15(9), 1970; https://doi.org/10.3390/v15091970
Submission received: 31 July 2023 / Revised: 11 September 2023 / Accepted: 14 September 2023 / Published: 21 September 2023
(This article belongs to the Collection SARS-CoV-2 and COVID-19)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

Since the identification of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 2019, the coronavirus disease 2019 (COVID-19) has become a major global public health threat. As of August 2023, there have been more than 770 million registered infections, posing a major socioeconomic burden. Reports from the World Health Organization (WHO) estimate that there was excess mortality in 2020 and 2021, with up to 16.6 million deaths attributed to COVID-19 [1,2]. Even though intense efforts led to the fast development of efficient prophylactic vaccines, treatment options for COVID-19 are still limited. Until recently, such treatments have primarily focused on alleviating symptoms or modulating the host’s immune system [3,4]. In contrast, there are only few approved drugs directly targeting SARS-CoV-2 proteins [5]. Such direct-acting antivirals include monoclonal antibodies that target the SARS-CoV-2 spike protein. At present, monoclonal antibodies have become inefficient and are no longer recommended due to currently circulating variants of concern [6,7]. At present, The United States Food and Drug Administration (FDA) has authorized three direct-acting antiviral therapeutics for treatment of COVID-19: (i) the polymerase inhibitor remdesivir (Veklury®, Gilead Sciences, Inc., Foster City, CA, USA), (ii) the polymerase inhibitor molnupiravir (Lagevrio®, Merck Sharp & Dohme LLC, Rahway, NJ, USA), and (iii) the main protease (Mpro) inhibitor nirmatrelvir (Paxlovid®, Pfizer, Inc., New York, NY, USA) [8]. Remdesivir, the first drug approved for COVID-19 treatment by both the FDA and the European Medicines Agency (EMA), is a nucleoside analog that targets the SARS-CoV-2 polymerase; it functions by template-dependent inhibition and by delay of RNA chain termination [9,10,11,12]. In clinical trials, treatment with remdesivir reduced the risk of hospitalization or death by 87% [13]. Nevertheless, the intravenous administration of remdesivir necessitates substantial resources [14]. The oral drug molnupiravir is a nucleoside analog that induces mutagenesis in the viral RNA genome. Clinical trials demonstrated a 31% reduction in the risk of hospitalization or death with molnupiravir treatment [15]. While it initially received emergency authorization from the FDA and EMA, it later failed to secure full authorization from the EMA due to its limited clinical efficacy [16,17,18,19]. The oral SARS-CoV-2 Mpro inhibitor nirmatrelvir has been approved for use by both the FDA and EMA and is currently the first-in-line COVID-19 treatment [20,21,22,23]. Nirmatrelvir has in clinical trials been shown to reduce the risk of COVID-19-associated hospitalization or death by 89% [24]. Further, nirmatrelvir is being considered for treatment of long-term effects of COVID-19 [25]. Additional Mpro inhibitors are in preclinical and clinical development, with some being in advanced clinical trials, such as ensitrelvir. Ensitrelvir, which has been approved in Japan, is now under Fast Track designation for COVID-19 treatment by the FDA [26,27,28]. Clinical trials demonstrated rapid viral clearance in patients with mild to moderate COVID-19 during treatment with ensitrelvir [29].
Because of the high mutation rate during RNA virus replication, it is expected that treatment of patients with direct-acting antivirals may give rise to treatment-resistant SARS-CoV-2 variants [30,31,32]. Further, there are case reports on patients failing treatment with nirmatrelvir [33,34,35]. Initial proof-of-concept studies have reported substitutions associated with resistance to currently authorized direct-acting antivirals including nirmatrelvir. Most of these studies were based on in vitro resistance selection and the study of naturally occurring polymorphisms in the drug targets [36,37,38,39,40,41,42,43,44]. In these studies, for nirmatrelvir, key resistance-associated substitutions (RAS) were found to be localized to Mpro position 166 [36,37,38,39,44]. Importantly, changes at Mpro position 166 were also observed in COVID-19 patients who were treated with nirmatrelvir, including one patient with treatment failure, highlighting the clinical relevance of the above-mentioned in vitro studies [45,46]. However, mechanisms leading to treatment failure remain to be studied, and a comprehensive overview of nirmatrelvir RAS is needed.
In previous drug-repurposing studies, we and others found that the hepatitis C virus (HCV) NS3 protease inhibitor boceprevir was active against SARS-CoV-2 by targeting Mpro [47,48,49]. However, the efficacy of boceprevir against SARS-CoV-2 is limited, and therefore, improved Mpro inhibitors have been developed. These efforts have led to the development of compounds, such as nirmatrelvir, that are structurally related to boceprevir but show improved activity against SARS-CoV-2 [50,51,52]. Here, we report boceprevir escape viruses with RAS in Mpro that had high fitness and conferred cross-resistance to nirmatrelvir in an infectious cell culture system. These findings will contribute to a more comprehensive overview of SARS-CoV-2 genetic correlates of resistance to nirmatrelvir, the first-in-line COVID-19 treatment.

2. Materials and Methods

2.1. Inhibitors, Cells, and Viruses

Boceprevir and nirmatrelvir were synthesized by Acme Bioscience (Palo Alto, CA, USA) and dissolved and stored as previously described [47].
All cell culture experiments were carried out in African green monkey kidney VeroE6 cells kindly provided by J. Dubuisson and human lung carcinoma A549-hACE2 cells (InvivoGen, Toulouse, France), which were maintained as previously described [36,47,53].
The virus stock of the so-called original SARS-CoV-2 was derived from the patient isolate SARS-CoV-2/human/Denmark/DK-AHH1/2020 (GenBank: MZ049597), representing an originally circulating variant with the D614G substitution in the spike protein, and was generated as described previously [54]. Stocks of the polyclonal boceprevir escape viruses were obtained as described in “Selection for SARS-CoV-2 Resistance to Boceprevir”. Stocks of the recombinant SARS-CoV-2 mutants were obtained as described in “Generation of Recombinant SARS-CoV-2 Mutants”. All infectious cell culture experiments were carried out under biosafety conditions in accordance with Danish regulations and with permission from the Danish authorities.

2.2. Selection for SARS-CoV-2 Resistance to Boceprevir

To select for SARS-CoV-2 resistance to boceprevir, the original SARS-CoV-2 was passaged in the presence of increasing concentrations of boceprevir as previously described for nirmatrelvir [36]. To initiate the passaging, VeroE6 cells were seeded in T25 flasks (Thermo Fisher Scientific, Roskilde, Denmark) at a density of 1 × 106 cells per flask, and both virus at 0.00002 multiplicity of infection (MOI) and specified concentrations of inhibitor were added the following day. Cell culture supernatants were collected every 48–72 h and stored at −80 °C; to maintain a subconfluent monolayer, cells were split, and the supernatant was replaced with fresh inhibitor containing medium. Upon splitting of cells, duplicate cultures were plated on chamber slides for immunostaining as described in “Immunostaining of Chamber Slides”. Five independent boceprevir escape experiments were carried out as specified. Selected supernatants were collected at the peak of infection and recovered viruses were subjected to next-generation sequencing (NGS). Polyclonal boceprevir escape virus 1 and 2 (BOC-EV1 and BOC-EV2) stocks were generated by inoculating 3 × 106 VeroE6 cells, plated the previous day in T80 flasks, with 15 µL supernatant derived from Day 74 from escape 1 or Day 57 from escape 2. Upon splitting of cells, duplicate cultures were plated on chamber slides for immunostaining as described in “Immunostaining of Chamber Slides”. Cell culture supernatants were harvested at peak of infection and stored at −80 °C. Polyclonal escape viruses were used for short-term concentration–response and longer-term treatments.

2.3. Generation of Recombinant SARS-CoV-2 Mutants

A bacterial artificial chromosome (BAC) based reverse genetics system [55] was used to generate SARS-CoV-2 recombinants harboring Mpro mutations. These mutants were engineered using infusion cloning and megaprimer-based cloning.
In vitro transcriptions and transfections using these constructs were carried out as previously described [36,55]. In brief, for in vitro transcriptions, the mMESSAGE mMACHINE T7 Transcription Kit (Thermo Fisher, Waltham, MA, USA) was used. The RNA transcripts were then transfected into VeroE6 cells using Lipofectamine 2000 (Thermo Fisher Scientific). A 1 µg amount of in vitro RNA transcripts was used for the L50F, A173V, and L50F + A173V recombinants, and 2 µg was used for the C160F, A191V, L50F + C160F + A173V, L50F + A173V + A191V and the original recombinants. Cell culture supernatants were collected on Day 2 and Day 3 post-transfection. A 250 µL volume of supernatant was then used to infect VeroE6 cells in T80 flasks to produce virus stocks of the recombinant SARS-CoV-2 mutants. Upon splitting of cells, duplicate cultures were plated on chamber slides for immunostaining as described in “Immunostaining of Chamber Slides”. Cell culture supernatants were harvested at the peak of infection and stored at −80 °C. Virus stocks were used for short-term concentration–response treatments, longer-term treatments, and viral fitness experiments.

2.4. Antiviral Short-Term Concentration-Response Treatments

Short-term concentration-response treatments were carried out as previously described [36,47,53]. In brief, on Day −1, VeroE6 or A549-hACE2 cells were seeded at 10,000 cells per well in 96-well plates (Thermo Fisher Scientific, Roskilde, Denmark). On Day 0, cells were incubated with original SARS-CoV-2 as a control, specified polyclonal SARS-CoV-2 escape viruses, or recombinant mutants. Following 60 min incubation, the specified serially diluted inhibitors were added. All dilutions were tested in 4 or 7 replicates. Each treatment experiment included a non-treated infected control in 8 or 14 replicates and a non-treated non-infected control in 8 or 12 replicates. The original SARS-CoV-2 was included in each treatment experiment as a reference. On Day 2 post-infection (DPI) and treatment, the cells were fixed and immunostained as described in “Immunostaining of 96-Well Plates”. The % residual infectivity was calculated as counts of individual treated infected wells related to mean counts of non-treated infected wells. The 50% effective concentration (EC50) values were determined using GraphPad Prism 8.0.0 and applying the formula Y = Top/(1 + 10(Log10EC50−X)∗HillSlope). Log10(EC50) and log10SEM(EC50) (SEM, standard error of the mean) derived from replicate experiments were used to determine p-values by Z-test, as previously described [56]. p-values below 0.0001 were considered significant.

2.5. Immunostaining of 96-Well Plates

Following short-term concentration–response treatments, 96-well plates were fixed with methanol and stained with first antibody SARS-CoV-2 spike chimeric monoclonal antibody (Sino Biological #40150-D004, Beijing, China) at 1:5000 dilution and second antibody F(ab’)2-Goat anti-human IgG Fc Cross-Adsorbed Secondary Antibody, HRP (Invitrogen #A24476, Carlsbad, CA, USA) or Goat F(ab’)2 Anti-Human IgG—Fc (HRP), preabsorbed (Abcam #ab98595, Cambridge, UK) at 1:2000 dilution. Next, infected cells were visualized by staining with DAB substrate, BrightDAB kit (Immunologic#BS04-110, Duiven, The Netherlands). The ImmunoSpot Series 5 UV Analyzer (CTL Europe GmbH, Bonn, Germany) was used to evaluate the number of single SARS-CoV-2 spike-positive cells per well. Representative images of stained 96-well plates are shown in a previous publication [57].

2.6. Longer-Term Antiviral Treatments

Longer-term treatments were carried out in VeroE6 cells. In brief, on Day −1, VeroE6 cells were seeded at 1 × 106 cells per flask in T25 flasks. On Day 0, specified polyclonal SARS-CoV-2 escape viruses, recombinant mutants, or original SARS-CoV-2 as a control were added at 0.00002 MOI together with the specified inhibitors. On Days 1, 3, and 5, cells were split, and supernatants were replaced with fresh inhibitor containing medium. Upon splitting of cells, cell culture supernatants were collected and stored at −80 °C for subsequent quantification of viral RNA, and duplicate cultures were plated on chamber slides for immunostaining as described in “Immunostaining of Chamber Slides”.

2.7. Immunostaining of Chamber Slides

During and following longer-term treatment assays, chamber slides were fixed with methanol and stained with first antibody SARS-CoV-2 spike chimeric monoclonal antibody (Sino Biological #40150-D004, Beijing, China) at 1:1000 dilution and second antibody Alexa Fluor 488 goat anti-human immunoglobulin G (IgG) (H + L) (Invitrogen #A-11013, Paisley, UK) at 1:500 dilution combined with Hoechst 33342 (Invitrogen, Paisley, UK) at 1:1000 dilution. Fluorescence microscopy (ZEISS Axio Vert.A1, Jena, Germany) was used to evaluate the percentages of SARS-CoV-2 spike-positive cells, using the following designations: 0% positive cells (no cells infected), single positive cells, and 10 to 90% positive cells, in steps of 10%.

2.8. Determination of SARS-CoV-2 RNA Titers

Assays were carried out as previously reported [36,47,58]. In brief, cell culture supernatants from longer-term treatments were extracted using Trizol LS Reagent (Life Technologies, Carlsbad, CA, USA) and chloroform (Sigma, Saint Louis, MI, USA). SARS-CoV-2 RNA was purified using the Zymo RNA Clean and Concentrator-5 kit (ZymoResearch, Irvine, CA, USA) for RT-qPCR according to the manufacturer’s protocol. Viral RNA was quantified by RT-qPCR on the LifeCycler 96 System (Roche) using the TaqMan Fast Virus 1-Step Master Mix (Thermo Fisher). For each sample, the RNA titer is given as a mean of two technical replicates. The lower limit of quantification (LLOQ) of the assay was determined as: (mean of RNA titers in non-infected control culture supernatants) + (3 standard deviations).

2.9. Determination of SARS-CoV-2 Infectivity Titers

Assays were carried out as previously reported [36,58]. In brief, 10,000 VeroE6 cells per well were seeded in 96-well plates on Day −1, and cells were infected with 10-fold serially diluted cell culture supernatants on Day 0 in 4 replicate cultures. On Day 3, cells were fixed and immunostained as described in “Immunostaining of 96-Well Plates”. SARS-CoV-2 infectivity titers were determined as 50% tissue culture-infectious dose per mL (TCID50/mL) using the Reed-Muench method. TCID50/mL titers are given as a means with SEM calculated from three replicate experiments with four technical replicates in each experiment. The lower limit of detection (LLOD) was 2 log10 TCID50/mL.

2.10. Evaluation of Genetic Stability of Recombinant SARS-CoV-2 Mutants

To evaluate their genetic stability, recombinant SARS-CoV-2 mutants were passaged four times in VeroE6 cells without inhibitor as previously described [36]. In brief, VeroE6 cells were seeded in T25 flasks at a density of 1 × 106 cells per flask on Day −1. Cells were infected with 250 µL supernatant derived at the peak of infection from the previous culture on Day 0. From Day 1, to maintain a subconfluent monolayer, cells were split, and the supernatant was replaced with fresh medium from Day 1. Upon splitting of cells, duplicate cultures were plated on chamber slides for immunostaining as described in “Immunostaining of Chamber Slides”. Viruses recovered from the fourth passage culture supernatants collected at the peak of infection were subjected to NGS analysis.

2.11. Next-Generation Sequencing of SARS-CoV-2 Genomes

The SARS-CoV-2 RNA was extracted from virus containing cell culture supernatant as described above; RT-PCR was used to generate five amplicons and the NEBNext Ultra II FS DNA Library Prep kit (New England BioLabs, Ipswich, MA, USA) was used for library preparations [54]. NGS analysis was done as previously described [36,47,53,54,55]. Linkage analysis of Mpro substitutions was done as previously described [32].

2.12. Analysis of Cell Viability by the MTS Assay in VeroE6 and A549-hACE2 Cells

To confirm that applied concentrations of inhibitors were not cytotoxic, cell viability was evaluated using the CellTiter 96 Aqueous One Solution Cell Proliferation Assay (Promega, Madison, WI, USA) as previously described [36,47,53,57]. Inhibitor concentrations were tested in 3 replicates including a non-treated control in 12 replicates.

2.13. Structural Analysis of the Effect of the L50F and A173V Substitutions in the Mpro-boceprevir Structure

An X-ray crystallography structure of the SARS-CoV-2 Mpro with boceprevir bound was obtained from the Protein Data Bank (PDB) with PDB entry 7k40 [59]. The PyMOL Molecular Graphics System 2.5.0. Schrödinger, LLC was used for manipulations, analysis, and visualization of the structure. More specifically, we carried out in silico mutagenesis for mutating A173 to V173 and chose the best-fitting rotamer from a library of backbone-dependent rotamers. We used the PyMOL distance tool for visualizing hydrogen bonds between boceprevir and amino acids in the Mpro structure, and the PyMOL plugin shows bumps for visualizing steric clashes in the structure.

3. Results

3.1. Viral Escape from Boceprevir Treatment

We previously characterized the efficacy of a panel of protease inhibitors originally developed for treatment of chronic HCV infection against SARS-CoV-2 in VeroE6 and A549-hACE2 cells [47]. Compared to the other tested protease inhibitors, boceprevir, structurally related to the current first-in-line SARS-CoV-2 Mpro inhibitor nirmatrelvir, had a favorable selectivity index (SI=half maximal cytotoxic concentration (CC50)/EC50). This facilitated SARS-CoV-2 escape experiments in VeroE6 cells to identify substitutions associated with Mpro inhibitor resistance (Supplementary Table S1). During such escape experiments, we observed that the original SARS-CoV-2 was suppressed by 3-fold EC50 boceprevir. However, boceprevir escape viruses could overcome up to 7-fold EC50 boceprevir. Higher boceprevir concentrations could not be applied due to cytotoxicity [47]. Sequence analysis revealed that polyclonal escape viruses from 4 of 5 independent escape experiments acquired a combination of the substitutions L50F and A173V in Mpro in virtually all viral genomes, while the fifth only acquired A173V in 78% of the genomes of the viral population (Supplementary Table S1). Thus, viral escape from boceprevir was linked to acquisition of Mpro substitutions warranting further investigation.

3.2. Sensitivity of Polyclonal Boceprevir Escape Viruses to Boceprevir and Nirmatrelvir

Compared to the original virus, selected polyclonal escape viruses with Mpro substitutions L50F and A173V, in the following termed BOC-EV1 and BOC-EV2, showed up to 4.7-fold increased EC50 (p < 0.0001) compared to the original virus in short-term concentration–response treatments with boceprevir (Figure 1, Supplementary Table S2). In addition, resistance was confirmed in longer-term treatments where BOC-EV1 and BOC-EV2 could overcome 5-fold EC50 of boceprevir, while the original virus was fully suppressed with no infected cells observed and viral RNA titers close to or below the LLOQ (Figure 2). Importantly, when testing these polyclonal escape viruses for nirmatrelvir cross-resistance, we found that they had up to 7.3-decreased nirmatrelvir sensitivity (p < 0.0001) compared to the original virus in VeroE6 and A549-hACE2 cells (Figure 1, Supplementary Table S2). The observed difference in nirmatrelvir potency between the two cell types can be attributed to the inherent expression of an efflux transporter by the VeroE6 cells, leading nirmatrelvir to be effectively transported out of the cells [60]. We additionally confirmed nirmatrelvir cross-resistance in longer-term treatments, where BOC-EV1 and BOC-EV2 could overcome up to 7.5- and 6.5-fold EC50 of nirmatrelvir, respectively (Figure 2), while the original virus was fully suppressed. Thus, polyclonal escape viruses were resistant to boceprevir and showed cross-resistance to nirmatrelvir, the first-in-line drug for treatment of COVID-19.

3.3. Sensitivity of Mutants with Engineered Mpro Substitutions to Boceprevir and Nirmatrelvir

The identified putative resistance substitutions L50F and A173V were characterized using a reverse genetics system that reflected the sequence of the original SARS-CoV-2 used in the escape experiments [55]. Thus, we engineered SARS-CoV-2 mutants harboring L50F and A173V, singly and in combination. In short-term concentration–response treatments, L50F and A173V singly conferred only small changes in boceprevir sensitivity (up to 2.0- and 1.4-fold decrease in sensitivity (p < 0.0001), respectively). The double mutant showed a 1.8-fold decreased boceprevir sensitivity in VeroE6 and a 2.6-fold decreased sensitivity in A549-hACE2 cells (p < 0.0001) (Figure 3, Supplementary Table S2). These substitutions had a somewhat greater impact on viral sensitivity to nirmatrelvir. L50F and A173V singly both conferred a decrease in nirmatrelvir sensitivity of up to 2.3-fold (p < 0.0001), while the double mutant showed a 3.5- and 2.1-fold decrease in sensitivity in VeroE6 and A549-hACE2 cells (p < 0.0001), respectively. Resistance of the L50F + A173V double mutant to nirmatrelvir was confirmed in longer-term treatments, as this mutant spread under treatment with 7.5-fold EC50 of nirmatrelvir while the original virus was fully suppressed (Figure 4). Thus, the L50F + A173V double mutant showed resistance to nirmatrelvir. Identification of L50F + A173V as nirmatrelvir RAS contributes to a comprehensive overview of RAS mediating resistance to this clinically relevant drug.

3.4. Fitness of Mutants with Engineered Mpro Substitutions

Next, we investigated the impact of the identified substitutions on viral fitness. We demonstrated high viral fitness for the engineered L50F + A173V double mutant and the L50F single mutant, as they showed high infectivity titers, comparable to those of the original virus, in both transfection and viral passage cultures (Figure 5). In contrast, the A173V single mutant showed markedly reduced fitness with viral infectivity titers close to or below the LLOD in transfection cultures. For viral passage cultures, recorded peak infectivity titers showed smaller differences. In addition, we carried out an analysis of genetic stability of these engineered mutants during a serial passage without drug pressure by sequencing analysis. The double mutant and the L50F mutant maintained the engineered substitutions in Mpro during four consecutive viral passages (Supplementary Table S3). In contrast, in the A173V mutant, the engineered substitution had partially reverted after four viral passages. None of the mutants acquired additional substitutions in Mpro during four viral passages (Supplementary Table S3). Moreover, in one independent experiment, the A173V mutant acquired L50F in 33% of viral genomes in the second viral passage. Thus, L50F compensated for the fitness cost of A173V. The L50F + A173V double mutant had high fitness, which could facilitate its emergence, spread, and persistence in populations.

3.5. Contribution of Additional Mpro Substitutions to Viral Resistance and Fitness

During selection of polyclonal escape viruses, sequencing analysis revealed that BOC-EV1 had acquired the additional Mpro substitutions C160F and A191V in 48% and 20% of viral genomes, respectively (Supplementary Table S1). Linkage analysis further revealed that C160F and A191V were not present on the same viral genomes, but individually in combination with L50F + A173V. To investigate the impact of these additional substitutions they were engineered singly and in combination with L50F + A173V. In short-term concentration–response treatments using boceprevir and nirmatrelvir, triple mutants L50F + A173V + C160F or L50F + A173V + A191V showed a decrease in susceptibility to both drugs that was comparable to that of the L50F + A173V double mutant (Figure 6). Similarly, in TCID50 assays, the triple mutants showed infectivity titers that were comparable to that of the L50F + A173V double mutant. When engineered singly, C160F or A191V did not show an impact on resistance or fitness of SARS-CoV-2. Thus, the Mpro substitutions C160F and A191V had no obvious impact on viral resistance or fitness using the applied assays.

3.6. Natural Occurrence of Identified Mpro Substitution

To assess the likelihood of emergence of the identified RAS we investigated if they are naturally occurring. Analysis of sequences extracted from the GISAID database prior to the widespread use of nirmatrelvir in the clinic showed that Mpro positions 50, 160, 173, and 191 were largely conserved but showed some degree of natural variation (Figure 7, Supplementary Table S4). Specifically, substitutions L50F, C160F, and A191V were found at comparable frequencies of 0.04%, 0.03%, and 0.09%, while substitution A173V was found at 0.002%. Thus, all identified Mpro substitutions including L50F and A173V were naturally occurring, even though at a low percentage. The natural occurrence of Mpro substitutions increases the likelihood of their spread in populations, considering a more widespread treatment with nirmatrelvir in the future.

3.7. Predicted Influence of Identified Mpro Substitutions on Mpro Structure

To provide a structural explanation for the observed resistance we carried out an analysis of the structural consequences of the identified RAS. Structural analysis revealed that in the Mpro structure, L50 and A173 are positioned on opposite sides of the S2 subsite of the active site defined by M49 and M165 [61] (Figure 8A). The backbone of A173 forms hydrogen bonds to the backbone of M165, and the sidechain is buried in a hydrophobic environment (Figure 8B). The A173V mutation introduces a bulkier sidechain that, as suggested by in silico mutagenesis, will lead to steric clashes with M165, L167, and F185 (Figure 8C). Thus, structural Mpro analysis suggested that A173V resulted in changes in the Mpro active site. Changes in the Mpro active site induced by A173V are likely perturbing its interaction with Mpro inhibitors, resulting in the observed resistance.

4. Discussion

In this study, we identified nirmatrelvir RAS L50F and A173V based on resistance selection during boceprevir escape experiments using an infectious SARS-CoV-2 cell culture system. While A173V was the main driver of resistance, L50F compensated for the fitness cost of A173V. Thus, the identified L50F + A173V escape viruses showed resistance and high fitness in our in vitro model. Importantly, this study refines the resistance profile of the first-in-line compound for treatment of COVID-19, nirmatrelvir. It is highly relevant to work towards a comprehensive overview of Mpro inhibitor RAS, which goes beyond the findings from initial proof-of-concept studies [36,37,38,39,44], as current COVID-19 treatments rely on the Mpro inhibitor nirmatrelvir and as additional protease inhibitors are being developed.
To our knowledge, by now, SARS-CoV-2 escape from boceprevir has not been reported. Even though boceprevir has limited clinical potential for treatment of COVID-19, it is relevant to study boceprevir escape because of the structural similarities between boceprevir and nirmatrelvir, and additional Mpro inhibitors in clinical development. Thus, as exemplified in this study, resistance profiles can potentially overlap among structurally related Mpro inhibitors, and therefore, boceprevir escape studies can help identify novel RAS of clinical relevance.
Due to its cytotoxicity, we could apply only relatively low concentrations of boceprevir. This likely explains why A173V was selected, which only conferred a relatively low degree of resistance compared to substitutions at position 166, reported by us and others to confer high resistance to nirmatrelvir [36,37,38,39,44]; the effect of RAS at position 166 on boceprevir resistance remains to be investigated in the future. Interestingly, the two polyclonal escape viruses showed somewhat higher resistance than the engineered L50F + A173V mutant. Such small discrepancies in phenotype for recombinant mutants versus polyclonal escape viruses have been observed for other viruses and inhibitors and might be due to substitutions outside the drug target contributing to resistance, either directly or indirectly via enhanced fitness [62,63].
Underlining the relevance of A173V as a nirmatrelvir RAS, a previous study identified A173V in SARS-CoV-2 selected during nirmatrelvir escape experiments in cell culture. Reverse genetics in an infectious culture system revealed that A173V conferred a small (1.8-fold) decrease in nirmatrelvir susceptibility in VeroE6 cells. While the combination of A173V with T21I or T21I + T304I resulted in a 3.1- and 15-fold decrease in nirmatrelvir susceptibility, respectively, the combination of A173V with L50F was not investigated [37]. Another study investigated naturally occurring Mpro polymorphisms in an infectious culture system. Here, the authors reported an 8.1-fold decrease in nirmatrelvir susceptibility for A173V and a 62.5-fold decrease for ΔP168 + A173V, while ΔP168 alone did not confer resistance [44]. In line with these findings, we described that culture-infectious SARS-CoV-2 escape viruses and engineered mutants with A173V have decreased nirmatrelvir susceptibility. In addition, we investigated the L50F + A173V double mutant, which has not been studied before by reverse genetics. Importantly, we confirmed that A173V together with L50F mediated resistance in longer-term treatments. We have previously shown that small EC50 differences in short-term treatment assays can translate into bigger differences in longer-term treatments for HCV protease inhibitor RAS, which had been associated with clinical resistance [64]. The maximum serum concentration (Cmax)/EC50 ratio for nirmatrelvir is 1.1 and 55 for the original SARS-CoV-2 in VeroE6 and A549-hACE2 cells, respectively, based on EC50 determined in this and a previous study [36]. The Cmax/EC50 ratio was reduced to 0.3 and 22 for the engineered L50F + A173V mutant and to 0.1 and 8.8 for BOC-EV1, the most resistant of the two polyclonal escape viruses, in VeroE6 and A549-hACE2 cells, respectively. While such a reduction could have an impact on clinical efficacy, the clinical relevance of these substitutions needs to be determined in future research. While certain substitutions in Mpro, including E166V, have been detected in patients following nirmatrelvir treatment, A173V has so far not been detected [21].
Together with our previous findings where we showed that L50F compensated for the fitness cost of E166V [36], our current data suggest that L50F has a broader fitness compensating effect for Mpro resistance substitutions. Importantly, we have shown that the pre-existence of naturally occurring L50F facilitated the selection of nirmatrelvir-resistant SARS-CoV-2 escape viruses with additional Mpro substitutions [36]. Thus, the background frequency of L50F and A173V, as recorded in the GISAID database as well as the high fitness of L50F + A173V, increase the risk of emergence and subsequent spread of Mpro inhibitor-resistant variants within populations. The emergence of Mpro inhibitor-resistant variants could potentially jeopardize the availability of effective treatment options for patients with severe COVID-19.
In one escape experiment, we identified Mpro substitutions C160F or A191V, selected in addition to L50F + A173V. Further reverse genetics analysis showed that C160F or A191V did not contribute to resistance in short-term treatment assays or viral fitness in TCID50 assays. These substitutions might be co-selected polymorphisms; however, A191V also occurred in the fifth escape experiment, even though on different genomes than A173V. It might be more likely that these substitutions slightly contribute to resistance or fitness to a level that cannot be detected in the applied assays.
L50F and A173V are located at structural elements that are important for the conformational flexibility of the active site, respectively, in the P2 helix and near the P3-P4 loop (Figure 8A). This conformational flexibility determines which substrates/inhibitors can bind to Mpro [48]. We hypothesize that the resistance observed for A173V against boceprevir and nirmatrelvir can be explained by the steric clashes (Figure 8C), which may induce changes in the conformational dynamics of the active site region, leading to a less favorable binding of boceprevir and nirmatrelvir. More specifically, we hypothesize that the steric clashes may change the interactions between the active site methionines (M49 and M165) and the inhibitors. This is supported for nirmatrelvir by molecular dynamics simulation (MDS) results of Mpro-nirmatrelvir showing that A173V increases the flexibility of M49 and the P2 helix [44]. The lower fitness of the A173V mutant could possibly be explained by A173V also leading to less favorable substrate binding. Considering the location of L50 in the P2 helix and its proximity to M49, the fitness compensatory effect of L50F might be explained by L50F + A173V compensating the conformational changes in the P2 helix induced by A173 and thereby improving substrate binding while maintaining reduced affinity for the inhibitors.
Future research should focus on further refinement of resistance profiles for Mpro inhibitors such as nirmatrelvir and ensitrelvir, also investigating cross-resistance. To this end, it will be relevant to carry out MDS analyses to assess the structural consequences of RAS on interactions of Mpro with relevant inhibitors and with its natural substrates.
A drawback of the study is that it is in vitro only, which does not fully reflect in vivo conditions, and thus, it will be important to further investigate the clinical relevance of identified RAS. Resistant Mpro variants have been shown to be transmissible in hamsters [65] and have already been found in patients. However, more systematic studies are needed to associate treatment failure with the identified RAS. Thus, for patients who experience treatment failure, it will be relevant to carry out sequence analysis to investigate if they are infected with SARS-CoV-2 variants that harbor putative RAS. Finally, it will be important to develop Mpro inhibitors with efficacy against SARS-CoV-2 variants with so far described nirmatrelvir RAS, by screening of drug libraries or by optimization of current Mpro inhibitors.
In conclusion, we report highly fit SARS-CoV-2 mutants with substitutions L50F and A173V in Mpro that confer decreased susceptibility to the oral SARS-CoV-2 Mpro inhibitor, nirmatrelvir. This study contributes to refining the resistance profile of nirmatrelvir, the first-in-line COVID-19 treatment, and the identified RAS could be included in resistance testing of emerging Mpro inhibitors. Characterization of Mpro inhibitor resistance profiles facilitates population monitoring to ensure the availability of safe and efficient treatments for patients suffering from COVID-19. Further, as Mpro is highly conserved between different coronaviruses, the definition of Mpro inhibitor resistance profiles contributes to pandemic preparedness.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/v15091970/s1, Supplementary Table S1: Non-synonymous mutations in boceprevir escape viruses localizing to sequences encoding Mpro; Supplementary Table S2: Potency of boceprevir and nirmatrelvir against original SARS-CoV-2 and SARS-CoV-2 mutants; Supplementary Table S3: Non-synonymous mutations in serially passaged SARS-CoV-2 mutants in the complete open reading frame (ORF); Supplementary Table S4: Naturally occurring substitutions in SARS-CoV-2 Mpro.

Author Contributions

K.A.G., Y.Z. and J.M.G. conceived this project. K.A.G., Y.Z. and J.M.G. designed the experiments. K.A.G., Y.Z., L.A.R., L.V.P., A.B., C.R.D.H. and A.O. carried out the experiments. All authors analyzed and interpreted the data (K.A.G., Y.Z., L.A.R., L.V.P., A.B., C.R.D.H., A.O., U.F., G.H.J.P., S.R., J.B. and J.M.G.). K.A.G., Y.Z., L.A.R. and J.M.G. prepared an initial manuscript draft. All authors contributed to the manuscript (K.A.G., Y.Z., L.A.R., L.V.P., A.B., C.R.D.H., A.O., U.F., G.H.J.P., S.R., J.B. and J.M.G.). J.M.G. supervised the study. All authors have read and agreed to the published version of the manuscript (K.A.G., Y.Z., L.A.R., L.V.P., A.B., C.R.D.H., A.O., U.F., G.H.J.P., S.R., J.B. and J.M.G.).

Funding

This research was funded by PhD stipends from the Candys Foundation (to K.A.G., A.O., J.B. and J.M.G.) and the China Scholarship Council (to Y.Z. and J.M.G.) and by grants from the Amager and Hvidovre Hospital Research Foundation (to C.R.D.H. and J.M.G.), the Danish Agency for Science and Higher Education (to S.R., J.B. and J.M.G.), the Independent Research Fund Denmark (to J.B.), the Innovation Fund Denmark (to J.M.G.), the Novo Nordisk Foundation including a Distinguished Investigator Grant (to J.B.), the Mauritzen la Fontaine Fonden (to J.B.), the Mauritzen La Fontaine Familiefond (to J.B.), the Region H Foundation (to J.B. and J.M.G.), the Toyota Foundation (to J.M.G.), and the Weimann Foundation (to U.F.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant data are included in the manuscript.

Acknowledgments

We would like to thank Bjarne Ørskov Lindhardt (Copenhagen University Hospital-Hvidovre) and Charlotte Menne Bonefeld (University of Copenhagen) for support and Lotte Mikkelsen and Anna-Louise Sørensen (Copenhagen University Hospital-Hvidovre) for laboratory assistance. We thank Jean Dubuisson and Sandrine Belouzard for providing VeroE6 cells.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization WHO Coronavirus Disease (COVID-19). Available online: https://covid19.who.int/ (accessed on 28 June 2023).
  2. Acosta, E. Global Estimates of Excess Deaths from COVID-19. Nature 2022, 613, 31–33. [Google Scholar] [CrossRef]
  3. Wang, Z.; Yang, L. The Therapeutic Potential of Natural Dietary Flavonoids against SARS-CoV-2 Infection. Nutrients 2023, 15, 3443. [Google Scholar] [CrossRef] [PubMed]
  4. Li, G.; Hilgenfeld, R.; Whitley, R.; De Clercq, E. Therapeutic Strategies for COVID-19: Progress and Lessons Learned. Nat. Rev. Drug Discov. 2023, 22, 449–475. [Google Scholar] [CrossRef]
  5. Yang, L.; Wang, Z. Bench-to-Bedside: Innovation of Small Molecule Anti-SARS-CoV-2 Drugs in China. Eur. J. Med. Chem. 2023, 257, 115503. [Google Scholar] [CrossRef] [PubMed]
  6. Cox, M.G.; Peacock, T.P.; Harvey, W.T.; Hughes, J.; Wright, D.W.; Willett, B.J.; Thomson, E.; Gupta, R.K.; Peacock, S.J.; Robertson, D.L.; et al. SARS-CoV-2 Variant Evasion of Monoclonal Antibodies Based on in Vitro Studies. Nat. Rev. Microbiol. 2022, 21, 112–124. [Google Scholar] [CrossRef] [PubMed]
  7. Focosi, D.; McConnell, S.; Casadevall, A.; Cappello, E.; Valdiserra, G.; Tuccori, M. Monoclonal Antibody Therapies against SARS-CoV-2. Lancet Infect. Dis. 2022, 22, e311–e326. [Google Scholar] [CrossRef]
  8. Side-by-Side Overview of Therapeutics Authorized or Approved for the Prevention of COVID-19 Infection or Treatment of Mild to Moderate COVID-19. Available online: https://covidpr.pregistry.com (accessed on 4 September 2023).
  9. Fact Sheet for Health Care Providers Emergency Use Authorization (Eua) of Veklury ® (Remdesivir). Available online: https://www.fda.gov/media/143189/download (accessed on 11 June 2023).
  10. Veklury|European Medicines Agency. Available online: https://www.ema.europa.eu/en/medicines/human/EPAR/veklury (accessed on 28 June 2023).
  11. FDA Approves First Treatment for COVID-19|FDA. Available online: https://www.fda.gov/news-events/press-announcements/fda-approves-first-treatment-covid-19 (accessed on 30 July 2023).
  12. Lamb, Y.N. Remdesivir: First Approval. Drugs 2020, 80, 1355–1363. [Google Scholar] [CrossRef]
  13. Gottlieb, R.L.; Vaca, C.E.; Paredes, R.; Mera, J.; Webb, B.J.; Perez, G.; Oguchi, G.; Ryan, P.; Nielsen, B.U.; Brown, M.; et al. Early Remdesivir to Prevent Progression to Severe COVID-19 in Outpatients. N. Engl. J. Med. 2022, 386, 305–315. [Google Scholar] [CrossRef]
  14. Wang, Z.; Yang, L.; Song, X.Q. Oral GS-441524 Derivatives: Next-Generation Inhibitors of SARS-CoV-2 RNA-dependent RNA Polymerase. Front. Immunol. 2022, 13, 1015355. [Google Scholar] [CrossRef]
  15. Jayk Bernal, A.; Gomes da Silva, M.M.; Musungaie, D.B.; Kovalchuk, E.; Gonzalez, A.; Delos Reyes, V.; Martín-Quirós, A.; Caraco, Y.; Williams-Diaz, A.; Brown, M.L.; et al. Molnupiravir for Oral Treatment of COVID-19 in Nonhospitalized Patients. N. Engl. J. Med. 2022, 386, 509–520. [Google Scholar] [CrossRef]
  16. Dyer, O. COVID-19: FDA Expert Panel Recommends Authorising Molnupiravir but Also Voices Concerns. BMJ 2021, 375, n2984. [Google Scholar] [CrossRef] [PubMed]
  17. Fact Sheet for Patients and Caregivers Emergency Use Authorization (EUA) of LAGEVRIOTM (Molnupiravir) Capsules For Coronavirus Disease 2019 (COVID-19). Available online: https://www.fda.gov/media/155055/download (accessed on 11 June 2023).
  18. Refusal of the Marketing Authorisation for Lagevrio (Molnupiravir). Available online: https://www.ema.europa.eu/en/documents/smop-initial/questions-answers-refusal-marketing-authorisation-lagevrio-molnupiravir_en.pdf (accessed on 28 June 2023).
  19. EMA Issues Advice on Use of Lagevrio (Molnupiravir) for the Treatment of COVID-19|European Medicines Agency. Available online: https://www.ema.europa.eu/en/news/ema-issues-advice-use-lagevrio-molnupiravir-treatment-covid-19 (accessed on 15 June 2023).
  20. COVID-19: EMA Recommends Conditional Marketing Authorisation for Paxlovid. Available online: https://www.ema.europa.eu/en/news/covid-19-ema-recommends-conditional-marketing-authorisation-paxlovid (accessed on 28 June 2023).
  21. Fact Sheet for Healthcare Providers: Emergency Use Authorization for Paxlovid. Available online: https://www.fda.gov/media/155050/download (accessed on 28 June 2023).
  22. FDA Approves First Oral Antiviral for Treatment of COVID-19 in Adults | FDA. Available online: https://www.fda.gov/news-events/press-announcements/fda-approves-first-oral-antiviral-treatment-covid-19-adults (accessed on 15 June 2023).
  23. Paxlovid|European Medicines Agency. Available online: https://www.ema.europa.eu/en/medicines/human/EPAR/paxlovid (accessed on 30 July 2023).
  24. Hammond, J.; Leister-Tebbe, H.; Gardner, A.; Abreu, P.; Bao, W.; Wisemandle, W.; Baniecki, M.; Hendrick, V.M.; Damle, B.; Simón-Campos, A.; et al. Oral Nirmatrelvir for High-Risk, Nonhospitalized Adults with COVID-19. N. Engl. J. Med. 2022, 386, 1397–1408. [Google Scholar] [CrossRef] [PubMed]
  25. Peluso, M.J.; Anglin, K.; Durstenfeld, M.S.; Martin, J.N.; Kelly, J.D.; Hsue, P.Y.; Henrich, T.J.; Deeks, S.G. Effect of Oral Nirmatrelvir on Long COVID Symptoms: 4 Cases and Rationale for Systematic Studies. Pathog. Immun. 2022, 7, 95. [Google Scholar] [CrossRef] [PubMed]
  26. Sasaki, M.; Tabata, K.; Kishimoto, M.; Itakura, Y.; Kobayashi, H.; Ariizumi, T.; Uemura, K.; Toba, S.; Kusakabe, S.; Maruyama, Y.; et al. S-217622, a SARS-CoV-2 Main Protease Inhibitor, Decreases Viral Load and Ameliorates COVID-19 Severity in Hamsters. Sci. Transl. Med. 2023, 15, eabq4064. [Google Scholar] [CrossRef]
  27. Xocova® (Ensitrelvir Fumaric Acid) Tablets 125mg Approved in Japan for the Treatment of SARS-CoV-2 Infection, under the Emergency Regulatory Approval System|News|Shionogi Co., Ltd. Available online: https://www.shionogi.com/global/en/news/2022/11/e20221122.html (accessed on 6 September 2023).
  28. Shionogi Receives, U.S. FDA Fast Track Designation for Ensitrelvir Fumaric Acid, an Investigational Oral Antiviral for COVID-19. Available online: https://www.shionogi.com/global/en/news/2023/04/20230404.html (accessed on 6 September 2023).
  29. Mukae, H.; Yotsuyanagi, H.; Ohmagari, N.; Doi, Y.; Sakaguchi, H.; Sonoyama, T.; Ichihashi, G.; Sanaki, T.; Baba, K.; Tsuge, Y.; et al. Efficacy and Safety of Ensitrelvir in Patients With Mild-to-Moderate Coronavirus Disease 2019: The Phase 2b Part of a Randomized, Placebo-Controlled, Phase 2/3 Study. Clin. Infect. Dis. 2023, 76, 1403–1411. [Google Scholar] [CrossRef]
  30. Hayden, F.G.; Belshe, R.B.; Clover, R.D.; Hay, A.J.; Oakes, M.G.; Soo, W. Emergence and Apparent Transmission of Rimantadine-Resistant Influenza A Virus in Families. N. Engl. J. Med. 1989, 321, 1696–1702. [Google Scholar] [CrossRef]
  31. Mason, S.; Devincenzo, J.P.; Toovey, S.; Wu, J.Z.; Whitley, R.J. Comparison of Antiviral Resistance across Acute and Chronic Viral Infections. Antiviral Res. 2018, 158, 103–112. [Google Scholar] [CrossRef]
  32. Jensen, S.B.; Fahnøe, U.; Pham, L.V.; Serre, S.; Tang, Q.; Ghanem, L.; Pedersen, M.S.; Ramirez, S.; Humes, D.; Pihl, A.F.; et al. Evolutionary Pathways to Persistence of Highly Fit and Resistant Hepatitis C Virus Protease Inhibitor Escape Variants. Hepatology 2019, 70, 771–787. [Google Scholar] [CrossRef]
  33. Anderson, A.S.; Caubel, P.; Rusnak, J.M. Nirmatrelvir–Ritonavir and Viral Load Rebound in COVID-19. N. Engl. J. Med. 2022, 387, 1047–1049. [Google Scholar] [CrossRef]
  34. Charness, M.E.; Gupta, K.; Stack, G.; Strymish, J.; Adams, E.; Lindy, D.C.; Mohri, H.; Ho, D.D. Rebound of SARS-CoV-2 Infection after Nirmatrelvir-Ritonavir Treatment. N. Engl. J. Med. 2022, 387, 1045–1047. [Google Scholar] [CrossRef]
  35. Wang, Y.; Chen, X.; Xiao, W.; Zhao, D.; Feng, L. Rapid COVID-19 Rebound in a Severe COVID-19 Patient during 20-Day Course of Paxlovid. J. Infect. 2022, 85, e134. [Google Scholar] [CrossRef]
  36. Zhou, Y.; Gammeltoft, K.A.; Ryberg, L.A.; Pham, L.V.; Tjørnelund, H.D.; Binderup, A.; Duarte Hernandez, C.R.; Fernandez-Antunez, C.; Offersgaard, A.; Fahnøe, U.; et al. Nirmatrelvir-Resistant SARS-CoV-2 Variants with High Fitness in an Infectious Cell Culture System. Sci. Adv. 2022, 8, eadd7197. [Google Scholar] [CrossRef]
  37. Iketani, S.; Mohri, H.; Culbertson, B.; Hong, S.J.; Duan, Y.; Luck, M.I.; Annavajhala, M.K.; Guo, Y.; Sheng, Z.; Uhlemann, A.C.; et al. Multiple Pathways for SARS-CoV-2 Resistance to Nirmatrelvir. Nature 2023, 613, 558–564. [Google Scholar] [CrossRef] [PubMed]
  38. Heilmann, E.; Costacurta, F.; Moghadasi, S.A.; Ye, C.; Pavan, M.; Bassani, D.; Volland, A.; Ascher, C.; Weiss, A.K.H.; Bante, D.; et al. SARS-CoV-2 3CLpro Mutations Selected in a VSV-Based System Confer Resistance to Nirmatrelvir, Ensitrelvir, and GC376. Sci. Transl. Med. 2023, 15, eabq7360. [Google Scholar] [CrossRef] [PubMed]
  39. Jochmans, D.; Liu, C.; Donckers, K.; Stoycheva, A.; Boland, S.; Stevens, S.K.; De Vita, C.; Vanmechelen, B.; Maes, P.; Trüeb, B.; et al. The Substitutions L50F, E166A, and L167F in SARS-CoV-2 3CLpro Are Selected by a Protease Inhibitor In Vitro and Confer Resistance to Nirmatrelvir. MBio 2023, 14, e0281522. [Google Scholar] [CrossRef] [PubMed]
  40. Kiso, M.; Yamayoshi, S.; Iida, S.; Furusawa, Y.; Hirata, Y.; Uraki, R.; Imai, M.; Suzuki, T.; Kawaoka, Y. In Vitro and in Vivo Characterization of SARS-CoV-2 Resistance to Ensitrelvir. Nat. Commun. 2023, 14, 4231. [Google Scholar] [CrossRef] [PubMed]
  41. Szemiel, A.M.; Merits, A.; Orton, R.J.; MacLean, O.A.; Pinto, R.M.; Wickenhagen, A.; Lieber, G.; Turnbull, M.L.; Wang, S.; Furnon, W.; et al. In Vitro Selection of Remdesivir Resistance Suggests Evolutionary Predictability of SARS-CoV-2. PLoS Pathog. 2021, 17, e1009929. [Google Scholar] [CrossRef]
  42. Checkmahomed, L.; Carbonneau, J.; Du Pont, V.; Riola, N.C.; Perry, J.K.; Li, J.; Paré, B.; Simpson, S.M.; Smith, M.A.; Porter, D.P.; et al. In Vitro Selection of Remdesivir-Resistant SARS-CoV-2 Demonstrates High Barrier to Resistance. Antimicrob. Agents Chemother. 2022, 66, e0019822. [Google Scholar] [CrossRef]
  43. Stevens, L.J.; Pruijssers, A.J.; Lee, H.W.; Gordon, C.J.; Tchesnokov, E.P.; Gribble, J.; George, A.S.; Hughes, T.M.; Lu, X.; Li, J.; et al. Mutations in the SARS-CoV-2 RNA-Dependent RNA Polymerase Confer Resistance to Remdesivir by Distinct Mechanisms. Sci. Transl. Med. 2022, 14, eabo0718. [Google Scholar] [CrossRef] [PubMed]
  44. Moghadasi, S.A.; Heilmann, E.; Khalil, A.M.; Nnabuife, C.; Kearns, F.L.; Ye, C.; Moraes, S.N.; Costacurta, F.; Esler, M.A.; Aihara, H.; et al. Transmissible SARS-CoV-2 Variants with Resistance to Clinical Protease Inhibitors. Sci. Adv. 2023, 9, eade8778. [Google Scholar] [CrossRef]
  45. Emergency Use Authorization (EMA) for Paxlovid. Available online: https://www.fda.gov/media/155194/download (accessed on 8 September 2023).
  46. Zuckerman, N.S.; Bucris, E.; Keidar-Friedman, D.; Amsalem, M.; Brosh-Nissimov, T. Nirmatrelvir Resistance—De Novo E166V/L50V Mutations in an Immunocompromised Patient Treated with Prolonged Nirmatrelvir/Ritonavir Monotherapy Leading to Clinical and Virological Treatment Failure—A Case Report. Clin. Infect. Dis. 2023, ciad494. [Google Scholar] [CrossRef]
  47. Gammeltoft, K.A.; Zhou, Y.; Duarte Hernandez, C.R.; Galli, A.; Offersgaard, A.; Costa, R.; Pham, L.V.; Fahnøe, U.; Feng, S.; Scheel, T.K.H.; et al. Hepatitis C Virus Protease Inhibitors Show Differential Efficacy and Interactions with Remdesivir for Treatment of SARS-CoV-2 in Vitro. Antimicrob. Agents Chemother. 2021, 65, e0268020. [Google Scholar] [CrossRef]
  48. Kneller, D.W.; Phillips, G.; O’Neill, H.M.; Jedrzejczak, R.; Stols, L.; Langan, P.; Joachimiak, A.; Coates, L.; Kovalevsky, A. Structural Plasticity of SARS-CoV-2 3CL Mpro Active Site Cavity Revealed by Room Temperature X-ray Crystallography. Nat. Commun. 2020, 11, 3202. [Google Scholar] [CrossRef] [PubMed]
  49. Ma, C.; Sacco, M.D.; Hurst, B.; Townsend, J.A.; Hu, Y.; Szeto, T.; Zhang, X.; Tarbet, B.; Marty, M.T.; Chen, Y.; et al. Boceprevir, GC-376, and Calpain Inhibitors II, XII Inhibit SARS-CoV-2 Viral Replication by Targeting the Viral Main Protease. Cell Res. 2020, 30, 678–692. [Google Scholar] [CrossRef]
  50. Kneller, D.W.; Li, H.; Phillips, G.; Weiss, K.L.; Zhang, Q.; Arnould, M.A.; Jonsson, C.B.; Surendranathan, S.; Parvathareddy, J.; Blakeley, M.P.; et al. Covalent Narlaprevir- and Boceprevir-Derived Hybrid Inhibitors of SARS-CoV-2 Main Protease. Nat. Commun. 2022, 13, 2268. [Google Scholar] [CrossRef] [PubMed]
  51. Göhl, M.; Zhang, L.; El Kilani, H.; Sun, X.; Zhang, K.; Brönstrup, M.; Hilgenfeld, R. From Repurposing to Redesign: Optimization of Boceprevir to Highly Potent Inhibitors of the SARS-CoV-2 Main Protease. Molecules 2022, 27, 4292. [Google Scholar] [CrossRef] [PubMed]
  52. Alugubelli, Y.R.; Geng, Z.Z.; Yang, K.S.; Shaabani, N.; Khatua, K.; Ma, X.R.; Vatansever, E.C.; Cho, C.C.; Ma, Y.; Xiao, J.; et al. A Systematic Exploration of Boceprevir-Based Main Protease Inhibitors as SARS-CoV-2 Antivirals. Eur. J. Med. Chem. 2022, 240, 114596. [Google Scholar] [CrossRef]
  53. Zhou, Y.; Gammeltoft, K.A.; Galli, A.; Offersgaard, A.; Fahnøe, U.; Ramirez, S.; Bukh, J.; Gottwein, J.M. Efficacy of Ion-Channel Inhibitors Amantadine, Memantine and Rimantadine for the Treatment of SARS-CoV-2 In Vitro. Viruses 2021, 13, 2082. [Google Scholar] [CrossRef]
  54. Ramirez, S.; Fernandez-Antunez, C.; Galli, A.; Underwood, A.; Pham, L.V.; Ryberg, L.A.; Feng, S.; Pedersen, M.S.; Mikkelsen, L.S.; Belouzard, S.; et al. Overcoming Culture Restriction for SARS-CoV-2 in Human Cells Facilitates the Screening of Compounds Inhibiting Viral Replication. Antimicrob. Agents Chemother. 2021, 65, e0009721. [Google Scholar] [CrossRef]
  55. Fahnøe, U.; Pham, L.V.; Fernandez-Antunez, C.; Costa, R.; Rivera-Rangel, L.R.; Galli, A.; Feng, S.; Mikkelsen, L.S.; Gottwein, J.M.; Scheel, T.K.H.; et al. Versatile SARS-CoV-2 Reverse-Genetics Systems for the Study of Antiviral Resistance and Replication. Viruses 2022, 14, 172. [Google Scholar] [CrossRef]
  56. Gottwein, J.M.; Scheel, T.K.H.; Jensen, T.B.; Ghanem, L.; Bukh, J. Differential Efficacy of Protease Inhibitors against HCV Genotypes 2a, 3a, 5a, and 6a NS3/4A Protease Recombinant Viruses. Gastroenterology 2011, 141, 1067–1079. [Google Scholar] [CrossRef] [PubMed]
  57. Zhou, Y.; Gilmore, K.; Ramirez, S.; Settles, E.; Gammeltoft, K.A.; Pham, L.V.; Fahnøe, U.; Feng, S.; Offersgaard, A.; Trimpert, J.; et al. In Vitro Efficacy of Artemisinin-Based Treatments against SARS-CoV-2. Sci. Rep. 2021, 11, 15471. [Google Scholar] [CrossRef] [PubMed]
  58. Offersgaard, A.; Hernandez, C.R.D.; Pihl, A.F.; Costa, R.; Venkatesan, N.P.; Lin, X.; Van Pham, L.; Feng, S.; Fahnøe, U.; Scheel, T.K.H.; et al. SARS-CoV-2 Production in a Scalable High Cell Density Bioreactor. Vaccines 2021, 9, 706. [Google Scholar] [CrossRef]
  59. Andi, B.; Kumaran, D.; Kreitler, D.F.; Soares, A.S.; Keereetaweep, J.; Jakoncic, J.; Lazo, E.O.; Shi, W.; Fuchs, M.R.; Sweet, R.M.; et al. Hepatitis C Virus NS3/4A Inhibitors and Other Drug-like Compounds as Covalent Binders of SARS-CoV-2 Main Protease. Sci. Rep. 2022, 12, 12197. [Google Scholar] [CrossRef] [PubMed]
  60. Owen, D.R.; Allerton, C.M.N.; Anderson, A.S.; Aschenbrenner, L.; Avery, M.; Berritt, S.; Boras, B.; Cardin, R.D.; Carlo, A.; Coffman, K.J.; et al. An Oral SARS-CoV-2 Mpro Inhibitor Clinical Candidate for the Treatment of COVID-19. Science 2021, 374, 1586–1593. [Google Scholar] [CrossRef]
  61. MacDonald, E.A.; Frey, G.; Namchuk, M.N.; Harrison, S.C.; Hinshaw, S.M.; Windsor, I.W. Recognition of Divergent Viral Substrates by the SARS-CoV-2 Main Protease. ACS Infect. Dis. 2021, 7, 2591–2595. [Google Scholar] [CrossRef]
  62. Serre, S.B.N.; Krarup, H.B.; Bukh, J.; Gottwein, J.M. Identification of Alpha Interferon-Induced Envelope Mutations of Hepatitis C Virus In Vitro Associated with Increased Viral Fitness and Interferon Resistance. J. Virol. 2013, 87, 12776–12793. [Google Scholar] [CrossRef]
  63. Sheldon, J.; Beach, N.M.; Moreno, E.; Gallego, I.; Pineiro, D.; Martinez-Salas, E.; Gregori, J.; Quer, J.; Esteban, J.I.; Rice, C.M.; et al. Increased Replicative Fitness Can Lead to Decreased Drug Sensitivity of Hepatitis C Virus. J. Virol. 2014, 88, 12098–12111. [Google Scholar] [CrossRef]
  64. Pham, L.V.; Jensen, S.B.; Fahnøe, U.; Pedersen, M.S.; Tang, Q.; Ghanem, L.; Ramirez, S.; Humes, D.; Serre, S.B.N.; Schønning, K.; et al. HCV Genotype 1-6 NS3 Residue 80 Substitutions Impact Protease Inhibitor Activity and Promote Viral Escape. J. Hepatol. 2019, 70, 388–397. [Google Scholar] [CrossRef] [PubMed]
  65. Abdelnabi, R.; Jochmans, D.; Donckers, K.; Trüeb, B.; Ebert, N.; Weynand, B.; Thiel, V.; Neyts, J. Nirmatrelvir-Resistant SARS-CoV-2 Is Efficiently Transmitted in Female Syrian Hamsters and Retains Partial Susceptibility to Treatment. Nat. Commun. 2023, 14, 2124. [Google Scholar] [CrossRef]
Viruses 15 01970 g001
Figure 1. Short-term concentration–response treatments of polyclonal SARS-CoV-2 boceprevir escape viruses with boceprevir and nirmatrelvir. VeroE6 or A549-hACE2 cells, that were infected with boceprevir escape virus 1 (BOC-EV1), boceprevir escape virus 2 (BOC-EV2), or original SARS-CoV-2 and that were treated with specified concentrations of boceprevir (BOC) or nirmatrelvir (NIR) were visualized by immunostaining of spike protein and counted automatically. Data points represent residual infectivity determined as count of infected cells in treated infected wells relative to the mean of counts in non-treated infected control wells. Data points are given as means of 4 or 7 replicates with SEM. Representative curves and 50% effective concentration (EC50) values from replicate experiments were determined in GraphPad Prism 8.0.0. In the figure, one representative curve is shown, while EC50 values are median EC50 values calculated based on 1 to 13 replicate experiments. Fold resistance values were calculated as EC50Escape Virus/EC50Original. Median EC50 values, fold resistance values, and p-values are summarized in Supplementary Table S2. For the original SARS-CoV-2 treated with nirmatrelvir, the curves have been presented in a previous publication, while median EC50 values have been calculated in a different manner [36].
Figure 1. Short-term concentration–response treatments of polyclonal SARS-CoV-2 boceprevir escape viruses with boceprevir and nirmatrelvir. VeroE6 or A549-hACE2 cells, that were infected with boceprevir escape virus 1 (BOC-EV1), boceprevir escape virus 2 (BOC-EV2), or original SARS-CoV-2 and that were treated with specified concentrations of boceprevir (BOC) or nirmatrelvir (NIR) were visualized by immunostaining of spike protein and counted automatically. Data points represent residual infectivity determined as count of infected cells in treated infected wells relative to the mean of counts in non-treated infected control wells. Data points are given as means of 4 or 7 replicates with SEM. Representative curves and 50% effective concentration (EC50) values from replicate experiments were determined in GraphPad Prism 8.0.0. In the figure, one representative curve is shown, while EC50 values are median EC50 values calculated based on 1 to 13 replicate experiments. Fold resistance values were calculated as EC50Escape Virus/EC50Original. Median EC50 values, fold resistance values, and p-values are summarized in Supplementary Table S2. For the original SARS-CoV-2 treated with nirmatrelvir, the curves have been presented in a previous publication, while median EC50 values have been calculated in a different manner [36].
Viruses 15 01970 g001
Viruses 15 01970 g002
Figure 2. Longer-term treatment of polyclonal escape viruses with boceprevir and nirmatrelvir. VeroE6 cells that were infected with the original SARS-CoV-2 and the specified polyclonal boceprevir escape viruses (BOC-EV1 and BOC-EV2) were treated longer term with the specified fold-EC50 of boceprevir (BOC) or nirmatrelvir (NIR). % SARS-CoV-2 infected culture cells were determined by immunostaining of spike protein relative to immunostaining of nuclei on the specified days post-infection (DPI). Viral RNA (vRNA) titers in culture supernatants on the specified DPI were determined by RT-qPCR. LLOQ, lower limit of quantification. nd, not determined as culture was terminated due to virus induced cell death.
Figure 2. Longer-term treatment of polyclonal escape viruses with boceprevir and nirmatrelvir. VeroE6 cells that were infected with the original SARS-CoV-2 and the specified polyclonal boceprevir escape viruses (BOC-EV1 and BOC-EV2) were treated longer term with the specified fold-EC50 of boceprevir (BOC) or nirmatrelvir (NIR). % SARS-CoV-2 infected culture cells were determined by immunostaining of spike protein relative to immunostaining of nuclei on the specified days post-infection (DPI). Viral RNA (vRNA) titers in culture supernatants on the specified DPI were determined by RT-qPCR. LLOQ, lower limit of quantification. nd, not determined as culture was terminated due to virus induced cell death.
Viruses 15 01970 g002
Viruses 15 01970 g003
Figure 3. Short-term concentration–response treatments of engineered SARS-CoV-2 mutants with boceprevir and nirmatrelvir. VeroE6 cells or A549-hACE2 cells that were infected with SARS-CoV-2 engineered with specified RAS or original SARS-CoV-2 and that were treated with specified concentrations of boceprevir (BOC) or nirmatrelvir (NIR) were visualized by immunostaining of spike protein and counted automatically. Data points represent residual infectivity determined as count of infected cells in treated infected wells relative to the mean of counts of non-treated infected control wells. Data points are given as means of 4 or 7 replicates with SEM. Representative curves and 50% effective concentration (EC50) values from replicate experiments were determined in GraphPad Prism 8.0.0. In the figure, one representative curve is shown, while EC50 values are median EC50 values calculated based on 1 to 13 replicate experiments. Fold resistance values were calculated as EC50Mutant/EC50Original. Median EC50 values, fold resistance values, and p-values are summarized in Supplementary Table S2. For the original SARS-CoV-2 and the L50F mutant treated with nirmatrelvir, the curves have been presented in a previous publication, while median EC50 values have been calculated in a different manner [36].
Figure 3. Short-term concentration–response treatments of engineered SARS-CoV-2 mutants with boceprevir and nirmatrelvir. VeroE6 cells or A549-hACE2 cells that were infected with SARS-CoV-2 engineered with specified RAS or original SARS-CoV-2 and that were treated with specified concentrations of boceprevir (BOC) or nirmatrelvir (NIR) were visualized by immunostaining of spike protein and counted automatically. Data points represent residual infectivity determined as count of infected cells in treated infected wells relative to the mean of counts of non-treated infected control wells. Data points are given as means of 4 or 7 replicates with SEM. Representative curves and 50% effective concentration (EC50) values from replicate experiments were determined in GraphPad Prism 8.0.0. In the figure, one representative curve is shown, while EC50 values are median EC50 values calculated based on 1 to 13 replicate experiments. Fold resistance values were calculated as EC50Mutant/EC50Original. Median EC50 values, fold resistance values, and p-values are summarized in Supplementary Table S2. For the original SARS-CoV-2 and the L50F mutant treated with nirmatrelvir, the curves have been presented in a previous publication, while median EC50 values have been calculated in a different manner [36].
Viruses 15 01970 g003
Viruses 15 01970 g004
Figure 4. Longer-term treatment of engineered SARS-CoV-2 mutants with nirmatrelvir. VeroE6 cells that were infected with the engineered L50F + A173V double mutant or the original SARS-CoV-2 were treated longer term with the specified fold-EC50 of nirmatrelvir (NIR). % SARS-CoV-2 infected culture cells were determined by immunostaining of spike relative to immunostaining of nuclei on the specified days post-infection (DPI). Viral RNA (vRNA) titers in culture supernatants on the specified DPI were determined by RT-qPCR. LLOQ, lower limit of quantification. For the original SARS-CoV-2, these data have been extracted from a previous publication [36].
Figure 4. Longer-term treatment of engineered SARS-CoV-2 mutants with nirmatrelvir. VeroE6 cells that were infected with the engineered L50F + A173V double mutant or the original SARS-CoV-2 were treated longer term with the specified fold-EC50 of nirmatrelvir (NIR). % SARS-CoV-2 infected culture cells were determined by immunostaining of spike relative to immunostaining of nuclei on the specified days post-infection (DPI). Viral RNA (vRNA) titers in culture supernatants on the specified DPI were determined by RT-qPCR. LLOQ, lower limit of quantification. For the original SARS-CoV-2, these data have been extracted from a previous publication [36].
Viruses 15 01970 g004
Viruses 15 01970 g005
Figure 5. Infectivity titers of engineered SARS-CoV-2 mutants. Viral infectivity titers of supernatants derived from VeroE6 cell cultures infected with SARS-CoV-2 mutants with specified engineered RAS or the original virus. Cultures were either transfected with RNA transcripts or infected with supernatant from first viral passage cultures. As an exception, for infection with the A173V mutant, supernatant from the transfection experiment was applied. Supernatants for determination of infectivity titers were obtained on Days 2 and 3 post-transfection, while virus stocks were supernatant derived at the peak of infection of the respective first or second passage cultures. Datapoints represent 50% tissue culture-infectious doses (TCID50) per mL and are given as means of 3 replicates with SEM. LLOD, lower limit of detection. For the original SARS-CoV-2 and L50F, these data have been published previously [36]. Infectivity titers shown in one graph are determined in the same experiment.
Figure 5. Infectivity titers of engineered SARS-CoV-2 mutants. Viral infectivity titers of supernatants derived from VeroE6 cell cultures infected with SARS-CoV-2 mutants with specified engineered RAS or the original virus. Cultures were either transfected with RNA transcripts or infected with supernatant from first viral passage cultures. As an exception, for infection with the A173V mutant, supernatant from the transfection experiment was applied. Supernatants for determination of infectivity titers were obtained on Days 2 and 3 post-transfection, while virus stocks were supernatant derived at the peak of infection of the respective first or second passage cultures. Datapoints represent 50% tissue culture-infectious doses (TCID50) per mL and are given as means of 3 replicates with SEM. LLOD, lower limit of detection. For the original SARS-CoV-2 and L50F, these data have been published previously [36]. Infectivity titers shown in one graph are determined in the same experiment.
Viruses 15 01970 g005
Viruses 15 01970 g006
Figure 6. Short-term concentration–response treatments and infectivity titers of engineered SARS-CoV-2 mutants with additional Mpro substitutions. (A) VeroE6 cells that were infected with SARS-CoV-2 engineered with specified RAS or original SARS-CoV-2 and that were treated with specified concentrations of boceprevir (BOC) or nirmatrelvir (NIR) were visualized by immunostaining of spike protein and counted automatically. Data points represent residual infectivity determined as count of infected cells in treated infected wells relative to the mean of counts of non-treated infected controls. Data points are given as means of 4 or 7 replicates with SEM. Representative curves and 50% effective concentration (EC50) values from 1 to 13 replicate experiments were determined in GraphPad Prism 8.0.0. In the figure, one representative curve is shown, while EC50 values are median EC50 values calculated based on replicate experiments. Fold resistance values were calculated as EC50Mutant/EC50Original. Median EC50 values, fold resistance values, and p-values are summarized in Supplementary Table S2. For the original SARS-CoV-2 treated with nirmatrelvir, the curves have been presented in a previous publication, while median EC50 values have been calculated in a different manner [36]. (B) Viral infectivity titers of supernatants derived at the peak of infection from VeroE6 cell cultures infected with SARS-CoV-2 mutants with specified engineered RAS. Cultures were infected with supernatant derived at the peak of infection from first viral passage cultures. Datapoints represent 50% tissue culture-infectious doses (TCID50) per mL and are given as means of 3 replicates with SEM. LLOD, lower limit of detection. For the original SARS-CoV-2, these data have been published previously [36]. Infectivity titers shown are determined in the same experiment.
Figure 6. Short-term concentration–response treatments and infectivity titers of engineered SARS-CoV-2 mutants with additional Mpro substitutions. (A) VeroE6 cells that were infected with SARS-CoV-2 engineered with specified RAS or original SARS-CoV-2 and that were treated with specified concentrations of boceprevir (BOC) or nirmatrelvir (NIR) were visualized by immunostaining of spike protein and counted automatically. Data points represent residual infectivity determined as count of infected cells in treated infected wells relative to the mean of counts of non-treated infected controls. Data points are given as means of 4 or 7 replicates with SEM. Representative curves and 50% effective concentration (EC50) values from 1 to 13 replicate experiments were determined in GraphPad Prism 8.0.0. In the figure, one representative curve is shown, while EC50 values are median EC50 values calculated based on replicate experiments. Fold resistance values were calculated as EC50Mutant/EC50Original. Median EC50 values, fold resistance values, and p-values are summarized in Supplementary Table S2. For the original SARS-CoV-2 treated with nirmatrelvir, the curves have been presented in a previous publication, while median EC50 values have been calculated in a different manner [36]. (B) Viral infectivity titers of supernatants derived at the peak of infection from VeroE6 cell cultures infected with SARS-CoV-2 mutants with specified engineered RAS. Cultures were infected with supernatant derived at the peak of infection from first viral passage cultures. Datapoints represent 50% tissue culture-infectious doses (TCID50) per mL and are given as means of 3 replicates with SEM. LLOD, lower limit of detection. For the original SARS-CoV-2, these data have been published previously [36]. Infectivity titers shown are determined in the same experiment.
Viruses 15 01970 g006
Viruses 15 01970 g007
Figure 7. Naturally occurring variations in Mpro at amino acid positions 50, 160, 173, and 191. Datapoints represent % of viruses identified in the GISAID database with specified amino acid substitutions at the specified Mpro positions. del, deletion.
Figure 7. Naturally occurring variations in Mpro at amino acid positions 50, 160, 173, and 191. Datapoints represent % of viruses identified in the GISAID database with specified amino acid substitutions at the specified Mpro positions. del, deletion.
Viruses 15 01970 g007
Viruses 15 01970 g008
Figure 8. Structural overview and analysis of L50F and A173V in the Mpro-boceprevir structure (PDB entry: 7k40 [59]). (A) Figure showing the position of L50 and A173 (limegreen) in the active site of the Mpro-boceprevir structure. The structural elements P2 helix (orange), P3-P4 loop (red), and P5 loop (pink) are highlighted, and boceprevir is shown in an aquamarine stick representation. The catalytic residues (H41 and C145) and the methionines (M49 and M165) are shown as sticks. L50 and A173 are positioned on opposite sides of boceprevir in proximity to M49 and M165, respectively. (B) Interactions of A173 in the Mpro structure. The backbone of A173 forms hydrogen bonds with M165 (yellow dashes), whereas the side chain forms hydrophobic interactions with M165, L167, and F185. (C) Modeling of the substitution A173V predicts steric clashes (red discs highlighted by black boxes) between the sidechain of V173 and the surrounding residues M165, L167, and F185.
Figure 8. Structural overview and analysis of L50F and A173V in the Mpro-boceprevir structure (PDB entry: 7k40 [59]). (A) Figure showing the position of L50 and A173 (limegreen) in the active site of the Mpro-boceprevir structure. The structural elements P2 helix (orange), P3-P4 loop (red), and P5 loop (pink) are highlighted, and boceprevir is shown in an aquamarine stick representation. The catalytic residues (H41 and C145) and the methionines (M49 and M165) are shown as sticks. L50 and A173 are positioned on opposite sides of boceprevir in proximity to M49 and M165, respectively. (B) Interactions of A173 in the Mpro structure. The backbone of A173 forms hydrogen bonds with M165 (yellow dashes), whereas the side chain forms hydrophobic interactions with M165, L167, and F185. (C) Modeling of the substitution A173V predicts steric clashes (red discs highlighted by black boxes) between the sidechain of V173 and the surrounding residues M165, L167, and F185.
Viruses 15 01970 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gammeltoft, K.A.; Zhou, Y.; Ryberg, L.A.; Pham, L.V.; Binderup, A.; Hernandez, C.R.D.; Offersgaard, A.; Fahnøe, U.; Peters, G.H.J.; Ramirez, S.; et al. Substitutions in SARS-CoV-2 Mpro Selected by Protease Inhibitor Boceprevir Confer Resistance to Nirmatrelvir. Viruses 2023, 15, 1970. https://doi.org/10.3390/v15091970

AMA Style

Gammeltoft KA, Zhou Y, Ryberg LA, Pham LV, Binderup A, Hernandez CRD, Offersgaard A, Fahnøe U, Peters GHJ, Ramirez S, et al. Substitutions in SARS-CoV-2 Mpro Selected by Protease Inhibitor Boceprevir Confer Resistance to Nirmatrelvir. Viruses. 2023; 15(9):1970. https://doi.org/10.3390/v15091970

Chicago/Turabian Style

Gammeltoft, Karen Anbro, Yuyong Zhou, Line Abildgaard Ryberg, Long V. Pham, Alekxander Binderup, Carlos Rene Duarte Hernandez, Anna Offersgaard, Ulrik Fahnøe, Günther Herbert Johannes Peters, Santseharay Ramirez, and et al. 2023. "Substitutions in SARS-CoV-2 Mpro Selected by Protease Inhibitor Boceprevir Confer Resistance to Nirmatrelvir" Viruses 15, no. 9: 1970. https://doi.org/10.3390/v15091970

APA Style

Gammeltoft, K. A., Zhou, Y., Ryberg, L. A., Pham, L. V., Binderup, A., Hernandez, C. R. D., Offersgaard, A., Fahnøe, U., Peters, G. H. J., Ramirez, S., Bukh, J., & Gottwein, J. M. (2023). Substitutions in SARS-CoV-2 Mpro Selected by Protease Inhibitor Boceprevir Confer Resistance to Nirmatrelvir. Viruses, 15(9), 1970. https://doi.org/10.3390/v15091970

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop