Usama Ramadan Abdelmohsenab,
Amgad Albohyc,
Basma S. Abdulrazikc,
Soad A. L. Bayoumid,
Lourin G. Malakd,
Iman S. A. Khallafd,
Gerhard Bringmann*e and
Salwa F. Farag*df
aDepartment of Pharmacognosy, Faculty of Pharmacy, Deraya University, 7 Universities Zone, 61111 New Minia City, Egypt
bDepartment of Pharmacognosy, Faculty of Pharmacy, Minia University, Minia 61519, Egypt. E-mail: Usama.ramadan@mu.edu.eg
cDepartment of Pharmaceutical Chemistry, Faculty of Pharmacy, The British University in Egypt (BUE), Suez Desert Road, El-Sherouk City, Cairo 11837, Egypt. E-mail: amgad.albohy@bue.edu.eg; basma.sabry@bue.edu.eg
dPharmacognosy Department, Faculty of Pharmacy, Assiut University, Assiut 71526, Egypt. E-mail: soad.bayoumi@pharm.aun.edu.eg; lourinmalak@aun.edu.eg; iman.khallaf@pharm.aun.edu.eg; farag_s@yahoo.com
eInstitute of Organic Chemistry, University of Würzburg, Am Hubland, 97074 Würzburg, Germany. E-mail: bringman@chemie.uni-wuerzburg.de
fPharmacognosy Department, College of Pharmacy, Taif University, Taif, Saudi Arabia
First published on 7th May 2021
COVID-19 is a global pandemic first identified in China, causing severe acute respiratory syndrome. One of the therapeutic strategies for combating viral infections is the search for viral spike proteins as attachment inhibitors among natural compounds using molecular docking. This review aims at shedding light on the antiviral potential of natural products belonging to the natural-products class of coumarins up to 2020. Moreover, all these compounds were filtered based on ADME analysis to determine their physicochemical properties, and the best 74 compounds were selected. Using virtual-screening methods, the selected compounds were investigated for potential inhibition of viral main protease (Mpro), viral methyltransferase (nsp16/10 complex), viral recognition binding domain (RBD) of S-protein, and human angiotensin-converting enzyme 2 (ACE2), which is the human receptor for viral S-protein targets, using molecular-docking studies. Promising potential results against SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) and methyltransferase (nsp16) are presented.
Computer-aided virtual screening of the natural compounds library, natural product activity and species source against type II transmembrane serine protease (TMPRSS2), the priming agent of SARS-CoV-2 and essential for the viral pathogenesis, revealed that compound NPC306344 showed significant interaction with the active-site residues of TMPRSS2, with a binding energy score of −14.69 kcal mol−1.7 Another in silico exploration in marine natural products revealed that some classes of compounds, such as phlorotannins, flavonoids, and pseudopeptides, can inhibit the SARS-CoV-2 virus.8 Epigallocatechin gallate, curcumin, apigenin, and chrysophanol were found to be anti-SARS-CoV-2 active, too.9 Acetylglucopetunidin, isoxanthohumol, and ellagic acid have been suggested as potential drug candidates for the development of anti-SARS-CoV-2 therapeutics.10 Although all these results from molecular docking studies need further in vitro and in vivo research to confirm the anti-SARS-CoV-2 potential, this encouraged us to analyze coumarins as a promising class of phenolic natural products.
Coumarins (2H-1-benzopyran-2-ones) are naturally occurring heterocyclic compounds, which consist of fused benzene and α-pyrone rings. They are derivatives of the lactone of O-hydroxycinnamic acid. The name coumarin is derived from the word coumarou, the French name for the tonka bean (Dipteryx odorata, Fabaceae), from which coumarin was first isolated in 1820. More than 1300 coumarins have meanwhile been isolated from natural sources including plants, bacteria, and fungi. Coumarins were reported in about 150 different plant species belonging to more than 30 different families such as Apiaceae, Rutaceae, Astraceae, Fabaceae, Moraceae, Guttiferae, Thymelaeaceae, Oleaceae, Calophyllaceae, Caprifoliaceae, and Nyctaginaceae.11 Natural coumarins are classified into simple coumarins, furanocoumarins (psoralene and angelicin types), pyranocoumarins (xanthyletin and seselin types), phenylcoumarins, and dicoumarins11,12 (Fig. 1).
Depending on the substitution pattern on the coumarin ring, numerous pharmacological activities of coumarins were reported. These activities are anticoagulant as for dicoumarol, anticonvulsant as for imperatorin and osthole, antioxidant as described for fraxin, antihyperglycemic as for fraxidin, antitubercular as in the case of scopoletin and umbelliferone,13 antifungal as for psoralene and imperatorin, anticancer as reported for bergapten, neuroprotective as in the case of osthole, anti-inflammatory as for esculetin, antibacterial as novobiocin and coumermycin, and antihypertensive as visnadine, while xanthotoxin is effective in the treatment of psoriasis.12,13 Furthermore, antiviral activities of coumarins have been reported.14–30 They include antihepatitis B virus as for nordentatin and clausarin,17 antihepatitis C virus as in the case of wedelolactone, glycycoumarin, and glycyrin,14,15 anti-HIV as for calanolides A, B, and C and for inophyllums A, B, C, E, and P, anti-human influenza virus as for eleutheroside B1,15,17 anti-dengue virus as for ramosin, myrsellinol, and myresellin,17 anti-herpes simplex HSV as in the case of rutamarin,14 anti-poliovirus as for isoobtusitin,31 and anti-chikungunya virus as for coumarins A and B from Mammea americana.15 Recently, in a preliminary screening 29 naturally occurring coumarins as potential SARS-CoV-2 replication inhibitors by docking technique have been investigated, out of which 17 compounds were found to bind to the active site through the interaction with the catalytic dyad, His41 and Cys145, along with other neighboring residues. Among all the investigated compounds, corymbocoumarin, methylgalbanate, and heraclenol, exhibited the best binding efficiency and could be considered as potential Mpro protease inhibitors.16 All these reported data about the antiviral activities of different coumarins encouraged us for investigation of coumarins docking studies against SARS-Cov-2.
Coumarins are characterized by being widely distributed, stable, soluble, low-molecular weight compounds, without adverse side effects and toxicity and with the feasibility of chemical modification to produce new semisynthetic derivatives, and they have been of great interest in the field of medicinal chemistry.17 Coumarins have received considerable attention as most promising candidates for anti-viral drugs, playing a role in targeting various cellular pathways that inhibit the growth and replication of viruses.17
Possible approaches for treating COVID-19 with natural coumarins could be interfering with the viral life cycle stages including binding of the virion with the respective receptor present at the cell surface. In this stage excess soluble forms of angiotensin-converting enzyme 2 (ACE2) are involved, thus the use of ACE2 inhibitors could be a possible strategy to treat COVID-19, followed by fusion and entry, transmembrane protease serine 2 can activate spike proteins and plays a significant role in the process of SARS-CoV-2 infection of host cells, viral genome replication, virus proteases as 3C-like protease responsible for cleaving viral peptides into functional units for virus replication and packaging in host cells, thus SARS-CoV-2 protease inhibitors will interfere with this stage. Transcription, translation stages involving RNA-dependent RNA polymerase and reverse transcriptase which are important enzymes of the coronavirus replication–transcription complex and virion assembly, then ultimately budding and release. Natural coumarin compounds have been found to act through a diverse set of the above-mentioned mechanisms.14,17,32
A series of 218 coumarins from different families (Fig. 2) including 80 representatives with previously reported antiviral activities (series 1: AV-1 to AV-80) (Fig. 3) and 138 general coumarins (series 2: C-1 to C-140) were investigated for their potential effects against the COVID-19 virus. Initially, the compounds were screened for their physicochemical properties and their possible oral availability using SwissADME and BOILED-Egg diagrams.34 Only compounds that showed oral availability without crossing the blood–brain barrier (BBB) were selected. The rationale behind this step was that there was particular interest in compounds that could be taken orally, which is essential for mass treatment. At the same time the selected compounds should not be able to cross the BBB to avoid any possible side effects or complications in the brain. Among the first series, 28 compounds were selected, which included antiviral (AV) compounds with the following numbers: 4, 5, 7, 9, 10, 11, 16, 18, 20, 21, 25, 29, 30, 32, 33, 34, 35, 36, 37, 38, 41, 52, 53, 64, 65, 66, 67, and 78 (Fig. 4). For the second series, 46 compounds were selected, which included C-type (coumarin) compounds with the following numbers: 6, 7, 8, 9, 11, 21, 33, 34, 36, 37, 40, 44, 45, 47, 48, 49, 50, 52, 53, 56, 59, 73, 74, 75, 82, 87, 88, 89, 93, 97, 98, 99, 101, 102, 103, 104, 105, 107, 111, 115, 116, 129, 133, 134, 135, and 136 (Fig. 5). For the structures, see the ESI.†
The next step was the investigation of potential effects of these compounds against COVID-19-related targets. Using docking studies, four SARS-CoV-2 targets were investigated. These targets included Mpro, viral methyltransferase (nsp16/10 complex), RBD of S-protein and human angiotensin-converting enzyme 2 ACE2, which is the human receptor for viral S-protein. The Protein Data Bank (PDB) codes for these proteins are 6LU7, 6W4H, 6M0J, and 6VW1, respectively. After the preparation of protein targets and ligands, the compounds were docked in the active sites of the four targets. Initially, the co-crystallized ligand, if available, was redocked in the active site to validate the docking procedure and then the RMSD between the docked and the co-crystallized ligand was calculated using the DockRMSD server.24 Docking of test compounds was done after RMSD was found to reach an accepted value (<2 Å). The results of the test compounds are shown in Table 1. All results better than −8.0 kcal mol−1 were shown in bold and the top-ranked hits are shown in italics.
# | 6LU7 | 6W4H, nsp16 | 6VW1 | 6M0J |
---|---|---|---|---|
No. | Mpro | Methyl-transferase | Human ACE2 | RBD-COV-2-S protein |
AV-4 | −6.9 | −7.4 | −5.9 | −6.2 |
AV-5 | −8.4 | −7.6 | −6.5 | −6.7 |
AV-7 | −7.5 | −8.2 | −7.0 | −6.7 |
AV-9 | −6.6 | −8.1 | −5.9 | −6.4 |
AV-10 | −6.8 | −7.9 | −5.9 | −6.3 |
AV-11 | −7.2 | −8.1 | −7.2 | −6.6 |
AV-16 | −7.5 | −8.8 | −6.1 | −6.9 |
AV-18 | −7.6 | −7.4 | −6.7 | −5.6 |
AV-20 | −7.5 | −7.8 | −6.1 | −6.6 |
AV-21 | −7.7 | −7.1 | −5.9 | −6.3 |
AV-25 | −7.1 | −7.8 | −6.1 | −6.0 |
AV-29 | −7.5 | −9.2 | −6.4 | −6.7 |
AV-30 | −7.6 | −7.8 | −6.9 | −6.8 |
AV-32 | −8.2 | −7.2 | −5.7 | −6.6 |
AV-33 | −7.9 | −7.2 | −5.7 | −6.9 |
AV-34 | −8.0 | −7.1 | −5.6 | −6.9 |
AV-35 | −7.0 | −7.6 | −8.0 | −6.2 |
AV-36 | −6.2 | −6.4 | −5.7 | −5.9 |
AV-37 | −8.2 | −7.3 | −5.4 | −7.3 |
AV-38 | −7.8 | −7.2 | −5.2 | −6.5 |
AV-41 | −7.6 | −7.8 | −5.9 | −6.7 |
AV-52 | −8.5 | −8.8 | −6.6 | −7.2 |
AV-53 | −8.8 | −7.9 | −6.8 | −7.3 |
AV-64 | −6.8 | −7.8 | −6.1 | −6.2 |
AV-65 | −6.7 | −7.3 | −6.4 | −5.9 |
AV-66 | −7.6 | −8.3 | −6.7 | −6.7 |
AV-67 | −7.9 | −8.3 | −6.4 | −6.8 |
AV-78 | −7.0 | −7.5 | −5.7 | −6.1 |
C-6 | −6.4 | −8.0 | −5.8 | −5.8 |
C-7 | −7.7 | −7.6 | −6.2 | −6.4 |
C-8 | −5.9 | −6.6 | −5.8 | −5.4 |
C-9 | −7.4 | −7.6 | −5.7 | −6.1 |
C-11 | −7.9 | −7.8 | −6.2 | −6.5 |
C-21 | −7.1 | −7.6 | −5.7 | −6.6 |
C-33 | −8.1 | −9.1 | −6.6 | −7.2 |
C-34 | −7.5 | −8.8 | −6.2 | −6.7 |
C-36 | −7.6 | −8.9 | −6.5 | −7.0 |
C-37 | −7.4 | −8.4 | −7.4 | −6.5 |
C-40 | −8.1 | −9.0 | −6.8 | −7.1 |
C-44 | −7.8 | −8.8 | −6.5 | −6.9 |
C-45 | −7.4 | −9.1 | −6.8 | −7.1 |
C-47 | −7.7 | −9.0 | −7.7 | −6.9 |
C-48 | −7.8 | −9.1 | −7.2 | −7.0 |
C-49 | −8.4 | −9.8 | −6.6 | −8.2 |
C-50 | −7.3 | −8.0 | −7.9 | −6.5 |
C-52 | −8.5 | −9.1 | −6.7 | −6.9 |
C-53 | −7.8 | −8.5 | −8.0 | −6.8 |
C-56 | −8.4 | −8.2 | −5.5 | −7.1 |
C-59 | −6.5 | −6.9 | −5.8 | −5.9 |
C-73 | −6.9 | −7.5 | −5.4 | −6.0 |
C-74 | −8.4 | −8.2 | −5.6 | −6.4 |
C-75 | −8.3 | −8.2 | −5.2 | −6.2 |
C-82 | −7.0 | −7.1 | −5.9 | −5.8 |
C-87 | −8.2 | −7.9 | −6.4 | −6.6 |
C-88 | −8.7 | −8.1 | −6.2 | −6.6 |
C-89 | −8.2 | −8.2 | −5.1 | −6.2 |
C-93 | −8.2 | −8.5 | −6.2 | −7.1 |
C-97 | −8.7 | −9.1 | −6.9 | −6.8 |
C-98 | −8.0 | −8.5 | −6.8 | −6.8 |
C-99 | −7.4 | −7.7 | −5.9 | −6.6 |
C-101 | −8.4 | −8.7 | −6.2 | −6.8 |
C-102 | −8.4 | −8.0 | −6.6 | −7.5 |
C-103 | −7.5 | −8.2 | −7.3 | −6.8 |
C-104 | −8.3 | −8.5 | −6.3 | −6.8 |
C-105 | −7.6 | −7.6 | −6.1 | −6.2 |
C-107 | −6.5 | −6.4 | −6.3 | −5.6 |
C-111 | −6.4 | −7.2 | −6.9 | −6.1 |
C-115 | −6.8 | −7.5 | −6.9 | −6.2 |
C-116 | −6.7 | −7.5 | −7.0 | −6.1 |
C-129 | −6.8 | −7.2 | −6.7 | −6.0 |
C-133 | −7.3 | −8.2 | −7.3 | −6.3 |
C-134 | −6.1 | −6.5 | −5.9 | −5.7 |
C-135 | −6.7 | −7.2 | −5.2 | −5.5 |
C-136 | −6.5 | −7.0 | −6.5 | −5.7 |
−7.9 (N3) | −8.2 (SAM) | — | — |
In general, the best docking scores (reaching −9.8 kcal mol−1) were seen with SARS-CoV-2 methyltransferase. Several of the compounds investigated here displayed better docking scores than the co-crystallized ligand for both the viral main protease and the methyltransferase. The other two targets, on the contrary, showed a lower number of hits and also lower docking scores.
Screening of the compounds against the viral main protease picked 22 structures with docking scores better than that of the co-crystallized ligand N3 (Table 1). Validation of the docking procedure for Mpro was reported earlier38,39. Among these compounds, the top three representatives were 6-O-β-D-glucopyranosyl-5-hydroxyangelicin (C-88), inophyllum G2 (AV-53) and bavacoumestan A (C-97) (Fig. 7), which showed docking scores between −8.7 and −8.8 kcal mol−1. All the three compounds were able to fit in the same active site as N3 (Fig. 7A, C, and D). The top-ranked compound, 6-O-β-D-glucopyranosyl-5-hydroxyangelicin (C-88), was picked as a representative example; its specific interactions with the viral main protease are shown in Fig. 7B.
Docking of the compounds in the active site of the viral methyltransferase showed the highest docking scores among all the tested targets (Table 1). Here the top compounds were daphnorin (C-49) and glycycoumarin (AV-29) (Fig. 6), which had docking scores of −9.8 and −9.2 kcal mol−1, respectively. Redocking of the internal ligand (S-adenosyl-L-methionine, SAM) was able to predict the binding pose with a very high accuracy (Fig. 8A). Daphnorin (C-49) was docked in the same position as SAM. The terminal coumarin ring overlapped with the purine ring of SAM, while the central coumarin ring overlapped with the SAM sugar moiety (Fig. 8B). In addition, daphnorin (C-49) was able to form several hydrogen bonds similar to those formed by SAM, along with extra hydrogen bonds. The interactions of daphnorin (C-49) with amino acids in the active site of the viral methyltransferase (6W4H) are shown in Fig. 8C. Another example of compounds with high affinity is glycycoumarin (AV-29), which is displayed in Fig. 8D, where the coumarin ring is also overlapped with the purine ring of SAM.
Fig. 7 Docking results of compounds in Mpro (6LU7). (A) C-88 (salmon) overlapped with N3 (blue, co-crystallized ligand) in the active site of the main protease. (B) Interactions of C-88. (C) AV-53 (green) overlapped with N3 (blue). (D) C-97 (gray) overlapped with N3 (blue). |
Finally, docking in the active site of viral S-protein and human ACE2 was more selective. Several compounds have shown low docking scores and only limited number of compounds have displayed docking scores above −8.0 kcal mol−1. With respect to RBD of S-protein (6M0J), daphnorin (C-49) was found to be the best-docked compound with a docking score of −8.2 kcal mol−1. The compound was found to bind at the surface that should interact with human ACE2 (Fig. 9A). Its sugar part forms two hydrogen bonds with Q493, which is known to be important in the interaction with ACE2,40 suggesting its potential role in preventing the recognition between S-protein and ACE2 (Fig. 9B). RBD of S-protein bound to ACE2 (6M0J) was found to align well with SARS-CoV-2 S trimer (7DK3) with RMSD of 1.47 suggesting similar interactions. Docking in human ACE2 provided similar results with only two compounds with a score of −8.0 kcal mol−1, which include isodispar B (AV-35) and daphnogirin (C-53). Both molecules were found to bind internally away from the interacting surface between S-protein and ACE2. Binding interaction of both compounds are shown in Fig. 9A and D.
It worth to mention that docking in compounds was proceeded by preparation steps which include stripping of water molecules. This step is essential to evacuate the active site for the docking step. Effect of water molecules on binding could be further investigated using other techniques such as molecular dynamics which is beyond the scope of this article.
ACE2 | Angiotensin-converting enzyme 2 |
ADME | Absorption, distribution, metabolism, and excretion |
BBB | Blood–brain barrier |
BOILED-Egg | Brain Or IntestinaL EstimateD permeation method |
Mpro | SARS-CoV-2 main protease |
Nsp | Non-structural protein |
PDB | Protein Data Bank |
RBD | Recognition binding domain |
RdRp | RNA-dependent RNA polymerase |
RMSD | Root-mean-square deviation |
SAM | S-Adenosyl-L-methionine |
TMPRSS2 | Type II transmembrane serine protease |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ra01989a |
This journal is © The Royal Society of Chemistry 2021 |