Protease Inhibitors: Developing evolution-resistant drugs for COVID-19
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.
References
-
Dynamic profiling of β-coronavirus 3CL Mpro protease ligand-binding sitesJournal of Chemical Information and Modeling 61:3058–3073.https://doi.org/10.1021/acs.jcim.1c00449
-
Deep mutational scanning: a new style of protein scienceNature Methods 11:801–807.https://doi.org/10.1038/nmeth.3027
-
Combating susceptibility to drug resistance: lessons from HIV-1 proteaseChemistry & Biology 11:1333–1338.https://doi.org/10.1016/j.chembiol.2004.08.010
Article and author information
Author details
Publication history
- Version of Record published: July 26, 2022 (version 1)
Copyright
© 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.
Metrics
-
- 474
- Page views
-
- 145
- Downloads
-
- 1
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
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)
Further reading
-
- Evolutionary Biology
The globally invasive mosquito subspecies Aedes aegypti aegypti is an effective vector of human arboviruses, in part because it specializes in biting humans and breeding in human habitats. Recent work suggests that specialization first arose as an adaptation to long, hot dry seasons in the West African Sahel, where Ae. aegypti relies on human-stored water for breeding. Here, we use whole-genome cross-coalescent analysis to date the emergence of human-specialist populationsand thus further probe the climate hypothesis. Importantly, we take advantage of the known migration of specialists out of Africa during the Atlantic Slave Trade to calibrate the coalescent clock and thus obtain a more precise estimate of the older evolutionary event than would otherwise be possible. We find that human-specialist mosquitoes diverged rapidly from ecological generalists approximately 5000 years ago, at the end of the African Humid Period—a time when the Sahara dried and water stored by humans became a uniquely stable, aquatic niche in the Sahel. We also use population genomic analyses to date a previously observed influx of human-specialist alleles into major West African cities. The characteristic length of tracts of human-specialist ancestry present on a generalist genetic background in Kumasi and Ouagadougou suggests the change in behavior occurred during rapid urbanization over the last 20–40 years. Taken together, we show that the timing and ecological context of two previously observed shifts towards human biting in Ae. aegypti differ; climate was likely the original driver, but urbanization has become increasingly important in recent decades.
-
- Genetics and Genomics
- Evolutionary Biology
Brain size and cortical folding have increased and decreased recurrently during mammalian evolution. Identifying genetic elements whose sequence or functional properties co-evolve with these traits can provide unique information on evolutionary and developmental mechanisms. A good candidate for such a comparative approach is TRNP1, as it controls proliferation of neural progenitors in mice and ferrets. Here, we investigate the contribution of both regulatory and coding sequences of TRNP1 to brain size and cortical folding in over 30 mammals. We find that the rate of TRNP1 protein evolution (ω) significantly correlates with brain size, slightly less with cortical folding and much less with body size. This brain correlation is stronger than for >95% of random control proteins. This co-evolution is likely affecting TRNP1 activity, as we find that TRNP1 from species with larger brains and more cortical folding induce higher proliferation rates in neural stem cells. Furthermore, we compare the activity of putative cis-regulatory elements (CREs) of TRNP1 in a massively parallel reporter assay and identify one CRE that likely co-evolves with cortical folding in Old World monkeys and apes. Our analyses indicate that coding and regulatory changes that increased TRNP1 activity were positively selected either as a cause or a consequence of increases in brain size and cortical folding. They also provide an example how phylogenetic approaches can inform biological mechanisms, especially when combined with molecular phenotypes across several species.