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 was summarized; Table S1. Among 101 different bioactives; 36 bioactives 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 enrichment analysis
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. Similarly, 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, 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 Figure 1. Similarly, the combined bioactives-proteins-pathways also reflected the anthrogallol to target the highest number of proteins. Likewise, TNFRSF1A was majorly targeted by the highest number of bioactives modulating the pathways in the cancer pathway (Figure 2).
3.3. Prediction of probable anti-viral activity
The anthraquinones were found to be anti-viral agents against Adenovirus, CMV, Hepatitis B, Hepatitis 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 Figure 3.
3.4. In silico molecular docking
Torososide B was predicted to have the highest binding affinity (–8.7 kcal/mol) with Papain-Like Protease with 9 hydrogen bond interactions with THR302, ASP303, TYR274, TYR265, ARG167, TYR269, ASP165. Further, Torososide B was predicted to possess the highest binding affinity (–9.3 kcal/mol) with Coronavirus Main Proteinase with 14 hydrogen bond interaction with LEU287, TYR237, THR199, ARG131, LYS137, LYS5, GLU290, ILE281, LEU282, PHE3. Similarly, 1,3,6-trihydroxy–2-methyl–9,10-anthraquinone–3-O-(6’-O-acetyl)-β-D-xylopyranosyl-(1->2)-β-D-glucopyranoside was predicted to possess the highest binding affinity (–8.7 kcal/mol) with spike protein with the highest number of hydrogen bond interactions i.e. 6 with ASP820, ILE816, ASP815, GLN825, MET703. The binding affinity of each compound with individual targets with hydrogen bond interactions and residues is summarized in Table S2. The interaction of Torososide B with Papain-Like Protease and Coronavirus Main Proteinase 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 Figure 4.