Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 23 September 2022
Sec. Inflammation Pharmacology

Network pharmacology and computer-aided drug design to explored potential targets of Lianhua Qingwen and Qingfei Paidu decoction for COVID-19

Liyuan LiLiyuan Li1Xiaoying WangXiaoying Wang1Xiao GuoXiao Guo1Yikun LiYikun Li1Qiuhang Song,
Qiuhang Song1,2*Aiying Li,,
Aiying Li1,2,3*
  • 1College of Basic Medicine, Hebei University of Chinese Medicine, Shijiazhuang, China
  • 2Hebei Key Laboratory of Chinese Medicine Research on Cardio-Cerebrovascular Disease, Hebei University of Chinese Medicine, Shijiazhuang, China
  • 3Hebei Higher Education Institute Applied Technology Research Center on TCM Formula Preparation, Shijiazhuang, China

Coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2, has spread globally, affecting people’s lives worldwide and hindering global development. Traditional Chinese Medicine (TCM) plays a unique role in preventing and treating COVID-19. Representative prescriptions for the COVID-19 treatment, Lianhua Qingwen (LHQW) and Qingfei Paidu Decoction (QFPD), effectively alleviate COVID-19 symptoms, delaying its progression and preventing its occurrence. Despite the extensive similarity in their therapeutic effects, the mechanisms and advantages of LHQW and QFPD in in treating mild-to-moderate COVID-19 remain elusive. To characterize the mechanisms of LHQW and QFPD in treating COVID-19, we used integrated network pharmacology and system biology to compare the LHQW and QFPD components, active compounds and their targets in Homo sapiens. LHQW and QFPD comprise 196 and 310 active compounds, some of which have identical targets. These targets are enriched in pathways associated with inflammation, immunity, apoptosis, oxidative stress, etc. However, the two TCM formulas also have specific active compounds and targets. In LHQW, arctiin, corymbosin, and aloe-emodin target neurological disease-related genes (GRM1 and GRM5), whereas in QFPD, isofucosterol, baicalein, nobiletin, oroxylin A, epiberberine, and piperlonguminine target immunity- and inflammation-related genes (mTOR and PLA2G4A). Our findings indicate that LHQW may be suitable for treating mild-to-moderate COVID-19 with nervous system symptoms. Moreover, QFPD may effectively regulate oxidative stress damage and inflammatory symptoms induced by SARS-CoV-2. These findings may provide references for the clinical application of LHQW and QFPD.

Introduction

Coronavirus disease 2019 (COVID-19) was first reported at the end of 2019. It then spread into a severe global public health pandemic, causing irreparable harm to human health. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative pathogen of COVID-19. Patients with mild-to-moderate COVID-19 often present cough, fever, muscle pain, and diarrhea, while patients with severe disease often present a strong inflammatory response and multiple organ failure. In fact, some patients exhibit central nervous system dysfunction and olfactory sensory disorder (Asadi-Pooya & Simani, 2020; Havervall et al., 2021; Nguyen et al., 2022). Moreover, recovered patients often experience sequelae, such as dyspnea and fatigue (Carfì et al., 2020). Various therapeutic strategies have been developed for treating COVID-19, including those derived from Traditional Chinese Medicine (TCM) (Salian et al., 2021). Among these TCM therapies, Lianhua Qingwen (LHQW) and Qingfei Paidu decoction (QFPD) in “Three Chinese Patent Medicines and Three TCM Prescriptions” were commonly used with high efficacy (Zhang S. et al., 2022).

Moreover, the components in LHQW capsules, namely forsythia, honeysuckle, isatis root, and herbal houttuynia, have properties that help heat-clearing, detoxification, and lung ventilation, which are significant in the treatment of influenza caused by the influenza A virus infection (Li L. C. et al., 2020; Gao et al., 2020; Zhuang et al., 2020; Chu et al., 2021). In addition, QFPD comprises four classic Chinese medicine prescriptions, namely Shegan Mahuang decoction, Maxing Shigan decoction, Xiaochaihu decoction, and Wuling powder. QFPD contains 21 main components, namely Belamcandae rhizoma, Alisma orientale, and Polyporus umbellatus, among others.

QFPD is widely used throughout China (Liu et al., 2020) for treating COVID-19 and has demonstrated efficacy in preventing disease progression, improving clinical symptoms, alleviating lung lesions, and reducing mortality rates, without inducing severe adverse effects (Zhang L. et al., 2022). Additionally, QFPD plays roles in inflammation, immunoregulation, neuroprotection, and lung injury reduction, among other physiological processes (Chen et al., 2020). Although both LHQW and QFPD can be used to treat mild-to-moderate COVID-19, their components are distinct, indicating that the similarity and specificity of the biological mechanisms underlying the effects of the two formulations warrant further investigation using contemporary platforms.

Previous studies demonstrated that patients with mild-and-moderate COVID-19 often show cough, inflammation, immune response, nervous system dysfunction, etc. To reveal the mechanism of LHQW and QFPD in treating these symptoms, we compared the LHQW and QFPD components, active compounds, and their targets in Homo sapiens. We found that LHQW and QFPD contain 196 and 310 active compounds respectively, they have common and specific targets. The common targets are enriched in pathways associated with inflammation, immunity, apoptosis, and oxidative stress. The special active compounds (Arctiin, Corymbosin, and Aloe-emodin) of LHQW target neurological disease-related genes (GRM1 and GRM5), and the special active compounds (isofucosterol, baicalein, nobiletin, oroxylin A, epiberberine, and piperlonguminine) of QFPD target immunity- and inflammation-related genes (mTOR and PLA2G4A). This finding suggested that LHQW may improve mild-to-moderate COVID-19 with nervous system symptoms, whereas QFPD may potentially regulate oxidative stress damage and inflammatory symptoms. These findings may provide potential references for the clinical application of LHQW and QFPD. Collectively, our study provides theoretical support for the effective administration of combinatorial TCM and allopathic medicine for the treatment of COVID-19.

Materials and methods

Selection and analysis of LHQW and QFPD active compounds

The composition, channel tropism, properties, and flavors of all LHQW and QFPD components were obtained from the Chinese Pharmacopoeia (2015 edition) (Commission, 2015). The active compounds with oral bioavailability (OB) ≥ 30% and drug-likeness (DL) ≥ 0.18 were derived from the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP, http://tcmspw.com/tcmsp.php) (Ru et al., 2014).

Acquisition of common and specific active compounds targets in LHQW and QFPD

Representative genes targeted by LHQW and QFPD active compounds were identified from the TCMSP database. To compare the similarity and specificity of LHQW and QFPD, we constructed a Venn diagram of target genes to determine their common and specific targets.

Construction of a protein-protein interaction (PPI) network and pathway enrichment analysis

The STRING (https://string-db.org/) (Szklarczyk et al., 2021) database was used to construct a PPI network of common targets, which was then topologically analyzed to identify hub targets. The targets in the PPI network were screened by “degree,” “between centrality,” and “closeness centrality,” and targets with scores greater than the median of the three parameters were selected to build a protein-protein interaction network (PIN). The molecular complex detection (MCODE) algorithm was used in the Cytoscape software to obtain the protein modules. In addition, the KOBAS (http://kobas.cbi.pku.edu.cn/) (Xie et al., 2011) database was used for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. The most enriched GO and KEGG pathways were extracted to analyze and compare the similarities between LHQW and QFPD.

Analyzing LHQW and QFPD specificity

The TCMSP database was used to identify active compounds that could target specific genes, as well as the TCM components that contained these active compounds. To determine the relationship between TCM components, active compounds, and specific targets, “TCM component–active compound–specific target” networks were constructed using Cytoscape software. The MCODE algorithm of the Cytoscape software was used to build the protein modules of specific genes in the networks.

Molecular docking analysis

The most significant target genes of LHQW and QFPD genes were selected, and the structures of proteins encoded by the target genes were obtained from the RCSB Protein Data Bank (RCSB PDB, https://www.rcsb.org/). The 3D structures of the active compounds were obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/). GOLD 5.2 and Molecular Operating Environment v2019.0102 (MOE) software were used for virtual docking, and ChemScore, ChemPLP and binding energy value were used in combination to evaluate the docking results. Finally, docking results were visualized using PyMol 1.8 and MOE.

Molecular dynamics simulation

After molecular docking, Amber 18 software was used to simulate the molecular dynamics of the targets and active compounds of the TCM preparation. Using Amber 99sb force parameters, the system was executed for 50 ns. Every 10 ps, trajectory data were saved with a 2-ps time step. The Gibbs free energy of the active compounds and proteins was calculated using the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method. The results were estimated using Gibbs free energy, root mean square deviation (RMSD), root mean square fluctuation (RMSF), and solvent-accessible surface area (SASA) metrics.

Results

Comparison of LHQW and QFPD composition, properties, flavors, and channel tropism

LHQW and QFPD components, properties, flavors, and channel tropism, are listed in Supplementary Table S1. Based on the TCM principle that most components in LHQW have “cold properties,” this formulation can be considered a “cold medicine”. In contrast, QFPD could be regarded as a ‘warm medicine’ as most of its components possess “warm properties.” Channel tropism is a critical TCM concept in which a drug can exhibit selective action in a specific part of the human body. Eight of the thirteen components in LHQW act through the lung channel, while the remaining five components act through the stomach channel. Nine of the twenty-one components in QFPD act through the lung channel, while nine belong to the spleen channel. Collectively, LHQW seemingly exerts a heat-clearing effect primarily through the lung channel, whereas QFPD exerts a dampness-eliminating effect primarily through the lung channel (Figure 1).

FIGURE 1
www.frontiersin.org

FIGURE 1. Main channel tropism of Lianhua Qingwen (LHQW) and Qingfei Paidu Decoction (QFPD) based on the TCM theory.

Comparison of LHQW and QFPD active components and their potential targets

LHQW contained 196 active compounds that targeted 249 genes, while QFPD contained 310 active compounds that targeted 261genes (Supplementary Table S2). The constructed Venn diagrams of the LHQW and QFPD active compounds and their targets identified 135 common compounds, while the other compounds were unique to either formula (Supplementary Figure S1A). Moreover, LHQW and QFPD had 226 common targets, whereas LHQW and QFPD targeted 23 and 35 specific targets, respectively (Supplementary Figure S1B). These data suggest that the two formulations may exert similar effects in treating mild-to-moderate COVID-19. The two prescriptions may contain the same components, or different components with the same active compounds, different active compounds with the same targets. However, they also have certain unique compounds and targets, thereby eliciting varying effects in the treatment of mild-to-moderate COVID-19.

PPI analysis of the common active compound targets in LHQW and QFPD

After entering the 226 common targets into the STRING database, we constructed a PPI network, and 49 hub targets were selected to build the PIN (Figure 2A). The results of the PIN topology analysis (Figure 2B) showed that it possessed the robustness and vulnerability topologies required for network research (γ = 0.404, γ < 3), and these targets could be considered hub targets. Following the biological functional clustering of targets in the PIN, three protein modules (score >4) were obtained (Figures 2C–E). These were believed to play an important role in network regulation and may be related to immune function and inflammatory responses. These results indicate that LHQW and QFPD may play comparable roles in treating mild-to-moderate COVID-19 by regulating the protein modules in the PIN.

FIGURE 2
www.frontiersin.org

FIGURE 2. Protein interaction network (PIN) construction and network topology analysis. (A) PIN was constructed using the common host targets of LHQW and QFPD. (B) PIN topology analysis was performed using the Network Analyzer of the Cytoscape software. (C–E) The significant module was identified from the PIN using the molecular complex detection (MCODE) method with a score >4.0; the red and yellow nodes represent hub genes, yellow represents the lower degree and red represents the higher degree.

GO and KEGG pathway enrichment analysis of common targets

The GO enrichment results for the 49 key common genes targeted by LHQW and QFPD (Figure 3A) showed that both formulations might have roles in inflammation, immunity, energy metabolism, and oxidative stress through various biological pathways, such as cytokine activity and responses to xenobiotic stimuli, nutrient levels, and oxidative stress. KEGG enrichment results further showed that the signaling pathways of the two formulations were primarily associated with various viral infections (hepatitis B, Epstein-Barr virus), inflammatory reactions (advanced glycation end-product [AGE]-receptor for AGE [RAGE] signaling pathway involved in diabetic complications, IL-17 signaling pathway), and immunoregulation (Toll-like receptor signaling pathway, PI3K-Akt signaling pathway). These results suggest that LHQW and QFPD may slow the progression of mild and moderate COVID-19 by regulating similar mechanisms, including inflammation, multiple viral infections, immunoregulation, energy metabolism, and oxidative stress (Figure 3B).

FIGURE 3
www.frontiersin.org

FIGURE 3. Gene Ontology (GO) and KEGG enrichment analysis of common genes of LHQW and QFPD targeted during the treatment of mild-to-moderate COVID-19. (A) GO analysis was performed for the top 10 GO terms identified; X-axis represents GO term categories, while the y-axis represents the enrichment score. (B) A KEGG analysis of the genes was performed, and the top 15 pathways are shown. The color of the dot represents the p-value. The dot size represents the number of target genes enriched in the pathway.

Analysis of LHQW and QFPD specificity

Table 1 lists the active compounds and their specific targets obtained from the TCMSP database. After importing the relationship between TCM components as well as active compounds and their specific targets into the Cytoscape software, the “TCM component–active compound–specific target” networks (Figures 4A–D) and three protein modules (score ≥4.44) of specific LHQW and QFPD targets were constructed (Figures 4C,D). These modules suggest that metabotropic glutamate receptors GRM1 and GRM5 play a vital role in conduction and neural network regulation.

TABLE 1
www.frontiersin.org

TABLE 1. Specific targets of the active compounds present in LHQW and QFPD.

FIGURE 4
www.frontiersin.org

FIGURE 4. “TCM component–active compound–specific target” networks and protein modules of LHQW and QFPD. The green diamond represents the TCM components, the orange triangle represents the active compounds, and the cyan fusiform represents the network targets. (A) Interconnection between TCM components, active compounds, and the specific targets of the active compounds in LHQW. (B) Interconnection between TCM components, active compounds, and their specific targets of the active compounds in QFPD. (C) The significant module identified from the specific targets of LHQW. (D) The significant module identified from the specific targets of QFPD.

Validation of specific target-chemical ingredient interactions by molecular docking

mTOR regulates the growth and metabolism of immune cells, as well as macrophage apoptosis (Zou et al., 2020). The therapeutic potential of mTOR inhibitors for COVID-19 has been reported (Fagone et al., 2020). Phospholipase A2 group IVA (PLA2G4A) is closely related to the arachidonic acid metabolic pathway, which can promote the degradation of phospholipids, and produce arachidonic acid, subsequently causing cytokine storm activation (Harizi et al., 2008). The “TCM component–active compound–specific target” networks analysis revealed that GRM1 and GRM5 could core specific proteins of LHQW, whereas mTOR and PLA2G4Awere potential core specific proteins of QFPD.

We next sought to determine which formulation was superior in treating mild-to-moderate COVID-19 with symptoms of depression (Asadi-Pooya & Simani, 2020; Havervall et al., 2021) or inflammation via targeting specific targets (GRM1, GRM5, mTOR and PLA2G4A). We performed molecular docking through GOLD and MOE to analyze the binding capacity between specific targets and the active compounds obtained from the ingredient–target network. ChemScore, ChemPLP, and binding energy values were applied to evaluate the docking results.

Arctiin exhibited better binding activity than the controls ([11C]LY-2428703 and ADX-50938) when docked with GRM1 and GRM5 (Supplementary Tables S3, S4). It primarily formed three hydrogen bonds, each with THR-815, CYS-746, and GLN-660 residues on GRM1 (Figure 5A), and six hydrogen bonds with SER-1809, GLY-624, SER-654, ASN-1747, TYR-1792, and THR-1735 residues on GRM5 (Figure 5G). According to the target-chemical interactions, arctiin interacts with THR-815, LEU-757, and ASN-760 residues on GRM1 (Figure 5D) and TRP-1785, SER-1809, ALA-1810, SER-658, GLY-624, SER-654, ILE-651, and VAL-1740 residues on GRM5 (Figure 5J). In addition, β-carotene exhibited better binding activity than the control ([11C]LY-2428703) when docked with GRM1 (Supplementary Table S4). Other chemical ingredients also formed hydrogen bonds and interacted with different residues on specific targets (Figure 5).

FIGURE 5
www.frontiersin.org

FIGURE 5. Molecular models of active compounds with specific key LHQW target genes. (A) Molecular model of arctiin, (B) β-carotene, (C) and aloe-emodin binding to GRM1. (G) Molecular model of arctiin, (H) corymbosin, (I) and aloe-emodin binding to GRM5. Yellow dashed lines represent H-bonds, thick pink sticks represent the active compounds, and thick green sticks represent residues in the protein-binding site. (D) Interactions between arctiin, (E) β-carotene, (F) aloe-emodin, and GRM1. (J) Interaction between arctiin, (K) corymbosin, (L) aloe-emodin, and GRM5. Blue dashed lines represent H-bonds between the active compounds and residue backbones, and green dashed lines represent H-bonds between the active compounds and residue sidechains. The dashed lines point to the acceptor.

These results suggest that LHQW may be effective in treating psychological disorder symptoms in mild-to-moderate COVID-19 by targeting neurological disease-related targets (GRM1 and GRM5) through arctiin, corymbosin, and aloe-emodin. Meanwhile, QFPD is associated with inflammation and functions by targeting immunity- and inflammation-related factors (mTOR and PLA2G4A) through isofucosterol, baicalein, nobiletin, oroxylinA, epiberberine, and piperlonguminine (Figure 6, Supplementary Tables S5, S6).

FIGURE 6
www.frontiersin.org

FIGURE 6. Molecular models of active compounds with specific QFPD target genes. (A) Molecular model of isofucosterol binding to mTOR. (B) Molecular model of baicalein, (C) nobiletin, (D) oroxylin A, (I) epiberberine, (G) piperlonguminine, (K) baicalein, (L) and isofucosterol binding to PLA2G4A. Yellow dashed lines represent H-bonds, thick pink sticks represent active compounds, and thick green sticks represent residues at the protein-binding site. (E) Interaction of isofucosterol, (F) baicalein, (G) nobiletin, and mTOR. (H) Interaction of oroxylin A, (M) epiberberine, (N) piperlonguminine, (O) baicalein, (P) isofucosterol, and PLA2G4A. Blue dashed lines represent H-bonds between the active compounds and residue backbones, and the green dashed lines indicate H-bonds between the active compounds and residue sidechains. The dashed lines point to the acceptor.

Validation of specific target-chemical ingredient interactions by molecular dynamic simulations

Amber18 was employed to analyze the binding capacity between specific targets (GRM1 and GRM5) and active compounds (arctiin and β-carotene) that demonstrated effects superior to the positive controls. The results of the target–ingredient binding (RMSD, RMSF, and SASA values) are shown in Figure 7 and Supplementary Figure S2, while those of the Gibbs free energy analyses are summarized in Table 2. RMSD can be used to characterize the stability of the active compounds that bind targets. The RMSD of arctiin binding to GRM1 is stable after 20 ns, with an RMSD value of approximately 0.9 Å at 20–35 ns (average RMSD = 0.69 Å; Figure 7A). β-carotene binding to GRM1 was stable at approximately 0.3 Å at 0–20 ns and approximately 0.5 Å at 25–50 ns, with an average RMSD value of 0.42 Å (Figure 7B). The RMSD of arctiin binding to GRM5 was stable after 5 ns, and the RMSD value was stable at approximately 0.7 Å at 5–25 ns and approximately 0.9 Å at 25–50 ns (average RMSD = 0.86 Å; Figure 7C). These data suggest that arctiin and β-carotene can bind stably with GRM1, and arctiin binds stably with GRM5. Hence, arctiin and β-carotene can stably bind to their targets, and may represent potential targeted drugs.

FIGURE 7
www.frontiersin.org

FIGURE 7. Root mean square deviation (RMSD) of molecular dynamic simulations. (A) RMSD of arctiin binding to GRM1. (B) RMSD of β-carotene binding to GRM1 (C) RMSD of arctiin binding to GRM5.

TABLE 2
www.frontiersin.org

TABLE 2. Gibbs free energy of arctiin and β-carotene with GRM1 and GRM5.

Discussion

The outbreak and spread of the novel coronavirus disease, COVID-19, turned into a global public health emergency with serious socio-economic repercussions (Shen et al., 2021). Effective treatment of COVID-19 with allopathic medicine has faced several challenges, including its applicability to children and other specific subgroups of patients, as well as its relatively high cost (Ansems et al., 2021; Agarwal et al., 2022). Meanwhile, treatment with TCM can alleviate symptoms and prevent fatal deterioration of COVID-19 patients. As such, TCM has been incorporated into various COVID-19 treatment guidelines in China, thereby playing an indispensable role in its management (Li Q. et al., 2020; Huang et al., 2021).

LHQW and QFPD have shown efficacy in the treatment of mild-to-moderate COVID-19 cases. Indeed, LHQW has been widely used in China (Li L. C. et al., 2020). Clinically, LHQW has exhibited numerous advantages in treating COVID-19 (Runfeng et al., 2020; Hu et al., 2021). LHQW inhibits virus transmission and regulates immune function in influenza virus infections. In addition, QFPD exerts therapeutic effects via immune regulation, anti-infection and anti-inflammatory properties, and multi-organ protection in patients with mild-to-moderate COVID-19. In fact, the overall treatment efficacy rate for COVID-19 reportedly exceeds 90% (Luo et al., 2020; Wang et al., 2021).

According to clinical studies, an abnormal concentration of various cytokines (such as TNF-α, IL-6, IL-8) and an excessive release of reactive oxygen species in COVID-19 could cause an overactive immune response in COVID-19 (Behrens & Koretzky, 2017; Huang et al., 2020; Pedersen & Ho, 2020). More specifically, the continuous production of inflammatory factors can cause cytokine storm, which aggravates lung tissue damage and results in oxidative stress damage (Bettelli et al., 2006; Dai et al., 2018; Erlich et al., 2020). Quercetin and kaempferol present in LHQW can inhibit the overexpression of cytokines and alleviate multiple organ injuries in COVID-19 (Zhang et al., 2021). Similarly, QFPD contains active compounds capable of alleviating inflammation and preventing the initiation of the cytokine storm in COVID-19 by regulating multiple pathways (Zhao et al., 2021). Our research supports these findings as we show that LHQW and QFPD can regulate the immune system, reduce inflammation, and inhibit oxidative stress damage via related pathways, including those involved in immunity and inflammation, apoptosis, and viral infection, ultimately alleviating COVID-19 symptoms.

Although there are similarities between the two formulations, LHQW and QFPD have specific targets and contain different active compounds. In the present study, the active compounds in LHQW exhibited an excellent ability to bind to specific targets, namely GRM1 and GRM5. GRM1 and GRM5 are type I metabotropic glutamate receptors involved in various mental illnesses, including depression, anxiety, and cognitive impairment by regulating the function of neurotransmitters such as glutamate, dopamine, and norepinephrine (Pittaluga, 2016; Borbely et al., 2022). The long-term administration of two antidepressants (imipramine and citalopram) in mice inhibited the transcription of ionotropic glutamate receptors in multiple brain regions, suggesting that glutamate receptors are associated with depression (Boyer et al., 1998). Meanwhile, QFPD may play a role in COVID-19 treatment by targeting mTOR and PLA2G4A. As a serine-threonine kinase (Huang et al., 2020), decreased expression of mTOR can interfere with viral replication and protein synthesis.

Drugs that target the mTOR of Middle East respiratory syndrome coronaviruses are currently available (Kindrachuk et al., 2015). The “SARS-COV-2–human PPI network” constructed in this study, indicates that mTOR may be a potential key target of COVID-19, and its inhibitors may elicit curative anti-SARS-COV-2 effects (Gordon et al., 2020). Meanwhile, PLA2G4A can participate in arachidonic acid synthesis and produce large amounts of inflammatory factors, leading to cytokine storm initiation (Harizi et al., 2008). Current clinical studies have shown that the expression levels of inflammatory factors (IL-2, IL-7, IFN-γ, and TNF-α) positively correlate with the severity of COVID-19 (Huang et al., 2020), suggesting that PLA2G4A may be a key target in the development of COVID-19.

Our findings indicate that LHQW and QFPD may treat COVID-19 via pathways associated with inflammation, immunity, apoptosis, oxidative stress, etc. Meanwhile, the unique active compounds (Arctiin, Corymbosin, and Aloe-emodin) in LHQW target neurological disease-related genes (GRM1 and GRM5), and the special active compounds (isofucosterol, baicalein, nobiletin, oroxylin A, epiberberine, and piperlonguminine) in QFPD target immunity- and inflammation-related genes (mTOR and PLA2G4A). These indicated that LHQW might be suitable for mild-to-moderate COVID-19 with nervous system symptoms (Asadi-Pooya & Simani, 2020; Havervall et al., 2021) and QFPD may have a potential effect on regulating oxidative stress damage and inflammatory symptoms. In addition, the active compounds identified in this study can be further modified and optimized as candidate drugs for mild-to-moderate COVID-19.

However, this study has the following limitations, the active compounds are from TCMSP, and the assessment of relationships between the different compounds was insufficient. Hence, further studies should be performed to verify these results in vitro and in vivo employing, genomics, metabolomics, and proteomics approaches. Our study may be the basis for LHQW and QFPD in treating mild-to-moderate COVID-19, providing candidate active compounds for the clinical application of LHQW and QFPD.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

QS and AL designed the experiments and wrote the manuscript; LL, XW, XG, and YL analyzed the data, and all of the authors reviewed the manuscript.

Funding

This project was supported by the Projects of Hebei Provincial Administration of Traditional Chinese Medicine (NO. 2022087); Doctoral Research Funding Project of Hebei University of Chinese Medicine (NO. BSZ2020011); Science and Technology Project of Hebei Education Department (NO. QN2021106).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.1013428/full#supplementary-material

References

Agarwal, A., Rochwerg, B., Lamontagne, F., Siemieniuk, R. A., Agoritsas, T., and Askie, L. A living WHO guideline on drugs for Covid-19. Pract. Guidel. BMJ (Clinical research ed.). 377 (2022) m3379. doi:10.1136/bmj.m3379

PubMed Abstract | CrossRef Full Text | Google Scholar

Ansems, K., Grundeis, F., Dahms, K., Mikolajewska, A., Thieme, V., Piechotta, V., et al. Remdesivir for the treatment of COVID-19. Cochrane Database Syst. Rev. 8 (2021) CD014962. doi:10.1056/NEJMoa2007764

PubMed Abstract | CrossRef Full Text | Google Scholar

Asadi-Pooya, A. A., and Simani, L., Central nervous system manifestations of COVID-19: A systematic review. J. Neurol. Sci. 413 (2020). 116832. doi:10.1016/j.jns.2020.116832

PubMed Abstract | CrossRef Full Text | Google Scholar

Behrens, E. M., and Koretzky, G. A., Review: Cytokine storm syndrome: Looking toward the precision medicine era. Arthritis Rheumatol. 69 (2017) 1135–1143. doi:10.1002/art.40071

PubMed Abstract | CrossRef Full Text | Google Scholar

Bettelli, E., Carrier, Y., Gao, W., Korn, T., Strom, T. B., Oukka, M., et al. Reciprocal developmental pathways for the generation of pathogenic effector TH17 and regulatory T cells. Nature 441 (2006) 235–238. doi:10.1038/nature04753

PubMed Abstract | CrossRef Full Text | Google Scholar

Borbely, E., Simon, M., Fuchs, E., Wiborg, O., Czeh, B., and Helyes, Z., Novel drug developmental strategies for treatment-resistant depression. Br. J. Pharmacol. 179 (2022) 1146–1186. doi:10.1111/bph.15753

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyer, P. A., Skolnick, P., and Fossom, L. H., Chronic administration of imipramine and citalopram alters the expression of NMDA receptor subunit mRNAs in mouse brain. A quantitative in situ hybridization study. J. Mol. Neurosci. 10 (1998) 219–233. doi:10.1007/BF02761776

PubMed Abstract | CrossRef Full Text | Google Scholar

Carfì, A., Bernabei, R., and Landi, F. Persistent symptoms in patients after acute COVID-19. Jama 324 (2020) 603–605. doi:10.1001/jama.2020.12603

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Wang, Y. K., Gao, Y., Hu, L. S., Yang, J. W., Wang, J. R., et al. Protection against COVID-19 injury by qingfei paidu decoction via anti-viral, anti-inflammatory activity and metabolic programming. Biomed. Pharmacother. 129 (2020) 110281. doi:10.1016/j.biopha.2020.110281

PubMed Abstract | CrossRef Full Text | Google Scholar

Chu, L., Huang, F., Zhang, M., Huang, B., and Wang, Y., Current status of traditional Chinese medicine for the treatment of COVID-19 in China. Chin. Med. 16 (2021) 63. doi:10.1186/s13020-021-00461-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Commission, C. P., Pharmacopoeia of the people’s Republic of China., China Medical Science Press, Beijing, China, 2015.

Google Scholar

Dai, J., Gu, L., Su, Y., Wang, Q., Zhao, Y., Chen, X., et al. Inhibition of curcumin on influenza A virus infection and influenzal pneumonia via oxidative stress, TLR2/4, p38/JNK MAPK and NF-κB pathways. Int. Immunopharmacol. 54 (2018) 177–187. doi:10.1016/j.intimp.2017.11.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Erlich, J. R., To, E. E., Liong, S., Brooks, R., Vlahos, R., O'Leary, J. J., et al. Targeting evolutionary conserved oxidative stress and immunometabolic pathways for the treatment of respiratory infectious diseases. Antioxid. Redox Signal. 32 (2020) 993–1013. doi:10.1089/ars.2020.8028

PubMed Abstract | CrossRef Full Text | Google Scholar

Fagone, P., Ciurleo, R., Lombardo, S. D., Iacobello, C., Palermo, C. I., Shoenfeld, Y., et al. Transcriptional landscape of SARS-CoV-2 infection dismantles pathogenic pathways activated by the virus, proposes unique sex-specific differences and predicts tailored therapeutic strategies. Autoimmun. Rev. 19 (2020) 102571. doi:10.1016/j.autrev.2020.102571

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, D., Niu, M., Wei, S. Z., Zhang, C. E., Zhou, Y. F., Yang, Z. W., et al. Identification of a pharmacological biomarker for the bioassay-based quality control of a thirteen-component TCM formula (Lianhua qingwen) used in treating influenza A virus (H1N1) infection. Front. Pharmacol. 11 (2020) 746. doi:10.3389/fphar.2020.00746

PubMed Abstract | CrossRef Full Text | Google Scholar

Gordon, D. E., Jang, G. M., Bouhaddou, M., Xu, J., Obernier, K., O'Meara, M. J., et al. A SARS-CoV-2-human protein-protein interaction map reveals drug targets and potential drug-repurposing. Nature. (2020).583. 459–468. doi:10.1038/s41586-020-2286-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Harizi, H., Corcuff, J. B., and Gualde, N., Arachidonic-acid-derived eicosanoids: Roles in biology and immunopathology. Trends Mol. Med. 14 (2008) 461–469. doi:10.1016/j.molmed.2008.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Havervall, S., Rosell, A., Phillipson, M., Mangsbo, S. M., Nilsson, P., Hober, S., et al. Symptoms and functional impairment assessed 8 Months after mild COVID-19 among health cae workers. Jama 325 (2021) 2015–2016. doi:10.1001/jama.2021.5612

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, K., Guan, W. J., Bi, Y., Zhang, W., Li, L., Zhang, B., et al. Efficacy and safety of lianhuaqingwen capsules, a repurposed Chinese herb, in patients with coronavirus disease 2019: A multicenter, prospective, randomized controlled trial [phytomedicine 85 (2021) 153242]. Phytomedicine 93 (2021) 153775. doi:10.1016/j.phymed.2021.153775

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Wang, Y., Li, X., Ren, L., Zhao, J., Hu, Y., et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 395 (2020) 497–506. doi:10.1016/S0140-6736(20)30183-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, K., Zhang, P., Zhang, Z., Youn, J. Y., Wang, C., Zhang, H., et al. Traditional Chinese Medicine (TCM) in the treatment of COVID-19 and other viral infections: Efficacies and mechanisms. Pharmacol. Ther. 225 (2021) 107843. doi:10.1016/j.pharmthera.2021.107843

PubMed Abstract | CrossRef Full Text | Google Scholar

Kindrachuk, J., Ork, B., Hart, B. J., Mazur, S., Holbrook, M. R., Frieman, M. B., et al. Antiviral potential of ERK/MAPK and PI3K/AKT/mTOR signaling modulation for Middle East respiratory syndrome coronavirus infection as identified by temporal kinome analysis. Antimicrob. Agents Chemother. 59 (2015) 1088–1099. doi:10.1128/AAC.03659-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L. C., Zhang, Z. H., Zhou, W. C., Chen, J., Jin, H. Q., Fang, H. M., et al. Lianhua Qingwen prescription for Coronavirus disease 2019 (COVID-19) treatment: Advances and prospects. Biomed. Pharmacother. 130 (2020a) 110641. doi:10.1016/j.biopha.2020.110641

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Wang, H., Li, X., Zheng, Y., Wei, Y., Zhang, P., et al. The role played by traditional Chinese medicine in preventing and treating COVID-19 in China. Front. Med. 14 (2020b) 681–688. doi:10.1007/s11684-020-0801-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, M., Gao, Y., Yuan, Y., Yang, K., Shi, S., Zhang, J., et al. Efficacy and safety of integrated traditional Chinese and western medicine for corona virus disease 2019 (COVID-19): A systematic review and meta-analysis. Pharmacol. Res. 158 (2020) 104896. doi:10.1016/j.phrs.2020.104896

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, E., Zhang, D., Luo, H., Liu, B., Zhao, K., Zhao, Y., et al. Treatment efficacy analysis of traditional Chinese medicine for novel coronavirus pneumonia (COVID-19): An empirical study from wuhan, hubei province, China. Chin. Med. 15 (2020) 34. doi:10.1186/s13020-020-00317-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, N. N., Hoang, V. T., Dao, T. L., Meddeb, L., Cortaredona, S., Lagier, J. C., et al. Long-term persistence of olfactory and gustatory disorders in COVID-19 patients. Front. Med. 9 (2022) 794550. doi:10.3389/fmed.2022.794550

CrossRef Full Text | Google Scholar

Pedersen, S. F., and Ho, Y. C., SARS-CoV-2: A storm is raging. J. Clin. Invest. 130 (2020) 2202–2205. doi:10.1172/JCI137647

PubMed Abstract | CrossRef Full Text | Google Scholar

Pittaluga, A., Presynaptic release-regulating mGlu1 receptors in central nervous system. Front. Pharmacol. 7 (2016) 295. doi:10.3389/fphar.2016.00295

PubMed Abstract | CrossRef Full Text | Google Scholar

Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. Tcmsp: A database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform. 6 (2014) 13. doi:10.1186/1758-2946-6-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Runfeng, L., Yunlong, H., Jicheng, H., Weiqi, P., Qinhai, M., Yongxia, S., et al. Lianhuaqingwen exerts anti-viral and anti-inflammatory activity against novel coronavirus (SARS-CoV-2). Pharmacol. Res. 156 (2020) 104761. doi:10.1016/j.phrs.2020.104761

PubMed Abstract | CrossRef Full Text | Google Scholar

Salian, V. S., Wright, J. A., Vedell, P. T., Nair, S., Li, C., Kandimalla, M., et al. COVID-19 transmission, current treatment, and future therapeutic strategies. Mol. Pharm. 18 (2021) 754–771. doi:10.1021/acs.molpharmaceut.0c00608

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, P., Li, J., Tu, S., Wu, Y., Peng, Y., Chen, G., et al. Positive effects of lianhuaqingwen granules in COVID-19 patients: A retrospective study of 248 cases. J. Ethnopharmacol. 278 (2021) 114220. doi:10.1016/j.jep.2021.114220

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 49 (2021) D605-D612. doi:10.1093/nar/gkaa1074

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Zhu, H., Li, M., Liu, Y., Lai, H., Yang, Q., et al. Efficacy and safety of qingfei paidu decoction for treating COVID-19: A systematic review and meta-analysis. Front. Pharmacol. 12 (2021) 688857. doi:10.3389/fphar.2021.688857

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, C., Mao, X., Huang, J., Ding, Y., Wu, J., Dong, S., et al. Kobas 2.0: A web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 39 (2011) W316–W322. doi:10.1093/nar/gkr483

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Gao, R., Zhou, Z., Tang, X., Lin, J., Wang, L., et al. A network pharmacology based approach for predicting active ingredients and potential mechanism of Lianhuaqingwen capsule in treating COVID-19. Int. J. Med. Sci. 18 (2021) 1866–1876. doi:10.7150/ijms.53685

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Yang, Z., Chen, Z. L., Li, Z. N., Yue, S. J., Li, J. J., et al. Efficacy and safety of "three Chinese patent medicines and three TCM prescriptions" for COVID-19: A systematic review and network meta-analysis. Evid. Based. Complement. Altern. Med. 2022 (2022a) 4654793. doi:10.1155/2022/4654793

CrossRef Full Text | Google Scholar

Zhang, L., Ma, Y., Shi, N., Tong, L., Liu, S., Ji, X., et al. Effect of qingfei paidu decoction combined with western medicine treatments for COVID-19: A systematic review and meta-analysis. Phytomedicine. 102 (2022b) 154166. doi:10.1016/j.phymed.2022.154166

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Tian, S., Lu, D., Yang, J., Zeng, H., Zhang, F., et al. Systems pharmacological study illustrates the immune regulation, anti-infection, anti-inflammation, and multi-organ protection mechanism of Qing-Fei-Pai-Du decoction in the treatment of COVID-19. Phytomedicine. 85 (2021) 153315. doi:10.1016/j.phymed.2020.153315

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhuang, W., Fan, Z., Chu, Y., Wang, H., Yang, Y., Wu, L., et al. Chinese patent medicines in the treatment of coronavirus disease 2019 (COVID-19) in China. Front. Pharmacol. 11 (2020) 1066. doi:10.3389/fphar.2020.01066

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, Z., Tao, T., Li, H., and Zhu, X., mTOR signaling pathway and mTOR inhibitors in cancer: progress and challenges. Cell Biosci. 10 (2020) 31. doi:10.1186/s13578-020-00396-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: COVID-19, traditional Chinese medicine, systems biology, GRM1/5, arctiin, β-carotene

Citation: Li L, Wang X, Guo X, Li Y, Song Q and Li A (2022) Network pharmacology and computer-aided drug design to explored potential targets of Lianhua Qingwen and Qingfei Paidu decoction for COVID-19. Front. Pharmacol. 13:1013428. doi: 10.3389/fphar.2022.1013428

Received: 07 August 2022; Accepted: 30 August 2022;
Published: 23 September 2022.

Edited by:

Pier Maria Fornasari, REGENHEALTHSOLUTIONS, Italy

Reviewed by:

Taif Shah, Kunming University of Science and Technology, China
Mithun Rudrapal, Rasiklal M. Dhariwal Institute of Pharmaceutical Education and Research, India

Copyright © 2022 Li, Wang, Guo, Li, Song and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Qiuhang Song, songqiuhang@163.com; Aiying Li, lay1963@126.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.