key: cord-0016920-rvblg2hc
authors: Zhang, Yunsen; Zhao, Zikuang; Chen, Huimin; Fu, Yutong; Wang, Wenxiang; Li, Qi; Li, Xuanhao; Wang, Xiaobo; Fan, Gang; Zhang, Yi
title: The Underlying Molecular Mechanisms Involved in Traditional Chinese Medicine Smilax china L. for the Treatment of Pelvic Inflammatory Disease
date: 2021-04-08
journal: Evid Based Complement Alternat Med
DOI: 10.1155/2021/5552532
sha: 8abd071dda99caf6178eb14d9edd4f5c4f0fbc86
doc_id: 16920
cord_uid: rvblg2hc

Smilax china L. (SCL) is extensively used in the treatment of pelvic inflammatory disease (PID). This study aimed to clarify the potential active ingredients of SCL and mechanisms on PID. SCL was widely distributed in Japan, South Korea, and China, which was traditionally considered heat-clearing, detoxicating, and dampness-eliminating medicine. Systems pharmacology revealed that 32 compounds in SCL may interact with 19 targets for immunoenhancement, antiapoptosis, anti-inflammation, and antioxidant activity of the PID model. Molecular docking revealed that isorhamnetin, moracin M, rutin, and oxyresveratrol may have higher binding potential with prostaglandin-endoperoxide synthase 2 (PTGS2), mitogen-activated protein kinase 1 (MAPK1), siderocalin (LCN2), tumor necrosis factor (TNF), and matrix metalloprotein-9 (MMP9), respectively. Molecular dynamics simulation showed that the binding modes of moracin M-MAPK1, rutin-TNF, and oxyresveratrol-MMP9 complexes were more stable, evidenced by relatively smaller fluctuations in root mean square deviation values. Conclusively, SCL may treat PID by inhibiting inflammatory factors, antitissue fibrosis, and microbial growth.

Pelvic inflammatory disease (PID), the infection and inflammation of the female upper genital tract, is a common cause of infertility, chronic pain, and ectopic pregnancy [1] . Diagnosis and management are challenging, largely due to a polymicrobial etiology that is not fully delineated [2] . Reportedly, it is estimated that 2.5 million American women aged between 18 and 44 have received a PID diagnosis in their lifetime [3] , and one in eight women with a history of PID encountered difficulties in getting pregnant [4] . PID treatment is mainly based on broad-spectrum antibiotic regimens and surgical treatment [2, 5] . Antibiotics are effective in lessening short-term morbidity but have no effects on long-term complications, due to the disease's complex mechanism and long-term process [6] . Although the incidence of PID has decreased because of screening for gonorrhea and chlamydia and the early intervention of broad-spectrum antibiotics, damage to the reproductive system caused by infection has not been ameliorated [7] . erefore, the therapeutic goal for the treatment of PID ought to include both short-term microbiological effects and long-term prevention of sequelae [6] . Also, the use of antibiotics is limited by the emergence of antibiotic resistance and PID without an identified pathogen. To inhibit progress, alleviate the long-term sequelae of PID, and avoid antibiotic resistance, it is often used in conjunction with traditional Chinese medicine (TCM) [8, 9] . Recent studies showed that the anti-inflammatory and immune mechanisms of PID were related to T cells, B cells, IG, cytokines (e.g., IL-6, TNF-α, and IL-1β), prostaglandin-endoperoxide synthase 2 (PTGS2), matrix metalloprotein-9 (MMP9), and TLRs signaling pathways [8, 10] .

From the theoretical perspective of TCM, the internal pathogenesis of PID is the disharmony of yin-yang (a general term for all kinds of pathological changes due to imbalance and incoordination of yin and yang) and insufficiency of healthy qi, and the external is that the dampness-heat and heat toxin invaded the thoroughfare and conception vessels, uterus, and uterine vessels. It caused blood blockage and abdominal pain ( Figure 1 ). Hence, the core excitation of PID onset is the blood stasis induced by dampness-heat (dampness-heat syndrome) [11, 12] . Smilax china L. (SCL), a Liliaceae plant, commonly known as "Baqia" (or "Jin Gang Teng"), is widely distributed in Asia ( Figure 2 ). From 1887 to 2020, it was found 7636 times ( Figure 2 (a)) according to GBIF online database (https:// www.gbif.org/species/5295472/metrics) and mainly distributed in Japan (3630 times) (Figure 2 Up to now, SCL has been included in Chinese pharmacopeia [13] with the effects of heat-clearing, detoxicating, and dampness-eliminating. It has still been widely used in TCM for the treatment of PID and formulated into granules, syrup, pills, and capsules, demonstrating a good curative effect [9, 14] . Based on the previous studies, steroid saponin, flavonoids, glycosides, and stilbenes are the principal chemical compounds in SCL and demonstrated anti-inflammatory effect via TLR-4-mediated signaling pathway [14] [15] [16] .

e flavonoid derivatives such as engeletin, isorhamnetin, and quercetin are the main constituents for the treatment of PID by inhibiting extracellular regulatory protein kinase and SMAD2/3 protein phosphorylation, thereby relieving the degree of fibrosis in the uterus via ERK1/2 and TGFβ-SMAD2/3 signaling pathways [17] . However, most of these studies showed the rough mechanisms of extract in SCL against PID, and the main active compounds and how these compounds interact with PIDrelated targets to interfere with relevant signaling pathways are still unclear.

TCM is designed to maintain the balance of the body's functions utilizing a lot of intricate compounds in herbs. Because multiple constituents may produce synergistic regulation on different targets, elucidating the mechanisms of TCM always takes lots of time and resources. ere is no doubt that systems pharmacology has recently emerged as a new field including physiology, genetics, biochemistry, and molecular simulation via integrating various research methods to investigate the complicated mechanisms of multiple compounds [18, 19] . Nowadays, systems pharmacology has been applied for revealing the pharmacological mechanism of TCM from the perspective of entirety. For instance, most of the ingredients from wellresearched herbs are carried out by molecular simulation, such as pharmacophore matching and inverse-docking to clarify the candidate targets which is available for researchers to further illustrate the integral mechanism of TCM [20] [21] [22] . Otherwise, the particle interaction, a role of the fundament of integral regulation, is the same concerned. us, static molecular docking and calculation of molecular mechanics-generalized Born surface area (MM-GBSA) free binding energy will provide a view of good binding pose and binding free energy to ensure that the complexes of compounds and targets possess enough energy to engender reaction of biochemistry [23, 24] . Finally, molecular dynamics (MD) simulation is utilized to evaluate the stability of protein-ligand complexes obtained from molecular docking using root mean square deviation (RMSD) and explore the noncovalent interaction between the active ingredients of herbs and the predicted targets, such as hydrogen bond (H-bond) and decomposition of molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) energy of amino acid residues [25] . Overall, a schematic representation of the workflow in this study is shown in Figure 4 , which makes it possible to systematically decode the active compounds and mechanism of TCM in the network. Hence, we evaluated the whole candidate targets of active compounds and provided a perspective of the integral mechanism via enriching the functions of targets and dissected the molecular mechanism of SCL in the treatment of PID using computational systems pharmacology.

SCL compounds were systematically listed as ligands from published paper mining [26, 27] , TCMID (http://www. megabionet.org/tcmid/) [28] , SymMap (https://www. symmap.org/detail/SMHB00008) [29] , and TCMSP databases (http://tcmspw.com/tcmsp.php) [30] . All compound structures from PubChem (https://pubchem.ncbi.nlm.nih. gov/) [31] were filtered by utilizing the "Lipinski rules" of the Molinspiration database (https://www.molinspiration.com/ cgi-bin/properties) [32] . In the field of drug discovery, the Lipinski rules were used to screen the compound database to eliminate molecules that were unsuitable for drug use, including n·OHNH ≤ 5, n·ON ≤ 10, MW ≤ 500, and miLogP ≤5. Compounds that met the Lipinski rules and others that did not but possessed good bioactivity were used in systems pharmacology and molecular docking [33, 34] . e 2D structures (.sdf format) of all compounds were generated by ChemBioOffice2014 [35] . Figure 1 : Effects of SCL on PID from view point of traditional Chinese medicine. From the theoretical perspective of TCM, the left and right circles represent the development process from physiological state of human (arrows of five phases are of the same size) to pathological state of PID patients (fire phase increased and water phase decreased). And the main interfering process of SCL against PID is between two circles. Above all, the internal pathogenesis of PID is the disharmony of yin-yang (a general term for all kinds of pathological changes due to imbalance and incoordination of yin and yang) and insufficiency of healthy qi (a collective designation for all normal functions of the human body and the abilities to maintain health, including the abilities of self-regulation and adaptation), and the external is that the dampness-heat and heat toxin invaded the thoroughfare and conception vessels, uterus, and uterine vessels. It caused blood blockage and abdominal pain. Hence, the core excitation of PID onset is the blood stasis induced by dampness-heat.

(a) J a p a n : 3 6 3 0 Evidence-Based Complementary and Alternative Medicine

Compounds that satisfied Lipinski rules were uploaded to SwissTargetPrediction (http://www.swisstargetprediction.ch/) [36] , PharmMapper (http://www.lilab-ecust.cn/pharmmapp er/) [37] , and SEA (http://sea.bkslab.org/) [38] to obtain preliminarily candidate targets. All the genes should be from "Homo sapiens" to clarify the function of critical targets and were proofread by the UniProt database (https://www. uniprot.org/) [39] .

Systems pharmacology was applied to analyze the interaction between SCL and PID and the selection of critical targets. To identify the intersection of targets between PID and SCL, the title of "pelvic inflammatory disease" was placed in GeneCards (https://www. genecards.org/) [40] , DisGeNET (https://www.disgenet. org/) [41] , DrugBank (https://www.drugbank.ca/) [42] , and existing research [43, 44] to obtain gene names of PID targets, which manually confirmed that each target had a clinical study in PID. Additionally, the intersection genes (compound-and disease-related targets) were used to perform annotation analysis of the obtained crossover genes by using Gene Ontology (GO) and the KEGG pathway analysis functions in the STRING platform (https://string-db.org/) [45] , with the intersection genes directly mapped to the pathway. Cytoscape [46] was used to visualize a network of "Compound-target-pathway (C-T-P network)". In this network, each compound, target, and pathway was indicated by nodes, and the interactions between each node were described by edges. e network was established to project an overview of the interactions among compounds, targets, and pathways.

To reveal the direct and indirect roles active compounds of SCL played in the pelvic inflammatory targets, the intersection targets were introduced to the STRING platform and a graphical network of PPI was generated [47] . In the network, each node represents all the proteins produced by a single, protein-coding gene locus, and edges represent proteinprotein associations which are meant to be specific and meaningful, i.e., proteins jointly contribute to a shared function; this does not necessarily mean they are physically binding each other. Ultimately, to clarify the interfering mechanism, compounds broke into PPI network targets with a high degree and relating to the critical pathway in the PPI network were selected to operate molecular docking. [48] , and the crucial targets were docked with active compounds. Schrödinger Glide was used to pretreat the 3D protein structure for docking, including adjusting the bond orders to ensure the stability of the chemical bonds between the atoms, adding the missing hydrogen atoms and amino acid residues, optimizing the orientation of amino acids and hydrogen atoms, optimizing the distribution of H-bonds, and removing water molecules and heterogeneous molecules. Finally, energy minimization with a force field OPLS-2005 was supplied.

At the end of pretreatment, the Receptor Grid Generation module was employed to select the ligand-binding cavity and generate a grid in protein. Next, the active compounds were introduced to Maestro and optimized by the liquid simulation of OPLS-2005 all-atoms force field in the LigPrep module, as well as combined into a ligand package. e docking accuracy was evaluated by standard precision (SP), as well as flexible docking. At this point, the docking preparation was completed, the scaling factor and partial charge cutoff of van der Waals radius scaling 1.0 and 0.25 were used to generate the grids on active sites, and ligand package was selected to perform molecular docking in the Ligand Docking module.

Next, the generated Glide G-score was used to assess the affinity between compounds and proteins. Moreover, crucial targets were docked with their self-ligands to set positive contrast, and their Glide G-scores were used to measure whether compounds possessed a good affinity to the protein and standardize the score of compounds to be visualized as a heat map by MeV [49] (step 1: normalized genes/rows; step 2: hierarchical clustering: average dot product-complete linkage clustering).

Docking. Good poses and good scores were obvious by SP docking, but what was the binding free energy of the docked complex was another problem. Docking results showed the active compounds did bind to the active site of the protein, but could this association last long enough to elicit any potential biological response, as biological response largely depends upon the binding free energy of the association. erefore, the docked complexes in SP mode were subjected to binding free energy calculation (MM-GBSA) using the Prime module of Maestro [50] . A total of 14 active compounds were selected for this analysis. Evidence-Based Complementary and Alternative Medicine 5 obtained by molecular docking. Otherwise, all atoms of the system in complexes are sanctioned to motion and interact for a fixed period (∼45 ns), and the trajectories of atoms and molecules are defined through Newton's equations of motion. In this study, the best poses of SP docked ligand/protein with an excellent binding free energy were employed to operate explicit solvent MD simulation. GROMACS [51] as a computational tool of MD simulation was used for this purpose.

Due to five complexes that would perform MD simulation, we had to protocol slightly various systems to reach our expectation as shown in Supplementary Table 1. In general, coordinates and charge of ligands were generated using PRODRG 2.5 [52] , and protein and SPC water model were described by the GROMOS96 43A1 force field [53] and defined in a solvent box with 8.5 nm × 8.5 nm × 8.5 nm. Na + and Clions were added to ensure the overall neutrality of the systems. Each MD simulation system was first relaxed to remove possible steric crashes by the steepest descent energy minimization algorithm and stopped minimization when the maximum force is <100.0 KJ/mol. In the second step, a 100 ps simulation was performed utilizing the canonical ensemble (NVT ensemble) using the modified Berendsen thermostat with a slowly ascending temperature from 0 K to 300 K, a fast temperature relaxation constant of 0.1 ps, and a temperature coupling of protein and ligand to prevent system bursting. Next, the periodic boundary condition was employed to produce the constant temperature and pressure (NPT) ensembles. e pressure was set at 1.0 bar and was controlled by the isotropic pressure scaling protocol applied in GROMACS. Moreover, no cutoff limit was used for electrostatic forces by employing the particle mesh Ewald (PME) algorithm. All bonds were constrained using the LINCS algorithm. en the simulation time for each system was 45 ns, and the trajectories of simulated systems were saved every 10 ps. Finally, disintegrated calculation of binding free energy was hired to qualitatively analyze the principal force of interaction between protein and ligand. All results were visualized by QtGrace (https://sourceforge.net/projects/ qtgrace/), VMD (http://www.ks.uiuc.edu/Research/vmd/), and LigPlot + (https://www.ebi.ac.uk/thornton-srv/ software/LigPlus/).

After filtering was operated for Lipinski rules, thirty-two in sixty-eight compounds were screened as active compounds in SCL and are listed in Supplementary Table 2 . Chemical formats of smile and SDF were generated to predict the potential targets. en 718 potential targets (some were duplicates) associated with compounds were screened from the Swis-sTargetPrediction, PharmMapper, and SEA server, respectively (Supplementary Table 3 ). Otherwise, eighty-six PID-related targets were retrieved from databases (Supplementary Table 4 ). en, to clarify the relationship between herb and disease, a network of compounds, targets, and pathways was visualized by Cytoscape. We found that a total of 32 compounds could act on 19 key targets and associate with 16 relevant GO annotations and 10 effective pathways ( Figure 5 ).

Analysis. 19 intersection targets were analyzed using the PPI network in the STRING platform ( Figure 6 ). Relevant parameters of the network were as follows: (1) Otherwise, red nodes (MMP9, tumor necrosis factor (TNF), interleukin-6 (IL-6), PTGS2, neutrophil gelatinase-associated lipocalin (LCN2), mitogen-activated protein kinase 1 (MAPK1), interleukin-2 (IL-2), and signal transducer and activator of transcription 3 (STAT3)) were used to highlight the IL-17 signaling pathway and 17 cell differentiation. All highlighted targets except IL-2, STAT3, and IL-6 were selected as core targets to illuminate the SCL interfering mechanism against PID utilizing molecular docking and MD simulation.

Energy. e heat map was employed to stick out the features of 32 active compounds as shown in Figure 7 (original data are presented in Supplementary Table 5 ). Compounds that had a high activity clustered together excellently, with a high affinity to PTGS2, LCN2, TNF, MAPK1, and MMP9. Rutin (10), isorhamnetin (18) , oxyresveratrol (30) , and moracin M (44) were found to occupy the top score, which exceeded or neared the original ligands in verified docking. Furthermore, compounds demonstrated a binding affinity to one or several targets. According to the distance metric of average dot product in Mev, 14 capital protein-ligand molecular interactions were analyzed in Table 1 . Dihydrokaempferide (62) Gramine (60) Dihydrokaempferol (12) Moracin M (44) Dihydroquercetin (48) Piceatannol (21) Kaempferide (58) Resveratrol (5) Oxyresveratrol (30) Polydatin (4) Rutin (10) Quercetin (17) Isorhamnetin (18) Kaempferol (11) Coumarin (55) Maackoline (43) Astilbin (3) Engeletin (1) Isoengeletin (2) Trans-2-hydroxycinnamic acid (46) Diosgenin (63) Sarsasapogenin (8) Butyl β-D-fructopyranoside (38) 3,5-Di-4-β-D-glu acid (42) Protocatechuic acid (40) Caffeic acid (51) Methylprotodioscin_qt (53) Oleanolic acid (50) Gallic acid (45) Vanillic acid (41) Syringic acid (39) Gentisic acid (47) 0.89265305 0.738021 0.0 LCN2 PTGS2 MAPK1 TNF MMP9 Figure 7 : Clustering heat map between compounds and PID-related targets. All compounds with a docking score of proteins were divided into 3 clusters from left to right considering the affinity. 14 active compounds on the right were operated MM-GBSA binding free energy. "Di": dimethoxy; "glu": glucopyranosyl cinnamic. [59] . In this study, five complexes were used to calculate the RMSD within 45 ns. Firstly, RMSD analysis and complexes in solvent depicted that PTGS2 (Figures 9(a) and 9(b)) started to stabilize after 10 ns, and it maintained 20 ns stability until 30 ns and slightly increased after 30 ns. Interestingly, a similar fluctuation was captured in the change of total H-bonds at 30 ns (Supplementary Figure 1A) , and the average number of H-bonds before 25 ns was 2.7 but increased to 3.5 during the last 15 ns. is fluctuation indicated that the stable growth of the H-bond and the transition of one stabilized configure to another. e ligand (Figure 9 (b)) started stabilizing after 10 ns and maintained to 45 ns. e average RMSD of protein and ligand was 0.312 and 0.082 nm, respectively.

Secondly, MAPK1-moracin M complex (Figures 9(c) and 9(d)) gained stability at around 10 ns, and the RMSD pattern of ligand suddenly increased at 27 ns and maintained to 45 ns. At the same time, the corresponding change appeared in the number of total H-bonds (Supplementary Figure 1B) . e average number of H-bonds before 27 ns was 3.7 and increased to 4.1 during the last 18 ns. e average RMSD of protein and ligand was 0.286 and 0.116 nm, respectively.

irdly, for the LCN2-rutin system (Figures 9(e) and 9(f )), stable protein-ligand interaction was also observed. e protein and rutin were both equilibrated at ∼11 ns with a slight fluctuation. e main fluctuation of protein was observed at 27∼30 ns and equilibrated during the last 15 ns. In the view of total H-bond number (Supplementary Figure 1C) , the H-bonds of complex initially increased before ∼22 ns (average number: 2.8), reduced by 0.7 at 22∼32 ns, and stabilized at 3.2 during the last 13 ns. e average RMSD of LCN2 and rutin was 0.247 and 0.164 nm, respectively. 

Fourthly, for multiple chain (four chains) system of TNF-rutin (Figures 9(g) and 9 (h)), it was time-consuming to equilibrate. erefore, a stable protein-ligand interaction was observed during the last 15 ns. e RMSD of TNF constantly rose before ∼26 ns. e ligand initially stabilized at ∼16 ns with a slight fluctuation and ultimately stabilized at 30 ns, and a sudden increase in the RMSD of ligand lasted 3 ns was captured. e total H-bond number was selected for analysis during the last 20 ns. A constant change and a large span of the H-bond number were recorded before 30 ns, which indicated an unstable H-bonding. en, the span of H-bond variation gradually reduced, and 2∼4 stable H-bonds were retained to the end. e average RMSD of TNF and rutin was 0.389 and 0.189 nm, respectively.

Finally, for the MMP9-oxyresveratrol system (Figures 9(i) and 9(j)), the RMSD analysis depicted that the protein achieved stability at around 7 ns and fluctuated slightly at ∼16 ns, and the ligand equilibrated at 1 ns and maintained until 40 ns with a slight fluctuation at ∼ 17 ns. Meanwhile, the change of the H-bond number gradually stabilized in the last 25 ns. e average RMSD of MMP9 and oxyresveratrol was 0.282 and 0.097 nm, respectively. Overall, the RMSD and total H-bond change analysis of five MD simulation complexes illustrated that the all above systems maintained a period of stability with ligands.

Hydrogen bonding is among the most essential parameters to understand the binding affinity of small molecules towards a biomacromolecule (e.g., protein). A large number of H-bonds present in between protein and small molecules signify a strong binding affinity. In this regard, H-bonds between the natural molecules and disease-related proteins were monitored over the MD simulation time. 

In the TNF-rutin complex, 25 different H-bonds were detected. On average, we observed 2.2 H-bonds, and the distance and angle were 2.94Å and 16.93°, respectively. Due to a constant change and large span of the H-bond number mentioned in "Section 3.4.1", the occupancy of amino acid residues including H-bond receptor (Val 123 main and Tyr 119 side) and donor (Ser 60 main) was generally lower.

As for the MMP9-oxyresveratrol complex, 20 different H-bonds were detected. ere were on average 4.3 H-bonds, and the relevant distance and angle were 2.78Å and 14.78°, respectively. ree H-bond receptors (e.g., Leu 104 main, His 230 main, and Pro 102 main) and one donor (Leu 234 main) were bound to oxyresveratrol stably. ese results hinted that the natural molecules interacted effectively towards the active site of PID-related proteins with a significant property of H-bonds.

To explore the interaction between proteins and their ligands, energy decomposition of each complex was performed by using g_mmpbsa [60] module of GROMACS software. A total of 100 snapshots were extracted from the stable and continuous trajectories for the free energy calculation. e binding free energy could be divided into van der Waals Interaction (∆Evdw), electrostatic energy (∆Eele), polar solvation interaction (∆Epol), and solvent-accessible surface area (∆Esasa). e results were listed in Table 4 and Figure 10 . To understand the energy contribution, ∆Bind (−146.535 ± 1.934 KJ/mol) of the PTGS2-isorhamnetin complex was decomposed to each amino acid residue. Figure 10 (a) hinted that the important residues including Val 523, Leu 352, Phe 518, and His 90 involved in the binding site of PTGS2 (Figure 8(a) ) showed a positive tendency for binding. As Figure 10(b) shows, Lys 145, Glu 62, and Asp 158 in MAPK1 showed a negative tendency to bind, and the latter two residues were within 4Å of ligand, but Glu 62 was the important residue which formed excellent H-bond interaction with the ligand. Otherwise, residues (e.g., Val 30 and Ile 75) that showed positive binding tendency were within 4Å to the ligand. e results showed that residues which had a positive binding tendency were closer to ligands (within 4Å) and even formed H-bonds and pi-pi stacking interaction, and the negative residues were apt to far away from the ligand as shown in Figures 10(c) -10(e). All the above results were consistent with our previous docking results.

Traditionally, in modern medicine, drugs were designed to target specific proteins relevant to the disease. However, herbs could possess even hundreds of compounds with multiple targets in TCM, which presented a significant obstacle to the exploration of drug mechanisms [61] . Fortunately, systems pharmacology, molecular docking, and MD simulation provided a holistic perspective to clarify the potential mechanisms for illustrating the integral regulation-based active compounds of TCM and conducting the consequent experiment [19] . erefore, the above methods were employed to reveal the intricate mechanisms of SCL in the treatment of PID.

According to the C-T-P network ( Figure 5 ), we proposed a simple inference as shown in Figure 11 (details are presented in Supplementary Figure 2) , with the IL-17 signaling pathways (the highest value of −log10 (FDR)) considered potentially efficacious utilizing in-depth excavation of the above network. IL-17 signaling pathway reportedly had dual regulatory roles in proinflammatory and host defense processes [62, 63] . On the one hand, excessive secretion of IL-17A and IL-17F from 17 cells can induce massive inflammatory factors including IL-6, IL-1β, and TNF-α. IL-17 signaling pathway can also synthesize prostaglandin E2 (PGE2) by inducing PTGS2, and the vasodilator effect of PGE2 also promoted inflammatory cells to enter the site of inflammation, so the activation of IL-17 signaling pathway had a strong proinflammatory effect [64] . In previous studies, SCL downregulated the expression of IL-6, interleukin-1 beta (IL-1β), TNF-α, IL-2, and PTGS2 in rats of the PID model, but its regulatory mechanisms involved remain unclear [17, 26] . However, combined with the results of this study, this may be the SCL inhibition of the IL-17 signaling pathway.

On the other hand, in the process of immunoregulation, the regulation of IL-17 signaling pathway can recruit neutrophil to the inflammatory region, releasing myeloperoxidase (MPO) [65] and inducing gene expression of LCN2 and matrix metalloproteinases (MMPs) [63, 66] . ese proteins make important impacts in host defense, which may contribute to alleviating PID symptoms by enhancing the immune function of the body. For instance, MMP9 was activated to promote embryo formation, wound healing, and transfer of inflammatory cells [66] [67] [68] . Interestingly, it was reported that SCL upregulated the expression of MMP-2 and MMP-9 in rats of the PID model and downregulated the MMPs inhibitor TIMP-1, thereby restoring the balance between MMPs and TIMP-1 and reducing tissue fibrosis during PID [17] . Collectively, these results demonstrated the multitarget regulation of compounds in SCL.

Although the above network provided a clear view of the integral regulation of SCL, there was still a barrier to validate the facticity of each research data and it was difficult to consider the critical effect of each result. Hence, how to select a crucial result (targets or pathways) to concern and validate was an issue worth pondering carefully for researchers. erefore, PPI analysis could play an essential role which helps us find the key targets and pathways and explore critical radioactive targets in the network. In this study, IL-17 signaling pathway and 17 cell differentiation-related targets were considered focuses in the PPI. As previously reported, PTGS2 is inducible and usually produces inflammatory prostaglandins, which mediate responses to physiological stress (infection and inflammation), stimulate chronic inflammation, and is a target for nonsteroidal anti-inflammatory drugs (NSAIDs) [56] . As shown in Table 1 and Figure 8 , the binding capacity of isorhamnetin to PTGS2 (Glide G-score: −9.757) was very close to rofecoxib (−9.800), suggesting that isorhamnetin may be a potential novel COX-2 inhibitor somewhat analogous to rofecoxib. Additionally, all the active compounds ( Figure 7 ) mostly belonged to flavonoids and stilbenes. Based on existing pharmacodynamic investigations, these flavonoids and stilbenes have achieved obvious anti-inflammatory effects and ameliorated fibrosis in PID animal models by inhibiting the synthesis or release of histamine, 5-hydroxytryptamine (5-HT), and PGE2, as well as enhancing the production of MMP9 in uteri [17, 69] . Particularly, engeletin (1), polydatin (4), and resveratrol (5) have inhibited the release of IL-6 and TNF-α [70] . Rutin (10) inhibited the release of IL-2 and TNF-α [26] .

ese results and existing experiments indicated that active ingredients in SCL, including isorhamnetin (18) , polydatin (4), oxyresveratrol (30) , and piceatannol (21) , could inhibit the activity of PTGS2 to decrease the synthesis of PGE2 [71] , which contributes to the restoration of PID by inflammatory inhibition. What is more, MD simulation of five complexes Table 1 ) whatever single chain (LCN2 and MMP9) and multiple chains (PTGS2, MAPK1, and TNF) could stabilize ultimately ( Figure 9 ). e fluctuation of various systems was significantly relevant to H-bonding change. As shown in Table 2 , the properties of H-bond in different systems were decoded. For the MAPK1-moracin M and MMP9-oxyresveratrol systems, the average H-bond numbers were greater than the LCN2rutin and TNF-rutin systems, but the detected H-bonds were lower than them. e structures of rutinose and flavonoid aglycones on rutin were easy to form H-binding interaction to surrounding amino acid residues in the protein. Otherwise, the H-bonding interaction and energy contribution of crucial residues were analyzed (Table 3 and Figure 10 ); the significant H-bond-related residues including Phe 518 on PTGS2, Glu 62 on MAPK1, Ser 68 on LCN2, Val 123 on TNF, and Leu 104 on MMP9 were of great importance to maintain the system's stability. e significant residues of energy contribution (e.g., Val 523 and Leu 352 on PTGS2, Ile 75 on MAPK1, Trp 79 on LCN2, Leu 55 on TNF, and Phe 110 on MMP9) indicated the van der Waals interaction and electrostatic energy were the same important. So far, a part of the results in network pharmacology has been validated and discussed via existing experiments, molecular docking, and MD simulation, which can yet be regarded as an effective method to elucidate the active compound of SCL and relevant mechanisms against PID. In our study, the protein LCN2 which can inhibit microorganisms by chelating the iron ions [72] was found to have an important impact and phenolic compounds (e.g., moracin M) were also considered as key compounds in treating PID.

Mechanistically, compared to previous research studies on the main constituents for the treatment of PID by relieving the degree of fibrosis in the uterus via ERK1/2 and TGFβ-SMAD2/ 3 signaling pathways, we reported the critical active compounds and relevant binding modes of IL-17 signaling pathway and 17 cell differentiation-related targets in the treatment of PID by inhibiting inflammatory factors, antitissue fibrosis, and microbial growth. However, some tough problems associated with disease targets and traditional herbs still bother investigators. e first thing is that the vague and a small number of valuable targets based on current research studies cannot orient to the whole disease. erefore, systems pharmacology and molecular simulation are used to reveal the partly network mechanisms and molecular actions rather than the whole disease network mechanisms. And similar to the current drug discovery strategies, compared with the number of disease targets, the quality of the targets (druggability and crystal reliability) is paid more heed by researchers. Secondly, as for molecular docking, due to the limitation of computational power and force field algorithm, a newly and generally applicable force field to improve the accuracy of molecular docking and fast methods of the binding free energy calculation for virtual screening [73] are urgent to develop. Enough computational accuracy provides researchers enough confidence to conduct subsequent experiments. Finally, the specific components in TCM are needed to extract and identify completely to enrich the efficacious material basis. It is comparatively easy to illustrate the molecular mechanism once obtaining reliable targets and adequate ingredients. All above, there is a shortcut to uncover the overall network mechanism of TCM against diseases through cross cooperation among pharmacology, medicinal chemistry, and computational chemistry.

Considering the promising treatment potential of SCL on PID, efforts are in demand to reveal the acting targets for SCL and the unclear mechanisms behind the therapeutic potentials. In this study, the computational systems pharmacology method was applied to explore the active ingredients of SCL and provided an integral view of the mechanism against PID. e principal 32 potent ingredients for the treatment of PID were uncovered to regulate 718 candidate targets. Furthermore, in the PPI and C-T-P network analysis, 8 of 19 PIDrelated targets were mapped to the IL-17-signaling pathway and 17 cell differentiation. We focused on five reported PID-related targets PTGS2, MAPK1, LCN2, TNF, and MMP9. e interactions between active compounds and PIDrelated targets were described with static and dynamic evaluations. A total of 14 active compounds, including rutin (−40.46 kcal/mol), isorhamnetin (−51.06 kcal/mol), oxyresveratrol (−56.71 kcal/mol), and moracin M (−37.01 kcal/ mol), showed greater binding force to the therapeutic targets. At the same time, the amino acid residues in the hydrophobic cavity which played an important role in the process of complexes were revealed to guide the design of relevant drugs.

Overall, active ingredients of SCL exhibited a strong affinity to therapeutic targets of PID, thereby contributing to decreasing inflammation, ameliorating fibrosis, and inhibiting or eliminating microorganisms via bidirectional regulation of the IL-17 signaling pathway. However, it was a time-consuming and risky process to draw this kind of conclusion. Analysis of the network had to be up against a problem on how to select the principal results to focus, which is a big challenge but a core in systems pharmacology.

erefore, valid validation of results (e.g., static molecular docking, MD simulation, animals, and biochemistry) is equally important. rough the analysis of this study, IL-17 pathway was found to probably play a critical role in the development and treatment of PID, but relevant research was lacking and incomplete. Hence, serum, integrated pharmacodynamics, and pharmacokinetics will be utilized to clarify the components in serum and relevant therapeutic mechanisms of SCL on PID.

Compound-target-pathway ∆Eele:

Electrostatic energy ∆Epol:

Polar solvation interaction ∆Esasa:

Solvent-accessible surface area ∆Evdw:

van der Waals interaction GO:

Gene Ontology 

e data used to support the findings of the study are included within Supplementary Materials.

e authors declare that they have no conflicts of interest. 

Epidemiology, pathogenesis and treatment of pelvic inflammatory disease

Recommendations and rationale for the treatment of pelvic inflammatory disease

Susceptibility of endometrial isolates recovered from women with clinical pelvic inflammatory disease or histological endometritis to antimicrobial agents

Self-reported infertility and associated pelvic inflammatory disease among women of reproductive age-national health and nutrition examination survey

Systemic evaluation and Meta-analysis of efficacy and safety of Kangfu Xiaoyan Suppository in the treatment of pelvic inflammatory disease

Pelvic inflammatory disease

Pelvic inflammatory disease

Clinical efficacy of kangfuyan capsule and antibiotic therapy for chronic pelvic inflammatory disease and influence on levels of inflammatory factors

Guidelines for clinical application of Chinese patent medicine in the treatment of sequelae of pelvic inflammatory diseases

Study on the Mechanism of Fuyanshu Capsule by Regulating the TLRs/MyD88 Pathway through KLF4 and A Systematic Review of Expelling Dampness and Dissipating Blood Stasis erapy of Chinese Patent Medicine for SPID

Chinese medicine treatment of pelvic inflammatory disease

World Health Organization, Regional Office for the Western Pacific. WHO International Standard Terminologies on Traditional Medicine in the Western Pacific Region, WHO Regional Office for the Western Pacific

Chinese Pharmacopoeia Commission, Pharmacopoeia of the People's republic of china, China Press of Traditional Chinese Medicine

Tremendous health benefits and clinical aspects of Smilax china

Research progress on chemical component and pharmacological mechanism of Chinese medicine Smilax China L

e flavonoid-enriched extract from the root of Smilax china L. inhibits inflammatory responses via the TLR-4-mediated signaling pathway

Protection of flavonoids from Smilax china L. rhizome on phenol mucilage-induced pelvic inflammation in rats by attenuating inflammation and fibrosis

Computational systems pharmacology analysis of cannabidiol: a combination of chemogenomics-knowledgebase network analysis and integrated in silico modeling and simulation

Insight into ginkgo biloba L. Extract on the improved spatial learning and memory by chemogenomics knowledgebase, molecular docking, molecular dynamics simulation, and bioassay validations

Systems pharmacology approach to investigate the mechanism of Oryeong-san formula for the treatment of hypertension

Network pharmacology interpretation of fuzheng-jiedu decoction against colorectal cancer

A network pharmacology approach to reveal the underlying mechanisms of artemisia annua on the treatment of hepatocellular carcinoma

Exploring the inhibitory activity of valproic acid against the HDAC family using an MMGBSA approach

e MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities

Exploring the interaction between human focal adhesion kinase and inhibitors: a molecular dynamic simulation and free energy calculations

Study on the effect of active ingredients of anti-inflammatory of Smilax china L

Study on the Anti-chronic Pelvic Inflammatory Disease Effective Fraction and its Composition of Smilax china L

TCMID: traditional Chinese medicine integrative database for herb molecular mechanism analysis

SymMap: an integrative database of traditional Chinese medicine enhanced by symptom mapping

TCMSP: a database of systems pharmacology for drug discovery from herbal medicines

PubChem substance and compound databases

Calculation of molecular properties and bioactivity score

Selecting the right compounds for screening: does Lipinski's Rule of 5 for pharmaceuticals apply to agrochemicals?

Synthesis of novel N-(1, 3-thiazol-2-yl) benzamide clubbed oxadiazole scaffolds: urease inhibition, Lipinski rule and molecular docking analyses

SwissTargetPrediction: a web server for target prediction of bioactive small molecules

Pharmmapper server: a web server for potential drug target identification using pharmacophore mapping approach

Relating protein pharmacology by ligand chemistry

UniProt, a hub for protein information

GeneCards Version 3: e Human Gene Integrator

DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants

DrugBank: a knowledgebase for drugs, drug actions and drug targets

ERK MAP kinase pathway induces the inflammatory factor il-8 following 1 chlamydia trachomatis infection

Shrimp anti-lipopolysaccharide factor (SALF), an antimicrobial peptide, inhibits proinflammatory cytokine expressions through the MAPK and NF-κB pathways in LPSinduced HeLa cells

STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets

Cytoscape: a software environment for integrated models of biomolecular interaction networks

Network pharmacology-based preventive effect of XZF on cutaneous toxicities induced by EGFR inhibitor

e RCSB protein data bank: redesigned web site and web services

Mev: multiexperiment viewer

Accurate prediction of the relative potencies of members of a series of kinase inhibitors using molecular docking and MM-GBSA scoring

GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers

PRODRG: a tool for high-throughput crystallography of protein-ligand complexes

GROMOS96 43a1 performance on the characterization of glycoprotein conformational ensembles through molecular dynamics simulations

Siderocalin (Lcn 2) also binds carboxymycobactins, potentially defending against mycobacterial infections through iron sequestration

Small-molecule inhibition of TNF

Crystal structure of rofecoxib bound to human cyclooxygenase-2

Identification of allosteric ERK2 inhibitors through in silico biased screening and competitive binding assay

Discovery of a highly selective chemical inhibitor of matrix metalloproteinase-9 (MMP-9) that allosterically inhibits zymogen activation

Targeting SARS-CoV-2: a systematic drug repurposing approach to identify promising inhibitors against 3C-like proteinase and 2′-Oribose methyltransferase

G_mmpbsa-A GRO-MACS tool for high-throughput MM-PBSA calculations

Material basis and molecular mechanisms of Dachengqi decoction in the treatment of acute pancreatitis based on network pharmacology

e activation and regulation of IL-17 receptor mediated signaling

IL-17 signaling in host defense and inflammatory diseases

Study on signaling and function of IL-17

Reaction of myeloperoxidase compound I with chloride, bromide, iodide, and thiocyanate

IL-17 family: cytokines, receptors and signaling

Agonist-biased signaling via matrix metalloproteinase-9 promotes extracellular matrix remodeling

Doxycycline-mediated inhibition of matrix metalloproteinases improves healing after rotator cuff repair

Screening of antiinflammatory active components of Smilax China L

Spectrum-effect relationship between UPLC fingerprint of Smilax china and anti-pelvic inflammation in rat

Anti-inflammatory and anti-nociceptive activities of Smilax China L. aqueous extract

Human metabolome-derived cofactors are required for the antibacterial activity of siderocalin in urine

Rigorous free energy simulations in virtual screening