Safae El Mazouria, Tarik Aanniza*, Jihane Touhtouha, Ilham Kandoussia, Mohammed Hakmia, Lahcen Belyamanib, Azeddine Ibrahimia, Mouna Ouadghiria
aMedical Biotechnology Laboratory (MedBiotech), Bioinova Research Center, Medical and Pharmacy School, Mohammed Vth University, Rabat, Morocco
bEmergency Department, Military Hospital Mohammed V, Medical and Pharmacy School, Mohammed Vth University, Rabat, Morocco
*Corresponding Author: Tarik Aanniz, Medical Biotechnology Laboratory (MedBiotech), Bioinova Research Center, Medical and Pharmacy School, Mohammed Vth University, Rabat, 10102, Morocco
Received: 25 February 2021; Accepted: 08 March 2021; Published: 18 March 2021
The coronavirus disease 19 (Covid-19) pandemic caused by the SARS-CoV-2 virus has become a humanitarian crisis. Considering the severity of the situation we have performed a virtual screening of anthraquinones derivative drugs and phytochemicals targeting simultaneously multiple essential proteins of SARS-CoV-2 namely Mpro, PLpro, RdRp and the spike.
Among the 9 screened anthraquinones derivative drugs, valrubicin, idarubicin, daunorubicin, doxorubicin, epirubicin and diacerein were the most potent inhibitors of SARS-CoV-2 Mpro, PLpro, RdRp and Spike simultaneously. Valrubicin has the best affinity towards the spike protein (-9.5 kcal/mol), RdRp (-8.2 kcal/mol) and PLpro (-7.9 kcal/mol) while idarubicin and doxorubicin were the most effective against Mpro (-8.3 kcal/mol). No toxicity measurements are required for these drugs since they were tested prior to their approval by the FDA. Of the 140 screened phytochemicals anthraquinones were the most potent candidates. Hypericin and rhein were able to bind to the active site of all four targets, while chrysophanol, aloesaponarin II, emodine, aloe-emodine, physcion and danthron simultaneously bound to the active site of SARS-CoV-2 Mpro, Spike and RdRp. Hypericin showed the best affinity towards the spike protein (-9.7 kcal/mol), RdRp (-10.2 kcal/mol) and PLpro (-7.8 kcal/mol), while chrysophanol was the most effective one against Mpro (-8.4 kcal/mol).
Our overall prediction findings indicate that anthraquinones may inhibit the activity of the four essential proteins of SARS-CoV-2 and those results can pave the way in drug discovery.
Anthraquinones; Virtual screening; Spike; RdRp; PLpro; Mpro
Anthraquinones articles; Virtual screening articles; Spike articles; RdRp articles; PLpro articles; Mpro articles
This work was carried out under national funding from the Moroccan Ministry of Higher Education and Scientific Research (COVID-19 program) to A.I. This work was also supported by a grant from the Moroccan Institute of Cancer Research and the PPR-1 program to A.I.
The world is facing a pandemic named Covid-19 caused by the novel coronavirus known as SARS-CoV-2 that emerged in the city of Wuhan, China by the end of 2019. Coronaviruses have repeatedly evolved during the past 1000 years [1]. They represent a diverse family of positive-sense RNA viruses capable of infecting the pulmonary system, the intestines, the liver and the nerve cells in human and animal hosts [2]. Almost a year after the first case detected, there are still no approved specific antiviral drugs against Covid-19. Notably, many antiviral, anti-malaria and anti-inflammatory molecules have been tested, with very limited success [3;4].
Even in the presence of approved vaccines, antivirals are needed in some specific circumstances. Many people are not allowed to be vaccinated such as pregnant women, immunocompromised and people with severe allergic reactions. Manufacturing and distributing a vaccine at large scale will present significant challenges. The use of monoclonal antibodies is a promising outlook and provide an effective and specific therapeutic option. Nevertheless, large-scale production is complex, time consuming and expensive, meaning poor countries might be priced out and only wealthy countries could afford them [5]. Monoclonal antibodies are administrated relatively at high doses and are difficult for generic-drug makers to duplicate.
The process of a new drug development is cost and time consuming, scientists are trying to solve the puzzle through drug repurposing, the most realistic in the present pandemic, or the use of compounds from natural sources, offering novel testable hypotheses [6;7]. Anthraquinones, structurally related to 9,10-dioxoanthracene, are an important privileged structure in medicinal chemistry for more than a millennium. They have been proven to serve as a suitable starting point for optimization leading to the drug development. Anthraquinones have been known and used as laxatives for centuries and currently are employed as antiarthritics, anti-inflammatory, antifungal, antiviral, antimalarial, antimicrobial, antiplatelet, antidiabetic, neuroprotective, antioxidants, antibacterial, for the treatment of multiple sclerosis, as immune boosters [8-10] and as anticancer agents. Marketed drugs containing anthraquinone scaffold include amrubicin, daunorubicin, diacerein, doxorubicin, epirubicin, idarubicin, mitoxantrone, and valrubicin [10].
Anthraquinones may produce potential damage to cells because of the close similarity with the toxic analogue, anthracene [11]. However, there is no clear-cut correlation, and the substitution pattern appears to be an important determinant and all serious side effects have typically been the result of abuse and chronic use [9]. Few studies indicated that some anthraquinones may inhibit coronavirus activity by targeting different steps of the infection process [12-14]. Anthraquinones are well known for their immunity boosting properties, making them valuable to tackle the Covid-19 [8]. Nevertheless, few studies attempted to identify simultaneous inhibitors of these targets.
Overall, these observations prompted us to investigate the potential effect of anthraquinones on SARS-CoV-2 infection. The present work concentrates on the effectiveness of anthraquinones as simultaneous multi-target binders to the Mpro, PLpro, RdRp and Spike of the SARS-CoV-2.
2.1. Data sources
A dataset of 9 anthraquinones derivative drugs and 140 phytochemicals were screened against 4 crucial SARS-CoV-2 targets namely Mpro, PLpro, RdRp and Spike proteins. Phytochemicals database was obtained based on Lipinski’s rule of five. Molecular docking was performed using Autodock Vina to predict the ability of small molecules to bind to the active site of the used targets to identify potent inhibitors.
2.2. Preparation of receptors
The three-dimensional crystal structure of SARS-CoV-2 Mpro (PDB ID: 6lu7), PLpro (PDB ID: 6xaa), RdRp (PDB ID: 7bv2) and the spike protein (PDB ID: 6lzg) were collected from the RCSB Protein Data Bank (PDB). The crystal structures were prepared according to Pitsillou et al. (2020) [15]. The solvation parameters, charges assignment, fragmental volumes, and protein optimization were checked [16-18]. Kollman united atom charges, solvation parameters and polar hydrogens were added to the targets using AutoDockTools [19].
2.3. Preparation of ligands
The ligand molecules were retrieved from PubChem database [20] and were converted to PDBQT format [21]. Gasteiger charges were assigned and then non-polar hydrogens were merged.
2.4. Virtual screening and molecular docking protocols
For all of the targets, a grid box (X: 20 Å, Y: 20 Å and Z: 20 Å) was generated around the active site residues. The center of the grid box for the spike protein was X: -26.84, Y: 18.32 and Z: -14.06, for Mpro was X: -26.28, Y: 12.59, Z: 58.96, for PLpro was X: -24.83, Y: 20.63, Z: -1.87 and for RdRp was X: 98.76, Y: 98.79, Z: 98.98. During the docking process using Autodock Vina, the protein was considered to be rigid and the ligand to be flexible [7;21]. The binding affinities of the drugs were measured in kcal/mol unit and sorted according to the higher negative values, which imply the best binding affinities [22].
2.5. Analysis and visualization
Based on calculated binding affinities, the ligands with the lowest binding energy (less than -7 kcal/mol) were chosen for the interactions analysis using PyMol for 3D visualization and BIOVIA Discovery Studio Visualizer for 2D visualization.
The rapid development and identification of efficient therapies against SARS-CoV-2 remains a major challenge. Repurposing of available drugs and the use of natural compounds always offer a great promise [7;8;23-25]. Currently, in silico studies are providing lots of preliminary data about potential drugs, which can be a great help in further in vitro and in vivo studies [15].
In this study, we screened 9 anthraquinones drugs and 120 phytochemicals against 4 targets of SARS-CoV-2 using structure based virtual screening approach towards four targets of the SARS-CoV-2 i.e., PLpro, Mpro, RdRp and the spike protein. A potent effect leading to more efficiency against SARS-CoV-2 could be the result of the multi-target inhibition by anthraquinones.
Anthraquinones derivative drugs repurposing
Considering -7 kcal/mol as a cutoff value, anthraquinones drugs i.e., valrubicin, idarubicin, daunorubicin, doxorubicin, epirubicin and diacerein exhibited the best binding affinity simultaneously towards all of the four targets. The binding energy of all the potent anthraquinones drugs are shown in Table 1 and their structures are shown in Figure 1. Valrubicin is the drug with the best affinity towards the spike protein (-9.5 kcal/mol), RdRp (-8.2 kcal/mol) and PLpro (-7.9 kcal/mol) while idarubicin and doxorubicin were the most effective against Mpro (-8.3 kcal/mol). Whereas, mitoxantrone and pixantrone showed low affinities towards all of the four targets. No toxicity measurements are required for these drugs since they were tested prior to their approval by the FDA.
Table 1: The binding energy and H-bonded interactions between anthraquinones derivative drugs and the SARS-CoV-2 targets
Table 2: The binding energy and H-bonded interactions between natural anthraquinones and the SARS-CoV-2 targets
Natural anthraquinones screening
Among the 140 phytochemicals screened, anthraquinones were the most potent inhibitors of the used targets. Anthraquinones namely rhein and hypericin also exhibited good binding affinities simultaneously towards all of the four targets, while chrysophanol, aloesaponarin II, emodin, aloe-emodin, danthron and physcion docked all the targets except PLpro (Figure 1; Table 1). Hypericin showed the best affinity towards the spike protein (-9.7 kcal/mol), RdRp (-10.2 kcal/mol) and PLpro (-7.8 kcal/mol), while chrysophanol exhibited the best affinity towards Mpro (-8.4 kcal/mol). Details of the binding energy of the best natural anthraquinones are tabulated in Table 2 and their structures are shown in Figure 1. The remaining phytochemicals showed low affinities. All of the prioritized natural anthraquinones are obeying Lipinski’s rule of five. The toxicity of anthraquinones is well documented and discussed [9;26].
Molecular interaction studies
All of the screened molecules occupied the active site of the selected targets. Researchers suggest that non-covalent interactions with the active site of the protein may lead to possible inhibition and blockage of the target.
Valrubicin, idarubicin, daunorubicin, doxorubicin and amrubicin have shown numerous interactions with the active site of the spike protein including 6 H-bonds, diacerein formed 5 H-bonds, while epirubicin formed 4 H-bonds. Idarubicin, doxorubicin and epirubicin formed 5 H-bonds with the active site of Mpro, valrubicin, daunorubicin, diacerein formed 2 H-bonds while the complex amrubicin-Mpro lacked hydrogen bonding interactions. Daunorubicin was interestingly locked into the binding site of RdRp, forming 7 H-bonds, while epirubicin formed 6 H-bonds, followed by idarubicin, valrubicin and amrubicin that formed 5 H-bonds, doxorubicin formed 4 H-bonds while the complex diacerein-Mpro was stabilized with 3 hydrogen bonding interactions. The interactions of diacerein, daunorubicin and epirubicin with the active site of PLpro were maintained by 4 H-bonds, valrubicin formed 3 H-bonds, idarubicin formed 2 H-bonds while doxorubicin and amrubicin formed one H-bond.
Details regarding polar interactions involving hydrogen-bonds are shown in the 3D interaction plots of Figure 2 and are tabulated in Table 1.
As regards to natural anthraquinones, 6 H-bonds can be seen between aloe-emodin and the active site of the spike protein, rhein, aloe-saponarin II, and emodin formed 4 H-bonds, hypericin and chrysophanol formed 3 H-bonds while physcion and danthron formed 2 H-bonds. The interactions between the active site of the Mpro and aloe-emodin and hypericin consisted of 5 H-bonds, chrysophanol and emodin formed 4 H-bonds, aloe-saponarin II formed 3 H-bonds, rhein and physcion formed 2 H-bonds while the complex danthron-Mpro was stabilized with one H-bond. Rhein, physcion and emodin formed 4 H-bonds with the active site of the RdRp, aloe-emodin, hypericin and chrysophanol formed 3 H-bonds, while aloe-saponarin II and danthron formed 2 H-bonds. 4 H-bonds were formed between rhein and the active site of the PLpro, while hypericin formed 3 H-bonds. The interactions between the remaining natural anthraquinones and PLpro were not considered as they show low affinities.
Details regarding polar interactions involving hydrogen-bonds are shown in the 3D interaction plots of Figure 3 and are tabulated in Table 2.
The use of drugs with known mechanisms of action has emerged as a valid tool to identify compounds that can contribute to contain Covid-19 [7;27-29]. Assessing evidences from molecular docking studies, it was clearly seen that anthraquinones derivative drugs i.e. valrubicin, idarubicin, daunorubicin, doxorubicin, epirubicin and diacerein exhibited good binding affinity towards PLpro, Mpro, RdRp and the spike simultaneously indicating their competency of inhibiting SARS-CoV-2.
Valrubicin, an antitumoral, was identified as potent inhibitor of SARS-CoV-2 Mpro [30]. Our study showed that valrubicin has a good affinity towards the spike protein as well as RdRp, PLpro and Mpro respectively, which allows us to suggest that valrubicin could be a promising inhibitor of SARS-CoV-2. In another study, idarubicin exhibited strong affinity towards endoribonuclease (EndoU) [31]. Here we showed that idarubicin can also bind to PLpro, Mpro, RdRp and the spike protein indicating its multi-target effect against SARS-CoV-2. Daunorubicin was reported to be a potential inhibitor of SARS-CoV-2 Mpro based on docking studies performed by Jimenez-Alberto et al. (2020) [32]. The current study showed that daunorubicin has also high affinity towards RdRp, Spike and PLpro. Al- Motawa et al. (2020) [33] reported that doxorubicin, an antitumor drug with historical antiviral activity, could be used to defeat SARS-CoV-2 by increasing cellular methylglyoxal (MG) concentration to virucidal levels. Relatively short-term treatment with doxorubicin increasing cellular MG may be beneficial in patients with Covid-19. Our finding brings additional supports and mechanisms to these results by proving that doxorubicin exhibits a very good affinity towards the four crucial SARS-CoV-2 proteins (PLpro, Mpro, RdRp and the spike protein). Epirubicin used for chemotherapy to treat breast cancer was identified as a potential inhibitor of SARS-CoV-2 Mpro in a study conducted by Khan et al. (2020) [28] which promotes our docking results showing that epirubicin has an inhibitory effect on SARS-CoV-2 Mpro, in addition to RdRp, Spike and PLpro. de Oliveira et al. (2020) [34] hypothesized that diacerein is a multi-target drug useful for Covid-19 treatment. By inhibiting the IL-1β production, this drug may have an impact on viral infection, disease onset, progression and outcome by controlling hyperinflammatory conditions [35, 36]. Amrubicin, used in the treatment of lung cancer was identified as an inhibitor of SARS-CoV-2 Mpro by Jimenez-Alberto et al. (2020) [32]. Moreover, in the present work we found that amrubicin has a good affinity towards Spike, RdRp and PLpro.
Many studies showed that anthraquinones like hypericin, rhein, emodin, aloe-emodin and chrysophanol could hinder the entry, replication and release of many viruses including coronaviruses. The antiviral activity of emodin and rhein has well been documented. Rhein from Rheum palmatum ethanol extract inhibited hepatitis B virus (HBV) replication [37]. Barnard et al. (1992) [38] reported that emodin, rhein and hypericin inhibited human cytomegalovirus (HCMV) while Chang et al. (2014) [39] reported the antiviral activity of aloe-emodin, chrysophanol, rhein, emodin and physcion against Japanese encephalitis virus (JEV). Lin et al. (2005) [40] showed also that aloe-emodin, emodin, and chrysophanol inhibited the SARS-CoV 3CLpro.
Interest in hypericin is ascribed to the discovery that it possesses extremely high activity towards enveloped viruses including Hepatitis C virus, herpes simplex virus, murine cytomegalovirus, sindbis virus, hepatitis, some leukemic viruses, and avian infectious bronchitis virus (IBV). Additionally, it has also been used in the clinical treatment of AIDS patients [41-47]. Balmeh et al. (2020) [48] showed that hypericin has a high affinity towards ACE2, GRP78, AT1R, and TMPRSS2 making it a promising compound to prevent the Covid-19 infection. Hypericin, could also bind the N-terminus and C-terminus of the SARS-CoV-2 Nsp14 [49]. Here we report that hypericin could effectively bind to the active site of SARS-CoV-2 Spike, RdRp Mpro and PLpro which corroborates its antiviral potential.
Rhein suppresses lung inflammatory injury induced by human respiratory syncytial virus, through inhibiting NLRP3 inflammasome activation via NF-κB pathway in mice [50]. Influenza A virus (IAV) adsorption and replication was inhibited by rhein through decreasing IAV-induced oxidative stress and downregulating the activation of the TLR4, Akt, p38, JNK MAPK, and NF- κB-mediated signaling pathways and the production of inflammatory cytokines and matrix metalloproteinases [50]. The inhibitory effect of this molecule on IL-6 and TNF-α has been reported [34]. Importantly, rhein inhibits the interaction between the SARS-CoV-2 S protein and ACE2 suggesting that rhein is a potential therapeutic agent for the treatment of SARS-CoV-2 [12]. Preincubation of rhein with biotinylated S protein inhibited S protein-ACE2 interaction. Savarino et al. (2007) [13] reported the effect of rhein on the activity of the human liver enzyme cathepsin B within endosomes which is a critical step in SARS-CoV-2 infection [14]. Rhein was also reported for exhibiting anti-coagulant effects [34]. In our study, we showed that rhein was also able to dock the active site of SARS-CoV-2 PLpro, Mpro, RdRP and the spike protein. All these findings make rhein advantageous as an inhibitor of SARS-CoV-2.
Different antiviral activities of chrysophanol towards poliovirus type 2 and 3, Coxsackie virus type A21 and B4, human rhinovirus type 2, Ross River virus and herpes simplex virus type 1, have been investigated [51]. Subbaiyan et al. (2020) [52] found that this compound has one of the highest binding affinities towards the viral spike protein which is consistent with our findings.
The antiviral properties of emodin have also been confirmed, mainly on enveloped viruses. Emodin protected against H1N1 [53]. Batista et al. (2019) [54] showed that emodin triggers a strong virucidal effect on ZIKA by reducing the hydrodynamic radius of virus particle in solution. Emodin could hinder the SARS-CoV by blocking both viral entry and release [12]. Luo et al. (2009) [55] showed that emodin and aloe-emodin blocked the interaction between S protein of SARS-CoV and ACE2. In another study, emodin inhibited the SARS-CoV 3a protein and S protein [56; 57]. Furthermore, here we report that emodin could inhibit PLpro, Mpro, RdRP and the spike protein which offers a potential therapeutic approach for SARS-CoV-2 [12].
Li et al. (2014) [58] reported the antiviral activity of aloe-emodin against influenza A virus by up-regulating galectin-3, and thioredoxin as well as down-regulating nucleoside diphosphate kinase A while Lin et al. (2008) [59] reported that it acted as an interferon-inducing agent with antiviral activity against Japanese encephalitis virus and enterovirus 71. In addition, here we show that aloe-emodin could bind perfectly to the active site of Mpro, RdRp and the spike protein making it a potent inhibitor of SARS-CoV-2.
Otherwise, little is known about the antiviral activity of physcion, danthron and aloesaponarin II. Here we report also that they are able to perfectly dock Mpro, and RdRp proteins suggesting them as potent inhibitors of SARS-CoV-2.
Safety issues of anthraquinones should definitely be considered while developing derivatives for medical use, but it should be borne in mind that toxicity is based on the entire structure of a molecule. Moreover, the cytotoxicity was found to be related to the dose and the duration of treatment [60;61]. Drug treatment of Covid-19 may be shorter than cancer chemotherapy: for example, median hospitalization time of patients surviving severe symptoms of Covid-19 was 28 days [62] while a typical course of cancer chemotherapy is 6 months or longer [60;61].
To sum up, the present study corroborates that anthraquinones are promising scaffolds for anti-SARS-CoV-2 antivirals development. Further validation by in vitro and in vivo studies are needed. If a high potency antiviral effect of these agents is found, low dose and short duration of treatment are expected to decrease the risk of adverse effects.
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
This work was carried out under national funding from the Moroccan Ministry of Higher Education and Scientific Research (COVID-19 program) to A.I. This work was also supported by a grant from the Moroccan Institute of Cancer Research and the PPR-1 program to A.I.