key: cord-033951-77tfhm5b
authors: Ma, Chunlong; Hu, Yanmei; Townsend, Julia Alma; Lagarias, Panagiotis I.; Marty, Michael Thomas; Kolocouris, Antonios; Wang, Jun
title: Ebselen, Disulfiram, Carmofur, PX-12, Tideglusib, and Shikonin Are Nonspecific Promiscuous SARS-CoV-2 Main Protease Inhibitors
date: 2020-10-09
journal: ACS Pharmacol Transl Sci
DOI: 10.1021/acsptsci.0c00130
sha: 
doc_id: 33951
cord_uid: 77tfhm5b

[Image: see text] Among the drug targets being investigated for SARS-CoV-2, the viral main protease (M(pro)) is one of the most extensively studied. M(pro) is a cysteine protease that hydrolyzes the viral polyprotein at more than 11 sites. It is highly conserved and has a unique substrate preference for glutamine in the P1 position. Therefore, M(pro) inhibitors are expected to have broad-spectrum antiviral activity and a high selectivity index. Structurally diverse compounds have been reported as M(pro) inhibitors. In this study, we investigated the mechanism of action of six previously reported M(pro) inhibitors, ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12, using a consortium of techniques including FRET-based enzymatic assay, thermal shift assay, native mass spectrometry, cellular antiviral assays, and molecular dynamics simulations. Collectively, the results showed that the inhibition of M(pro) by these six compounds is nonspecific and that the inhibition is abolished or greatly reduced with the addition of reducing reagent 1,4-dithiothreitol (DTT). Without DTT, these six compounds inhibit not only M(pro) but also a panel of viral cysteine proteases including SARS-CoV-2 papain-like protease and 2A(pro) and 3C(pro) from enterovirus A71 (EV-A71) and EV-D68. However, none of the compounds inhibits the viral replication of EV-A71 or EV-D68, suggesting that the enzymatic inhibition potency IC(50) values obtained in the absence of DTT cannot be used to faithfully predict their cellular antiviral activity. Overall, we provide compelling evidence suggesting that these six compounds are nonspecific SARS-CoV-2 M(pro) inhibitors and urge the scientific community to be stringent with hit validation.

A new coronavirus, SARS-CoV-2, started to circulate among humans in late 2019 and quickly evolved to be a global pandemic. As of August 31, 2020, there have been more than 6 million positive cases with over 183 000 deaths in the United States alone. The origin of SARS-CoV-2 is still under investigation, and the closest strain is the bat coronavirus RaTG13, which shares 96% sequence similarity with SARS-CoV-2. 1 It is unknown whether SARS-CoV-2 transmitted directly from bats to humans or through an intermediate host. 2 There are currently no vaccines or antiviral drugs available for SARS-CoV-2. Encouraging progress has been made in vaccine development, and several RNA-, DNA-, and adenovirus-based vaccine candidates are now in phase III clinical trials. 3 For small-molecule antivirals, remdesivir was granted emergency use authorization in the United States.

SARS-CoV-2 is an enveloped, positive-sense, single-stranded RNA virus that belongs to the betacoronavirus genera, which also includes SARS-CoV, MERS-CoV, HCoV-OC43, and HCoV-HKU1. SARS-CoV-2 shares ∼80% sequence similarity with SARS-CoV. As such, many of the reported antivirals against SARS-CoV-2 were originally developed for either SARS-CoV or other related coronaviruses. 4 SARS-CoV-2 infects ACE2-expressing cells and enters the cell through either direct cell surface fusion or the endosomal pathway. 5 For direct cell surface fusion, the host membrane protease TMPRSS2 cleaves the viral spike protein, triggering viral membrane fusion with the host cell membrane. 6 For endosomal entry, cathepsin L mediates the cleavage of viral spike protein. 7 Once the viral RNA is released in the cytoplasm, it undergoes translation into viral polyproteins pp1a and pp1ab, which are subsequently cleaved by two viral proteases, the main protease (M pro ), also called 3-chymotrypsin-like protease (3CL pro ), and the papain-like protease (PL pro ). The released viral proteins can then assemble to form the viral polymerase RdRp complex to catalyze the replication of viral RNA. Finally, progeny virions are released from the infected cells through exocytosis and are ready for the next round of infection.

Compounds that interfere with any step in the viral life cycle can theoretically inhibit viral replication. Among the list of drug targets pursued as SARS-CoV-2 small-molecule antivirals, the viral polymerase RdRp and the protease M pro are the most extensively studied. The RdRp inhibitor remdesivir received emergency use authorization in the United States. It showed broad-spectrum antiviral activity against multiple coronaviruses in cell culture, including SARS-CoV-2, SARS-CoV, and MERS-CoV, and it also had in vivo efficacy in a SARS-CoV infection mouse model. 8 EIDD-2801, an RdRp inhibitor, is another promising drug candidate with broad-spectrum antiviral activity. 9 M pro is a cysteine protease that cleaves the viral polyprotein at more than 11 sites. It has a unique substrate preference of glutamine at the P1 position, while no host protease is known to have such a preference. 10 As such, the most potent M pro inhibitors such as GC376 and N3 all contain a 2-pyrrolidone substitution in the P1 position as a mimetic of the glutamine in the substrate. Several crystal structures of M pro in complex with inhibitors have been solved, showing that the pyrrolidone forms multiple hydrogen bonds with the His163 and Glu166 side chains and the main chain of Phe140. 11 −15 In addition to the classic pyrrolidone-containing M pro inhibitors, several noncanonical M pro inhibitors have also been reported with both enzymatic inhibition and cellular antiviral activity. 11, 12 In this study, we aim to validate six previously reported M pro inhibitors: ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12 ( Figure 1 ). 12 Among these six compounds, ebselen is a clinical candidate with antiinflammatory and antioxidant activities. In preclinical studies, ebselen was reported to react with cysteine residues from completely unrelated proteins including the C-terminal domain of the HIV-1 capsid, 16 Mycobacterium tuberculosis transpeptidase LdtMt2, 17 glutamate dehydrogenase, 18 Clostridium dif f icile toxins TcdA and TcdB, 19 Mycobacterium tuberculosis (Mtb) antigen 85C enzyme, 20 hepatitis C virus NS3 helicase, 21 plant cysteine protease papain, 22 glutathione Stransferases, 22 and many others. However, the inhibition of papain by ebselen was abolished by the addition of reducing regents including glutathione (GSH), 2-mercaptoethanol, and sodium borohydride. 22 Ebselen was also reported to induce protein unfolding for the insulin-degrading enzyme. 23 Ebselen analogues were synthesized and were found to inhibit both SARS-CoV-2 M pro and PL pro . 24 Disulfiram inhibits a panel of diverse enzymes including methyltransferase, 25 urease, 26 and kinase 27 through reaction with the cysteine residues. A study also showed that disulfiram inhibits PL pro from SARS-CoV and MERS-CoV with IC 50 values of 14.2 and 22.7 μM, respectively. 28 However, the inhibition was completely lost in the presence of 5 mM β-mercaptoethanol (IC 50 > 300 μM). 28 Carmofur inhibits human acid ceramidase by covalently modifying the catalytic C143 residue. 29 PX-12 inhibits tubulin polymerization through cysteine oxidation. 30 Tideglusib is an irreversible inhibitor of glycogen synthase kinase-3β (GSK-3β). 31 Ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin were recently reported as SARS-CoV-2 M pro inhibitors with IC 50 values ranging from 0.67 to 21.39 μM in the FRET-based enzymatic assay. 12 Among the six compounds, ebselen inhibited SARS-CoV-2 replication with an EC 50 value of 4.67 ± 0.80 μM in plaque reduction assay. Disulfiram was able to reduce the viral replication by ∼30% in the viral titer reduction assay at 10 μM, while tideglusib, carmofur, and PX-12 had no significant antiviral effect. Intriguingly, in a follow up study, carmofur was shown to inhibit SARS-CoV-2 viral replication with an EC 50 of 24.30 μM, and the X-ray crystal structure of M pro with carmofur was solved (PDB: 7BUY). 32 The description of the enzymatic assay did not specify whether the reducing reagent dithiothreitol (DTT) was added or not. It is standard practice to add DTT or another reducing reagent such as glutathione (GSH) or β-mercaptoethanol (β-ME) in the enzymatic assay of cysteine protease to ensure that the enzyme is in the active form by reducing the catalytic cysteine residue as well as preventing nonspecific thiol reactive compounds from covalently modifying the catalytic cysteine.

The questions we are trying to address in this study are whether the inhibition of M pro by these compounds is specific and whether their enzymatic inhibition potency IC 50 values can be used to faithfully predict the cellular antiviral activity. In other words, do the reported IC 50 values of ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12 against SARS-CoV-2 M pro reflect specific enzymatic inhibition, or are they due to nonspecific inactivation of the enzyme? We tested these compounds against a panel of related and unrelated viral cysteine proteases, the SARS-CoV-2 PL pro , and the 2A protease (2A pro ) and 3C protease (3C pro ) from EV-A71 and EV-D68, in ACS Pharmacology & Translational Science pubs.acs.org/ptsci Article a consortium of assays with or without DTT. Collectively, our results showed that in the absence of DTT, ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12 nonspecifically inhibit all six viral cysteine proteases including SARS-CoV-2 M pro . However, despite their potent inhibition of enzymatic activity of 2A pro and 3C pro from EV-A71 and EV-D68 in the FRET assay in the absence of DTT, none of the compounds showed cellular antiviral activity against EV-A71 and EV-D68. Therefore, it can be concluded that the enzymatic inhibition potency of cysteine protease inhibitors measured in the absence of DTT cannot be used to predict the cellular antiviral activity. Overall, although there is an immediate need for SARS-CoV-2 antivirals, the scientific community needs to be cautious about the nonspecific effect of promiscuous compounds, and secondary assays should be performed at the early stage to triage hits that lack translational potential.

Inhibition of SARS-CoV-2 M pro , PL pro , EV-A71 and EV-D68 2A pro , and 3C pro by Ebselen, Disulfiram, Carmofur, PX-12, Tideglusib, and Shikonin Is DTT-Dependent. To dissect the effect of DTT on the enzymatic inhibition of SARS-CoV-2 M pro by ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin, we performed dose−response titration in the FRET-based enzymatic assay with and without 4 mM DTT. A known M pro inhibitor, GC376, was included as a positive control. It was found that all compounds inhibited M pro in the absence of DTT (Figure 2A , red curves), and the IC 50 values are generally in agreement with previous published results, 12 except for those of shikonin and PX-12, which showed more than 10-fold difference. However, none of the compounds showed potent inhibition against M pro in the presence of 4 mM DTT (IC 50 > 25 μM) (Figure 2A , black curves). Carmofur showed weak inhibition with an IC 50 value of 28.2 ± 9.5 μM in Table 1 ). These results suggest that the claimed inhibition of M pro by these six compounds might not be target-specific. To test this hypothesis, we next tested these six compounds against five other viral cysteine proteases, among which SARS-CoV-2 PL pro , EV-A71 2A pro , and EV-D68 2A pro have no sequence similarity with M pro , while EV-A71 3C pro and EV-D68 3C pro share similar chymotrypsin-like folding with M pro , despite only showing 16.3 and 19.6% sequence similarities with M pro . GRL0617 was included as a positive control for SARS-CoV-2 PL pro , 33 and telaprevir was included as a positive control for EV-A71 2A pro and EV-D68 2A pro . 34 GC376 was used as a positive control for both EV-A71 3C pro and EV-D68 3C pro . If the inhibition of M pro by ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin is specific, in either the absence or the presence of DTT, then these six compounds are expected to show little or no inhibition against the unrelated SARS-CoV-2 PL pro , EV-A71 2A pro , and EV-D68 2A pro . For SARS-CoV-2 PL pro , we observed results similar to those for SARS-CoV-2 M pro : All compounds displayed a potent inhibitory effect in the absence of DTT, while little or no inhibition was observed in the presence of DTT ( Figure 2B ). Shikonin was less potent against PL pro than against M pro in the absence of DTT with IC 50 values of 1.5 and 55.3 μM, respectively ( Table 1 ). The potency of shikonin increased about 2-fold, and the IC 50 values were 55.3 and 28.2 μM with and without DTT, respectively ( Table 1) . As expected, the noncovalent SARS-CoV-2 PL pro inhibitor GRL0617 35 inhibits SARS-CoV-2 PL pro similarly in the present or the absence of 4 mM DTT ( Figure 2B , last column). For EV-A71 and EV-D68 2A pro and 3C pro , we observed trends similar to that of M pro : All six compounds showed potent enzymatic inhibition in the absence of DTT, and the inhibition was abolished with the addition of DTT (Figures 2C−F). In contrast, there is no significant shift of the potency with and without DTT for GC376 in inhibiting EV-A71 3C pro ( Figure 2D ) and EV-D68 3C pro ( Figure 2F ) and for telaprevir in inhibiting EV-A71 2A pro ( Figure 2C ) and EV-D68 2A pro ( Figure 2E ).

Next, we tested whether another reducing agent, GSH, could also abolish the inhibitory effect of these promiscuous inhibitors against SARS-CoV-2 M pro . Ebselen and disulfiram were chosen as representative examples. In the absence of 4 mM DTT or 1 mM GSH, ebselen and disulfiram completely inhibit M pro enzymatic activity at 20 μM ( Figure 3 ); however, no inhibition was observed for ebselen and disulfiram when either 4 mM DTT or 1 mM GSH was present in the reaction buffer ( Figure 3 ). In contrast, the inhibition by the positive control GC376 was not affected by the reducing agent DTT or GSH.

Collectively, the enzymatic assay results suggest that ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin are promiscuous cysteine protease inhibitors that inhibit not only M pro but also five other related and unrelated viral cysteine proteases including SARS-CoV-2 PL pro and EV-A71 and EV-D68 2A pro and 3C pro in the absence of DTT, and the 

The values shown in bold were reported in reference 12 N.T. = not tested. inhibition is abolished with the addition of a reducing reagent, either DTT or GSH.

Ebselen, Disulfiram, Carmofur, PX-12, Tideglusib, and Shikonin Did Not Bind to SARS-CoV-2 M pro in the Presence of DTT in the Thermal Shift Assay. To investigate whether ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin bind directly to M pro or other related and unrelated cysteine proteases, we performed a thermal shift binding assay. In the thermal shift binding assay, a temperature gradient is applied to denature a protein in the presence of a fluorescence dye. When the protein unfolds, the hydrophobic region is exposed to the fluorescence dye, and an increased fluoresce signal is observed. Specific binding of a small molecule to the native state of a protein usually stabilizes the protein, leading to a shift of the melting temperature (ΔT m ). 34,36−38 Here we measured the T m change upon addition of these six compounds in the presence or absence of 4 mM DTT against six viral cysteine proteases including M pro . All compounds were tested at 40 μM except shikonin, which was tested at 10 μM as it quenches the SYPRO orange dye fluorescence signal at 40 μM. Compounds were preincubated with 3 μM protease in its corresponding enzymatic reaction buffer at 30°C for 30 min with or without 4 mM DTT, then a 20−90°C temperature gradient was applied and T m was calculated. As expected, for the positive controls in the presence of 4 mM DTT, a significant T m increase was observed in the binding of GC376 to SARS-CoV-2 M pro ( Figure 4A ), EV-A71 3C pro ( Figure 4D ), and EV-D68− 3C pro ( Figure 4F ). Binding of GRL0617 to SARS-CoV-2 PL pro also led to significant stabilization ( Figure 4B Taken together, when 4 mM DTT was present in the assay buffer, there was no binding between the six viral cysteine proteases and the six compounds ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin; without DTT, these compounds appear to nonspecifically bind to these proteases, leading to destabilization.

Ebselen, Disulfiram, PX-12, Tideglusib, and Shikonin Did Not Bind to SARS-CoV-2 M pro , while Carmofur Showed Binding in the Presence of DTT in the Native Mass Spectrometry Binding Assay. To corroborate the results from the thermal shift binding assay, we next performed native mass spectrometry (MS)-based binding assays. Native MS analysis revealed that SARS-CoV-2 M pro forms a dimer that was measured to have a mass of 67 595 Da ( Figure 5A ). A small abundance of monomer was measured with a mass of 33 796 Da, but the intact dimer was the predominant signal (data not shown). The addition of GC376 revealed that up to two ligands bound per dimer (Figure 5B) , suggesting a binding ratio of one drug per monomer, which is consistent with its mechanism of action revealed by X-ray crystallography. 11 The addition of 4 mM DTT shifted the equilibrium to one ligand per dimer ( Figure 5C ). When carmofur was added to SARS-CoV-2 M pro , it bound up to three ligands per dimer, with the most abundant signal being that for the two bound per dimer ( Figure 5J ). When 4 mM DTT was added, it disrupted the ligand binding of carmofur, with the most abundant signal being that for the dimer without ligand bound ( Figure 5K ). Nevertheless, signals corresponding to one ligand per dimer and two ligands per dimer could still be detected, suggesting that carmofur has moderate binding toward SARS-CoV-2 M pro Figure 5K ). This result is consistent with our enzymatic assay result, which showed that carmofur inhibits M pro with an IC 50 of 28.2 ± 9.5 μM in the presence of 4 mM DTT (Figure 2A) . Taken together, the inhibition of M pro by carmofur has certain degree of specificity, although the potency is moderate. Similarly, disulfiram, ebselen, and PX-12 bound up to three, four, and five ligands per dimer respectively ( Figure 5F ,D,N) in the absence of DTT, and the addition of 4 mM DTT completely disrupted this ligand binding ( Figure 5G ,E,O), and only the dimer signal without ligand was detected. The complete disruption of this ligand binding with the addition of DTT suggests that these compounds bind nonspecifically to SARS-CoV-2 M pro . Shikonin was found to bind up to four ligands per dimer to SARS-CoV-2 M pro ( Figure 5L ), and this binding was completely disrupted upon addition of 4 mM DTT ( Figure  5M ). Tideglusib did not bind to M pro in either the absence or the presence of 4 mM DTT at concentrations ranging from 10 to 40 μM ( Figure 5H ,I).

The observation that ebselen, disulfiram, carmofur, shikonin, and PX-12 can bind to M pro with more than two ligands per dimer in the absence of DTT indicates that these compounds might not only modify the catalytic cysteine C145 but also possibly bind to allosteric sites or covalently modify other cysteine residues on M pro .

Molecular Dynamics (MD) Simulations of the Binding of M pro to GC376, Carmofur, and Ebselen. We performed MD simulations to compare the stability of the binding interactions identified in the X-ray structures of SARS-CoV-2 M pro in complex with GC376 and carmofur. The MD simulations of M pro with ebselen were carried using the highest scored docking pose.

Our previous study showed that when GC376 binds to the SARS-CoV-2 M pro , the covalent thioketal adduct can adapt both the S-and R-configuration. 11 The MD simulations showed that the complexes formed with GC376 in either the Sor R-configuration did not deviate from the starting X-ray structures with an RMSD in protein and ligand positions from the X-ray structure smaller than ca. 2 Å for the protein and smaller than ca. 2.3 Å for the ligand (Figure 6C ,F). The MD simulations verified stabilizing interactions observed in the Xray structure, which remain stable inside the binding cavity as shown in the frequency interaction plot in Figure 6A ,D. From the protein−ligand contact plots of the GC376 in the Sconfiguration ( Figure 6A,B) , it is shown that it forms multiple hydrogen bonds, i.e., (a) between the thiohemiketal P1 hydroxyl group and C145 peptidic NH; (b) between 2pyrrolidone's NH at the P1 and the side chain of E166 and between peptidic NH at P1 with H164 side chain imidazole; (c) between carbamate CO, NH, and benzyloxy oxygen at P2 with peptidic NH at M165, Q189 side chain CO, and Q189 side chain NH, respectively. In the X-ray structure, the pyrrolidone's CO forms a hydrogen bond with H163 side chain imidazole; the MD simulation plot in Figure 6A represents an average description of the interactions from an ensemble and not from a single snapshot in the X-ray structure. The polar 2-pyrrolidone group is oriented toward the solventexposed S1 pocket, while isobutyl group at P2 position is oriented comfortably toward the hydrophobic S2 site formed by H41, M49, and M169. The benzyloxy group facing toward L167 moves freely in the area between Q189, A191, Q192, L167, and P168 ( Figure 6B ). In the R-configuration, GC376 is stabilized through hydrogen bonding interaction with C145, H164, E166, and Q189, as well as G143 and H41, but the frequency of hydrogen bond interactions is less when pyrrolidone side chain is buried in the S1 pocket ( Figure  6D,E) .

In contrast to the numerous stabilizing interactions of GC376 in both the S-and R-configurations, carmofur interacts mainly with C145 through van der Waals and hydrogen bonding interactions ( Figure 6G,H) . In carmofur, the MD simulations show that the RMSD in protein and ligand positions from the X-ray structure are both ca. 3 Å for the protein and the ligand. The hexyl side chain is oriented inside the binding cavity, from S1 to S3, but the interactions of the drugs cannot be specific for M pro without directing hydrogen bond and van der Waals interactions complementary to the cavity. Ebselen interacts through hydrogen bonding interactions mainly with E166 and van der Waals interactions with M49 and M165 in the hydrophobic S2 binding area. Similar to carmofur, these interactions are not adequate to effectively trap the small drug inside the wide binding area of M pro ( Figure  6J,K) , and the drug can rotate around the phenyl−CO bond, resulting in a high RMSD of ca. 5.4 Å ( Figure 6L ). Compared to GC376, the protein−ligand contact plots of carmofur and ebselen suggest that these two compounds bind to SARS-CoV-2 M pro with reduced affinity, corroborating with their nonspecific inhibition mechanism.

Ebselen, Disulfiram, Carmofur, PX-12, Tideglusib, and Shikonin Had No Cellular Antiviral Activity against EV-A71 and EV-D68. If the enzymatic inhibition potency IC 50 values obtained in the absence of DTT can be used to faithfully predict the cellular antiviral activity, then one would expect all six compounds, ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin, will have potent antiviral activity against EV-A71 and EV-D68. To test this hypothesis, the cellular antiviral activity of ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin against EV-A71 and EV-D68 viruses were tested in RD cells using the viral cytopathic effect (CPE) assay. 34 GC376 and telaprevir were included as positive controls as 3C pro and 2A pro inhibitors. GC376 inhibited EV-A71 and EV-D68 with EC 50 values of 0.2 and 0.9 μM, respectively ( Table 2) . Telaprevir inhibited EV-D68 with an EC 50 of 0.4 μM (Table 2 ). However, none of the six compounds showed antiviral activity against either EV-A71 or EV-D68 at the highest nontoxic drug concentration (Table 2 ).

To combat the COVID-19 pandemic, researchers around the globe are racing to come up with effective countermeasures. Promising progress has been made in developing vaccines and antiviral drugs. Antivirals are necessary complements of vaccines and are needed for postinfection treatment. Among the viral proteins under investigation as drug targets for SARS-CoV-2, the viral protein RdRp is the most extensively studied, which was followed by the viral protein M pro . 39 Antiviral drug discovery targeting M pro started with the initial efforts of developing of inhibitors against rhinovirus 3C protease (3C pro ). Rhinovirus 3C pro , enterovirus 3C pro , human norovirus 3CL protease, and coronavirus 3CL protease (M pro ) all share the same substrate preference for glutamate at the P1 position, suggesting that 3C pro or 3CL pro inhibitors are promising drug candidates for broad-spectrum antivirals. Over the past few decades, significant progress has been made in designing 3C pro or 3CL pro inhibitors. Rupintrivir (AG7088) and AG7404 are prominent examples of human rhinovirus 3C pro inhibitors that have been evaluated in clinical trials for the treatment of rhinovirus infection. For coronaviruses, GC376 is one of the most advanced lead compounds. It showed broad-spectrum in vitro antiviral activity against SARS-CoV and MERS-CoV, and in vivo antiviral activity in cats infected with feline infectious peritonitis virus. 40,41 Given the sequence similarity between SARS-CoV-2 and SARS-CoV M pro , it became apparent that existing 3C pro or 3CL pro inhibitors might be active against SARS-CoV-2 M pro . Indeed, the most potent M pro inhibitors reported so far such as N3, 13a, 13b, and GC376 all contain the pyrrolidone substitution in the P1 position, with variations in the reactive warhead and P2, P3, and P4 substitutions. 11−14 Interestingly, six compounds, ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin, that share no structural similarity with GC376 were claimed as novel SARS-CoV-2 M pro inhibitors. 12 MS/MS analysis revealed that ebselen, PX-12, and carmofur were able to covalently modify the catalytic cysteine C145 of SARS-CoV-2 M pro . 12 In line with these documented polypharmacology of ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin, we are interested in validating these compounds against SARS-CoV-2 M pro inhibition. Our enzymatic assay results showed that the inhibition of SARS-CoV-2 M pro by these six compounds is dependent on the reducing reagent DTT. In the absence of DTT, all six compounds ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin showed potent inhibition against not only M pro but also two related viral proteases the EV-A71 and EV-D68 3C pro , as well as three unrelated viral proteases, SARS-CoV-2 PL pro and EV-A71 and EV-D68 2A pro (Figure 2 ). However, upon addition of 4 mM DTT, the broad-spectrum enzymatic inhibition of these compounds was largely diminished, except for carmofur and tideglusib, which had weak inhibition against M pro and PL pro with IC 50 values of 28.2 and 30.4 μM, respectively ( Figure 2 ; Table 1 ). In line with the enzymatic assay results, thermal shift binding assay and native MS assay showed that ebselen, disulfiram, PX-12, tideglusib, and shikonin did not bind to M pro in the presence of DTT, while carmofur could still bind to M pro with the addition of DTT (Figures 4 and 5) . These results suggest that the inhibition of M pro by carmofur has certain specificity, although the potency is relatively weak. In contrast, the inhibitory effect and binding of control compounds GC376 against M pro and EV-A71 and EV-D68 3C pro , GRL0618 against PL pro , and telaprevir against EV-A71 and EV-D68 2A pro were not affected by the addition of DTT (Figures. 2, 4, and 5 ). MD simulations provided additional evidence showing that the drug-bound M pro complex is more stable for specific inhibitor GC376 than for promiscuous compounds ebselen and carmofur ( Figure 6 ). Furthermore, it is generally assumed that for specific inhibitors, the enzymatic inhibition potency IC 50 value could be used to predict the cellular antiviral activity. However, despite their apparent inhibition of the EV-A71 and EV-D68 2A pro and 3C pro in the absence of DTT (Table 1) , none of the six compounds, ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin, showed cellular antiviral activity against EV-A71 and EV-D68 in the CPE assay (Table 2 ). Therefore, caution should be taken when interpreting the enzymatic assay inhibition IC 50 values of cysteine proteases obtained in the absence of reducing reagent such as DTT or GSH. In the absence of DTT, the apparent inhibition might be due to either alkylation or oxidation of the cysteine residue by reactive compounds. To rule out such a nonspecific effect, reducing reagents such as DTT, β-ME, or GSH should be added to the enzymatic buffer. Specific cysteine protease inhibitors should not show significant IC 50 shift upon addition of reducing reagent in both the enzymatic assay and the binding assay. Moreover, counter-screening against unrelated cysteine proteases should also be performed as a secondary assay to confirm the specificity. Although these promiscuous compounds such as ebselen have been frequently highlighted as promising drug candidates, 42, 43 the scientific community should be cautious in interpreting the pharmacology of these compounds and be aware of their nonspecific effects.

Cell Lines and Viruses. Human rhabdomyosarcoma (RD) cells were maintained in Dulbecco's modified Eagle's medium (DMEM), supplemented with 10% fetal bovine serum (FBS) and 1% penicillin−streptomycin antibiotics. Cells were kept at 37°C in a 5% CO 2 atmosphere. EV-D68 strain US/MO/14− 18947 (ATCC NR-49129) was purchased from ATCC and amplified in RD cells prior to infection assays. EV-A71 strain 5865/SIN/000009 was obtained from Dr. Chan at the Department of Medical Microbiology, Faculty of Medicine, University of Malaya. 44 Protein Expression and Purification. SARS-CoV-2 M pro . SARS-CoV-2 M pro gene from strain BetaCoV/Wuhan/ WIV04/2019 in the pET29a(+) vector with E. coli codon optimization was ordered from GenScript (Piscataway, NJ). The expression and purification of SARS-CoV-2 M pro was described as previously. 11 SARS-CoV-2 PL pro . SARS-CoV-2 papain-like protease (PL pro ) gene (ORF 1ab 1564−1876) from strain BetaCoV/ Wuhan/WIV04/2019 with E. coli codon optimization was ordered from GenScript in the pET28b(+) vector. The pET28b(+) plasmid with SARS-CoV-2 PL pro gene was transformed into BL21(DE3) cells with kanamycin selection. A single colony was picked to inoculate 10 mL of LB media and was cultured 37°C overnight. This 10 mL culture was added to 1 L of LB media and grown to around OD 600 of 0.8. This culture was cooled on ice for 15 min, then induced with 0.5 mM IPTG. Induced cultures were incubated at 18°C for an additional 24 h and then harvested and lysed the same way as SARS-CoV-2 M pro protein. 11 The supernatant was incubated with Ni-NTA resin for overnight at 4°C on a rotator. The Ni-NTA resin was thoroughly washed with 30 mM imidazole in wash buffer (50 mM Tris [pH 7.5], 150 mM NaCl, 2 mM DTT), and PL pro protein was eluted from Ni-NTA with 300 mM imidazole. Eluted PL pro was dialyzed against 100-fold volume dialysis buffer (50 mM Tris [pH 7.5], 150 mM NaCl, 2 mM DTT) in a 10 000 kDa molecular weight cutoff dialysis tubing.

EV-A71 2A pro . The EV-A71 2A pro gene from strain EV-A71/ 7D3 (genbank accession number MF973167) with E. coli codon optimization was ordered from GenScript in the pET28b(+) vector. The expression and purification of EV-A71 2A pro is same as that for SARS-CoV-2 PL pro described in the above section.

EV-A71 3C pro . The EV-A71 3C pro gene from strain EV-A71/ 7D3 (genbank accession number MF973167) with E. coli codon optimization was ordered from GenScript in the pET28b(+) vector. The expression and purification of EV-A71 3C pro is same as that for SARS-CoV-2 PL pro described in the above section.

EV-D68 2A pro . The EV-D68 2A pro gene from strain US/KY/ 14−18953 with E. coli codon optimization was ordered from GenScript in the pET28b(+) vector. The expression and purification of EV-D68 2A pro was described as previously. 34 EV-D68 3C pro . The EV-AD68 3C pro gene from strain US/ KY/14−18953 with E. coli codon optimization was ordered from GenScript in the pET28b(+) vector. The expression and purification of EV-D68 3C pro is same as that for SARS-CoV-2 PL pro described in the above section.

FRET Substrate Peptide Synthesis. The FRET-based peptide substrates used for the enzymatic assay are shown below:

• SARS-CoV-2 M pro substrate: Dabcyl-KTSAVLQ/ SGFRKME-Edans The synthesis of SARS-CoV-2 M pro , PL pro , EV-A71 2A pro , EV-D68 2A pro ,and EV-D68 3C pro substrates were described previously. 11, 34 Enzymatic Assays. The IC 50 values of the testing compounds against various SARS-CoV-2, EV-A71, and EV-D68 proteases in the presence or in the absence of 4 mM DTT were measured with a common protocol as the following: First, 100 μL of protease (SARS-CoV-2 M pro at 100 nM; SARS-CoV-2 PL pro at 200 nM; EV-A71 2A pro at 3 μM; EV-A71 3C pro at 2 μM; EV-D68 2A pro at 1 μM; or EV-D68 3C pro at 100 nM) was incubated with various concentrations of testing inhibitors at 30°C for 30 min in its reaction buffer in a 96-well plate, and then the reaction was initiated by adding FRET substrate (SARS-CoV-2 M pro and PL pro substrates at 10 μM; EV-A71 and EV-D68 substrates at 20 μM). The reaction was monitored for 2 h, and the initial velocity was calculated using the data from the first 15 min by linear regression. The IC 50 was calculated by plotting the initial velocity against various concentrations of testing inhibitor by using a four parameters dose−response curve in Prism (v8.0) software. The reaction buffers used were as follows:

• SARS-CoV-2 M pro reaction buffer: 20 mM HEPES, pH 6.5, 120 mM NaCl, 0.4 mM EDTA, and 20% glycerol • SARS-CoV-2 PL pro reaction buffer: 50 mM HEPES, pH7.5, 0.01% triton X-100

• EV-A71 2A pro reaction buffer: 50 mM Tris pH 7.0, 150 mM NaCl, 10% glyceol • EV-A71 3C pro reaction buffer: 50 mM Tris pH 7.0, 150 mM NaCl, 1 mM EDTA, 10% glycerol • EV-D68 2A pro reaction buffer: same as EV-A71 2A pro reaction buffer • EV-D68 3C pro reaction buffer: same as EV-A71 3C pro reaction buffer Thermal Shift Binding Assay (TSA). The thermal shift binding assay (TSA) was carried out using a Thermal Fisher QuantStudio 5 Real-Time PCR System as described previously. 34, 36 Briefly, 3 μM protease in its enzymatic reaction buffer (see the "Enzymatic Assays" section for the reaction buffer components) in the presence of 4 mM DTT or in the absence of DTT was incubated with testing compounds at 30°C for 30 min in a 96-well PCR plate. SYPRO orange dye (1×) was added, and the fluorescence of the well was monitored under a temperature gradient range from 20 to 90°C with 0.05°C /s incremental step. The melting temperature (T m ) was calculated as the mid log of the transition phase from the native to the denatured protein using a Boltzmann model (Protein Thermal Shift Software v1.3).

Native Mass Spectrometry. The native MS binding assay of SARS-CoV-2 M pro was carried out using previously described methods. 11 Briefly, purified SARS-CoV-2 M pro was buffer exchanged into 0.2 M ammonium acetate (pH 6.8) at a protein concentration of 6 μM. Each of the ligands tested (GC376, ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12) was diluted to 200 and 100 μM in ethanol. The compounds were then titrated into the protein sample to give a final drug concentration of 10, 20, or 40 μM. For the ligand binding studies containing dithiothreitol (DTT), a 40 mM stock of DTT was dissolved in water. A final concentration of 4 mM DTT was added to each of those samples. For the ligand binding studies without DTT added, an equal volume of nanopore water was added to the samples in place of DTT. The final concentration of protein in each of the samples was 4.9 μM. Each sample contained 4.5 μL of protein, 0.5 μL of ligand, and 0.5 μL of DTT or water. The samples were mixed and incubated at room temperature for 30 min prior to analysis.

Native MS was performed as previously described using Q-Exactive HF quadrupole-Orbitrap mass spectrometer with the Ultra-High Mass Range (UHMR) research modifications (Thermo Fisher Scientific). All of the samples were ionized in positive ion mode using 0.9 kV capillary voltage with the temperature set to 200°C. The resolution of the instrument was set to 15 000 for all samples except for samples containing the compound Jun8−38−3, for which the resolution was set to 30 000. The trapping gas pressure within the instrument was set to 3.50 V of source fragmentation was applied for each of the samples to aid in desolvation of the sample. All samples were analyzed between a 500−15 000 m/z range. All of the data were deconvolved and analyzed using UniDec. 45 Cytopathic Effect Assay (CPE). The EC 50 and CC 50 values for the protease inhibitors investigated in this study were measured using RD cells as described previously. 36 Briefly, RD cells were seeded and grown overnight to ∼90% confluence in a 96-well plate at 37°C and 5% CO 2 . For EV-D68 virus infection, cells were washed with PBS saline and infected with virus diluted in DMEM medium with 2% FBS and 30 mM MgCl 2 . Viruses were incubated with cells for 1 h at 33°C followed by addition of various concentrations of testing protease inhibitors in DMEM medium with 30 mM MgCl 2 . For EV-A71 virus infection, the procedures are identical to those for EV-D68 virus, except that 30 mM MgCl 2 was omitted in all the media, and viruses were infected and incubated at 37°C instead of 33°C. Three days after infection, cells were stained with 66 μg/mL of neutral red dye for 2 h, and neutral red uptake was measured at an absorbance at 540 nM. CC 50 was measured similarly but in the absence of viral infection.

The proximal origin of SARS-CoV-2

A Snapshot of the Global Race for Vaccines Targeting SARS-CoV-2 and the COVID-19 Pandemic

Research and Development on Therapeutic Agents and Vaccines for COVID-19 and Related Human Coronavirus Diseases

Cell entry mechanisms of SARS-CoV-2

SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor

Cathepsin Lselective inhibitors: A potentially promising treatment for COVID-19 patients

Broadspectrum antiviral GS-5734 inhibits both epidemic and zoonotic coronaviruses

An orally bioavailable broad-spectrum antiviral inhibits SARS-CoV-2 in human airway epithelial cell cultures and multiple coronaviruses in mice

Substrate specificity profiling of SARS-CoV-2 M pro protease provides basis for anti-COVID-19 drug design. bioRxiv (Biochemistry)

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

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

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved alpha-ketoamide inhibitors

) Structure and inhibition of the SARS-CoV-2 main protease reveals strategy for developing dual inhibitors against M pro and cathepsin L. bioRxiv (Biophysics)

Targeting the Mycobacterium tuberculosis transpeptidase LdtMt2 with cysteine-reactive inhibitors including ebselen

Ebselen induces reactive oxygen species (ROS)-mediated cytotoxicity in Saccharomyces cerevisiae with inhibition of glutamate dehydrogenase being a target

A small-molecule antivirulence agent for treating Clostridium difficile infection

Ebselen inhibits hepatitis C virus NS3 helicase binding to nucleic acid and prevents viral replication

Interaction of ebselen with glutathione S-transferase and papain in vitro

Identification of ebselen as a potent inhibitor of insulin degrading enzyme by a drug repurposing screening

2020) Ebselen derivatives are very potent dual inhibitors of SARS-CoV-2 proteases -PL pro and M pro in in vitro studies. bioRxiv (Biochemistry)

Disulfiram is a direct and potent inhibitor of human O 6 -methylguanine-DNA methyltransferase (MGMT) in brain tumor cells and mouse brain and markedly increases the alkylating DNA damage

Inhibition of Urease by Disulfiram, an FDA-Approved Thiol Reagent Used in Humans

Structural basis for inactivation of Giardia lamblia carbamate kinase by disulfiram

Disulfiram can inhibit MERS and SARS coronavirus papain-like proteases via different modes

Molecular Mechanism of Inhibition of Acid Ceramidase by Carmofur

-methylpropyl)-dithio]-1H-imidazole inhibits tubulin polymerization through cysteine oxidation

Evidence for Irreversible Inhibition of Glycogen Synthase Kinase-3β by Tideglusib

Structural basis for the inhibition of SARS-CoV-2 main protease by antineoplastic drug carmofur

A noncovalent class of papain-like protease/deubiquitinase inhibitors blocks SARS virus replication

Validating Enterovirus D68−2A pro as an Antiviral Drug Target and the Discovery of Telaprevir as a Potent D68−2A pro Inhibitor

Characterization and Noncovalent Inhibition of the Deubiquitinase and deISGylase Activity of SARS-CoV-2 Papain-Like Protease

A Novel Capsid Binding Inhibitor Displays Potent Antiviral Activity against Enterovirus D68

Pharmacological Characterization of the Mechanism of Action of R523062, a Promising Antiviral for Enterovirus D68

Development of broad-spectrum enterovirus antivirals based on quinoline scaffold

Potential Therapeutic Agents and Associated Bioassay Data for COVID-19 and Related Human Coronavirus Infections

Reversal of the Progression of Fatal Coronavirus Infection in Cats by a Broad-Spectrum Coronavirus Protease Inhibitor

Potential therapeutic use of ebselen for COVID-19 and other respiratory viral infections

Molecular characterization of ebselen binding activity to SARS-CoV-2 main protease

Enterovirus A71 DNA-Launched Infectious Clone as a Robust Reverse Genetic Tool

Bayesian deconvolution of mass and ion mobility spectra: from binary interactions to polydisperse ensembles

The authors declare no competing financial interest.