The Pharmacological Mechanisms of Mori Folium and Eucommiae Cortex Extracts against Immunosuppression from Network Pharmacology Perspectives

1. Introduction

The immune system plays an important role in anti-infection, anti-tumors as it acts as a body’s defense system by protecting our body cells, tissues, and organs from invading infections by harmful microorganisms and other disease-causing microbes [1,2]. Immunomodulation is the interaction between immune cells and immune molecules in the immune system, as well as with other systems such as the neuroendocrine system, to maintain the body at the most appropriate level. It involves two major mechanisms: immune-stimulation and immunosuppression [3], where the immunosuppression refers to the inhibition of the immune response and is becoming increasingly common, especially when patients received transplanted organs or bone marrow [2].

Morus alba L., belonging to the Moraceae family, is one of the most valuable plants and rich in natural products, whose pharmaceutical name is Mori Folium (MF) [4,5]. Mori Folium contains a variety of active ingredients, such as polysaccharides [6], flavonoids [7], alkaloids [8], etc. The comprehensive effects of these active ingredients reflect the pharmacological effects of MF, including antidiabetic, anti-inflammatory, antibacterial, cardiovascular and cardioprotective, hypolipidemic, antioxidant, and antiatherogenic abilities [9,10].

Eucommiae Cortex (EC) is derived from the dried bark of Eucommia ulmoides Oilv., where active compounds are separated and analyzed by modern chemical methods, and a total of 112 compounds are identified, mainly including lignans, iridoids, phenolics, steroid and terpenoids and flavonoids [11]. Currently, many studies have shown that the pharmacological effects of EC are antihypertensive, hypolipidemic, anti-obesity, antidiabetic, neuroprotective, antioxidative, antifatigue, anti-aging, antitumor, anti-inflammatory and enhancing immune-function [11].

Network pharmacology is a new field combining traditional Chinese medicine (TCM) and network pharmacology, which is based on the theory of system biology and the viewpoint of network pharmacology. Network pharmacology facilitate the under-standing of the compatibility law of prescriptions, identification of TCM medicine, prediction of disease-related targets and the action mechanism of TCM [12]. In this study, the network pharmacology, molecular docking and molecular dynamics simulations methods were used to explore the pharmacological and molecular mechanisms of the anti- immunosuppression activity of Mori Folium and Eucommiae Cortex (MFEC) extracts. And the underlying mechanism of MFEC would provide a new insight into the screening of potential bioactivity and facilitate the development of drugs for anti-immunosuppression treatment from the active compounds.

2. Materials and Methods

2.1. Animal experiments

The animals were housed in the animal house of Anhui Agricultural University, with commercial standard diets and water ad libitum. And among the experimented animals, the half is males and the other half is females. Besides, in current study, all animal experiments completely comply with the ARRIVE guidelines and conducted in accordance with the U.K. Animals (Scientific Procedures) Act, 1986 and associated guidelines, EU Directive 2010/63/EU for animal experiments, or the National Research Council's Guide for the Care and Use of Laboratory Animals.

In this study, fifty mice were randomly divided into five groups, with each group of ten mice: control group, model group, MF group, EC group and MFEC group. After one week accommodation, on the 8th-10th day, except the control group, the other groups were intraperitoneally injected with cyclophosphamide at the rate of 80 mg/kg once a day. The control group was injected with the same amount of normal saline. On the 11th-17th day, MF group, EC group and MFEC group were fed with 200 mg/kg MF, EC, MFEC respectively. The other groups were fed with the same amount of normal saline. All the experiments pertaining to animals comply with the commonly-accepted “3R”, according to the guideline of Anhui Agricultural University about the protection of animals. The spleen of mice was collected and weighted (g), and the spleen index of mice was calculated as follows:

The spleen index (%)=spleen weight (g)/body mass of mice (g) x100%

Additionally, the IF-2, IF-6 and TNF-α level in serum was determined according to the ELISA kits for mouse (Shanghai Jianglai Biotechnology Co., Ltd., Shanghai, China), and 100 µL of negative control, standard or diluted serum sample to be tested were added to the microtiter plate to make duplicate wells. The microtiter plate was covered with film, placed horizontally, and incubated at room temperature for 60 min. After discarding the liquid, wash the plate 4 times with washing solution, dry the plate after each operation. Add 100 μL of enzyme-labeled antibody to each well, add the mem-brane, and incubate at room temperature in the dark for 30 min. Similarly, after dis-carding the liquid, wash it with washing solution for 4 times, and dry the plate after each operation. Add 100 μL of TMB substrate solution to each well, and incubate at room temperature in the dark for 10 min, then add 100 μL of stop solution to each well to stop the reaction. The optical density value (OD value) of each well was measured at the wavelength of 450 nm by microplate reader.

2.2. Database construction of active ingredients and potential targets

Chemical ingredients of MFEC were obtained from Traditional Chinese Medicine Database and Analysis Platform (TCMSP) (https://old.tcmsp-e.com/tcmsp.php). The active compounds were selected according to the pharmacokinetics parameters including absorption, distribution, metabolism, and excretion. Basically, the criteria for the related parameters were set follows: oral bioavailability (OB), not less than 30% and drug-likeness (DL), not less than 0.18.

The potential targets of the active ingredients were obtained from TCMSP, TCM-ID (http://119.3.41.228:8000/tcmid/), PubChem (https://pubchem.ncbi.nlm.nih.gov/) and Swiss Target Prediction database (http://www.swisstargetprediction.ch/). Also, all the targets were integrated and the duplicates were removed.

2.3. Acquisition and screening of immunosuppression-associated targets

GeneCards database (https://www.genecards.org/) was used to extract the immunosuppression-related targets, and the median value of relevance score was used to screen the obtained targets. Only the targets with inference score no less than 9.33 were included. Besides, the Uniprot database (https://www.uniprot.org/) was used to convert the obtained disease-related targets and the potential targets of the active ingredients obtained in Section 2.2 into the gene symbol formats. By merging all acquired genes, all targets related to immunosuppression were collected to establish a gene-library of anti-immunosuppression targets. The Venny 2.1 tool (https://bioinfogp.cnb.csic.es/tools/venny/) was used to determine the number of intersection genes between MFEC potential targets and disease-related genes, and to plot a Venn diagram.

2.4. Network construction and topological analysis

The network construction and topological analysis were investigated using String 11.5 database (https://string-db.org/) and Cytoscape 3.8.1 software. Intersection ingredients-targets network was constructed to understand the associations between active ingredients and intersection targets of MFEC. Network generation and visualization were performed using Cytoscape software. In order to construct PPI network, common targets were imported into String 11.5 database to generate a network and the topological analysis was conducted using Cytoscape software, where the degree centrality, closeness centrality, as well as the betweenness centrality were obtained to evaluate the central properties of nodes in the network.

2.5. GO and KEGG pathway enrichment analysis

To further elucidate pharmacological mechanisms of MFEC in immunosuppression treatment, Metascape database (https://metascape.org/), a web-based portal where covers a comprehensive gene list annotation, was used to perform Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The results were analysed and sorted according to the adjusted P value. Additionally, the target-pathway network was established using Cytoscape software.

2.6. Molecular docking verification

The molecular docking was performed to predict the binding affinity and interactions between the targets with the active ingredients. Based on the topological parameters of the ingredients-targets and PPI networks, ingredients and targets with more than 2-fold the median value of degree were selected as candidate ligands and recep-tors for molecular docking. The software ChemBio3D 19.0 was used to calculate minimizing energy and the ligand structures were exported by MOL 2 format. The 3D structures of selected receptors were screened and downloaded from the RCSB PDB database (https://www.rcsb.org/,). The Sybyl-x 2.1.1 software could perform the molecular docking of active compounds with putative targets. All the softwares and databases used in this study were described in supplementary information Table S4.

2.7. Molecular dynamics simulations

The above part of molecular docking results (arachidonic acid-HSP90AA1; Dehydrodieugenol-HSP90AA1; icosa-11,14,17-trienoic acid methyl ester-HSP90AA1; (-)- Tabernemontanine-AR; beta-carotene-AR; beta-sitosterol-AR) were performed by Gromacs 2020.6 software for molecular dynamics simulations. The simulation system was adjusted to a sodium chloride solution at 37℃ in order to replicate the actual human in vivo environment. Considering the active ingredient of MF and EC as multi-carbon ring skeleton structures, Charmm36 force field and TIP3P water model was chosen for MDs. During the MDs, the relevant hydrogen bonds are constrained by the Linear Constraint Solver (LINCS) algorithm with the integration step of 2 fs, while the non-bond interaction cut-off value is set to 10 Å and updated every 10 steps. And the electrostatic interactions are calculated using the Particle-mesh Ewald (PME) with a cut-off value of 1.2 nm. In order to optimize the original conformation of the protein in the solvent, the protein receptor-ligand small molecule complex was pre-equilibrated for 100 ps prior to simulations, followed by canonical ensemble (NVT) equilibration for 100 ps using a modified Berendsen temperature coupling algorithm with the coupling time constant of 0.1 ps, which allowed the complex to be warmed up to 310 K with the solvent system. Subsequently, the solvent and complex were pressure equilibrated, while the pressure was increased to 1 bar using a Berendsen constant pressure with constant-pressure and constant-temperature (NPT) equilibration of 100 ps. Finally, MDs of the complex were performed for 50 ns.

3. Results

3.1. Mice spleen index, IL-2, IL-6, and TNF-α level in serum

It can be seen from Figure 1 that compared with the control group, the spleen index of the mice in the model group decreased. After treatment with MF, EC and MFEC, the spleen index of the mice increased compared with the model group, and the spleen index of the MFEC group increased more obviously, but statistical analysis showed no significant difference. Comparing the levels of TNF-α in all groups, the content of TNF-α in the mf group was slightly lower than that in the control group, but there was no significant change in the other groups. Taken together, these data indicate that the drug may have a certain therapeutic effect on the immunosuppression of mice.

3.2. Database construction of active ingredients and potential targets

Based on the criteria of the OB≥ 30% and DL≥ 0.18, a total of 53 active ingredients (25 in MF, 24 in EC, 4 common ingredients, in supplementary materials Table S1) were retrieved in TCMSP database. All the potential targets of active ingredients were found in TCMSP, TCM-ID, Swiss Target Prediction and Pubchem database. After removing duplicates, 423 potential targets were screened out of the two traditional Chinese medicines.

3.3. Immune suppression-related targets

The immunosuppression-associated targets were retrieved in the GeneCards database. As a result, a total of 11461 disease-related genes were obtained. Based on 3 times screenings of the median value, 1435 disease-related genes were selected. Furthermore, 204 common targets (in supplementary materials Table S2) between MFEC targets and disease-related targets were filtered as the critical targets for further study (Figure 2A).

3.4. Intersection compound-targets network

The compound-target interaction network was visualized in Figure 2B with 259 nodes and 665 edges by Cytoscape. In the network, nodes represent the screened drug, active ingredients and genes, while the connections between the nodes represent the interactions between these biological analyses. In the current study, active compounds with degree 2-fold higher than the median value were selected to perform molecular docking. These 11 active compounds were shown in supplementary materials Table S3.

3.5. PPI network

To analyze the protein-protein interaction, the 204 common targets of ingredients and diseases were imported into the String database and Cytoscape to show complex interaction between proteins encoded by these targets. The PPI network covered 203 nodes (the protein IGHG1 was not identified) and 4897 edges (Figure 2C). Among the 203 proteins, AKT1, AR, CASP3, HSP90AA1, JUN, MAPK14, MMP2, PTSG2, TNF, were hub proteins finally selected as putative targets to conduct molecular docking.

3.6. GO and KEGG enrichment analysis

3.6.1. GO enrichment

GO enrichment analysis was performed using Metascape software on 204 common targets of immunosuppression and Mori Folium and Eucommiae Cortex in terms of biological process (BP), cellular composition (CC) and molecular function (MF)[13]. As seen in Figure 3A, Biological process mainly included positive regulation of vitamin D biosynthetic process, positive regulation of calcidiol 1-monooxygenase activity, response to carbon monoxide, activation of cysteine-type endopeptidase activity involved in apoptotic signaling pathway and positive regulation of apoptotic process involved in morphogenesis. The cellular composition mainly enriched in Bcl-2 family protein complex, cyclin/CDK positive transcription elongation factor complex and nuclear cyclin-dependent protein kinase holoenzyme complex. The molecular function mainly involved phosphatidylinositol-3,4-bisphosphate 5-kinase activity, BH3 domain binding, ErbB-3 class receptor binding, 1-phosphatidylinositol-4-phosphate 3-kinase activity, phosphatidylinositol-4,5-bisphosphate 3-kinase activity.

3.6.2. KEGG enrichment

KEGG pathway analysis could indicate possible signal pathway of 204 proteins, and there were 362 pathways been obtained totally. Figure 3B and supplementary material Table S5 showed specific information about the underlying mechanisms involved in MFEC to treat immune suppression. The top three enriched KEGG pathways were: pathways in cancer, AGE-RAGE signaling pathway in diabetic complications, fluid shear stress and atherosclerosis.

3.7. Molecular docking

9 proteins were selected according to the PPI network and the structure of these nine proteins were downloaded from the PDB database. Sybyl-X software was used to remove water and excess ligands in these 9 proteins and then hydrogenated. The last step in pre-processing the protein was to generate binding pockets. These 9 proteins were docked with 11 active components (obtained in 3.4), the related results were shown in Table 1, Figure 4A,B. In Table 1, we could see that among all 11 active compounds, iristectorigenin A, quercetin, acid methyl ester, arachidonic acid might be an important active compound. Among all 9 protein targets, iristectorigenin A had better docking results with AKT1 and AR, and there were three docking sites with each target. Quercetin and CASP3 docking have five binding sites, which is the most in all dockaing. Acid methyl ester binds to HSPAA1 and TNF in a docking site, respectively. Arachidonic acid has the best docking effect with four of the screened targets and has a high score. In addition to two binding sites with PTGS2, arachidonic acid has one docking site with JUN, MAPK14 and MMP2.

3.8. Molecular dynamics simulations

To further validate the binding stability of drug active ingredients with AR and HSP90AA1 in vivo, the analysis of Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation(RMSF), Radius of Gyration (Rg) and Hydrogen bond (Hbond) were obtained by Gromacs 2020.6 software. RMSD was used to analyze the binding stability of the AR and HSP90AA1 to the corresponding ligand receptor in the MDs. Figure 5A showed that the (-)-Tabernemontanine-AR complex reached dynamic equilibrium within a short time (10 ns) and the RMSD value remained around 0.12 nm until the end of the simulation, indicating that AR bound well with (-)-Tabernemontanine and can form a stable complex. Meanwhile, the binding of AR with beta-sitosterol and beta-carotene reached the equilibrium at 10-15 ns, and its RMSD fluctuation value was less than 0.3 nm, which suggested that was stable drug-protein complex binding. Figure 5B suggested that arachidonic acid-HSP90AA1 and icosa-11,14,17-trienoic acid methyl este-HSP90AA1 always had large fluctuations in the RSMD curves during the MDs, indicating that the active ingredient of the drug could not bind stably to the protein target. RMSF was used to reflect the fluctuations and the degree of motion drasticity of protein residues throughout the MDs. Figure 5C,D showed that AR and HSP90AA1 residues had greater flexibility and adaptability when key drug active ingredients were bound to hub protein target.

Rg indicated the tightness of the overall structure of the protein [14], as well as characterized the changes in peptide chain relaxation of the protein during the simulation. Figure 6A indicated that the complexes formed by AR bounds to corresponding three drug active molecules had more stable radius of gyration and they gradually decreased in the MDs, which suggested that the proteins gradually converge and the complex structures become stable after binding to the small molecules. The Rg curves of HSP90AA1 showed some fluctuations (Figure 6B). As was seen in Figure 6C, beta-sitosterol-AR complex formed hydrogen bonds only in part of the time period, and (-)-Tabernemontanine-AR complex basically remained more than 1 hydrogen bonds during the simulation. The interaction between beta-carotene with AR never formed hydrogen bond (Figure 6C) because beta-carotene-AR complex structure haven’t hydrogen bond and the interaction be dominated by hydrophobic interactions. Figure 6D showed that the number of hydrogen bonds of arachidonic acid-HSP90AA1 complex were low throughout the MDs. Figure 6D showed that the number of hydrogen bonds formed by the icosa-11,14,17-trienoic acid methyl ester-HSP90AA1 complex had 1-2 hydrogen bonds. Dehydrodieugenol-HSP90AA1 complex maintained one hydrogen bond for most of the time throughout the MDs, which indicates that Dehydrodieugenol is more stable in binding to HSP90AA1 combined with RMSD analysis.

4. Discussion

In recent decades, TCM has attracted worldwide attention due to its exact curative effects, relatively less low toxicity and low cost. Due to the complexity components of TCM and the various biological systems where they were involved, how to elucidate the mechanism of action has become a challenge [12,15]. Network pharmacology takes the main active ingredients and target proteins of TCM as nodes, and edges to represent their interactions by generating an interaction network, to elucidate the mechanism of action of TCM prescriptions from the molecular level [16].

In the present study, the network pharmacology is used to explore the possibility of MFEC in the treatment of immunosuppression. The obtained relevant targets were used to dig out the putative related pathways, and to speculate the possible mechanism of MFEC treatment of immunosuppression. Also, molecular docking was used to verify the key components and targets of MFEC involved in the treatment of immunosuppression.

During the active component data-mining of MFEC, a total of 53 active components were collected. Finally, 11 critical components were selected according to the PPI and topological parameters. The 11 components were icosa-11,14,17-trienoic acid methyl ester, iristectorigenin A, arachidonic acid, tetramethoxyluteolin in MF; dehydrodieugenol, (-)-tabernemontanine, (9R)-6'-methoxycinchonan-9-ol in EC, and their common active ingredients quercetin, kaempferol, beta-sitosterol, beta-carotene. These 11 compounds mainly belong to flavonoids, alkaloids, phenolics, terpenes, phytosterols and organic acids that play pharmacological effects in MFEC. Both flavonoids and phenolic active ingredients in MFEC belong to polyphenolic compounds and they contain a large number of aromatic hydroxyl groups, which account for their anti-inflammatory, anti-glycemic pharmacological effects, and prevent oxidative dam-age and cell death [11,17,18]. In addition, studies have shown that kaempferol can alleviate neuronal damage caused by kainic acid-induced epilepsy in mice, restore the proliferation of T lymphocytes to a certain extent, reduce the phagocytic function of macrophages, and inhibit the apoptosis of thymocytes [19]. Besides, kaempferol has certain immunomodulatory effects [19]. To explore the protective effect of organic ac-ids on inflammatory injury in acute tracheobronchitis, LPS wad used to establish respiratory inflammation in mice [20]. Results showed that organic acids could treat acute tracheobronchitis by regulating the TLR4/NF-κB signaling pathway, indicating that organic acids have certain anti-inflammatory effects.

The pharmacological effects of alkaloids are very extensive, and can be used for anticancer, antimalarial, antiviral, antihypertensive, antispasmodic, antiarrhythmic, analgesic, antibacterial and antidiabetic, as well as central nervous system stimulants, muscle relaxants, with vasodilatory properties [21]. The beta-sterols in MFEC belong to the class of phytosterols that have been shown to have cholesterol-lowering, anti-cancer, anti-atherosclerotic, anti-inflammatory and antioxidant effects [22]. Be-ta-carotene belongs to the group of carotenoids in terpenoids and has been proven to be the most abundant pigment and fat-soluble antioxidant in nature [23]. In addition, studies have demonstrated that the addition of fully oxidized beta-carotene to the feed can enhance the immunity and performance of sows [24].

Comprehensive analysis of GO enrichment results showed that the active com-ponents in MFEC may act on organisms through nucleus, cytoplasm, membrane, cyto-sol and organelle membrane, exert the molecular functions of phosphatidylinosi-tol-3-kinase (PI3K) family protein kinases, BH3 domain binding, nitric oxide and other synthase regulators, vascular endothelial growth factor (VEGF) activator receptors, RNA polymerase 2 transcription coactivator, and so on, participate in the reaction to carbon monoxide and iron ions, positively regulate vitamin D biosynthesis, monooxygenase activity and cell apoptosis in vivo, negative regulation of myosin light chain phosphatase and other series of biological processes, so as to achieve the goal of therapeutic immunosuppression. KEGG enrichment analysis showed that the main signaling pathways involved in the common targets of MFEC and immunosuppression included pathway in cancer, advanced glycation end products and their receptors (AG-Es-RAGE) signaling pathway in diabetic complications, PI3K-Akt signaling pathway, tumor necrosis factor (TNF) signaling pathway, fluid shear stress and atherosclerosis, kapos sarcoma-associated herpes virus infection, hepatitis B and C, prostate and pancreatic cancer. AGEs are heterogeneous glycation products of proteins, lipids and nucleotides, and their receptors are called RAGE, which are multiligand transmembrane receptors of the immunoglobulin superfamily [25]. When the body's redox homeostasis is disrupted, ROS accumulate, resulting in oxidative stress. Studies have shown that some active components in TCM could improve the level of intracellular oxidative stress, reduce the content of AGEs and ROS, and down-regulate the level of oxidative stress by inhibiting the AGEs-RAGE signaling pathway [26]. From the CC analysis, PI3K is a cytoplasmic lipid kinase, which is a heterodimer composed of a regulatory subunit p85 and a catalytic subunit p110. And the active PI3K phosphorylates to generate the second messenger phosphatidylinositol 3,4,5-triphosphate (PIP3), which further induces the phosphorylation of AKt, and thus participates in the regulation of various life processes such as growth, apoptosis, and oxidative stress [27]. TNF is one of the most well-studied cytokines in the immune system, which can regulate various processes such as cell communication, differentiation and death. These regulatory functions are related to various diseases of the body, including the body's autoimmunity [28].

The nine protein targets selected for molecular docking in this study were jointly selected by combining the analysis results of compounds-targets network and PPI network. AKT1, one of the AKT isoforms, is an oncogene that is ubiquitous in neurons as an important downstream substrate in PI3K-Akt signaling pathway [29]. AKT1 also plays a key role in the normal development of the nervous system and memory formation when participating in the regulation of cell survival and growth [30]. When AKT1 is mutated, it can increase the risk of schizophrenia [31]. The androgen receptor (AR) is a ligand-activated nuclear receptor that interacts with the estrogen receptor (ER), glucocorticoid receptor (GR), progesterone receptor (PR) and mineralocorticoid receptor (MR), both belong to the type I nuclear receptor subfamily [32]. AR coordinates the expression of androgen-regulated transcriptomes in the nucleus and is critical for prostate development, homeostasis, and carcinogenesis. Studies have shown that AR outside the nucleus can form a complex with Akt under androgen stimulation, and induce the phosphorylation and activation of Akt [32,33]. CASP3 is one of the final effector proteins in the apoptotic response. Tian et al. [34]have shown that under hypoxia and nutrient deprivation conditions, the activity of CASP3 in human nucleus pulposus-derived mesenchymal hepatocytes can be enhanced to promote cell apoptosis. HSP90AA1 is the α isoform of the heat shock protein (HSP) family with a molecular weight of 90 kDa. HSP90 has the function of maintaining protein stability and is involved in the arrangement and maintenance of almost the entire cytoskeleton [35]. Additionally, HSP90AA1 can promote autophagy through the PI3K/Akt/mTOR signaling pathway, and inhibit apoptosis through the JNK/P38 signaling pathway to im-prove the drug resistance of osteosarcoma cells [17]. JUN is an important component of the transcriptional activation protein complex AP1 (Activation Protein-1, AP1), which is widely involved in many biological processes such as cell proliferation, differentiation and apoptosis. JUN can form dimers by combining with other components of the AP1 complex to jointly regulate the expression of downstream target genes [36]. MAPK14 belongs to the mitogen-activated protein kinase (MAPK) family of proteins, enhancing the formation of tumor platelet aggregates that interact with lung endothelium to form lung metastases [37]. MMP2 belongs to the family of matrix metalloproteinases (MMPs), which are known as inflammatory mediators [38]. In the process of atherosclerosis (AS), MMP2 can degrade a variety of collagen and basement mem-brane components, lyse the collagen fibers of AS plaques, and reduce the thickness of the fibrous cap at the plaque. high important evaluation factor [39]. PTGS2, also known as cyclooxygenase 2 (COX-2), is a pro-inflammatory factor that is normally not expressed in most tissues. Its expression is induced by a variety of stimuli, such as epi-dermal growth factor (EGF), interleukin-1 (IL-1), and TNF, and inflammatory responses may also induce the production of PTGS2 [40,41]. TNF mainly refers to TNF-α, as the first inflammatory factor in the inflammatory response, TNF can activate lymphocytes and neutrophils, and play a crucial role in initiating and expanding the inflammatory cascade [42].

Molecular docking revealed drug-target and protein-target interactions from the molecular perspective, and MDs have been widely used to study the binding stability between proteins and molecules and evaluate the structural features of protein-ligand systems [43]. Therefore, in this study, we adopted a combination of molecular docking and molecular dynamics simulation to screen the key targets of ML and EC for enhancing immune function and to verify the stability of active ingredient binding to protein targets. The molecular docking results suggests that iristectorigenin A, quercetin, acid methyl ester, arachidonic acid exhibited great binding affinity to 9 hub targets of immunosuppression. From the analysis of the results of RMSD, RMSF, RG and Hbond of MDs, it is clear that (-)-Tabernemontanine-AR, beta-sitosterol-AR and Dehydrodieugenol-HSP90AA1 complexes bind stably. Molecular docking and further MDs validation finally identified AR, MAPK14, TNF and HSP90AA1 as the core targets for immunosuppression.

5. Conclusions

Mori Folium and Eucommiae Cortex extracts have been demonstrated to be effective for the treatment of immunosuppression. However, the underlying mechanism is unclear. In the current study, network pharmacology, molecular docking, molecular dynamics simulations and experiment validation were conducted to decipher the pharmacological mechanisms of MFEC extracts. Our results showed that the interactions were in consistent with the published literatures and could be employed to the treatment for immunosuppression. Generally, network pharmacology analysis is applied to investigate the pharmacological mechanism of TCM. We reckon that it is a good alternative to explore mechanism by combing network pharmacology analysis, molecular docking and molecular dynamics simulation multiple pathways work together to enhance the body’s immune function and improve immune suppression. This study provides a theoretical basis for in-depth exploration of the pharmacological effects of MFEC and its mechanism of action in the treatment of immunosuppression.