Evidence-Based Complementary and Alternative Medicine

Evidence-Based Complementary and Alternative Medicine / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 8406127 | https://doi.org/10.1155/2021/8406127

Zhen-Zhen Wang, Yuan Jia, Kamal D. Srivastava, Weihua Huang, Raj Tiwari, Anna Nowak-Wegrzyn, Jan Geliebter, Mingsan Miao, Xiu-Min Li, "Systems Pharmacology and In Silico Docking Analysis Uncover Association of CA2, PPARG, RXRA, and VDR with the Mechanisms Underlying the Shi Zhen Tea Formula Effect on Eczema", Evidence-Based Complementary and Alternative Medicine, vol. 2021, Article ID 8406127, 17 pages, 2021. https://doi.org/10.1155/2021/8406127

Systems Pharmacology and In Silico Docking Analysis Uncover Association of CA2, PPARG, RXRA, and VDR with the Mechanisms Underlying the Shi Zhen Tea Formula Effect on Eczema

Academic Editor: Nasiara Karim
Received02 Jul 2020
Revised18 Mar 2021
Accepted27 Mar 2021
Published17 May 2021

Abstract

Eczema is a complex chronic inflammatory skin disease impacted by environmental factors, infections, immune disorders, and deficiencies in skin barrier function. Shi Zhen Tea (SZT), derived from traditional Chinese medicine Xiao-Feng-San, has shown to be an effective integrative therapy for treating skin lesions, itching, and sleeping loss, and it facilitates reduction of topical steroid and antihistamine use in pediatric and adult patients with severe eczema. Yet, its active compounds and therapeutic mechanisms have not been elucidated. In this study, we sought to investigate the active compounds and molecular mechanisms of SZT in treating eczema using systems pharmacology and in silico docking analysis. SZT is composed of 4 medicinal herbs, Baizhu (Atractylodis macrocephalae rhizome), Jingjie (Schizonepetae herba), Kushen (Sophorae flavescentis radix), and Niubangzi (Arctii fructus). We first identified 51 active compounds from SZT and their 81 potential molecular targets by high-throughput computational analysis, from which we identified 4 major pathways including Th17 cell differentiation, metabolic pathways, pathways in cancer, and the PI3K-Akt signaling pathway. Through network analysis of the compound-target pathway, we identified hub molecular targets within these pathways including carbonic anhydrase II (CA2), peroxisome proliferator activated receptor γ (PPAR γ), retinoid X receptor α (RXRA), and vitamin D receptor (VDR). We further identified top 5 compounds including cynarine, stigmasterin, kushenol, β-sitosterol, and (24S)-24-propylcholesta-5-ene-3β-ol as putative key active compounds on the basis of their molecular docking scores with identified hub target proteins. Our study provides an insight into the therapeutic mechanism underlying multiscale benefits of SZT for eczema and paves the way for developing new and potentially more effective eczema therapies.

1. Introduction

Eczema, characterized by itchy, scaly, erythematous, and oozing skin, is the most common chronic inflammatory skin disorder, affecting 10–20% of children and 1–3% of adults [1, 2]. In addition, allergic diseases such as allergic rhinitis, food allergy, and asthma are frequent comorbidities in moderate-to-severe eczema patients. Eczema severely impacts the life quality of patients and exerts tremendous financial and social burden [3]. In the United States, direct economic costs for the treatment of atopic eczema reached approximately $3.8 billion per year in 2015 [4]. Eczema was ranked first in 2010 global “skin disease burden,” which takes into account severe health loss, including psychological, social, and financial consequences [5]. The pathogenesis of eczema is complicated and is associated with genetic and environmental factors [6]. Severe eczema often fails to respond adequately to therapies specifically targeting a single cytokine or inflammatory mediator. To date, there is no effective cure for eczema, although some treatments can temporarily relieve skin symptoms and reduce inflammation [7]. Thus, the development of new treatments to effectively regulate the immune system and rebuild the skin barrier is an unmet need.

Traditional Chinese Medicine (TCM) has demonstrated significant success in treating complex diseases by regulating multiple pathways working to return dysfunctional organ systems back to normal [8]. Over the past 20 years, multiple studies have reported the efficacy of TCM for the treatment of eczema without serious adverse effects [9, 10]. For example, composite poria granules have shown clinical efficacy in eczema treatment [11]. Recently, Run-Zao-Zhi-Yang capsules (RZZYC) have been reported to provide symptom relief, maintain long-term remission, and improve quality of life, with fewer recurrences [12]. In addition, the well-known famous TCM formula Xiao-Feng-San (XFS) has been used for the treatment of skin disorders since the ancient Ming Dynasty in China [13]. It was reported to be one of the top 5 most commonly prescribed herbal formulas for eczema (total prescriptions = 381,282) in Taiwan [14]. We developed Shi Zhen Tea (SZT), a derivative of XFS that has been successfully used in the US as integrative eczema therapy [8]. It is comprised of four medicinal herbs, Sophorae flavescentis radix (Kushen, K), Schizonepetae herba (Jingjie, J), Arctii fructus (Niubangzi, N), and Atractylodis macrocephalae rhizoma (Baizhu, B), which display strong ability of dispersing cold, dispelling wind, and dehumidification. It showed efficacy in treating moderate-to-severe eczema in both children and adults [15]. Disease severity scores decreased in all patients, and most patients stopped taking steroids or antihistamines after the 3-month treatment. However, as with most traditional medicines, its mechanisms of action are elusive, due to presence of multiple, complex compounds, and their metabolites. By combining discovery of new actionable disease targets and isolation of active compounds, the systems pharmacology approach together with high-throughput computational analysis provides a powerful tool to study mechanisms underlying TCM function. A growing body of evidence [1618] suggests that the use of systems pharmacology can lead to an integrated understanding of drug action on the human interactome which can be used for drug discovery and pharmacological analysis of multiple agents on complex diseases along with prediction of adverse effects. Recently, certain classical Chinese prescriptions, such as Huo-Xiang-Zheng-Qi formula [19], Dan-Shen formula [20], Xiao-Zheng-Fang formula [21], and Bu-Fei-Yi-Shen formula [22], as well as antiasthma herbal medicine intervention (ASHMI), a formula previously developed by us [23], have been analyzed and evaluated by systems pharmacology.

The current study seeks to uncover the therapeutic mechanisms’ underlying efficacy of SZT on eczema using systems pharmacology, through combined evaluation of absorption, distribution, metabolism, and excretion (ADME), compound feature mapping, drug target mining, target enrichment network and pathway analyses, and in silico molecular docking. The whole workflow of the study is shown in Figure 1. Results of these analyses will contribute to the understanding of how SZT formula functions, on the molecular level, in the treatment of eczema, providing the rationale and tools to improve the formulation for better TCM treatment efficacy.

2. Materials and Methods

2.1. Building of the Compound Database

The SZT formula is composed of four Chinese herbs: Sophorae flavescentis radix (Kushen, K), Schizonepetae herba (Jingjie, J), Arctii fructus (Niubangzi, N), and Atractylodis macrocephalae rhizoma (Baizhu, B). The compounds of these four herbs were collected from the TCM Systems Pharmacology (TCMSP) database [24].

2.2. Active Compound Prediction by Oral Bioavailability and Drug-Likeness Evaluation

Oral bioavailability (OB), one of the vital properties of drugs, is the percentage that is absorbed into the system from an orally administered dose of drugs. The compounds with OB ≥ 30% (20% variability making threshold up to 50%) in TCMSP were selected as active compounds considering the average of OB (30% with 10%–50% variability) in clinical studies [23].

Drug-likeness (DL) is the similarity of a given compound with the physiochemical or/and structural properties of existing drugs, which have been used to evaluate a drug’s potential early in the process of discovery. A Tanimoto coefficient was developed to predict the DL value of compounds [25]. The formula is as follows:where A represents the molecular properties of herb compounds and B represents the average molecular descriptor of molecules from the DrugBank database (drugbank.ca). DL ≥ 0.18 was determined as the criterion to select potential active compounds, where 0.18 was the mean value of DL index in DrugBank [26].

2.3. Compound Feature Mapping

The pharmacological activities of a compound are directly related to their physicochemical properties. The physicochemical parameters included in Lipinkski’s Rule of Five, including molecular weight, lipohydro partition coefficient, hydrogen donor and acceptor, and rotatable bonds, have been used to evaluate druggability of compounds [27]. Thus, estimation of relevant physicochemical properties of active compounds may provide information about their drug likeness. Thus, the compound feature mapping method [23] was applied to visualize all active compounds from SZT according to the following six physicochemical properties, molecular weight (MW), hydrogen bond donors (nHDon), hydrogen bond acceptor (nHAcc), octanol/water partition coefficient (AlogP), topological polar surface area (TPSA), and number of rotatable bonds (RBN). All six properties were analyzed by principal component analysis (PCA), where six-dimensional data were reduced to two-dimensional data with acceptable former descriptor variation in order to best display distribution of active ingredients in chemical space.

2.4. Target Fishing and Classification

First, a simplified molecular-input line-entry system (SMILES) of potential active compounds was prepared by ALOGPS 2.1 program (vcclab.org/web/alogps/). Target predictions were conducted using the SMILES of potential compounds through several published databases, including TCMSP [24], HitPick (mips.helmholtz-muenchen.de/hitpick/), Swiss Target Prediction [28] (swisstargetprediction.ch/), Similarity Ensemble Approach (SEA, sea.bkslab.org/) [29], PubChem (pubchem.ncbi.nlm.nih.gov/) [30], and DrugBank [31]. All the predicted targets obtained from these databases were mined in prevailing databases as follows to select targets for eczema: therapeutic target database (TTD, db.idrblab.net/ttd/) [32], genetic association database (GAD, geneticassociationdb.nih.gov/), DisGeNet (disgenet.org/) [33], and Open Targets Platform (targetvalidation.org/) [34]. Selected targets were finally mapped to the UniProt Database (uniprot.org/) [35] for normalization.

2.5. Gene Ontology and Pathway Analysis

Target enrichment gene ontology (GO) and pathway analyses lead us to gain mechanistic insight from the molecular level into biological function level. GO was introduced by mapping targets to the DAVID database (david.ncifcrf.gov/) [36]. The GO biological process terms with false discovery rate (FDR) < 0.01 were selected. Pathways were obtained by mapping targets to KOBAS 3.0 (kobas.cbi.pku.edu.cn) [37]. The significant pathways with FDR <0.01 were selected.

2.6. Network Construction and Analysis

With predicted targets of compounds and significant pathways, we built two biological networks using Cytoscape (v3.2.1). The compound-target-disease (C-T-D) network, containing active compounds, their related targets, and disease, provides general information about the pharmacological mechanism of SZT formula in the molecular level. The principal pathways, extracted with its related targets and active compounds, link the pathways to candidate compounds through targets and further interweave a compound-target-pathway-disease (C-T-P-D) network. The properties of these two networks were validated by NetworkAnalyzer [38], a plugin of Cytoscape.

2.7. Molecular Docking

To explore the binding modes and offer more insights into the interaction between molecular targets and compounds, molecular docking was performed on hub targets and all related bioactive compounds by AutoDock Vina [39]. Protein crystal structures with excellent resolution were downloaded from RCSB protein data bank (rcsb.org/) [40]. The structures of the ligand were directly downloaded from PubChem (pubchem.ncbi.nlm.nih.gov/) [30] without further optimization. Proteins and compounds were prepared by AutoDockTools (v1.5.6) [41]. The three-dimensional molecular graphics were prepared by the PyMOL system [42] (pymol.org) and Discovery Studio [43]. Generally, all hydrogens and Gasteriger charges were added to each molecule. Docking areas and Autogrid parameters were set based on the binding pockets of proteins. The lowest binding energy was selected and illustrated as the best binding conformation.

3. Results and Discussion

3.1. Active Compounds Obtained from SZT

In total, 441 compounds from the four herbs of SZT were obtained from TCM Systems Pharmacology (TCMSP), of which 159 were from Jingjie, 113 from Kushen, 55 from Baizhu, and 114 from Niubangzi. Such complexity of SZT makes it difficult to explore its dominant compounds and uncover the biological mechanism of formula. Thus, it is necessary to determine potentially active ingredients’ database from the SZT formula. Many compounds found in TCM, while being biologically active, may lack the absorption and bioavailability necessary to overcome barriers during the oral administration process and reach the targets, making it difficult to investigate the mechanisms of herbs. Poor oral bioavailability (OB) may result in low efficacy and cause a new drug to fail in clinical trials. High OB is a critical parameter to evaluate the possibility of a molecule becoming a drug [44]. We used oral bioavailability (OB) and drug-likeness (DL) as criteria to select the compounds as follows. OB ≥ 30% and DL ≥ 0.18 were used as threshold to narrow down the compound scope, based on the theory that drugs require proper physicochemical properties to cross various barriers and play their pharmacological effects. This threshold was determined on account of not only extracting sufficient information from formula but also keeping the balance of active compound numbers from different herbs.

As a result, 51 active compounds from 441 were selected for further analysis. Abbreviation, chemical name, molecular structure, property, and species of each active compound are listed in Table 1, which included flavonoid (25.5%), alkaloid (25.5%), phenylpropanoid (15.7%), terpene (13.7%), steroid (11.7%), and others (7.9%). Selected compounds with satisfactory OB and DL value were assumed as potential active compounds. Among them, 25 compounds (OB ≥ 50% and DL ≥ 0.18) are from 113 compounds of herb Kushen. For example, matrine (K15, OB = 63.77% and DL = 0.25), one of the major tetracycloquinolizindine alkaloids, was reported to exhibit anti-inflammatory [45], antiviral [46], and antiallergic properties [47]. Flavone norkurarinol (K2, OB = 51.28% and DL = 0.64) showed various biological activities such as antioxidant, antibacterial, anti-influenza, and anti-inflammatory activities [48]. For herb Jingjie, 11 of 159 compounds met the filter criteria (OB ≥ 30% and DL ≥ 0.18). Luteolin (J1, OB = 36.16% and DL = 0.25), a flavone present in many herbs, has been proven to inhibit inflammation responses [49]. Quercetin (J2, OB = 46.43% and DL = 0.28), a well-known free-radical scavenger, can stimulate the immune system by decreasing the production of proinflammatory cytokines, suppressing IL-4 production, improving Th1/Th2 balance, and restraining the formation of antigen-specific IgE antibody [50]. In addition, phytosterols, such as β-sitosterol (J3, OB = 36.91% and DL = 0.75), sitosterol (J4, OB = 36.91% and DL = 0.75), and stigmasterol (J5, OB = 43.83% and DL = 0.76) have been the subject of increased interest for potential treatment of atopic dermatitis (AD) to replace steroid therapy [51]. For herb Niubangzi, 8 active compounds from 114 were identified to show OB ≥ 30% and DL ≥ 0.18. Especially, arctiin (N2, OB = 34.45% and DL = 0.84), one of the major active ingredients, can be biotransformed to its more potent aglycone arctigenin by intestinal microbiota [52]. Cynarin (N8, OB = 31.75% and DL = 0.68), an ester of quinic acid and two molecules of caffeic acid, have been proved to display good antioxidant, antiradical, and anticholinergic effects [53]. Moreover, β-sitosterol (J3, N3, OB = 36.91% and DL = 0.75) and supraene (J6, N5, OB = 33.55% and DL = 0.42) are shared active compounds by both Niubangzi and Jingjie. For herb Baizhu, 7 active compounds from 55 ingredients were selected to meet the criterion (OB ≥ 30% and DL ≥ 0.18). Among them, α-amyrin (B3, OB = 39.51% and DL = 0.76), a pentacyclic triterpene, shows potential anti-inflammatory effects [54]. In particular, it significantly suppresses the scratching behavior in a mouse model of pruritus by inhibiting degranulation of mast cell [55], beneficial for skin lesions in eczema.


HerbComp. IDComp. nameMol. Wt.O.B.D.L.StructureSp.

KushenK1kushenolJ_qt286.350.860.24F
KushenK2Norkurarinol442.5551.280.64F
KushenK3KushenolJ580.5951.390.74F
KushenK4(+)-Lupanine248.4152.710.24A
KushenK5Norartocarpetin286.2554.930.24F
KushenK6Sophranol264.4155.420.28A
KushenK7(+)-Lehmannine246.3958.340.25A
KushenK8(+)-Allomatrine248.4158.870.25A
KushenK9Sophoridine248.4160.070.25A
KushenK10Leachianone,g356.460.970.4F
KushenK11Isosophocarpine246.3961.570.25A
KushenK12Anagyrine244.3762.010.24A
KushenK13Leontalbinine246.3962.080.25A
KushenK14Lehmanine246.3962.230.25A
KushenK15Matrine248.4163.770.25A
KushenK16Sophocarpine246.3964.260.25A
KushenK1713,14-Dehydrosophoridine246.3965.340.25A
KushenK18Inermin284.2865.830.54P
KushenK19cis-Dihydroquercetin304.2766.440.27F
KushenK20Isomatrine248.4168.680.25A
KushenK21Formononetin268.2869.670.21F
KushenK22K22256.2771.120.18F
KushenK23Inermine284.2875.180.54P
KushenK24Phaseolin322.3878.20.73P
KushenK25Glyceollin338.3897.270.76P
JingjieJ1Luteolin286.2536.160.25F
JingjieJ2Quercetin302.2546.430.28F
JingjieJ3beta-Sitosterol414.7936.910.75S
JingjieJ4Sitosterol414.7936.910.75S
JingjieJ5Stigmasterol412.7743.830.76S
JingjieJ6Supraene410.833.550.42T
JingjieJ7Diosmetin300.2831.140.27F
JingjieJ8Campest-5-en-3-beta-ol400.7637.580.71S
JingjieJ9J9302.347.740.27F
JingjieJ10Schizonepetoside B330.4231.020.28T
JingjieJ11Schkuhrin I420.554.450.52O
NiubangziN1Neoarctin A742.8839.990.27P
NiubangziN2Arctiin534.6134.450.84P
NiubangziN3beta-Sitosterol414.7936.910.75S
NiubangziN4Kaempferol286.2541.880.24F
NiubangziN5Supraene410.833.550.42T
NiubangziN6beta-Carotene536.9637.180.58T
NiubangziN7N7386.4852.30.48P
NiubangziN8Cynarin(e)516.4931.760.68P
BaizhuB1B1276.4135.950.21T
BaizhuB2B2428.8236.230.78S
BaizhuB3α-Amyrin426.839.510.76T
BaizhuB43β-Acetoxyatractylone274.3954.070.22T
BaizhuB5B5355.4460.310.31O
BaizhuB6B6312.3962.40.22O
BaizhuB7B7356.4563.370.3O

K22: (2R)-7-hydroxy-2-(4-hydroxyphenyl)chroman-4-one; J9: 5,7-dihydroxy-2-(3-hydroxy-4-methoxyphenyl)chroman-4-one; N7: (3R,4R)-3,4-bis[(3,4-dimethoxyphenyl)methyl]oxolan-2-one; B1: 8β-ethoxy atractylenolide III; B2: (24S)-24-propylcholesta-5-ene-3beta-ol; B5: 14-acetyl-12-senecioyl-2E,8E,10E-atractylentriol; B6: 12-senecioyl-2E,8E,10E-atractylentriol; B7: 14-acetyl-12-senecioyl-2E,8Z,10E-atractylentriol.
3.2. Compound Feature Mapping Based on Physicochemical Parameters

The drug-like physicochemical properties of active compounds were investigated by analyzing six common drug-associated physicochemical parameters including molecular weight (MW), hydrogen bond donors (nHDon), hydrogen bond acceptor (nHAcc), octanol/water partition coefficient (AlogP), topological polar surface area (TPSA), and number of rotatable bonds (RBN). Principal component analysis (PCA) was applied to display the distribution of active compounds in chemical space considering all six physicochemical parameters. The first principle component (PC1) displays the direction of maximal variance which takes 51.92% of total variance, and PC2 represents the second largest variations that account for 37.95% (Figure 2). These two principle components account for 89.87% of the total variances, indicating that they sufficiently represent the six physicochemical features of these compounds. Figure 2 illustrates the relatedness distribution of the active compounds. Most of the active compounds from SZT are clustered in the red circled area, indicating they have similar physicochemical features. It is noted that the binding of compounds to their targets are up to the structure and property of compounds. The result implies that the targets of these active compounds, derived from different herbs but with similar physicochemical properties, may be associated and overlapping.

3.3. Potential Target Fishing

Establishing compound-target interactions and enriching those in signaling pathways become increasingly necessary for explicating the mechanism of drug action and their underlying pharmacological effects. The targets of 51 active compounds were fished by several popular databases, including HitPick, Swiss Target Prediction, Similarity Ensemble Approach (SEA), PubChem, and DrugBank. The integrated targets were further mined in targets selected from the therapeutic target database (TTD), genetic association database (GAD), DisGeNet, and Open Targets Platform with terms “eczema” and “atopic dermatitis.” Finally, a total of 81 potential targets were screened out for 51 active compounds with 327 interactions. Most compounds were found to interact with more than one target for eczema, which implies comprehensive regulation and extensive pharmacological actions of SZT. For example, phaseolin (K24) from Kushen interacts with 15 potential targets related to eczema, while kaempferol (N4) from Niubangzi links to 12 potential targets. Among these targets, vascular endothelial growth factor A (VEGFA) [56], epidermal growth factor receptor (EGFR) [57], and interleukin 2 (IL-2) [58] are believed to play crucial roles during the pathogenesis of AD. β-sitosterol (N3, J3) from both Niubangzi and Jingjie targets 11 proteins, among which peroxisome-proliferator-activated receptor (PPAR) family and retinoic-acid-receptor-related orphan receptors (RORs) might be the main regulated targets of β-sitosterol [59]. RORs regulate the function of specific subsets of T cells and innate lymphoid cells, which are key drivers of inflammatory disease in barrier tissues [60]. Moreover, the retinoid X receptor (RXR) is an intriguing and essential member of nuclear receptors activated by 9-cis retinoic acid to regulate cell differentiation, metabolism, and cell apoptosis [61]. The strong reduction of retinoid signaling and retinoid concentration in affected as well as nonaffected skin of individuals with AD has been discovered, compared to healthy individuals [62]. A selective activator of RAR or RXR may influence the retinoid signaling in skin, and RORA and RORC have been reported to be a potential critical therapeutic target for treating skin diseases of inflammatory etiology [63].

C-T-D network construction is used to display the active compounds and potential targets of SZT for eczema treatment.

3.3.1. Network Analysis

Active compound and target data extracted above were employed to construct a C-T-D network containing 129 nodes and 406 edges, involving 47 active compounds and 81 potential targets. Four active compounds K11, K25, B5, and B6 shared no targets with eczema and were, thus, excluded. The network demonstrates multi-interactions between compounds and targets (Figure 3). The diamonds and circles represent active compounds from SZT and their potential targets, respectively. Of 81 targets, Kushen, Jingjie, Niubangzi, and Baizhu recognized 52, 35, 43, and 19 targets, respectively. To evaluate the importance of each node in this network, the degree of each node was defined as the number of edges connected with this node, representing the influence of the node. The network average degree of node was 6.3.

3.3.2. Target Analysis

Most of the targets from our network construction are closely associated with the pathogenesis of eczema. We believe that targets with highest degree are the most vital ones regulated by SZT, thus chosen as hub targets for further analysis. Interestingly, carbonic anhydrase II (CA2) had the largest degree, hit by 17 compounds (5 of Kushen, 3 of Niubangzi, 7 of Jingjie, and 1 of Baizhu). CA2 monitors pH regulation, water transport, and hydration of CO2, which is related to diverse diseases such as glaucoma, tumor, epilepsy, and diabetes [64]. CA2 inhibitors, such as acetazolamide and brinzolamide, have been developed and established to treat glaucoma, seizure disorder, and acute mountain sickness [65]. It was reported that CA2 was upregulated by Th2 cytokines in lesion skin of AD patients, which caused a pH rise in skin [66]. The higher surface pH may affect the integrity and cohesion of skin and compromise barrier function of skin in AD [67]. Thus, the CA2 inhibitor may relieve a series of symptoms in dermatitis triggered by itching [68]. Moreover, VEGFA was targeted by 12 active compounds (6 of Kushen, 4 of Jingjie, and 2 of Niubangzi). VEGFA displays a key role in vascular permeability, vasodilation, and angiogenesis. VEGFA can be secreted by various immunocytes in inflamed skin, such as mast cells, macrophages, eosinophils, basophils, and Th17 cells, which appear critical not only in lymphatic vessel expansion but also in antigen clearance and inflammation resolution through enhancement of lymphangiogenesis [56]. Meanwhile, 12 compounds from SZT were found to interact with PPARG. PPARG belongs to a subfamily of nuclear hormone receptors and plays a distinct physiological role in regulating expression of genes involved in cellular proliferation, specific components of the Th2 inflammatory pathway, and maintenance of the skin barrier [69].

Gene ontology (GO) and pathway analysis reveal potential regulation of SZT in the immune, inflammatory, and metabolic processes.

Utilizing the identified targets as an enriched gene set, we employed gene-set enrichment analysis to find GO biological process terms in the DAVID database and pathways in the KOBAS database. The top 15 biological process GO terms with False Discovery Rate (FDR) < 0.01 are ranked by enrichment score (–logFDR) in Figure 4(a). FDR is to conceptualize the rate of type I errors in null hypothesis testing when conducting multiple comparison [70]. Most targets are found to be closely related to several biological processes, such as the steroid-hormone-mediated signaling pathway, inflammatory response, retinoic acid receptor signaling pathway, intracellular receptor signaling pathway, and cellular response to lipopolysaccharide. Moreover, most biological process terms listed have been shown to be highly associated with the immunity, inflammation, and/or eczema. In affected or nonaffected skin of patients with AD, the retinoid transport, synthesis, concentration, and signaling were reported to be strongly decreased, suggesting the intrinsic influence of the retinoid signaling pathway to the pathological process of AD [71]. Particularly, the steroid-hormone-mediated signaling pathway with hub targets, such as PPARG, RXRA, RORC, and VDR, is an important regulatory pathway in immunity and inflammation [72]. Inflammatory response, positive regulation of the nitric oxide biosynthetic process, and response to lipopolysaccharide are strongly associated with the inflammatory reaction of the body. The GO results imply that SZT might treat eczema by integral regulating the immune and inflammatory biological processes.

A total of 137 pathways with FDR <0.01 were obtained in KOBAS. The top 15 most significant pathways are listed in Figure 4(b). Most pathways are strongly associated with the immune process, metabolic process, and complex pathways in combination of both. For instance, Th17 cell differentiation, inflammatory bowel disease, tuberculosis, kaposi-sarcoma-associated herpesvirus infection, and human cytomegalovirus infection are more inclined to the immune response induced by internal and/or external stimuli, while the AGE-RAGE signaling pathway in diabetic complications, metabolic pathways, and adipocytokine signaling pathway are more prone to the mediation of the metabolic process in the body. Moreover, both pathways in cancer and the PI3K-Akt signaling pathway are complex and are related to both immune and metabolic processes. Compared with biological process GO analysis, KEGG pathway enrichment provided more detailed information about the signal transmission during biological processes. Thus, key pathways were selected to establish further the C-T-P-D network to determine the relatively critical compounds and hub targets.

A C-T-P-D network was established to reveal the relationships between compounds, targets, pathways and disease.

We speculate that hub targets with high degree in the network are the most promising targets regulated by SZT. Thus, based on the number of hub targets contained in these pathways, four key pathways were selected, including the Th17 cell differentiation pathway (P1), pathways in cancer (P2), metabolic pathways (P3), and PI3K-Akt signaling pathway (P4). Then, the compounds and targets involved with these four pathways were chosen to construct the C-T-P network, as shown and listed in Figure 5(a). Diamonds, circles, triangles, and hexagons represent bioactive compounds, corresponding targets, pathways, and disease eczema, respectively. Node color from green to red and node size are proportional to its degree in the C-T-P-D network. The network average degree of node was 7.6, with the degree defined as the number of connections. Moreover, detailed information including the involved hub targets, related pathways, and potential ligands is displayed in Figure 5(b). The imbalance of Th1 and Th2 in the pathogenesis of eczema has been widely investigated [73]. Recently, the crucial role of Th17 cell in allergic contact dermatitis has also been demonstrated by several studies. Increasing Th17 activation was confirmed in Asian and pediatric AD [74, 75]. Once activated, Th17 cells have the capacity to produce IL-17A, IL-17F, IL-22, and IL-26. In particular, IL-17 leads to the expression of proinflammatory cytokines, responsible for eosinophil- and neutrophil-mediated inflammation [76]. We found that 16 targets were involved in the Th17 cell differentiation pathway, including hub targets IL-2, RORC, RORA, and RXRA. Among them, IL-2 is released by activated Th1 cells to regulate inflammation and tissue damage. The regulatory roles of RORs and RXRs in the function of lymphoid cell have already been discussed previously [60, 61]. Pathways in cancer, including the MAPK signaling pathway, Jak-STAT signaling pathway, calcium signaling pathway, estrogen signaling pathway, and VEGF signaling pathway, were all identified in the pathway analysis. The hub targets, such as Prostaglandin-Endoperoxide Synthase 2 (PTGS2), EGFR, RXRA, PPARG, VEGFA, and IL-2, were all included in this complex pathway, pathway in cancer. PTGS2, also known as cyclooxygenase-2 (COX-2), is currently the widely used anti-inflammation target of nonsteroidal anti-inflammatory drugs. The inhibition of PTGS2 results in preferable anti-inflammatory effects [77]. In addition to immune function, the association between AD and metabolic syndrome has attracted more attention [78]. Metabolic syndrome includes abdominal obesity, hypertension, insulin resistance, and dyslipidemia [79]. Collective data indicate that central obesity is positively associated with AD, which is likely linked to low-grade systemic inflammation induced by adipocytes (leptin and adiponectin) generated by visceral adipose tissue [80]. CA2, VDR, PTGST1, and PTGST2 are all included in metabolic pathways. The PI3K-Akt signaling pathway is involved in a wide variety of cellular processes, including cellular growth, migration, and proliferation, especially the process of inflammation and immune response. It was reported that the expression of PI3K and Akt in blood was significantly higher in AD patients than in the healthy controls, which might be associated with the intrinsic activation of T-cell and cytokine secretion [81]. Moreover, specific inhibition of PI3K expression led to significant inhibition of T-cell proliferation and secretion of cytokines, such as IL-6 and IL-10. The PI3K-Akt signaling pathway has been reported to be associated with LPS-induced acute inflammatory responses and autophagy regulation in immune response [19].

3.4. Molecular Docking Analysis Investigates the Binding Activities

To investigate binding between potential targets and active compounds, in silico molecular docking was used to calculate the binding energy and evaluate the binding mode between the compound and its target. From hub targets, CA2 (P3), PPARG (P2), and RXRA (P1, P2, and P4) were selected to represent the four different pathways. Additionally, activation of VDR causes marked induction of skin barrier genes and antimicrobial peptide genes in lesion skin, ameliorating allergen-triggered eczema in the murine AD model [82]. Also, it is a potentially novel therapeutic target for treatment of eczema. Thus, molecular docking was conducted for CA2, PPARG, RXRA, and VDR and their interacting active compounds using AutoDock Vina, considering proteins as rigid molecules. For CA2, promising ligands are docking into receptor-binding-domain as inhibitors. Also, for PPARG, RXRA, and VDR, promising ligands are calculated as their agonists. The ligands and their molecular dockings were demonstrated in Figure 6. The results revealed that all compounds as predicted ligands show moderate to strong binding affinity with related proteins. In particular, compounds cynarine (N8), stigmasterol (J5), (24S)-24-propylcholesta-5-ene-3β-ol (B2), kushenol (K1), and β-sitosterol (N3) had the strongest binding affinities to their receptors. The binding modes of complexes CA2-N8 (−9.1 kcal/mol), PPARG-J5 (−9.3 kcal/mol), PPARG-B2 (−9.2 kcal/mol), RXRA-K1 (−8.1 kcal/mol), and VDR-N3 (−10.6 kcal/mol) are shown in Figure 7.

For CA2, 15 ligands were docked into the CA2 active site and the cynarine (N8) showed highest binding affinity. The hydrophobic pocket formed by Val121, Val143, and Leu198 of CA2 were also observed in the binding of CA2-N8. The hydrogen bond formed between the ligand and Thr199, Thr200, Asn67, and Asp72 further increases the stability of the ligand in the binding site (Figures 7(a) and 7(b)). For PPARG, both stigmasterol (J5) and (24S)-24-propylcholesta-5-ene-3β-ol (B2) were well fitted into the binding cavity of the protein by hydrophobic interaction. The binding pocket of PPARG-stigmasterol (J5) shown in Figures 7(c) and 7(d) indicates that nine amino acid residues from PPARG interacted with J5 by hydrophobic interactions, including Leu330, Ile326, Val339, Arg288, Cyc285, Ile341, Leu270, and Ile262. Both stigmasterol and (24S)-24-propylcholesta-5-ene-3β-ol are phytosterols with similar chemical structure. The binding energy of PPARG-B2 (Figures 7(e) and 7(f)) is −9.2 kcal/mol, very close to that of PPARG-stigmasterol (J5) (−9.3 kcal/mol). Moreover, the binding of PPARG-N3 (β-sitosterol) has previously been studied [83] with similar affinity. For RXRA, docking ligands into the ligand-binding domain was used to evaluate the possibility of these compounds acting as agonists. The strongest binding existed between RXRA and kushenol (K1). Hydrogen binding (Cys432), π-π stacking (Phe313), and hydrophobic interactions (Ile268, Ala271, Ala272, and Leu309) formed between protein residues and ligand kushenol (K1) were illustrated in Figures 7(g) and 7(h). VDR is a member of nuclear receptor superfamily of transcriptional regulators, which can be activated by its specific ligand vitamin D3. The docking result indicated that binding conformation of β-sitosterol (N3) fitted well in binding domain of VDR crystallographical structure forming stable complex VDR-N3 (β-sitosterol). All 14 amino acid residues showed hydrophobic interactions with ligand N3 (β-sitosterol) (Figures 7(i) and 7(j)). Their forceful interaction might attribute to the similar structure of vitamin D3 and β-sitosterol. Inspiringly, it has been reported that β-sitosterol promoted the immune function of vitamin D3 [84], suggesting the anti-inflammation and immune enhancement function of β-sitosterol. The abovementioned strong interactions further imply the functioning route of SZT compounds through proteins and their corresponding biological activities.

4. Conclusions

Eczema is a chronic inflammatory condition causing pruritus, erythema, dryness, scaling, and vesiculopapular rash that negatively impact the quality of life of the affected patients [1]. With no definite cure for eczema, common medicines, such as oral antihistamines, topical corticosteroids, emollients, and even immune-modulating agents, can not only relieve symptoms but also sometimes cause severe side effects [85]. Significantly, these agents do not target upstream molecular dysregulation that is central to eczema. There is an urgent need to develop new treatments to not only relieve symptoms but also address the underlying molecular mechanisms that drive the disease to potentially cure eczema with fewer side effects. TCM has a long and well-known history in Asia treating diseases in a multicompound-multitarget modality [86]. TCM as a complementary or alternative medicine is gradually becoming more widely accepted throughout the world, and molecular basis of its activity is only recently being elucidated [87]. XFS is an acclaimed TCM formula in China documented in “Orthodox Manual of External Medicine” [13], which eliminates itching in the body. We modified the formula to retain four pivotal herbs to develop SZT that has excellent efficacy in treatment of eczema. However, molecular mechanisms of SZT remain poorly defined. Recently, systems pharmacology based on network analysis has provided a new method to uncover the mechanisms of complex formula in TCM [1618]. We explored the molecular mechanisms of SZT for eczema treatment using a systems pharmacology approach, including ADME evaluation, herb feature mapping, drug target mining, network and pathway analyses, and in silico molecular docking.

In summary, we found 51 active compounds out of 441 compounds that were retrieved based on OB and DL screening with satisfactory pharmacokinetic properties. We further studied these active compounds for their physicochemical properties and drug-likeness. Findings of drug target mining supported that active compounds targeted single or multiple proteins involved in the pathological process or treatment of eczema. The C-T-D network analysis, GO, and pathway analyses conducted using DAVAD and KOBAS databases allowed us to elucidate the pharmacological effects of SZT from a holistic perspective. Our results imply that SZT efficacy in eczema may be attributed to immune and metabolic functions via regulation of multiple pathways and networks. Among these, four key pathways stood out, including Th17 cell differentiation, pathways in cancer, metabolic pathways, and PI3K-Akt signaling pathway. Moreover, the C-T-P-D network highlighted the critical active compounds. Finally, molecular docking studies between active compounds and proteins CA2, PPARG, RXRA, and VDR confirmed moderate to strong binding affinities. These data further support the potential biological activities of these active compounds. Overall, the formula affects a regulation of the immune system and metabolic processes, instead of a single target or a single pathway. In vitro and in vivo studies are currently in progress to investigate, confirm, and expand knowledge into the molecular alterations induced by SZT active compounds. In conclusion, systems pharmacology was applied in this work to shed light on the molecular mechanism and pharmacological regulation of SZT for eczema treatment. Our results will help further optimize SZT and facilitate development of more effective drugs to cure eczema.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

ZZW, YL, WH, RT, ANW, JG, and MM have no financial conflicts of interest to disclose. KDS receives salary from General Nutraceutical Technology LLC. XML received research support from the National Institutes of Health (NIH)/National Center for Complementary and Alternative Medicine (NCCAM), Food Allergy Research and Education (FARE) and Winston Wolkoff Integrative Medicine Fund for Allergies and Wellness, the Parker Foundation, and Henan University of Chinese Medicine; received consulting fees from Food Allergy Research and Education (FARE), Johnson & Johnson Pharmaceutical Research & Development, LLC, and Bayer Global Health LLC; royalties from UpToDate; travel expenses from the National Center for Complementary and Alternative Medicine (NCCAM) and FARE, Henan University of Chinese Medicine, Harvard TCM summit, China Allergy Society Conference, and North American TCM Symposium; shares US patents PCT/US05/008417, PCT/US 10,500,169, PCP/14/762,416, and PCP 14/762,416; is a member of Herbs Springs, LLC, General Nutraceutical Technology LLC, and Health Freedom LLC; and takes compensation from practice at Integrative Health and Acupuncture PC.

Acknowledgments

This research was funded by the Henan University of Chinese Medicine, visiting fellowship grant to ZZW (CTIM-01-2019), and Scientific Research Foundation (KYQD 2019-05). This research was also supported by the Henan Province Scientific and Technological Project (no. 202102310472) and Study of Integrative Medicine forImmunology and Wellness, New York Medical College. The authors thank Henry Ehrlich for critical reading of the manuscript; Wei Zhou, Lihua Cao and Hongjuan He, Nan Yang, Anish R Maskey, Ke Li, and Ibrahim A Musa for helpful discussion.

References

  1. N. Novak, “New insights into the mechanism and management of allergic diseases: atopic dermatitis,” Allergy, vol. 64, no. 2, pp. 265–275, 2009. View at: Publisher Site | Google Scholar
  2. T. E. Shaw, G. P. Currie, C. W. Koudelka, and E. L. Simpson, “Eczema prevalence in the United States: data from the 2003 national survey of children’s health,” Journal of Investigative Dermatology, vol. 131, no. 1, pp. 67–73, 2011. View at: Publisher Site | Google Scholar
  3. J. I. Silverberg, “Health care utilization, patient costs, and access to care in US adults with eczema,” JAMA Dermatology, vol. 151, no. 7, pp. 743–752, 2015. View at: Publisher Site | Google Scholar
  4. P. D. Arkwright, C. Motala, H. Subramanian, J. Spergel, L. C. Schneider, and A. Wollenberg, “Management of difficult-to-treat atopic dermatitis,” The Journal of Allergy and Clinical Immunology: In Practice, vol. 1, no. 2, pp. 142–151, 2013. View at: Publisher Site | Google Scholar
  5. R. J. Hay, N. E. Johns, H. C. Williams et al., “The global burden of skin disease in 2010: an analysis of the prevalence and impact of skin conditions,” Journal of Investigative Dermatology, vol. 134, no. 6, pp. 1527–1534, 2014. View at: Publisher Site | Google Scholar
  6. C. J. Apfelbacher, A. Jupiter, B. Carter, E. Weisshaar, and M. M. Böhmer, “Oral H1 antihistamines as ‘add‐on’ therapy to topical treatment for eczema,” Cochrane Database of Systematic Reviews, vol. 4, 2016. View at: Google Scholar
  7. R. Nowicki, M. Trzeciak, A. Wilkowska et al., “Special paper atopic dermatitis: current treatment guidelines. Statement of the experts of the dermatological section, polish society of allergology, and the allergology section, polish society of dermatology,” Advances in Dermatology and Allergology, vol. 4, no. 4, pp. 239–249, 2015. View at: Publisher Site | Google Scholar
  8. X.-M. Li and H. Ehrlich, Traditional Chinese Medicine, Western Science, and the Fight against Allergic Disease, World Scientific Publishing, Singapore, 1st edition, 2016.
  9. K. L. E. Hon, T. F. Leung, P. C. Ng et al., “Efficacy and tolerability of a Chinese herbal medicine concoction for treatment of atopic dermatitis: a randomized, double-blind, placebo-controlled study,” British Journal of Dermatology, vol. 157, no. 2, pp. 357–363, 2007. View at: Publisher Site | Google Scholar
  10. S. Gu, A. W. H. Yang, C. C. L. Xue et al., “Chinese herbal medicine for atopic eczema,” Cochrane Database of Systematic Reviews, vol. 9, 2013. View at: Google Scholar
  11. X.-R. Lin, C.-X. Tu, X.-M. Meng, C.-M. Yang, M.-Y. Gao, and L. Gu, “Studies on treating eczema by Chinese herbal medicine with anti-type IV allergic activity,” Chinese Journal of Integrated Traditional and Western Medicine, vol. 7, no. 1, pp. 7–11, 2001. View at: Google Scholar
  12. D. Huang, K. Chen, F.-R. Zhang et al., “Efficacy and safety of Run Zao Zhi Yang capsule on chronic eczema: a multiple-center, randomized, double-blind, placebo-controlled clinical study,” Journal of Dermatological Treatment, vol. 30, no. 7, pp. 677–684, 2019. View at: Publisher Site | Google Scholar
  13. Y.-H. Lin, Y.-C. Chen, S. Hu, H.-Y. Chen, J.-L. Chen, and S.-H. Yang, “Identifying core herbal treatments for urticaria using Taiwan’s nationwide prescription database,” Journal of Ethnopharmacology, vol. 148, no. 2, pp. 556–562, 2013. View at: Publisher Site | Google Scholar
  14. H.-Y. Chen, Y.-H. Lin, S. Hu, S.-h. Yang, J.-l. Chen, and Y.-C. Chen, “Identifying Chinese herbal medicine network for eczema: implications from a nationwide prescription database,” Evidence-Based Complementary and Alternative Medicine, vol. 2015, Article ID 347164, 9 pages, 2015. View at: Publisher Site | Google Scholar
  15. E. Thanik, J. A. Wisniewski, A. Nowak-Wegrzyn, H. Sampson, and X.-M. Li, “Effect of traditional Chinese medicine on skin lesions and quality of life in patients with moderate to severe eczema,” Annals of Allergy, Asthma & Immunology, vol. 121, no. 1, pp. 135-136, 2018. View at: Publisher Site | Google Scholar
  16. J. Hansen and R. Iyengar, “Computation as the mechanistic bridge between precision medicine and systems therapeutics,” Clinical Pharmacology & Therapeutics, vol. 93, no. 1, pp. 117–128, 2013. View at: Publisher Site | Google Scholar
  17. S. Zhao and R. Iyengar, “Systems pharmacology: network analysis to identify multiscale mechanisms of drug action,” Annual Review of Pharmacology and Toxicology, vol. 52, no. 1, pp. 505–521, 2012. View at: Publisher Site | Google Scholar
  18. M. Danhof, “Systems pharmacology—towards the modeling of network interactions,” European Journal of Pharmaceutical Sciences, vol. 94, pp. 4–14, 2016. View at: Publisher Site | Google Scholar
  19. M. Zhao, Y. Chen, C. Wang et al., “Systems pharmacology dissection of multi-scale mechanisms of action of Huo-Xiang-Zheng-Qi formula for the treatment of gastrointestinal diseases,” Frontiers in Pharmacology, vol. 9, p. 1448, 2019. View at: Publisher Site | Google Scholar
  20. X. Li, X. Xu, J. Wang et al., “A system-level investigation into the mechanisms of Chinese Traditional Medicine: compound Danshen Formula for cardiovascular disease treatment,” PLoS One, vol. 7, no. 9, Article ID e43918, 2012. View at: Publisher Site | Google Scholar
  21. Y. Liu, H. Xing, X. Jiang, Y. Chen, M. Huang, and S. Yu, “Network pharmacology-based preventive effect of XZF on cutaneous toxicities induced by EGFR inhibitor,” Biomedicine & Pharmacotherapy, vol. 123, Article ID 109755, 2020. View at: Publisher Site | Google Scholar
  22. J. Li, P. Zhao, Y. Li, Y. Tian, and Y. Wang, “Systems pharmacology-based dissection of mechanisms of Chinese medicinal formula Bufei Yishen as an effective treatment for chronic obstructive pulmonary disease,” Scientific Reports, vol. 5, no. 1, Article ID 15290, 2015. View at: Publisher Site | Google Scholar
  23. W. Zhou, Z. Chen, W. Li et al., “Systems pharmacology uncovers the mechanisms of anti-asthma herbal medicine intervention (ASHMI) for the prevention of asthma,” Journal of Functional Foods, vol. 52, pp. 611–619, 2019. View at: Publisher Site | Google Scholar
  24. J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, p. 13, 2014. View at: Publisher Site | Google Scholar
  25. Y. Yamanishi, M. Kotera, M. Kanehisa, and S. Goto, “Drug-target interaction prediction from chemical, genomic and pharmacological data in an integrated framework,” Bioinformatics, vol. 26, no. 12, pp. i246–i254, 2010. View at: Publisher Site | Google Scholar
  26. H. Liu, J. Wang, W. Zhou, Y. Wang, and L. Yang, “Systems approaches and polypharmacology for drug discovery from herbal medicines: an example using licorice,” Journal of Ethnopharmacology, vol. 146, no. 3, pp. 773–793, 2013. View at: Publisher Site | Google Scholar
  27. C. A. Lipinski, F. Lombardo, B. W. Dominy, and P. J. Feeney, “Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings,” Advanced Drug Delivery Reviews, vol. 23, no. 1–3, pp. 3–25, 1997. View at: Publisher Site | Google Scholar
  28. A. Daina, O. Michielin, and V. Zoete, “SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules,” Nucleic Acids Research, vol. 47, no. W1, pp. W357–W364, 2019. View at: Publisher Site | Google Scholar
  29. M. J. Keiser, B. L. Roth, B. N. Armbruster, P. Ernsberger, J. J. Irwin, and B. K. Shoichet, “Relating protein pharmacology by ligand chemistry,” Nature Biotechnology, vol. 25, no. 2, pp. 197–206, 2007. View at: Publisher Site | Google Scholar
  30. S. Kim, J. Chen, T. Cheng et al., “PubChem 2019 update: improved access to chemical data,” Nucleic Acids Research, vol. 47, no. D1, pp. D1102–D1109, 2018. View at: Publisher Site | Google Scholar
  31. D. S. Wishart, C. Knox, A. C. Guo et al., “DrugBank: a knowledgebase for drugs, drug actions and drug targets,” Nucleic Acids Research, vol. 36, no. suppl_1, pp. D901–D906, 2007. View at: Publisher Site | Google Scholar
  32. Y. Wang, S. Zhang, F. Li et al., “Therapeutic target database 2020: enriched resource for facilitating research and early development of targeted therapeutics,” Nucleic Acids Research, vol. 48, no. D1, pp. D1031–D1041, 2019. View at: Google Scholar
  33. J. Piñero, J. M. Ramírez-Anguita, J. Saüch-Pitarch et al., “The DisGeNET knowledge platform for disease genomics: 2019 update,” Nucleic Acids Research, vol. 48, no. D1, pp. D845–D855, 2019. View at: Google Scholar
  34. D. Carvalho-Silva, A. Pierleoni, M. Pignatelli et al., “Open targets platform: new developments and updates two years on,” Nucleic Acids Research, vol. 47, no. D1, pp. D1056–D1065, 2018. View at: Publisher Site | Google Scholar
  35. T. U. Consortium, “UniProt: a worldwide hub of protein knowledge,” Nucleic Acids Research, vol. 47, no. D1, pp. D506–D515, 2018. View at: Google Scholar
  36. D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists,” Nucleic Acids Research, vol. 37, no. 1, pp. 1–13, 2008. View at: Publisher Site | Google Scholar
  37. C. Xie, X. Mao, J. Huang et al., “Kobas 2.0: a web server for annotation and identification of enriched pathways and diseases,” Nucleic Acids Research, vol. 39, no. suppl_2, pp. W316–W322, 2011. View at: Publisher Site | Google Scholar
  38. Y. Assenov, F. Ramírez, S.-E. Schelhorn, T. Lengauer, and M. Albrecht, “Computing topological parameters of biological networks,” Bioinformatics, vol. 24, no. 2, pp. 282–284, 2007. View at: Publisher Site | Google Scholar
  39. O. Trott and A. J. Olson, “AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading,” Journal of Computational Chemistry, vol. 31, no. 2, pp. 455–461, 2010. View at: Publisher Site | Google Scholar
  40. H. M. Berman, J. Westbrook, Z. Feng et al., “The protein data bank,” Nucleic Acids Research, vol. 28, no. 1, pp. 235–242, 2000. View at: Publisher Site | Google Scholar
  41. M. F. Sanner, “Python: a programming language for software integration and development,” Journal of Molecular Graphics & Modelling, vol. 17, no. 1, pp. 57–61, 1999. View at: Google Scholar
  42. W. L. DeLano, “Pymol: an open-source molecular graphics tool,” CCP4 Newsletter on Protein Crystallography, vol. 40, pp. 82–92, 2002. View at: Google Scholar
  43. DassaultSystèmesBIOVIAD, “Discovery studio,” 2020. View at: Google Scholar
  44. D. F. Veber, S. R. Johnson, H.-Y. Cheng, B. R. Smith, K. W. Ward, and K. D. Kopple, “Molecular properties that influence the oral bioavailability of drug candidates,” Journal of Medicinal Chemistry, vol. 45, no. 12, pp. 2615–2623, 2002. View at: Publisher Site | Google Scholar
  45. J. R. Guzman, J. S. Koo, J. R. Goldsmith, M. Mühlbauer, A. Narula, and C. Jobin, “Oxymatrine prevents NF-κB nuclear translocation and ameliorates acute intestinal inflammation,” Scientific Reports, vol. 3, no. 1, p. 1629, 2013. View at: Publisher Site | Google Scholar
  46. Y. Yang, J. Xiu, X. Zhang et al., “Antiviral effect of matrine against human enterovirus 71,” Molecules, vol. 17, no. 9, pp. 10370–10376, 2012. View at: Publisher Site | Google Scholar
  47. X.-M. Li and L. Brown, “Efficacy and mechanisms of action of traditional Chinese medicines for treating asthma and allergy,” Journal of Allergy and Clinical Immunology, vol. 123, no. 2, pp. 297–306, 2009. View at: Publisher Site | Google Scholar
  48. H.-M. Oh, S. W. Lee, M. H. Park et al., “Norkurarinol inhibits toll-like receptor 3 (TLR3)-mediated pro-inflammatory signaling pathway and rotavirus replication,” Journal of Pharmacological Sciences, vol. 118, no. 2, pp. 161–170, 2012. View at: Publisher Site | Google Scholar
  49. C. M. Park and Y.-S. Song, “Luteolin and luteolin-7-O-glucoside inhibit lipopolysaccharide-induced inflammatory responses through modulation of NF-κB/AP-1/PI3K-Akt signaling cascades in RAW 264.7 cells,” Nutrition Research and Practice, vol. 7, no. 6, pp. 423–429, 2013. View at: Publisher Site | Google Scholar
  50. J. Mlcek, T. Jurikova, S. Skrovankova, and J. Sochor, “Quercetin and its anti-allergic immune response,” Molecules, vol. 21, no. 5, p. 623, 2016. View at: Publisher Site | Google Scholar
  51. S.-J. Kim, “The ameliorative effect of β-sitosterol on DNCB-induced atopic dermatitis in mice,” Biomedical Science Letters, vol. 23, no. 4, pp. 303–309, 2017. View at: Publisher Site | Google Scholar
  52. Q. Gao, M. Yang, and Z. Zuo, “Overview of the anti-inflammatory effects, pharmacokinetic properties and clinical efficacies of arctigenin and arctiin from Arctium lappa L.,” Acta Pharmacologica Sinica, vol. 39, no. 5, pp. 787–801, 2018. View at: Publisher Site | Google Scholar
  53. M. Topal, H. Gocer, F. Topal et al., “Antioxidant, antiradical, and anticholinergic properties of cynarin purified from the Illyrian thistle (Onopordum illyricum L.),” Journal of Enzyme Inhibition and Medicinal Chemistry, vol. 31, no. 2, pp. 266–275, 2016. View at: Publisher Site | Google Scholar
  54. S. A. Holanda Pinto, L. M. S. Pinto, G. M. A. Cunha, M. H. Chaves, F. A. Santos, and V. S. Rao, “Anti-inflammatory effect of α, β-Amyrin, a pentacyclic triterpene from Protium heptaphyllum in rat model of acute periodontitis,” Inflammopharmacology, vol. 16, no. 1, pp. 48–52, 2008. View at: Publisher Site | Google Scholar
  55. F. A. Oliveira, R. C. P. Lima-Junior, W. M. Cordeiro et al., “Pentacyclic triterpenoids, α,β-amyrins, suppress the scratching behavior in a mouse model of pruritus,” Pharmacology Biochemistry and Behavior, vol. 78, no. 4, pp. 719–725, 2004. View at: Publisher Site | Google Scholar
  56. A. Genovese, A. Detoraki, F. Granata, M. R. Galdiero, G. Spadaro, and G. Marone, “Angiogenesis, lymphangiogenesis and atopic dermatitis,” New Trends in Allergy and Atopic Eczema, vol. 96, pp. 50–60, 2012. View at: Publisher Site | Google Scholar
  57. A. Sääf, A. Pivarcsi, M. C. G. Winge et al., “Characterization of EGFR and ErbB2 expression in atopic dermatitis patients,” Archives of Dermatological Research, vol. 304, no. 10, pp. 773–780, 2012. View at: Publisher Site | Google Scholar
  58. K.-H. Hsieh, C.-C. Chou, and S.-F. Huang, “Interleukin 2 therapy in severe atopic dermatitis,” Journal of Clinical Immunology, vol. 11, no. 1, pp. 22–28, 1991. View at: Publisher Site | Google Scholar
  59. S. Ramalingam, M. Packirisamy, M. Karuppiah et al., “Effect of β-sitosterol on glucose homeostasis by sensitization of insulin resistance via enhanced protein expression of PPRγ and glucose transporter 4 in high fat diet and streptozotocin-induced diabetic rats,” Cytotechnology, vol. 72, no. 3, pp. 357–366, 2020. View at: Publisher Site | Google Scholar
  60. J. Dai, M.-K. Choo, J. M. Park, and D. E. Fisher, “Topical ROR inverse agonists suppress inflammation in mouse models of atopic dermatitis and acute irritant dermatitis,” Journal of Investigative Dermatology, vol. 137, no. 12, pp. 2523–2531, 2017. View at: Publisher Site | Google Scholar
  61. M. I. Dawson and Z. Xia, “The retinoid X receptors and their ligands,” Biochimica et Biophysica Acta (BBA)—Molecular and Cell Biology of Lipids, vol. 1821, no. 1, pp. 21–56, 2012. View at: Publisher Site | Google Scholar
  62. J. Mihály, J. Gericke, G. Aydemir et al., “Reduced retinoid signaling in the skin after systemic retinoid-X receptor ligand treatment in mice with potential relevance for skin disorders,” Dermatology, vol. 225, no. 4, pp. 304–311, 2012. View at: Publisher Site | Google Scholar
  63. N. Malhotra, J. M. Leyva-Castillo, U. Jadhav et al., “RORα-expressing T regulatory cells restrain allergic skin inflammation,” Science Immunology, vol. 3, no. 21, Article ID eaao6923, 2018. View at: Publisher Site | Google Scholar
  64. M. Mboge, B. Mahon, R. McKenna, and S. Frost, “Carbonic anhydrases: role in pH control and cancer,” Metabolites, vol. 8, no. 1, p. 19, 2018. View at: Publisher Site | Google Scholar
  65. C. T. Supuran, A. Scozzafava, and A. Casini, “Carbonic anhydrase inhibitors,” Medicinal Research Reviews, vol. 23, no. 2, pp. 146–189, 2003. View at: Publisher Site | Google Scholar
  66. M. Kamsteeg, P. L. J. M. Zeeuwen, G. J. de Jongh et al., “Increased expression of carbonic anhydrase II (CA II) in lesional skin of atopic dermatitis: regulation by Th2 cytokines,” Journal of Investigative Dermatology, vol. 127, no. 7, pp. 1786–1789, 2007. View at: Publisher Site | Google Scholar
  67. J.-P. Hachem, D. Crumrine, J. Fluhr, B. E. Brown, K. R. Feingold, and P. M. Elias, “pH directly regulates epidermal permeability barrier homeostasis, and stratum corneum integrity/cohesion,” Journal of Investigative Dermatology, vol. 121, no. 2, pp. 345–353, 2003. View at: Publisher Site | Google Scholar
  68. J. Schalkwijk, M. Kamsteeg, G. J. D. Jongh, and P. L. J. M. Zeeuwen, New Use of Carbonic Anhydrase Inhibitors, University of Saskatchewan, Nijmegen, Netherlands, 2007.
  69. P. Sertznig and J. Reichrath, “Peroxisome proliferator-activated receptors (PPARs) in dermatology,” Dermato-Endocrinology, vol. 3, no. 3, pp. 130–135, 2011. View at: Publisher Site | Google Scholar
  70. Y. Benjamini and Y. Hochberg, “Controlling the false discovery rate: a practical and powerful approach to multiple testing,” Journal of the Royal Statistical Society: Series B (Methodological), vol. 57, no. 1, pp. 289–300, 1995. View at: Publisher Site | Google Scholar
  71. J. Mihály, A. Gamlieli, M. Worm, and R. Rühl, “Decreased retinoid concentration and retinoid signalling pathways in human atopic dermatitis,” Experimental Dermatology, vol. 20, no. 4, pp. 326–330, 2011. View at: Publisher Site | Google Scholar
  72. L. I. McKay and J. A. Cidlowski, “Molecular control of immune/inflammatory responses: interactions between nuclear factor-κb and steroid receptor-signaling pathways,” Endocrine Reviews, vol. 20, no. 4, pp. 435–459, 1999. View at: Publisher Site | Google Scholar
  73. K. Eyerich and N. Novak, “Immunology of atopic eczema: overcoming the Th1/Th2 paradigm,” Allergy, vol. 68, no. 8, pp. 974–982, 2013. View at: Publisher Site | Google Scholar
  74. S. Noda, M. Suárez-Fariñas, B. Ungar et al., “The Asian atopic dermatitis phenotype combines features of atopic dermatitis and psoriasis with increased TH17 polarization,” Journal of Allergy and Clinical Immunology, vol. 136, no. 5, pp. 1254–1264, 2015. View at: Publisher Site | Google Scholar
  75. P. M. Brunner, A. Israel, N. Zhang et al., “Early-onset pediatric atopic dermatitis is characterized by TH2/TH17/TH22-centered inflammation and lipid alterations,” Journal of Allergy and Clinical Immunology, vol. 141, no. 6, pp. 2094–2106, 2018. View at: Publisher Site | Google Scholar
  76. P. Rerknimitr, A. Otsuka, C. Nakashima, and K. Kabashima, “The etiopathogenesis of atopic dermatitis: barrier disruption, immunological derangement, and pruritus,” Inflammation and Regeneration, vol. 37, no. 1, p. 14, 2017. View at: Publisher Site | Google Scholar
  77. C. S. Williams, M. Mann, and R. N. DuBois, “The role of cyclooxygenases in inflammation, cancer, and development,” Oncogene, vol. 18, no. 55, pp. 7908–7916, 1999. View at: Publisher Site | Google Scholar
  78. J. Lee, H. Jung, K. Han et al., “Association between metabolic syndrome and atopic dermatitis in Korean adults,” Acta Dermato Venereologica, vol. 97, no. 1, pp. 77–80, 2017. View at: Publisher Site | Google Scholar
  79. E. Reisin and M. A. Alpert, “Definition of the metabolic syndrome: current proposals and controversies,” The American Journal of the Medical Sciences, vol. 330, no. 6, pp. 269–272, 2005. View at: Publisher Site | Google Scholar
  80. Z. Ali, C. S. Ulrik, T. Agner, and S. F. Thomsen, “Association between atopic dermatitis and the metabolic syndrome: a systematic review,” Dermatology, vol. 234, no. 3-4, pp. 79–85, 2018. View at: Publisher Site | Google Scholar
  81. D. Q. Xiao, P. H. Xu, H. Deng et al., “The characterization and clinical significance of PI3K/Akt signaling pathway activation in the peripheral T cells of pediatric patients with atopic dermatitis,” International Journal of Clinical and Experimental Medicine, vol. 10, pp. 2904–2910, 2017. View at: Google Scholar
  82. B. Hartmann, R. Riedel, K. Jörß et al., “Vitamin D receptor activation improves allergen-triggered eczema in mice,” Journal of Investigative Dermatology, vol. 132, no. 2, pp. 330–336, 2012. View at: Publisher Site | Google Scholar
  83. H. A. Ahmed and I. Y. Alkali, “In silico molecular docking studies of some phytochemicals against peroxisome-proliferator activated receptor gamma (PPAR-γ),” GSC Biological and Pharmaceutical Sciences, vol. 5, no. 2, pp. 1–5, 2018. View at: Publisher Site | Google Scholar
  84. L. Alappat, M. Valerio, and A. B. Awad, “Effect of vitamin D and β-sitosterol on immune function of macrophages,” International Immunopharmacology, vol. 10, no. 11, pp. 1390–1396, 2010. View at: Publisher Site | Google Scholar
  85. M. Megna, M. Napolitano, C. Patruno et al., “Systemic treatment of adult atopic dermatitis: a review,” Dermatology and Therapy, vol. 7, no. 1, pp. 1–23, 2017. View at: Publisher Site | Google Scholar
  86. X. Zhou, S. W. Seto, D. Chang et al., “Synergistic effects of Chinese herbal medicine: a comprehensive review of methodology and current research,” Frontiers in Pharmacology, vol. 7, p. 201, 2016. View at: Publisher Site | Google Scholar
  87. J. Dong, “The relationship between traditional Chinese medicine and modern medicine,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 153148, 10 pages, 2013. View at: Publisher Site | Google Scholar

Copyright © 2021 Zhen-Zhen 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.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views450
Downloads507
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.