Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 17 April 2020
Sec. Inflammation Pharmacology
This article is part of the Research Topic Multiple Pharmacological and Toxicological Approaches to Explore Efficacy and Safety of Plant-derived Active Ingredients View all 8 articles

Kunxian Capsule for Rheumatoid Arthritis: Inhibition of Inflammatory Network and Reducing Adverse Reactions Through Drug Matching

Yujun Tang&#x;Yujun TangYi Zhang&#x;Yi ZhangLin LiLin LiZhijun XieZhijun XieChengping Wen*Chengping Wen*Lin Huang*Lin Huang*
  • College of Basic Medical Science, Zhejiang Chinese Medical University, Hangzhou, China

Tripterygium wilfordii Hook.f and Tripterygium hypoglaucum (H.Lév.) Hutch is effective herbs to prevent aggravation of Rheumatoid arthritis (RA). However, both of them show severe side effects in the reproductive system and other systems. Kunxian Capsule (KX), a Traditional Chinese Medicine (TCM) patent prescription, comprised of 4 herbs, including H.Lév. Hutch, is reported to be an available prescription in treating RA with fewer side effects as compares to Tripterygium tablets. To reveal the pharmacological mechanism of KX in RA treatment and side effect alleviation, we collected related information of KX from open-access databases and performed various analyses. 1354 targets were identified in KX. These targets were enriched in the calcium signaling pathway, cAMP signaling pathway, cGMP-PKG signaling pathway and PI3K-AKT signaling pathway, forming biological functions, such as cofactor binding, coenzyme binding, etc. These pathways or functions mostly affect cell cycle, differentiation, and maturation of Th17 cells, macrophage, and synovial fibroblast. These targets also act on the IL-17 signaling pathway, Th17 cell differentiation signaling pathway and TNF signaling pathway, which is related to inflammation response inhibition. Next, a disease network was constructed, which indicated IMPDH2, MTHFD1 are the key genes answering for the side effects of H.Lév. Hutch. The side effect–related genes lead to the negative regulation of nucleic acid, which could be restored by the rest 3 herbs through some positive amino acid metabolism. In conclusion, KX is a relatively safe alternative approach in RA intervention.

Introduction

Rheumatoid arthritis (RA) is a chronic autoimmune disease characterized by symmetry and peripheral polyarthritis. Epidemiology shows that the prevalence of RA is 0.5% to 1%, varying by geography and population (Silman and Pearson, 2002; Alamanos et al., 2006). Chronic inflammatory arthritis often leads to joint damage and disability as the disease progressed.

Currently, treatment for RA includes non-steroidal anti-inflammatory drugs, glucocorticoids, immunosuppressants, and alternative therapy (such as Chinese medicine). KunXian Capsule (KX) is a Traditional Chinese Medicine (TCM) patent prescription which is widely used in RA treatment. There are four main herbs in KX, namely Tripterygium hypoglaucum (H.Lév.) Hutch (Kun Ming Shan Hai Tang), Epimedium brevicornu Maxim (Yin Yang Huo), Cuscuta chinensis Lam (Tu Si Zi), and Lycium barbarum L (Gou Qi Zi). Previous research has demonstrated that KX can down-regulate chemokines, inhibit inflammation response and has analgesic and immunosuppressive effects of the whole-course intervention (Jia and Liu, 2016). However, the molecular mechanism still needs to be revealed. A systematic review of 528 patients in 7 RCTs showed that the total effective rate of KX combined with methotrexate group is higher than methotrexate alone group in RA treatment (Zhou et al., 2016), the advantage lies in controlling systemic inflammation and reducing the number of swollen joints and morning stiffness. As meta-analysis reported, the adverse reactions of KX, including reproductive toxicity, liver damage, renal dysfunction, leukopenia, gastrointestinal reactions, but the adverse reaction rate is significantly lower than the Tripterygium tablet (Liu et al., 2017; Cao et al., 2018). Thus, we suspect that the drug combination of herbs in KX may work on reducing side effects.

In 2008, Hopkins proposed the concept of “network pharmacology”. Network pharmacology could provide a new strategy for drug development by analyzing the intervention of drugs on disease networks (Hopkins, 2008). It is difficult to make a detailed and comprehensive study of TCM compounds due to multi-component and multi-target characterization. Network pharmacology or system biology provides new methods for the understanding of complex Chinese medicine pharmacological mechanisms (Hopkins, 2008). At present, many databases made a significant contribution to the development of systems biology, such as the Bioinformatics Analysis Tool for Molecular mechanism of TCM (BATMAN-TCM) developed by Liu (Liu et al., 2016). BATMAN-TCM was applied in the study exploring the mechanism of Curcuma longa L (E’Zhu) in the treatment of breast cancer (Kong et al., 2017). In this article, we try to explore the pharmacological mechanism of KX in treating RA and explain how the drugs combination of herbs works on reducing adverse reactions. At the same time, we predicted potential therapeutic targets that guide in-depth research.

Results

The Biological Function of KX

1385 targets of KX were collected from TCM-related databases, and 784 of them come from H.Lév. Hutch. 28635 RA-related genes were gathered from CTD. The number of overlapped targets between KX and RA is 1354, which were defined as anti-RA targets. 50 key genes were identified by CytoHubba, which were defined as hub anti-RA targets. The hub anti-RA targets included AGT, GNG2, ANXA1, etc. (Supplementary Materials 1). Gene Ontology (GO) analysis result showed that the main functions of anti-RA targets were: cofactor binding, coenzyme binding, organic acid binding, carboxyl binding, vitamin binding, etc. (Figure 1A) while the hub anti-RA targets were enriched in G protein related biological function. Such as G protein-coupled amine receptor activity, neurotransmitter receptor activity, catecholamine binding, G protein-coupled amine receptor binding, adrenergic receptor activity, G protein coupling Neurotransmitter activity, G-protein coupled peptide activity, G protein beta/gamma sub-complex complex (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1 Gene ontology (GO) analysis of Kunxian Capsule. (A) The GO analysis of all anti-RA targets of Kunxian Capsule. (B) The GO analysis of hub anti-RA targets of Kunxian Capsule.

We also applied KEGG pathway enrichment analysis for anti-RA targets and hub anti-RA targets. The analysis indicated that they are overlapped in the calcium signaling pathway, cAMP signaling pathway, neural ligand interaction, and cGMP-PKG signaling pathway (Figures 2A, B). Besides, KX is associated with various inflammatory pathways and cell cycle pathways, such as IL-17 signaling pathway, TH17 cell differentiation, TNF signaling pathway, PI3K-AKT signaling pathway. (Supplementary Figures S1–S6).

FIGURE 2
www.frontiersin.org

Figure 2 Kyoto encyclopedia of genes and genomes (KEGG) analysis of Kunxian capsule. (A) The KEGG analysis of all anti-RA targets of Kunxian capsule. (B) The KEGG analysis of hub anti-RA targets of Kunxian capsule.

Protein-protein interaction (PPI) network was carried out to identify the key targets of KX in RA treatment. The whole network contains 1348 nodes, 7164 edges. For better analysis in adverse reactions and presentation, we select the nodes whose degree is above 10 to construct a PPI network (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3 Protein protein interaction network. The yellow nodes represent targets in three herbs only, the blue nodes represent targets in Tripterygium hypoglaucum hutch only, the pink nodes represent the shared targets.

Mechanism of Adverse Reactions

Network Construction of Adverse Reactions

From the systematic reviews mentioned above (Zhou et al., 2016), we identified 4 adverse reactions: reproductive toxicity, liver damage, renal dysfunction, leukopenia. The adverse responses involve a total of 136 genes, 9 of which were overlapped with the targets of H.Lév. Hutch, i.e., IMPDH2, MTHFD1, PKD2, PTGIS, SLC25A13, AGTR1, NR0B1, AR, and AGT (Figure 4). IMPDH2 is involved in cellular guanine metabolism, which could play a critical role in DNA and RNA synthesis (Dalal et al., 2009). MTHFD1 could encode a protein that possesses 5,10-methenyltetrahydrofolate cyclohydrolase, 5,10-methylenetetrahydrofolate dehydrogenase, and 10-formyltetrahydrofolate synthetase activities. Each of these enzymatic activities catalyzes one of three sequential reactions in the transformation of 1-carbon derivatives of tetrahydrofolate (THF), which are substrates for de novo purine syntheses (Pietrzik et al., 2010; Scaglione and Panzavolta, 2014). To summarize, MTHFD1 and IMPDH2 are involved in the transcriptional synthesis of RNA and DNA. Indeed, IMPDH2 is the target of Mycophenolate mofetil (MMF) (Dalal et al., 2009), and the MTHFD1, similar to methotrexate (MTX), is associated with the metabolism of THF (de Jonge et al., 2005). It is not surprising that all the side effect, such as sperm reduction and leukopenia of MMF and MTX, overlap KX’s because they shared the same mechanism. Here we consider they are instead side effects than disease treating targets because they did not enrich in main therapeutic mechanism pathways in KX.

FIGURE 4
www.frontiersin.org

Figure 4 Disease network diagram. The red diamond represents some of the side effects caused by Tripterygium hypoglaucum hutch, and the green dots represent the disease-related genes. The direct linkage between genes was shown in PPI data of Biogrid.

We then analyze the rest herbs (Epimedium brevicornu Maxim, Cuscuta chinensis Lam, and Lycium barbarum L) of KX, which pointed out that they may help in reducing the side effect of H.Lév. Hutch. The analysis of the GO enrichment of rest herbs shows that the primary function involved in animal organ development, positive regulation of the cellular biosynthetic process, positive regulation of protein metabolic process, positive regulation of gene expression, etc. (Table 1). The KEGG enrichment analysis of the rest herbs shows that they have a regulatory effect on amino acid metabolism pathways (Table 2) and amino acids can be transformed into nucleic acid via specific metabolic pathways. To sum up, the rest herbs produce a positive regulation of the nucleic acid and amino acid biosynthetic process, positive regulation of transcription effects, resulting in consequences of alleviating adverse reactions caused by H.Lév. Hutch.

TABLE 1
www.frontiersin.org

Table 1 GO analysis of the rest of 3 herbs.

TABLE 2
www.frontiersin.org

Table 2 KEGG analysis of the rest of three herbs.

Validation From Published Biological Evaluation

KX Induces Apoptosis Through PI3K/AKT/mTOR Signaling Pathway

The PI3K/AKT/mTOR pathway has a pivotal role in mediating many cellular responses, including protection from apoptosis. H.Lév. Hutch is one of the ingredients of KX, which belongs to the Tripterygium wilfordii genus. Active ingredients include Tripterine (TP), Celastrol (Cel), wilforgine, etc. We made a review of the Tripterygium wilfordii genus herbs or their active ingredients to validate our hypothesis. The KX capsule or its active ingredients could promote varieties of cell apoptosis, including T cells, macrophages, dendritic cells, and fibroblast-like synoviocytes (FLS). Experimental evidences of KX capsule or its active ingredients for apoptosis were summarized in Table 3.

TABLE 3
www.frontiersin.org

Table 3 Experimental evidence of KX capsule or its active ingredients for apoptosis.

The KX Capsule or Its Ingredients Inhibit the Inflammatory Network of Th17, Macrophage, and FLS

The KX capsule or its ingredients could target at NF-κB, AP-1, IKK, HIF-1α, resulting in varieties of phenotype changes, such as inhibition of inflammatory cytokines released by macrophages, differentiation of T cell, invasion of FLS. Experimental evidence of the KX capsule or its active ingredients for inhibiting inflammatory network were summarized in Table 4.

TABLE 4
www.frontiersin.org

Table 4 KX capsule or its active ingredients for inhibiting inflammatory network.

The Relationship Between Amino Acid Metabolism and Toxicity

Nucleotide synthesis, a biological process that needs the exitance of cofactor tetrahydrofolate, is required in cell survival. Methotrexate (MTX), the main treatment in RA, inhibits the enzyme dihydrofolate reductase and then lead to tetrahydrofolate exhausted. The depletion of tetrahydrofolate causes cell death by suppressing DNA and RNA production (Kanarek et al., 2018). Toxicants, such as deoxynivalenol, could cause amino acids and nucleotides deficiency through repressing a variety of biosynthesis of amino acids and nucleotides, finally lead to hepatotoxicity, immunotoxicity and other toxicity (Ji et al., 2016; Geng et al., 2019; Tang et al., 2019). Furthermore, disruption of the TCA cycle and glutamate metabolism plays an essential role in the pathogenesis of amyotrophic lateral sclerosis (ALS), a disease that is closely related to the death of the neurocytes (Veyrat-Durebex et al., 2016). For cancer cell proliferation, these cells up-regulate serine synthesis, one-carbon (folate) metabolism, and the glycine cleavage system (SOG pathway) to fulfill the biosynthetic requirements of purines, ATP and NADPH. By contrast, MTX could inhibit the proliferation of cancer cells via disrupting the pathways (Tedeschi et al., 2013). TP (one of the active ingredients of H.Lév. Hutch) increases the level of cysteinylglycine and cysteine conjugates derived from NADPH-independent metabolism and promotes hepatotoxicity, which could be ameliorated by the supplement of glutathione (Du et al., 2014). Another active ingredient of H.Lév. Hutch, Cel suppresses tryptophan catabolism, shows its cytotoxic effect against colon cancer cells (Qi et al., 2018).

In contrast, positive regulation of the amino acid metabolism or supplement of some amino acids may help in cell survival and detoxication. MTX induces an increase in intestinal permeability. After MTX treatment, the Caco-2 cells barrier function is decreased, which was assessed by transepithelial electrical resistance, but could be reversed by the supplement of glutamine, glutamate, arginine, and leucine (Beutheu et al., 2013). Similarly, MTX-induced gastrointestinal toxicity in rats could be alleviated by exogenous glutamine and arginine (Gulgun et al., 2010) and supplement of branched-chain amino acids are beneficial effects on neuronal survival and axon regeneration (Park et al., 2017). Epimedium brevicornu Maxim, Cuscuta chinensis Lam, and Lycium barbarum L play a pivotal role in amino acid or nucleotide metabolism. Lycium barbarum L contains a variety of amino acids. Up to now, 18 amino acids have been identified, including 8 essential amino acids. This is the only plant that has been reported to contain taurine (Luo et al., 2004; Tian and Wang, 2006; Hsieh et al., 2018). Supplementation with Lycium barbarum polysaccharides (LBP) stimulates the growth and metabolic activity via increasing amino acid and nucleotide metabolism, short-chain fatty acid-related metabolism (Wang et al., 2019b). Significant metabolic disorders were observed in rats with a prolonged injection of hydrocortisone. Epimedium rebalances the multiple metabolic pathways involved in amino acid and phospholipid metabolism, for example, up-regulate taurine, glycine, and β-glucose levels (Huang et al., 2013; Pan et al., 2016). By far, there is insufficient data to reveal the three herbs in ameliorating H.Lév. Hutch-related side effects via exact molecule mechanism or amino acid/nucleotide metabolism. Indeed, the three herbs alleviate toxicity in different ways. Total flavonoids of Epimedium protect the reproductive system against cyclophosphamide-induced toxicity by inhibiting oxidative stress (Yuan et al., 2014). By contrast, Cuscuta chinensis Lam detoxicates cyclophosphamide-induced toxicity through regulating cytokines GM-CSF and TNF-α (Raju et al., 2015). Herein, for endeavoring to illustrate the coexisting mechanism of causing and alleviating side effects, we supposed amino acid and nucleotide metabolism approach.

Discussion

KX can act on the calcium signaling pathway, cAMP signaling pathway, cGMP-PKG signaling pathway, and G protein-coupled receptors (GPCRs). All these pathways are involved in cell signal transmission and could have crosstalk with PI3K/AKT/mTOR signaling pathway, which is used to achieve the proliferation, differentiation, and apoptosis of immune cells in RA (Shu et al., 2017; Feng and Qiu, 2018; Wang et al., 2019a). In physiological states, regulation of PI3K by GPCRs mainly occurs through Gβγ, which directly binds to the p110 catalytic subunit and protein kinase A (PKA) (Wang et al., 2019a). PI3K phosphorylates AKT, which has many downstream cascades, including mTOR, a transcription factor encodes genes of cell proliferation and survival (Keppler-Noreuil et al., 2016). Second, JAK-STAT signaling can interconnect with the PI3K/AKT/mTOR pathway (Rawlings et al., 2004). Given that JAKs are associated with cytokine receptors, such as cytokines IFNγ, IL-2, IL-4, and IL-10 (Rodig et al., 1998), it may exist a signaling transmission from cytokines receptors to mTOR. Lastly, cAMP, a second messenger, could be affected by the calcium signaling pathway, could activate PI3K via intermediate molecule such as PKA. The molecules mentioned above were high degree targets in KX, once they were inhibited, resulting in the decrease of proliferation, differentiation, and apoptosis of immune cells in RA (Figure 5A).

FIGURE 5
www.frontiersin.org

Figure 5 Mechanism of Kunxian capsule in treating RA. The red nodes represent targets in Kunxian capsule. (A) Through intermediate protein Gβγ, GPCRs regulated the activation of PI3Ks. The phosphorylation of PI3K activated AKT which lead to downstream cascades, such as controlling cell proliferation, protein translation. JAK-STAT signaling is related with the proliferation of the inflammatory cells. Lastly, cAMP, a second messenger affected by calcium signaling pathway, could activate PI3K via intermediate molecules, such as PKA and mTOR, which are high degree targets in KX. KX inhibits the three pathways and decrease cell proliferation and promote cell apoptosis in RA. (B) Th17 cells could be activated by IL-1β (through mTOR-ROR) or IL-6 (through STAT3-HIF1-ROR). Subsequently, the activated TH17 cells released IL-17, which activated the downstream effector cells, such as monocytes, macrophages, and synovial fibroblast. These cells could release inflammatory cytokines which could promoted cartilage damage. KX could target at the critical checkpoints in the pathways, downregulating the inflammatory cells proliferation and alleviating cartilage damage.

At the same time, KX can alleviate the inflammatory response and reduce the disease activities by regulating multiple key inflammatory signaling pathways in RA, such as IL-17 signaling pathway, Th17 cell differentiation signaling pathway, and TNF signaling pathway (Lubberts, 2015). In physiological status, Th17 cells could be activated by IL-1β (through mTOR-ROR) or IL-6 (through STAT3-HIF1-ROR). Subsequently, IL-17 released by Th17 stimulated the downstream effector cells, such as monocytes, macrophages, and synovial fibroblasts. These cells could release inflammatory cytokines, which could promote cartilage damage. KX could target at the critical checkpoints in the pathways, downregulating the inflammatory cells proliferation and alleviating cartilage damage (Figure 5B).

In terms of adverse reactions, it was found that H.Lév. Hutch can act on IMPDH2 and MTHFD1.IMPDH2 encodes the rate-limiting enzyme (inosine 5′-monophosphate dehydrogenase) in the de novo guanine nucleotide biosynthesis. It is thus involved in maintaining cellular guanine deoxy- and ribonucleotide pools, which needed for DNA and RNA synthesis (Dalal et al., 2009). The functions of MTHFD1 have been described previously. Tetrahydrofolate is an essential coenzyme in the one-carbon transferase system. The coenzyme is produced by folic acid catalyzed by folate reductase to form dihydrofolate. Dihydrofolate is then catalyzed by dihydrofolate reductase to form tetrahydrofolate. Tetrahydrofolic acid is a carrier of the one-carbon group, which can participate in the purine’s synthesis, and promote the formation of healthy blood cells (Watkins et al., 2011).

IMPDH2 and MTHFD1 are involved in nucleic acid metabolism and are responsible for perm reduction and leukopenia. Targets of ASS1, GOT1, NOS2 in rest herbs form a positive regulation of amino acid. Subsequently, amino acid could be changed to purines via the specific metabolic pathway, which compensating purines metabolic disrupted by H.Lév. Hutch. Besides, the rest herbs take effect in the transition of purines to nucleic acids (Figure 6). Interestingly, they enriched in drug metabolism-cytochrome P450 signaling pathway, indicating that they may help in promoting drug metabolism.

FIGURE 6
www.frontiersin.org

Figure 6 The mechanism of side effect caused by Tripterygium hypoglaucum hutch and the mechanism of alleviation compensated by the rest of herbs. Tripterygium hypoglaucum hutch may inhibit the activity of IMPDH2 which encodes the rate-limiting enzyme (inosine 5′-monophosphate dehydrogenase) in the de novo guanine nucleotide biosynthesis. As a result, the metabolism of nucleic acid was inhibited because of the downregulation of guanine. Besides, Tripterygium hypoglaucum hutch may inhibit the activity of MTHFD1 which encodes a protein that possesses the folic acid metabolism related enzymes activities, disrupting the purines converted to nucleic acid. The rest herbs could act on ASS1, GOT1, NOS2 et al. in the amino acid metabolism pathways, including glycine, serine, and threonine and form a positive regulation of amino acid. After that amino acid could be changed to purines via specific metabolic pathway, which compensating purines metabolic disrupted by Tripterygium hypoglaucum hutch. Besides, the rest herbs could play roles in positive regulation of gene expression, positive regulation of transcription, DNA-templated, positive regulation of RNA biosynthetic process, helping in the transition of purines to nucleic acid.

Zhang (Zhang et al., 2014) has validated that H.Lév. Hutch, one herb in KX, could decrease the IL-17 level in serum in the CIA rats, but the exact mechanism remains unknown. Although some targets have been validated by different researchers using active ingredients, there may be synergism, antagonism, potentiation to show us a different result. Some of the targets with high degrees on these two pathways, such as mTOR, NF-κB, can be experimentally verified in the future. We speculate that KX can inhibit the targets mentioned above, decreasing TH17 cell differentiation and reducing IL-17 levels. It also needs to verify that whether KX could disrupt the PI3K-AKT-mTOR pathway and, thus, promote the apoptosis of related inflammatory cells. Interestingly, KEGG enrichment analysis showed that KX also has neuroactive ligand-receptor interaction, which may have central analgesic effects.

The mechanism prediction of KX via bioinformatics is limited because the minority components of the herbs are analyzed in our study, which may make little sense in clinical practice. However, we believe that the systematic approach to infer therapeutic mechanisms of KX for RA treatment provides new insights into the understanding of this TCM formula and provides a new perspective in TCM patent prescription studies. What is more, we provide new insights into the understanding of drug matching of TCM to alleviate side effects.

Conclusion

KX treats RA via inhibiting the inflammatory network and inducing apoptosis through PI3K/AKT/mTOR signaling pathway. H.Lév. Hutch, one of the ingredients of KX, could cause some side effects, such as sperm reduction and leukopenia. The rest 3 herbs (Epimedium brevicornu Maxim, Cuscuta chinensis Lam, and Lycium barbarum L) could play a role of detoxication through positive regulation of the amino acid metabolism or supplement of some amino acids. In conclusion, KX is a relatively safe alternative approach in RA treatment.

Materials and Methods

Target Prediction of Kunxian Capsule

There are four main components of KX, H.Lév. Hutch (Kun Ming Shan Hai Tang), Epimedium brevicornu Maxim (Yin Yang Huo), Cuscuta chinensis Lam (Tu Si Zi), and Lycium barbarum L (Gou Qi Zi). The component and targets of Epimedium brevicornu Maxim (Yin Yang Huo), Cuscuta chinensis Lam (Tu Si Zi), and Lycium barbarum L (Gou Qi Zi) were searched at BATMAN-TCM. The Cut-off value is set to 20. Since the BATMAN-TCM (http://bionet.ncpsb.org/batman-tcm/) database does not record H.Lév. Hutch component data, a literature search is required. Finally, we used the data of a published master’s thesis in Chinese (Li, 2019). The thesis summarized the components of H.Lév. Hutch, which has been identified by now. Then the compound was sorted and converted to PubChem_CID in PubChem (https://pubchem.ncbi.nlm.nih.gov/).

Screening of Related Genes and Targets

CTD is a publicly available database that contains robust disease genes. “Rheumatoid arthritis” was entered into the CTD database to obtain any possible disease target genes.

OMIM is a comprehensive, authoritative compendium of human genes and genetic phenotypes database, which can be used to establish disease networks via pathological disorder. The OMIM database was retrieved by PubMed, and the retrieved relevant disease genes were recorded.

Visualization and Functional Analysis

Drug-Target-Disease Interaction Network and Adverse Reaction Disease Network Construction

Targets of the Kunxian capsule and RA-related genes were input into STRING (https://string-db.org/) to attain protein interaction networks, with a minimum required interaction score 0.9. Cytoscape 3.7.0 (https://cytoscape.org/) software was applied for PPI visualization. The adverse reaction disease network was constructed using the gene retrieved in OMIM with a minimum required interaction score of 0.4.

Screening of Hub Genes

CytoHubba, a novel Cytoscape plugin that can be used to rank node features in the network. CytoHubba provides 11 topological analysis methods, including Degree, Edge Percolated Component, Maximum Neighborhood Component, Density of Maximum Neighborhood Component, Maximal Clique Centrality, and six shortest path centralities (Bottleneck, EcCentricity, Closeness, Radiality, Betweenness, Stress). Among the 11 methods, the MCC PPI network has better performance in predicting the accuracy of key targets (Shu et al., 2017). The MCC algorithm is used to construct the PPI network and select the top 50 key targets.

GO Enrichment and KEGG Pathway Analysis

Bioconductor’s “clusterProfiler,” “pathview,” and other functional packages (Hueber et al., 2010; Li et al., 2012; Jiao et al., 2018) in the R software were utilized for GO enrichment, KEGG pathway analysis, and plotting. The cutoff values and P value are both set to 0.05, and the top 20 results are selected for drawing.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author Contributions

YT: data analyses, figure preparation, manuscript preparation, and study initiation. YZ: interpretation, data analyses, and manuscript preparation/submission. LL: data analyses, and manuscript preparation. ZX: data interpretation, figure preparation. CW: critically reviewed the manuscript, project funding. LH: critically reviewed the manuscript and study initiation. All authors approved the final version of the manuscript.

Funding

Research was funded by National Key R&D Program of China (2018YFC1705500).

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.

Acknowledgments

We would like to thank Yujie Tang who majors in the arts. She helped us to arrange elements position and colors of Figures 5 and 6.

Supplementary Material

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

Abbreviations

RA, Rheumatoid arthritis; KX, KunXian capsule; (H.Lév.) Hutch, Tripterygium hypoglaucum Hutch; TCM, Traditional Chinese medicine; BATMAN-TCM, Bioinformatics analysis tool for molecular mechanism of TCM; AGT, Angiotensinogen; GNG2, G protein subunit gamma 2; ANXA1, Annexin A1; GO, Gene ontology; TNF, Tumor necrosis factor; PI3K, Phosphoinositide 3-kinase; AKT (also known as PKB), Protein kinase B; PPI, Protein-protein interaction; IMPDH1, Inosine monophosphate dehydrogenase 1; MTHFD1, Methylenetetrahydrofolate dehydrogenase, cyclohydrolase, and formyltetrahydrofolate synthetase 1; PKD2, Polycystin 2; PTGIS, Prostaglandin I2 synthase; SLC25A13, Solute carrier family 25 member 13; AGTR1, Angiotensin II receptor type 1; NR0B1, Nuclear receptor subfamily 0 group B member 1; AR, Androgen receptor; THF, Tetrahydrofolate; MMF, Mycophenolate mofetil; MTX, Methotrexate; mTOR, Mechanistic target of rapamycin kinase; TP, Tripterine; Cel, Celastrol; FLS, Fibroblast-like synoviocytes; NF-κB, Nuclear factor kappa B; AP-1, Activator protein 1; IKK, IκB kinase; HIF-1α, Hypoxia-inducible factor-1α; TCA cycle, Tricarboxylic acid cycle; ALS, Amyotrophic lateral sclerosis; LBP, Lycium barbarum polysaccharides; GM-CSF, Granulocyte-macrophage colony-stimulating factor; GPCRs, G protein-coupled receptors; PKA, Protein kinase A; JAK, Janus kinase; STAT, Signal transducer and activator of transcription; IFNγ, Interferon gamma; IL, Interleukin; Th, T helper cell; ASS1, Argininosuccinate synthase 1; GOT1, Glutamic-oxaloacetic transaminase 1; NOS2, Nitric oxide synthase 2; CIA, Collagen-induced arthritis.

References

Alamanos, Y., Voulgari, P. V., Drosos, A. A. (2006). Incidence and Prevalence of Rheumatoid Arthritis, Based on the 1987 American College of Rheumatology Criteria: A Systematic Review. Semin. Arthritis Rheumatism 36 (3), 182–188. doi: 10.1016/j.semarthrit.2006.08.006

CrossRef Full Text | Google Scholar

Bao, X., Cui, J., Wu, Y., Han, X., Gao, C., Hua, Z., et al. (2006). The roles of endogenous reactive oxygen species and nitric oxide in triptolide-induced apoptotic cell death in macrophages. J. Mol. Med. 85 (1), 85. doi: 10.1007/s00109-006-0113-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Beutheu, S., Ghouzali, I., Galas, L., Dechelotte, P., Coeffier, M. (2013). Glutamine and arginine improve permeability and tight junction protein expression in methotrexate-treated Caco-2 cells. Clin. Nutr. 32 (5), 863–869. doi: 10.1016/j.clnu.2013.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, Y., Yun, N., Zhou, A. (2018). Meta-analysis of ADR Induced by Tripterysium Glycosides Tablet. China Pharm. 29 (01), 125–130. doi: 10.6039/j.issn.1001-0408.2018.01.31

CrossRef Full Text | Google Scholar

Cascão, R., Vidal, B., Raquel, H., Neves-Costa, A., Figueiredo, N., Gupta, V., et al. (2012). Effective treatment of rat adjuvant-induced arthritis by celastrol. Autoimmun. Rev. 11 (12), 856–862. doi: 10.1016/j.autrev.2012.02.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Dalal, P., Grafals, M., Chhabra, D., Gallon, L. (2009). Mycophenolate mofetil: safety and efficacy in the prophylaxis of acute kidney transplantation rejection. Ther. Clin. Risk Manag 5 (1), 139–149. doi: 10.2147/tcrm.s3068

PubMed Abstract | CrossRef Full Text | Google Scholar

de Jonge, R., Hooijberg, J. H., van Zelst, B. D., Jansen, G., van Zantwijk, C. H., Kaspers, G. J., et al. (2005). Effect of polymorphisms in folate-related genes on in vitro methotrexate sensitivity in pediatric acute lymphoblastic leukemia. Blood 106 (2), 717–720. doi: 10.1182/blood-2004-12-4941

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, F., Liu, Z., Li, X., Xing, J. (2014). Metabolic pathways leading to detoxification of triptolide, a major active component of the herbal medicine Tripterygium wilfordii. J. Appl. Toxicol. 34 (8), 878–884. doi: 10.1002/jat.2906

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, F. B., Qiu, H. Y. (2018). Effects of Artesunate on chondrocyte proliferation, apoptosis and autophagy through the PI3K/AKT/mTOR signaling pathway in rat models with rheumatoid arthritis. BioMed. Pharmacother. 102, 1209–1220. doi: 10.1016/j.biopha.2018.03.142

PubMed Abstract | CrossRef Full Text | Google Scholar

Geng, N., Ren, X., Gong, Y., Zhang, H., Wang, F., Xing, L., et al. (2019). Integration of metabolomics and transcriptomics reveals short-chain chlorinated paraffin-induced hepatotoxicity in male Sprague-Dawley rat. Environ. Int. 133, 105231. doi: 10.1016/j.envint.2019.105231

PubMed Abstract | CrossRef Full Text | Google Scholar

Gulgun, M., Karaoglu, A., Kesik, V., Kurt, B., Erdem, O., Tok, D., et al. (2010). Effect of proanthocyanidin, arginine and glutamine supplementation on methotrexate-induced gastrointestinal toxicity in rats. Methods Find Exp. Clin. Pharmacol. 32 (9), 657–661. doi: 10.1358/mf.2010.32.9.1516694

PubMed Abstract | CrossRef Full Text | Google Scholar

Ho, L. J., Chang, W. L., Chen, A., Chao, P., Lai, J. H. (2013). Differential immunomodulatory effects by Tripterygium wilfordii Hook f-derived refined extract PG27 and its purified component PG490 (triptolide) in human peripheral blood T cells: potential therapeutics for arthritis and possible mechanisms explaining in part Chinese herbal theory “Junn-Chenn-Zuou-SS”. J. Transl. Med. 11, 294. doi: 10.1186/1479-5876-11-294

PubMed Abstract | CrossRef Full Text | Google Scholar

Hopkins, A. L. (2008). Network pharmacology: the next paradigm in drug discovery. Nat. Chem. Biol. 4 (11), 682–690. doi: 10.1038/nchembio.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsieh, L. Y., Chan, H. H., Kuo, P. C., Hung, H. Y., Li, Y. C., Kuo, C. L., et al. (2018). A feasible and practical (1)H NMR analytical method for the quality control and quantification of bioactive principles in Lycii Fructus. J. Food Drug Anal. 26 (3), 1105–1112. doi: 10.1016/j.jfda.2018.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D., Yang, J., Lu, X., Deng, Y., Xiong, Z., Li, F. (2013). An integrated plasma and urinary metabonomic study using UHPLC–MS: Intervention effects of Epimedium koreanum on ‘Kidney-Yang Deficiency syndrome’ rats. J. Pharmaceutical Biomed. Anal. 76, 200–206. doi: 10.1016/j.jpba.2012.12.022

CrossRef Full Text | Google Scholar

Hueber, W., Patel, D. D., Dryja, T., Wright, A. M., Koroleva, I., Bruin, G., et al. (2010). Effects of AIN457, a fully human antibody to interleukin-17A, on psoriasis, rheumatoid arthritis, and uveitis. Sci. Transl. Med. 2 (52), 52ra72. doi: 10.1126/scitranslmed.3001107

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, J., Sun, J., Pi, F., Zhang, S., Sun, C., Wang, X., et al. (2016). GC-TOF/MS-based metabolomics approach to study the cellular immunotoxicity of deoxynivalenol on murine macrophage ANA-1 cells. Chemico-Biological Interact. 256, 94–101. doi: 10.1016/j.cbi.2016.06.017

CrossRef Full Text | Google Scholar

Jia, X., Liu, P. (2016). Research progress of Kunxian capsule in the treatment of rheumatoid arthritis. World Latest Med. Inf. 25–26. doi: 10.3969/j.issn.1671-3141.2016.16.015

CrossRef Full Text | Google Scholar

Jiao, Y., Ding, H., Huang, S., Liu, Y., Sun, X., Wei, W., et al. (2018). Bcl-XL and Mcl-1 upregulation by calreticulin promotes apoptosis resistance of fibroblast-like synoviocytes via activation of PI3K/Akt and STAT3 pathways in rheumatoid arthritis. Clin. Exp. Rheumatol 36 (5), 841–849.

PubMed Abstract | Google Scholar

Kanarek, N., Keys, H. R., Cantor, J. R., Lewis, C. A., Chan, S. H., Kunchok, T., et al. (2018). Histidine catabolism is a major determinant of methotrexate sensitivity. Nature 559 (7715), 632–636. doi: 10.1038/s41586-018-0316-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Keppler-Noreuil, K. M., Parker, V. E., Darling, T. N., Martinez-Agosto, J. A. (2016). Somatic overgrowth disorders of the PI3K/AKT/mTOR pathway & therapeutic strategies. Am. J. Med. Genet. C Semin. Med. Genet. 172 (4), 402–421. doi: 10.1002/ajmg.c.31531

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, Q., Ma, Y., Yu, J., Chen, X. (2017). Predicted molecular targets and pathways for germacrone, curdione, and furanodiene in the treatment of breast cancer using a bioinformatics approach. Sci. Rep. 7 (1), 15543. doi: 10.1038/s41598-017-15812-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, P. P., Liu, D. D., Liu, Y. J., Song, S. S., Wang, Q. T., Chang, Y., et al. (2012). BAFF/BAFF-R involved in antibodies production of rats with collagen-induced arthritis via PI3K-Akt-mTOR signaling and the regulation of paeoniflorin. J. Ethnopharmacol 141 (1), 290–300. doi: 10.1016/j.jep.2012.02.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G. Q., Liu, D., Zhang, Y., Qian, Y. Y., Zhu, Y. D., Guo, S. Y., et al. (2013). Anti-invasive effects of celastrol in hypoxia-induced fibroblast-like synoviocyte through suppressing of HIF-1alpha/CXCR4 signaling pathway. Int. Immunopharmacol 17 (4), 1028–1036. doi: 10.1016/j.intimp.2013.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. (2019). Chemical components of Tripterygium hypoglaucum (Lévl.) Hutchins. from Guizhou Province. dissertation (Guiyang: Guizhou University).

Google Scholar

Liu, Z., Guo, F., Wang, Y., Li, C., Zhang, X., Li, H., et al. (2016). BATMAN-TCM: a Bioinformatics Analysis Tool for Molecular mechANism of Traditional Chinese Medicine. Sci. Rep. 6 (1), 21146. doi: 10.1038/srep21146

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Li, N., Sun, J., Han, Y., Wu, C., Lin, C. (2017). Analysis of Adverse Reaction of Common Chinese Herbs with Antirheumatic Activity. Liaoning J. Tradit. Chin. Med. 44 (10), 2138–2140. doi: 10.13192/j.issn.1000-1719.2017.10.041

CrossRef Full Text | Google Scholar

Liu, Q. (2004). The Effects of triptolide on Regulating of Dendritic Cel Functions and Inhibiting Tumor Cell Growth. dissertation (Wuhan: Hubei Medical University).

Google Scholar

Lubberts, E. (2015). The IL-23-IL-17 axis in inflammatory arthritis. Nat. Rev. Rheumatol 11 (7), 415–429. doi: 10.1038/nrrheum.2015.53

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, Q., Cai, Y., Yan, J., Sun, M., Corke, H. (2004). Hypoglycemic and hypolipidemic effects and antioxidant activity of fruit extracts from Lycium barbarum. Life Sci. 76 (2), 137–149. doi: 10.1016/j.lfs.2004.04.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, S., Chen, A., Han, Z., Wang, Y., Lu, X., Yang, Y. (2016). 1H NMR-based metabonomic study on the effects of Epimedium on glucocorticoid-induced osteoporosis. J. Chromatography B 1038, 118–126. doi: 10.1016/j.jchromb.2016.10.015

CrossRef Full Text | Google Scholar

Park, K. K., Luo, X., Mooney, S. J., Yungher, B. J., Belin, S., Wang, C., et al. (2017). Retinal ganglion cell survival and axon regeneration after optic nerve injury in naked mole-rats. J. Comp. Neurol. 525 (2), 380–388. doi: 10.1002/cne.24070

PubMed Abstract | CrossRef Full Text | Google Scholar

Pietrzik, K., Bailey, L., Shane, B. (2010). Folic acid and L-5-methyltetrahydrofolate: comparison of clinical pharmacokinetics and pharmacodynamics. Clin. Pharmacokinet 49 (8), 535–548. doi: 10.2165/11532990-000000000-00000

PubMed Abstract | CrossRef Full Text | Google Scholar

Ping, Q., Zhou, Y., Zhang, S., Cao, J., Xu, L., Fang, G., et al. (2015). Study on effects of Tripterygium wilfordii polycoride in resisting macrophage inflammation and regulating inflammation via TLR4/NF-κB. China J. Chin. materia Med. 40 (16), 3256–3261.

Google Scholar

Qi, Y., Wang, R., Zhao, L., Lv, L., Zhou, F., Zhang, T., et al. (2018). Celastrol Suppresses Tryptophan Catabolism in Human Colon Cancer Cells as Revealed by Metabolic Profiling and Targeted Metabolite Analysis. Biol. Pharmaceutical Bull. 41 (8), 1243–1250. doi: 10.1248/bpb.b18-00171

CrossRef Full Text | Google Scholar

Raju, N., Sakthivel, K. M., Kannan, N., Vinod Prabhu, V., Guruvayoorappan, C. (2015). Cuscuta chinensis Ameliorates Immunosuppression and Urotoxic Effect of Cyclophosphamide by Regulating Cytokines - GM-CSF and TNF-Alpha. Appl. Biochem. Biotechnol. 176 (3), 742–757. doi: 10.1007/s12010-015-1608-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rawlings, J. S., Rosler, K. M., Harrison, D. A. (2004). The JAK/STAT signaling pathway. J. Cell Sci. 117 (8). 1281 doi: 10.1242/jcs.00963

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodig, S. J., Meraz, M. A., White, J. M., Lampe, P. A., Riley, J. K., Arthur, C. D., et al. (1998). Disruption of the Jak1 gene demonstrates obligatory and nonredundant roles of the Jaks in cytokine-induced biologic responses. Cell 93 (3), 373–383. doi: 10.1016/s0092-8674(00)81166-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Scaglione, F., Panzavolta, G. (2014). Folate, folic acid and 5-methyltetrahydrofolate are not the same thing. Xenobiotica 44 (5), 480–488. doi: 10.3109/00498254.2013.845705

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, J., Zhang, F., Zhang, L., Wei, W. (2017). G protein coupled receptors signaling pathways implicate in inflammatory and immune response of rheumatoid arthritis. Inflammation Res. 66 (5), 379–387. doi: 10.1007/s00011-016-1011-5

CrossRef Full Text | Google Scholar

Silman, A. J., Pearson, J. E. (2002). Epidemiology and genetics of rheumatoid arthritis. Arthritis Res. Ther. 4 (3), S265. doi: 10.1186/ar578

CrossRef Full Text | Google Scholar

Su, J., Cheng, J., Sun, H. X., Diao, Z. Y., Zhen, X., Yang, J., et al. (2014). Tripterygium glycosides impairs the proliferation of granulosa cells and decreases the reproductive outcomes in female rats. Birth Defects Res. B Dev. Reprod. Toxicol. 101 (3), 283–291. doi: 10.1002/bdrb.21111

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Z., Li, Y., Jiang, Y., Cheng, J., Xu, S., Zhang, J. (2019). Cellular metabolomics reveals glutamate and pyrimidine metabolism pathway alterations induced by BDE-47 in human neuroblastoma SK-N-SH cells. Ecotoxicol. Environ. Saf. 182, 109427. doi: 10.1016/j.ecoenv.2019.109427

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedeschi, P. M., Markert, E. K., Gounder, M., Lin, H., Dvorzhinski, D., Dolfi, S. C., et al. (2013). Contribution of serine, folate and glycine metabolism to the ATP, NADPH and purine requirements of cancer cells. Cell Death Dis. 4 (10), e877–e877. doi: 10.1038/cddis.2013.393

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, M., Wang, M. (2006). Studies on extraction, isolation and composition of Lycium barbarum polysaccharides. China J. Chin. Materia Med. 31 (19), 1603–1607.

Google Scholar

Tong, Z., Cheng, L., Song, J., Wang, M., Yuan, J., Li, X., et al. (2018). Therapeutic effects of Caesalpinia minax Hance on complete Freund’s adjuvant (CFA)-induced arthritis and the anti-inflammatory activity of cassane diterpenes as main active components. J. Ethnopharmacol. 226, 90–96. doi: 10.1016/j.jep.2018.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Veyrat-Durebex, C., Corcia, P., Piver, E., Devos, D., Dangoumau, A., Gouel, F., et al. (2016). Disruption of TCA Cycle and Glutamate Metabolism Identified by Metabolomics in an In Vitro Model of Amyotrophic Lateral Sclerosis. Mol. Neurobiol. 53 (10), 6910–6924. doi: 10.1007/s12035-015-9567-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Zuo, S., Liu, Z., Ji, X., Yao, Z., Wang, X. (2018). Study on the efficacy and mechanism of triptolide on treating TNF transgenic mice with rheumatoid arthritis. BioMed. Pharmacother. 106, 813–820. doi: 10.1016/j.biopha.2018.07.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D., Zhou, W., Chen, J., Wei, W. (2019a). Upstream regulators of phosphoinositide 3-kinase and their role in diseases. J. Cell Physiol. 234, 14460–14472. doi: 10.1002/jcp.28215

CrossRef Full Text | Google Scholar

Wang, M., Xie, Z., Li, L., Chen, Y., Li, Y., Wang, Y., et al. (2019b). Supplementation with compound polysaccharides contributes to the development and metabolic activity of young rat intestinal microbiota. Food Funct. 10 (5), 2658–2675. doi: 10.1039/C8FO02565G

PubMed Abstract | CrossRef Full Text | Google Scholar

Watkins, D., Schwartzentruber, J. A., Ganesh, J., Orange, J. S., Kaplan, B. S., Nunez, L. D., et al. (2011). Novel inborn error of folate metabolism: identification by exome capture and sequencing of mutations in the MTHFD1 gene in a single proband. J. Med. Genet. 48 (9), 590–592. doi: 10.1136/jmedgenet-2011-100286

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, V. K. W., Qiu, C., Xu, S. W., Law, B. Y. K., Zeng, W., Wang, H., et al. (2019). Ca(2+) signaling plays a role in celastrol-mediated suppression of synovial fibroblasts of rheumatoid arthritis patients and experimental arthritis in rats. Br. J. Pharmacol. 176 (16), 2922–2944. doi: 10.1111/bph.14718

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, D., Wang, H., He, H., Jia, L., He, Y., Wang, T., et al. (2014). Protective Effects of Total Flavonoids from Epimedium on the Male Mouse Reproductive System Against Cyclophosphamide-Induced Oxidative Injury by Up-Regulating the Expressions of SOD3 and GPX1. Phytotherapy Res. 28 (1), 88–97. doi: 10.1002/ptr.4956

CrossRef Full Text | Google Scholar

Zhang, F., Zhou, H., Cui, D., Yan, Y., Dong, C. (2014). Effects of THH on IL- 6, IL- 17 and IFN-γ in Collagen Induced Arthritis Rat. Chin. Foreign Med. Res. 12 (13), 138–139. doi: 10.14033/j.cnki.cfmr.2014.13.006

CrossRef Full Text | Google Scholar

Zhou, J., Xiao, W., Wu, R. (2016). Effectiveness and Safety of Kunxian Capsule for Rueumatoid Arthritis:A Systematic Review. J. Liaoning Univ. Traditional Chin. Med. 18 (10), 122–126. doi: 10.13194/j.issn.1673-842x.2016.10.037

CrossRef Full Text | Google Scholar

Keywords: Kunxian Capsule, rheumatoid arthritis, pharmacological mechanism, side effects alleviation, drug matching, bioinformatics

Citation: Tang Y, Zhang Y, Li L, Xie Z, Wen C and Huang L (2020) Kunxian Capsule for Rheumatoid Arthritis: Inhibition of Inflammatory Network and Reducing Adverse Reactions Through Drug Matching. Front. Pharmacol. 11:485. doi: 10.3389/fphar.2020.00485

Received: 03 June 2019; Accepted: 27 March 2020;
Published: 17 April 2020.

Edited by:

Claudio Ferrante, Università degli Studi G. d’Annunzio Chieti e Pescara, Italy

Reviewed by:

Lorenzo Di Cesare Mannelli, University of Florence, Italy
Sengul Uysal, Selçuk University, Turkey

Copyright © 2020 Tang, Zhang, Li, Xie, Wen and Huang. 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: Chengping Wen, wengcp@163.com; Lin Huang, huanglin@zcmu.edu.cn

These authors have contributed equally to this work

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.