Protease Inhibitors: Developing evolution-resistant drugs for COVID-19

Analyzing how mutations affect the main protease of SARS-CoV-2 may help researchers develop drugs that are effective against current and future variants of the virus.
  1. Daniel M Weinreich  Is a corresponding author
  1. Department of Ecology, Evolution and Organismal Biology, Brown University, United States
  2. Center for Computational Molecular Biology, Brown University, United States

As of mid-July 2022, more than 12 billion doses of vaccines against SARS-CoV-2, the virus responsible for COVID-19, have been administered to over 5 billion individuals (WHO, 2022). However, mutations that allow the virus to evade vaccines are spreading globally, leading to a need for boosters and updated vaccines (Callaway, 2021; Vogel, 2022). This immunity ‘arms race’ illustrates how efforts to control the COVID-19 pandemic could benefit from additional pharmaceutical approaches.

One of these approaches is the use of inhibitor drugs to block the action of viral proteases, the enzymes that cleave polyproteins encoded in the viral genome to yield several functional proteins. If these proteases are inactivated, the virus cannot synthesize the proteins it requires to reproduce, so inhibiting these enzymes may be an effective way of treating viral diseases. For example, protease inhibitors are key components of highly active antiretroviral therapy (HAART), which renders HIV a chronic disease.

In the case of SARS-CoV-2, many inhibitor candidates have been identified for its main protease, which is called Mpro (Jin et al., 2020; Narayanan et al., 2022). One of these inhibitors, with the trade name Paxlovid, has been authorized for emergency use by the United States Federal Drugs Administration (FDA), and pharmacists have been allowed to prescribe it since July 2022 (Owen et al., 2021; FDA, 2022). But can Mpro inhibitors be designed to avoid the rapid obsolescence we have observed in vaccines as a result of viral evolution? Now, in eLife, Julia Flynn, Daniel Bolon and colleagues at the University of Massachusetts (UMass) Chan Medical School and the Novartis Institutes for Biomedical Research report the results of experiments assessing the effects of different mutations on the proteolytic activity of MPro (Flynn et al., 2022).

Using a method called deep mutational scanning (Fowler and Fields, 2014), Flynn et al. generated a ‘library’ that contained almost every single missense mutation of the Mpro enzyme relative to the sequence found in the original Wuhan isolate (Wu et al., 2020). Missense mutations are changes to a single nucleotide in the sequence that alter the resulting amino acid. Assessing the activity of the MPro variants that result from each conceivable missense mutation revealed how they each affect the activity of the enzyme. To do this, Flynn et al. placed the library into yeast cells, and assessed how well each variant of the protease worked in three different environments.

The results showed that the proteolytic activity of each variant was correlated across the three assays, implying that the assays capture the same properties of MPro. Flynn et al. also found a correlation between the results of their assays and previous measurements of the catalytic rates for individual variants. Additionally, the vast majority of the 290 missense mutations of MPro observed most frequently (at least 100 times) in COVID-19 patients exhibited activity comparable to that of the wild type. This further confirmed that Flynn et al.’s assay captures the biologically relevant function of the enzyme. Finally, the team found that all of the mutation-intolerant residues in MPro (the sites at which at least 17 out of the 20 alternative amino acids block the protein’s activity) are highly conserved in other coronaviruses. This suggests that inhibitors designed against the current SARS-CoV-2 form of the enzyme will likely be active against future, emergent outbreaks.

Most protease inhibitors act by competitively occupying the enzyme’s binding pocket to the exclusion of its native substrate. Therefore, Flynn et al. decided to further analyze missense mutations in the vicinity of the binding pocket that nevertheless preserve enzyme activity, reasoning that these mutations might represent opportunities for the virus to evolve resistance to inhibitors. Unfortunately, they found that mutations at many of those residues preserve the proteolytic ability of Mpro (Cho et al., 2021).

Similarly, and perhaps of more immediate concern, Flynn et al. identified three missense mutations that reduced the binding stability to Paxlovid by at least 1 kcal/mol, corresponding to an approximately five-fold reduction in the drug’s binding affinity, while maintaining proteolytic activity. These represent likely resistance mutations; indeed, a subset of these mutations have now been observed in laboratory populations of SARS-CoV-2 challenged with this version of Paxlovid, as well as in currently circulating isolates (Service, 2022). It will be important to monitor the frequency of these three mutations in patients being treated with Paxlovid as its use becomes more widespread.

But the potential implications of the findings of Flynn et al. go further than providing insight into prospects for resistance evolution against existing Mpro inhibitors. Some years ago, the notion of the ‘substrate envelope’ of a viral protease – the overlapping volume occupied by all of its native substrates – was introduced (King et al., 2004). The inspiration was that the substrate specificity of the HIV protease seemed to be determined by the shape of the substrate rather than by a specific amino acid sequence. Because resistance mutations must allow an enzyme to distinguish between its inhibitors and its substrates, they might work by introducing atomic overlaps or disfavored electrostatics in the region where the inhibitor protrudes beyond the substrate envelope. In a recent companion study, researchers at UMass and Novartis (including many authors from Flynn et al.) have reported a high-resolution structure of the SARS-CoV-2 Mpro substrate envelope (Shaqra et al., 2022).

Together, these two articles provide clear guidance for the development of evolution-resistant protease inhibitors against SARS-CoV-2. To find clinically durable candidates, compounds should be screened first to identify those that lie entirely within the enzyme’s substrate envelope. Of these, care should be taken to select those that interact with one or more mutation-intolerant residues. This will make it more difficult for the protease to evolve the ability to discriminate between its native substrates and the drug without blocking its activity. Of course, Darwinian evolution is a tremendously powerful force, and resistance has evolved against every antimicrobial class ever deployed in the clinic. Nevertheless, the functional and structural resolution now available for Mpro-inhibitors offers new optimism 30 months into the COVID-19 pandemic.


  1. Website
    1. WHO
    (2022) WHO Coronavirus (COVID-19) Dashboard
    Accessed July 6, 2022.

Article and author information

Author details

  1. Daniel M Weinreich

    Daniel M Weinreich is in the Department of Ecology, Evolution and Organismal Biology, and the Center for Computational Molecular Biology, Brown University, Providence, United States

    For correspondence
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1424-7541

Publication history

  1. Version of Record published: July 26, 2022 (version 1)


© 2022, Weinreich

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.


  • 571
    Page views
  • 152
  • 1

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Daniel M Weinreich
Protease Inhibitors: Developing evolution-resistant drugs for COVID-19
eLife 11:e81334.

Further reading

    1. Developmental Biology
    2. Evolutionary Biology
    Kwi Shan Seah, Vinodkumar Saranathan
    Research Article

    The study of color patterns in the animal integument is a fundamental question in biology, with many lepidopteran species being exemplary models in this endeavor due to their relative simplicity and elegance. While significant advances have been made in unraveling the cellular and molecular basis of lepidopteran pigmentary coloration, the morphogenesis of wing scale nanostructures involved in structural color production is not well understood. Contemporary research on this topic largely focuses on a few nymphalid model taxa (e.g., Bicyclus, Heliconius), despite an overwhelming diversity in the hierarchical nanostructural organization of lepidopteran wing scales. Here, we present a time-resolved, comparative developmental study of hierarchical scale nanostructures in Parides eurimedes and five other papilionid species. Our results uphold the putative conserved role of F-actin bundles in acting as spacers between developing ridges, as previously documented in several nymphalid species. Interestingly, while ridges are developing in P. eurimedes, plasma membrane manifests irregular mesh-like crossribs characteristic of Papilionidae, which delineate the accretion of cuticle into rows of planar disks in between ridges. Once the ridges have grown, disintegrating F-actin bundles appear to reorganize into a network that supports the invagination of plasma membrane underlying the disks, subsequently forming an extruded honeycomb lattice. Our results uncover a previously undocumented role for F-actin in the morphogenesis of complex wing scale nanostructures, likely specific to Papilionidae.

    1. Evolutionary Biology
    Hironori Funabiki, Isabel E Wassing ... Thomas Carroll
    Research Article

    5-Methylcytosine (5mC) and DNA methyltransferases (DNMTs) are broadly conserved in eukaryotes but are also frequently lost during evolution. The mammalian SNF2 family ATPase HELLS and its plant ortholog DDM1 are critical for maintaining 5mC. Mutations in HELLS, its activator CDCA7, and the de novo DNA methyltransferase DNMT3B, cause immunodeficiency-centromeric instability-facial anomalies (ICF) syndrome, a genetic disorder associated with the loss of DNA methylation. We here examine the coevolution of CDCA7, HELLS and DNMTs. While DNMT3, the maintenance DNA methyltransferase DNMT1, HELLS, and CDCA7 are all highly conserved in vertebrates and green plants, they are frequently co-lost in other evolutionary clades. The presence-absence patterns of these genes are not random; almost all CDCA7 harboring eukaryote species also have HELLS and DNMT1 (or another maintenance methyltransferase, DNMT5). Coevolution of presence-absence patterns (CoPAP) analysis in Ecdysozoa further indicates coevolutionary linkages among CDCA7, HELLS, DNMT1 and its activator UHRF1. We hypothesize that CDCA7 becomes dispensable in species that lost HELLS or DNA methylation, and/or the loss of CDCA7 triggers the replacement of DNA methylation by other chromatin regulation mechanisms. Our study suggests that a unique specialized role of CDCA7 in HELLS-dependent DNA methylation maintenance is broadly inherited from the last eukaryotic common ancestor.