Jing Wu1, Bo Feng2,3, Li-Xin Gao1,3, Chun Zhang1, Jia Li2,3,4, Da-Jun Xiang5, Yi Zang3, Wen-Long Wang1. 1. School of Life Sciences and Health Engineering, Jiangnan University, Wuxi 214122, China. 2. School of Life Science and Biopharmaceutics, Shenyang Pharmaceutical University, Shenyang 110016, China. 3. State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China. 4. Zhongshan Institute for Drug Discovery, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Zhongshan Tsuihang New District, Zhongshan 528400, China. 5. Xishan People's Hospital of Wuxi City, Wuxi 214105, China.
Abstract
The COVID-19 pandemic caused by SARS-CoV-2 is a global burden on human health and economy. The 3-Chymotrypsin-like cysteine protease (3CLpro) becomes an attractive target for SARS-CoV-2 due to its important role in viral replication. We synthesized a series of 8H-indeno[1,2-d]thiazole derivatives and evaluated their biochemical activities against SARS-CoV-2 3CLpro. Among them, the representative compound 7a displayed inhibitory activity with an IC50 of 1.28 ± 0.17 μM against SARS-CoV-2 3CLpro. Molecular docking of 7a against 3CLpro was performed and the binding mode was rationalized. These preliminary results provide a unique prototype for the development of novel inhibitors against SARS-CoV-2 3CLpro.
The COVID-19 pandemic caused by SARS-CoV-2 is a global burden on human health and economy. The 3-Chymotrypsin-like cysteine protease (3CLpro) becomes an attractive target for SARS-CoV-2 due to its important role in viral replication. We synthesized a series of 8H-indeno[1,2-d]thiazole derivatives and evaluated their biochemical activities against SARS-CoV-2 3CLpro. Among them, the representative compound 7a displayed inhibitory activity with an IC50 of 1.28 ± 0.17 μM against SARS-CoV-2 3CLpro. Molecular docking of 7a against 3CLpro was performed and the binding mode was rationalized. These preliminary results provide a unique prototype for the development of novel inhibitors against SARS-CoV-2 3CLpro.
Entities:
Keywords:
COVID-19; Mpro inhibitors; drug design and synthesis; structure-activity relationships (SAR)
The global pandemic of coronavirus disease (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posted major challenges to public health systems and the economy worldwide [1,2,3,4,5]. There have been 434 million confirmed cases of COVID-19 worldwide as of the end of February 2022, and almost 6 million deaths have been reported [6]. Although multiple effective vaccines against COVID-19 are available, reinfections and breakthrough infections are frequently reported [7,8]. In addition, the virus is continuing to evolve, and a new variant named Omicron enables the virus to evade the immune protective barrier due to a large number of mutations in the receptor binding sites [9,10,11]. Therefore, it is urgent to develop effective drugs and specific treatments for people who are infected by COVID-19 with severe symptoms.3CLpro (also called Mpro) plays an essential role during replication and transcription of SARS-CoV-2 and has been regarded as an attractive target for treating COVID-19 and other coronavirus-caused diseases [12,13,14]. The development of 3CLpro inhibitors has attracted much attention from medicinal chemists and the pharmaceutical industry. The collective efforts culminated in the recent approval of Paxlovid (nirmatrelvir) by FDA for the treatment of SARS-CoV-2 [15]. As shown in Figure 1, Most known 3CLpro inhibitors are peptidomimetic inhibitors containing a warhead of Michael acceptor, such as nirmatrelvir with nitrile [16], YH-53 with benzothiazolyl ketone [17], compound 1 with α-ketoamide [18], and compound 2 with aldehyde [19]. Others are nonpeptidic inhibitors including covalent and noncovalent inhibitors. Covalent inhibitors, such as Carmofur, Shikonin [20], and 3 [21], are identified by high-throughput screening. Noncovalent inhibitor CCF0058981 [22] and flavonoid analogs (baicalin, baicalein, and 4′-O-Methylscutellarein) [23,24] were obtained through structure-based optimization and from traditional Chinese medicines, respectively.
Figure 1
SARS-CoV-2 3CLpro inhibitors.
In pursuit of novel 3CLpro inhibitors, we identified 8H-indeno[1,2-d]thiazole derivative 4 as a novel SARS-CoV-2 3CLpro inhibitor (IC50 = 6.42 ± 0.90 μM) through high-throughput screening of our compound collection (Figure 2). This result provided us with an opportunity to explore novel small molecule inhibitors against SARS-CoV-2 3CLpro. Herein, we designed and synthesized a series of 8H-indeno[1,2-d]thiazole derivatives, evaluated their inhibitory activities against SARS-CoV-2 3CLpro, and elucidated the SARs. Selected compound 7a was subjected to molecular docking to predict the binding mode with SARS-CoV-2 3CLpro.
Figure 2
Structure of 8H-indeno[1,2-d]thiazole derivatives.
2. Results and Discussion
2.1. Design and Synthesis of 8H-Indeno[1,2-d]thiazole Derivatives
Based on the structure of compound 4, 14 new 8H-indeno[1,2-d]thiazole derivatives (compounds 7a–7l, and 10a–10b) (shown in Scheme 1 and Scheme 2) were designed and synthesized through a two-step synthesis from the appropriate ketone and thiourea [25,26,27,28]. Adjusting the methoxy group of compound 4 from position 5 to position 6 afforded compound 7a. Considering the effects of steric hindrance and electron withdrawing, compounds 7b–7e were synthesized by substitution of the methoxy group for the butoxy, isobutoxy, and methyl groups and for the chlorine atom. After replacing the 3,5-dimethoxybenzamido moiety in compound 7a with 3,4,5-trimethoxybenzamido, 3,5-diacetoxybenzamido, 3-methoxybenzamido, 3-fluorobenzamido, thiophene-2-carboxamido, and 4-chlorobenzamido, compounds 7f–7k were obtained. To evaluate the effect of ring expansion, compound 7l was synthesized. Finally, ring opening analogues 10a and 10b were synthesized to elucidate the effect of the central ring on the inhibition of 3CLpro.
Scheme 1
(a) thiourea, bromine, ethanol, 100 °C, 5–6 h; (b) aromatic acid, HATU, DIPEA, DMF, r t, 2–3 h, 25–50%.
Scheme 2
(a) thiourea, iodine, 110 °C, 10 h; (b) aromatic acid, HATU, DIPEA, DMF, r t, 2–3 h, 35–40%.
2.2. SARS-CoV-2 3CLpro Inhibitory Activities and Structure-Activity Relationships
All synthesized compounds were evaluated for inhibitory activity against SARS-CoV-2 3CLpro using PF-07321332 as positive control [29,30,31], and the results were detailed in Table 1. We initially prepared 7a from the commercially available compound 5a by the route outlined in Scheme 1. We noticed that compound 7a with 6-methoxy group on the phenyl ring exhibited inhibitory activity against SARS-CoV-2 3CLpro with 1.28 ± 0.17 μM, about five times more potent than compound 4 with 5-methoxy group on the phenyl ring. The result indicated that the position of the methoxy group on the phenyl ring significantly affected inhibitory activities against SARS-CoV-2 3CLpro. To explore the SAR of this seemingly important position, methoxy group on compound 7a was replaced by butoxy (7b), isobutoxy (7c), methyl groups (7d), and chlorine atom (7e); the inhibitory activities of the corresponding compounds 7b–7e were completely abolished. These results demonstrated that the effect of steric hindrance at this position was detrimental to inhibitory activities. The SAR of R was explored next. Replacement of the 3,5-dimethoxybenzamido moiety with 3,4,5-trimethoxybenzamido moiety, 3,5-diacetoxybenzamido moiety, 3-methoxybenzamido moiety, 3-fluorobenzamido moiety, thiophene-2-carboxamido moiety, and 4-chlorobenzamido moiety led to compounds 7f, 7g, 7h, 7i, 7j, and 7k, respectively. The inhibitory activity of compounds 7f and 7g dropped significantly, while compound 7h almost maintained its inhibitory activities. These results indicated that the extra steric hindrance had negative impact on the inhibitory activities. Compared to compound 7h, the inhibitory activities of compounds 7i–7k diminished; these results indicated that introduction of an electron-withdrawing group or heterocyclic ring on the scaffold of 8H-indeno[1,2-d]thiazole took negative roles for inhibitory activities. Expanding the five-membered ring on compound 7a to a six-membered ring led to compound 7l, which unfortunately did not show any inhibitory activity against SARS-CoV-2 3CLpro. Opening the five-membered ring on compound 7a resulted in compounds 10a and 10b, which also lost inhibitory activities. These results indicated that the five-membered ring on compound 7a is important for the inhibitory activity against SARS-CoV-2 3CLpro.
Table 1
Inhibitory activities of target compounds against SARS-CoV-2 3CLpro.
Compd.
R1
R3
n
SARS-CoV-2 3CLpro
Inhibition (%) at 20 μM
IC50 (μM)
7a
methoxy
1
89.5 ± 2.0
1.28 ± 0.17
7b
butoxy
1
0.5 ± 4.9
>20
7c
isobutoxy
1
−3.1 ± 1.7
>20
7d
methyl
1
21.7 ± 2.2
>20
7e
chloro
1
27.2 ± 5.3
>20
7f
methoxy
1
5.0 ± 5.6
>20
7g
methoxy
1
32.6 ± 6.8
>20
7h
methoxy
1
72.5 ± 6.1
2.86 ± 0.11
7i
methoxy
1
20.3 ± 4.7
>20
7j
methoxy
1
31.9 ± 18.2
>20
7k
methoxy
1
1.5 ± 4.5
>20
7l
methoxy
2
−13.1 ± 1.7
>20
10a
-
-
1.9 ± 2.1
>20
10b
-
-
1.8 ± 3.5
>20
PF-07321332(nirmatrelvir)
99.5 ± 0.1
0.012 ± 0.001
2.3. Predicting Binding Mode of with 3CLpro
To explore the interaction mode between small molecule 7a and 3CLpro (PDB code: 6M2N) [23], we carried out molecular docking by applying AutoDock 4.2 program [31,32,33,34]. Figure 3a showed that 7a docked well into the binding pockets S1 and S2 of 3CLpro, in which the S1, S2 sites play a key role in substrate recognition [35]. As illustrated in Figure 3b, the indene moiety of compound 7a buried deeply into the hydrophobic S2 subsite with π-electrons with Arg188 and hydrophobic interaction with Met165; the 3,5-dimethoxybenzamido moiety of compound 7a formed strong H-bonds with Asn142, Glu166 on S1 subsite, while compounds 4 and 7h escaped from S1 subsite, as shown in Supplementary Materials Figures S1 and S2.
Figure 3
(a) surf representation of the compound 7a (bonds representation) in the 3CLpro S1 (red), S2 (blue), S1′ (green), S4 (orange) binding pocket; (b) the docking results of 7a and 3CLpro (PDB code: 6M2N, active residues in 3.0 Å range around 7a).
3. Materials and Methods
3.1. Chemistry
All chemical reagents are reagent grade and used as purchased. 1H NMR (400 MHz) spectra were recorded on a Bruker AVIII 400 MHz spectrometer (Bruker, Billerica, MA, USA). The chemical shifts were reported in parts per million (ppm) using the 2.50 signal of DMSO (1H NMR) and the 39.52 signal of DMSO (13C NMR) as internal standards. ESI Mass spectra (MS) were obtained on a SHIMADZU 2020 Liquid Chromatograph Mass Spectrometer (SHIMADZU, Kyoto, Japan).
3.1.1. General Procedure for the Synthesis of Compounds 7a–7k (Exemplified by 7a)
To a solution of 5f (528.2 mg, 3.0 mmol) in dry ethanol (10 mL) were added thiourea (456.7 mg, 6.0 mmol) and bromine (0.2 mL, 3.3 mmol) at room temperature. The reaction solution was stirred at 100 °C for 5–6 h. At the end of the reaction, the solvent was evaporated and aqueous ammonium hydroxide (25%) was added to the residue. The precipitated solid 6f was collected without purification for the next step. The mixture of 6f (255.2 mg, 1.1 mmol), 3,5-dimethoxybenzoic acid (182.1 mg, 1.0 mmol), HATU (380.2 mg, 1.0 mmol), and DIPEA (0.5 mL 3.0 mmol) in DMF (6 mL) was stirred at room temperature for 2 h. The reaction mixture was quenched with water. The aqueous layer was extracted with EtOAc (20 mL × 2). The combined organic layers were dried over Na2SO4. The residue was purified by column chromatography on silica gel (eluting with DCM) to afford compound 7l (103.0 mg, yield 26%) as a white solid.1H NMR (400 MHz, DMSO-d6) δ 12.66 (s, 1H), 7.66 (dd, J = 8.4, 2.0 Hz, 1H), 7.32 (t, J = 2.0 Hz, 2H), 6.88 (s, 1H), 6.85 (dd, J = 8.4, 2.4 Hz, 1H), 6.73 (d, J = 2.4 Hz, 1H), 3.83 (s, 6H), 3.77 (s, 3H), 3.00–2.91 (m, 4H) ppm. 13C NMR (100 MHz, DMSO-d6) δ 164.31, 160.45, 158.42, 156.44, 143.66, 136.68, 133.90, 124.23, 123.33, 121.55, 114.08, 111.82, 105.75, 105.06, 55.59, 55.09, 28.65, 20.74 ppm. MS (ESI): m/z calcd for C21H21N2O4S [M + H]+ 397.12, found 396.95.
3.1.3. General Procedure of Synthesis of 10a–10b (Exemplified by 10a)
A mixture of 8a (10.0 mmol, 1.0 equiv), thiourea (20.0 mmol, 2.0 equiv), and iodine (10.0 mmol, 1.0 equiv) was stirred at 110 °C for 10 h. After the reaction was completed, the residue was triturated with MTBE and adjusted to pH 9–10 with 25% ammonium hydroxide. The precipitated solid was collected and washed with EtOAc (30 mL × 2) and NaHCO3 (15 mL × 2) aqueous solution. The separated organic layer dried over Na2SO4 and evaporated to dryness to afford crude product 9a. The mixture of 9a (3.3 mmol, 1.1 equiv), aromatic acid (3.0 mmol, 1.0 equiv), HATU (3.0 mmol, 1.0 equiv), and DIPEA (9.0 mmol, 3.0 equiv) in DMF (20 mL) was stirred at room temperature for 2 h. Then the reaction mixture was quenched with water. The aqueous layer was extracted with EtOAc (30 mL × 2). The combined organic layers were dried over Na2SO4. The residue was purified by column chromatography on silica gel (eluting with DCM) to afford compound 10a as a white solid (406.7 mg, yield 35%). 1H NMR (400 MHz, DMSO-d6) δ 12.67 (s, 1H), 7.57 (d, J = 8.4 Hz, 1H), 7.32 (d, J = 2.0 Hz, 2H), 7.21 (s, 1H), 6.86 (d, J = 2.4 Hz, 1H), 6.83 (dd, J =8.4, 2.8 Hz, 1H), 6.74 (t, J = 2.4 Hz, 1H), 3.83 (s, 6H), 3.77 (s, 3H), 2.43 (s, 3H) ppm. 13C NMR (100 MHz, DMSO-d6) δ 164.74, 160.45, 158.60, 157.76, 149.01, 136.97, 134.22, 130.74, 127.31, 115.96, 111.22, 110.14, 105.80, 104.90, 55.58, 55.04, 21.26 ppm. MS (ESI): m/z calcd for C20H21N2O4S [M + H]+ 385.12, found 385.20.3,5-dimethoxy-N-(4-(4-methoxy-3-methylphenyl)thiazol-2-yl)benzamide (10b), eluting with DCM, yield = 40%; 1H NMR (400 MHz, DMSO-d6) δ 12.70 (s, 1H), 7.77 (d, J = 2.4 Hz, 1H), 7.75 (s, 1H), 7.49 (s, 1H), 7.33 (d, J = 2.4 Hz, 2H), 6.99 (d, J = 8.8 Hz, 1H), 6.74 (t, J = 2.4 Hz, 1H), 3.84 (s, 6H), 3.82 (s, 3H), 2.20 (s, 3H) ppm. 13C NMR (100 MHz, DMSO-d6) δ 164.59, 160.46, 158.26, 157.15, 149.27, 133.86, 128.07, 126.72, 125.68, 124.62, 110.36, 106.38, 105.79, 105.09, 55.59, 55.31, 16.21 ppm. MS (ESI): m/z calcd for C20H21N2O4S [M + H]+ 385.12, found 385.25.
3.2. Molecule Docking
The protease structure, SARS-CoV-2 3CLpro enzyme (PDB code: 6M2N) with 2.2 Å, was obtained from the the Protein Data Bank at the RCSB site (http://www.rcsb.org (accessed on 6 March 2022)). The molecule docking used the Lamarckian genetic algorithm local search method and the semiempirical free energy calculation method in the AutoDock 4.2 program. Additionally, the charge was added by Kollman in AutoDock 4.2, The docking methold was employed on rigid receptor conformation, all the rotatable torsional bonds of compound 7a were set free, the size of grid box was set at to 10.4 nm × 12.6 nm × 11.0 nm points with a 0.0375 nm spacing and grid center (−33.798 −46.566 39.065), and the other parameters were maintained at their default settings.
3.3. Enzymatic Activity and Inhibition Assays
The enzyme activity and inhibition assays of SARS-CoV-2 3CLpro have been described previously [20,36]. Briefly, the recombinant SARS-CoV-2 3CLpro (40 nM at a final concentration) was mixed with each compound in 50 μL of assay buffer (20 mM Tris, pH 7.3, 150 mM NaCl, 1% Glycerol, 0.01% Tween-20) and incubated for 10 min. The reaction was initiated by adding the fluorogenic substrate MCA-AVLQSGFRK (DNP) K (GL Biochem, Shanghai, China), with a final concentration of 40 μM. After that, the fluorescence signal at 320 nm (excitation)/405 nm (emission) was immediately measured by continuous 10 points for 5 min with an EnVision multimode plate reader (Perkin Elmer, Waltham, MA, USA). The initial velocity was measured when the protease reaction was proceeding in a linear fashion; plots of fluorescence units versus time were fitted with linear regression to determine initial velocity. Plots of initial velocity versus inhibitor concentration were fitted using a four-parameter concentration–response model in GraphPad Prism 8 to calculate the IC50 values. All data are shown as mean ± SD, n = 3 biological replicates.
4. Conclusions
In summary, we synthesized a series of 8H-Indeno[1,2-d]thiazole derivatives and evaluated their biochemical activities against SARS-CoV-2 3CLpro. Among them, the representative compound 7a displayed inhibitory activity with an IC50 of 1.28 ± 0.17 μM against SARS-CoV-2 3CLpro. Molecular docking elucidated that 7a was well-docked into the binding pockets S1 and S2 of 3CLpro. These preliminary results could provide a possible opportunity for the development of novel inhibitors against SARS-CoV-2 3CLpro with optimal potency and improved pharmacological properties.
Authors: Vicky Mody; Shashidharamurthy Taval; Joanna Ho; Savannah Wills; Ahmed Mawri; Latasha Lawson; Maximilian C C J C Ebert; Guillaume M Fortin; Srujana Rayalam Journal: Commun Biol Date: 2021-01-20