Asmaa I. Owis*ab,
Marwa S. El-Hawaryac,
Dalia El Amira,
Hesham Refaatd,
Eman Alaaeldinde,
Omar M. Alyf,
Mahmoud A. Elrehanygh and
Mohamed S. Kamelci
aDepartment of Pharmacognosy, Faculty of Pharmacy, Beni-Suef University, Beni-Suef, Egypt
bDepartment of Pharmacognosy, Faculty of Pharmacy, Heliopolis University for Sustainable Development, Cairo, Egypt
cDepartment of Pharmacognosy, Faculty of Pharmacy, Deraya University, Minia, Egypt
dDepartment of Pharmaceutics, Faculty of Pharmacy, Deraya University, Minia, Egypt
eDepartment of Pharmaceutics, Faculty of Pharmacy, Minia University, Minia, Egypt
fDepartment of Medicinal Chemistry, Faculty of Pharmacy, Minia University, Minia, Egypt
gDepartment of Biochemistry and Molecular Biology, Faculty of Pharmacy, Deraya University, Minia, Egypt
hDepartment of Biochemistry and Molecular Biology, Faculty of Medicine, Minia University, Minia, Egypt
iDepartment of Pharmacognosy, Faculty of Pharmacy, Minia University, Minia, Egypt
First published on 12th April 2021
Several studies are now underway as a worldwide response for the containment of the COVID-19 outbreak; unfortunately, none of them have resulted in an effective treatment. Salvadora persica L. (Salvadoraceae), commonly known as meswak, is one of the popular plants used by Muslims as an oral hygiene tool. It is documented that the meswak possesses antiviral activity, but no report discusses its use for coronavirus treatment. Herein, a mixture of 11 flavonoids prepared from the aqueous plant extract and its liposomal formulation were shown to inhibit SARS-CoV-2 in an in vitro A549 cell line culture and a RT-PCR test almost as well as the FDA-approved anti-COVID-19 agent, remdesivir. Encapsulation within liposomal formulation led to a highly significant increase in the percentage of inhibition of viral replication from 38.09 ± 0.83 to 85.56 ± 1.12% in a flavonoid mixture and its liposomal preparation, respectively, and this figure approached that obtained for remdesivir (91.20 ± 1.71%). Preliminary tests were also performed, including a total flavonoid assay, a molecular docking study, a 3CL-protease inhibition assay and a cytotoxicity study. It was worthy to find a cheap, readily available, safe natural source for promising anti-SARS-CoV-2 agents, that leak their phytochemicals into the aqueous saliva during regular use as a brushing agent.
A few protocols are being used to treat laboratory confirmed or suspected COVID-19 patients hospitalized with severe disease, for example, chloroquine derivatives,11 convalescent plasma,12 azithromycin13 and remdesivir, with limited known information about its safety on human cells.14 Natural products are considered as viable lead candidates that have a highly significant role in drug discovery. Salvadora persica L. (Salvadoraceae) or chewing stick, commonly known in Arabic as meswak, is one of the popular plants used by Muslims as an oral hygiene tool. Several scientific studies have worked on the pharmacological profile of S. persica, with a wealth of published literature now available. Besides the effects on oral health, its biological activities have been found to include antimicrobial, hypoglycemic, antioxidant, anticonvulsant, antiulcer, sedative, analgesic, antiosteoporosis, anti-inflammatory, hypolipidemic, antitumor and antiviral effects,15,16 but there are no reports dealing with coronavirus. A literature survey also revealed the presence of a variety of phytochemicals in S. persica extracts including alkaloids, tannins, glycosides, saponins, essential oils and flavonoids.15 Previous research on an aqueous extract of S. persica led to the isolation and identification of 11 flavonoid glycosides.17 Interestingly, 10 of them exhibited significant binding stability at the N3 binding site, found through a molecular docking study against the COVID-19 main protease (Mpro), suggesting their potential effect on blocking viral replication.14 Therefore, the objectives of this study were to prepare a flavonoid-rich fraction (FRF) from an aqueous extract of the aerial parts of S. persica, and investigate its anti-COVID-19 activity using the preliminary steps as a molecular docking study on the contact surface of hACE2-COVID 19 spike protein complex, a 3CL-protease inhibition assay, and cytotoxic activity. The final objective was to prepare a liposomal formulation (FRF-Lip) from this fraction to increase its viral inhibitory action and assess the potential ability of both fractions to affect human SARS-CoV-2 using an in vitro assay with remdesivir as reference. As far as is known, this is the first practical study dealing with the effect of S. persica on COVID-19.
Fig. 2 The alpha sphere (A) of the main binding pocket between hACE2 (B) and the COVID-19 spike protein (C). |
The interaction energies of compounds 1–11 revealed that most of the compounds showed a prominent interaction with both hACE2 and the viral spike protein (Table 1, Fig. 3 and 4). The most promising compound, with the highest binding energy (−9.4799 kcal per mole), was compound 2. Compounds 5, 1, 3, 4 and 7, with binding energies of −7.6404, −7.3869, −7.3769 and −7.2929 kcal per mole, respectively, showed higher binding energies than hesperidin. Whereas compounds 9, 8, 10, 11 and 6 showed moderate binding energies which were less than that of hesperidin (−6.6629 kcal per kmole). The six compounds (1–6) with the highest binding energies were able to bind to the spike protein with similar hydrogen bonding with Asp 405 or Asp 406, and Lys 403. Compounds 1, 2 and 5 were like hesperidin and were able to bind to hACE2 with hydrogen bonding with Glu 37. Compounds 1, 2, 6 and 7 had similar binding with His 34, either by hydrogen bonding or H–pi bonding. It is worth noting that all the compounds showed many other types of hydrogen bonding in addition to pi–pi and other hydrophobic bonds with the spike protein and hACE2. Interestingly, compound 2 could be used as a starting point for developing more potent derivatives with better therapeutic properties. Flavonols are the most studied, and documented, compounds tested against SARS-CoV-2.24. Pandey et al. reported that kaempferol and isorhamnetin interacted with the S2 domain of the S protein with high binding affinities when compared to hydroxychloroquine (positive control).25 Mouffouk et al. proved that the presence of ortho di-hydroxyl groups in kaempferol created stronger hydrogen bonds and a stable hACE2-COVID 19 spike protein complex than those of isorhamnetin.24 These data were consistent with the results obtained in the research reported here.
Compound no. | dG kcal per mole | Receptor | |
---|---|---|---|
Amino acid/type of bonding/distance (Å)/binding energy (kcal per mole) | |||
hACE2 | Spike protein | ||
1 | −7.3869 | GLU 37/H–donor/3.06/−0.8 | ASP 405/H–donor/3.07/−1.3 |
HIS 34/H–acceptor/2.96/−1.9 | LYS 403/H–acceptor/2.97/−5.2 | ||
ARG 408/pi–cation/4.08/−1.3 | |||
ARG 408/pi–cation/4.06/−0.6 | |||
2 | −9.4799 | GLU 37/H–donor/2.68/−3.1 | ASP 405/H–donor/3.01/−0.8 |
ASP 30/H–donor/2.91/−3.8 | GLY 496/H–acceptor/2.76/−1.3 | ||
HIS 34/H–pi/3.85/−0.9 | LYS 403/H–acceptor/2.90/−5.0 | ||
ALA 387/pi–H/4.39/−0.6 | |||
3 | −7.3769 | ASP 405/H–donor/3.01/−2.3 | |
ASP 406/H–donor/2.76/−4.0 | |||
ASP 406/H–donor/3.11/−1.9 | |||
ASP 406/H–donor/3.26/−1.0 | |||
LYS 403/H–acceptor/2.96/−4.1 | |||
LYS 403/H–acceptor/3.34/−0.9 | |||
4 | −7.2929 | ASP 30/H–donor/3.14/−1.1 | ASP 405/H–donor/2.76/−3.0 |
LYS 403/H–acceptor/3.39/−0.6 | |||
5 | −7.6404 | GLU 37/H–donor/3.08/−1.1 | ASP 406/H–donor/2.67/−1.5 |
GLU 37/H–donor/3.49/−0.8 | LYS 403/H–acceptor/2.76/−7.6 | ||
ARG 393/H–acceptor/3.06/−1.4 | |||
ASN 33/H–acceptor/2.91/−2.7 | |||
6 | −5.2061 | HIS 34/H–donor/2.98/−1.7 | ARG 408/H–acceptor/3.10/−2.5 |
ARG 393/H–acceptor/2.87/−3.3 | |||
HIS 34/pi–pi/3.90/−0.0 | |||
7 | −7.1293 | HIS 34/H–donor/2.78/−1.2 | TYR 453/H–donor/2.64/−1.0 |
ARG 393/H–acceptor/3.15/−1.7 | LYS 403/H–acceptor/3.53/−1.0 | ||
LYS 353/H–acceptor/2.97/−0.8 | ARG 408/pi–cation/3.85/−1.4 | ||
HIS 34/H–pi/4.04/−1.8 | |||
8 | −5.9453 | ASP 30/H–donor/2.99/−3.1 | LYS 403/H–acceptor/2.91/−5.7 |
ASN 33/H–acceptor/2.97/−2.5 | |||
9 | −6.2471 | ASP 30/H–donor/3.38/−0.8 | ASP 406/H–donor/2.94/−3.6 |
LYS 403/H–acceptor/2.90/−3.3 | |||
10 | −5.7829 | ARG 408/H–acceptor/3.38/−0.9 | |
ARG 408/H–acceptor/3.47/−0.8 | |||
11 | −5.3156 | SER 494/H–donor/2.80/−1.0 | |
GLN 409/H–acceptor/2.75/−2.4 | |||
ARG 408/H–acceptor/2.83/−4.1 | |||
Hesperidin | −6.6629 | GLU 37/H–donor/2.75/−2.3 | GLN 409/H–acceptor/2.76 |
Fig. 3 A 3D representation of the docking poses of compound 2 on the COVID-19 spike binding site with hACE-2. |
Fig. 4 A 2D representation of the docking of compound 2 on the COVID-19 spike binding site with hACE-2. |
Studying the physicochemical and ADME properties revealed that all the compounds had a reasonable topological polar surface area (TPSA) between 20 Å2 and 130 Å2 except for the reference compound hesperidin (Table 2). Entries of the Moriguchi logarithm of the partition coefficient (MLOGP) < 4.15 were of reasonable lipophilicity to be taken orally. Compounds 5, 7, 9 and 2 were of reasonable lipophilicity. The predicted gastrointestinal (GI) absorption and blood–brain barrier (BBB) permeability showed that compounds 9, 2, 4, 3 and 1 were expected to be highly absorbed from the GI tract whereas compounds 9, 2, 4, and 1 can pass the BBB (Table 2). It is worth noting that all the tested compounds showed no pan assay interference compounds (PAINS) alerts. Fortunately, the possible deleterious effects on the CNS of COVID-19 infected patients26–31 could be treated or prevented by compounds 2, 1 and 4 as they were the top-scoring compounds with the capability of crossing the BBB.
Compound no. | TPSA Å2 | logPo/w (MLOGP) | GI absorption | BBB permeant | PAINS alert |
---|---|---|---|---|---|
a TPSA: topological polar surface area, MLOGP: Moriguchi logarithm of the partition coefficient, BBB: blood–brain barrier, PAINS: pan assay interference compounds. | |||||
1 | 37.30 | 4.38 | High | Yes | 0 |
2 | 54.37 | 3.59 | High | Yes | 0 |
3 | 37.30 | 4.75 | High | No | 0 |
4 | 37.30 | 4.29 | High | Yes | 0 |
5 | 110.02 | 3.99 | Low | No | 0 |
6 | 69.56 | 5.61 | Low | No | 0 |
7 | 110.02 | 3.88 | Low | No | 0 |
8 | 89.79 | 4.98 | Low | No | 0 |
9 | 37.30 | 3.81 | High | Yes | 0 |
10 | 110.02 | 4.71 | Low | No | 0 |
11 | 110.02 | 4.67 | Low | No | 0 |
Hesperidin | 234.29 | −3.04 | Low | No | 0 |
The virtual protective and Mpro inhibitory role14 of the meswak flavonoids against SARS-CoV-2 was proved using in vitro experiments. Instead of using individual flavonoids, a mixture of 11 flavonoids (FRF) was used in this analysis, based on the synergism capability of natural products, which is thought to have more powerful antiviral activity.32 The FRF were successfully prepared from the aqueous extract of S. persica leaves and stems. The identity of the obtained fractions was confirmed by co-spot-TLC with previously isolated and characterized flavonoids.14,17 Total flavonoid content was determined on an aqueous extract and FRF with contents of 9.8 ± 1.17, and 7.4 ± 1.78 mg Ru per g sample, respectively. This confirmed that the FRFs contained most of the flavonoid content present in the aqueous extract.
According to virtual docking results, the FRFs had an inhibitory effect on both the Mpro enzyme14 and spike protein-ACE-2. The 3CL protease inhibition assay was performed as a preliminary practical evaluation. The FRFs exhibited a significant inhibitory action against 3CL-protease with an IC50 = 8.59 ± 0.3 μg mL−1 whereas that of the standard (tannic acid) was 2.1 ± 0.2 μM.33 This result was consistent with the stabilized binding of all the FRF-compounds (1–11) at the N3-binding site of Mpro as described in a previous docking study.14 The glycosylated flavonols with at least one rhamnose unit in their structure at the carbon 3 position can bind with a higher affinity to the active site of the Mpro than to the corresponding aglycones.34
In an attempt to enhance the potential activity of FRF as an anti-SARS-CoV-2 agent, a nano-sized liposomal formulation (FRF-Lip) was prepared. This hypothesis was based upon the fact that nano-sized particles have a high surface area and possess the ability to introduce both water-soluble and lipophilic components, such as flavonoids, through the cellular membrane into the cells.20 The adopted spraying technique successfully produced a uniformly nanosized FRF-loaded liposomal formulation with a size of 43.7 ± 1.5 nm and a PDI of 0.32. The image obtained with TEM showed the vesicular multilayered structure of the prepared FRF-Lip as shown in Fig. 5 and the % EE y was 47.8 ± 2.1%.
Safety and possible secondary effects were considered as major requirements during the search for a new antiviral lead. The natural FRF obtained from meswak, a plant used regularly by Muslims, was screened for its cellular toxicity effect, and for the determination of the appropriate concentration using the in vitro inhibition assay. The estimated concentration associated with 50% cytotoxicity (CC50) was found to be 24.5 ± 1.9 μg mL−1.
The significant results of the effects of FRFs on 3CL-protease highlighted the need to investigate their behavior, using cell culture, on the human virus and to compare the obtained results with that of the FRF-Lip. The prospective anti-human coronavirus effect of FRF and FRF-Lip was investigated by in vitro RT-PCR testing which specifically detects SARS-CoV-2 RNA quantitatively. Encapsulation of FRF within the liposomal formulation led to a high significant increase in the percentage of inhibition in viral replication from 38.09 ± 0.83% to 85.56 ± 1.12% in FRF and FRF-Lip, respectively (p < 0.0001). This could be attributed to the higher surface area and better dispersibility of the liposomal formulation which was able to prevent precipitation and aggregation of the lipophilic components. Hence, it could optimize the expected poor cellular uptake of the hydrophilic components, and enhance the delivery of such challenging natural products.20 It was useful to find a natural product having the ability to inhibit SARS-CoV-2 replication in a value very close to that of the FDA-approved anti-COVID-19: remdesivir (91.20 ± 1.71%).
The 11 compounds in the FRF are flavonols having a 3,4′-dihydroxy-2-phenylchromen-4-one backbone, namely kaempferol and its 3′-methoxylated derivative: isorhamnetin (Fig. 1). They contain O-sugar groups attached at position 3 of the carbon ring. These 3-O-glycosylated kaempferol analogs were previously documented to have antiviral activity generally and anti-SARS-CoV activity specifically.35
In addition to all the previous results together, the COVID-19 disease is characterized by an excessive inflammatory response leading to a major lung illness and consequent possible mortality. Anti-inflammatory agents are likely to be effective against the subsequent elevated cytokine levels, such as γ-interferon, typically associated with COVID-19 patients.36 The use of S. persica aqueous extract and its flavonoids as anti-inflammatory natural agents37 is not only likely to reduce SARS-CoV-2 viral infectivity, but also likely to reduce host inflammatory response. Moreover, S. persica (meswak) itself, as a natural famous plant used by Muslims, has the advantages of being a cheap, readily available, and safe drug source. Interestingly, human watery saliva can act as an aqueous extractor for its various innate phytochemicals during regular traditional use for oral hygiene by a mechanical cleaning process.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ra00142f |
This journal is © The Royal Society of Chemistry 2021 |