Anthraquinone Derivatives as an Immune Booster and their Therapeutic Option Against COVID-19

Pukar Khanal ·B. M. Patil ·Jagdish Chand ·Yasmin Naaz

Abstract

Keywords 3CLpro·Anthroquine derivatives·Coronavirus·COVID-19·Immune boost

1 Introduction

Presently, CoV Disease (COVID-19) has lead to millions of death throughout the world beginning from the December of 2019 [1].Further, the risk of getting affected with COVID-19 is supplemented in subjects with lower immunity primarily special subjects like pediatrics/geriatrics and the patients suffering from infectious and non-infectious diseases [2].Although approaches are being made to prevent this virus to get spread via social distancing and so on, boosting of immunity in subjects could play an important role in inhibiting the transmission of the virus and its invasion into the body.Although investigations are undergoing to develop the vaccine against COVID-19, it may still take time as the drug discovery process is much complicated.Hence, it is important to identify any alternative approach as prophylaxis against COVID-19 which can be implemented via the utilization of immune boosters from natural sources.Three targets of novel coronavirus i.e.3C-like protease (3CLpro),papain-like protease (PLpro), and spike protein [3—6] are being targeted by multiple investigators to identify the new hit molecule for the management of COVID-19.Further,inflammation and cell necrosis are contributing factors in worsening the COVID-19 pathogenesis which suggests identifying the molecule with anti-viral, anti-inflammatory and immune-boosting properties.

Anthroquinolines are the group of compounds from multiple folk medicines likeSennaspecies which are utilized in Ayurvedic system of medicines and Traditional Chinese Medicines for the management of various infectious and non-infectious diseases [7, 8].Further, anthraquinone derivatives are also reported for anti-viral property [9], antiinflammatory effi cacy [10], and as immune booster [11].Hence, in the COVID-19 infection, it may be beneficial if the bioactives are identified with an immune boost, anti-inflammatory, and anti-viral properties which can be demonstrated via the concept of network pharmacology or polypharmacological approach.Hence, based on the above theme, of anti-viral/anti-inflammatory/immune boosting reports of anthraquinones we attempted to screen the multiple anthraquinone derivatives as an immune booster and anti-viral effi cacy using in silico molecular docking and other system biology tools.

2 Materials and methods

2.1 Bioactives and their Druglikeness Score

The reported phytoconstituents under the phytochemistry of anthraquinone were retrieved from available literature/Chemical Entities of Biological Interest (ChEBI) records(https://www.ebi.ac.uk/chebi/).All the compounds were then predicted for the druglikeness score by querying the SMILES of each molecule in MolSoft (https://molso ft.com/mprop/).

2.2 Target Prediction and their Enrichment Analysis to Assess Immune-boosting Effi cacy

Anthraquinone derivatives with positive druglikeness scores were queried in DIGEP-Pred [12] to identify “Proteins based targets” (up-regulated/downregulated proteins) at the probable activity of 0.5.The list of regulated proteins was queried in STRING [13] to identify the biological process, cellular function, and molecular processes of combined gene-set.Further, the probably modulated pathways were also identified concerning the Kyoto Encyclopedia of Genes and Genomes database.Network among the bioactives, their targets, and modulated pathways was constructed using the Cytoscape [14] version 3.5.1.Any duplicates interconnection of two nodes was eliminated to avoid the appearance of a false hit.The whole network was analyzed using the“network analyzer” tool based on node size and count representing the edge count as “low values to small size” and“low values to bright colors” respectively as explained previously [15].

2.3 Prediction of Probable Anti-viral Activity

Anti-viral activity of each compound was predicted by querying the SMILES in Prediction of Activity Spectra for Substances [16] at the pharmacological activity (Pa) >pharmacological inactivity (Pi) and retrieving the probable biological spectrum for keyword “anti-viral”.Records were queried for their probable pharmacological spectrum against multiple viruses like Adenovirus, Cytomegalovirus (CMV), Hepatitis B, Hepatitis C, Hepatitis, Herpes,Human immunodeficiency virus (HIV), Influenza A, Influenza, Parainfluenza, Picornavirus, Poxvirus, Rhinovirus, and Trachoma.

2.4 In Silico Molecular Docking

2.4.1 Preparation of Ligand Molecules

All the 3D.sdf format of ligand molecules with positive drug-likeness scores was retrieved from PubChem database(https://pubch em.ncbi.nlm.nih.gov/) or structures were drawn in ChemSketch (https://www.acdla bs.com/resou rces/freew are/chems ketch/) as applicable.The ligand molecules were converted into.pdb format using Discovery studio 2019[17].All the bioactives were minimized using MMFF94 forcefield [18] using conjugate gradients as an optimization algorithm.After the minimization of energy, all ligand molecules were converted into.pdbqt format.

2.4.2 Preparation of Macromolecules

Three target proteins of COVID-19 i.e.3clpro (PDB:6LU7),PLpro (PDB:4M0W), and spike proteins (homology modeled target, accession number:AVP78042.1 as query sequence and PDB:6VSB as a template using SWISSMODEL [19]) were selected.The retrieved proteins from Research Collaboratory for Structural Bioinformatics database were in complex with other heteroatoms which were removed using Discovery studio 2019 and saved in.pdb format.

2.4.3 Ligand-protein Docking

The ligand molecules were docked with protein molecules using autodock4 [20] by setting 8 exhaustiveness as default to obtain 10 different poses of ligand molecules.After docking the pose of ligand scoring lowest binding energy was selected to visualize the ligand—protein interaction in Discovery Studio 2019 as explained previously [21, 22].

3 Results

3.1 Bioactives and their Druglikeness Score

The complete datasheet of 101 compounds including their name, ChEBI ID, molecular formula/weight, and synonym including phytochemistry for the retrieved compounds were summarized (Table S1).Among 101 different compounds,36 were identified with positive druglikeness scores.Among them, laccaic acid A scored highest druglikeness score i.e.0.85 with molecular weight 537.09, 12 hydrogen bond acceptor, 8 hydrogen bond donors, and 2.88 MolLogP.The details of the druglikeness score of each compound are summarized in Table 1.

3.2 Target Prediction and their Enrichment Analysis to Assess Immune-Boosting Effi cacy

Among the compounds with positive druglikeness score,anthragallol was predicted to modulate the highest number of genes i.e.25.Similarly, human carbonyl reductase 1 (CBR1) was targeted by the highest number of bioactives i.e.33.Further, the enrichment analysis identified the modulation of 54 different pathways in which pathways in cancer was majorly modulated by regulating 12 genes(AR, CASP8, CDK4, CTNNB1, EPAS1, HMOX1, KLK3,MMP2, NFE2L2, NOS2, RAC1, and RARA) under the background of 515 proteins at the false discovery rate of 7.71E−08.Table 2 summarizes the gene enrichment analysis of the modulated gene set along with modulated pathways with their respective gene codes.The protein—protein interaction of modulated proteins is presented in Fig.1.Similarly,the combined bioactives-proteins-pathways also reflected the anthrogallol to target the highest number of proteins.Likewise, TNFRSF1A and pathways in cancer were majorly targeted/modulated protein and pathways respectively (Fig.2).

3.3 Prediction of Probable Anti-viral Activity

The anthraquinones were found to be anti-viral agents against adenovirus, CMV, hepatitis B and C, hepatitis, herpes, HIV, influenza A, influenza, parainfluenza, picornavirus, poxvirus, rhinovirus, and trachoma.Among them, the majority of the compounds were active against herpes virus i.e.13.28%.The overall activity of compounds against multiple viruses is summarized in Fig.3.

3.4 In Silico Molecular Docking

Torososide B was predicted to have the highest binding affi nity (−8.7 kcal/mol) with PLpro with 9 hydrogen bond interactions via THR302, ASP303, TYR274, TYR265,ARG167, TYR269, ASP165.Further, Torososide B was predicted to possess the highest binding affi nity (−9.3 kcal/mol) with 3CLpro with 14 hydrogen bond interactions with LEU287, TYR237, THR199, ARG131, LYS137,LYS5, GLU290, ILE281, LEU282, and PHE3.Similarly,1,3,6-trihydroxy-2-methyl-9,10-anthraquinone-3-O-(6′-Oacetyl)-β-D-xylopyranosyl-(1—>2)-β-D-glucopyranoside was predicted to possess the highest binding affi nity (−8.7 kcal/mol) with spike protein with the highest number of hydrogen bond interactions i.e.6 with ASP820, ILE816, ASP815,GLN825, and MET703.The binding affi nity of each compound with individual targets with the number of hydrogen bond interactions and residues is summarized in Table S2.The interaction of Torososide B with PLpro and 3CLpro and 1,3,6-trihydroxy-2-methyl-9,10-anthraquinone-3-O-(6′-O-acetyl)-β-D-xylopyranosyl-(1—>2)-β-D-glucopyranoside with spike protein is presented in Fig.4.

4 Discussion

During the COVID-19 infection, severe necrosis and inflammation lead to defects in the supply of necessary nutrients and oxygen into the cells which are more terrible in the subjects with compromised immunity.Hence, in the present study, we investigated multiple anthraquinone derivatives from various traditional medicines to act against COVID-19 targets i.e.3CLpro, PLpro, and spike protein, and their combined immune-boosting effi cacy.Initially, we calculated the druglikeness score of each molecule based on the“Lipinksi’s Rule of Five” [23] as the majority of the plantbased medicines are utilized via the oral route which identified 36 different compounds with positive druglikeness score and considered to get absorbed orally (Table 1) which were contemplated for further study.

The conventional drug discovery process utilizes the concept of “single drug-single protein-single disease” [24]which may not be applicable in the management of infectious diseases.This is due to the affi nity of the pathogens(viruses/bacteria) to affect the multiple homeostatic functions of protein molecules.It means multiple proteins from the pathogens are involved to generate this effect.Hence, this can be managed via the utilization of modified drug development process “multi compound-multi protein-disease”interaction in which multiple bioactives regulate multiple proteins [25] which can also be taken as a basic key of boosting the immune system.Hence, in the present study, the combined synergistic phenomena of anthraquinones were investigated rather than a single bioactive molecule to identify multiple pathways that are directly or indirectly involved in the immune system.

Gene set enrichment analysis identified multiple pathways like p53 signaling pathway [26], PI3K-Akt signaling pathway [27], Rap1 signaling pathway [28], NF-kappa B signaling pathway [29] which are directly involved in the boosting the immune system.Similarly, some other pathways like pathways in cancer, PPAR signaling pathway, colorectal cancer, chemical carcinogenesis, estrogen signaling pathway were also identified which reflects the potency of anthraquinones to be beneficial in the subjects which are suffering from these pathways associated diseases like cancer.Further,pathways like p53 signaling pathways, PI3K-Akt, Wnt signaling pathways are also associated with diseases like diabetes and obesity where the immune system is compromised.Hence, regulation of these pathways could be beneficial inmanaging the diseases from which they are suffering, and boosting the immune system will also act as prophylaxis against COVID-19.Additionally, the enrichment analysis also identified the modulation of multiple pathways that are associated with a pathogenic infection like viral myocarditis and tuberculosis which reflects the potency of anthraquinone derivatives to manage the infectious diseases.Further, herbal medicines rich in anthraquinones also possess the anti-viral potency against various viruses.Hence, we attempted to identify the probable anti-viral activity of the anthraquinones with positive druglikeness score which identified their effi cacy against multiple viruses like the rhino, influenza,herpes trachoma, pox, and CMV.

Table 1 (continued)

3CLpro alters the ubiquitin system to incorporate the viral polypeptides and deregulates the homeostatic task of functional proteins [30] which was majorly targeted by torososide B.Further, PLpro alters the function of protein phosphatase 1A and protein phosphatase 1B into the replicase proteins to adjust viral life cycle [31] which was majorly inhibited by torososide B.Similarly, spike protein utilizes angiotensin-converting enzyme 2 (ACE-2)as a receptor to enter inside the host cell [32, 33] which was chiefly regulated by modulated by 1,3,6-trihydroxy-2-methyl-9,10-anthraquinone-3-O-(6′-O-acetyl)-β-Dxylopyranosyl-(1→ 2)-β-D-glucopyranoside.These results reflect the probability of the anthraquinone derivatives to act as the anti-viral against COVID-19.However, as the time proceeds, it is to be understood that the binding affi nity of probable lead hit molecules may get altered due to mutation in the possible protein targets and the inhibitory function may not occur as predicted.

5 Conclusion

The present study utilized in silico molecular docking tools to identify the binding affi nity of previously recorded anthraquinones derivatives against 3clpro, PLpro, and spike protein which identified Torososide B and 1,3,6-trihydroxy-2-methyl-9,10-anthraquinone-3-O-(6′-O-acetyl)-β-Dxylopyranosyl-(1→ 2)-β-D-glucopyranoside as a lead hits.Similarly, the combined synergies of the network identified the modulation of multiple pathways involved in the immune system like p53, chemokine, and PI3K-Akt signaling pathways.Additionally, anthraquinone derivatives were also identified as the modulators of the disease pathways where the immune system is compromised like diabetes and obesity.All these results suggest the probable therapeutic option of the utilization of anthraquinones as an immune booster and anti-viral against novel coronavirus by acting on the three targets as investigated.However, the present findings are only based on the computer simulations which need to be further validated using well designed experimental protocols.

Table 2 Enrichment analysis of modulated proteins by the reported anthraquinone derivatives

Table 2 (continued)

AcknowledgementsAll the authors are thankful to Principal KLE College of Pharmacy, Belagavi, KAHER Belagavi for his support for this completion of this work.

Author ContributionsPukar Khanal developed the protocol, performed the work, and drafted the manuscript.Prof.B.M.Patil reviewed and finalized the manuscript draft.Jagdish Chand and Yasmin Naaz had an equal contribution to mining the database and assisting in the work performed.

FundingNot available.

Compliance with Ethical Standards

Conflict of interestThere is no conflict of interest.

Research Involving Human and Animal ParticipantsThis manuscript doesn’t include any animal or human studies.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creat iveco mmons.org/licen ses/by/4.0/.