Next Article in Journal
Determination of N-Acetyl-l-cysteine Ethyl Ester (NACET) by Flow Injection Analysis and Spectrophotometric Detection Using Different Thiol-Sensitive Ligands
Previous Article in Journal
Electron-Deficient Acetylenes as Three-Modal Adjuvants in SNH Reaction of Pyridinoids with Phosphorus Nucleophiles
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 26
Issue 22
10.3390/molecules26226825
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Antiviral Activities of Halogenated Emodin Derivatives against Human Coronavirus NL63

by
Monika Horvat
1,
Martina Avbelj
2,
María Beatriz Durán-Alonso
3,
Mihailo Banjanac
4,
Hrvoje Petković
2,* and
Jernej Iskra
1,*
1
Faculty of Chemistry and Chemical Technology, University of Ljubljana, Večna pot 113, 1000 Ljubljana, Slovenia
2
Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, 1000 Ljubljana, Slovenia
3
Unit of Excellence, Institute of Biology and Molecular Genetics (IBGM), University of Valladolid-CSIC, 47003 Valladolid, Spain
4
Fidelta d.o.o., Prilaz baruna Filipovića 29, 10000 Zagreb, Croatia
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(22), 6825; https://doi.org/10.3390/molecules26226825
Submission received: 25 October 2021 / Revised: 6 November 2021 / Accepted: 8 November 2021 / Published: 11 November 2021
Browse Figures
Versions Notes

Abstract

:
The current COVID-19 outbreak has highlighted the need for the development of new vaccines and drugs to combat Severe Acute Respiratory Syndrome-Coronavirus-2 (SARS-CoV-2). Recently, various drugs have been proposed as potentially effective against COVID-19, such as remdesivir, infliximab and imatinib. Natural plants have been used as an alternative source of drugs for thousands of years, and some of them are effective for the treatment of various viral diseases. Emodin (1,3,8-trihydroxy-6-methylanthracene-9,10-dione) is a biologically active anthraquinone with antiviral activity that is found in various plants. We studied the selectivity of electrophilic aromatic substitution reactions on an emodin core (halogenation, nitration and sulfonation), which resulted in a library of emodin derivatives. The main aim of this work was to carry out an initial evaluation of the potential to improve the activity of emodin against human coronavirus NL63 (HCoV-NL63) and also to generate a set of initial SAR guidelines. We have prepared emodin derivatives which displayed significant anti-HCoV-NL63 activity. We observed that halogenation of emodin can improve its antiviral activity. The most active compound in this study was the iodinated emodin analogue E_3I, whose anti-HCoV-NL63 activity was comparable to that of remdesivir. Evaluation of the emodin analogues also revealed some unwanted toxicity to Vero cells. Since new synthetic routes are now available that allow modification of the emodin structure, it is reasonable to expect that analogues with significantly improved anti-HCoV-NL63 activity and lowered toxicity may thus be generated.

1. Introduction

Severe Acute Respiratory Syndrome-Coronavirus-2 (SARS-CoV-2) is a family of enveloped positive-sense RNA viruses that cause life-threatening respiratory infections and severe pneumonia in humans [1,2]. Coronavirus (CoV) entry into host cells (pulmonary and parabronchial epithelial cells) is mediated by spike protein, which is responsible for binding to receptors ACE-2 and mediating virus–host membrane fusion [3,4,5]. The development of effective antiviral drugs with a broad spectrum of activity has been hampered by viral diversity and the ability of SARS-CoV to mutate rapidly, even during an epidemic. It is therefore very important to develop antiviral drugs that effectively and safely inhibit the spread of SARS-CoV, or at least significantly alleviate the symptoms of SARS-CoV infection. In particular, the development of simple, small compounds that can be produced and administered inexpensively would be of great importance. Recently, several potential repurposed drugs against COVID-19 (SARS-CoV-2 virus) have been found, such as remdesivir, infliximab and imatinib. Remdesivir has potent antiviral activity and has already been approved for urgent use [6,7]. However, to curb the spread of infection, it is important to identify new drug-leads that are more broadly effective against CoV.
Nature is not only a source of emerging mutant viruses, but also a reservoir of natural products that play a crucial role in drug development. Emodin (1,3,8-trihydroxy-6-methylanthracene-9,10-dione), a potent natural bioactive anthraquinone, is found in various plants, lichens and molds, such as Cassia obtusifolia and Cassia occidentalis, Rhamnus orbiculatus, Aloe vera, Japanese knotweed, Polygonum multiflorum, Rheum palmatum, Scutellaria baicalensis and Rumex chalepensis [8,9,10,11,12,13]. Emodin is known for its anti-oxidant, anti-ulcerogenic, anti-bacterial, anti-fibrotic, anti-inflammatory, anti-cardiovascular, anti-viral and anti-cancer activities [14,15,16,17,18,19]. It has demonstrated antitumor activity against various cancers such as leukemia, squamous cell carcinoma of human tongue, lung cancer, gallbladder cancer, breast cancer, colon cancer and others [9,14,15,20,21]. Modified emodin compounds have therefore shown relevant pharmacological activity [22,23,24,25,26,27]. Of particular interest and promise are the results obtained with halogenated derivatives of emodin. In 2014, Huang and co-workers found that halogenated emodin derivatives can exert a potent inhibitory activity on bacterial topoisomerase I and DNA gyrase. The best results were obtained with 2,4-diiodoemodin [28]. In 2017, the research group led by Sukhatme and Sun reported the structure–activity relationship (SAR) of emodin and emodin derivatives as ATP citrate lyase (ACL) inhibitors. Halogenated emodin analogues (2-iodoemodin, 2-chloroemodin, 4-chloroemodin and 2,4-dibromoemodin) showed significantly increased activity [29]. Later, Tansakul’s group demonstrated that the hydroxyl and methyl groups were crucial for anti-MRSA (anti-methicillin-resistant Staphylococcus aureus) activity. All compounds containing two halogenated atoms (I, Br or Cl) at positions 2 and 4 were active against MRSA. The best results were obtained in the presence of an iodine atom (2,4-diiodoemodin) [30]. 4-chloroemodin was found to significantly inhibit the growth of gram-positive bacteria, especially that of common drug-resistant MRSA and VRE (vancomycin-resistant enterococci) isolates, through a dual antibacterial mechanism that interacts with the bacterial cell membrane and DNA [31].
In addition, emerging evidence suggests that emodin displays broad spectrum antiviral activities against herpes simplex viruses (HSV-1 and HSV-2) [32,33], hepatitis B virus (HBV) [34,35], Japanese encephalitis virus (JEV) [13], Human cytomegalovirus [36], Influenza A [37], Zika virus [38], Coxsackie B virus [39,40], Poliovirus [41], Cypridine herpesvirus 3 (CyHV-3) [42] and in a number of viral diseases. Through its antiviral activity, emodin can also prevent or reduce SARS-CoV infection [43,44,45,46,47,48]. Hsiang and co-workers reported that emodin can block the interaction of SARS-CoV spike protein with ACE-2 and infectivity of spike protein-pseudotyped retrovirus on Vero E6 cells [49]. In 2011, Schwarz and co-workers showed that emodin can inhibit the 3a ion channel of coronavirus as well as the release of SARS-CoV from infected cells [19]. Singha Roy and Das performed a blind molecular docking analysis of natural anthraquinones against SARS-CoV-2 main protease Mpro. The results suggest that natural emodins may prove to be effective inhibitors of COVID-19 by binding to the catalytic dyad, HIS41 and CYS145, through non-covalent forces near the active site [50].
The aim of the current study was to synthesize a series of emodin analogues and investigate their activities against human coronavirus NL63 (HCoV-NL63) [51]. For this work, we focused on the HCoV-NL63 virus, which causes mild to moderate upper respiratory tract infections in children, severe lower respiratory tract infection, croup and bronchiolitis [52]. HCoV-NL63 and SARS-CoV-2 both belong to the group of common human coronaviruses and they both use the ACE2 receptor to infect cells. HCoV-NL63 is thus a low pathogenic common coronavirus that may be used to study SARS under BSL2 conditions. This is the reason why HCoV-NL63 has been suggested as a suitable surrogate virus for studying SARS-CoV-2 [53].
Herein, we report the strategies followed for the selective introduction of NO2, SO3H and halogen atoms into the anthraquinone ring of emodin. Taking into account the potency of halogenated emodins, our work focused on the selective decoration of emodin, containing different halogen atoms and substitution patterns (Figure 1) and the evaluation of their antiviral activity against HCoV-NL63. Our results demonstrate that the presence of different functional groups in the emodin scaffold has a significant impact on their anti-HCoV-NL63 activity.

2. Results and Discussion

2.1. Synthesis of Emodin Derivatives

The natural product emodin E_H provides an entry point for the introduction of various functional groups on aromatic rings by electrophilic aromatic substitution. Our aim was to create a library of emodin derivatives by modifying the aromatic ring through halogenation, nitration, amination and sulfonation.
First, we investigated the selective halogenation of emodin E_H by classical reagents (N-chlorosuccinimide—NCS, N-bromosuccinimide—NBS, N-iodoosuccinimide—NIS) and by an alternative method—oxidative halogenation with hydrogen peroxide as oxidant. 2-Iodoemodin E_I was prepared according to the published methods [28] using I2 and NaHCO3 in 73% yield (Table 1, entry 1) and in the presence of NIS reagent in 84% yield (Table 1, entry 6). Oxidative iodination with I2 and 30% hydrogen peroxide in 2-MeTHF also allowed the selective synthesis of the same product with a better yield of 93% (Table 1, entry 16). For the oxidative iodination, 2 equivalents of iodine and 2.5 equivalents of hydrogen peroxide were used. While conducting the experimental work, it was observed that the reaction does not proceed in the absence of the oxidant H2O2 (entry 11) and that a higher amount of iodine is required for the quantitative conversion to E_I (entries 12–15). We also found that the amount of hydrogen peroxide has a minimal effect on the increase in conversion. A larger number of iodine substituents could not be introduced into the aromatic core of emodin using oxidative iodination. To introduce additional iodine atoms, I2/NaHCO3 or NIS had to be used. Nevertheless, the diiodinated emodin E_2I could not be selectively prepared, regardless of the amount of iodine or NIS used. 2,4,7-Triiodoemodin E_3I was selectively formed by applying the classical method (I2, NaHCO3) [28] in 81% yield or using NIS reagent in 79% yield as a brown-orange powder (Table 1, entries 3 and 10). We also attempted to synthesize the fully iodinated emodin E_4I, but this compound could not be prepared regardless of the reaction conditions. We tried increasing the temperature and adding the activator H2SO4 to the reagent NIS. In all cases, the triiodinated product appeared selectively. The positions of iodine on the aromatic rings in compounds E_I and E_3I were confirmed by 2D NMR spectroscopy (Figures S6–S13, Supplementary Material).
Next, bromination was carried out using NBS and by oxidative halogenation using H2O2/HBr. Neither of these methods could be effective for selective synthesis of monobromo-derivative E_Br, as shown by the results in Table 2. Although a reaction with 1 equivalent of NBS was carried out in an ice bath, both E_Br and E_2Br were formed simultaneously (Table 2, entries 1 and 2). The same results were observed for oxidative bromination with HBr and H2O2. Regardless of the amounts of HBr and hydrogen peroxide used, a mixture of the two products was always formed. Using 2D NMR spectroscopy, we found that in the case of E_Br, bromine binds to site 2 in emodin, either when the NBS reagent or HBr/H2O2 is used.
E_2Br was selectively prepared according to a published procedure [29] with a slight modification consisting of treating emodin E_H with NBS in THF at 0 °C (88% yield). Due to the high reactivity of emodin, the bromination required a relatively short reaction time (30 min) and a low temperature (0 °C). The same product was obtained by the oxidative halogenation method (HBr, H2O2) in 2,2,2 trifluoroethanol (TFE) in 91% yield (Table 2, entry 14). Unfortunately, the oxidative method, although more environmentally friendly, did not allow the introduction of more bromine atoms. Despite the higher amount of HBr and the use of the activating solvent TFE, the reaction stopped at the dibrominated product. The use of an NBS reagent allowed a greater number of bromine atoms to be introduced, but we encountered problems with the selectivity in the preparation of E_3Br. Regardless of the temperature at which the reaction was carried out and the amount of reagent used, E_2Br and/or E_4Br were also formed (Table 2, entries 4–6). Selectively, we prepared the orange-colored product 2,4,6,8-tetrabromo-1,3,5-trihydroxy-7-methylanthracene-9,10-dione E_4Br (Table 2, entry 7).
Chlorination of emodin gave similar results to bromination. The monosubstituted product could not be prepared selectively with the NCS reagent or oxidatively with HCl/H2O2, since the dichloro derivative E_2Cl was also formed. It was observed that the regioselectivity of chlorination to E_Cl depends on the method used—HCl/H2O2 or NCS—as determined by conducting 1D and 2D NMR spectroscopy on the crude reaction mixture consisting of E_H, E Cl or/and E_2Cl. The results showed that when NCS was used, the Cl was bound at position 2 (E_Cl-2), whereas when HCl/H2O2 was used, it was bound at position 4 (E_Cl-4) (Figures S1–S3). Good yields of the disubstituted emodin 2,4-dichloroemodin E_2Cl were obtained selectively with NCS [30] using H2SO4 as catalyst in a 3-h process at room temperature (Table 3, entry 3). When emodin was treated with 6 equivalents of NCS in the presence of H2SO4 at a reflux temperature, no formation of three- or four-chloroemodin occurred. Only two chlorine atoms could be introduced into emodin. Oxidative chlorination was proved to be an unsuitable procedure for the selective synthesis of chlorine products. Despite the combination of different ratios of HCl and H2O2, emodin was always converted into a mixture of different products that could not be separated (Table 3, entries 5–9).
Next, we investigated the selective nitration of E_H. We attempted to selectively prepare nitroemodin with one, two or three NO2 functional groups, but none of the conditions tested led to successful results. The results showed that despite the use of small amounts of nitric and sulfuric acids, a mixture of various nitrated products E_X formed that could not be separated (Table 4, entries 1–7). It was concluded that the nitration reaction is not selective towards a particular product despite the absence of H2SO4. We prepared the tetranitroemodin E_4NO2 using 6.2 equiv. HNO3 and 10 equiv. H2SO4 in 91% yield (Table 4, entry 8), reducing excess reagents by an order of magnitude according to a published method [24] (Table 4, entry 9).
In addition, we also prepared amino-substituted emodin E_NH2 by the method described in an earlier publication (Figure 2) [24].
Sulfonation was carried out with oleum at room temperature. After 24 h, the reaction mixture was analyzed by NMR; three different products were observed. According to the NMR spectra obtained from the crude reaction mixture, trisulfonated, disulfonated and monosulfonated emodin were formed. However, the composition of the mixture changed with time and reverse sulfonation occurred. Hence, following the work-up procedure, the crude reaction mixture was stirred in MeCN/hexane and after 3 h the product was completely converted to a stable monosubstituted E_SO3H, which was isolated as a brown solid in 76% yield (Figure 3).
1,3,8-Trimethoxy-6-methyl-9,10-anthraquinone E_OCH3 was obtained in good yield (81%) by applying a previously described reaction procedure (Me2SO4 and base K2CO3) [30]. Brominated trimethoxyemodin E_Br_OCH3 was synthesized from tetrabrominated emodin E_4Br by the same method. After purification, the product E_4Br_OCH3 was obtained in 87% yield.

2.2. Antiviral Activity of Emodin Analogues

The primary objective of this study was to determine the potential of a library of 11 emodin analogues (Figure 4) to combat the cytopathic effects of HCoV-NL63. The non-transformed primate kidney Vero cell line was used for this work. An additional series of studies was performed to evaluate the effects of each compound on the viability of Vero cells.

2.2.1. Evaluation of the Viability of Vero Cells in the Presence of Emodin and Emodin Analogues

Published reports indicate that emodin displays a broad spectrum of activities, including cytotoxicity [14]. Therefore, in addition to evaluating the activity of emodin and the emodin analogues against HCoV-NL63, it was important to test the potential cytotoxicity of these compounds. Viability assays were performed using Vero cells; IC50 curves for each compound and other additional data are shown in Figure 5 and in Table S1 (raw data in Table S4).
Interestingly, only E_I showed higher toxicity to Vero cells compared to emodin (E_H). E_4Br, E_Cl and E_OMe were less toxic than emodin, while virtually no toxicity was observed with E_SO3H, E_NH2 and E_Br_OMe. However, it is important to consider that although all emodin analogues were soluble in DMSO at 50 mM concentrations, some analogues exhibited a tendency to precipitate when added to the cell culture medium; compounds containing methoxy groups were the most difficult to dissolve in DMSO and were therefore probably the least soluble in aqueous solutions. It is therefore possible that these compounds gave lower cytotoxicity indexes due to their poor solubility in the cell culture medium.

2.2.2. Anti-Viral Activity of Emodin and Emodin Derivatives

In another series of studies, we evaluated the potential of emodin and emodin derivatives to protect Vero cells against the cytopathic effects induced by HCoV-NL63. Dead cell protease activity was used as a measure of impaired cell viability, as described in the Methods section; chloroquine and remdesivir were used as controls. The results obtained are presented in Figure 6 and Table 5. Additional data can be found in Supplementary Materials (Table S2; raw data in Table S5). Three compounds, E_OMe, E_SO3H and E_Br_OMe, did not exhibit significant anti-viral activity; as described above, this could be due to their lower solubility in aqueous solutions. Emodin and the emodin analogues E_4NO2 and E_I impaired Vero cell viability and had anti-viral effects at very similar concentrations. Five other compounds, E_3I, E_4Br, E_2Br, E_2Cl and E_NH2, showed anti-viral activity at concentrations lower than those at which they impaired Vero cell viability; among these, E_NH2 exhibited the least toxicity to Vero cells. On the other hand, E_3I was the compound that demonstrated the strongest anti-HCoV-NL63 activity, along with the largest difference (about 10-fold) between the IC50 value for anti-HCoV-NL63 activity and that for cytotoxicity, thus providing the largest therapeutic window. Nevertheless, the IC50 value for cytotoxicity of E_3I was still considerably high under the applied experimental conditions, namely identical to that of E_H (Table S1).
As mentioned earlier, chloroquine and remdesivir were used as controls, since these drugs are being considered as potentially effective against COVID-19. Interestingly, the anti-HCoV-NL63 activity of some of the emodin analogues was much higher than that of chloroquine. On the other hand, the emodin analogues E_3I and E_2Br displayed anti-HCoV-NL63 activity comparable to that of remdesivir; all the data obtained on the activities of chloroquine and remdesivir are presented in Supplementary Materials (Table S3 and Figure S5; raw data in Tables S6 and S7).

3. Materials and Methods

Emodin was purchased from Fluorochem Ltd. (Glossop, UK). All other reagents and solvents were of reagent-grade quality and were obtained from commercial suppliers Honeywell (Seelze, Germany) and Sigma-Aldrich (Taufkirchen, Germany). TLC was performed on Merck-60-F254 plates (Merck, Darmstadt, Germany) using mixtures of EtOAc:EtOH (10:1), CH2Cl2:EtOH (100:1) or EtOAc:MeOH (20:1). Crude emodin preparations were purified by column chromatography on silica gel (63–200 µm, 70–230 mesh ASTM; Honeywell, Seelze, Germany). The isolated compounds were characterized by 1H, 13C NMR spectra, HRMS and IR analysis. 1H and 13C NMR spectra were recorded on Bruker Avance III 500 instruments (Bruker, Billerica, MA, USA). IR spectra were recorded on Bruker Alpha II FTIR Instrument (Bruker, Billerica, MA, USA). HR-MS were recorded on LC MS system Agilent 6224 Accurate Mass TOF LC/MS (Agilent Technologies, Santa Clara, CA, USA).

3.1. Compound Synthesis and Structure Confirmation

1,3,8-trihydroxy-2-iodo-6-methylanthracene-9,10-dione (E_I) [28]. Iodine (254 mg, 1.0 mmol) and 30% H2O2 (255 µL, 2.5 mmol) were added to a stirred solution of emodin (135 mg, 0.5 mmol) in 2-MeTHF (5 mL). The reaction mixture was stirred at room temperature for 24 h. The reaction was monitored by TLC (CH2Cl2:EtOH = 100:1). After the reaction was complete, the reaction mixture was washed with NaHSO3 and extracted with dichloromethane (3 × 30 mL). The organic layer was washed with water (1 × 30 mL), dried over anhydrous Na2SO4 and evaporated under a vacuum. The crude reaction product was washed with hexane (5 mL) and acetonitrile (5 mL) to remove soluble impurities. The product was dried in vacuum to provide the product (184 mg, 93%) as an orange solid. 1H NMR (500 MHz, DMSO-d6, 25 °C): δ = 13.07 (s, 1H, OH), 12.26 (s, 1H, OH), 11.83 (s, 1H, OH), 7.47 (s, 1H, ArH), 7.21 (s, 1H, ArH), 7.16 (s, 1H, ArH), 2.41 (s, 3H, CH3) ppm. 13C NMR (126 MHz, DMSO-d6, 25 °C): δ = 189.4, 181.2, 165.2, 163.5, 161.4, 148.6, 134.2, 132.7, 124.3, 120.6, 113.1, 108.5, 106.8, 83.2, 21.6 ppm. IR: 3339, 1665, 1615, 1473, 1379, 1263, 1168, 949 cm−1. HRMS (ESI−): m/z calcd for C15H9IO5 394.9422 [M-H]−, found: 394.9428 [M-H]−.
1,3,8-trihydroxy-2,4,7-triiodo-6-methylanthracene-9,10-dione (E_3I) [28]. Iodine (1.3 g, 5.0 mmol) was added at a temperature of 0 °C to a stirred solution of emodin (135 mg, 0.5 mmol) in THF (13 mL) and water (13 mL). NaHCO3 (3.5 g, 42.0 mmol) was then added in a stepwise manner. The reaction mixture was stirred for 24 h at room temperature. The reaction was monitored by TLC (CH2Cl2:EtOH = 100:1). After completion of the reaction, the mixture was extracted with dichloromethane (3 × 30 mL). The organic layer was washed with water (1 × 30 mL), dried over anhydrous Na2SO4 and evaporated under vacuum. The crude reaction product was washed with hexane (5 mL) and acetonitrile (3 × 5 mL) to remove soluble impurities. The product was dried in vacuum to provide the product (262 mg, 81%) as an orange–brown solid. 1H NMR (500 MHz, DMSO-d6, 25 °C): δ = 13.61 (s, 1H, OH), 12.60 (s, 1H, OH), 7.47 (s, 1H, ArH), 2.48 (s, 3H, CH3) ppm. 13C NMR (126 MHz, DMSO-d6, 25 °C): δ = 187.2, 180.2, 165.1, 163.8, 159.5, 151.8, 133.3, 132.0, 120.4, 111.8, 110.0, 101.3, 84.1, 82.5, 29.3 ppm. IR: 3359, 1614, 1377, 1236, 1111, 1043 cm−1. HRMS (ESI−): m/z calcd for C15H7I3O5 646.7355 [M-H]−, found: 646.7380 [M-H]−.
2,4-dibromo-1,3,8-trihydroxy-6-methylanthracene-9,10-dione (E_2Br) [28]. HBr (48%, 2.0 mmol) and 30% H2O2 (255 µL, 2.5 mmol) were added to a stirred solution of emodin (135 mg, 0.5 mmol) in TFE (5 mL). The reaction mixture was stirred at room temperature for 24 h. The reaction was monitored by TLC (CH2Cl2:EtOH = 100:1). After completion of the reaction, the mixture was washed with NaHSO3 and extracted with CH2Cl2 (3 × 30 mL). The organic layer was washed with water (1 × 30 mL), dried over anhydrous Na2SO4 and evaporated under a vacuum. The crude reaction product was washed with hexane (5 mL) and acetonitrile (5 mL) to remove soluble impurities. The product was dried in a vacuum to provide the product (194.7 mg, 91%) as an orange solid. 1H NMR (500 MHz, DMSO-d6, 25 °C): δ = 13. 89 (s, 1H, OH), 12.10 (s, 1H, OH), 7.42 (d, J = 0.6 Hz, 1H, ArH), 7.09 (d, J = 0.6 Hz, 1H, ArH), 2.40 (s, 3H, CH3) ppm. 13C NMR (126 MHz, DMSO-d6, 25 °C): δ = 188.8, 181.2, 161.1, 160.7, 149.2, 134.0, 130.5, 123.9, 121.1, 113.2, 110.6, 106.6, 105.7, 67.5, 22.1 ppm. IR: 3346, 1664, 1376, 1289, 1250, 1213, 1114 cm−1. HRMS (ESI−): m/z calcd for C15H8Br2O5 424.8666 [M-H]−, found: 424.8667 [M-H]−.
2,4,5,7-tetrabromo-1,3,8-trihydroxy-6-methylanthracene-9,10-dione (E_4Br) [54]. N-Bromosuccinimide (890 mg, 5.0 mmol) was added to a solution of emodin (270 mg, 1.0 mmol) in THF (5 mL) and stirred for 24 h at room temperature. After completion of the reaction, the mixture was extracted with ethyl acetate (3 × 30 mL). The organic layer was dried over anhydrous Na2SO4 and the solvent was evaporated under a vacuum. The crude reaction product was washed with hexane (5 mL) and acetonitrile (3 × 5 mL) to remove soluble impurities. The product was dried in a vacuum to provide the product (486 mg, 83%) as an orange solid. 1H NMR (500 MHz, DMSO-d6, 25 °C): δ = 2.72 (s, 3H, CH3) ppm. 13C NMR (126 MHz, DMSO-d6, 25 °C): δ = 185.4, 183.0, 162.5, 159.5, 156.3, 147.3, 133.5, 132.5, 119.7, 115.6, 113.4, 107.9, 106.8, 103.8, 25.6 ppm. IR: 3302, 1679, 1621, 1566, 1379, 1223, 1161, 1129, 1049, 805, 760. IR: 3404, 1666, 1617, 1368, 1315, 1231, 1205, 1123, 1046 cm−1. HRMS (ESI−): m/z calcd for C15H6Br4O5 580.6876 [M-H]−, found: 580.6886 [M-H]−.
2,4-dichloro-1,3,8-trihydroxy-6-methylanthracene-9,10-dione (E_2Cl) [30]. N-Chlorosuccinimide (333.8 mg, 2.5 mmol) and conc. H2SO4 (0.5 mL) were added to a solution of emodin (270 mg, 1.0 mmol) in THF (5 mL) and stirred for 3 h at room temperature. After the reaction was complete, the mixture was extracted with ethyl acetate (3 × 30 mL). The organic layer was dried over anhydrous Na2SO4 and evaporated under vacuum. The crude reaction product was washed with hexane (5 mL) and acetonitrile (3 × 5 mL) to remove soluble impurities. The product was dried in vacuum to provide the product (274.7 mg, 81%) as a yellow solid. 1H NMR (500 MHz, DMSO-d6, 25 °C): δ = 13.57 (s, 1H, OH), 11.70 (s, 1H, OH), 7.47 (s, 1H, ArH), 7.16 (s, 1H, ArH), 2.42 (s, 3H, CH3) ppm. 13C NMR (126 MHz, DMSO-d6, 25 °C): δ = 185.0, 182.0, 165.6, 160.7, 160.2, 147.2, 133.9, 127.5, 123.4, 122.3, 120.3, 113.5, 112.6, 104.2, 22.0 ppm. IR: 3312, 1622, 1536, 1379, 1199, 1162, 1100 cm−1. HRMS (ESI−): m/z calcd for C15H8Cl2O5 336.9676 [M-H]−, found: 336.9686 [M-H]−.
1,3,8-trihydroxy-6-methyl-2,4,5,7-tetranitroanthracene-9,10-dione (E_4NO2) [24]. HNO3 (64%, 0.5 mL) was added dropwise over 10 min onto a mixture of emodin (270 mg, 1.0 mmol) and concentrated H2SO4 (96%, 5 mL) at 0 °C. The reaction mixture was stirred at 0 °C for 30 min. The mixture was poured onto ice and extracted with ethyl acetate (3 × 30 mL). The organic layer was dried over anhydrous Na2SO4 and the solvent was evaporated under vacuum. The crude reaction was purified by column chromatography using an ethyl acetate/ethanol (10/1) mobile phase. The solvent was evaporated in vacuum to provide the product (408 mg, 91%) as an orange–red solid. 1H NMR (500 MHz, DMSO-d6, 25 °C): δ = 2.21 (s, 3H, CH3) ppm. 13C NMR (126 MHz, DMSO-d6, 25 °C): δ = 183.8, 177.3, 161.0, 156.6, 151.9, 143.2, 140.4, 139.9, 134.8, 128.8, 124.3, 122.4, 117.0, 96.7, 12.3 ppm. IR: 3351, 1638, 1538, 1369, 1170 cm−1. HRMS (ESI+): m/z calcd for C15H6N4O13 451.0004 [M+H]+, found: 451.0002 [M+H]+.
4-((3-aminopropyl)amino)-1,3,8-trihydroxy-6-methylanthracene-9,10-dione (E_NH2) [24]. (Diacetoxyiodo)benzene (354 mg, 1.1 mmol) was added to a solution of emodin (270 mg, 1.0 mmol) in 1,3-diaminopropane (40 mL). The reaction mixture was stirred for 24 h at room temperature. After the reaction was complete, the mixture was poured into cold water (200 mL), 10 M HCl (60 mL), neutralized with saturated NaHCO3 (200 mL), and extracted with ethyl acetate (3 × 75 mL). The organic layer was dried over anhydrous Na2SO4 and evaporated under vacuum. The crude reaction product was purified by column chromatography using ethyl acetate/methanol (20/1) mobile phase. The solvent was evaporated in vacuum to provide the product (243 mg, 71%) as a violet solid. 1H NMR (500 MHz, DMSO-d6, 25 °C): δ = 14.92 (s, 1H, OH), 13.61 (s, 1H, OH), 12.34 (s, 1H, OH), 7.50 (s, 1H, ArH), 6.81 (s, 1H, ArH), 5.53 (s, 1H, ArH), 4.17 (q, J = 6.4 Hz, 2H, CH2), 2.89 (t, J = 6.9 Hz, 2H, CH2), 2.36 (s, 3H, CH3), 1.93–1.87 (m, 2H, CH2) ppm. 13C NMR (126 MHz, DMSO-d6, 25 °C): δ = 179.2, 176.6, 173.1, 168.4, 160.9, 149.7, 142.7, 135.2, 119.5, 117.5, 116.5, 107.4, 105.6, 100.4, 41.3, 36.9, 29.6, 22.1 ppm. IR: 3195, 1738, 1576, 1502, 1401, 1355 cm−1. HRMS: m/z calcd for C18H18N2O5 343.1288 [M+H]+, found: 343.1290 [M+H]+.
1,3,8-trihydroxy-6-methyl-9,10-dioxo-9,10-dihydroanthracene-5-sulfonic acid (E_SO3H). Oleum (8 mL) was added dropwise onto emodin (270 mg, 1.0 mmol) for over 15 min at room temperature. The reaction mixture was stirred for 24 h at room temperature. The reaction mixture was poured onto ice and extracted with ethyl acetate (3 × 30 mL). The organic layer was dried over anhydrous Na2SO4 and the solvent was evaporated under a vacuum. The crude reaction product was stirred in MeCN and hexane for 3 h and was then purified by column chromatography using ethyl acetate/ethanol (10/1) mobile phase. The solvent was evaporated in vacuum to provide the product (266 mg, 76%) as a brown solid. 1H NMR (500 MHz, DMSO-d6, 25 °C): δ = 11.40 (s, 1H, OH), 11.16 (s, 1H, OH), 7.35 (d, J = 1.9 Hz, 1H, ArH), 7.14 (d, J = 1.9 Hz, 1H, ArH), 6.75 (s, 1H, ArH), 2.38 (s, 3H, CH3) ppm. 13C NMR (126 MHz, DMSO-d6, 25 °C): δ = 158.2, 151.1, 146.8, 139.7, 134.7, 133.5, 126.7, 121.4, 112.1, 110.6, 109.4, 108.9, 105.0, 96.7, 14.6 ppm. IR: 3451, 1738, 1600, 1426, 1367, 1202, 1094 cm−1. HRMS (ESI−): m/z calcd for C15H10O8S 349.0024 [M-H]−, found: 349.0032 [M-H]−.
1,3,8-trimethoxy-6-methylanthracene-9,10-dione (E_OCH3) [30]. Potassium carbonate (415 mg, 3.0 mmol) was added to a solution of emodin (100 mg, 0.37 mmol) in acetone (7 mL). Then dimethyl sulfate (285 µL, 3.0 mmol) was added slowly and the reaction mixture was stirred at reflux for 24 h. The reaction mixture was allowed to cool to room temperature. After cooling to room temperature, the solvent was evaporated. Then water (5 mL) and acetone (5 mL) were added to the reaction mixture under stirring for 15 min. The product was filtrated off, washed with water, and dried in vacuum to provide the product (94 mg, 81%) as a yellow–white solid. 1H NMR (500 MHz, Chloroform-d, 25 °C): δ = 7.64 (s, 1H, ArH), 7.32 (s, 1H, ArH), 7.09 (s, 1H, ArH), 6.76 (s, 1H, ArH), 3.98 (s, 3H, OCH3), 3.96 (s, 3H, OCH3), 3.95 (s, 3H, OCH3), 2.47 (s, 3H, CH3) ppm. 13C NMR (126 MHz, Chloroform-d, 25 °C): δ = 184.4, 181.8, 163.7, 161.7, 159.8, 144.6, 136.4, 134.4, 121.5, 119.6, 119.0, 118.4, 105.3, 101.9, 56.5, 56.5, 55.9, 22.1 ppm. IR: 1657, 1599, 1322, 1241, 1021, 946, 910 cm−1. HRMS (ESI+): m/z calcd for C18H16O5 313.1071 [M+H]+, found: 313.1077 [M+H]+.
2,4,5,7-tetrabromo-1,3,8-trimethoxy-6-methylanthracene-9,10-dione (E_Br_OCH3). Potassium carbonate (415 mg, 3.0 mmol) was added to a solution of brominated emodin 5 (216 mg, 0.37 mmol) in acetone (7 mL). Then dimethyl sulfate (285 µL, 3.0 mmol) was added slowly and the reaction mixture was heated to reflux for 24 h. The reaction mixture was allowed to cool to room temperature. After cooling to room temperature, the solvent was evaporated. Then water (5 mL) and acetone (5 mL) were added to the reaction mixture with stirring for 15 min. The product was filtrated off, washed with water and dried in a vacuum to provide the product (202 mg, 87 %) as a light pink solid. 1H NMR (500 MHz, Chloroform-d, 25 °C): δ = 4.02 (s, 3H, OCH3), 4.01 (s, 3H, OCH3), 3.98 (s, 3H, OCH3), 2.78 (s, 3H, CH3) ppm. 13C NMR (126 MHz, Chloroform-d, 25 °C): δ = 184.1, 179.5, 160.1, 156.9, 155.1, 147.0, 135.8, 135.0, 128.7, 128.0, 126.7, 121.7, 117.7, 112.2, 63.2, 63.1, 61.1, 25.9 ppm. IR: 2153, 2036, 1695, 1367, 1322, 1216, 991 cm−1. HRMS: m/z calcd for C18H12Br4O5 624.7491 [M+H]+, found: 624.7487 [M+H]+.

3.2. Evaluation of Antiviral Activity

Compound preparation. For testing purposes, all compounds were dissolved in DMSO to a final concentration of 50 mM. With these stock solutions, mother plates were prepared in DMSO and generated stocks for testing 8 point–dose responses. 1:1 dilutions were prepared, starting from the original 50 mM stock. Namely, 30 µL of stock solution was added to the first column in the plate. 15 µL of DMSO was added to all the other wells; stocks to test successive doses were produced by transferring 15 µL of compound solution from the preceding column to the next. Each mother plate held the test solutions corresponding to two compounds. Chloroquine and Remdesivir were prepared in three-fold dilutions, chloroquine starting at a 30 µM concentration and Remdesivir at 10 µM.
Cell culturing and viral infection. Vero cells (Cercopithecus aethiops lung epithelial cells, ATCC, CRL-81), were seeded in EMEM medium containing 10% FBS in inner wells of 96-well white plates, at a density of 20000 cells/well. The following day, the plating medium was removed and replaced with 100 µL of fresh EMEM medium containing 2% FBS; a 100 nL solution of each compound was then also added to the wells using a Mosquito pipetting device (TTPlabtech, Cambridge, England). Cells were kept in this medium for 4 days at 33 °C in 5% CO2 before evaluating the cytotoxicity of each compound.
When conducting experiments to evaluate the anti-viral activity of the emodin compounds, the culture medium was changed at 24 h post-plating to 50 µL of fresh EMEM medium containing 2% FBS; 100 nL solution of each test compound was also added to the wells, as described above and an additional 50 µL of the same medium containing a 1:5 dilution of a HCoV-NL63 viral stock was added (Human Coronavirus, Strain NL63, FR-304, IRR (International Reagent Resource). Cultures were maintained in this medium for 4 days at 33 °C, in 5% CO2; the anti-viral activity of each compound was then measured.

3.3. Evaluation of Cytotoxic and Anti-Viral Activities

In order to evaluate any potential cytotoxicity of each of the test compounds, 50 µL/well of Cell Titer-Glo Luminescent reagent was added to each culture following a 4 day-incubation in the presence of these compounds. Luminescence was measured after 5 min of incubation using a Spectra Max i3 instrument (Molecular Devices, San Jose, CA, USA); the obtained values were considered proportional to cellular ATP content. Cytotoxic activity of the tested compounds on Vero cells was expressed as the percentage of live cells compared to that in control cultures that had not been exposed to any compound and had been assigned live cell rates of 100%.
Assays to evaluate anti-viral activity were carried out at 4 days post-infection by adding 50 µL/well of Cyto Tox-Fluor reagent to cultures priorly infected with HCoV-NL63 virus. Plates were incubated for 1 h and dead-cell protease activity was assayed by measuring the fluorescent product that had formed (480/520), using a Spectra Max i3 instrument. Anti-viral activity of each compound was presented as the percentage of live cells compared to that in control Vero cultures that had been infected with the HCoV-NL63 virus but had not been exposed to any of the test compounds; these latter cultures were assigned live cell rates of 100%. All IC50 values were calculated using GraphPad Prism 8 software (GraphPad Software, LLC, San Diego, CA, USA).
Assays were considered valid when the signal to background ratio (signal from cells infected with the virus versus signal from uninfected cells) was higher than 2 and Z’ was higher than 0.35.

4. Conclusions

Emodin is a natural anthraquinone commonly found in plants. Numerous studies (reviewed by Dong and coworkers [14]) show that emodin has a wide spectrum of pharmacological properties. Among these, it is now clear that it also has anti-viral activity. However, emodin can also cause toxicity, such as hepatotoxicity and nephrotoxicity. In accordance with these data, our study demonstrated the toxicity of emodin to Vero cells, a kidney cell line. The main objective of this work was to apply a series of synthetic approaches to the preparation of emodin analogues and to subsequently conduct structure-activity relation (SAR) studies on them, which could shed light on how to improve the anti-HCoV-NL63 activity and to reduce the cytotoxicity of these compounds. Our work, in agreement with the published literature reporting improved activity of halogenated emodin derivatives, showed that halogenation of emodin resulted in enhanced anti-HCoV-NL63 activity; however, cytotoxicity was still relatively high. On the other hand, the introduction of charged groups, such as SO3H and NH2, resulted in a significant decrease in undesirable cytotoxicity, although this also resulted in decreased anti-HCoV-NL63 activity. In addition, the introduction of methoxy functionalities into emodin did not seem to improve its anti-HCoV-NL63 activity; this was probably related to the lower solubility of these compounds in aqueous solutions.
In summary, we applied various synthetic routes to prepare a limited number of emodin analogues that allowed us to perform some initial studies aimed at improving the anti-HCoV-NL63 activity of emodin. Of particular importance was the establishment of SAR guidelines related to anti-HCoV-NL63 activity. The availability of new synthetic routes for the diversification of an emodin structure should allow the generation of new emodin analogues with significantly improved properties, such as increased solubility and, most importantly, stronger anti-HCoV-NL63 activity and reduced toxicity, so that a much wider therapeutic window can realistically be generated.

5. Patents

Patent pending (LU500249).

Supplementary Materials

The following are available online. Figure S1: Chlorination of emodin E_H with NCS and HCl/H2O2, Figure S2: Characterization of E_Cl-2 with 1H NMR and 2D (HSQC, HMBC) NMR in DMSO, 500 MHz, Figure S3: Characterization of E_Cl-4 with 1H and 2D (HSQC, HMBC) NMR in DMSO, 500 MHz, Figure S4: Following a desulfonation of trisubstituted emodin E_3SO3H to monosubstituted E_SO3H in MeCN/hexane at room temperature by 1H NMR (bottom spectra: crude reaction mixture, middle spectra: after 1 h, top spectra: after 3 h), Figure S5: IC50 curves for anti-HCoV-NL63 effects of standard compounds chloroquine and Remdesivir, Figure S6: 1H NMR spectrum of compound E_I in DMSO, 500 MHz, Figure S7: 13C NMR spectrum of compound E_I in DMSO, 500 MHz, Figure S8: 2D HSQC NMR spectrum of compound E_I in DMSO, 500 MHz, Figure S9: 2D HMBC NMR spectrum of compound E_I in DMSO, 500 MHz, Figure S10: 1H NMR spectrum of compound E_3I in DMSO, 500 MHz, Figure S11: 13C NMR spectrum of compound E_3I in DMSO, 500 MHz, Figure S12: 2D HSQC NMR spectrum of compound E_3I in DMSO, 500 MHz, Figure S13: 2D HMBC NMR spectrum of compound E_3I in DMSO, 500 MHz, Figure S14: 1H NMR spectrum of compound E_2Br in DMSO, 500 MHz, Figure S15: 13C NMR spectrum of compound E_2Br in DMSO, 500 MHz, Figure S16: 2D HSQC NMR spectrum of compound E_2Br in DMSO, 500 MHz, Figure S17: 2D HMBC NMR spectrum of compound E_2Br in DMSO, 500 MHz, Figure S18: 1H NMR spectrum of compound E_4Br in DMSO, 500 MHz, Figure S19: 13C NMR spectrum of compound E_4Br in DMSO, 500 MHz, Figure S20: 1H NMR spectrum of compound E_2Cl in DMSO, 500 MHz, Figure S21: 13C NMR spectrum of compound E_2Cl in DMSO, 500 MHz, Figure S22: 2D HSQC NMR spectrum of compound E_2Cl in DMSO, 500 MHz, Figure S23: 2D HMBC NMR spectrum of compound E_2Cl in DMSO, 500 MHz, Figure S24: 1H NMR spectrum of compound E_4NO2 in DMSO, 500 MHz, Figure S25: 13C NMR spectrum of compound E_4NO2 in DMSO, 500 MHz, Figure S26: 1H NMR spectrum of compound E_NH2 in DMSO, 500 MHz, Figure S27: 13C NMR spectrum of compound E_NH2 in DMSO, 500 MHz, Figure S28: 1H NMR spectrum of compound E_SO3H in DMSO, 500 MHz, Figure S29: 13C NMR spectrum of compound E_SO3H in DMSO, 500 MHz, Figure S30: 2D HSQC NMR spectrum of compound E_SO3H in DMSO, 500 MHz, Figure S31: 2D HMBC NMR spectrum of compound E_SO3H in DMSO, 500 MHz, Figure S32: 2D NOESY NMR spectrum of compound E_SO3H in DMSO, 500 MHz, Figure S33: 1H NMR spectrum of compound E_OCH3 in CDCl3, 500 MHz, Figure S34: 13C NMR spectrum of compound E_OCH3 in CDCl3, 500 MHz, Figure S35: 1H NMR spectrum of compound E_Br_OCH3 in CDCl3, 500 MHz, Figure S36: 13C NMR spectrum of compound E_Br_OCH3 in CDCl3, 500 MHz, Table S1: IC50 (µM) values indicating the effect of emodin and emodin analogues on Vero cell viability, Table S2: IC50 (µM) values for anti HCoV-NL63 effects of emodin and emodin analogs, Table S3: IC50 (µM) values for anti HCoV-NL63 effects of chloroquine and Remdesivir, Table S4: Percentage of live cells following exposure to various concentrations of emodin and emodin analogues, Table S5: Percentage of live cells following viral infection and exposure to various concentrations of emodin and emodin analogues, Table S6: Percentage of live cells following viral infection and exposure to various concentration of chloroquine, Table S7: Percentage of live cells following viral infection and exposure to various concentration of Remdesivir.

Author Contributions

Conceptualization: J.I. and H.P.; Data curation: M.H., M.A., M.B.D.-A., M.B., H.P. and J.I.; Investigation: M.H., M.A. and M.B.; Methodology: M.H., M.A., M.B.D.-A., M.B., H.P. and J.I.; Project administration: H.P. and J.I.; Supervision: H.P. and J.I.; Visualization: M.H., M.B.D.-A., H.P. and J.I.; Writing—original draft: M.H., M.B.D.-A., H.P. and J.I.; Writing—review and editing: M.H., M.B.D.-A., M.B., H.P. and J.I.; Funding acquisition: H.P. and J.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Javna Agencija za Raziskovalno Dejanvnost RS, grant number P1-0134 and P4–0116.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data present in this study are available in Supplementary Materials.

Acknowledgments

We thank the Slovenian Research Agency (P1-0134) and the APPLAUSE project for financial support to J.I. The APPLAUSE project is co-financed by the European Regional Development Fund through the Urban Innovative Actions (UIA) initiative. We also thank the Slovenian Research Agency (P4–0116) for financial support to H.P. The authors are grateful to the staff of The Centre for Research Infrastructure at the Faculty of Chemistry and Chemical Technology (IC UL FCCT).

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are available from the authors.

References

  1. Xia, S.; Zhu, Y.; Liu, M.; Lan, Q.; Xu, W.; Wu, Y.; Ying, T.; Liu, S.; Shi, Z.; Jiang, S.; et al. Fusion mechanism of 2019-nCoV and fusion inhibitors targeting HR1 domain in spike protein. Cell. Mol. Immunol. 2020, 17, 765–767. [Google Scholar] [CrossRef]
  2. Xia, S.; Liu, M.; Wang, C.; Xu, W.; Lan, Q.; Feng, S.; Qi, F.; Bao, L.; Du, L.; Liu, S.; et al. Inhibition of SARS-CoV-2 (previously 2019-nCoV) infection by a highly potent pan-coronavirus fusion inhibitor targeting its spike protein that harbors a high capacity to mediate membrane fusion. Cell Res. 2020, 30, 343–355. [Google Scholar] [CrossRef] [Green Version]
  3. Attia, Y.A.; El-Saadony, M.T.; Swelum, A.A.; Qattan, S.Y.A.; Al-qurashi, A.D.; Asiry, K.A.; Shafi, M.E.; Elbestawy, A.R.; Gado, A.R.; Khafaga, A.F.; et al. COVID-19: Pathogenesis, advances in treatment and vaccine development and environmental impact-an updated review. Environ. Sci. Pollut. Res. 2021, 24, 22241–22264. [Google Scholar] [CrossRef] [PubMed]
  4. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.-H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271–280. [Google Scholar] [CrossRef] [PubMed]
  5. Walls, A.C.; Park, Y.-J.; Tortorici, M.A.; Wall, A.; McGuire, A.T.; Veesler, D. Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell 2020, 181, 281–292. [Google Scholar] [CrossRef] [PubMed]
  6. Awadasseid, A.; Wu, Y.; Tanaka, Y.; Zhang, W. Effective drugs used to combat SARS-CoV-2 infection and the current status of vaccines. Biomed. Pharmacother. 2021, 137, 111330. [Google Scholar] [CrossRef] [PubMed]
  7. Pruijssers, A.J.; George, A.S.; Schäfer, A.; Leist, S.R.; Gralinksi, L.E.; Dinnon, K.H.; Yount, B.L.; Agostini, M.L.; Stevens, L.J.; Chappell, J.D.; et al. Remdesivir Inhibits SARS-CoV-2 in Human Lung Cells and Chimeric SARS-CoV Expressing the SARS-CoV-2 RNA Polymerase in Mice. Cell Rep. 2020, 32, 107940. [Google Scholar] [CrossRef]
  8. Izhaki, I. Emodin—A secondary metabolite with multiple ecological functions in higher plants. New Phytol. 2002, 155, 205–217. [Google Scholar] [CrossRef] [Green Version]
  9. Hsu, S.-C.; Chung, J.-G. Anticancer potential of emodin. BioMedicine 2012, 2, 108–116. [Google Scholar] [CrossRef]
  10. Benova, B.; Adam, M.; Pavlikova, P.; Fischer, J. Supercritical fluid extraction of piceid, resveratrol and emodin from Japanese knotweed. J. Supercrit. Fluids 2010, 51, 325–330. [Google Scholar] [CrossRef]
  11. Bashtanova, U.B.; Beckett, K.P.; Flowers, T.J. Review: Physiological Approaches to the Improvement of Chemical Control of Japanese Knotweed (Fallopia japonica). Weed Sci. 2009, 57, 584–592. [Google Scholar] [CrossRef]
  12. Chen, H.; Tuck, T.; Ji, X.; Zhou, X.; Kelly, G.; Cuerrier, A.; Zhang, J. Quality Assessment of Japanese Knotweed (Fallopia japonica) Grown on Prince Edward Island as a Source of Resveratrol. J. Agric. Food Chem. 2013, 61, 6383–6392. [Google Scholar] [CrossRef] [PubMed]
  13. Chang, S.-J.; Huang, S.-H.; Lin, Y.-J.; Tsou, Y.-Y.; Lin, C.-W. Antiviral activity of Rheum palmatum methanol extract and chrysophanol against Japanese encephalitis virus. Arch. Pharm. Res. 2014, 37, 1117–1123. [Google Scholar] [CrossRef]
  14. Dong, X.; Fu, J.; Yin, X.; Cao, S.; Li, X.; Lin, L.; Huyiligeqi; Ni, J. Emodin: A Review of Its Pharmacology, Toxicity and Pharmacokinetics. Phytother. Res. 2016, 30, 1207–1218. [Google Scholar] [CrossRef] [PubMed]
  15. Tuli, H.S.; Aggarwal, V.; Tuorkey, M.; Aggarwal, D.; Parashar, N.C.; Varol, M.; Savla, R.; Kaur, G.; Mittal, S.; Sak, K. Emodin: A metabolite that exhibits anti-neoplastic activities by modulating multiple oncogenic targets. Toxicol. In Vitro 2021, 73, 105142. [Google Scholar] [CrossRef] [PubMed]
  16. Shrimali, D.; Shanmugam, M.K.; Kumar, A.P.; Zhang, J.; Tan, B.K.H.; Ahn, K.S.; Sethi, G. Targeted abrogation of diverse signal transduction cascades by emodin for the treatment of inflammatory disorders and cancer. Cancer Lett. 2013, 341, 139–149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Li, Q.; Gao, J.; Pang, X.; Chen, A.; Wang, Y. Molecular Mechanisms of Action of Emodin: As an Anti-Cardiovascular Disease Drug. Front. Pharmacol. 2020, 11, 559607. [Google Scholar] [CrossRef] [PubMed]
  18. Cui, Y.; Chen, L.-J.; Huang, T.; Ying, J.-Q.; Li, J. The pharmacology, toxicology and therapeutic potential of anthraquinone derivative emodin. Chin. J. Nat. Med. 2020, 18, 425–435. [Google Scholar] [CrossRef]
  19. Schwarz, S.; Wang, K.; Yu, W.; Sun, B.; Schwarz, W. Emodin inhibits current through SARS-associated coronavirus 3a protein. Antivir. Res. 2011, 90, 64–69. [Google Scholar] [CrossRef] [PubMed]
  20. Yu, T.; Guo, M.; Wen, S.; Zhao, R.; Wang, J.; Sun, Y.; Liu, Q.; Zhou, H. Poly(ethylene glycol) dimethyl ether mediated oxidative scission of aromatic olefins to carbonyl compounds by molecular oxygen. RSC Adv. 2021, 11, 13848–13852. [Google Scholar] [CrossRef]
  21. Li, B.; Zhao, X.; Zhang, L.; Cheng, W. Emodin Interferes with AKT1-Mediated DNA Damage and Decreases Resistance of Breast Cancer Cells to Doxorubicin. Front. Oncol. 2021, 10, 58853. [Google Scholar] [CrossRef]
  22. Paudel, P.; Shrestha, S.; Park, S.E.; Seong, S.H.; Fauzi, F.M.; Jung, H.A.; Choi, J.S. Emodin Derivatives as Multi-Target-Directed Ligands Inhibiting Monoamine Oxidase and Antagonizing Vasopressin V1A Receptors. ACS Omega 2020, 5, 26720–26731. [Google Scholar] [CrossRef]
  23. Paudel, P.; Seong, S.H.; Fauzi, F.M.; Bender, A.; Jung, H.A.; Choi, J.S. Establishing GPCR Targets of hMAO Active Anthraquinones from Cassia obtusifolia Linn Seeds Using In Silico and In Vitro Methods. ACS Omega 2020, 5, 7705–7715. [Google Scholar] [CrossRef] [Green Version]
  24. Teich, L.; Daub, K.S.; Krügel, V.; Nissler, L.; Gebhardt, R.; Eger, K. Synthesis and biological evaluation of new derivatives of emodin. Bioorg. Med. Chem. 2004, 12, 5961–5971. [Google Scholar] [CrossRef]
  25. Wang, X.P.; Xu, W.F. Facile Synthesis of Emodin Derivatives as Potential MMPIs. Bull. Korean Chem. Soc. 2005, 26, 1923–1924. [Google Scholar] [CrossRef]
  26. Zhu, X.; Chen, Q.; Yang, Y.; Ai, X.; Chen, S.; Song, Y. Synthesis and anti-inflammatory effects of novel emodin derivatives bearing azole moieties. Arch. Pharm. 2020, 353, 1900264. [Google Scholar] [CrossRef] [PubMed]
  27. Yang, K.; Jin, M.-J.; Quan, Z.-S.; Piao, H.-R. Design and Synthesis of Novel Anti-Proliferative Emodin Derivatives and Studies on Their Cell Cycle Arrest, Apoptosis Pathway and Migration. Molecules 2019, 24, 884. [Google Scholar] [CrossRef] [Green Version]
  28. Duan, F.; Li, X.; Cai, S.; Xin, G.; Wang, Y.; Du, D.; He, S.; Huang, B.; Guo, X.; Zhao, H.; et al. Haloemodin as Novel Antibacterial Agent Inhibiting DNA Gyrase and Bacterial Topoisomerase I. J. Med. Chem. 2014, 57, 3707–3714. [Google Scholar] [CrossRef] [PubMed]
  29. Koerner, S.K.; Hanai, J.-I.; Bai, S.; Jernigan, F.E.; Oki, M.; Komaba, C.; Shuto, E.; Sukhatme, V.P.; Sun, L. Design and synthesis of emodin derivatives as novel inhibitors of ATP-citrate lyase. Eur. J. Med. Chem. 2017, 126, 920–928. [Google Scholar] [CrossRef] [PubMed]
  30. Chalothorn, T.; Rukachaisirikul, V.; Phongpaichit, S.; Pannara, S.; Tansakul, C. Synthesis and antibacterial activity of emodin and its derivatives against methicillin-resistant Staphylococcus aureus. Tetrahedron Lett. 2019, 60, 151004. [Google Scholar] [CrossRef]
  31. Duan, F.; Xin, G.; Niu, H.; Huang, W. Chlorinated emodin as a natural antibacterial agent against drug-resistant bacteria through dual influence on bacterial cell membranes and DNA. Sci. Rep. 2017, 7, 12721. [Google Scholar] [CrossRef] [Green Version]
  32. Xiong, H.-R.; Luo, J.; Hou, W.; Xiao, H.; Yang, Z.-Q. The effect of emodin, an anthraquinone derivative extracted from the roots of Rheum tanguticum, against herpes simplex virus in vitro and in vivo. J. Ethnopharmacol. 2011, 133, 718–723. [Google Scholar] [CrossRef] [PubMed]
  33. Hsiang, C.Y.; Ho, T.Y. Emodin is a novel alkaline nuclease inhibitor that suppresses herpes simplex virus type 1 yields in cell cultures. Br. J. Pharmacol. 2008, 155, 227–235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Shuangsuo, D.; Zhengguo, Z.; Yunru, C.; Xin, Z.; Baofeng, W.; Lichao, Y.; Yan’an, C. Inhibition of the replication of hepatitis B virus in vitro by emodin. Med. Sci. Monit. 2006, 12, Br302–Br306. [Google Scholar]
  35. Parvez, M.K.; Al-Dosari, M.S.; Alam, P.; Rehman, M.; Alajmi, M.F.; Alqahtani, A.S. The anti-hepatitis B virus therapeutic potential of anthraquinones derived from Aloe vera. Phytother. Res. 2019, 33, 2960–2970. [Google Scholar] [CrossRef] [PubMed]
  36. Barnard, D.L.; Huffman, J.H.; Morris, J.L.B.; Wood, S.G.; Hughes, B.G.; Sidwell, R.W. Evaluation of the antiviral activity of anthraquinones, anthrones and anthraquinone derivatives against human cytomegalovirus. Antivir. Res. 1992, 17, 63–77. [Google Scholar] [CrossRef]
  37. Dai, J.-P.; Wang, Q.-W.; Su, Y.; Gu, L.-M.; Zhao, Y.; Chen, X.-X.; Chen, C.; Li, W.-Z.; Wang, G.-F.; Li, K.-S. Emodin Inhibition of Influenza A Virus Replication and Influenza Viral Pneumonia via the Nrf2, TLR4, p38/JNK and NF-kappaB Pathways. Molecules 2017, 22, 1754. [Google Scholar] [CrossRef] [Green Version]
  38. Batista, M.N.; Braga, A.C.S.; Campos, G.R.F.; Souza, M.M.; Matos, R.P.A.; Lopes, T.Z.; Candido, N.M.; Lima, M.L.D.; Machado, F.C.; Andrade, S.T.Q.d.; et al. Natural Products Isolated from Oriental Medicinal Herbs Inactivate Zika Virus. Viruses 2019, 11, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Zhang, H.M.; Wang, F.; Qiu, Y.; Ye, X.; Hanson, P.; Shen, H.; Yang, D. Emodin inhibits coxsackievirus B3 replication via multiple signalling cascades leading to suppression of translation. Biochem. J. 2016, 473, 473–485. [Google Scholar] [CrossRef] [PubMed]
  40. Ding, Y.; Xu, J.; Cheng, L.-B.; Huang, Y.-Q.; Wang, Y.-Q.; Li, H.; Li, Y.; Ji, J.-Y.; Zhang, J.-H.; Zhao, L. Effect of Emodin on Coxsackievirus B3m-Mediated Encephalitis in Hand, Foot, and Mouth Disease by Inhibiting Toll-Like Receptor 3 Pathway In Vitro and In Vivo. J. Infect. Dis. 2020, 222, 443–455. [Google Scholar] [CrossRef] [PubMed]
  41. Semple, S.J.; Pyke, S.M.; Reynolds, G.D.; Flower, R.L.P. In vitro antiviral activity of the anthraquinone chrysophanic acid against poliovirus. Antiv. Res. 2001, 49, 169–178. [Google Scholar] [CrossRef]
  42. Wang, Z.; Zheng, N.; Liang, J.; Wang, Q.; Zu, X.; Wang, H.; Yuan, H.; Zhang, R.; Guo, S.; Liu, Y.; et al. Emodin resists to Cyprinid herpesvirus 3 replication via the pathways of Nrf2/Keap1-ARE and NF-κB in the ornamental koi carp (Cyprinus carpio haematopterus). Comp. Biochem. Physiol. Part C: Toxicol. Pharmacol. 2021, 246, 109023. [Google Scholar] [CrossRef] [PubMed]
  43. Adhikari, B.; Marasini, B.P.; Rayamajhee, B.; Bhattarai, B.R.; Lamichhane, G.; Khadayat, K.; Adhikari, A.; Khanal, S.; Parajuli, N. Potential roles of medicinal plants for the treatment of viral diseases focusing on COVID-19: A review. Phytother. Res. 2021, 35, 1298–1312. [Google Scholar] [CrossRef] [PubMed]
  44. Vougogiannopoulou, K.; Corona, A.; Tramontano, E.; Alexis, M.N.; Skaltsounis, A.-L. Natural and Nature-Derived Products Targeting Human Coronaviruses. Molecules 2021, 26, 448. [Google Scholar] [CrossRef]
  45. Remali, J.; Aizat, W.M. A Review on Plant Bioactive Compounds and Their Modes of Action Against Coronavirus Infection. Front. Pharmacol. 2021, 11, 589044. [Google Scholar] [CrossRef] [PubMed]
  46. Illian, D.N.; Siregar, E.S.; Sumaiyah, S.; Utomo, A.R.; Nuryawan, A.; Basyuni, M. Potential compounds from several Indonesian plants to prevent SARS-CoV-2 infection: A mini-review of SARS-CoV-2 therapeutic targets. Heliyon 2021, 7, e06001. [Google Scholar] [CrossRef] [PubMed]
  47. Zhou, Y.; Hou, Y.; Shen, J.; Huang, Y.; Martin, W.; Cheng, F. Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2. Cell Discov. 2020, 6, 14. [Google Scholar] [CrossRef] [Green Version]
  48. Islam, M.T.; Sarkar, C.; El-Kersh, D.M.; Jamaddar, S.; Uddin, S.J.; Shilpi, J.A.; Mubarak, M.S. Natural products and their derivatives against coronavirus: A review of the non-clinical and pre-clinical data. Phytother. Res. 2020, 34, 2471–2492. [Google Scholar] [CrossRef]
  49. Ho, T.-Y.; Wu, S.-L.; Chen, J.-C.; Li, C.-C.; Hsiang, C.-Y. Emodin blocks the SARS coronavirus spike protein and angiotensin-converting enzyme 2 interaction. Antivir. Res. 2007, 74, 92–101. [Google Scholar] [CrossRef] [PubMed]
  50. Das, S.; Singha Roy, A. Naturally Occurring Anthraquinones as Potential Inhibitors of SARS-CoV-2 Main Protease: A Molecular Docking Study. In ChemRxiv; Cambridge Open Engage: Cambridge, UK, 2020. [Google Scholar] [CrossRef]
  51. Singh, J.; Pandit, P.; McArthur, A.G.; Banerjee, A.; Mossman, K. Evolutionary trajectory of SARS-CoV-2 and emerging variants. Virol. J. 2021, 18, 166. [Google Scholar] [CrossRef] [PubMed]
  52. Abdul-Rasool, S.; Fielding, B.C. Understanding Human Coronavirus HCoV-NL63. Open Virol. J. 2010, 4, 76–84. [Google Scholar] [CrossRef] [PubMed]
  53. Chakraborty, A.; Diwan, A. NL63: A Better Surrogate Virus for studying SARS-CoV-2. Integr. Mol. Med. 2020, 7, 1–9. [Google Scholar] [CrossRef]
  54. Muzychkina, R.A.; Pribytkova, L.N. Bromination of emodin. Chem. Nat. Compd. 1990, 26, 524–527. [Google Scholar] [CrossRef]
Molecules 26 06825 g001 550
Figure 1. Derivatization of emodin.
Figure 1. Derivatization of emodin.
Molecules 26 06825 g001
Molecules 26 06825 g002 550
Figure 2. Synthesis of amino emodin E_NH2.
Figure 2. Synthesis of amino emodin E_NH2.
Molecules 26 06825 g002
Molecules 26 06825 g003 550
Figure 3. Synthesis of E_SO3H.
Figure 3. Synthesis of E_SO3H.
Molecules 26 06825 g003
Molecules 26 06825 g004 550
Figure 4. Emodin and emodin derivatives prepared in the course of this study.
Figure 4. Emodin and emodin derivatives prepared in the course of this study.
Molecules 26 06825 g004
Molecules 26 06825 g005 550
Figure 5. IC50 curves reflecting the effect of the tested compounds on Vero cell viability.
Figure 5. IC50 curves reflecting the effect of the tested compounds on Vero cell viability.
Molecules 26 06825 g005
Molecules 26 06825 g006 550
Figure 6. IC50 curves representing the anti-HCoV-NL63 effects of emodin and emodin derivatives.
Figure 6. IC50 curves representing the anti-HCoV-NL63 effects of emodin and emodin derivatives.
Molecules 26 06825 g006
Table
Table 1. Iodination of emodin E_H with I2/NaHCO3, NIS or I2/H2O2.
Table 1. Iodination of emodin E_H with I2/NaHCO3, NIS or I2/H2O2.
Molecules 26 06825 i001
EntryReagent (equiv.)SolventConditionsRelative Distribution a
E_HE_IE_2IE_3I
1I2 (2), NaHCO3THF/H2Ort., 1 h-100 (73%)--
2I2 (4), NaHCO3THF/H2Ort., 1 h-100--
3I2 (10), NaHCO3THF/H2Ort., 24 h---100 (81%)
4I2 (10), NaHCO3THF/H2O60 °C, 24 h---100
5NIS (1)THFrt., 24 h397--
6NIS (1.3)THFrt., 3 h-100 (84%)--
7NIS (1.3)THFrt., 24 h-937-
8NIS (2)THFrt., 24 h-8614-
9NIS (4)THFrt., 24 h--5644
10NIS (4)THF60 °C, 24 h---100 (79%)
11I2 (1)2-MeTHFrt., 24 h100---
12I2 (0.5), H2O2 (8)2-MeTHFrt., 24 h5644--
13I2 (1), H2O2 (8)2-MeTHFrt., 24 h2179--
14I2 (1.5), H2O2 (2.5)2-MeTHFrt., 24 h1684--
15I2 (1.5), H2O2 (4)2-MeTHFrt., 24 h1387--
16I2 (2), H2O2 (2.5)2-MeTHFrt., 24 h-100 (93%)--
17I2 (8), H2O2 (8)2-MeTHFrt., 24 h-100--
Reaction conditions: Emodin (0.1 mmol), reagent (I2 (0.2–1.0 mmol), NIS (N-iodosuccinimide) (0.1–0.4 mmol), H2O2 (30%, 0.25–0.8 mmol)), solvent (1 mL), a Conversion to product was determined by 1H NMR.
Table
Table 2. Bromination of emodin E_H with NBS and HBr/H2O2.
Table 2. Bromination of emodin E_H with NBS and HBr/H2O2.
Molecules 26 06825 i002
EntryReagent (Equiv.)SolventConditionsRelative Distribution a
E_HE_BrE_2BrE_3BrE_4Br
1NBS (1)THF0 °C, 10 min344026--
2NBS (1.5)THF0 °C, 15 min174340--
3NBS (2.2)THF0 °C, 30 min--100 (88%)--
4NBS (3)THF0 °C, 24 h--6337-
5NBS (3)THFrt., 24 h--2773-
6NBS (4)THFrt., 24 h---6040
7NBS (5)THFrt., 24 h----100 (83%)
8HBr (1), H2O2 (2.5)TFE0 °C, 2 h2971---
9HBr (1.3), H2O2 (5)TFE0 °C, 2 h67312--
10HBr (2.5), H2O2 (5)TFE0 °C, 2 h88111--
11HBr (1), H2O2 (5)TFErt., 24 h97912--
12HBr (2), H2O2 (5)TFErt., 24 h-2872--
13HBr (2.4), H2O2 (5)TFErt., 24 h-1486--
14HBr (4), H2O2 (5)TFErt., 24 h--100 (91%)--
Reaction conditions: Emodin (0.1 mmol), reagent (NBS (N-bromosuccinimide) (0.1–0.5 mmol), HBr (48%, 0.1–0.4 mmol), H2O2 (30%, 0.25–0.5 mmol)), solvent (1 mL). a Conversion to product was determined by 1H NMR.
Table
Table 3. Chlorination of emodin E_H with NCS and HCl/H2O2.
Table 3. Chlorination of emodin E_H with NCS and HCl/H2O2.
Molecules 26 06825 i003
EntryReagent (Equiv.)SolventConditionsRelative Distribution a
E_HE_ClE_2ClE_3Cl
1NCS (1)THFrt., 24 h100---
2NCS (2)THFrt., 24 h49429-
3NCS b (2.5)THFrt., 3 h--100 (81%)-
4NCS b (6)THF60 °C, 24 h--100-
5HCl (1), H2O2 (2)TFErt., 24 h56368-
6HCl (2), H2O2 (5)TFErt., 24 h106822-
7HCl (3), H2O2 (5)TFErt., 24 h-443917
8HCl (4), H2O2 (5)TFErt., 24 h-413722
9HCl (5), H2O2 (10)TFErt., 24 h-175033
Reaction conditions: Emodin (0.1 mmol), reagent (NCS (N-chlorosuccinimide) (0.1–0.6 mmol), HCl (37%, 0.1–0.5 mmol), H2O2 (30%, 0.2–1.0 mmol)), solvent (1 mL), a Conversion to product was determined by 1H NMR, b H2SO4.
Table
Table 4. Nitration of emodin with HNO3 and H2SO4.
Table 4. Nitration of emodin with HNO3 and H2SO4.
EntryReagent (Equiv.) aConditionsConv. b
E_XE_4NO2
1HNO3 (1), H2SO4 (4), MeCN0 °C, 1 h2%-
2HNO3 (2.5), H2SO4 (10)0 °C, 1 h78%22%
3HNO3 (2.5), H2SO4 (20)0 °C, 1 h79%21%
4HNO3 (3), H2SO4 (10)0 °C, 1 h59%41%
5HNO3 (3), H2SO4 (20)0 °C, 1 h27%73%
6HNO3 (4.4), H2SO4 (20)0 °C, 1 h10%90%
7HNO3 (6.2), H2SO4 (5)0 °C, 0.25 h13%87%
8HNO3 (6.2), H2SO4 (10)0 °C, 0.5 h-100 (91%)
9HNO3 (63), H2SO4 (80)0 °C, 1 h; rt, 4 h-100 % [24]
a Emodin (0.1 mmol), reagent (HNO3 (0.1–6.3 mmol), H2SO4 (0.4–8.0 mmol)), MeCN/no solvent, b Conversion to product was determined by 1H NMR relative to E_H. E_X is a mixture of mono-, di- and trinitroemodin.
Table
Table 5. IC50 (µM) values corresponding to the anti HCoV-NL63 effects (AV) and the effects on Vero cell viability (expressed as IC50 values, CV) of each of the tested compounds.
Table 5. IC50 (µM) values corresponding to the anti HCoV-NL63 effects (AV) and the effects on Vero cell viability (expressed as IC50 values, CV) of each of the tested compounds.
IC50 (µM)
Exp.CT
E_H2.54.9
E_4NO26.16.1
E_I1.33.6
E_3I0.54.9
E_4Br1.77.2
E_2Br1.05.4
E_2Cl1.17.5
E_OMe>508.7
E_SO3H22.0>50
E_NH26.341.8
E_Br_OMe>50>50
Remdesivir0.61
Chloroquine19.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Horvat, M.; Avbelj, M.; Durán-Alonso, M.B.; Banjanac, M.; Petković, H.; Iskra, J. Antiviral Activities of Halogenated Emodin Derivatives against Human Coronavirus NL63. Molecules 2021, 26, 6825. https://doi.org/10.3390/molecules26226825

AMA Style

Horvat M, Avbelj M, Durán-Alonso MB, Banjanac M, Petković H, Iskra J. Antiviral Activities of Halogenated Emodin Derivatives against Human Coronavirus NL63. Molecules. 2021; 26(22):6825. https://doi.org/10.3390/molecules26226825

Chicago/Turabian Style

Horvat, Monika, Martina Avbelj, María Beatriz Durán-Alonso, Mihailo Banjanac, Hrvoje Petković, and Jernej Iskra. 2021. "Antiviral Activities of Halogenated Emodin Derivatives against Human Coronavirus NL63" Molecules 26, no. 22: 6825. https://doi.org/10.3390/molecules26226825

APA Style

Horvat, M., Avbelj, M., Durán-Alonso, M. B., Banjanac, M., Petković, H., & Iskra, J. (2021). Antiviral Activities of Halogenated Emodin Derivatives against Human Coronavirus NL63. Molecules, 26(22), 6825. https://doi.org/10.3390/molecules26226825

Article Metrics

Back to TopTop