Transmissible SARS-CoV-2 variants with resistance to clinical protease inhibitors

Efficacy of late-onset antiviral treatment in immunocompromised hosts with persistent SARS-CoV-2 infection

Immunocompromised people are at high risk of prolonged severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and progression to severe coronavirus disease 2019 (COVID-19). However, the efficacy of late-onset direct-acting antiviral (DAA) therapy with therapeutics in clinical use and experimental drugs to mitigate persistent viral replication is unclear. In this study, we employed an immunocompromised mouse model, which supports prolonged replication of SARS-CoV-2 to explore late-onset treatment options. Tandem immuno-depletion of CD4+ and CD8+ T cells in C57BL/6 mice followed by infection with SARS-CoV-2 variant of concern (VOC) beta B.1.351 resulted in prolonged infection with virus replication for 5 weeks after inoculation. Early-onset treatment with nirmatrelvir/ritonavir (paxlovid) or molnupiravir was only moderately efficacious, whereas the experimental therapeutic 4'-fluorouridine (4'-FlU, EIDD-2749) significantly reduced virus load in the upper and lower respiratory compartments 4 days postinfection (dpi). All antivirals significantly lowered virus burden in a 7-day treatment regimen initiated 14 dpi, but paxlovid-treated animals experienced rebound virus replication in the upper respiratory tract 7 days after treatment end. Viral RNA was detectable 28 dpi in paxlovid-treated animals, albeit not in the molnupiravir or 4'-FlU groups, when treatment was initiated 14 dpi and continued for 14 days. Low-level virus replication continued 35 dpi in animals receiving vehicle but had ceased in all treatment groups. These data indicate that late-onset DAA therapy significantly shortens the duration of persistent virus replication in an immunocompromised host, which may have implications for clinical use of antiviral therapeutics to alleviate the risk of progression to severe disease in highly vulnerable patients.

IMPORTANCE

Four years after the onset of the global coronavirus disease 2019 (COVID-19) pandemic, the immunocompromised are at greatest risk of developing life-threatening severe disease. However, specific treatment plans for this most vulnerable patient group have not yet been developed. Employing a CD4+ and CD8+ T cell-depleted immunocompromised mouse model of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, we explored therapeutic options of persistent infections with standard-of-care paxlovid, molnupiravir, and the experimental therapeutic 4'-fluorouridine (4'-FlU). Late-onset treatment initiated 14 days after infection was efficacious, but only 4'-FlU was rapidly sterilizing. No treatment-experienced viral variants with reduced susceptibility to the drugs emerged, albeit virus replication rebounded in animals of the paxlovid group after treatment end. This study supports the use of direct-acting antivirals (DAAs) for late-onset management of persistent SARS-CoV-2 infection in immunocompromised hosts. However, treatment courses likely require to be extended for maximal therapeutic benefit, calling for appropriately powered clinical trials to meet the specific needs of this patient group.

Efficacy of late-onset antiviral treatment in immune-compromised hosts with persistent SARS-CoV-2 infection

The immunocompromised are at high risk of prolonged SARS-CoV-2 infection and progression to severe COVID-19. However, efficacy of late-onset direct-acting antiviral (DAA) therapy with therapeutics in clinical use and experimental drugs to mitigate persistent viral replication is unclear. In this study, we employed an immunocompromised mouse model, which supports prolonged replication of SARS-CoV-2 to explore late-onset treatment options. Tandem immuno-depletion of CD4 + and CD8 + T cells in C57BL/6 mice followed by infection with SARS-CoV-2 variant of concern (VOC) beta B.1.351 resulted in prolonged infection with virus replication for five weeks after inoculation. Early-onset treatment with nirmatrelvir/ritonavir (paxlovid) or molnupiravir was only moderately efficacious, whereas the experimental therapeutic 4'-fluorourdine (4'-FlU, EIDD-2749) significantly reduced virus load in upper and lower respiratory compartments four days post infection (dpi). All antivirals significantly lowered virus burden in a 7-day treatment regimen initiated 14 dpi, but paxlovid-treated animals experienced rebound virus replication in the upper respiratory tract seven days after treatment end. Viral RNA was detectable 28 dpi in paxlovid-treated animals, albeit not in the molnupiravir or 4'-FlU groups, when treatment was initiated 14 dpi and continued for 14 days. Low-level virus replication continued 35 dpi in animals receiving vehicle but had ceased in all treatment groups. These data indicate that late-onset DAA therapy significantly shortens the duration of persistent virus replication in an immunocompromised host, which may have implications for clinical use of antiviral therapeutics to alleviate the risk of progression to severe disease in highly vulnerable patients.

Importance

Four years after the onset of the global COVID-19 pandemic, the immunocompromised are at greatest risk of developing life-threatening severe disease. However, specific treatment plans for this most vulnerable patient group have not yet been developed. Employing a CD4 + and CD8 + T cell-depleted immunocompromised mouse model of SARS-CoV-2 infection, we explored therapeutic options of persistent infections with standard-of-care paxlovid, molnupiravir, and the experimental therapeutic 4'-FlU. Late-onset treatment initiated 14 days after infection was efficacious, but only 4'-FlU was rapidly sterilizing. No treatment-experienced viral variants with reduced susceptibility to the drugs emerged, albeit virus replication rebounded in animals of the paxlovid group after treatment end. This study supports the use of direct-acting antivirals for late-onset management of persistent SARS-CoV-2 infection in immunocompromised hosts. However, treatment courses likely require to be extended for maximal therapeutic benefit, calling for appropriately powered clinical trials to meet the specific needs of this patient group.

A low-background, fluorescent assay to evaluate inhibitors of diverse viral proteases

Multiple coronaviruses (CoVs) can cause respiratory diseases in humans. While prophylactic vaccines designed to prevent infection are available for severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), incomplete vaccine efficacy, vaccine hesitancy, and the threat of other pathogenic CoVs for which vaccines do not exist have highlighted the need for effective antiviral therapies. While antiviral compounds targeting the viral polymerase and protease are already in clinical use, their sensitivity to potential resistance mutations as well as their breadth against the full range of human and preemergent CoVs remain incompletely defined. To begin to fill that gap in knowledge, we report here the development of an improved, noninfectious, cell-based fluorescent assay with high sensitivity and low background that reports on the activity of viral proteases, which are key drug targets. We demonstrate that the assay is compatible with not only the SARS-CoV-2 Mpro protein but also orthologues from a range of human and nonhuman CoVs as well as clinically reported SARS-CoV-2 drug-resistant Mpro variants. We then use this assay to define the breadth of activity of two clinically used protease inhibitors, nirmatrelvir and ensitrelvir. Continued use of this assay will help define the strengths and limitations of current therapies and may also facilitate the development of next-generation protease inhibitors that are broadly active against both currently circulating and preemergent CoVs. IMPORTANCE Coronaviruses (CoVs) are important human pathogens with the ability to cause global pandemics. Working in concert with vaccines, antivirals specifically limit viral disease in people who are actively infected. Antiviral compounds that target CoV proteases are already in clinical use; their efficacy against variant proteases and preemergent zoonotic CoVs, however, remains incompletely defined. Here, we report an improved, noninfectious, and highly sensitive fluorescent method of defining the sensitivity of CoV proteases to small molecule inhibitors. We use this approach to assay the activity of current antiviral therapies against clinically reported SARS-CoV-2 protease mutants and a panel of highly diverse CoV proteases. Additionally, we show this system is adaptable to other structurally nonrelated viral proteases. In the future, this assay can be used to not only better define the strengths and limitations of current therapies but also help develop new, broadly acting inhibitors that more broadly target viral families.

In vitro and in vivo characterization of SARS-CoV-2 strains resistant to nirmatrelvir

Nirmatrelvir, an oral antiviral agent that targets a SARS-CoV-2 main protease (3CLpro), is clinically useful against infection with SARS-CoV-2 including its omicron variants. Since most omicron subvariants have reduced sensitivity to many monoclonal antibody therapies, potential SARS-CoV-2 resistance to nirmatrelvir is a major public health concern. Several amino acid substitutions have been identified as being responsible for reduced susceptibility to nirmatrelvir. Among them, we selected L50F/E166V and L50F/E166A/L167F in the 3CLpro because these combinations of substitutions are unlikely to affect virus fitness. We prepared and characterized delta variants possessing Nsp5-L50F/E166V and Nsp5-L50F/E166A/L167F. Both mutant viruses showed decreased susceptibility to nirmatrelvir and their growth in VeroE6/TMPRSS2 cells was delayed. Both mutant viruses showed attenuated phenotypes in a male hamster infection model, maintained airborne transmissibility, and were outcompeted by wild-type virus in co-infection experiments in the absence of nirmatrelvir, but less so in the presence of the drug. These results suggest that viruses possessing Nsp5-L50F/E166V and Nsp5-L50F/E166A/L167F do not become dominant in nature. However, it is important to closely monitor the emergence of nirmatrelvir-resistant SARS-CoV-2 variants because resistant viruses with additional compensatory mutations could emerge, outcompete the wild-type virus, and become dominant.

Fitness effects of mutations to SARS-CoV-2 proteins

Knowledge of the fitness effects of mutations to SARS-CoV-2 can inform assessment of new variants, design of therapeutics resistant to escape, and understanding of the functions of viral proteins. However, experimentally measuring effects of mutations is challenging: we lack tractable lab assays for many SARS-CoV-2 proteins, and comprehensive deep mutational scanning has been applied to only two SARS-CoV-2 proteins. Here we develop an approach that leverages millions of publicly available SARS-CoV-2 sequences to estimate effects of mutations. We first calculate how many independent occurrences of each mutation are expected to be observed along the SARS-CoV-2 phylogeny in the absence of selection. We then compare these expected observations to the actual observations to estimate the effect of each mutation. These estimates correlate well with deep mutational scanning measurements. For most genes, synonymous mutations are nearly neutral, stop-codon mutations are deleterious, and amino-acid mutations have a range of effects. However, some viral accessory proteins are under little to no selection. We provide interactive visualizations of effects of mutations to all SARS-CoV-2 proteins (https://jbloomlab.github.io/SARS2-mut-fitness/). The framework we describe is applicable to any virus for which the number of available sequences is sufficiently large that many independent occurrences of each neutral mutation are observed.

A single inactivating amino acid change in the SARS-CoV-2 NSP3 Mac1 domain attenuates viral replication and pathogenesis in vivo

Despite unprecedented efforts, our therapeutic arsenal against SARS-CoV-2 remains limited. The conserved macrodomain 1 (Mac1) in NSP3 is an enzyme exhibiting ADP-ribosylhydrolase activity and a possible drug target. To determine the therapeutic potential of Mac1 inhibition, we generated recombinant viruses and replicons encoding a catalytically inactive NSP3 Mac1 domain by mutating a critical asparagine in the active site. While substitution to alanine (N40A) reduced catalytic activity by ~10-fold, mutations to aspartic acid (N40D) reduced activity by ~100-fold relative to wildtype. Importantly, the N40A mutation rendered Mac1 unstable in vitro and lowered expression levels in bacterial and mammalian cells. When incorporated into SARS-CoV-2 molecular clones, the N40D mutant only modestly affected viral fitness in immortalized cell lines, but reduced viral replication in human airway organoids by 10-fold. In mice, N40D replicated at >1000-fold lower levels compared to the wildtype virus while inducing a robust interferon response; all animals infected with the mutant virus survived infection and showed no signs of lung pathology. Our data validate the SARS-CoV-2 NSP3 Mac1 domain as a critical viral pathogenesis factor and a promising target to develop antivirals.

Molecular Epidemiology of SARS-CoV-2 Omicron Sub-Lineages Isolated from Turkish Patients Infected with COVID-19

Early detection and characterization of new variants and their impacts enable improved genomic surveillance. This study aims to evaluate the subvariant distribution of Omicron strains isolated from Turkish cases to determine the rate of antiviral resistance of RdRp and 3CLpro inhibitors. The Stanford University Coronavirus Antiviral & Resistance Database online tool was used for variant analyses of the strains uploaded to GISAID as Omicron (n = 20.959) between January 2021 and February,2023. Out of 288 different Omicron subvariants, B.1, BA.1, BA.2, BA.4, BE.1, BF.1, BM.1, BN.1, BQ.1, CK.1, CL.1, and XBB.1 were the main determined subvariants, and BA.1 (34.7%), BA.2 (30.8%), and BA.5 (23.6%) were reported most frequently. RdRp and 3CLPro-related resistance mutations were determined in n = 150, 0.72% sequences, while the rates of resistance against RdRp and 3CLpro inhibitors were reported at 0.1% and 0.6%, respectively. Mutations that were previously associated with a reduced susceptibility to remdesivir, nirmatrelvir/r, and ensitrelvir were most frequently detected in BA.2 (51.3%). The mutations detected at the highest rate were A449A/D/G/V (10.5%), T21I (10%), and L50L/F/I/V (6%). Our findings suggest that continuous monitoring of variants, due to the diversity of Omicron lineages, is necessary for global risk assessment. Although drug-resistant mutations do not pose a threat, the tracking of drug mutations will be necessary due to variant heterogenicity.

Rapid assembly of SARS-CoV-2 genomes reveals attenuation of the Omicron BA.1 variant through NSP6

Although the SARS-CoV-2 Omicron variant (BA.1) spread rapidly across the world and effectively evaded immune responses, its viral fitness in cell and animal models was reduced. The precise nature of this attenuation remains unknown as generating replication-competent viral genomes is challenging because of the length of the viral genome (~30 kb). Here, we present a plasmid-based viral genome assembly and rescue strategy (pGLUE) that constructs complete infectious viruses or noninfectious subgenomic replicons in a single ligation reaction with >80% efficiency. Fully sequenced replicons and infectious viral stocks can be generated in 1 and 3 weeks, respectively. By testing a series of naturally occurring viruses as well as Delta-Omicron chimeric replicons, we show that Omicron nonstructural protein 6 harbors critical attenuating mutations, which dampen viral RNA replication and reduce lipid droplet consumption. Thus, pGLUE overcomes remaining barriers to broadly study SARS-CoV-2 replication and reveals deficits in nonstructural protein function underlying Omicron attenuation.

The impact of SARS-CoV-2 3CL protease mutations on nirmatrelvir inhibitory efficiency. Computational insights into potential resistance mechanisms

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.

Alkyne Derivatives of SARS-CoV-2 Main Protease Inhibitors Including Nirmatrelvir Inhibit by Reacting Covalently with the Nucleophilic Cysteine

Nirmatrelvir (PF-07321332) is a nitrile-bearing small-molecule inhibitor that, in combination with ritonavir, is used to treat infections by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Nirmatrelvir interrupts the viral life cycle by inhibiting the SARS-CoV-2 main protease (Mpro), which is essential for processing viral polyproteins into functional nonstructural proteins. We report studies which reveal that derivatives of nirmatrelvir and other Mpro inhibitors with a nonactivated terminal alkyne group positioned similarly to the electrophilic nitrile of nirmatrelvir can efficiently inhibit isolated Mpro and SARS-CoV-2 replication in cells. Mass spectrometric and crystallographic evidence shows that the alkyne derivatives inhibit Mpro by apparent irreversible covalent reactions with the active site cysteine (Cys145), while the analogous nitriles react reversibly. The results highlight the potential for irreversible covalent inhibition of Mpro and other nucleophilic cysteine proteases by alkynes, which, in contrast to nitriles, can be functionalized at their terminal position to optimize inhibition and selectivity, as well as pharmacodynamic and pharmacokinetic properties.

Rapid assembly of SARS-CoV-2 genomes reveals attenuation of the Omicron BA.1 variant through NSP6

Although the SARS-CoV-2 Omicron variant (BA.1) spread rapidly across the world and effectively evaded immune responses, its viral fitness in cell and animal models was reduced. The precise nature of this attenuation remains unknown as generating replication-competent viral genomes is challenging because of the length of the viral genome (30kb). Here, we designed a plasmid-based viral genome assembly and resc ue strategy (pGLUE) that constructs complete infectious viruses or noninfectious subgenomic replicons in a single ligation reaction with >80% efficiency. Fully sequenced replicons and infectious viral stocks can be generated in 1 and 3 weeks, respectively. By testing a series of naturally occurring viruses as well as Delta-Omicron chimeric replicons, we show that Omicron nonstructural protein 6 harbors critical attenuating mutations, which dampen viral RNA replication and reduce lipid droplet consumption. Thus, pGLUE overcomes remaining barriers to broadly study SARS-CoV-2 replication and reveals deficits in nonstructural protein function underlying Omicron attenuation.

A strategy for evaluating potential antiviral resistance to small molecule drugs and application to SARS-CoV-2

Alterations in viral fitness cannot be inferred from only mutagenesis studies of an isolated viral protein. To-date, no systematic analysis has been performed to identify mutations that improve virus fitness and reduce drug efficacy. We present a generic strategy to evaluate which viral mutations might diminish drug efficacy and applied it to assess how SARS-CoV-2 evolution may affect the efficacy of current approved/candidate small-molecule antivirals for Mpro, PLpro, and RdRp. For each drug target, we determined the drug-interacting virus residues from available structures and the selection pressure of the virus residues from the SARS-CoV-2 genomes. This enabled the identification of promising drug target regions and small-molecule antivirals that the virus can develop resistance. Our strategy of utilizing sequence and structural information from genomic sequence and protein structure databanks can rapidly assess the fitness of any emerging virus variants and can aid antiviral drug design for future pathogens.

Genomic Tracking of SARS-CoV-2 Variants in Myanmar

In December 2019, the COVID-19 disease started in Wuhan, China. The WHO declared a pandemic on 12 March 2020, and the disease started in Myanmar on 23 March 2020. In December 2020, different variants were brought worldwide, threatening global health. To counter those threats, Myanmar started the COVID-19 variant surveillance program in late 2020. Whole genome sequencing was done six times between January 2021 and March 2022. Among them, 83 samples with a PCR threshold cycle of less than 25 were chosen. Then, we used MiSeq FGx for sequencing and Illumina DRAGEN COVIDSeq pipeline, command line interface, GISAID, and MEGA version 7 for data analysis. In January 2021, no variant was detected. The second run, during the rise of cases in June 2021, showed Alpha, Delta, and Kappa variants. The third and the fourth runs in August and December showed only a Delta variant. Omicron and Delta variants were detected during the fifth run in January 2022. The sixth run in March 2022 showed only Omicron BA.2. Amino acid mutation at the receptor binding domain of Spike glycoprotein started since the second run coupling with high transmission, recurrence, and vaccine escape. We also found the mutation at the primer targets used in current RT-PCR platforms, but there was no mutation at the existing antiviral drug targets. The occurrence of multiple variants and mutations claimed vigilance at ports of entry and preparedness for effective control measures. Genomic surveillance with the observation of evolutionary data is required to predict imminent threats of the current disease and diagnose emerging infectious diseases.

Small molecules in the treatment of COVID-19

The outbreak of COVID-19 has become a global crisis, and brought severe disruptions to societies and economies. Until now, effective therapeutics against COVID-19 are in high demand. Along with our improved understanding of the structure, function, and pathogenic process of SARS-CoV-2, many small molecules with potential anti-COVID-19 effects have been developed. So far, several antiviral strategies were explored. Besides directly inhibition of viral proteins such as RdRp and Mpro, interference of host enzymes including ACE2 and proteases, and blocking relevant immunoregulatory pathways represented by JAK/STAT, BTK, NF-κB, and NLRP3 pathways, are regarded feasible in drug development. The development of small molecules to treat COVID-19 has been achieved by several strategies, including computer-aided lead compound design and screening, natural product discovery, drug repurposing, and combination therapy. Several small molecules representative by remdesivir and paxlovid have been proved or authorized emergency use in many countries. And many candidates have entered clinical-trial stage. Nevertheless, due to the epidemiological features and variability issues of SARS-CoV-2, it is necessary to continue exploring novel strategies against COVID-19. This review discusses the current findings in the development of small molecules for COVID-19 treatment. Moreover, their detailed mechanism of action, chemical structures, and preclinical and clinical efficacies are discussed.

Sierra SARS-CoV-2 sequence and antiviral resistance analysis program

INTRODUCTION

Although most laboratories are capable of employing established protocols to perform full-genome SARS-CoV-2 sequencing, many are unable to assess sequence quality, select appropriate mutation-detection thresholds, or report on the potential clinical significance of mutations in the targets of antiviral therapy METHODS: We describe the technical aspects and benchmark the performance of Sierra SARS-CoV-2, a program designed to perform these functions on user-submitted FASTQ and FASTA sequence files and lists of Spike mutations. Sierra SARS-CoV-2 indicates which sequences contain an unexpectedly large number of unusual mutations and which mutations are associated with reduced susceptibility to clinical stage mAbs, the RdRP inhibitor remdesivir, or the Mpro inhibitor nirmatrelvir RESULTS: To assess the performance of Sierra SARS-CoV-2 on FASTQ files, we applied it to 600 representative FASTQ sequences and compared the results to the COVID-19 EDGE program. To assess its performance on FASTA files, we applied it to nearly one million representative FASTA sequences and compared the results to the GISAID mutation annotation. To assess its performance on mutations lists, we applied it to 13,578 distinct Spike RBD mutation patterns and showed that exactly or partially matching annotations were available for 88% of patterns CONCLUSION: Sierra SARS-CoV-2 leverages previously published data to improve the quality control of submitted viral genomic data and to provide functional annotation on the impact of mutations in the targets of antiviral SARS-CoV-2 therapy. The program can be found at https://covdb.stanford.edu/sierra/sars2/ and its source code at https://github.com/hivdb/sierra-sars2.

Paxlovid-like nirmatrelvir/ritonavir fails to block SARS-CoV-2 transmission in ferrets

Despite the continued spread of SARS-CoV-2 and emergence of variants of concern (VOC) that are capable of escaping preexisting immunity, therapeutic options are underutilized. In addition to preventing severe disease in high-risk patients, antivirals may contribute to interrupting transmission chains. The FDA has granted emergency use authorizations for two oral drugs, molnupiravir and paxlovid. Initial clinical trials suggested an efficacy advantage of paxlovid, giving it a standard-of-care-like status in the United States. However, recent retrospective clinical studies suggested a more comparable efficacy of both drugs in preventing complicated disease and case-fatalities in older adults. For a direct efficacy comparison under controlled conditions, we assessed potency of both drugs against SARS-CoV-2 in two relevant animal models; the Roborovski dwarf hamster model for severe COVID-19 in high-risk patients and the ferret model of upper respiratory tract disease and transmission. After infection of dwarf hamsters with VOC omicron, paxlovid and molnupiravir were efficacious in mitigating severe disease and preventing death. However, a pharmacokinetics-confirmed human equivalent dose of paxlovid did not significantly reduce shed SARS-CoV-2 titers in ferrets and failed to block virus transmission to untreated direct-contact ferrets, whereas transmission was fully suppressed in a group of animals treated with a human-equivalent dose of molnupiravir. Prophylactic administration of molnupiravir to uninfected ferrets in direct contact with infected animals blocked productive SARS-CoV-2 transmission, whereas all contacts treated with prophylactic paxlovid became infected. These data confirm retrospective reports of similar therapeutic benefit of both drugs for older adults, and reveal that treatment with molnupiravir, but not paxlovid, may be suitable to reduce the risk of SARS-CoV-2 transmission.