key: cord-0951892-hgnx9zwl
authors: Wang, Zhengyuan; Zhao, Yao; Wang, Qingxing; Xing, Yangfei; Feng, Lu; Kong, Juan; Peng, Chao; Zhang, Leike; Yang, Haitao; Lu, Min
title: Identification of proteasome and caspase inhibitors targeting SARS-CoV-2 M(pro)
date: 2021-06-01
journal: Signal Transduct Target Ther
DOI: 10.1038/s41392-021-00639-8
sha: 9f308ca5ae8aab74b4c3654fad38c966ec10af5a
doc_id: 951892
cord_uid: hgnx9zwl

nan

Since the beginning of 2020, the Coronavirus (CoV) Disease 2019 (COVID-19) pandemic has posed formidable challenges to public health security. The main protease (M pro , 3CL pro ) of CoVs plays essential roles in viral replication, making them attractive targets for antiviral drug development [1] [2] [3] . Dozens of SARS-CoV-2 M pro inhibitors have been reported with some entering clinical trials, but none is approved for COVID-19 treatment to date [1] [2] [3] . In this study, we discovered that the proteasome inhibitor MG132 and caspase inhibitors such as Z-VAD(OMe)-FMK are effective SARS-CoV-2 M pro inhibitors.

We recently identified arsenic trioxide as an effective mutant p53 rescue compound functioning by increasing protein thermostability 4 . Here we used protein thermostability as a readout to screen SARS-CoV-2 M pro -thermostabilizing compounds from a library of 4198 chemical entities containing US Food and Drug Administration (FDA)-approved drugs and clinical-stage or known-target compounds (Supplementary  Table S1 ). In the differential scanning fluorimetry (DSF) assay, Z-VAD(OMe)-FMK, MG132, boceprevir, thermopsine, and baicalein were identified as the top five SARS-CoV-2 M pro -thermostabilizing compounds, which increased the melting temperature (T m ) of M pro by at least 2°C (Fig. 1a , Supplementary  Fig. S1a , and Supplementary Table S1 ). The destabilizing hits ( Fig. 1a and Supplementary Fig. S1b) were not pursued further because of their potential promiscuity of binding or other undesirable properties. The top five M pro -stabilizing compounds were validated in concentration titration, whereby Z-VAD(OMe)-FMK was consistently the most potent M pro -stabilizing compound ( Supplementary Fig. S1c, d) .

We next predicted the M pro -binding modes of the five hits and performed a structure-activity relationship (SAR) analysis. The substrate-binding pocket of M pro , containing four subsites (S1′, S1, S2, and S4), is highly conserved among all CoV M pro homologs ( Supplementary Fig. S2a, upper panel) . We previously designed a series of M pro inhibitors harboring a classic core structure, in which the four subgroups (warhead, R1, R2, and R3) were assigned to occupy the four subsites of M pro , respectively ( Supplementary Fig. S2a , lower panel) 2 , as exemplified by N3 1 , 11a 2 , GC376 3 , and MI-23 5 ( Supplementary Fig. S2b ). Among the top five hits, baicalein and thermopsine are natural products with low molecular weight ( Supplementary Fig. S1d) , and thus they are unlikely to compactly occupy all four subsites. Boceprevir has been reported to function based on the classic core structure 3 , and its M pro -binding mechanism has been elucidated 5 . The remaining two hits, MG132 and Z-VAD(OMe)-FMK, structurally harbor a classic core structure and were thus proposed to bind to the active site of M pro in a conventional binding mode ( Supplementary Fig. S2c) . Notably, Z-VAD(OMe)-FMK lacks an R2 subgroup in the proposed binding mode ( Supplementary Fig. S2c ). To test the proposed binding mode, we collected 34 commercially available compounds sharing high structural similarity with MG132 and Z-VAD(OMe)-FMK ( Supplementary Fig. S2d ). Among them, three caspase inhibitors and three proteasome inhibitors detectably stabilized M pro (Fig. 1b, (3) small-sized R2 is apparently associated with high potency (comparison among the first 6 bars). In summary, the three most potent Z-VAD(OMe)-FMK, Z-DEVD-FMK, and Z-IETD-FMK contain a core structure with an FMK warhead and a small-sized R2.

To validate that M pro is a direct target of MG132 on the atomic level, we solved a high-resolution crystal structure of M pro -MG132 (Supplementary Table S2 ; 1.7 Å resolution). The two M pro molecules formed a symmetry homodimer ( Supplementary Fig.  S3a ). As proposed in Fig. S2c , MG132 binds to the active site based on the classic core structure ( Supplementary Fig. S3b) , whereby an aldehyde, side chain of Leu, side chain of the second Leu, and benzyloxycarbonyl (Cbz) act as the warhead, R1, R2, and R3, respectively (Fig. 1c) . Mass spectrometry (MS) did not reveal an obvious molecular weight increase of M pro upon incubation with MG132 ( Supplementary Fig. S3c ; ebselen was used as a control), which is presumably due to the highly dynamic and reversible bond between the MG132 aldehyde and the targeted cysteine. The mode of MG132 binding to SARS-CoV-2 M pro significantly differs from that of its binding to the proteasomal 20S subunit because of the different shape of the binding pocket ( Supplementary Fig. S3d ). Nevertheless, it is similar to those of the reported classical core-structure-based M pro inhibitors (Fig.  1d) . Notably, MG132 subgroups R1, R2, and R3 are all derived from hydrophobic Leu and Cbz, which undergo extensive hydrophobic interactions with M pro (Fig. 1e) . The core structure of MG132 forms four hydrogen bonds with M pro (comparable with the reported inhibitors), whereby the R1, R2, and R3 subgroups form only one hydrogen bond with M pro (Supplementary Fig. S3e) .

We next solved a high-resolution crystal structure of the M pro -Z-VAD(OMe)-FMK complex (Supplementary Table S2 ; 1.8 Å resolution), wherein the asymmetric unit contains one molecule ( Supplementary Fig. S4a ). Z-VAD(OMe)-FMK binds to the active site of M pro (Supplementary Fig. S4b) , however, in an unexpected binding mode (Fig. 1f) . The FMK group acts as a previously unreported Cys145-binding warhead, which was confirmed by MS ( Supplementary Fig. S4c , d; Z-DEVD-FMK and Z-IETD-FMK with the FMK group were also confirmed to covalently bind to Cys145).

FMK is frequently used as a warhead to tether cysteines of caspases. However, it has not been reported to be used in SARS-CoV-2 M pro inhibitors, to our knowledge. The binding mode of M pro -Z-VAD(OMe)-FMK differs significantly from that of caspase-1-Z-VAD-FMK, whereby the Z-VAD(OMe)-FMK molecule linearly occupies the long narrow pocket of caspase-1 ( Supplementary   Fig. S4e ). Due to R3 re-orientation, the binding mode of Z-VAD (OMe)-FMK to M pro is completely different from that of previously reported M pro inhibitors (Fig. 1g) . The four features of M pro -Z-VAD (OMe)-FMK binding (Fig. 1h) -FMK acting as the warhead, absence of R2, re-orientation of flexible hydrophobic R3, and compact occupation of the S2-S4 joint site-may contribute to the Letter high potency of Z-VAD(OMe)-FMK in thermostabilizing SARS-CoV-2 M pro . During preparation of this manuscript, a crystal structure of M pro -Z-VAD(OMe)-FMK (PDB ID: 7C8B) was released in the PDB database by an independent group, confirming the observed unconventional binding mode.

The inhibitory activities of the three caspase inhibitors and MG132 against SARS-CoV-2 M pro were determined in vitro using a fluorescence resonance energy transfer (FRET)-based assay. The three caspase inhibitors exhibited potent M pro inhibitory effects with half-maximal inhibitory concentrations (IC 50 ) in the nanomolar and low micromolar range (Fig. 1i, 0.59-2.80 μM) , while MG132 displayed a higher IC 50 (3.91 μM). Generally speaking, the three caspase inhibitors were comparable in their IC 50 values with the reported rationally designed M pro inhibitors ( Supplementary Fig.  S5a, 0.03-30.0 μM) .

The antiviral activities of these four compounds against SARS-CoV-2 were next determined in Vero cells. At 24 h following SARS-CoV-2 infection, the three caspase inhibitors displayed potent half-maximal effective concentrations (EC 50 ) in the nanomolar and low micromolar range in Vero cells (Fig. 1j, 0.64-1.88 μM) . We noticed that Z-VAD(OMe)-FMK exhibited higher potency in enzymatic activity assay, whereas lower potency in antiviral activity assay when compared to the other two caspase inhibitors. It may be associated with potential 'off-target' effect of Z-VAD (OMe)-FMK when inhibiting M pro in the cell-based antiviral assay, for example, promiscuously binding to multiple caspases and proteins of host cells. None of the tested caspase inhibitors caused cytotoxicity in Vero cells, exhibiting half-maximal cytotoxic concentrations (CC 50 ) > 300 μM (Fig. 1j) . MG132 is cytotoxic to this cell line 3 , and its EC 50 could not be reliably determined (Fig. 1j) . The extended cytotoxicity studies consistently suggested that these three caspase inhibitors were relatively non-toxic to cells (Supplementary Fig. S5b ). Compared to the anti-SARS-CoV-2 activities of the reported M pro inhibitors measured in the same Vero cells, the non-toxic caspase inhibitors were superior to those of N3 1 , Boceprevir 3 , and GC-376 3 , comparable with the six most potent M pro inhibitors (MI-09/12/14/28/30/31) reported by Qiao and colleagues 5 , but less potent than our previously optimized 11a/b 2 (Supplementary Fig. S5a) .

In summary, we identified MG132, Z-VAD(OMe)-FMK and its structural analogs as direct M pro -binding small molecules with antiviral activity against SARS-CoV-2. MG132 is widely recognized as a CoV inhibitor, with several reported and proposed targets. Our study provides evidence at single-atom resolution that M pro is a direct target of MG132. In contrast to the reported rationally designed M pro inhibitors, Z-VAD(OMe)-FMK binds to the active site of M pro in an unconventional binding mode, which may contribute to the observed high potency of Z-VAD(OMe)-FMK in inhibiting SARS-CoV-2. Our studies provide compelling structural evidence to support an alternative strategy for the design of potent inhibitors against SARS-CoV-2 M pro .

The datasets generated in this study are available from the corresponding authors upon reasonable request.

Atomic coordinates and structure factors of the reported crystal structures have been deposited with the worldwide Protein Data Bank (https://www.rcsb.org) under the PDB codes 7CUT and 7CUU. i The hydrolytic activity of SARS-CoV-2 M pro was measured in the presence of increasing concentrations of the indicated compounds. The IC 50 values were determined based on dose-response curves using nonlinear regression. The data represent mean ± SEM from 3 experiments. The y-axis of the graphs represents mean % inhibition of hydrolytic activity. j Antiviral activity against SARS-CoV-2 and cytotoxicity of the indicated compounds in Vero cells. Vero cells were treated with a concentration series of the indicated compounds and infected with SARS-CoV-2 (MOI = 0.01). At 24 h post-infection, the virions in the culture supernatant were harvested, followed by determination of the viral RNA copy number by qRT-PCR. The cytotoxicity of the indicated compounds in Vero cells was determined using the CCK8 assay. The y-axes of the graphs represent mean % reduction of the virion yield or mean % cell viability.

Zhengyuan Wang 1 , Yao Zhao 2 , Qingxing Wang 3 , Yangfei Xing 1 , Lu Feng 2 , Juan Kong 2 , Chao Peng 4 , Leike Zhang 3 , Haitao Yang 2 and Min Lu 1

Structure of M(pro) from SARS-CoV-2 and discovery of its inhibitors

Structure-based design of antiviral drug candidates targeting the SARS-CoV-2 main protease

GC-376, and calpain inhibitors II, XII inhibit SARS-CoV-2 viral replication by targeting the viral main protease

Arsenic trioxide rescues structural p53 mutations through a cryptic allosteric site

SARS-CoV-2 M(pro) inhibitors with antiviral activity in a transgenic mouse model

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.