- Open Access
Network pharmacology-based predictions of active components and pharmacological mechanisms of Artemisia annua L. for the treatment of the novel Corona virus disease 2019 (COVID-19)
BMC Complementary Medicine and Therapies volume 22, Article number: 56 (2022)
Novel Corona Virus Disease 2019 (COVID-19) is closely associated with cytokines storms. The Chinese medicinal herb Artemisia annua L. (A. annua) has been traditionally used to control many inflammatory diseases, such as malaria and rheumatoid arthritis. We performed network analysis and employed molecular docking and network analysis to elucidate active components or targets and the underlying mechanisms of A. annua for the treatment of COVID-19.
Active components of A. annua were identified through the TCMSP database according to their oral bioavailability (OB) and drug-likeness (DL). Moreover, target genes associated with COVID-19 were mined from GeneCards, OMIM, and TTD. A compound-target (C-T) network was constructed to predict the relationship of active components with the targets. A Compound-disease-target (C-D-T) network has been built to reveal the direct therapeutic target for COVID-19. Molecular docking, molecular dynamics simulation studies (MD), and MM-GBSA binding free energy calculations were used to the closest molecules and targets between A. annua and COVID-19.
In our network, GO, and KEGG analysis indicated that A. annua acted in response to COVID-19 by regulating inflammatory response, proliferation, differentiation, and apoptosis. The molecular docking results manifested excellent results to verify the binding capacity between the hub components and hub targets in COVID-19. MD and MM-GBSA data showed quercetin to be the more effective candidate against the virus by target MAPK1, and kaempferol to be the other more effective candidate against the virus by target TP53. We identified A. annua’s potentially active compounds and targets associated with them that act against COVID-19.
These findings suggest that A. annua may prevent and inhibit the inflammatory processes related to COVID-19.
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) first appeared in Wuhan, China, in December 2019 and then spread rapidly across China, causing large-scale respiratory infections in the population. So far, we are still unable to stop the transmission of SARS-CoV-2 and the Corona Virus Disease 2019 (COVID-19) continues to threaten human health, causing several deaths daily worldwide. Currently, there are several effective vaccines against SARS-CoV-2, such as mRNA-1273, ChAdOx1 nCoV-19, CoronaVac, and BNT162b2. Although the vaccination program prevented thousands of deaths from COVID-19, it was accompanied by new vaccine-related complications, including vaccine-induced immune thrombocytopenia and thrombosis (VITT) , Bell’s palsy , and hypersensitivity myocarditis . Therefore, exploration of more effective agents, with a lower incidence of side effects, against COVID-19 is still an urgent requirement. Recently, several novel therapeutic methods have been hypothesized for combating the COVID-19 crisis, including allium sativum derived carbon dots, silver nanoparticles, decorin, bilirubin nanomedicine, nanoceria, electric stimulation, and other nanotechnologies, which are justified to be potential theranostic agents for the management of COVID-19 [4,5,6,7,8,9,10,11]. Here, we hypothesized that Artemisia annua L. (A. annua) could be a promising theranostic agent for COVID-19, so we made many justifications for this hypothesis through data mining.
SARS-CoV-2 can infect the airway cells and result in excessive inflammation leading to a cytokine storm. Pathologically, platelet and endothelial dysfunction are the essential components of COVID-19 infection . Cytokine storm and the elevated levels of circulating cytokines are also associated with various other infectious and immune-mediated conditions as seen in COVID-19 .
A. annua, which is known to possess antiviral abilities [14, 15], is suggested by Haq et al.  to be of clinical importance in this epidemic and make a case for more elaborated clinical trials and experimental studies to explore its effects on the SARS-CoV-2 virus. Furthermore, some components of A. annua have been reported to have effective antiviral properties and immunosuppressive effects in vivo . Artemisinin, an antimalarial lactone derived from A. annua, is reported to decrease the infiltration of immunomodulatory cells and inflammatory cytokines in vivo . The use of Artemisinin against different respiratory diseases has also been investigated in lung cancer models and inflammatory-driven respiratory disorders . Isohamnetin, an active component of A. annua, an in vivo study showed it can inhibit SARS-CoV-2 spike pseudotyped virus entering ACE2h cells . Quercetin, belonging to the active component of A. annua, a comprehensive review summarized it seems to protect against SARS-CoV-2 through different mechanisms of action . A previous stiduy indicated that another active component of A. annua named kaempferol may potentially interact with the SARS-CoV-2 main protease 3CLpro . These studies above demonstrated the anti-SARS-CoV-2 activity of A. annua.
Given these findings, there is still a lack of direct evidence proving that A. annua may play an important role in COVID-19 management. Therefore, it is essential to investigate the potential active components of A. annua against COVID-19 and figure out their mechanisms of action. Network pharmacology, which holds a systematic and holistic view in understanding the nature of drugs, has been endorsed by various British pharmacologists as a mean to represent datasets and reveal the nature of the interactions between several nodes [23, 24]. These nodes represent small molecules like genes and proteins.
Our study aims at exploring the relevant roles of the potentially active components of A. annua for the treatment of COVID-19 using network pharmacology tools. Through the complex network analysis of drug-target-disease signaling pathway, we lay the foundation for further study on A. annua and A.annua-based drugs to treat COVID-19, and provide a platform to show the effectiveness of Chinese medicine against this epidemic. The process of this research is shown in Fig. 1.
Active compounds screening
The Traditional Chinese Medicines for Systems Pharmacology Database and Analysis Platform (TCMSP, available online: https://www.tcmsp-e.com), a pharmacology platform built for herbal medicines, was used to gather all the chemical components of A. annua . Potential active components were defined as ingredients with OB ≥ 30% and DL ≥ 0.18. We further identified effective targets of each compound by DrugBank Database (https://www.drugbank.ca) .
Compound-target network construction
Same targets generated by the same-compound prediction were screened and removed to avoid repetition and were then entered into UniProt Database (https://www.uniprot.org) . The search for the gene symbols of the active compounds in A. annua was refined by selecting species as ‘Homo sapiens’, so as to exclude non-human, irregular, and repetitive targets. Compounds and targets data were finally introduced into Cytoscape 3.7.2 software to build a compound-target (C-T) network.
COVID-19 related targets collection
COVID-19 related targets were obtained by GeneCards Database (https://www.genecards.org) . We searched for the keywords Novel Coronavirus Pneumonia, COVID-19, SARS-COV-2, and calculated the correlation between the targets and COVID-19 through the Gifts algorithm. Targets with a score of 30 were selected as the potential COVID-19 targets according to their correlation ranking and the standard gene names of the targets were also figured out. Additionally, OMIM (https://omim.org)  and TTD databases (http://db.idrblab.net/ttd)  were used to obtain more targets.
Compound-target-disease network and protein-protein interaction (PPI) construction
The compound-target-disease (C-T-D) network was constructed on the direct therapeutic targets of A. annua that acts on COVID-19. The common intersection for A. annua and COVID-19 was obtained by Venny online analysis tool , and the selected common targets were uploaded to String 11.0 Database (https://string-db.org, version 11.0) to construct the PPI network. The organismic selection was for ‘Homo sapiens’, and the minimum required interaction score was set to ‘medium confidence (>0.400)’. The topological importance of nodes in the C-T-D network was evaluated by the degree, betweenness centrality, and closeness centrality . A node was considered a hub when the degree of the node was more than the twofold average degree value . Subsequently, three topological parameters were computed to pick up big hubs, in which the degree, betweenness centrality, and closeness centrality were greater than the median value of all hubs in the network .
Gene ontology (GO) enrichment analysis
GO functional enrichment analysis, such as biological process (BP), cellular component (CC) and molecular function (MF) of target proteins of A. annua acting on COVID-19, was performed using David 6.8.0 Database (https://david.ncifcrf.gov) . An adjusted P-value of 0.05 was considered to identify the enriched terms.
Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis
KEGG Database (https://www.kegg.jp/kegg/) [36,37,38] was used to perform KEGG pathway enrichment analysis of target proteins of A. annua acting on COVID-19. KEGG Mapper queried the location of common targets between A.annua and SARS-CoV-2 and shared targets in the key pathway.
iGemdock molecular docking
The structures of the hub targets were downloaded from PDB (http://www.pdb.org)  and saved as PDB files, and the structures of the hub compounds were obtained from TCMSP and saved as mol2 files. Before docking, polar hydrogens were added to target proteins, and water molecules and original ligands were removed from the target proteins using Discovery Studio (Version 2016). Then, all the removed ligands were also saved as mol2 files for further docking. Further, iGEMDOCK (http://gemdock.life.nctu.edu.tw/dock, Version 2.1)  was used to complete Molecular docking simulation. The default parameters were set to ‘standard docking’ and ‘docking accuracy settings’. In general, the binding affinity is determined by the binding free energy, which is relatively low when the conformation of the binding is stable.
GlideXP molecular docking
The composite system of Isorhamnetin, Quercetin, Kaempferol, and each protein using XGlide module in the Schrödinger2018–1 software. The 3D structure of the small molecule was downloaded from Pubchem. For the protein structure, we downloaded it from the PDB database. The PDB IDs were 3 W55 (MAP1), 5 T89 (VEGFA), 6MXY (TP53), 6SFJ (MAP8), 6ZR5 (MAP14). Subsequently, the Protein Preparation Wizard in Schrödinger was used to optimize each protein separately, including removing non-ligand molecules and water molecules, adding hydrogen atoms, and using the OPLS2005 force field to optimize the structure remove intermolecular collisions. We finally docked the prepared files through the XGlide docking module, used a sitemap to predict the docking site, and the standard precision (SP) algorithm for molecular docking.
Molecular dynamics simulations
The small molecule-protein complex system obtained by docking was used as the initial structure to perform all-atom molecular dynamics simulations. The charge of the small molecule was calculated by the antechamber module and the Hartree–Fock (HF) SCF/6-31G* of the Gaussian 09 software [41, 42]. Small molecules and proteins were described using GAFF2 small molecule force field and ff14SB protein force field, respectively [43, 44]. The four systems all used the LEAP module to add hydrogen atoms to the system, add a truncated octahedral TIP3P solvent box five at a distance of 8 Å from the system , and add Na+/Cl− to the system to balance the system charge.
The molecular dynamics simulation was performed with AMBER 18 software . Before the simulation, energy optimization was performed on the system, including the steepest descent method with 2500 steps and the conjugate gradient method with 2500 steps. After the system energy optimization was completed, the system’s temperature was increased by 200 ps under a fixed volume and a constant heating rate so that the temperature of the system slowly rose from 0 K to 298.15 K. Under the condition that the system maintains a temperature of 298.15 K, a 500 ps NVT (isothermal and isopyknic) system simulation was performed to distribute the solvent molecules in the solvent box uniformly. Then, in the case of NPT (isothermal and isobaric), a 500 ps equilibrium simulation was performed on the entire system. Finally, under periodic boundary conditions, the four composite systems were subjected to 4 ns NPT ((isothermal and isobaric) system simulations. In the simulation, the cutoff distance of the non-bond was set to 10 Å, the Particle mesh Ewald (PME) method was used to calculate the long-range electrostatic interaction , the SHAKE method was used to limit the length of the hydrogen atom bond , and the Langevin algorithm was used for temperature control , the collision frequency γ was set to 2 ps− 1. The system pressure was 1 atm, the integration step was 2 fs, and the trajectory was saved every 4 ps for subsequent stability analysis and binding free energy calculation.
Binding free energy calculations
We calculated the free energy of binding between proteins and ligands in all systems using the MM-GBSA method [50,51,52,53]. In this study, the MD trajectory of 2–4 ns was used for calculation, and the specific formula is as follows:
In formula (1), ΔEinternal, ΔEVDW, and ΔEelec represents internal energy, van der Waals interaction, and electrostatic interaction, respectively. The internal energy includes bond energy (Ebond), angular energy (Eangle), and torsion energy (Etorsion), and ΔGGB and ΔGSA collectively referred to as solvation free energy (ΔGGB is the free energy of polar solvation, and ΔGSA is the free energy of non-polar solvation). We used the GB model (igb = 2) developed by researchers such as Nguyen  to calculate ΔGGB. ΔGSA was calculated based on the product of surface tension (γ) and solvent accessible surface area (SA), ΔGSA =0.0072 × ΔSASA . We ignored Entropy changes due to high computational resources and low precision [50, 51].
Active compounds selection
19 active compounds found in A. annua were screened by TCMSP platform under the standard conditions for a drug, particularly OB ≥ 30% and DL ≥ 0.18, including eupatin, isorhamnetin, sitosterol, etc. The OB and DL values for each ingredient are shown in Table 1.
Compound-target network and compound-target-disease network
We got 19 compound nodes and 208 target nodes (Fig. 2a). We obtained 871 COVID-19 targets by merging and deleting the duplicate values. From the 871 COVID-19 target genes retrieved, the C-T-D network comprised of 88 nodes and 229 edges. The details are presented in Fig. 2b.
Protein-protein interaction (PPI) construction
The screened active compounds of A. annua were intersected with the COVID-19 related targets, and the Venn analysis tool was used to obtain 71 common targets (Fig. 3a). The PPI network diagram contained 71 nodes and 1071 edges (Fig. 3b). 36 key nodes were selected as hubs according to the used standard (degree value > mean value) with the cutoff of twofold median value of degree (degree> 28.65). The nodes obtained were VEGFA, IL6, TP53, TNF, MAPK1, MAPK8, CASP3, MAPK14, IL10, CXCL8, CCL2, FOS, PTGS2, RELA, EGFR, IL1B, IL4, CASP8, PTEN, ICAM1, IL2, ESR1, BCL2L1, IFNG, PPARG, NOS3, HMOX1, SERPINE1, STAT1, TGFB1, NOS2, CXCL10, CRP, MPO, CDKN2A and AP. VEGFA was the node with the highest degree value as 60 (Fig. 3c). With a value of 57, IL6 stood as the second highest node. TP53, TNF, and MAPK1 had the same degree value of 54, IL10 had a value of 46, IL4 a value of 40, and IL2 a value of 39. The higher the node value of the target, the more important role it played in the network regulation and hence was likely to be a key target for A. annua to treat COVID-19.
GO enrichment analysis
GO function enrichment analysis was carried out with the help of DAVID database. 545 GO entries were obtained and 429 biological processes (BP). Mainly a negative regulation of transcription was observed. Other observations were related to response to heat, positive regulation of fibroblast proliferation, defense response to Gram-negative bacterium, cell growth, negative regulation of epithelial cell proliferation, response to mechanical stimulus, removal of superoxide radicals, PERK-mediated unfolded protein response, and a positive regulation of macrophage differentiation. GO analysis also highlighted 44 cellular components (CC), mainly extracellular exosome, intracellular, endoplasmic reticulum, endoplasmic reticulum membrane, protein complex, neuron projection, axon, lysosome, extracellular matrix, and platelet alpha granule lumen. 72 molecular functions (MF), mainly transcription factor activity, sequence-specific DNA binding, protein kinase activity, protein serine kinase activity, transcription regulatory region DNA binding, kinase activity, RNA polymerase II core promoter proximal region sequence-specific DNA binding, calmodulin binding, and protein complex binding, were seen. 58 GO entries were finally screened out, as shown in Fig. 4.
KEGG pathway enrichment analysis
To clarify the role of A. annua targets in the signaling pathway, analysis of KEGG pathway enrichment was performed using the KEGG database. The analysis results showed 221 pathways of A. annua related to the COVID-19. The first 20 KEGG pathways were finally screened out by gene count and shown in Fig. 5. Then the COVID-19 pathway and the location of SARS-CoV-2 targets and overlapping genes of enriched pathways are listed in Fig. 6.
iGemdock molecular docking analysis
The hubs namely quercetin, isorhamnetin and kaempferol, with above-average values for degree (> 28.65), betweenness (> 0.0128) and closeness (> 0.619), were simultaneously regarded as the hub compounds. The details are presented in Table 2.
We considered the median value of all pairs of binding energy as a strong binding efficacy of less than − 99.49 KJ/mol. The docking simulation results indicated that isorhamnetin and quercetin exhibited excellent binding affinity, and the details on the binding energy are presented in Fig. 7. The binding energy between the removed ligands and their target proteins was calculated in order to further investigate the binding affinity. Our results showed that the binding affinity of isorhamnetin-TP53 and quercetin-TP53 were far more superior than that of the removed original ligands and targets. Thus, these results indicated that isorhamnetin and quercetin were the potential therapeutic active compounds and TP53 was the potential therapetic target.
GlideXP molecular docking analysis
To further enhance the study’s credibility, we used the GlideXP high-precision docking algorithm to obtain the combination of the above seven hub targets and the three hub active components. The active site parameters of each target were calculated and are listed in Table 3. The lower docking affinity reflects the stronger binding ability between components and targets. As shown in Table 3 and Fig. 8, The combination with a value lower than − 8.00 kcal/mol binds more tightly, suggesting that MPK1, MPK8, and MPK14 are essential potential targets. Therefore, A.annua may improve COVID-19 by regulating the activity of these proteins. Then we used molecular dynamics simulation to study more in-depth.
MD simulation analysis
Based on the binding modes of three natural products and four potential targets, we tried to further explore the stability of these 21 complexes through molecular dynamics simulation technology to clarify which complexes are more reliable. The results are shown in Fig. 9, which shows the root mean square fluctuations (RMSD) of 21 kinds of compounds in the MD process and the stability of the compound with the fluctuations. In contrast, most of the stablest fluctuations are the complex system composed of MAPK1, consistent with those mentioned above GlideXP molecular docking results. Figure 9 shows the fluctuation of the ligand during the complex dynamic simulation. The more violent the fluctuation, the more unstable it is. During dynamic simulation, the RMSD of each molecule in the corresponding protein changed with the simulation time to show the stability of the numerical reaction molecule at the active site. The lower the fluctuation, the more stable the binding.
In addition to investigating the stability between the three compounds and the four targets, we also adopted the high-precision binding free energy calculation method MM-GBSA to study the binding free energy of these complex systems more directly.
Table 4 and Fig. 10 show the specific results. We obtained the binding free energy of 21 systems through MM-GBSA calculation and decomposed the energy term. Consistent with the above calculations, the highest binding free energy combination was the hybrid system of MAPK1/Quercetin, and its binding free energy was − 36.7 kcal/mol. The higher combination was TP53/Kaempferol, and its binding free energy was − 34.33 kcal/mol. Our analysis showed that these high binding free energy combinations were mainly contributed by van der Waals forces and electrostatic interactions. In summary, we believe that the MAPK1/Quercetin, TP53/Kaempferol complex has a high degree of credibility and is necessary for further experimental verification.
Binding mode analysis
Complexes such as MAPK8/Isorhamnetin, MAPK8/Quercetin, MAPK8/Kaempferol, and MAPK14/Kaempferol were good interactions during GlideXP molecular docking binding free energy. They all displayed negative binding free energy values excluded from further analysis. MAPK1/Quercetin has good interaction during GlideXP molecular docking and favorable binding free energy. Although Complex TP53/Kaempferol did not display good interaction during GlideXP molecular docking, it was found to have good interaction during iGemdockimolecular docking and showed good binding free energy values. Therefore, we selected MAPK1/Quercetin and TP53/Kaempferol complexes to verify their binding mode. Figure 11 shows that the binding mode of MAPK1/Quercetin and TP53/Kaempferol complexes combine well. The formation of hydrogen bonds between quercetin and the active site residues of MAPK1 involved residues ASP106, MET108, LYS114, ASN154, and ASP167 (Fig. 11a). The formation of hydrogen bonds between kaempferol and the active site residues of TP53 involved residues MET1584 and SER1503 (Fig. 11b). Hence, our findings confirmed that A.annua might treat COVID-19 by regulating the activity of MAPK1 and TP53.
Besides A. annua’s wide usage as an antimalarial drug, there are many clinical and experimental studies showing that A. annua also has antiviral and anti-inflammatory properties. Mathieu Gendrot et al. conducted an inhibition test of SARS-CoV-2 replication in Vero E6 cells with 5 artemisinin-based combination therapies (ACTs) in vitro. In their study, mefloquine-artesunate exerted the highest antiviral activity with 72.1–18.3% inhibition at the expected maximum blood concentration (Cmax) for each ACT drug, at doses commonly administered in malaria treatment. Additionally, all the other combinations include artesunate-amodiaquine, artemether-lumefantrine, artesunate-pyronaridine, or dihydroartemisinin-piperaquine, showed antiviral inhibition in the same ranges (27.1 to 34.1%) . Therefore, we speculated that A. annua may have a significant effect on COVID-19 as well. In our study, we elucidated the relationship between A. annua and COVID-19 using a network pharmacology strategy. We screened 19 potential active ingredients and 208 potential targets of A. annua and obtained 71 common targets for A. annua and COVID-19. We constructed networks such as compound-target network, PPI network of A. annua and COVID-19, and performed GO and KEGG enrichment analysis to comprehensively understand the mechanisms behind A. annua’s capabilities in treating COVID-19.
In addition to the overall efficacy in clearing heat and detoxifying properties, each compound of A. annua can also play different roles. As shown by the C-T network analysis, quercetin, isorhamnetin and kaempferol compounds play a crucial role in manifesting treatment properties of A. annua in treating COVID-19. Quercetin can also prevent the virus from sticking to cells by interacting with proteins on the surface of the virus  and by inhibiting inflammatory cytokines . Isorhamnetin owns strong antiviral potency, which can block cytoplasmic lysosome acidification and reduce virus-induced ROS generation . Kaempferol is known to significantly suppress cell-autonomous immunity through down-regulation of p38 and JNK . Interestingly, its derivatives are also considered to have antiviral properties against the 3a channel protein of coronavirus, a channel that is related to the viral release mechanism . Artemisinin and its derivations were beneficial for the management of viral infection, and the antiviral effect is associated with enhanced type I interferon response of the host [62, 63]. Thus Quercetin, isorhamnetin, kaempferol, and artemisinin could be the main potentially active compounds of A. annua that are employed in the management of COVID-19 by inhibiting inflammatory mediators.
PPI network analysis showed that the degrees and betweenness centrality of 36 key genes were greater than the average. The further topology analysis of the C-T-D network revealed that 7 key target proteins include key genes for the pathological mechanisms of COVID-19, mainly the immune and the blood coagulation mechanism. For instance, vascular endothelial growth factor-A (VEGF-A) is a crucial regulator of angiogenesis, a process that forms new blood vessels from pre-existing vessels [64, 65]. Clinical data showed that several COVID-19 patients are presented with increasing levels of angiogenesis and endotheliopathy markers [66, 67]. Additionally, brain inflammation caused by SARS-CoV-2 was considered to be related to VEGFA, which can facilitate the accumulation of inflammatory cells and regulate the angiopoietins II level . TNF is an intensely studied proinflammatory cytokine, which can trigger either cell survival or cell death by regulating a complex inflammatory network [69, 70], and its hyperproduction is related to the pathogenesis of COVID-19 as well. The caspase family regulates crucial biological functions, such as cell death in apoptosis and pyroptosis, as well as performs non-cell death functions in inflammation, dendrite trimming, and cell differentiation . A spatially restricted activation of caspase-3 occurs in maturing megakaryocytes to promote proplatelet maturation and platelet shedding in the bloodstream .
TP53 is a transcription factor that functions towards the regulation of important cellular activities, such as cell cycle, senescence, and apoptosis, which can suppress inflammation in a plethora of human tissues  and thus is often mutated in certain malignancies. The increased MAPK activation may influence the release of many proinflammatory cytokines  and induce platelet aggregation in COVID-19 patients . Based on the above details, we find that A. annua can play a therapeutic role in the entire process of COVID-19 management thus, meeting the need for the systematic treatment of COVID-19 and futuristic studies, it would be worthwhile to explore the role of VEGFA, CASP3, TP53, TNF, MAPK1, MAPK8 and MAPK14 in the mechanisms manifested by A. annua. Through molecular docking simulation and computing binding free energy, we validated that 3 hub compounds quercetin, isorhamnetin, and kaempferol, particularly quercetin and isorhamnetin, had an excellent binding affinity with the 7 hub targets. The MD and MM-GBSA calculation data further indicated that MPK1/Quercetin and TP53/Kaempferol possessed the highest binding free energy, which demonstrated the repurposing possibilities of quercetin and kaempferol based on their binding activity with multiple COVID-19 targets, and supported their ability to function as anti-SARS-CoV-2.
The GO enrichment analysis revealed that the treatment of COVID-19 with A. annua was mainly related to BP, CC, and MF. BP includes negative regulation of transcription, response to heat, and positive regulation of fibroblast proliferation. CC includes extracellular exosome, intracellular, and endoplasmic reticulum. MF includes transcription factor activity, sequence-specific DNA binding, and protein kinase activity. The KEGG pathways of A. annua are mainly related to the roles of immune regulation. For instance, Gonadotropin-releasing hormone (GnRH) analogs can lead to pro-inflammatory changes in T lymphocytes . TRP channels regulate fundamental biological processes throughout the body, and the dysfunction of these channels has been causally linked to a number of disease states  that regulate the innate immune cell function in lung inflammation . The airway wall is an essential controller of inflammatory, immune and regenerative responses to viruses, and Asthma is a T lymphocyte-controlled disease of the airway wall caused by inflammation, thus a strong link between Asthma and infection with coronaviruses is established [79, 80]. Based on the GO and KEGG results, as mentioned above, it could be postulated that A. annua acts on COVID-19 mainly by regulating inflammatory responses, transcription, and proliferation.
Antimalarial drugs, for which lung concentration data are available, are found to be 10 to 160-fold more concentrated in the lungs than in blood . Artemisinin-type compound possesses an antiviral ability . Through the clinical observation, we primarily speculate that A. annua related drugs such as artemisinin piperaquine have considerable command on controlling the processes leading to lung inflammation.
The use of A. annua for COVID-19 could be mainly linked to its anti-inflammatory and immune regulatory properties. Its multiple active compounds aim for multiple targets, which were mapped to different pathways embodied in the complex integrated network mechanisms of multi-component, multi-targeted, and multi-channel regulated treatment of COVID-19. Although our study provides theoretical support and scientific evidence for A. annua’s capabilities in fighting against COVID-19, it has some limitations in terms of lacking specific experimental verification. Thus, our study puts foundation for further screening of these compounds and verification of the observations made with animal or cell experiments to clarify the main regulatory targets of A. annua.
In summary, in our study, we have used a network pharmacology strategy to predict the main active compounds and key targets of A. annua for the treatment of COVID-19 and speculated the potential mechanisms from multiple approaches and perspectives.
Our study initially found A. annua to have anti-inflammatory and immune regulatory properties. Further, we analyzed the potential active compounds and targets of A. annua for COVID-19 using the network pharmacology methods. The screened results show that A. annua can prevent and treat COVID-19 through multiple components, targets, and pathways that need further validation through in vitro or in vivo studies.
Availability of data and materials
The datasets generated and analysed during the current study are available in the Zenodo repository (https://zenodo.org/record/5565195#.YWXPgvkzY2x).
- A. annua :
Artemisia annua L
Artemisinin-based combination therapies
Corona Virus Disease 2019
Kyoto Encyclopedia of Genes and Genomes
Molecular dynamics simulation studies
Root mean square deviation
Severe acute respiratory syndrome coronavirus 2
Traditional Chinese Medicines for Systems Pharmacology Database and Analysis Platform
Transient receptor potential
COVID-19 rapid guideline: vaccine-induced immune thrombocytopenia and thrombosis (VITT). NICE guideline[NG200]. 2021. https://www.nice.org.uk/guidance/ng20.
Wan EYF, Chui CSL, Lai FTT, et al. Bell's palsy following vaccination with mRNA (BNT162b2) and inactivated (CoronaVac) SARS-CoV-2 vaccines: a case series and nested case-control study. Lancet Infect Dis. 2021. https://doi.org/10.1016/S1473-3099(21)00451-5.
Kounis NG, Mplani V, Koniari I, Velissaris D. Hypersensitivity myocarditis and COVID-19 vaccines. Kardiol Pol. 2021. https://doi.org/10.33963/KP.a2021.0166.
Kalkal A, Allawadhi P, Pradhan R, Khurana A, Bharani KK, Packirisamy G. Allium sativum derived carbon dots as a potential theranostic agent to combat the COVID-19 crisis. Sens Int. 2021;2:100102. https://doi.org/10.1016/j.sintl.2021.100102.
Allawadhi P, Singh V, Khurana A, et al. Silver nanoparticle based multifunctional approach for combating COVID-19. Sens Int. 2021;2:100101. https://doi.org/10.1016/j.sintl.2021.100101.
Allawadhi P, Singh V, Khurana I, et al. Decorin as a possible strategy for the amelioration of COVID-19. Med Hypotheses. 2021;152:110612. https://doi.org/10.1016/j.mehy.2021.110612.
Singh V, Allawadhi P, Khurana A, Banothu AK, Bharani KK. Critical neurological features of COVID-19: role of imaging methods and biosensors for effective diagnosis. Sens Int. 2021;2:100098. https://doi.org/10.1016/j.sintl.2021.100098.
Khurana I, Allawadhi P, Khurana A, et al. Can bilirubin nanomedicine become a hope for the management of COVID-19? Med Hypotheses. 2021;149:110534. https://doi.org/10.1016/j.mehy.2021.110534.
Khurana A, Allawadhi P, Khurana I, et al. Role of nanotechnology behind the success of mRNA vaccines for COVID-19. Nano Today. 2021;38:101142. https://doi.org/10.1016/j.nantod.2021.101142.
Allawadhi P, Khurana A, Allwadhi S, Joshi K, Packirisamy G, Bharani KK. Nanoceria as a possible agent for the management of COVID-19. Nano Today. 2020;35:100982. https://doi.org/10.1016/j.nantod.2020.100982.
Allawadhi P, Khurana A, Allwadhi S, et al. Potential of electric stimulation for the management of COVID-19. Med Hypotheses. 2020;144:110259. https://doi.org/10.1016/j.mehy.2020.110259.
Gu SX, Tyagi T, Jain K, et al. Thrombocytopathy and endotheliopathy: crucial contributors to COVID-19 thromboinflammation. Nat Rev Cardiol. 2021;18:194–209. https://doi.org/10.1038/s41569-020-00469-1.
Mangalmurti N, Hunter CA. Cytokine storms: understanding COVID-19. Immunity. 2020;53:19–25. https://doi.org/10.1016/j.immuni.2020.06.017.
Adhikari B, Marasini BP, Rayamajhee B, et al. Potential roles of medicinal plants for the treatment of viral diseases focusing on COVID-19: a review. Phytother Res. 2021;35:1298–312. https://doi.org/10.1002/ptr.6893.
Liu X, Cao J, Huang G, Zhao Q, Shen J. Biological activities of artemisinin derivatives beyond malaria. Curr Top Med Chem. 2019;19:205–22. https://doi.org/10.2174/1568026619666190122144217.
Haq FU, Roman M, Ahmad K, et al. Artemisia annua: trials are needed for COVID-19. Phytother Res. 2020;34:2423–4. https://doi.org/10.1002/ptr.6733.
Ho WE, Peh HY, Chan TK, Wong WS. Artemisinins: pharmacological actions beyond anti-malarial. Pharmacol Ther. 2014;142:126–39. https://doi.org/10.1016/j.pharmthera.2013.12.001.
Lee AS, Hur HJ, Sung MJ. The effect of artemisinin on inflammation-associated Lymphangiogenesis in experimental acute colitis. Int J Mol Sci. 2020;21. https://doi.org/10.3390/ijms21218068.
Cheong DHJ, Tan DWS, Wong FWS, Tran T. Anti-malarial drug, artemisinin and its derivatives for the treatment of respiratory diseases. Pharmacol Res. 2020;158:104901. https://doi.org/10.1016/j.phrs.2020.104901.
Zhan Y, Ta W, Tang W, et al. Potential antiviral activity of isorhamnetin against SARS-CoV-2 spike pseudotyped virus in vitro. Drug Dev Res. 2021. https://doi.org/10.1002/ddr.21815.
Di Petrillo A, Orru G, Fais A, Fantini MC. Quercetin and its derivates as antiviral potentials: a comprehensive review. Phytother Res. 2021. https://doi.org/10.1002/ptr.7309.
Khan A, Heng W, Wang Y, et al. In silico and in vitro evaluation of kaempferol as a potential inhibitor of the SARS-CoV-2 main protease (3CLpro). Phytother Res. 2021;35:2841–5. https://doi.org/10.1002/ptr.6998.
Hopkins AL. Network pharmacology. Nat Biotechnol. 2007;25:1110–1. https://doi.org/10.1038/nbt1007-1110.
Berger SI, Iyengar R. Network analyses in systems pharmacology. Bioinformatics. 2009;25:2466–72. https://doi.org/10.1093/bioinformatics/btp465.
Ru J, Li P, Wang J, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014;6:13. https://doi.org/10.1186/1758-2946-6-13.
Wishart DS, Feunang YD, Guo AC, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018;46:D1074–82. https://doi.org/10.1093/nar/gkx1037.
UniProt C. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 2019;47:D506–15. https://doi.org/10.1093/nar/gky1049.
Stelzer G, Rosen N, Plaschkes I, et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr Protoc Bioinformatics. 2016;54:1 30 31–31 30 33. https://doi.org/10.1002/cpbi.5.
Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. OMIM.org: online Mendelian inheritance in man (OMIM(R)), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2015;43(D789–798). https://doi.org/10.1093/nar/gku1205.
Wang Y, Zhang S, Li F, et al. Therapeutic target database 2020: enriched resource for facilitating research and early development of targeted therapeutics. Nucleic Acids Res. 2020;48:D1031–41. https://doi.org/10.1093/nar/gkz981.
Oliveros JC. VENNY. An interactive tool for comparing lists with Venn Diagrams; 2007.
Ashtiani M, Salehzadeh-Yazdi A, Razaghi-Moghadam Z, et al. A systematic survey of centrality measures for protein-protein interaction networks. BMC Syst Biol. 2018;12:80. https://doi.org/10.1186/s12918-018-0598-2.
Zhang Y, Li Y, Mao X, et al. Thyroid hormone synthesis: a potential target of a Chinese herbal formula Haizao Yuhu decoction acting on iodine-deficient goiter. Oncotarget. 2016;7:51699–712. https://doi.org/10.18632/oncotarget.10329.
Yu G, Zhang Y, Ren W, et al. Network pharmacology-based identification of key pharmacological pathways of Yin-Huang-Qing-Fei capsule acting on chronic bronchitis. Int J Chron Obstruct Pulmon Dis. 2017;12:85–94. https://doi.org/10.2147/COPD.S121079.
Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13. https://doi.org/10.1093/nar/gkn923.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. https://doi.org/10.1093/nar/28.1.27.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51. https://doi.org/10.1002/pro.3715.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51. https://doi.org/10.1093/nar/gkaa970.
Burley SK, Berman HM, Bhikadiya C, et al. RCSB protein data Bank: biological macromolecular structures enabling research and education in fundamental biology, biomedicine, biotechnology and energy. Nucleic Acids Res. 2019;47:D464–74. https://doi.org/10.1093/nar/gky1004.
Yang JM, Chen CC. GEMDOCK: a generic evolutionary method for molecular docking. Proteins. 2004;55:288–304. https://doi.org/10.1002/prot.20035.
Frisch M, Trucks G, Schlegel H, et al. Gaussian 09 (Revision D.01). Wallingford CT: Gaussian Inc.; 2009.
Wang J, Wang W, Kollman PA, Case DA. Antechamber, an accessory software package for molecular mechanical calculations. J Am Chem Soc. 2001;222:U403.
Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA. Development and testing of a general amber force field. J Comput Chem. 2004;25:1157–74. https://doi.org/10.1002/jcc.20035.
Maier JA, Martinez C, Kasavajhala K, Wickstrom L, Hauser KE, Simmerling C. ff14SB: improving the accuracy of protein side chain and backbone parameters from ff99SB. J Chem Theory Comput. 2015;11:3696–713. https://doi.org/10.1021/acs.jctc.5b00255.
Mark P, Nilsson L. Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K. J Phys Chem A. 2001;105(43):9954–60. https://doi.org/10.1021/jp003020w.
Salomon-Ferrer R, Case DA, Walker RC. An overview of the Amber biomolecular simulation package. Wiley Interdiscip Rev Comput Mo Sci. 2013;3(2):198–210. https://doi.org/10.1002/wcms.1121.
Sagui C, Darden TA. Molecular dynamics simulations of biomolecules: long-range electrostatic effects. Annu Rev Biophys Biomol Struct. 1999;28:155–79. https://doi.org/10.1146/annurev.biophys.28.1.155.
Kräutler V, Van Gunsteren WF, Hünenberger PH. A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations. J Comput Chem. 2001;22(5):501–8.
Larini L, Mannella R, Leporini D. Langevin stabilization of molecular-dynamics simulations of polymers by means of quasisymplectic algorithms. J Chem Phys. 2007;126:104101. https://doi.org/10.1063/1.2464095.
Chen Y, Zheng Y, Fong P, Mao S, Wang Q. The application of the MM/GBSA method in the binding pose prediction of FGFR inhibitors. Phys Chem Chem Phys. 2020;22:9656–63. https://doi.org/10.1039/d0cp00831a.
Hou T, Wang J, Li Y, Wang W. Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations. J Chem Inf Model. 2011;51:69–82. https://doi.org/10.1021/ci100275a.
Genheden S, Ryde U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin Drug Discov. 2015;10:449–61. https://doi.org/10.1517/17460441.2015.1032936.
Rastelli G, Del Rio A, Degliesposti G, Sgobba M. Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA. J Comput Chem. 2010;31:797–810. https://doi.org/10.1002/jcc.21372.
Nguyen H, Roe DR, Simmerling C. Improved generalized born solvent model parameters for protein simulations. J Chem Theory Comput. 2013;9:2020–34. https://doi.org/10.1021/ct3010485.
Weiser J, Shenkin PS, Still WC. Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO). J Comput Chem. 1999;20(2):217–30. https://doi.org/10.1002/(SICI)1096-987X(19990130)20:2<217::AID-JCC4>3.0.CO;2-A.
Gendrot M, Duflot I, Boxberger M, et al. Antimalarial artemisinin-based combination therapies (ACT) and COVID-19 in Africa: in vitro inhibition of SARS-CoV-2 replication by mefloquine-artesunate. Int J Infect Dis. 2020;99:437–40. https://doi.org/10.1016/j.ijid.2020.08.032.
Lopes BRP, da Costa MF, Genova Ribeiro A, et al. Quercetin pentaacetate inhibits in vitro human respiratory syncytial virus adhesion. Virus Res. 2020;276:197805. https://doi.org/10.1016/j.virusres.2019.197805.
Jasso-Miranda C, Herrera-Camacho I, Flores-Mendoza LK, et al. Antiviral and immunomodulatory effects of polyphenols on macrophages infected with dengue virus serotypes 2 and 3 enhanced or not with antibodies. Infect Drug Resist. 2019;12:1833–52. https://doi.org/10.2147/IDR.S210890.
Abdal Dayem A, Choi HY, Kim YB, Cho SG. Antiviral effect of methylated flavonol isorhamnetin against influenza. PLoS One. 2015;10:e0121610. https://doi.org/10.1371/journal.pone.0121610.
Dong W, Wei X, Zhang F, et al. A dual character of flavonoids in influenza a virus replication and spread through modulating cell-autonomous immunity by MAPK signaling pathways. Sci Rep. 2014;4:7237. https://doi.org/10.1038/srep07237.
Schwarz S, Sauter D, Wang K, et al. Kaempferol derivatives as antiviral drugs against the 3a channel protein of coronavirus. Planta Med. 2014;80:177–82. https://doi.org/10.1055/s-0033-1360277.
Wang X, Zheng B, Ashraf U, et al. Artemisinin inhibits the replication of flaviviruses by promoting the type I interferon production. Antivir Res. 2020;179:104810. https://doi.org/10.1016/j.antiviral.2020.104810.
D'Alessandro S, Scaccabarozzi D, Signorini L, et al. The Use of Antimalarial Drugs against Viral Infection. Microorganisms. 2020;8(1):85. https://doi.org/10.3390/microorganisms8010085.
Folkman J, Klagsbrun M. Angiogenic factors. Science. 1987;235:442–7. https://doi.org/10.1126/science.2432664.
Ferrara N, Gerber HP, LeCouter J. The biology of VEGF and its receptors. Nat Med. 2003;9:669–76. https://doi.org/10.1038/nm0603-669.
Meizlish ML, Pine AB, Goshua G, et al. Circulating markers of angiogenesis and Endotheliopathy in COVID-19. medRxiv. 2020. https://doi.org/10.1101/2020.06.29.20140376.
Ackermann M, Verleden SE, Kuehnel M, et al. Pulmonary vascular Endothelialitis, thrombosis, and angiogenesis in Covid-19. N Engl J Med. 2020;383:120–8. https://doi.org/10.1056/NEJMoa2015432.
Yin XX, Zheng XR, Peng W, Wu ML, Mao XY. Vascular endothelial growth factor (VEGF) as a vital target for brain inflammation during the COVID-19 outbreak. ACS Chem Neurosci. 2020;11:1704–5. https://doi.org/10.1021/acschemneuro.0c00294.
Blaser H, Dostert C, Mak TW, Brenner D. TNF and ROS crosstalk in inflammation. Trends Cell Biol. 2016;26:249–61. https://doi.org/10.1016/j.tcb.2015.12.002.
Brenner D, Blaser H, Mak TW. Regulation of tumour necrosis factor signalling: live or let die. Nat Rev Immunol. 2015;15:362–74. https://doi.org/10.1038/nri3834.
Julien O, Wells JA. Caspases and their substrates. Cell Death Differ. 2017;24:1380–9. https://doi.org/10.1038/cdd.2017.44.
Shalini S, Dorstyn L, Dawar S, Kumar S. Old, new and emerging functions of caspases. Cell Death Differ. 2015;22:526–39. https://doi.org/10.1038/cdd.2014.216.
Barabutis N, Schally AV, Siejka A. P53, GHRH, inflammation and cancer. EBioMedicine. 2018;37:557–62. https://doi.org/10.1016/j.ebiom.2018.10.034.
Grimes JM, Grimes KV. p38 MAPK inhibition: a promising therapeutic approach for COVID-19. J Mol Cell Cardiol. 2020;144:63–5. https://doi.org/10.1016/j.yjmcc.2020.05.007.
Manne BK, Denorme F, Middleton EA, et al. Platelet gene expression and function in patients with COVID-19. Blood. 2020;136:1317–29. https://doi.org/10.1182/blood.2020007214.
Sung N, Salazar Garcia MD, Dambaeva S, et al. Gonadotropin-releasing hormone analogues lead to pro-inflammatory changes in T lymphocytes. Am J Reprod Immunol. 2016;76:50–8. https://doi.org/10.1111/aji.12522.
Gotz V, Qiao S, Beck A, Boehm U. Transient receptor potential (TRP) channel function in the reproductive axis. Cell Calcium. 2017;67:138–47. https://doi.org/10.1016/j.ceca.2017.04.004.
Scheraga RG, Southern BD, Grove LM, Olman MA. The role of TRPV4 in regulating innate immune cell function in lung inflammation. Front Immunol. 2020;11:1211. https://doi.org/10.3389/fimmu.2020.01211.
Edwards MR, Bartlett NW, Hussell T, Openshaw P, Johnston SL. The microbiology of asthma. Nat Rev Microbiol. 2012;10:459–71. https://doi.org/10.1038/nrmicro2801.
Lambrecht BN, Hammad H. The airway epithelium in asthma. Nat Med. 2012;18:684–92. https://doi.org/10.1038/nm.2737.
Efferth T. Beyond malaria. The inhibition of viruses by artemisinin-type compounds. Biotechnol Adv. 2018;36:1730–7. https://doi.org/10.1016/j.biotechadv.2018.01.001.
This work was supported by Natural Science Foundation of China [Grant Number 81873218 and 82074301], Science and Technology Project of Guangdong Province [Grant Number 2020A0505090009 and 2021A0505030060], Science and Technology Project of Guangzhou City of Guangdong Province [2022-06-01-11-3016-0028], Special Research project on COVID-19 Prevention and Control of The Education Department of Guangdong Province [Grant Number 2020KZDZX1055].
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Tang, Y., Li, X., Yuan, Y. et al. Network pharmacology-based predictions of active components and pharmacological mechanisms of Artemisia annua L. for the treatment of the novel Corona virus disease 2019 (COVID-19). BMC Complement Med Ther 22, 56 (2022). https://doi.org/10.1186/s12906-022-03523-2
- Artemisia annua L.
- Network pharmacology
- Molecular docking
- Molecule dynamics