Umbelliferone and eriodictyol suppress the cellular entry of SARS-CoV-2
Cell & Bioscience volume 13, Article number: 118 (2023)
Artemisia argyi (A. argyi), also called Chinese mugwort, has been widely used to control pandemic diseases for thousands of years since ancient China due to its anti-microbial infection, anti-allergy, and anti-inflammation activities. Therefore, the potential of A. argyi and its constituents in reducing the infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was investigated in this study.
Among the phytochemicals in A. argyi, eriodictyol and umbelliferone were identified to target transmembrane serine protease 2 (TMPRSS2) and angiotensin-converting enzyme 2 (ACE2) proteins, the essential factors for the cellular entry of SARS-CoV-2, in both FRET-based enzymatic assays and molecular docking analyses. These two ingredients of A. argyi suppressed the infection of ACE2-expressed HEK-293 T cells with lentiviral-based pseudo-particles (Vpp) expressing wild-type and variants of SARS-CoV-2 spike (S) protein (SARS-CoV-2 S-Vpp) via interrupting the interaction between S protein and cellular receptor ACE2 and reducing the expressions of ACE2 and TMPRSS2. Oral administration with umbelliferone efficiently prevented the SARS-CoV-2 S-Vpp-induced inflammation in the lung tissues of BALB/c mice.
Eriodictyol and umbelliferone, the phytochemicals of Artemisia argyi, potentially suppress the cellular entry of SARS-CoV-2 by preventing the protein binding activity of the S protein to ACE2.
The global pandemic of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a serious health-threatening issue in most countries . The genomic evolution of SARS-CoV-2 creates diverse mutant strains, which become more transmittable and resistant to immune attack and anti-virus treatments . SARS-CoV-2 variants were classified into variants of interest (VOIs) and variants of concern (VOCs) by the World Health Organization (WHO) according to their transmission rate, disease severity, therapeutic response, and healthy problem . For example, the Omicron variant, identified as the lineage B.1.1.529 in South Africa in November 2021, has been viewed as one of the circulating and high lethal VOIs associated with its rapid human-to-human transmission, increased risk of reinfection, and resistance to vaccines [4,5,6].
SARS-CoV-2 belongs to the family of highly diverse coronaviruses (CoVs) constituted by positive-sense and single-strand RNA . For cellular infection, SARS-CoV-2 utilizes various factors for the entry of host cells, translation of RNA genome, and replication cycle, ultimately releasing viral progeny. The cellular access of SARS-CoV-2 relies on its receptor engagement of spike (S) protein with the specific host receptor angiotensin-converting enzyme2 (ACE2) . The receptor binding domain (RDB) of the S1 subunit of S protein directly mediates viral attachment to ACE2, while S2 is cleaved by the membrane proteases and exposed the fusion peptide (FP) for ACE2 engagement and membrane fusion [9, 10]. Transmembrane serine protease 2 (TMPRSS2) and cathepsin L are typically necessary to present the FP for membrane fusion with host cells [9, 11]. In addition, the S protein can be cleaved into S1 and S2 subunits by the proprotein convertase furin for the anchorage at the cellular membrane during the viral biosynthesis in the virus-producing cells, facilitating the virus entry [12, 13]. SARS-CoV-2 next undergoes clathrin-mediated endocytosis to gain access into cells . Following virus entry, the viruses release, replicate, and transcribe their genomic RNA following cellular entry to produce essential structural proteins that compose new and complete viral particles [14, 15]. Since SARS-CoV-2 has evolved into diverse variants, the therapeutic effectiveness of anti-virus treatments is limited and scarce. It is imperative to explore the potential protective treatments or integrated approaches against COVID-19.
Artemisia argyi (A. argyi), also named Chinese mugwort, belongs to the family of Asteraceae. This herb has been used to control pandemic diseases by burning dried herbs for thousands of years since ancient China  and still is a renowned traditional Chinese medicine (TCM) and food currently used in the Far East. A. argyi contains a broad-spectrum profile of bioactive compounds, including flavonoids, terpenoids, and caffeoylquinic acids, and contributes to anti-microbial, anti-allergic, anti-diabetic, and anti-inflammatory activities by modulating the immune system or defending oxidative stress [17, 18]. In addition, this TCM was also shown to suppress the proliferation of tumor cells with low cytotoxicity to normal cells [19, 20]. Recently, A. argyi has been implied as a potential TCM against severe acute respiratory syndrome (SARS), middle east respiratory syndrome (MERS), and COVID-19 [21, 22]. Nevertheless, its activity against the infection with SARS-CoV2 has not been demonstrated yet, and the active phytochemicals and underlying mechanisms also remain unclear.
Eriodictyol and umbelliferone, two constituents in Artemisia argyi, suppressed the enzymatic activities involved in the cellular entry of SARS-CoV-2
To identify the bioactive compounds with the potential against the enzymatic activities for S protein priming or viral replication, 14 known ingredients of A. argyi  were subjected to fluorescence resonance energy transfer (FRET)-based enzymatic activity assays, which principle was illustrated in Fig. 1A . As shown in Table 1, the enzymatic activity of TMPRSS2 was dramatically attenuated by eriodictyol (100.0% ± 0.3%), umbelliferone (74.1% ± 3.5%), and 13-Oxo-9E,11E-octadecadienoic acid (13-Oxo-ODE) (76% ± 4.7%) at 60 μM. However, none of these compounds at 60 μM suppressed more than 50% of 3CLpro (Mpro) activity (Table 1), which is the main protease in cleaving coronavirus polyprotein for the replication of SARS-CoV . 13-Oxo-ODE is a metabolite from de-hydrogenation of 13-hydroxyoctadecadienoic acid (13-HODE), an oxidized product of linoleic acid catalyzed by 12/15-lipoxygenase (LOX) [24, 25], and has been reported as an endogenous ligand for nuclear hormone receptor peroxisome proliferation-activated receptor gamma (PPARγ) involved in inflammatory bowel disease [26, 27]. But the clinical relevance of 13-Oxo-ODE as a potential contributor to nonalcoholic fatty liver disease  led us to focus on eriodictyol and umbelliferone for their anti-corona virus activity. Eriodictyol not only repressed the enzymatic activity of TMPRSS2 (Fig. 1B) but also abolished furin activity (Fig. 1C) in a dose-dependent manner. Following the priming of the S protein by TMPRSS2 or furin, the cleaved S protein interacts with the human cellular receptor ACE2 . Another FRET assay, which principle was illustrated in Fig. 1D , further showed the inhibitory effect of eriodictyol on the interaction between S protein and ACE2 in a dose-dependent manner (Fig. 1E). Taken together, eriodictyol and umbelliferone in A. argyi potentially suppressed the enzymatic activities of TMPRSS2 or furin to reduce the S1/S2 priming of S protein and interfered the interaction between S protein and ACE2.
Molecular docking analysis of eriodictyol and umbelliferone with SARS-CoV-2 proteins
Since eriodictyol and umbelliferone suppressed the enzymatic activity of TMPRSS2 and S protein/ACE2 interaction, the binding modes of these two compounds in the catalytic pockets of TMPRSS2 and the interface between the S protein and ACE2 were performed with Discovery Studio. The results indicated that eriodictyol (Fig. 2A) and umbelliferone (Fig. 2B) showed binding affinity to the catalytic pockets of TMPRSS2, and the changes in energy value were − 13.05 kcal/mol (eriodictyol/TMPRSS2) and − 14.82 kcal/mol (umbelliferone/TMPRSS2), respectively (Table 2). Furthermore, the interface between the S protein and human receptor ACE2 was also disrupted by these two compounds (Figs. 2C, D) with the free energy decreases of − 33.478 and -36.539 kcal/mol for eriodictyol/S protein-ACE2 and umbelliferone/S protein-ACE2 complexes, respectively (Table 2). Based on these results, eriodictyol and umbelliferone could be the phytochemicals in A. argy capable of targeting the potential residues (Table 2) of multiple proteins associated with the cellular entry of SARS-CoV-2.
A. argyi and its active compounds downregulate ACE2 and TMPRSS2 expressions in lung epithelial cells
The cytotoxic activities of A. argyi, eriodictyol, and umbelliferone were next examined in Beas 2B lung epithelial cells using MTT assays. Cytotoxicity concentration 50% (CC50) values of A. argyi, eriodictyol, and umbelliferone were 4471 μg/ml (Additional file 1: Fig. S1A), 212.6 μM (Additional file 1: Fig. S1B), and 302.4 μM (Additional file 1: Fig. S1C), respectively. Based on the level of CC50, we further addressed whether the expressions of ACE2 and TMPRSS2 are also affected by these potential agents in preventing the cellular entry of SARS-CoV-2. A. Argyi, eriodictyol, and umbelliferone significantly suppressed the and mRNA expressions of ACE2 and TMPRSS2 (Fig. 3A–D). While A. argyi suppressed ACE2 protein expression at 100 μg/ml (Fig. 3E) and eriodictyol reduced TMPRSS2 protein level at 50 μM (Fig. 3F), umbelliferone showed more potent inhibitory effects on both ACE2 and TMPRSS2 protein expressions in a dose-dependent manner (Fig. 3G). These results suggest that A. argyi, eriodictyol, and umbelliferone can inhibit SARS-CoV-2 infection by mainly repressing the expressions of cellular ACE2 and TMPRSS2 with different specificity and potency.
A. argyi, eriodictyol, and umbelliferone broadly suppress the cell entry of SARS-CoV-2 S-Vpp variants
The above data suggest that eriodictyol and umbelliferone could prevent the cellular entry of SARS-CoV-2 through multiple mechanisms, including inhibitions of cellular TMPRSS2 activity, cellular ACE2 and TMPRSS2 expressions, and S protein/ACE2 interaction. Therefore, we further assessed the efficacy of A. argyi extracts and its phytochemicals eriodictyol and umbelliferone in blocking the cellular entry of SARS-CoV-2 S protein-pseudo-type lentiviral particles (SARS-CoV-2 S-Vpp). HEK-293 T cells stably expressing ACE2 protein were employed to test the cellular entry of SARS-CoV-2 Vpp with S protein variants following the pretreatments with tested drugs. The results displayed that the infections with most SARS-CoV-2 S-Vpp variants, including wild type (WT) (Fig. 4A), B.1.351 (Beta) (Fig. 4C), Lineage P1 (Gamma) (Fig. 4D), and B1.617.2 (Delta) (Fig. 4E) were dramatically repressed by the treatments with A. argyi, eriodictyol, and umbelliferone in a dose-dependent manner without affecting the infection with VSVG control (Additional file 2: Fig. S2) and the cell viability. However, the infections with B1.1.7 (Alpha) (Fig. 4B) and B.1.429 (Epsilon) (Fig. 4F) were not affected by A. argyi extracts. The inhibitory concentration 50% (IC50) values of A. argyi, eriodictyol, and umbelliferone in antiviral activity are shown in Table 3. The IC50 values of A. argyi, eriodictyol, and umbelliferone in the detected viruses, except B1.1.7 (Alpha) and B.1.429 (Epsilon), were around 332 ~ 1534 μg/ml, 88 ~ 126 μM, and 67 ~ 100 μM, respectively. These findings suggested that A. argyi can inhibit infections with SARS-CoV-2.
Omicron (B.1.1.529), first reported in November 2021, is a highly mutated variant of SARS-CoV-2 and has become the predominant viral variant of concern (VOC) worldwide. The inhibitory effects of tested compounds on the cellular entry of SARS-CoV-2 S-Vpp of Omicron variants were next examined in HEK-293 T cells expressing ACE2 proteins. The results showed that umbelliferone significantly and dose-dependently suppressed the cellular entry of all tested Omicron variants BA.1 (Fig. 5A), BA.1.1 (Fig. 5B), BA.2 (Fig. 5C), and XBB.1 (Fig. 5D). A. argyi also significantly reduced the infections with BA.1, BA1.1, and XBB.1 variants. However, eriodictyol only slightly interfered with the infection with XBB.1 but not other Omicron variants. Our findings suggested that umbelliferone can be a potential agent more selective against the variants of SARS-CoV-2 Omicron with IC50 values of around 50 ~ 147 μM (Table 3).
Umbelliferone effectively suppressed the lung inflammation of mice induced by SARS-CoV-2 S-Vpp of WT or Omicron variants
Since umbelliferone showed promising activity against the cellular entry of SARS-CoV-2 S-Vpp in vitro, its inhibitory effects on SARS-CoV-2-associated inflammation in lung tissues were addressed in a mice model. As illustrated in Fig. 6A, BALB/c mice expressing human ACE2 were pre-treated with umbelliferone 70 mg/kg for 2 days, followed by challenge with WT, BA.2, or XBB.1 variants of SARS-CoV-2 S-Vpp and continuous treatment with umbelliferone for another 7 days. Compared to the infection of SARS-CoV-2 S-Vpp WT (Additional file 3: Fig. S3A), infections with BA.2 (Fig. 6B) and XBB.1 (Additional file 4: Fig. S4A) variants more obviously elicited inflammation as evidenced by the formation of granuloma (indicated with the arrow) in H&E staining and the inductions of inflammatory gene expressions, including Tnf, Infg, Il1b, Il6, Il10, and Tgfb1 in RT-qPCR analysis (Additional file 5: Figs. S5A–F). IHC staining with anti-CD3 and TNF-α antibodies further supported the SARS-CoV-2 S-Vpp-induced infiltration of immune cells and inflammation (Fig. 6B and Additional files 3, 4: Figs. S3A and S4A). The increases in the mRNA level of alpha-smooth muscle actin (α-SMA), a marker for fibrogenesis , further revealed the tendency of lung tissues to form lung fibrosis by the infection with the pseudo viral particles (Additional file 5: Fig. S5G). Similar to the effect in vitro (Fig. 3), pretreatment with umbelliferone suppressed the protein and mRNA expression of TMPRSS2 in lung tissues (Figs. 6C, D, and Additional files 3, 4: Figs. S3B, C and S4B, C). Also, it reduced the granuloma formation, CD3-positive immune cell infiltration, and TNF-α expression (Fig. 6B and Additional file 3: Figs. S3A and Additional file 4: S4A). The decreases in the mRNA expressions of inflammatory genes (Fig. 6E–J and Additional files 3, 4: Figs. S3D–I and S4D–I) and α-SMA expressions (Fig. 6K and Additional files 3, 4: Figs. S3J and S4J) further supported the potential of umbelliferone in preventing infection with SARS-CoV-2 variants in vivo.
Eriodictyol decreases the replication of SARS-CoV-2 by enzymatic inhibition of RdRp
After the entry and uncoating of SARS-CoV-2, 3-chymotrypsin-like protease (3CLpro) and papain-like protease (PLpro) are the virus-encoded cysteine proteases for cleaving the viral polyprotein and then generating nonstructural proteins during viral replication [30, 31]. RNA-dependent RNA polymerase (RdRp; also known as nsp12) is a viral replicase for synthesizing complementary RNA during the transcriptional cycle of SARS-CoV-2 . Therefore, we further examined the inhibitory effects of tested compounds on the enzymatic activity of 3CLpro, PLpro, and, RdRp by performing FRET-base enzymatic activity assays and found that eriodictyol dramatically reduced the activities of these three viral enzymes (Additional file 6: Figs. S6A–C). However, the 50% inhibition concentrations of eriodictyol for 3CLpro and PLpro were over 240 μM. These results displayed that A. argyi and its phytochemicals may interfere with the cellular entry and RNA replication of SARS-CoV-2 at relatively lower doses and showed inhibitory effects on the viral protease activities at higher doses.
Until the end of 2022, the COVID-19 pandemic has caused more than 60 million confirmed cases and 6 million deaths worldwide. The emergency of SARS-CoV-2 variants, which not only facilitate viral replication, transmission, and immune escape but also weaken the protective ability of recently developed vaccines [2, 33], generates more concerns about the prevention of SARS-CoV-2 infection. TCMs have been explored as potential therapeutics against COVID-19 and can improve the efficacy of standard treatments while diminishing disease deterioration [34, 35]. For example, Taiwan Chingguan Yihau (NRICM101) has been shown to disrupt virus invasion and host inflammation in patients with SARS-CoV2 infections [36, 37]. Our data showed that treatment with A. argyi, a well-known herbal medicine used in the Far East, can be a new strategy to prevent infections with multiple variants of SARS-CoV2 by suppressing their cellular entry as well as viral replication via targeting cellular proteins TMPRSS2 and ACE2 and viral protein RdRp, respectively (Fig. 7 and Additional file 6: Fig. S6D).
The binding affinity of jaceosidin and eupatilin, two phytochemicals in A. argyi, to SARS-CoV-2 3CLpro has been predicted in molecular docking simulation . In the present study, however, several flavonoids, including 5,7,3'-Trihydroxy-6,4',5'-trimethoxyflavone, β-Rhamnocitrin, eriodictyol, eupatilin, hispidulin, and jaceosidin containing the basic skeleton of bezno-γ-pyrene  did not show the inhibitory effect on the activity of SARS-CoV-2 3CLpro (Table 1), suggesting that the flavonoids in A. argyi mainly target the cellular entry of SARS-CoV-2 to decrease the deterioration from COVID-19.
Although the cellular entry of SARS-CoV-2. B1.1.7 (Alpha) also depends on the binding to human receptor ACE2, the transmissibility of B1.1.7 is lower than that of other classified SARS-CoV-2 mutants . The S1-receptor binding domain (RBD) of B1.1.7 may maintain an unfavorable conformation before ACE2 engagement, causing its weak binding affinity to human ACE2 [41, 42]. It explains the fewer effects of A. argyi, eriodictyol, and umbelliferone on the cellular entry of SARS-CoV-2. B1.1.7 (Alpha) (Fig. 4B). Similar to our findings, many other S protein-targeting agents, such as Echinacea purpurea  and Montelukast , showed efficacy in preventing infections with different SARS-CoV-2 mutants except B1.1.7.
SARS-CoV-2 lineage B.1.1.529, named Omicron strain, was documented as a VOC on 26th November 2021 by WHO and became the significant viral strain in the COVID-19 pandemic. Unlike other SARS-CoV-2 variants, Omicron is marked by many mutations across the entire genome, including its spike glycoprotein gene , contributing to immune escape, an attenuating ability of vaccines, and a higher frequency of reinfections [5, 6]. Unlike the alpha strain, these high mutations of the Omicron S protein are associated with its increased binding efficacy to host receptor ACE2 . Our data also showed that infections with Omicron variants also caused a higher frequency of inflammation in lung tissues (Fig. 6 and Additional files 3, 4, 5: Figs. S3, S4, S5). While other SARS-CoV-2 variants require the S protein priming by the host transmembrane protein TMPRSS2 for cellular entry, the infection of Omicron variants can be accomplished by an endocytic route in a TMPRSS2-independent manner [47, 48]. Our results found that both A. argyi and umbelliferone diminished the protein level of ACE2 (Fig. 4), accounting for their activity in blocking the cellular entry of Omicron variants (Figs. 5, 6 and Additional file 4: Fig. S4). Therefore, umbelliferone would be the critical phytochemical in A. argyi capable of reducing the infection with Omicron variants.
In general, this study explored the potent antiviral activity of A. argyi and its phytochemicals (eriodictyol and umbelliferone) against multiple SARS-CoV-2 variants with the IC50 values ranging from 332 to 1534 µg/ml and 67 to 126 µM, respectively. Mechanistically, A. argyi, eriodictyol, and umbelliferone repress the enzymatic activity of TMPRSS2 and furin for priming S protein, impeding the interaction between S protein and ACE2, and inhibit the RdRp-mediated viral replication (Fig. 7 and Additional file 6: Fig. S6D). In addition, A. argyi and umbelliferone also showed specific activity against Omicron variants due to their inhibitory effect on ACE2 protein expression.
Cell lines culture
Beas 2B epithelial cells, which were established from normal human bronchus, were grown in Dulbecco’s Modified Eagle Medium (DMEM) with low glucose (1 g/L) and sodium pyruvate (Gibco), and ACE2-expressing HEK-293 T cells were cultured in DMEM/Nutrient Mixture (F-12) medium (Gibco). The culture mediums for these two cell lines contained 10% fetal bovine serum (FBS; Gibco) and 1% penicillin–streptomycin (P/S; Thermo Fisher Scientific), and both cell lines were incubated at 37 ℃ in a humidified 5% CO2/95% air incubator.
Preparation of A. argyi extraction and pure compounds
The powder of A. argyi obtained from dry plants was dissolved and heated in sterilized H2O at 37 ℃ for 30 min. Eriodictyol (Cayman), umbelliferone (Sigma-Aldrich), and other ingredients of A. argyi listed in Table 1 were dissolved in DMSO.
Measurement of cell viability
The cell viability was measured in MTT (Sigma-Aldrich) and cell counting kit-8 (CCK-8; Dojindo) according to the manufacturer’s protocols. Beas 2B cells (5000 cells/well) in a 96-well plate were treated with A. argyi solution, eriodictyol, or umbelliferone in a dose-dependent manner for 2 days and then subjected to MTT assays. ACE2-expressing HEK-293 T cells were treated with A. argyi solution, eriodictyol, or umbelliferone for 2 days followed by infections with the variants of SARS-CoV-2 pseudoviruses, and their viabilities were then detected in CCK-8 assays.
Infection with pseudo-typed lentiviral particles (Vpp) of SARS-CoV-2 S protein mutants
The virus particle pseudo-typed (Vpp) of SARS-CoV-2 S protein mutants and the vesicular stomatitis-G (VSV-G) control bearing luciferase were obtained from RNA Technology Platform and Gene Manipulation Core, Academia Sinica in Taiwan. ACE2-expressing HEK-293 T cells were pre-treated with A. argyi, eriodictyol, and umbelliferone at the indicated concentrations for 2 days and then were infected with the pseudo lentivirus of VSVG and SARS-CoV-2 variants for 1 day. After determining the cell viability in cell counting kit-8 assay (CCK-8; Dojindo), the luciferase intensity of the lysates prepared with One-Glo™ Luciferase assay buffer (Promega) was measured in Luminescence Plate Reader.
Molecular docking simulation
The structures of proteins and compounds in this study were retrieved from Protein Data Bank (PDB, https://www.rcsb.org/) and were applied to molecular docking calculation (BIOVIA Discovery Studio) to simulate the binding efficacy of tested compounds to S protein and receptor ACE2 as well as TMPRSS2.
Fluorescence resonance energy transfer (FRET)-based enzymatic activity assay
In the examination of TMPRSS2 and Furin protease activities in the cleavage of SARS-CoV-2 spike protein, the recombinant proteins (15 μg/ml) of TMPRSS2 and Furin were pre-incubated in the assay buffer (25 mM Tris 8.0, 150 mM NaCl) with/without eriodictyol at room temperature for 30 min. Additionally, Camostat (CAM)  and naphthofluorescein (NAP)  were employed as the positive control against the activity of TMPRSS2 and furin, respectively. After adding 20 μM of the fluorescent protein substrates, the reaction of substrate cleavage was monitored continuously for 6 h by detecting mNeonGreen fluorescence (excitation: 506 nm/emission: 536 nm) using Synergy™ H1 hybrid multi-mode microplate reader (BioTek Instruments, Inc.) [23, 51]. The first 1 h of the reaction was used to calculate the initial velocity (V0). The initial velocity with each compound was calculated and normalized to DMSO control (as illustrated in Fig. 1A).
In the analysis of the enzymatic activities associated with SARS-CoV-2 replication, including 3C-like protease (3CLpro), papain-like protease (PLpro), and RNA-dependent RNA polymerase (RdRp), the tested compounds were first incubated with the recombinant proteins of these proteases at room temperature. Additionally, GC-376 (GC) , GRL-0617 (GRL) , and GS-443902 (GS, an active metabolite from remdesivir) , and cepharanthine , were employed as positive controls against the activity of 3CLpro, PLpro, RdRp, and S/ACE2 interaction, respectively. After 30 min of incubation, the fluorescent protein substrate was added to the reaction. The fluorescent signal was recorded by detecting the excitation wavelength of 505 nm and emission wavelength of 536 nm .
Time-resolved fluorescence resonance energy transfer (TR-FRET) assay
The interruption effects of tested compounds on the interaction between SARS-CoV-2 Spike S1 and human ACE2 proteins were measured by using TR-FRET assays according to the manufacturer’s protocol (BPS Bioscience, Inc.) [23, 51]. Briefly, the recombinant proteins of ACE2 and SARS-CoV-2 Spike S1 were incubated with or without the tested compounds at the indicated concentrations at room temperature for 1 h. TR-FRET signals were recorded by detecting the emission at a wavelength of 620 or 665 nm with the excitation at a wavelength of 340 nm (Fig. 1D).
In vivo infection
The 8-week-old female BALB/c mice (5 for each group) were purchased from LASCO. For pseudovirus infection, mice were first treated with 2.5 × 108 PFU of Adv-hACE2 via intranasal administration. After Adv-hACE2 infections for 7 days, mice were intranasally delivered 105 PFU of VSVG or SARS-CoV-2 pseudo-lentiviruses of wild-type spike protein, BA.2, and XBB.1 for 7 days . Before 2 days of infections with SARS-CoV-2 pseudo-lentiviruses, mice were orally fed with water, 70 mg/kg umbelliferone, or 10 mg/kg cepharanthine for 9 days (as illustrated in Fig. 5A). At the indicated harvest time point, mice were euthanized by CO2 inhalation, and lungs were collected for analyzing histopathology, TMPRSS2 expression, and the mRNA levels of inflammation cytokines in immunohistochemistry, western blotting, and RT-qPCR analyses, respectively. All the procedures of the animal experiments in this study were approved by the Institutional Animal Care and Use Committee of China Medical University (CMUIACUC-2021-87) in accordance with NIH guidelines.
Western blot analysis
Protein lysates from mouse lungs or Beas 2B lung cells prepared in RIPA buffer with protease and phosphatase inhibitors were separated in SDS-PAGE and then transferred to PVDF membranes. The membranes were blocked in 5% milk in TBST buffer (TBS with 0.1% Tween 20) and were incubated with primary antibodies against TMPRSS2 (ProteinTech), ACE2 (Genetex), or β-actin (Sigma-Aldrich) at 4 ℃ for overnight followed by the incubation with HRP-conjugated second antibody at room temperature for 1 h. After washing with TBST buffer, the immunoreactive signals were visualized using enhanced chemiluminescence with enhanced chemiluminescence (ECL) reagent.
Quantitative reverse transcription PCR (RT-qPCR)
Total RNA from mouse lungs or Beas 2B lung cells was extracted with TRIzol reagent (Invitrogen). The first-strand cDNA synthesis was started in a 20 μl mixture with 500 ng of mRNA, dNTP, Oligo-dT, and distilled DEPC H2O at 65 ℃ for 5 min, and was added with 5X First-Strand Buffer, 0.5 M DTT, 200 U/μl M-MLV reverse transcriptase (Invitrogen) for the incubation at 37 ℃, 25 ℃, and 37 ℃ for 2, 10, 50 min, respectively, followed by another incubation at 70 ℃ for 15 min. The cDNA products were subjected to real-time PCR analyses with SYBR Green Master Mix to detect the gene expressions of ACE2, TMPRSS2, inflammatory cytokines, and α-SMA with specific primer sets (Additional file 7: Table S1).
Immunohistochemistry (IHC) staining
The sections of paraffin-embedded lung tissue (5 μm) were dewaxed by xylene and rehydrated by 100%, 95%, 80%, and 75% of ethanol. These sections were incubated overnight with indicated primary antibodies and then were washed to remove unbound primary antibodies. These sections are then stained with a polymer HRP-conjugated secondary rabbit antibody for 30 min followed by a reaction with diaminobenzidine for 2 min. TissueFAXS was used to capture the images and perform the quantitative analysis.
Data are shown as the mean ± standard error of the mean (SEM). A two-tailed t-test was used for comparisons. A p-value < 0.05 was considered statistically significant.
Availability of data and materials
All of performed and analyzed results in this main article are available from the corresponding authors upon request.
- A.argyi :
Traditional Chinese medicine
Severe acute respiratory syndrome coronavirus 2
Proteins transmembrane serine protease 2
Angiotensin-converting enzyme 2
- S protein:
Coronavirus disease 2019
Variant of interests
Variant of concerns
World Health Organization
Middle east respiratory syndrome
Peroxisome proliferation-activated receptor gamma
Cytotoxicity concentration 50%
- SARS-CoV-2 S-Vpp:
SARS-CoV-2 S protein-pseudotyped lentiviral particles
Inhibitory concentration 50%
- 3CLpro :
- PLpro :
RNA-dependent RNA polymeras
Alpha-smooth muscle actin
Receptor binding domain
Fluorescence resonance energy transfer
Time-resolved fluorescence resonance energy transfer
Del Rio C, Malani PN. COVID-19-New insights on a rapidly changing epidemic. JAMA. 2020;323(14):1339–40.
Harvey WT, Carabelli AM, Jackson B, Gupta RK, Thomson EC, Harrison EM, Ludden C, Reeve R, Rambaut A, Consortium C-GU, et al. SARS-CoV-2 variants spike mutations and immune escape. Nat Rev Microbiol. 2021;19(7):409–24.
Coronaviridae Study Group of the International Committee on Taxonomy of V. The species severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nat Microbiol. 2020;5(4):536–44.
Karim SSA, Karim QA. Omicron SARS-CoV-2 variant: a new chapter in the COVID-19 pandemic. Lancet. 2021;398(10317):2126–8.
Pulliam JRC, van Schalkwyk C, Govender N, von Gottberg A, Cohen C, Groome MJ, Dushoff J, Mlisana K, Moultrie H. Increased risk of SARS-CoV-2 reinfection associated with emergence of omicron in South Africa. Science. 2022;376(6593):eabn4947.
Willett BJ, Grove J, MacLean OA, Wilkie C, De Lorenzo G, Furnon W, Cantoni D, Scott S, Logan N, Ashraf S, et al. SARS-CoV-2 omicron is an immune escape variant with an altered cell entry pathway. Nat Microbiol. 2022;7(8):1161–79.
Gao Y, Yan L, Huang Y, Liu F, Zhao Y, Cao L, Wang T, Sun Q, Ming Z, Zhang L, et al. Structure of the RNA-dependent RNA polymerase from COVID-19 virus. Science. 2020;368(6492):779–82.
Wang Q, Zhang Y, Wu L, Niu S, Song C, Zhang Z, Lu G, Qiao C, Hu Y, Yuen KY, et al. Structural and functional basis of SARS-CoV-2 entry by using human ACE2. Cell. 2020;181(4):894-904.e899.
Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S, Schiergens TS, Herrler G, Wu NH, 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(2):271-280.e278.
Wrapp D, Wang N, Corbett KS, Goldsmith JA, Hsieh CL, Abiona O, Graham BS, McLellan JS. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science. 2020;367(6483):1260–3.
Bayati A, Kumar R, Francis V, McPherson PS. SARS-CoV-2 infects cells after viral entry via clathrin-mediated endocytosis. J Biol Chem. 2021;296:100306.
Hoffmann M, Kleine-Weber H, Pöhlmann S. A multibasic cleavage site in the spike protein of SARS-CoV-2 Is essential for infection of human lung cells. Mol Cell. 2020;78(4):779-784.e775.
Shang J, Wan Y, Luo C, Ye G, Geng Q, Auerbach A, Li F. Cell entry mechanisms of SARS-CoV-2. Proc Natl Acad Sci U S A. 2020;117(21):11727–34.
Dömling A, Gao L. Chemistry and biology of SARS-CoV-2. Chem. 2020;6(6):1283–95.
V’Kovski P, Kratzel A, Steiner S, Stalder H, Thiel V. Coronavirus biology and replication: implications for SARS-CoV-2. Nat Rev Microbiol. 2021;19(3):155–70.
Luo DY, Yan ZT, Che LR, Zhu JJ, Chen B. Repellency and insecticidal activity of seven Mugwort (Artemisia argyi) essential oils against the malaria vector Anopheles sinensis. Sci Rep. 2022;12(1):5337.
Bisht D, Kumar D, Kumar D, Dua K, Chellappan DK. Phytochemistry and pharmacological activity of the genus Artemisia. Arch Pharm Res. 2021;44(5):439–74.
Nigam M, Atanassova M, Mishra AP, Pezzani R, Devkota HP, Plygun S, Salehi B, Setzer WN, Sharifi-Rad J. Bioactive compounds and health benefits of Artemisia Species. Nat Prod Commun. 2019;14(7):1934578X1985035.
Lang SJ, Schmiech M, Hafner S, Paetz C, Steinborn C, Huber R, Gaafary ME, Werner K, Schmidt CQ, Syrovets T, et al. Antitumor activity of an Artemisia annua herbal preparation and identification of active ingredients. Phytomedicine. 2019;62:152962.
Tin AS, Sundar SN, Tran KQ, Park AH, Poindexter KM, Firestone GL. Antiproliferative effects of artemisinin on human breast cancer cells requires the downregulated expression of the E2F1 transcription factor and loss of E2F1-target cell cycle genes. Anticancer Drugs. 2012;23(4):370–9.
Ren W, Liang P, Ma Y, Sun Q, Pu Q, Dong L, Luo G, Mazhar M, Liu J, Wang R, et al. Research progress of traditional Chinese medicine against COVID-19. Biomed Pharmacother. 2021;137:111310.
Zhao Z, Li Y, Zhou L, Zhou X, Xie B, Zhang W, Sun J. Prevention and treatment of COVID-19 using traditional Chinese Medicine: a review. Phytomedicine. 2021;85:153308.
Cheng FJ, Huynh TK, Yang CS, Hu DW, Shen YC, Tu CY, Wu YC, Tang CH, Huang WC, Chen Y, et al. Hesperidin is a potential inhibitor against SARS-CoV-2 infection. Nutrients. 2021;13(8):2800.
Earles SM, Bronstein JC, Winner DL, Bull AW. Metabolism of oxidized linoleic acid: characterization of 13-hydroxyoctadecadienoic acid dehydrogenase activity from rat colonic tissue. Biochim Biophys Acta. 1991;1081(2):174–80.
Reinaud O, Delaforge M, Boucher JL, Rocchiccioli F, Mansuy D. Oxidative metabolism of linoleic acid by human leukocytes. Biochem Biophys Res Commun. 1989;161(2):883–91.
Altmann R, Hausmann M, Spöttl T, Gruber M, Bull AW, Menzel K, Vogl D, Herfarth H, Schölmerich J, Falk W, et al. 13-Oxo-ODE is an endogenous ligand for PPARgamma in human colonic epithelial cells. Biochem Pharmacol. 2007;74(4):612–22.
Dubuquoy L, Rousseaux C, Thuru X, Peyrin-Biroulet L, Romano O, Chavatte P, Chamaillard M, Desreumaux P. PPARgamma as a new therapeutic target in inflammatory bowel diseases. Gut. 2006;55(9):1341–9.
Walls AC, Park YJ, Tortorici MA, Wall A, McGuire AT, Veesler D. Structure, function, and antigenicity of the SARS-CoV-2 Spike glycoprotein. Cell. 2020;181(2):281-292.e286.
Guo L, Bai S, Ding S, Zhao L, Xu S, Wang X. PM2.5 Exposure induces lung injury and fibrosis by regulating ferroptosis via TGF-beta signaling. Dis Markers. 2022;2022:7098463.
Báez-Santos YM, St John SE, Mesecar AD. The SARS-coronavirus papain-like protease: structure, function and inhibition by designed antiviral compounds. Antiviral Res. 2015;115:21–38.
Grum-Tokars V, Ratia K, Begaye A, Baker SC, Mesecar AD. Evaluating the 3C-like protease activity of SARS-coronavirus: recommendations for standardized assays for drug discovery. Virus Res. 2008;133(1):63–73.
Kirchdoerfer RN, Ward AB. Structure of the SARS-CoV nsp12 polymerase bound to nsp7 and nsp8 co-factors. Nat Commun. 2019;10(1):2342.
Krause PR, Fleming TR, Longini IM, Peto R, Briand S, Heymann DL, Beral V, Snape MD, Rees H, Ropero AM, et al. SARS-CoV-2 variants and vaccines. N Engl J Med. 2021;385(2):179–86.
Alam S, Sarker MMR, Afrin S, Richi FT, Zhao C, Zhou JR, Mohamed IN. Traditional herbal medicines, bioactive metabolites, and plant products against COVID-19: update on clinical trials and mechanism of actions. Front Pharmacol. 2021;12:671498.
Wang H, Xu B, Zhang Y, Duan Y, Gao R, He H, Li X, Li J. Efficacy and safety of traditional Chinese medicine in coronavirus disease 2019 (COVID-19): a systematic review and meta-analysis. Front Pharmacol. 2021;12:609213.
Tsai KC, Huang YC, Liaw CC, Tsai CI, Chiou CT, Lin CJ, Wei WC, Lin SJ, Tseng YH, Yeh KM, et al. A traditional Chinese medicine formula NRICM101 to target COVID-19 through multiple pathways: a bedside-to-bench study. Biomed Pharmacother. 2021;133:111037.
Tseng YH, Lin SJ, Hou SM, Wang CH, Cheng SP, Tseng KY, Lee MY, Lee SM, Huang YC, Lin CJ, et al. Curbing COVID-19 progression and mortality with traditional Chinese medicine among hospitalized patients with COVID-19: A propensity score-matched analysis. Pharmacol Res. 2022;184:106412.
Lu Y, Zhang B, Wang N, Li M, Xi N. Investigation of major flavonoids from Artemisia argyi as a potential COVID-19 drug: molecular docking and DFT calculations. Crystals. 2022;12(7):990.
Song X, Wen X, He J, Zhao H, Li S, Wang M. Phytochemical components and biological activities of Artemisia argyi. Journal of Functional Foods. 2019;52:648–62.
Cojocaru C, Cojocaru E, Turcanu AM, Zaharia DC. Clinical challenges of SARS-CoV-2 variants (Review). Exp Ther Med. 2022;23(6):416.
Socher E, Conrad M, Heger L, Paulsen F, Sticht H, Zunke F, Arnold P. Mutations in the B117 SARS-CoV-2 spike protein reduce receptor-binding affinity and induce a flexible link to the fusion peptide. Biomedicines. 2021;9(5):525.
Niu Z, Zhang Z, Gao X, Du P, Lu J, Yan B, Wang C, Zheng Y, Huang H, Sun Q. N501Y mutation imparts cross-species transmission of SARS-CoV-2 to mice by enhancing receptor binding. Signal Transduct Target Ther. 2021;6(1):284.
Vimalanathan S, Shehata M, Sadasivam K, Delbue S, Dolci M, Pariani E, D’Alessandro S, Pleschka S. Broad antiviral effects of Echinacea purpurea against SARS-CoV-2 variants of concern and potential mechanism of action. Microorganisms. 2022;10(11):2145.
Luedemann M, Stadler D, Cheng CC, Protzer U, Knolle PA, Donakonda S. Montelukast is a dual-purpose inhibitor of SARS-CoV-2 infection and virus-induced IL-6 expression identified by structure-based drug repurposing. Comput Struct Biotechnol J. 2022;20:799–811.
Viana R, Moyo S, Amoako DG, Tegally H, Scheepers C, Althaus CL, Anyaneji UJ, Bester PA, Boni MF, Chand M, et al. Rapid epidemic expansion of the SARS-CoV-2 omicron variant in southern Africa. Nature. 2022;603(7902):679–86.
Mannar D, Saville JW, Zhu X, Srivastava SS, Berezuk AM, Tuttle KS, Marquez AC, Sekirov I, Subramaniam S. SARS-CoV-2 omicron variant: antibody evasion and cryo-EM structure of spike protein-ACE2 complex. Science. 2022;375(6582):760–4.
Peacock TP, Brown JC, Zhou J, Thakur N, Newman J, Kugathasan R, Sukhova K, Kaforou M, Bailey D, Barclay WS. The SARS-CoV-2 variant, Omicron, shows rapid replication in human primary nasal epithelial cultures and efficiently uses the endosomal route of entry. bioRxiv doi: 10.1101/2021.12.31.474653
Zhao H, Lu L, Peng Z, Chen LL, Meng X, Zhang C, Ip JD, Chan WM, Chu AW, Chan KH, et al. SARS-CoV-2 Omicron variant shows less efficient replication and fusion activity when compared with delta variant in TMPRSS2-expressed cells. Emerg Microbes Infect. 2022;11(1):277–83.
Hu Y, Ma C, Szeto T, Hurst B, Tarbet B, Wang J. Boceprevir, calpain Inhibitors II and XII, and GC-376 have broad-spectrum antiviral activity against coronaviruses. ACS Infect Dis. 2021;7(3):586–97.
Cheng YW, Chao TL, Li CL, Chiu MF, Kao HC, Wang SH, Pang YH, Lin CH, Tsai YM, Lee WH, et al. Furin inhibitors block SARS-CoV-2 spike protein cleavage to suppress virus production and cytopathic effects. Cell Rep. 2020;33(2):108254.
Chen Y, Huang WC, Yang CS, Cheng FJ, Chiu YF, Chen HF, Huynh TK, Huang CF, Chen CH, Wang HC, et al. Screening strategy of TMPRSS2 inhibitors by FRET-based enzymatic activity for TMPRSS2-based cancer and COVID-19 treatment. Am J Cancer Res. 2021;11(3):827–36.
Fu Z, Huang B, Tang J, Liu S, Liu M, Ye Y, Liu Z, Xiong Y, Zhu W, Cao D, et al. The complex structure of GRL0617 and SARS-CoV-2 PLpro reveals a hot spot for antiviral drug discovery. Nat Commun. 2021;12(1):488.
Cox RM, Wolf JD, Lieber CM, Sourimant J, Lin MJ, Babusis D, DuPont V, Chan J, Barrett KT, Lye D, et al. Oral prodrug of remdesivir parent GS-441524 is efficacious against SARS-CoV-2 in ferrets. Nat Commun. 2021;12(1):6415.
Fan H, He ST, Han P, Hong B, Liu K, Li M, Wang S, Tong Y. Cepharanthine: a promising old drug against SARS-CoV-2. Adv Biol. 2022;6(12):e2200148.
Chen Y, Yang WH, Chen HF, Huang LM, Gao JY, Lin CW, Wang YC, Yang CS, Liu YL, Hou MH, et al. Tafenoquine and its derivatives as inhibitors for the severe acute respiratory syndrome coronavirus 2. J Biol Chem. 2022;298(3):101658.
Tseng SH, Lam B, Kung YJ, Lin J, Liu L, Tsai YC, Ferrall L, Roden RBS, Wu TC, Hung CF. A novel pseudovirus-based mouse model of SARS-CoV-2 infection to test COVID-19 interventions. J Biomed Sci. 2021;28(1):34.
Parts of experiments and data analysis were performed with the use of the Medical Research Core Facilities, Office of Research & Development at China Medical University, Taichung, Taiwan, ROC.
This study was supported by grants from the Ministry of Science Technology, Taiwan (MOST 110-2320-B-039-056-MY3, MOST 111-2320-B-039-016-MY3, MOST 111-2320-B-039-027-MY3, MOST 111-2811-B-039-028, and MOST 111-2320-B-039-036) and China Medical University Hospital (CMUH-DMR-111-025, CMUH-DMR-110-026, CMUH-DMR-110-028, and CMUHCH-DMR-110-006). This work was also supported by the Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Ministry of Education (MOE) in Taiwan.
Ethics approval and consent to participate
Consent for publication
The authors have declared no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The cytotoxicity of A. argyi, eriodictyol, and umbelliferone is illustrated in lung epithelial cells. Beas 2B cells were treated with A. argyi A eriodictyol B, and umbelliferone C at the indicated concentrations for 2 days. Then the cell viability was detected using MTT assays. The cytotoxicity concentration 50%value was calculated and displayed. Data was shown as mean±SEM from three independent experiments with triplicates.
A. argyi, eriodictyol, and umbelliferone did not affect the cellular entry of VSVG-pseudotyped lentivirus.HEK-293T cells expressing ACE2 were infected with VSVG-pseudotyped lentivirus encoding luciferase after treatments with A. argyi, eriodictyol, or umbelliferone at the indicated concentrations for 2 days, and the luciferase activities were subsequently measured to determine the infection rate with the pseudoviruses. Additionally, the cell viability also was examined in CCK-8 assays. Camostat, naphthofluorescein, and cepharanthinewere used as the positive controls for the inhibitions of TMPRSS2 and furin activities and the binding between S protein and ACE2, respectively. Data were shown as mean±SEM from three independent experiments with triplicates. *** p<0.001.
Umbelliferone represses the wild-type SARS-CoV-2-associated pulmonary inflammation in vivo.The mice infected with VSVG or wild-type SARS-CoV-2 pseudoviruses were sacrificed after oral administration with sterilized water, 70 mg/kg umbelliferone, and 10 mg/kg cepharanthinefor 9 days, and lungs were subsequently harvested for the examinations of histopathological changes in H&E staining A lymphocyte infiltration and TNF-α protein expression in IHC staining A TMPRSS2 protein in Western blot analysis B and the mRNA levels of indicated genes in RT-qPCR analysis C–J respectively. * p < 0.05, ** p < 0.01, *** p <0.001, and **** p < 0.0001.
Umbelliferone represses Omicron BA.2-associated pulmonary inflammation in vivo.The mice infected with VSVG or Omicron BA.2 pseudoviruses were sacrificed after oral administration with sterilized water, 70 mg/kg umbelliferone, and 10 mg/kg cepharanthinefor 9 days, and lungs were subsequently harvested for the examinations of histopathological changes in H&E staining A lymphocyte infiltration and TNF-α protein expression in IHC staining A TMPRSS2 protein in Western blot analysis B and the mRNA levels of indicated genes in RT-qPCR analysis C–J respectively. * p < 0.05, ** p < 0.01, and *** p <0.001. The granuloma formation is marked with a circle.
Omicron variants strongly elicited the expressions of inflammatory cytokines and fibrosis markers in vivo.. The mice infected with wild-type SARS-CoV-2 and Omicron variantspseudo lentiviruses were sacrificed, and total RNA extracted from lung tissues were subjected to RT-qPCR analyses for the mRNA levels of inflammatory cytokines, including Tnf A Ifng B Il1b C Il6 D Il10 E and Tgfb1 F and lung fibrosis marker Acta2 G. * p < 0.05, ** p < 0.01, and **** p < 0.0001.
Eriodictyol suppresses the activity of viral replication enzymes. A FRET-based enzymatic activity assay was performed to address the activity of RdRp A 3CLpro B and PLpro C under treatment with eriodictyol at the indicated concentrations. GS-443902, GC-376, and GRL-0617were used as the positive controls for the inhibitions of RdRp, 3CLpro, and PLpro activities, respectively. Data were shown as mean±SEM from three independent experiments with triplicates. ** p < 0.01 and **** p < 0.0001. D Schematics of the working model of eriodictyol in inhibiting viral replication by antagonizing RdRp activity. The working model was created by bioRENDER.
The primers for determining gene expressions in this study.
About this article
Cite this article
Cheng, FJ., Ho, CY., Li, TS. et al. Umbelliferone and eriodictyol suppress the cellular entry of SARS-CoV-2. Cell Biosci 13, 118 (2023). https://doi.org/10.1186/s13578-023-01070-y