Dear Editor,
The ongoing SARS-CoV-2 pandemic continues to be a significant threat to global health. First reported in November 2021, the Omicron variant (B.1.1.529) is more transmissible and can evade immunity better than previous SARS-CoV-2 variants, fueling an unprecedented surge in cases. To produce functional proteins from its polyprotein, SARS-CoV-2 relies on the cysteine proteases Nsp3/papain-like protease (PLpro) and Nsp5/main protease (Mpro)/3C-like protease to cleave at three and more than 11 sites, respectively.1 Therefore, Mpro and PLpro inhibitors are considered to be one of the most promising SARS-CoV-2 antivirals. On December 22, 2021, the Food and Drug Administration (FDA) issued an Emergency Use Authorization (EUA) for PAXLOVID, a ritonavir-boosted formulation of nirmatrelvir. Nirmatrelvir is a first-in-class orally bioavailable SARS-CoV-2 Mpro inhibitor.2 Thus, the scientific community must vigilantly monitor potential mechanisms of drug resistance, especially because SARS-CoV-2 is naïve to Mpro inhibitors. Mutations have been well identified in variants to this point.3 Notably, Omicron Mpro (OMpro) harbors a single mutation—P132H. In this study, we characterized the enzymatic activity, drug inhibition, and structure of OMpro while evaluating the past and future implications of Mpro mutations.
Using an established Förster resonance energy transfer (FRET) assay to monitor proteolytic activity,4 we show that OMpro and wild-type (WT) Mpro have equivalent affinity and catalytic constants for its substrate (Fig. 1a, b). However, subsequent melting temperature experiments reveal OMpro has a melting temperature (Tm) of 53.6 ± 0.1 °C; lower than the WT Mpro Tm of 56.2 ± 0.2 by 2.6 °C (Fig. 1b). Despite its lower melting temperature, OMpro and WT Mpro hydrolyze their substrates at comparable rates for up to 24 h at 37 °C (Fig. 1c). Further biochemical analysis suggests that OMpro and Mpro are equally susceptible to covalent inhibitors such as GC-376, PF-07321332 (nirmatrelvir), and PF-00835231 (Fig. 1d). Interestingly, when these inhibitors are incubated with WT Mpro or OMpro, their melting temperatures converge, despite the different melting temperature of their apo forms (Fig. 1e). This suggests that inhibitors may stabilize OMpro to a greater extent than WT Mpro. These described biochemical/biophysical results are in line with the equivalent antiviral potency of these molecules against WT and Omicron and the lack of observable differences in the active site.5,6,7,8
In parallel to our study, Ullrich et al. also characterized six Mpro mutants identified from the circulating SARS-CoV-2 variants.9 The catalytic efficiency (kcat/Km) values for WT and Omicron Mpro were 0.016 and 0.023 s−1 µM−1, similar to the values obtained in our study (Fig. 1b). In addition, the Pfizer team reported that PF-07321332 (nirmatrelvir) retained a potent inhibitory constant (Ki) against Omicron P132H Mpro (Ki = 0.635 nM), similar to WT (Ki = 0.933 nM).10 Our independent study further confirmed these results.
The crystal structure of OMpro in complex with GC-376 was determined at 2.05 Å resolution in the monoclinic space group I2 with a Rwork/Rfree of 0.179/0.219 (Supplementary information, Table S1; PDB ID 7TOB). The unit cell is the same as many previously solved WT Mpro crystal structures where a = 45.19 Å, b = 52.99 Å, c = 113.01 Å and α = 90.00°, β = 100.50°, γ = 90.00°. Overall, OMpro has a nearly identical structure to WT Mpro (Fig. 1f).4 The most pronounced differences involve the area around the mutation associated with OMpro, P132H. Found 22 Å from the catalytic cysteine Cys145, P132H lies between the catalytic domain and the dimerization domain, and thus, it does not impart any direct structural changes to the active site.
Clear electron density shows that His132 forms π-stacking interactions with the sidechain of Glu240 (Fig. 1g). As the imidazole sidechain of histidine is protonated near physiological pH, the interaction with Glu240 may be further strengthened through electrostatic interactions. Additionally, there is a newly formed water-mediated hydrogen bond with Glu240. We also find that Glu240 reorients itself toward the core to accommodate the bulkier His sidechain, where it forms a new hydrogen bond with Thr198 (Fig. 1h). Consequently, the Thr198 sidechain rotates ~90°, placing its hydroxyl group in close distance to Glu240. As a result of these conformational changes, certain portions of the enzyme appear to move very slightly, ~0.5 Å (Fig. 1f). We speculate that the lower thermal stability of apo OMpro may be due to minor residue adjustments to accommodate the bulkier His132 sidechain, ultimately destabilizing its structure. Because the mutation occurs at the interface between the dimerization domain and the catalytic domain, these movements can also affect intramolecular packing. Residue 132 is also located at the turn between two β-sheets, a position that naturally favors the cyclic sidechain of proline.
Although the P132H mutation does not appear to reduce enzymatic activity and inhibitor binding (Fig. 1a–d), the decrease in thermal stability (Fig. 1b) indicates that protein flexibility may be greater, which plays an important role in enzyme evolution, especially to broaden substrate profile or to alter ligand binding. Mpro can recognize a wide range of peptide substrates, although its P1 position preferentially binds to glutamine. Future studies on whether P132H and other mutations can influence enzymatic activity for larger substrate libraries or other known ligands will be necessary. In addition, the effect of the Mpro P132H mutant on the viral replication remains elusive, and further works need to generate recombinant viruses to study the viral replication kinetics and tropism.
Extensive sequencing of SARS-CoV-2 isolates has provided unprecedented insights into the mutations that occur in the viral RNA genome.3 Based on the annotations provided through CoVsurver enabled by GISAID (www.gisaid.org/epiflu-applications/covsurver-mutations-app), mutations of Nsp5 appear mostly stochastic (Supplementary information, Fig. S1); however, several positions have a disproportionately large number of mutations. Of the top 25 most common mutants (Fig. 1i), three are found on P132: P132H (489,444 occurrences; EPI_ISL_8931050), P132L (8813 occurrences; EPI_ISL_8768027), and P132S (7452 occurrences; EPI_ISL_8925342). L89F, K90R, and K88R are the second, third, and seventh, most common mutations, suggesting this β-sheet is also a hotspot. Notably, none of the 25 most common mutants involve residues in the active site or at the dimerization interface (Fig. 1j). With several exceptions, including P132H, the resulting amino acid is often similar in size and physicochemical properties, such as K→R. However, SARS-CoV-2 has not yet encountered Mpro antivirals. If nothing else, SARS-CoV-2 has taught us that widespread proliferation, low fidelity genome synthesis, and selective pressure might quickly produce drug-resistant phenotypes. Thus, it is important to monitor future mutations and their associated biochemical properties to anticipate future drug-resistance.
Data availability
The X-ray crystal structure of the Omicron Mpro P132H mutant in complex with GC376 was deposited in PDB with the code 7TOB.
References
Cannalire, R. et al. J. Med. Chem. https://doi.org/10.1021/acs.jmedchem.0c01140 (2020).
Owen, D. R. et al. Science 374, 1586–1593 (2021).
WHO. https://www.who.int/en/activities/tracking-SARS-CoV-2-variants Accessed on 19 February 2022. (2022).
Ma, C. et al. Cell Res. 30, 678–692 (2020).
Dabrowska, A. et al. bioRxiv https://doi.org/10.1101/2021.12.21.473268 (2021).
Rai, D. K. et al. bioRxiv https://doi.org/10.1101/2022.01.17.476644 (2022).
Vangeel, L. et al. Antiviral Res. 198, 105252 (2022).
Rosales, R. et al. bioRxiv https://doi.org/10.1101/2022.01.17.476685 (2022).
Ullrich, S. et al. Boorg. Med. Chem. Lett. https://doi.org/10.1016/j.bmcl.2022.128629 (2022).
Greasley, S. E. et al. bioRxiv https://doi.org/10.1101/2022.01.17.476556 (2022).
Acknowledgements
X-ray data were collected at the X-ray facility supported by the Spallation Neutron Source (Oak Ridge National Laboratory), a DOE Office of Science User Facility. This research was supported by the National Institutes of Health (NIH) (grants AI147325, AI157046, and AI158775) to J.W. We thank Eric M Lewandowski for reviewing this paper.
Author information
Authors and Affiliations
Contributions
M.D.S., J.W., and Y.C. conceived the experiments, interpreted the data, and composed the paper. Y.H. performed the biochemical assays and melting temperature experiments. Proteins were cloned and purified by X.Z. Crystal structures were solved and refined by M.D.S. Crystals were grown by M.V.G. X-ray datasets were prepared and collected by M.T.K and F.M.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Supplementary information
Rights and permissions
About this article
Cite this article
Sacco, M.D., Hu, Y., Gongora, M.V. et al. The P132H mutation in the main protease of Omicron SARS-CoV-2 decreases thermal stability without compromising catalysis or small-molecule drug inhibition. Cell Res 32, 498–500 (2022). https://doi.org/10.1038/s41422-022-00640-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41422-022-00640-y
This article is cited by
-
Discovery and evaluation of active compounds from Xuanfei Baidu formula against COVID-19 via SARS-CoV-2 Mpro
Chinese Medicine (2023)
-
Multiple pathways for SARS-CoV-2 resistance to nirmatrelvir
Nature (2023)
-
SARS-CoV-2: analysis of the effects of mutations in non-structural proteins
Archives of Virology (2023)
-
Rapid assembly of SARS-CoV-2 genomes reveals attenuation of the Omicron BA.1 variant through NSP6
Nature Communications (2023)
-
Ginkgolic acid and anacardic acid are reversible inhibitors of SARS-CoV-2 3-chymotrypsin-like protease
Cell & Bioscience (2022)