Investigation of the Potential Mechanism Governing the Effect of the <i>Shen Zhu San</i> on COVID-19 by Network Pharmacology

Wang, Yuxuan; Ru, Yuhua; Zhuo, Guowei; Sheng, Maozheng; Wang, Shuangqiu; Ma, Jiarui; Zhou, Chongyi; Sun, Xiaohe; Zeng, Yanqi; Zhang, Ya; Li, Hui; Lu, Zhigang; Wu, Depei; Wu, Minghua

doi:https://doi.org/10.1155/2020/8468303

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Results Discussion Conclusion Methods Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 8468303 | https://doi.org/10.1155/2020/8468303

Investigation of the Potential Mechanism Governing the Effect of the Shen Zhu San on COVID-19 by Network Pharmacology

Yuxuan Wang,^1,2,3Yuhua Ru,⁴Guowei Zhuo,⁵Maozheng Sheng,⁶Shuangqiu Wang,⁷Jiarui Ma,⁶Chongyi Zhou,⁷Xiaohe Sun,⁷Yanqi Zeng,⁵Ya Zhang,⁵Hui Li,¹and Zhigang Lu² et al.

Academic Editor: Mario Ledda

Received01 Mar 2020

Revised30 Jun 2020

Accepted22 Sept 2020

Published09 Nov 2020

Abstract

Background. Since December 2019, coronavirus disease 2019 (COVID-19) due to SARS-CoV-2 infection has emerged in Wuhan and rapidly spread throughout China and even to other countries. Combined therapy with modern medicine and traditional Chinese medicine has been proposed, in which Shen Zhu San (SZS) was regarded as one of the basic prescriptions. Methods. Network pharmacological approaches along with candidate compound screening, target prediction, target tissue location, protein-protein interaction network, gene ontology (GO), KEGG enrichment analyses, and gene microarray analyses were applied. Results. A total of 627 targets of the 116 active ingredients of SZS were identified. Targets in immune cells and tissues were much more abundant than those in other tissues. A total of 597 targets were enriched in the GO biological cellular process, while 153 signaling pathways were enriched according to the KEGG analysis. A total of 450 SARS-related targets were integrated and intersected with the targets of SZS to identify 40 common targets that were significantly enriched in five immune function aspects of the immune system process during GO analysis. Several inflammation-related pathways were found to be significantly enriched throughout the study. Conclusions. The therapeutic mechanisms of the effects of SZS on COVID-19 potentially involve four effects: suppressing cytokine storms, protecting the pulmonary alveolar-capillary barrier, regulating the immune response, and mediating cell death and survival.

1. Introduction

In early December 2019, a series of pneumonia cases of unknown cause were first confirmed in Wuhan, Hubei Province, China [1]. According to the World Health Organization (WHO), by June 30, 2020, there have been over ten million cases and five hundred thousand deaths in 215 countries and regions worldwide. Unfortunately, the number is still increasing. Efforts made by the Chinese Government to curb outbreaks included initiating first measures to isolate Wuhan, which were extended to the whole Hubei province, standing 35 million residents during the Chinese Spring Festival [2]. Meanwhile, 30 provinces of China had launched a first-level response to the major public health emergencies [3]. Because of the severity of this outbreak, the WHO declared it a Public Health Emergency of International Concern on January 30, 2020 [4], and then increased assessment of the risk of spread and impact of COVID-19 to very high at a global level on February 28, 2020 [5].

Chinese scientists rapidly isolated and confirmed a novel type of coronavirus named SARS-CoV-2 by the International Committee on the Taxonomy of Viruses on February 11, 2020 [6, 7]. Coronaviruses are positive-sense single-stranded RNA viruses, which belong to the family Coronaviridae. SARS-CoV-2 is the seventh member of the coronavirus family and the third coronavirus that causes human fatal illness after SARS-CoV and MERS-CoV [8]. SARS-CoV, which caused outbreaks of severe acute respiratory syndrome (SARS) in Guangdong Province in China in 2002, shares 79.5% genetic sequence similarity with SARS-CoV-2 as well as the same cell entry receptor [8], and the amino acid sequence identity between the SARS-CoV-2 and the SARS-CoV S-proteins is 76.47% [9]. Some preliminary studies suggested that SARS-CoV-2 was similar to SARS-CoV to some extent, based on full-length genome phylogenetic analyses [8, 10] and the putative similarity of cell entry mechanism and human cell receptor utilization [8, 11, 12]. The disease caused by the novel coronavirus was named coronavirus disease 2019 (COVID-19) by the WHO [7]. The clinical signs and symptoms have been reported mainly as viral pneumonia with fever and dyspnea, a few even developed acute respiratory distress syndrome (ARDS), with chest radiographs showing bilateral pulmonary infiltration [1, 13–15].

As there are no specific antiviral therapies for COVID-19 and the main treatments are supportive at present, the Proposed Diagnosis and Treatment (5^thedition) issued by China’s National Health Commission integrated therapy with traditional Chinese medicine (TCM) and modern medicine [16]. To date, TCM intervention has been officially specified as an integrant therapeutic strategy in several provinces in China for patients who develop COVID-19 symptoms. The experts in TCM who were assigned by China’s National Health Commission supplied a series of therapies to treat the COVID-19. Among these therapies, a traditional formula called Shen Zhu San (SZS, in Chinese: 神术散) was regarded as one of the basic prescriptions for the treatment of COVID-19 [17–19].

SZS, which is a thousand-year-old TCM formula, was first recorded in an ancient Chinese medical book called “Yang’s Hereditary Medical Formulary” and was later included in the “Prescriptions of Peaceful Benevolent Dispensary”, the official medical book of the Song Dynasty. SZS is mainly composed of seven medicinal herbs: Atractylodis Rhizoma, Angelicae Dahuricae Radix, Asari Radix Et Rhizoma, Glycyrrhizae Radix Et Rhizoma, Notopterygii Rhizoma Et Radix, Ligustici Rhizoma Et Radix, and Chuanxiong Rhizoma [20, 21]. In the long course of history, SZS has been widely used throughout China and has been proven to have therapeutic effects on epidemic infectious diseases. Modern toxicological studies demonstrated that the herbs in SZS have a high margin of drug safety [22–29]. According to ancient records and modern application, SZS is good for treating respiratory tract diseases with fever, shortness of breath, cough, nasal obstruction, headache, dizziness, body aches, and hematochezia caused by the wind-cold damp pathogen, for which the symptoms and TCM pathogenesis are highly consistent with those COVID-19. As a new pharmacological discipline, integrating network pharmacology emphasizes the concept of a “multicomponent, multitarget therapeutic network” and highlights the overall concept of TCM [30]. Network pharmacology offers a novel concept for understanding the multitargeted mechanism of the treatment of complex diseases with TCM [31]. In this study, we explored the potential mechanism of SZS in COVID-19 treatment by utilizing integrating network pharmacological approaches.

2. Results

2.1. Potential Mechanism of SZS in COVID-19 Treatment

2.1.1. Identification of Active Compounds and Targets of SZS

Using the TCMSP and SymMap databases, the chemical ingredients were identified based on the criteria of OB ≥ 30% and DL ≥ 0.18: 9 in Atractylodis Rhizoma, 18 in Angelicae Dahuricae Radix, 6 in Asari Radix Et Rhizoma, 71 in Glycyrrhizae Radix Et Rhizoma, 12 in Notopterygii Rhizoma Et Radix, 1 in Ligustici Rhizoma Et Radix, and 7 in Chuanxiong Rhizoma. There were 116 chemical ingredients chosen for further investigation after the removal of duplicates. A total of 2219 human component-target interactions were generated from the 116 active ingredients of SZS. Duplicates of the validated and predicted component targets were eliminated, and then the 627 targets of the active ingredients of SZS were screened (Table S1), among which the main compounds and targets are shown in Table 1.

2.1.2. Compound-Target Network and Target Tissue Location

The 627 targets of the active ingredients of SZS mentioned above were used to construct a component-target network classified by tissue. Genes are commonly expressed in multiple tissues and cell types. In this study, as shown in Figure 1, the targets were located in certain tissues where they had the highest mRNA expression level according to the BioGPS database [32, 33], because we were interested in the expression patterns of the targets in specialized cells and tissues, including immune cells and tissues, lung, brain, heart, kidney, and liver. The mRNA-specific expression patterns of targets offer significant clues for understanding the underlying therapeutic effects of SZS. Specifically, there were 207 targets located in immune cells and tissues: 28 in the lung, 126 in the brain, 34 in the heart, 20 in the kidney, and 42 in the liver. The nodes with different colors represent the corresponding relationships of targets and organizations, as illustrated in the lower-left corner of Figure 1. The number and the total degree of all the nodes of component targets in immune cells and tissues (marked in yellow) were much greater than those in other tissues, which indicated that SZS might exert therapeutic effect by regulating immune function.

2.1.3. Analysis of Targets in PPI Network

A PPI network was constructed to analyze and understand the mechanisms of the effects of SZS according to the protein-protein interactions. A total of 627 targets were inputted into the STRING database and the target proteins were rejected independently of the PPI network (Figure 2(a)) [34]. There were 201 nodes (representing active proteins) and 292 edges (representing the interaction between the active proteins and other proteins) in the network. All the nodes represented queried proteins and the first layer of interactors, and the protein-protein associations represented as edges were implied to be specific and meaningful. The results of the network topology analyses were as follows: the average degree of the nodes was 2.9054, and there were 81 nodes with a degree higher than the average degree. The average betweenness centrality of nodes was 0.0562, and there were 34 nodes with a higher than average betweenness centrality. The significance of key proteins was analyzed based on the comprehensive assessment of the degree, betweenness centrality, and closeness centrality of the nodes exported from the STRING database. According to a threshold degree value of ≥6, 31 critical nodes were further identified and are shown in bubble charts (Figure 2(b)).

(a)

(b)

Figure 2

The PPI network for SZS and hub target analysis. (a) The nodes indicate proteins, and the edges represent protein-protein associations. The cyan edges represent interactions from the curated database, and the purple edges were experimentally determined. The green, red, and blue edges represent the interactions that were determined according to the gene neighborhood, fusions, and co-occurrence, respectively. (b) The centrality of targets was evaluated according to the degree centrality, betweenness centrality, and closeness centrality, which exhibited variation in terms of the y-axis values, colors of nodes, and sizes of nodes.

2.1.4. GO Functional and KEGG Pathway Enrichment Analysis

As a bioinformatics analysis tool, GO defines the input genes by describing their function and the relationships between concepts [35]. After collection of data by ClueGO and CluePedia, GO annotation and enrichment of the SZS target protein genes were performed for three concepts of cell composition, molecular function, and biological process. The GO enrichment analysis of 627 potential targets (Figure 3(a)) showed that the main BP terms enriched in the second GO class were cellular process, biological regulation, and response to the stimulus. For further exploration of the changes in these biological functions, we selected the top 20 BP in the third GO class to generate bubble charts according to the value (Figure 3(b)). The results showed that the mechanism of SZS was mainly related to the regulation of biological quality, oxygen-containing compound response, drug response, chemical response, organic substance response, toxic substance response, chemical stimulus cellular response, ion transport, organic cyclic compound response, and nitrogen compound response. Among these processes, 597 genes were enriched in the BP involved in the cellular process, 531 genes were enriched in processes of biological regulation, and 515 genes were enriched in stimulus response.

(a)

(b)

(c)

Figure 3

GO functional and KEGG enrichment analysis. (a) The second GO class enrichment statistics for 627 targets. The x-axis represents the GO terms: red for BP terms, gray for MF terms, and blue for CC terms. The y-axis shows the number of enriched genes. (b) The bubble diagram shows the top 20 BP terms on the y-axis, while the annotated gene counts for the BP terms are presented on the x-axis. (c) KEGG pathway annotation. The x-axis represents the number of genes in the given classification. The details of the GO and KEGG analyses are presented in Tables S2 and S3, respectively.

There were a total of 153 significant signaling pathways identified by the KEGG enrichment analysis. The bar plot of the KEGG pathway annotation is presented in Figure 3(c). Environmental information processing was mainly enriched in terms of signal transduction and signaling molecules and interactions. Human diseases mainly involved cancers, infectious diseases, and endocrine and metabolic diseases. The data indicated that the targets of SZS were widely involved in regulating metabolism and cell growth and death. Furthermore, there were a few targets involved in genetic information processing. The specific pathways are illustrated in the Discussion section in detail.

2.2. Understanding SZS by Analogy: Investigation of the Potential Mechanism of the Effects of SZS in SARS-CoV Pneumonia Treatment

2.2.1. Screening of SARS-Related Targets and Determination of Common Targets between SZS and SARS

To date, there have been few data available describing the molecular mechanisms underlying COVID-19 so far. To better understand the potential therapeutic effect of SZS on COVID-19, we investigated the potential mechanism of the effects of SZS in SARS-CoV pneumonia treatment. Eighteen years ago, there was an outbreak of SARS-CoV pneumonia in Guangdong, China. As one of the few lethal coronaviruses, SARS-CoV-2 shares 79.5% of its sequence with SARS-CoV and has a similar receptor-binding domain structure to that of SARS-CoV [8, 10], which is the reason that the International Committee on the Taxonomy of Viruses recognized the novel coronavirus as closely related to SARS-CoV. The International Committee on the Taxonomy of Viruses decided to formally name the novel coronavirus SARS-CoV-2 [6]. Regardless of the unclear pathogenesis, the novel coronavirus has caused clusters of fatal pneumonia with a clinical presentation greatly resembling SARS, including the development of ARDS and cytokine storm-caused lethality [1, 13, 14, 36, 37]. Moreover, front-line medical workers announced that the abnormalities in laboratory results and possible pathophysiology of COVID-19, such as cellular immune deficiency, coagulation activation, and multiple organ injury, were similar to those of SARS-CoV pneumonia [14, 15, 36]. Thus, investigating the underlying mechanisms of the effects of SZS in relieving SARS-CoV pneumonia may partly benefit the current fight against COVID-19.

Hence, we integrated totally 450 disease targets related to SARS from the MalaCards database and GEO gene microarray analysis as described in Methods 5.7. The differentially expressed genes are presented in a heatmap (Figure 4(a)) and a volcano plot (Figure 4(b)). As illustrated by the Venn diagram (Figure 4(c)), the disease targets and the component targets intersected to identify 40 common targets, which implies the importance of understanding the role of SZS in COVID-19 treatment by analogy to that in SARS-CoV pneumonia treatment.

(a)

(b)

(c)

In addition, the compound-target networks as well as the relationships among the herb-compound-target-pathway interactions were constructed to better comprehend the complicated molecular mechanisms (Figure 5).

2.2.2. Analysis of Cotargets in PPI Network

Forty cotargets were inputted into the STRING database according to the criteria of a confidence level greater than 0.9 and the rejection of target proteins excluded by the PPI network (Figure 6(a)). There were 33 nodes and 133 edges in the network. The results of the network topology analyses were as follows: the average degree of the nodes was 8.0606, and there were 15 nodes with an average degree higher than average. The average betweenness centrality of the nodes was 0.0324, and there were 9 nodes with higher than average betweenness centrality. The significance of the critical proteins was analyzed based on the degree, betweenness centrality, and closeness centrality of the node exported from STRING database. According to the criterion of a degree ≥6, 21 critical nodes were identified and are presented in bubble charts (Figure 6(b)).

(a)

(b)

2.2.3. GO Functional Enrichment Analysis

The most significantly enriched GO terms of the 40 cotargets are shown according to BP (red bar), MF (gray bar), and CC (blue bar) (Figure 7(a)). The main BP terms were cellular process, regulation of biological process, and biological regulation. Furthermore, the 40 cotargets were found to be enriched in the immune system process by utilizing the ClueGO and CluePedia plugins in Cytoscape (Figure 7(b)) [38–40]. The enrichment of the immune system process was determined to a kappa score of ≥0.4 and a value set ≤0.01. The enriched GO terms were represented as nodes, whose sizes represented the significance of the enrichment of the term. The pie chart shows the proportion of each group associated with the 40 target genes (Figure 7(b)). There were 57 terms and 322 edges contained in the GO enrichment map. The terms were significantly enriched in terms of five immune functions: the regulation of the adaptive immune response, lymphocyte activation involved in the immune response, alpha-beta T-cell differentiation, positive regulation of myeloid cell differentiation, and regulation of the interferon-gamma-mediated signaling pathway.

(a)

(b)

Figure 7

GO functional enrichment analysis and pathway mapping. (a) The second level GO enrichment statistics for 40 cotargets. (b) The cotarget genes were mapped for the immune system process GO terms by utilizing Cytoscape equipped with the ClueGO and CluePedia plugins. Each node represented a GO term, and its size represented the significance. An edge indicates the existence of common genes: a finer line indicates a smaller overlap. The different functional groups of GO terms were reflected by different node colors and are shown on the pie chart. . The detailed statistics are provided in Table S4.

2.2.4. KEGG Pathway Enrichment Analysis

KEGG pathway analysis was conducted for the purpose of predicting the potential functions of the 40 cotargets. According to a kappa score of ≥0.4, 166 terms were found to be connected by 1247 edges (Figure 8(a)). The most significant KEGG terms for the target genes included the IL-17 signaling pathway, hepatitis C, Toll-like receptor signaling pathway, Yersinia infection, Epstein-Barr virus infection, C-type lectin receptor signaling pathway, toxoplasmosis, NOD-like receptor signaling pathway, influenza A, malaria, Measles, and inflammatory bowel disease. Beyond that, we noticed that several inflammation-related terms were also significantly enriched, such as the JAK-STAT signaling pathway, Th17 cell differentiation, NF-kappa B signaling pathway, RIG-I-like receptor signaling pathway, and HIF-1 signaling pathway (Figure 8(b)). These results were in accordance with those of the KEGG enrichment analysis for all the targets of SZS, suggesting that interventions in the inflammatory process might be the crucial mechanism involved in the treatment by SZS of the 2019 novel coronavirus or SARS-CoV pneumonia.

(a)

(b)

Figure 8

KEGG enrichment analysis and pathway mapping. (a) Functionally grouped network of the enriched categories was generated for the target genes using ClueGO and CluePedia plugins. Pathway terms were represented as nodes, and the node size represents the significance of the enrichment of the term. The pie chart shows the proportion of each group associated with 40 cotargets. (b) The bars show the number of genes related to the pathway terms. . The detailed statistics are provided in Table S5.

3. Discussion

3.1. SZS Is a Traditional TCM Formula Which Has Been Adopted to Treat COVID-19 in China

The manifestations of COVID-19 mimic those of SARS-CoV in terms of symptomatology, laboratory features, and chest imaging manifestations [1, 13–15], while there were still characteristics useful for differential diagnosis [14, 15, 36]. The most distinct symptoms of hospitalized COVID-19 patients were subsequent respiratory insufficiency caused by viral pneumonia, even followed by ARDS or acute respiratory failure [1, 13–15]. Combined with life support, TCM treatment has effectively reduced the severity and enhance the recovery of the COVID-19 in China [17, 18]. TCM already showed therapeutic effects on SARS during 2002–2003 [41–46]. Several TCM prescriptions have already been formulated for COVID-19 treatment, such as Ma Xing Shi Gan Decoction, Qing-Fei-Pai-Du-Tang, Lian Hua Qing Wen Capsule, and Shuang Huang Lian [47, 48], among which a classic formula named Shen Zhu San was regarded as one of the basic prescriptions [17–19, 49]. It is worth noting that different TCM prescriptions are suitable for COVID-19 cases with different symptoms. For example, Chinese doctors applied Ma Xing Shi Gan Decoction to treat COVID-19 patients who have a high-grade fever, while SZS is suitable for COVID-19 cases with a slight fever and body aches. SZS contains seven herbs that are made into a decoction that has been used to treat epidemic infectious diseases since the Song Dynasty in China. Some of the herbs and the compounds in SZS have been proven to have antiviral or anti-inflammatory influences in vitro or in vivo [50–55]. The database-based approaches are effective and seem to become prevalent in investigating the mechanism of TCM in treating COVID-19 [56, 57]. Our previous research utilized such approaches to systematically reveal the material basis and underlying biological mechanisms of Ma Xing Shi Gan Decoction in the treatment of COVID-19 [49]. In the present study, we utilized integrating network pharmacological approaches to investigate the potential mechanism of SZS involved in treating COVID-19.

The clinical manifestations, radiological evidence, and laboratory result abnormalities all hinted that the pulmonary pathological features were closely related to the inability of the lung to exchange gas, which is due to viral infection and persistent immune system responses, consequently resulting in downstream airway obstruction and alveolar structure loss [14, 15, 36, 58, 59]. Targeting the immune pathways that promote the inflammatory signals and contribute to diffuse alveolar injury may help to maintain vitality and earn a valuable time window for clearing the virus. Simultaneously, TCM therapy, such as SZS, could serve as an important adjuvant therapy [58, 60, 61]. Our data indicated that immune cells and tissues were the main targets of SZS (Figure 1), and the biological effect of SZS was widely involved in the regulation of the immune process, as illustrated in Figure 7(b).

3.2. SZS May Relieve COVID-19 through Suppressing Cytokine Storm

Resulting from virus-induced abnormal immune activation primarily, cytokine storm is an acute, lethal, and systemic complication connected with the severity of COVID-19 [1, 62]. Abnormal concentrations of several associated cytokines like TNF, IL-10, and IL-6 were found in COVID-19 patients [1, 63, 64]. As shown in Figure 9, SZS was likely to regulate inflammatory factor production through different mechanisms. Specifically, the TNF pathway was enriched in the KEGG enrichment analysis (Table S3), and TNF was a targeted protein in the PPI network (Figure 2(a)). TNF, a famous and perhaps the most intensely studied proinflammatory cytokine, is considered to play a prominent role in the cytokine storm and has been identified as a central cytokine in acute viral diseases, including those caused by the influenza virus, dengue virus, and Ebola virus [65, 66]. Besides, TNF is induced by secondary mediator cascades in respiratory epithelial cells after influenza virus infection, as opposed to a range of proinflammatory cytokines directly induced in respiratory epithelial cells [67, 68], which further indicates that TNF may act as a critical effector in amplifying and broadening the proinflammatory response and then escalating the cytokine storm [68, 69]. IL-10 is another core component of the cytokine system to prevent immunopathology during inflammatory responses [70]. Through the extracellular signal-regulated kinase/IL-10 axis, IL-10 is induced and stimulates STAT3, one of the 31 hub targets (Figure 2(b)), to inhibit the NF-κB signaling pathway in the early phase of infection. Proinflammatory cytokines such as TNF are downregulated, and the cytokine storm can thereby be alleviated. IL-10 dysregulation tends to injure the host because of either excessive pathogen proliferation or cytokine storm [71]. Although the specific effects of IL-6 are mixed, IL-6 has a predominantly anti-inflammatory effect to protect the host from diverse infections and tissue injuries [72]. However, high levels of IL-6 may exert disproportionate effects like activating the coagulation pathway and vascular endothelial cells (ECs), thereby generating acute cytokine storms [62].

Several ingredients of SZS, such as quercetin, kaempferol, and wogonin, have been proven to reduce TNF and IL-6, thus suppressing inflammatory activity [65, 73–76]. What is more, quercetin could also regulate the cytokine balance by increasing serum IL-10 [77]. From these results and literature review, we came to the assumption that by decreasing TNF, quercetin, kaempferol, and wogonin could inhibit the dysfunction of ECs, extravasation of cytokines and inflammatory cells, and injury of the parenchymal cell. The ingredients in SZS may decrease tissue factor (TF) (Table S1) activated by excessive IL-6, thus decreasing the incidence of a lethal cytokine storm through suppressing the coagulation pathway. Besides, quercetin may also protect the host from cytokine storms by activating STAT3 to inhibit the NF-κB/TNF-α axis. To be brief, quercetin, kaempferol, and wogonin in SZS potentially treat COVID-19 by suppressing cytokine storm.

3.3. Understanding the Role of SZS in COVID-19 by Analogy to SARS Treatment: SZS Could Participate in Regulating Immune Response and Mediating Cell Death and Survival

In the present study, due to the similarities between SARS-CoV pneumonia and COVID-19 from aspects of the etiology, clinical manifestation, and laboratory findings, we also attempted to understand SZS by analogy to the mechanisms of SZS in relieving SARS-CoV pneumonia. As illustrated in Figure 4, SZS shared 40 common targets associated with SARS. The acute phase of SARS is related to a severe reduction in the number of T cells in the blood [78].

Jun and Fos are contained in the 21 hub cotargets (Figure 6(b)) as well as the enriched GO terms of the immune system process (Figure 7(b)). c-Jun NH2-terminal kinases (JNKs) take part in cellular proliferation, differentiation, and death, as well as the response to proinflammatory molecules. In particular, they play a critical role in T helper cell (Th) activation and maintenance of Th1/Th2 polarization [79]. Fos regulates transcription of several cytokine genes affected by SARS-CoV infection, including the proinflammatory cytokines involved in the cytokine storm during SARS-CoV pathogenesis [80]. AP-1 consists of members of Jun and Fos families, and the activity of AP-1 is regulated by JNK and extracellular signal-regulated kinase (ERK), of which the activation induces transcription of c-Fos and succeeding AP-1 activation. Through the activation of JNK and ERK signaling cascades, SARS induces an increase in the activity of AP-1, which regulates the transcription of many cytokine genes affected in SARS-CoV infection [80]. As Jun was targeted by wogonin, beta-sitosterol, kaempferol, and formononetin while Fos was targeted by quercetin, we suggest that the immunomodulatory effects of SZS may be mediated by Jun and Fos to protect against cytokine storm for COVID-19 treatment.

In addition to the potential therapeutic mechanism mentioned above, several targets of SZS, such as STAT3 and MAPK, were closely related to apoptosis (Figures 2 and 3). The transcriptional activity of STAT3 is influenced by viral infection and is relevant to the antiapoptotic activity of cells, and MAPK is reported to promote both cell death and survival [81, 82]. Moreover, Akt was the target with the highest degree among the 31 hub genes (Figure 2(b)), and the Akt signaling pathway is involved in cellular prosurvival signaling, which has been proven to be involved in the development of SARS-CoV infection [83]. The potential impact of SZS on cell death and survival during COVID-19 deserves further exploration.

3.4. SZS May Treat COVID-19 via Improving Respiratory Function

As indicated in the pathological results of COVID-19 [84], the diffuse alveolar damage with fibrin rich hyaline membranes and multinucleated giant cells, triggered by virus infection and inflammatory dysregulation, was the leading cause of diminished respiratory function like dyspnea and ARDS [84, 85]. Resulting from alveolar endothelial and epithelial barriers disruption, elevated vascular permeability ultimately leads to the dysfunction of the pulmonary alveolar-capillary barrier [86, 87]. After virus infection, the activation of the phosphoinositide 3-kinase (PI3K)/Akt/eNOS signaling pathway leads to an increase in histamine and NO, which increases microvascular permeability [88]. Sesamin is a lignan that exhibits antioxidative and anti-inflammatory effects by blocking the Toll-like receptor 4 signaling pathway and protecting human ECs in the lung [89]. Sesamin seems to alleviate histamine and NO-induced leakage by potently inhibiting Akt and p38 (MAPK) activities [90]. Moreover, we noticed that the histamine H1/2/3 receptors were simultaneously targeted by several compounds in SZS, including alloisoimperatorin, cnidilin, prangenidin, prangenin, and caribine, which probably hints at the potential role of SZS in regulating histamine-mediated microvascular permeability. Owing to the excessive fibronectin inhibited by the Serpine1 driven tissue-type plasminogen activators (PLAT), the elevated vascular permeability stimulates inflammatory cell migration and proliferation and recruits neutrophils to the lung, thus leading to severe lung injury [91]. Several compounds in SZS targeted Serpine1 and PLAT, such as quercetin, inflacoumarin A, 4,9-dimethoxy-1-vinyl-β-carboline, and kanzonol F (Table 1). Quercetin has been proven to increase the expression of PLAT while suppressing Serpine1 [92, 93]. Meanwhile, several active ingredients of SZS also targeted the Serpine1 and PLAT, including inflacoumarin A, 4,9-dimethoxy-1-vinyl-β-carboline, and kanzonol F (Table 1). In summary, SZS may play a protective role in maintaining the alveolar endothelial and epithelial barriers from disruption, relieving pulmonary capillary leakage, inhibiting the inflammatory environment of lung, and ultimately improving respiratory function.

There were some limitations in the present work. Since the latest update time of TCMSP and SymMap was 2014 and 2018, respectively, there may be some omissions of the components of SZS. Also, the dose-effect relationship of each component in SZS is not involved in these databases [94]. Such limitations remind us of treating the network pharmacological results more cautiously. At present, nevertheless, the network pharmacological approaches are still beneficial to research the mechanism of TCM in treating COVID-19.

4. Conclusion

In the present study, we investigated the underlying mechanism of the effects of SZS in treating the COVID-19 by utilizing the methods of integrating network pharmacology. The therapeutic effects potentially involve suppressing cytokine storms, protecting the pulmonary alveolar-capillary barrier, regulating the immune response, and mediating cell death and survival. Further experiments are needed to validate the specific molecular mechanisms of the effects of SZS on COVID-19.

5. Methods

5.1. Constructing Database of Candidate Compounds

The constituents of the compounds in SZS were retrieved from the online public databases Traditional Chinese Medicine Systems Pharmacology (TCMSP) (http://lsp.nwu.edu.cn/tcmsp.php) and SymMap (https://www.symmap.org/) [95–98]. Each candidate’s druggability was analyzed according to its oral bioavailability (OB) and drug-likeness (DL). As an important indicator to objectively evaluate the intrinsic quality of drugs, OB is the degree and speed of the absorption of drugs into the circulatory system. The higher the OB of the compound is, the more likely the compound is to be applied clinically. DL refers to the sum of the pharmacokinetic and safety properties derived from the interactions of the physicochemical properties and structural factors, including solubility, permeability, and stability. The molecules with OB ≥ 30% and DL ≥ 0.18 were considered to exhibit relatively good pharmacokinetic properties and were identified as candidate compounds for analysis [99, 100]. All of the compound information was standardized according to “Canonical SMILES” based on the PubChem database (https://pubchem.ncbi.nlm.nih.gov/).

5.2. Predicting Targets of Active Compounds

The validated targets, which were confirmed by experiments, of the active ingredients of SZS were obtained from the TCMSP database [101]. Then, the potential target prediction of these active ingredients was performed by utilizing ChemMapper (http://www.lilab-ecust.cn/chemmapper/index.html), a web server for predicting potential drug targets based on the 3D structure similarity [102–105]. The similarity score is calculated from the molecular 3D similarity between the query and hit compounds and is scaled to [0, 2], while only five thresholds (0.8, 1.0, 1.2, 1.5, and 1.8, with 1.2 as the default value) can be set. All the compounds of which the score is below the threshold will not be exhibited. The closer the similarity score is to 2, the more likely a pharmacological link between the molecules is. In the present research, the similarity score threshold was set as the default value of 1.2 compared with drugs in the DrugBank database. Table S1 listed the targets of the active ingredients of SZS. As illustrated in the last column of Table S1, “validated” meant that the interactions of ingredients and targets were validated by the TCMSP database, while the others, which were shown in the form of prediction scores, were obtained from the ChemMapper. The prediction score is a standard score normalized to [0, 1]. The closer the prediction score is to 1, the higher confidence between the ingredient and the target is. Commonly, to screen the potential targets of the active compounds, the thresholds of prediction scores are set above 0, indicating good confidence [103–107]. To be specific, the prediction scores were outputted with three decimal places. If the score were still 0.000, which meant a very low level of confidence, the target would be identified as meaningless and removed. Moreover, duplicates of the validated and predicted targets were eliminated. In this way, the potential and rational targets of SZS could be taken into further analyses as far as possible, especially when the pathogenesis and therapeutic mechanism of COVID-19 are still not clear. All the targets obtained above were standardized according to their gene names and UniProt IDs by searching the UniProtKB (https://www.uniprot.org/) database with “Homo sapiens” as the species [108].

5.3. Constructing the Target Tissue Location Network

BioGPS is a database for querying and organizing gene annotation resources. It provides gene expression data for cells or tissues obtained by microarray analysis [32]. The target-organ location network was constructed with the use of the dataset: GeneAtlas U133 A, gcrma (http://biogps.org/#goto=welcome). First, the mRNA expression patterns of each target gene were determined in 84 different organ tissues. Second, the average values for each gene were analyzed. Third, the genes were located in the relevant organ tissues where the mRNA expression level was higher than the mean. Finally, a target-organ location network was constructed using Cytoscape 3.7.2 (https://cytoscape.org/), a well-known tool for network pharmacology research, to visualize the biological pathways and intermolecular interaction networks, among others. Furthermore, this supplied a basic set of features for data integration, analysis, and visualization for complex network analysis [109].

5.4. Constructing the Protein-Protein Interaction (PPI) Network

Since the chances of proteins performing their assigned functions individually are small, proteins involved in the biochemical process in the same cell tend to form macromolecular complexes through interactions to perform biological functions. Therefore, the exploration of protein interactions and their interaction networks are vital to understanding cellular organization, bioprocesses, and functions. In order to better understand protein interactions systematically, the associated targets were inputted into Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, Version: 11.0) (https://string-db.org/) to obtain the relevant information on protein interaction, for which the genes were determined as nodes and the interactions as lines in a network [34]. STRING is an online database of functional protein association networks, providing associations between proteins according to experimentally determined data, literature mining, databases, and gene associations like neighborhood, fusions, co-occurrence, coexpression, or protein homology [110]. To ensure high interaction confidence, the minimum score was set to the highest confidence value of 0.9 to include the broad scope of protein interactions. Besides, disconnected proteins in the network were excluded. Finally, the protein interactions were determined according to exported statistics. The degree of a node indicates the number of routes to which the node is connected.

5.5. Gene Ontology (GO) Enrichment Analysis

Gene ontology (GO) is an international standardized gene functional classification system that offers a dynamically updated controlled vocabulary and strictly defined concepts to comprehensively describe the properties of genes and their products in any organism [35]. GO has three ontologies: molecular function (MF), cellular component (CC), and biological process (BP). The basic unit of GO is the GO term. Each GO term belongs to a type of ontology. GO enrichment analysis provides all GO terms that are significantly enriched in targets compared to the genome background and filters out the targets that correspond to biological functions. Firstly, all targets were mapped to GO terms in the gene ontology database (http://www.geneontology.org/), the gene numbers were calculated for every term, and significantly enriched GO terms in targets compared to the genome background were defined by a hypergeometric test. The calculated value was then subjected to FDR correction. In this study, GO terms with an FDR ≤0.01, as the threshold, and the data were collected by ClueGO and CluePedia (Cytoscape plugins) [111–113]. Moreover, the targets were mapped for biological process (GO:0002376) and immune system process to explore the underlying impact on immune responses by utilizing Cytoscape equipped with the ClueGO and CluePedia plugins [38–40].

5.6. Pathway Enrichment Analysis

Genes usually interact with each other to play roles in certain biological functions. The pathway-based analysis helps to further understand gene biological functions. KEGG is the major public pathway-related database [114]. Pathway enrichment analysis identified significantly enriched metabolic pathways or signal transduction pathways involving targets compared with the whole genome background in the KEGG pathway database (http://www.genome.jp/kegg/). The calculated value was subjected to FDR correction by taking FDR ≤0.01 as a threshold and the data were collected by the ClueGO and CluePedia plugins [115, 116]. Moreover, R software version 3.6.1 (http://www.r-project.org) with several R packages, including clusterProfiler, org.Hs.eg.db, enrichplot, and ggplot2, was applied to draw graphs during the GO and pathway enrichment analyses, such as bar plots, bubble diagrams, and a KEGG pathway map [117]. The R packages are available on Bioconductor (https://www.bioconductor.org/) [118].

5.7. Screening of Genes Related to SARS

Genes related to SARS were derived from MalaCards (https://www.malacards.org/pages/info), a database of human genes and genetic disorders. The proteins acquired from MalaCards were used as hub proteins and submitted to STRING and Human Protein Reference Database (HPRD) (http://www.hprd.org/) to determine the proteins interacting with these hub proteins. HPRD is a resource for experimentally derived information about the human proteome including protein-protein interactions, posttranslational modifications, and tissue expression [119, 120]. Meanwhile, we identified the gene expression data regarding SARS from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/), an international public repository that archives and freely distributes high-throughput gene expression and other functional genomic datasets. The GEO dataset numbered GSE1739 is comprised of gene expression data obtained using the GPL201 Affymetrix Human HG-Focus Target Array from 14 samples, including 10 peripheral blood samples from adult patients with SARS and 4 controls. The workflow has been wrapped into the R packages “limma” and “edgeR” to identify the differentially expressed genes (DEGs) between the case and control samples [118]. In the current study, genes with an adjusted value <0.05 and |LogFC|>1.0 were regarded as DEGs. Meanwhile, the R package “pheatmap” was used to draw the heatmap plot and volcano plot, with the up- and downregulated genes presented.

5.8. Construction and Analysis of Visualization Networks

In order to understand the complex relationship between active compounds and potential targets, the compound-target networks as well as the relationships among the compound-target-disease interactions were established with the software Cytoscape-v3.7.2 to reveal the molecular mechanisms. This network was composed of nodes and edges. Nodes represent molecules (compounds, targets, pathways, herbs, or diseases), and edges indicate intermolecular interactions, namely, the connections between nodes.

Abbreviations

AJ:	Adherens junctions
ARDS:	Acute respiratory distress syndrome
BP:	Biological process
CC:	Cellular component
COVID-19:	Coronavirus disease 2019
DEGs:	Differentially expressed genes
DL:	Drug-likeness
ECM:	Extracellular matrix
GEO:	Gene Expression Omnibus
GO:	Gene ontology
HPRD:	Human Protein Reference Database
JNK:	c-Jun NH2-terminal kinases
MF:	Molecular function
MMPs:	Matrix metalloproteinases
OB:	Oral bioavailability
PI3K:	Phosphoinositide 3-kinase
PLAT:	Plasminogen activators
PPI:	Protein-protein interaction
ROS:	Reactive oxygen species
SARS:	Severe acute respiratory syndrome
SZS:	Shen Zhu San
TCM:	Traditional Chinese medicine
TCMSP:	Traditional Chinese Medicine Systems Pharmacology
Th:	T helper cell
TJ:	Tight junctions
TNF-α:	Tumor necrosis factor-α
USA:	The United States of America
VCAM-1:	Vascular cell adhesion molecule 1
WHO:	World Health Organization.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request..

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

W. M. H., W. D. P., and L. Z. G. conceptualized the study; W. Y. X. and R. Y. H. developed the methodology; Z. G. W. was responsible for software; M. J. R. and W. S. Q. performed the formal analysis; W. Y. X. and M. J. R. were involved in investigation; Z. C. Y., S. X. H., and Z. Y. Q. were responsible for resources; L. Z. G. was involved in data curation; W. Y. X., R. Y. H., M. J. R., W. S. Q., and Z. C. Y. were involved in writing and preparing the original draft; S. M. Z., W. S. Q., L. H., and Z. Y. were involved in reviewing and editing; W. M. H. and W. D. P. were involved in supervision; W. Y. X. was involved in project administration; W. M. H., W. D. P., and L. Z. G were responsible for funding acquisition. All authors have read and agreed to the published version of the manuscript. Yuxuan Wang, Yuhua Ru, and Guowei Zhuo contributed equally to this work.

Acknowledgments

The authors gratefully appreciate all the people fighting with the COVID-19 epidemics in the world. This research was funded by grants from the Evidence-Based Capacity Building Project (2019XZZX-NB007), Jiangsu Province Administration Project (2017ZX09304021), Jiangsu Medical Outstanding Talents Project (JCRCA2016002), and Jiangsu Provincial Key Medical Center (YXZXA2016002) to W. D. P. and the National Natural Science Foundation of China (81673412) to L.Z.G.

Supplementary Materials

Table S1: a database of component targets . Table S2: the statistics of GO functional enrichment of 627 targets. Table S3: the statistics of KEGG enrichment of 627 targets. Table S4: the statistics of GO functional enrichment of 40 cotargets in the immune system process. Table S5: the statistics of KEGG enrichment of 40 cotargets. (Supplementary Materials)

References

C. Huang, Y. Wang, X. Li et al., “Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China,” The Lancet, vol. 395, no. 10223, pp. 497–506, 2020.
View at: Publisher Site | Google Scholar
F.-S. Wang and C. Zhang, “What to do next to control the 2019-nCoV epidemic?” The Lancet, vol. 395, no. 10222, pp. 391–393, 2020.
View at: Publisher Site | Google Scholar
The State Council Information Office PRC, “SCIO briefing on the joint response to the new coronavirus-related pneumonia outbreak,” 2020, http://english.scio.gov.cn/pressroom/2020-01/29/content_75657018.htm.
View at: Google Scholar
World Health Organization, “WHO director-general’s statement on ihr emergency committee on novel coronavirus (2019-nCoV),” 2020, https://www.who.int/dg/speeches/detail/who-director-general-s-statement-on-ihr-emergency-committee-on-novel-coronavirus-(2019-ncov).
View at: Google Scholar
T. N. Lerner, L. Ye, and K. Deisseroth, “Communication in neural circuits: tools, opportunities, and challenges,” Cell, vol. 164, no. 6, pp. 1136–1150, 2016.
View at: Publisher Site | Google Scholar
A. Gorbalenya, S. Baker, R. Baric et al., Severe Acute Respiratory Syndrome-Related Coronavirus: the Species and its Viruses, a Statement of the Coronavirus Study Group, 2020.
Global Biodefense, “WHO names disease caused by new coronavirus: COVID-19. 2020,” https://globalbiodefense.com/2020/02/11/who-names-disease-caused-by-new-coronavirus-covid-19/.
View at: Google Scholar
P. Zhou, “A pneumonia outbreak associated with a new coronavirus of probable bat origin,” Nature, vol. 579, 2020.
View at: Google Scholar
H. Wang, X. Li, T. Li et al., “The genetic sequence, origin, and diagnosis of SARS-CoV-2,” European Journal Of Clinical Microbiology & Infectious Diseases: Official Publication of the European Society of Clinical Microbiology, vol. 39, no. 9, pp. 1629–1635, 2020.
View at: Publisher Site | Google Scholar
R. Lu, X. Zhao, J. Li et al., “Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding,” The Lancet, vol. 395, no. 10224, pp. 565–574, 2020.
View at: Publisher Site | Google Scholar
M. Letko, A. Marzi, and V. Munster, “Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses,” Nature Microbiology, vol. 5, 2020.
View at: Google Scholar
S. F. Ahmed, A. A. Quadeer, and M. R. McKay, “Preliminary identification of potential vaccine targets for the COVID-19 coronavirus (SARS-CoV-2) based on SARS-CoV immunological studies,” Viruses, vol. 12, no. 3, 2020.
View at: Publisher Site | Google Scholar
Q. Li, X. Guan, P. Wu et al., “Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia,” The New England Journal of Medicine, vol. 382, no. 13, pp. 1199–1207, 2020.
View at: Google Scholar
N. Chen, Epidemiological and Clinical Characteristics of 99 Cases of 2019 Novel Coronavirus Pneumonia in Wuhan, China: A Descriptive Study, Lancet, London, UK, 2020.
D. Wang, “Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in wuhan, China,” JAMA, vol. 395, 2020.
View at: Google Scholar
L. Lin and T. S. Li, “Interpretation of guidelines for the diagnosis and treatment of novel coronavirus (2019-nCoV) infection by the National Health Commission (trial version 5),” Zhonghua Yi Xue Za Zhi, vol. 100, no. 0, p. E001, 2020.
View at: Google Scholar
Y. G. Wang, “XinXing GuanZhuangBingDu (2019-nCoV) FeiYan ZhongYi LinChuangTeZheng yu BianZhengZhiLiao ChuTan,” Journal of Traditional Chinese Medicine, vol. 4, pp. 281–285, 2020, http://kns.cnki.net/kcms/detail/11.2166.R.20200129.1258.002.html.
View at: Google Scholar
W. K. Zheng, “ZhongYiYao FangZhi XinXing GuanZhuangBingDuGanRan de FeiYan GeDi ZhenLiao FangAn ZongHe FenXi,” Journal of Traditional Chinese Medicine, vol. 4, pp. 277–280, 2020, http://kns.cnki.net/kcms/detail/11.2166.r.20200206.1113.002.html.
View at: Google Scholar
P. Song, L Zhao, X Li et al., “Interpretation of the Traditional Chinese Medicine portion of the diagnosis and treatment protocol for corona virus disease 2019 (Trial Version 7),” Journal of Traditional Chinese Medicine = Chung I Tsa Chih Ying Wen pan, vol. 40, no. 3, pp. 497–508, 2020.
View at: Publisher Site | Google Scholar
H. Peng, Zhongyi Fangji Dacidian, vol. 7, Beijing: People’s Health Publishing House, Beijing, China, 1996.
Chinese Pharmacopoeia Commission, Pharmacopoeia of the People’s Republic of China 2020 Edition, vol. 1, China Medical Science Press, Beijing, China, 2020.
R. A. Isbrucker and G. A. Burdock, “Risk and safety assessment on the consumption of Licorice root (Glycyrrhiza sp.), its extract and powder as a food ingredient, with emphasis on the pharmacology and toxicology of glycyrrhizin,” Regulatory Toxicology and Pharmacology, vol. 46, no. 3, pp. 167–192, 2006.
View at: Publisher Site | Google Scholar
J. T. Azietaku, H. Ma, X.-A. Yu et al., “A review of the ethnopharmacology, phytochemistry and pharmacology of Notopterygium incisum,” Journal of Ethnopharmacology, vol. 202, pp. 241–255, 2017.
View at: Publisher Site | Google Scholar
Z. Chen, C. Zhang, F. Gao et al., “A systematic review on the rhizome of Ligusticum chuanxiong Hort. (Chuanxiong),” Food and Chemical Toxicology, vol. 119, pp. 309–325, 2018.
View at: Publisher Site | Google Scholar
M. Ramalingam and S.-J. Kim, “Phytochemical, toxicological and pharmacological studies of asiasari radix et rhizoma: a review,” Tropical Journal of Pharmaceutical Research, vol. 14, no. 3, pp. 545–554, 2015.
View at: Publisher Site | Google Scholar
T. Ramesh, K. Lee, H.-W. Lee, and S.-J. Kim, “Subacute toxicological evaluation ofAsiasari radixmethanol extract,” Drug and Chemical Toxicology, vol. 32, no. 3, pp. 243–251, 2009.
View at: Publisher Site | Google Scholar
K. Na-Bangchang, T. Plengsuriyakarn, and J. Karbwang, “Research and development of atractylodes lancea (thunb) dc. as a promising candidate for cholangiocarcinoma chemotherapeutics,” Evidence-Based Complementary and Alternative Medicine, vol. 2017, Article ID 5929234, 16 pages, 2017.
View at: Publisher Site | Google Scholar
Y. Wu, “Research advancement on pharmacological effect of radix Angelicae Dahuricae,” Lishizhen Medicine and Materia Medica Research (In Chinese), vol. 20, no. 3, pp. 625–627, 2009.
View at: Google Scholar
M. L. Fernandez and S. Vega-López, “Efficacy and safety of sitosterol in the management of blood cholesterol levels,” Cardiovascular Drug Reviews, vol. 23, no. 1, pp. 57–70, 2005.
View at: Publisher Site | Google Scholar
E. de Wit, N. van Doremalen, D. Falzarano, and V. J. Munster, “SARS and MERS: recent insights into emerging coronaviruses,” Nature Reviews Microbiology, vol. 14, no. 8, pp. 523–534, 2016.
View at: Publisher Site | Google Scholar
B. Liu, Z. Song, J. Yu, P. Li, Y. Tang, and J. Ge, “The atherosclerosis-ameliorating effects and molecular mechanisms of BuYangHuanWu decoction,” Biomedicine & Pharmacotherapy, vol. 123, p. 109664, 2020.
View at: Publisher Site | Google Scholar
C. Wu, X. Jin, G. Tsueng, C. Afrasiabi, and A. I. Su, “BioGPS: building your own mash-up of gene annotations and expression profiles,” Nucleic Acids Research, vol. 44, no. D1, pp. D313–D316, 2016.
View at: Publisher Site | Google Scholar
Z. Song, “Systems pharmacological approach to investigate the mechanism of for Alzheimer’s disease,” Evidence-Based Complementary and Alternative Medicine, vol. 2018, Article ID 5194016, 20 pages, 2018.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–D613, 2019.
View at: Publisher Site | Google Scholar
M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene Ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.
View at: Publisher Site | Google Scholar
J. F. W. Chan, A Familial Cluster of Pneumonia Associated with the 2019 Novel Coronavirus Indicating Person-To-Person Transmission: A Study of a Family Cluster, Lancet, London, UK, 2020.
N. Zhu, D. Zhang, W. Wang et al., “A novel coronavirus from patients with pneumonia in China,” The New England Journal of Medicine, vol. 382, no. 8, pp. 727–733, 2020.
View at: Google Scholar
V. L. Hegde, S. Tomar, A. Jackson et al., “Distinct MicroRNA expression profile and targeted biological pathways in functional myeloid-derived suppressor cells induced by Δ9-tetrahydrocannabinolin vivo,” Journal of Biological Chemistry, vol. 288, no. 52, pp. 36810–36826, 2013.
View at: Publisher Site | Google Scholar
P. Dong, “Identification of key genes and pathways in triple-negative breast cancer by integrated bioinformatics analysis,” BioMed Research International, vol. 2018, Article ID 2760918, 10 pages, 2018.
View at: Publisher Site | Google Scholar
S. Puente-Marin, “In Silico functional networks identified in fish nucleated red blood cells by means of transcriptomic and proteomic profiling,” Genes, vol. 9, no. 4, 2018.
View at: Publisher Site | Google Scholar
Z. Chen and T. Nakamura, “Statistical evidence for the usefulness of Chinese medicine in the treatment of SARS,” Phytotherapy Research, vol. 18, no. 7, pp. 592–594, 2004.
View at: Publisher Site | Google Scholar
Y. Luo, C.-Z. Wang, J. Hesse-Fong, J.-G. Lin, and C.-S. Yuan, “Application of Chinese medicine in acute and critical medical conditions,” The American Journal of Chinese Medicine, vol. 47, no. 6, pp. 1223–1235, 2019.
View at: Publisher Site | Google Scholar
M.-M. Zhang, X.-M. Liu, and L. He, “Effect of integrated traditional Chinese and Western medicine on SARS: a review of clinical evidence,” World Journal of Gastroenterology, vol. 10, no. 23, pp. 3500–3505, 2004.
View at: Publisher Site | Google Scholar
J. L. McKimm-Breschkin, S. Jiang, D. S. Hui, J. H. Beigel, E. A. Govorkova, and N. Lee, “Prevention and treatment of respiratory viral infections: presentations on antivirals, traditional therapies and host-directed interventions at the 5th ISIRV antiviral group conference,” Antiviral Research, vol. 149, pp. 118–142, 2018.
View at: Publisher Site | Google Scholar
P.-C. Leung, “The efficacy of Chinese medicine for SARS: a review of Chinese publications after the crisis,” The American Journal of Chinese Medicine, vol. 35, no. 4, pp. 575–581, 2007.
View at: Publisher Site | Google Scholar
J. Liu, E. Manheimer, Y. Shi, and C. Gluud, “Chinese herbal medicine for severe acute respiratory syndrome: a systematic review and meta-analysis,” The Journal of Alternative and Complementary Medicine, vol. 10, no. 6, pp. 1041–1051, 2004.
View at: Publisher Site | Google Scholar
Y. Yang, M. S. Islam, J. Wang, Y. Li, and X. Chen, “Traditional Chinese medicine in the treatment of patients infected with 2019-new coronavirus (SARS-CoV-2): a review and perspective,” International Journal of Biological Sciences, vol. 16, no. 10, pp. 1708–1717, 2020.
View at: Publisher Site | Google Scholar
L. L. D. Zhong, W. C. Lam, W. Yang et al., “Potential targets for treatment of coronavirus disease 2019 (COVID-19): a review of Qing-Fei-Pai-du-tang and its major herbs,” The American Journal of Chinese Medicine, vol. 48, no. 5, pp. 1051–1071, 2020.
View at: Publisher Site | Google Scholar
Y. X. Wang, J. R Ma, S. Q Wang et al., “Utilizing integrating network pharmacological approaches to investigate the potential mechanism of Ma Xing Shi Gan Decoction in treating COVID-19,” European Review for Medical and Pharmacological Sciences, vol. 24, no. 6, pp. 3360–3384, 2020.
View at: Publisher Site | Google Scholar
Y. Cheng, J.-Y. Mai, T.-L. Hou, J. Ping, and J.-J. Chen, “Antiviral activities of atractylon from Atractylodis rhizoma,” Molecular Medicine Reports, vol. 14, no. 4, pp. 3704–3710, 2016.
View at: Publisher Site | Google Scholar
J. Cinatl, B. Morgenstern, G. Bauer, P. Chandra, H. Rabenau, and H. Doerr, “Glycyrrhizin, an active component of liquorice roots, and replication of SARS-associated coronavirus,” The Lancet, vol. 361, no. 9374, pp. 2045-2046, 2003.
View at: Publisher Site | Google Scholar
L. Wang, R. Yang, B. Yuan, Y. Liu, and C. Liu, “The antiviral and antimicrobial activities of licorice, a widely-used Chinese herb,” Acta Pharmaceutica Sinica B, vol. 5, no. 4, pp. 310–315, 2015.
View at: Publisher Site | Google Scholar
Y. S. Eng, “Unraveling the molecular mechanism of traditional chinese medicine: formulas against acute airway viral infections as examples,” Molecules, vol. 24, no. 19, 2019.
View at: Google Scholar
T. Li and T. Peng, “Traditional Chinese herbal medicine as a source of molecules with antiviral activity,” Antiviral Research, vol. 97, no. 1, pp. 1–9, 2013.
View at: Publisher Site | Google Scholar
C.-L. Chen and D.-D. Zhang, “Anti-inflammatory effects of 81 Chinese herb extracts and their correlation with the characteristics of traditional Chinese medicine,” Evidence-based Complementary and Alternative Medicine, vol. 2014, Article ID 985176, 2014.
View at: Google Scholar
X. Ren, X.-X. Shao, X.-X. Li et al., “Identifying potential treatments of COVID-19 from Traditional Chinese Medicine (TCM) by using a data-driven approach,” Journal of Ethnopharmacology, vol. 258, p. 112932, 2020.
View at: Publisher Site | Google Scholar
B. Yang, H. Wang, C. Lu et al., “A network pharmacology study to uncover the multiple molecular mechanism of the Chinese patent medicine toujiequwen granules in the treatment of corona virus disease (COVID-19),” Research Square, vol. 1, 2020.
View at: Google Scholar
A. C. Kalil and P. G. Thomas, “Influenza virus-related critical illness: pathophysiology and epidemiology,” Critical Care (London, England), vol. 23, no. 1, p. 258, 2019.
View at: Publisher Site | Google Scholar
G. M.-K. Tse, “Pulmonary pathological features in coronavirus associated severe acute respiratory syndrome (SARS),” Journal of Clinical Pathology, vol. 57, no. 3, pp. 260–265, 2004.
View at: Publisher Site | Google Scholar
H. Chen and Q. Du, “Potential Natural Compounds for Preventing SARS-CoV-2 (2019-nCoV) Infection,” 2020.
View at: Google Scholar
J. E. M. Clark D Russell and J. K. Baillie, Clinical Evidence Does Not Support Corticosteroid Treatment for 2019-nCoV Lung Injury, Lancet, London, UK, 2020.
T. Tanaka, M. Narazaki, and T. Kishimoto, “Immunotherapeutic implications of IL-6 blockade for cytokine storm,” Immunotherapy, vol. 8, no. 8, pp. 959–970, 2016.
View at: Publisher Site | Google Scholar
Y. Fang, H. Zhang, Y. Xu, J. Xie, P. Pang, and W. Ji, “CT manifestations of two cases of 2019 novel Coronavirus (2019-nCoV) pneumonia,” Radiology, vol. 295, no. 1, pp. 208-209, 2020.
View at: Publisher Site | Google Scholar
M. Bassetti, A. Vena, and D. R. Giacobbe, “The novel Chinese Coronavirus (2019-nCoV) infections: Challenges for fighting the storm,” European journal of clinical investigation, vol. 50, no. 3, Article ID e13209, 2020.
View at: Google Scholar
H. K. Kim, H. R. Park, J. S. Lee, T. S. Chung, H. Y. Chung, and J. Chung, “Down-regulation of iNOS and TNF-α expression by kaempferol via NF-κB inactivation in aged rat gingival tissues,” Biogerontology, vol. 8, no. 4, pp. 399–408, 2007.
View at: Publisher Site | Google Scholar
J. R. Tisoncik, M. J. Korth, C. P. Simmons, J. Farrar, T. R. Martin, and M. G. Katze, “Into the eye of the cytokine storm,” Microbiology and Molecular Biology Reviews, vol. 76, no. 1, pp. 16–32, 2012.
View at: Publisher Site | Google Scholar
M. C. Chan, “Proinflammatory cytokine responses induced by influenza A (H5N1) viruses in primary human alveolar and bronchial epithelial cells,” Respiratory Research, vol. 6, no. 1, p. 135, 2005.
View at: Publisher Site | Google Scholar
M. M. Levy, “SCCM/ESICM/ACCP/ATS/SIS International Sepsis Definitions Conference,” Critical Care Medicine, vol. 31, no. 4, pp. 1250–1256, 2001.
View at: Google Scholar
P. Osterlund, J Pirhonen, N Ikonen et al., “Pandemic H1N1 2009 influenza A virus induces weak cytokine responses in human macrophages and dendritic cells and is highly sensitive to the antiviral actions of interferons,” Journal of Virology, vol. 84, no. 3, pp. 1414–1422, 2010.
View at: Publisher Site | Google Scholar
J. M. Rojas, M. Avia, V. Martín, and N. Sevilla, “IL-10: a multifunctional cytokine in viral infections,” Journal of Immunology research, vol. 2017, Article ID 6104054, 14 pages, 2017.
View at: Publisher Site | Google Scholar
M.-H. Tsai, C.-H. Chang, R.-K. Tsai et al., “Cross-regulation of proinflammatory cytokines by interleukin-10 and miR-155 in Orientia tsutsugamushi-infected human macrophages prevents cytokine storm,” Journal of Investigative Dermatology, vol. 136, no. 7, pp. 1398–1407, 2016.
View at: Publisher Site | Google Scholar
W. Y. Park, R. B. Goodman, K. P. Steinberg et al., “Cytokine balance in the lungs of patients with acute respiratory distress syndrome,” American Journal of Respiratory and Critical Care Medicine, vol. 164, no. 10, pp. 1896–1903, 2001.
View at: Publisher Site | Google Scholar
A. B. Granado-Serrano, M. Á. Martín, L. Bravo, L. Goya, and S. Ramos, “Quercetin attenuates TNF-induced inflammation in hepatic cells by inhibiting the NF-κB pathway,” Nutrition and Cancer, vol. 64, no. 4, pp. 588–598, 2012.
View at: Publisher Site | Google Scholar
H. Peng, “Role of quercetin on inhibiting TGF-beta1 and TNF-alpha in lung of rats with silicosis,” in Proceedings of the International Conference on Material and Environmental Engineering (ICMAEE 2014), Atlantis Press, Jiujiang, China, 2014.
View at: Google Scholar
J. Liu, X Li, Y Yue, J Li, T He, and Y He, “The inhibitory effect of quercetin on IL-6 production by LPS-stimulated neutrophils,” Cellular & Molecular Immunology, vol. 2, no. 6, pp. 455–460, 2005.
View at: Google Scholar
D. Yi, “Neuroprotective effects of ferulic acid against glutamate-induced neurotoxicity in pc12 cells,” Nutrients, vol. 6, 2008.
View at: Google Scholar
B. Jihong, “The effects of quercetin on serum IL-18 and IL-10 of nonalcoholic fatty liver disease in rats,” no. 1, p. 18, 2013.
View at: Google Scholar
R. Channappanavar, J. Zhao, and S. Perlman, “T cell-mediated immune response to respiratory coronaviruses,” Immunol Res, vol. 59, no. 1-3, pp. 118–128, 2014.
View at: Publisher Site | Google Scholar
N. Arbour, D. Naniche, D. Homann, R. J. Davis, R. A. Flavell, and M. B. A. Oldstone, “c-Jun NH2-terminal kinase (JNK)1 and JNK2 signaling pathways have divergent roles in CD8+ T cell-mediated antiviral immunity,” Journal of Experimental Medicine, vol. 195, no. 7, pp. 801–810, 2002.
View at: Publisher Site | Google Scholar
B. Varshney and S. K. Lal, “SARS-CoV accessory protein 3b induces AP-1 transcriptional activity through activation of JNK and ERK pathways,” Biochemistry, vol. 50, no. 24, pp. 5419–5425, 2011.
View at: Publisher Site | Google Scholar
T. Mizutani, “Tyrosine dephosphorylation of STAT3 in SARS coronavirus-infected Vero E6 cells,” FEBS Letters, vol. 577, no. 1-2, pp. 187–192, 2004.
View at: Publisher Site | Google Scholar
T. Mizutani, “Signal transduction in SARS-CoV-infected cells,” Annals of the New York Academy of Sciences, vol. 1102, no. 1, pp. 86–95, 2007.
View at: Publisher Site | Google Scholar
T. Mizutani, “JNK and PI3k/Akt signaling pathways are required for establishing persistent SARS-CoV infection in Vero E6 cells,” Biochimica et biophysica acta, vol. 1741, no. 1-2, pp. 4–10, 2005.
View at: Publisher Site | Google Scholar
Z. Xu, L. Shi, Y. Wang et al., “Pathological findings of COVID-19 associated with acute respiratory distress syndrome,” The Lancet Respiratory Medicine, vol. 8, no. 4, pp. 420–422, 2020.
View at: Publisher Site | Google Scholar
R. J. Mason, “Pathogenesis of COVID-19 from a cell biology perspective,” European Respiratory Journal, vol. 55, no. 4, 2020.
View at: Publisher Site | Google Scholar
M. A. Matthay, L. B. Ware, and G. A. Zimmerman, “The acute respiratory distress syndrome,” Journal of Clinical Investigation, vol. 122, no. 8, pp. 2731–2740, 2012.
View at: Publisher Site | Google Scholar
M. A. Matthay, “Acute respiratory distress syndrome,” Nature Reviews. Disease Primers, vol. 5, no. 1, p. 18, 2019.
View at: Publisher Site | Google Scholar
A. Di Lorenzo, C. Fernandez-Hernando, G. Cirino, and W. C. Sessa, “Akt1 is critical for acute inflammation and histamine-mediated vascular leakage,” Proceedings of the National Academy of Sciences, vol. 106, no. 34, pp. 14552–14557, 2009.
View at: Publisher Site | Google Scholar
X. Bai, “Sesamin enhances Nrf2-mediated protective defense against oxidative stress and inflammation in colitis via akt and erk activation,” Oxidative Medicine and Cellular Longevity, vol. 2019, Article ID 2432416, 20 pages, 2019.
View at: Publisher Site | Google Scholar
L. C. Li, H. M. Piao, M. Y. Zheng, Z. H. Lin, G. Li, and G. H. Yan, “Sesamin attenuates mast cell-mediated allergic responses by suppressing the activation of p38 and nuclear factor-κB,” Molecular Medicine Reports, vol. 13, no. 1, pp. 536–542, 2016.
View at: Publisher Site | Google Scholar
L. E. Gralinski, “Mechanisms of severe acute respiratory syndrome coronavirus-induced acute lung injury,” mBio, vol. 4, no. 4, 2013.
View at: Publisher Site | Google Scholar
W. J. Lan and X. X. Zheng, “Activity of Ginkgo biloba extract and quercetin on thrombomodulin expression and tissue-type plasminogen activator secretion by human umbilical vein endothelial cells,” Biomedical and Environmental Sciences : BES, vol. 19, no. 4, pp. 249–253, 2006.
View at: Google Scholar
N. C. Olave, M. H. Grenett, M. Cadeiras, H. E. Grenett, and P. J. Higgins, “Upstream stimulatory factor-2 mediates quercetin-induced suppression of PAI-1 gene expression in human endothelial cells,” Journal of Cellular Biochemistry, vol. 111, no. 3, pp. 720–726, 2010.
View at: Publisher Site | Google Scholar
Y.-X. Huang, W.-X. Wang, S. Zhang, Y.-P. Tang, and S.-J. Yue, “The database-based strategy may overstate the potential effects of traditional Chinese medicine against COVID-19,” Pharmacological Research, vol. 159, p. 105046, 2020.
View at: Publisher Site | Google Scholar
J. Ru, “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, p. 13, 2014.
View at: Publisher Site | Google Scholar
P. Wang, “Intermodule coupling analysis of Huang-Lian-Jie-Du decoction on stroke,” Front Pharmacology, vol. 10, p. 1288, 2019.
View at: Publisher Site | Google Scholar
L. Chen, Y. Cao, H. Zhang et al., “Network pharmacology-based strategy for predicting active ingredients and potential targets of Yangxinshi tablet for treating heart failure,” Journal of Ethnopharmacology, vol. 219, pp. 359–368, 2018.
View at: Publisher Site | Google Scholar
Y. Wu, F. Zhang, K. Yang et al., “SymMap: an integrative database of traditional Chinese medicine enhanced by symptom mapping,” Nucleic Acids Research, vol. 47, no. D1, pp. D1110–d1117, 2019.
View at: Publisher Site | Google Scholar
W. Guo, “Integrating network pharmacology and pharmacological evaluation for deciphering the action mechanism of herbal formula Zuojin pill in suppressing hepatocellular carcinoma,” Frontier Pharmacology, vol. 10, p. 1185, 2019.
View at: Publisher Site | Google Scholar
X. W. Zhang, “Dissection of pharmacological mechanism of Chinese herbal medicine Yihuo Huatan formula on chronic obstructive pulmonary disease: a systems pharmacology-based study,” Science Reports, vol. 9, no. 1, p. 13431, 2019.
View at: Publisher Site | Google Scholar
L. Song, Y Luo, X Wang et al., “Exploring the active mechanism of berberine against HCC by systematic pharmacology and experimental validation,” Molecular Medicine Reports, vol. 20, no. 5, pp. 4654–4664, 2019.
View at: Publisher Site | Google Scholar
L. Gao, J. Hao, Y.-Y. Niu et al., “Network pharmacology dissection of multiscale mechanisms of herbal medicines in stage IV gastric adenocarcinoma treatment,” Medicine, vol. 95, no. 35, p. e4389, 2016.
View at: Publisher Site | Google Scholar
M. Yu, “Identify the key active ingredients and pharmacological mechanisms of compound XiongShao capsule in treating diabetic peripheral neuropathy by network pharmacology approach,” Evidence-Based Complementary and Alternative Medicine, vol. 2019, Article ID 5801591, 11 pages, 2019.
View at: Publisher Site | Google Scholar
L. Gao, “Molecular targets of Chinese herbs: a clinical study of hepatoma based on network pharmacology,” Scientific Reports, vol. 6, p. 24944, 2016.
View at: Publisher Site | Google Scholar
J. Gong, C. Cai, X. Liu et al., “ChemMapper: a versatile web server for exploring pharmacology and chemical structure association based on molecular 3D similarity method,” Bioinformatics, vol. 29, no. 14, pp. 1827–1829, 2013.
View at: Publisher Site | Google Scholar
L. Gao, “Uncovering the anticancer mechanism of Compound Kushen Injection against HCC by integrating quantitative analysis, network analysis and experimental validation,” Scientific Reports, vol. 8, no. 1, p. 624, 2018.
View at: Publisher Site | Google Scholar
Y. Aisa, K. Yunusi, Q. Chen, and N. Mi, “Systematic understanding of the potential targets and pharmacological mechanisms of acteoside by network pharmacology approach,” Medicinal Chemistry Research, vol. 29, no. 4, pp. 793–801, 2020.
View at: Publisher Site | Google Scholar
T. UniProt Consortium, “UniProt: the universal protein knowledge base,” Nucleic Acid, vol. 46, no. 5, p. 2699, 2018.
View at: Google Scholar
D. Otasek, “Cytoscape Automation: empowering workflow-based network analysis,” Genome Biology, vol. 20, no. 1, p. 185, 2019.
View at: Publisher Site | Google Scholar
C. Ctortecka, V. Palve, B. M. Kuenzi et al., “Functional proteomics and deep network interrogation reveal a complex mechanism of action of midostaurin in lung cancer cells,” Molecular & Cellular Proteomics, vol. 17, no. 12, pp. 2434–2447, 2018.
View at: Publisher Site | Google Scholar
G. Bindea, B. Mlecnik, H. Hackl et al., “ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks,” Bioinformatics, vol. 25, no. 8, pp. 1091–1093, 2009.
View at: Publisher Site | Google Scholar
B. Mlecnik, J. Galon, and G. Bindea, “Automated exploration of gene ontology term and pathway networks with ClueGO-REST,” Bioinformatics, vol. 35, no. 19, pp. 3864–3866, 2019.
View at: Publisher Site | Google Scholar
B. Mlecnik, J. Galon, and G. Bindea, “Comprehensive functional analysis of large lists of genes and proteins,” Journal of Proteomics, vol. 171, pp. 2–10, 2018.
View at: Publisher Site | Google Scholar
M. Kanehisa and S. Goto, “KEGG: Kyoto Encyclopedia of Genes and Genomes,” Nucleic Acids Research, vol. 28, no. 1, pp. 27–30, 2000.
View at: Publisher Site | Google Scholar
J. Sun, “To unveil the molecular mechanisms of Qi and blood through systems biology-based investigation into Si-Jun-Zi-Tang and Si-Wu-tang formulae,” Scientific Reports, vol. 6, p. 34328, 2016.
View at: Publisher Site | Google Scholar
Y. Wang, “Systems pharmacology to investigate the interaction of berberine and other drugs in treating polycystic ovary syndrome,” Scientific Reports, vol. 6, Article ID 28089, 2016.
View at: Publisher Site | Google Scholar
G. Yu, L.-G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “Limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, p. e47, 2015.
View at: Publisher Site | Google Scholar
T. S. Keshava Prasad, R. Goel, K. Kandasamy et al., “Human Protein Reference Database--2009 update,” Nucleic Acids Research, vol. 37, pp. D767–D772, 2009.
View at: Publisher Site | Google Scholar
G. Xie, “A network pharmacology analysis to explore the effect of on traumatic brain injury,” BioMed Research International, vol. 2018, Article ID 3951783, 13 pages, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Yuxuan Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1003

Downloads

1017

Citations