4.1. Compound library preparation
All solvents were distilled before use. Commercially available reagents were used without further purification. Selenium powder (100 mesh) (Sigma-Aldrich, Saint Louis, MO, USA) used for Na2Se2 and Li2Se2 preparation had a purity of ≥ 99.5%. Freshly distilled MeCN was redistilled twice over P2O5 before preparation of ebselen and its derivatives. MeOH was distilled slowly over a mixture of LAH and CaH2 before hydrogenation of ebselen. CH2Cl2 (DCM) was distilled over P2O5 before preparation of ebselen analogues 25–34. Triethyl amine (Et3N) (POCh, Gliwice, Poland), distilled over NaOH, was stored over NaOH pellets. Anhydrous sodium carbonate (Na2CO3) (POCh, Gliwice, Poland) was ground in a mortar before use. The intermediates, 2-(chloroseleno)benzoyl chloride and bis[(2-chlorocarbony)phenyl] diselenide were prepared from anthranilic acid and elemental selenium via the formation of 2,2’-dicarboxydiphenyl diselenide – a key intermediate, according to the literature procedure [57, 58, 60, 61]. Preparative column chromatography was performed on Merck Si60 silica gel (63–200 µm). Analytical TLC was performed on PET foils precoated with silica gel (Merck silica gel, 60 F254) (Sigma-Aldrich, Saint Louis, MO, USA), and were visualized with light (λmax = 254 nm), or by staining with iodine steam. Melting points were determined on an Electrothermal IA 91100 digital melting-point apparatus using the standard open capillary method. IR spectra (4000–400 cm− 1) were recorded in KBr plates on a Perkin-Elmer 2000 FT-IR spectrometer or on a Fourier transform, Bruker VERTEX 70V spectrometer using diamond ATR accessory. Absorption maxima are reported in wavenumbers (cm− 1). 1H-NMR and 13C-NMR spectra (300.1, 399.8, 600.6 MHz and 75.48, 100.5, 151.0 MHz, respectively) were recorded on a Bruker DRX 300 (Bruker, Rheinstetten, Germany), Jeol 400yh (Jeol, Tokyo, Japan) and Bruker Avance II 600 (Bruker, Poznań, Poland) instruments. NMR spectra recorded in CHCl3-d1 and DMSO-d6 were referenced to the respective residual 1H or 13C signals of the solvents, and chemical shifts (δ) are given in parts per million (ppm), and coupling constants (J) are in Hz. 19F-NMR and 77Se-NMR (376.2 and 76.24 MHz, respectively) were collected on Jeol 400yh instrument. High-resolution mass spectra were collected using electrospray ionization on a Waters LCT Premier XE TOF instrument.
The literature procedure was adapted for the preparation of dilithium diselenide [59], diaryl diselenides 24 [71], 28 [57, 60], 31 [57], 32–33 [60], 36 [57, 61], 38 [57], and 2-(chloroseleno)benzoyl chloride (37) [57]. Purity and homogeneity of known compounds were confirmed by measuring their m.p. for ebselen [57], 1–2 [55], 3 [54], 7–8 [57], 9 [72], 11 [73], 16 [57], 19–22 [60], 23 [57], 24 [74], 26 [74], 28 [60], 29 [74], 31 [57], 32–33 [60], 36 [73], 37 [75], and 38 [73], or FT-IR spectra for ebselen [54], 3 [54], 10 [63], 11 [73], 16 [57], 31 [57], 1H- and/or 13C-NMR spectra for 1–2 [55], 7 [56, 57], 8 [60], 10 [73], 11 [73], 16 [57], 19–22 [60], 23 [57], 24 [71], 28 [60], 31 [57], 32 [60], and 77Se-NMR spectrum for 10 [63], and HRMS for 16 [57], and comparing them with literature data. All new 13, 30, 34, uncharacterized, 4 [76], 5 [74], 6 [77], 9 [72], 12 [78], 14 [78], 15 [79], 17–18 [80], 25 [79], 27 [79], and spectroscopically uncharacterized 26, 29 [74] selenium species were fully characterized. The hydrogen and carbon atom positions in the 1H-NMR and 13C-NMR spectra were supported by the dept-135 or COSY experiments and by 2D-NMR map analysis of the Heteronuclear Multiple-Quantum Correlation (HMQC) and Heteronuclear Multiple Bond Correlation (HMBC), Nuclear Overhauser Enhancement Spectroscopy (NOESY) if measured. See Supporting Information for detailed synthesis protocols and compounds spectroscopic characterization.
4.2. SARS-CoV-2 PLpro preparation
SARS-CoV-2 PLpro was prepared as described [26]. In brief, pGEX6P-1-SARS-CoV-2PLpro was transformed into BL21 (DE3) codon-plus E. coli cells and induced with 0.1 mM IPTG and 0.1 mM ZnSO4 at 18°C overnight. GST-fusion SARS-CoV-2 PLpro was purified using a standard protocol. The fusion protein was cleaved using GST-PreScission protease at 4°C overnight followed with desalting and passing through fresh glutathione beads to remove cleaved GST and GST-PreScission protease. The sample was further purified using Superdex 200 pg size-exclusion columns (GE) equilibrated with 20 mM Tris-Cl pH 8.0, 40 mM NaCl and 2 mM DTT. The peak fractions were pooled and concentrated to ~ 10 mg/ml and snap frozen in liquid nitrogen for later use.
4.3. SARS-CoV-2 Mpro preparation
SARS-CoV-2 Mpro was recombinantly produced as described16. Briefly, the gene of the Mpro was cloned into the PGEX-6p-1 vector, which has a Nsp4-Nsp5 and a PreScission cleavage site at the N- and C-termini, respectively, to generate the authentic target protein. The gene of the target protein was expressed in the E. coli of the BL21-Gold (DE3) (Novagen) strain. The recombinantly produced Mpro was purified by employing HisTrap FF (GE Healthcare) and ion-exchange chromatography (Q FF, GE Healthcare). Finally, the high-purity target protein was subjected to a buffer exchange (20 mM Tris, 150 mM NaCl, 1 mM EDTA, 1 mM DTT, pH 7.8) for further experiments.
4.4. SARS-CoV-2 nsp14 preparation
SARS-CoV-2 mRNA cap guanine N7-methyltransferase nsp14 was prepared as described previously [41]. Briefly, the nsp14 gene was cloned into the pET28 SUMO expression vector. Nsp14 was overexpressed in BL21 (DE3) RIL E. coli (Invitrogen), as a fusion protein with His-tagged SUMO. The fusion protein was purified using HisTrap FFTM column (Cytiva), followed by loading on HiTrap 26/10 Desalting column (Cytiva). To remove the N-terminal tag (6xHis-Sumo) Sumo protease (MCLAB) was added and then the nsp14 protein again purified on HisTrap FFTM column. Flow-through fractions containing nsp14 were collected and separated from N-terminal tag (6xHis-Sumo), and His-tagged Sumo protease. The flow-through fraction was further finally purified on a Superdex 75 pg HiLoad 26/600 gel filtration column (Cytiva). Fractions containing nsp14 were concentrated to 30 µM, flash frozen and stored at − 80°C in a buffer containing 50 mM HEPES (pH 8.0), 100 mM NaCl, 1 mM DTT, 10% glycerol.
4.5. Inhibitor screening
Evaluation of the compound library for inhibitors of SARS-CoV-2 PLpro and SARS-CoV-2 Mpro was carried out in Corning 96-wells plates. For PLpro, 1 µL of each compound in DMSO solution was added to the wells. Next, 79 µL of enzyme preincubated for 10 min at 37°C in assay buffer (50 mM Tris, 5 mM NaCl, 0.075% BSA, pH 7.5) was added to each well. The enzyme was incubated with the compounds at 37°C for 30 min. Next, 20 µL Ac-LRGG-ACC substrate in assay buffer was added to the wells. Final concentrations were: 100 nM enzyme, 10 µM substrate and 1 µM tested compounds. In the assay for Mpro, 1 µL of each compound in DMSO solution was added to the wells. Next, 79 µL of enzyme in assay buffer (50 mM Tris, 1 mM EDTA, pH 7.3) [81] was added to each well and the plate was incubated at room temperature for 2 min. Next, 20 µL of QS1 substrate in assay buffer was added to the wells. Final concentrations were: 100 nM enzyme, 50 µM substrate, and 100 nM or 1 µM tested compounds. Measurements were carried out at 37°C using a Molecular Devices Spectramax Gemini XPS spectrofluorometer. ACC fluorophore release was monitored for 30 min (λex = 355 nm, λem = 460 nm). For the further analysis, the linear range of the progress curves was used. Measurements were performed at least in duplicate. Results were presented as mean values of relative enzyme inhibition (%, compared to the control measurement without inhibitor) with standard deviations. During the assays, the DMSO concentration in the wells was < 2%.
4.6. IC50 determination
To determine IC50, the relative activity of the investigated proteases was assessed in at least 11 different concentrations of selected inhibitors. Initial compound concentrations were found experimentally. Serial dilutions of inhibitors in assay buffers (described above) were prepared in 96-well plates (20 µL of each dilution in wells). For SARS-CoV-2 PLpro, 60 µL enzyme preincubated for 10 min at 37°C in assay buffer was added to the wells. The enzyme was incubated with inhibitors for 30 min at 37°C. Next, 20 µL substrate (Ac-LRGG-ACC) in assay buffer was added to the wells. Final concentrations were 100 nM enzyme and 10 µM substrate. For SARS-CoV-2 Mpro, 60 µL enzyme was added with no preincubation. The enzyme was incubated with inhibitor for 2 min at room temperature. Next, 20 µL of substrate (QS1) in the assay buffer was added to the wells. Final concentrations were 100 nM for the enzyme and 50 µM for the substrate. Measurements were carried out at 37°C using a Molecular Devices Spectramax Gemini XPS spectrofluorometer. ACC fluorophore release was monitored for 30 min (λex = 355 nm, λem = 460 nm). IC50 values were determined with GraphPad Prism software using non-linear regression (dose-response – Inhibition equation) and presented as relative enzyme activity vs. inhibitor concentration. Measurements were performed at least in triplicate. Results are presented as mean values with standard deviations. During the assays, the DMSO concentration in wells was < 2%. See Supplementary Information for IC50 graphs.
4.6. IC50 determination with nsp14 N7-MTase
To determine IC50 parameters of ebselen analogs towards the nsp14 enzyme, we used the previously described Py-FLINT assay designed for N7-MTase activity studies [41, 64]. The Py-FLINT probe (1 µM) was incubated with SAM cosubstrate (20 µM), nsp14 (40 nM), and an inhibitor (half-log dilutions logCinh <-2.5;2>). Point fluorescence measurements (λex = 345 nm, λem = 378 nm) were carried out in 96-well black, non-binding assay plates at 30°C. Initial rates V were calculated by fitting a linear curve to the first 10 points (10 minutes). To the obtained dependences V(Cinh) the following four-parameter dose-response equation was fitted:
$$\frac{V}{{V}_{0}}=A1+\frac{A2-A1}{1+{\left(\frac{{C}_{inh}}{{IC}_{50}}\right)}^{p}} , \left(1\right)$$
where A1 and A2 are the bottom and top asymptotes, respectively; Cinh the inhibitor concentration; p is the Hill coefficient, and V/V0 is the ratio of the initial reaction rate with the inhibitor to that without the inhibitor. For curve fitting and IC50 calculations we used GraphPad Prism software.
4.7. Anti-SARS-CoV-2 and cytotoxicity assays in Vero E6 cells
The anti-SARS-CoV-2 activity was measured by determining the extent to which the compounds inhibited the virus-induced cytopathic effect (CPE) and reduced SARS-CoV-2 RNA in Vero E6 cells (ECACC 85020206). For the CPE-based assay, two-fold fold serial dilutions of compounds were added in triplicate in a 384-well plate with 5,000 Vero E6 cells in DMEM medium with 2% FBS, 100 U of penicillin/ml, and 100 µg of streptomycin/ml (all Merck). After 1 h incubation, SARS-CoV-2 (strain hCoV-19/Czech Republic/NRL_6632_2/2020 was isolated in a biosafety level 3 laboratory from nasopharyngeal swab by inoculating Vero CCL81 cells [ECACC 84113001]) was added at multiplicity of infection 0.05 IU/ml. Following three days incubation at 37°C in 5% CO2, the cell viability was determined by addition of XTT solution (Sigma-Aldrich) for 4 h and the absorbance was measured using EnVision plate reader (Perkin Elmer). Drug concentrations required to reduce viral cytopathic effect by 50% (EC50) were calculated using nonlinear regression from plots of percentage cell viability versus log10 drug concentration using GraphPad Prism v.9.0.0 Software. For RNA reduction-based assay, two-fold fold serial dilutions of compounds were added in triplicate in 96-well plate with 20,000 Vero cells plated day before in the same medium as above. After 1 h incubation, SARS-CoV-2 was added at multiplicity of infection 0.05 IU/cell. After 2 h, virus was removed and new compound was added to the cells. Cells were incubated for two days, then the medium was used as a template in RT-qPCR (Multiplex RT-PCR for COVID-19, Diana Biotechnologies, Czech Republic). Compound concentrations required to reduce SARS-CoV-2 RNA copy number by 50% (EC50) were calculated from plots of percentage of RNA copy number versus log10 drug concentration as above.
Cytotoxicity was evaluated by incubating two-fold serial dilutions of each compound with Vero E6 cells. Following three days incubation at 37°C in 5% CO2, the cell viability was determined by addition of XTT solution as above. The compound concentrations resulting in 50% reduction of absorbance (CC50) were calculated from plots of percentage of absorbance versus log10 drug concentration as above.