Abstract
Objective
To identify blood-absorbed components of Ninglung Deshee (ND) using liquid chromatography-mass spectrometry (LC-MS), and to investigate the mechanism of ND for treating depression using network pharmacology and molecular docking methods.
Methods
Ultra-high-performance liquid chromatography coupled with quadrupole electrostatic field orbitrap tandem high-resolution mass spectrometry (UHPLC-Q-Orbitrap HRMS) was used to identify the components of ND in the blood. Potential target sites of these chemical components were identified using the SwissADME database, while disease-related targets for depression were selected from the GeneCards database. Intersecting targets between the two databases were then identified. Cytoscape software was used to construct a “component-target-disease” network to screen the core components of ND relevant to the treatment of depression. The STRING platform was used to construct a protein-protein interaction (PPI) network of the intersecting targets to identify the core targets of ND for treating depression. Gene ontology (GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed using the DAVID platform. Schrödinger software was used to determine molecular docking of the core components and core target proteins.
Results
Seventy-six blood components associated with 747 drug targets were identified,. Of the 2,000 depression-related targets, 212 overlapped with those of ND. A total of 16 core targets and 16 core components were identified. The GO function and KEGG pathway enrichment analyses indicated that treatment of depression with ND’s may be related closely to the mitogen-activated protein kinase (MAPK) signaling and Prion disease signaling pathways. The molecular docking results showed that piperine had the highest binding energy with tumor necrosis factor (TNF).
Conclusion
The core components of ND such as piperine, liquiritigenin, and piperanine may act on targets including TNF, mitogen-activated protein kinase 3(MAPK3), estrogen receptor 1 (ESR1, and interleukin-6(IL-6), resulting in regulation of signaling pathways such as MAPK, TNF, and Prion disease involved in the treatment of depression. Accordingly, TNF as a therapeutic target and signaling pathway for ND to treat depression warrants further study, while other potential mechanisms and effector molecules of ND for treating depression require further detailed investigation.
1. Introduction
Depression is a major mental disorder 1 that is manifested clinically by a persistent low mood and loss of pleasure. The disorder is often accompanied by cognitive impairments such as lack of concentration and slow thinking, and physical discomforts such as sleep disturbances and decreased energy. 2 In severe cases, depression can lead to self-harm or suicidal thoughts and behaviors. 3 According to the Global Burden of Disease Assessment Report released by the World Health Organization (WHO), depression had already become the second leading cause of disease burden and disability by 2020, and was expected to become the leading cause of global disease burden by 2030. 4
Currently, the drugs used to treat depression mainly include tricyclic antidepressants, monoamine oxidase inhibitors, tetracyclic antidepressants, serotonin and norepinephrine reuptake inhibitors, norepinephrine, and dopamine reuptake inhibitors. In addition, newly discovered target drugs have been introduced in recent years, such as neurosteroids, indoleamine 2,3-dioxygenase inhibitors, and hormone drugs. However, these classic antidepressants generally have adverse effects such as high dependency, a potential for addiction, and easy development of resistance.5,6 Therefore, from the perspective of traditional medicine, there is significant practical value to develop antidepressants that are safe, reliable, and have minimal side effects and consistent quality.
ND is a modified traditional Tibetan medicine, based on the “Srogzin-11” formula. 7 The medicine was modified and formulated by Professor Tsering Namgyal (pinyin: CaiRang Nanjia) according to Tibetan medical theory and years of clinical experience. The formula consists of ten medicinal ingredients, including nutmeg [Myristica fragrans Houtt.], agarwood [Aquilaria sinensis (Lour.) Spreng.], Aucklandiae radix [Dolomiaea souliei (Franch.), C.Shih, (Rhodiola rosea L.), (Licorice [Glycyrrhiza uralensis Fisch. ex DC.]), (Terminalia chebula Retz.), (Fructus chaerospondiatis [Choerospondias axillaris (Roxb.) B.L.Burtt & A.W.Hill]), (Corydalis impatiens (Pall.) Fisch. ex DC.), (caraway [Carum carvi L.]), and (Schisandra chinensis (Turcz.) Baill.). According to Tibetan medicine theory, both ND and Srogzin-11 are used to treat“heart-lung”disease. In Tibetan medicine, the“lung”is a core concept fundamental to sustaining life, with its imbalance manifesting as a range of complex physical and mental disorders. 8 The syndrome described as“heart-lung” disease in Tibetan medicine closely resembles depression, and it is generally considered that depression falls within the scope of“heart-lung” disease.9,10 Based on this concept, both ND and Srogzin-11 are used to treat depression. Pharmacological studies have also confirmed the antidepressant effects of Srogzin-11.11-14 However, the mechanism of action and pharmacodynamic basis of ND remains unclear.
In order to obtain further information on ND the current study used UHPLC-Q-Orbitrap HRMS technology to analyze the constituents present in serum after administration of ND. These constituents were then used to investigate the pharmacodynamic material basis of the beneficial effects of ND on depression. 15 This involved constructing a “drug-target-disease-pathway” network to investigate the multi-component, multi-target, and multi-pathway characteristics of ND and elucidate the mechanisms of action of the absorbed components. This involved predicting the core components and targets of ND and also verifying the binding ability of these core components to the core targets using molecular docking. Taken together, these analyses provided a reference for understanding the pharmacodynamic material basis and mechanisms of action of ND.
2. Materials and Methods
2.1. Reagents and Materials
The study used the following equipment: Q Exactive Plus Orbitrap high-resolution liquid chromatography-mass spectrometer (Thermo Fisher); U3000 ultra-high-performance liquid chromatograph with automatic sampler (Thermo Fisher); Vortex-2 Genie vortex mixer (Scientific Industries); 5810R low-temperature centrifuge (Eppendorf); WD-9415C ultrasonic cleaner (Beijing 61 Instrument Factory); and a Waters ACQUITY UPLC HSS T3 C18 column 2.1 mm × 100 mm, 1.8 μm (Waters).
ND was provided by Tongren Drangsong Tibetan Medicinal Bath Hospital. Methanol, acetonitrile, and formic acid were purchased from Thermo Fisher, ammonium acetate from Sigma, and β-glucuronidase from Shanghai Yuanye Bio-Technology Co., Ltd. Deionized water was prepared using a Milli Q Advantage A10 ultrapure water system. Laboratory hardware, such as 1.5 mL centrifuge tubes, 15 mL centrifuge tubes, and sample vials were all purchased from Axygen.
2.2. Animals and Drug Administration
Specific pathogen-free (SPF)-grade female Sprague-Dawley rats (weighing 180-220 g) were provided by SPF (Beijing) Biotechnology Co., Ltd. (Beijing, China, approval number: SCHK-2024-0001) and housed in an animal room (24 ± 2°C, relative humidity 60 ± 5%) with a 12-h light/dark cycle. The rats were acclimated for one week on standard feed and water before being divided randomly into control and experimental groups (n=6 per group). ND was administered orally at 1.67 g/mL, 1 mL/100 g body weight, twice daily for three consecutive days. The control group received the same volume of physiological saline.
2.3. Collection and Preparation of Drug-Containing Serum
The day before the last day of administration, the rats were fasted for 12 h, but had free access to water. Blood samples were collected from the orbital sinus at 0.5, 1, and 2 h after the last administration. Blood collection at 0.5 and 1 h (0.2 - 0.4 mL each time) was performed quickly by skilled personnel without the use of anesthetics in order to avoid chemical interference in the analysis of blood components. For the final 2-h blood collection (large volume, followed by euthanasia), the rats were anesthetized by an intraperitoneal injection of pentobarbital sodium (30 mg/kg) to reduce stress and ensure the safety of the procedures. The collected blood samples were placed in centrifuge tubes and centrifuged at 3500 rpm for 20 min at 4°C. The supernatant serum was transferred to cryovials, labeled, and stored at -80°C for future use.
To measure the concentration of ND 150 μL of pH 5.0 sodium acetate buffer, 100 μL of the serum sample, and 10 μL of β-glucuronidase (concentration 1000 U/mL) were added sequentially into a 5 mL centrifuge tube, and then incubated with shaking at 37°C for 3 h. β-glucuronidase hydrolyze glucuronide conjugates produced by phase II metabolism in the body, and convert polyhydroxy compounds into their free forms. This ensures that the test results accurately reflect the parent drug concentration in the serum.
After incubation, a 100 μL aliquot of the above mixture was added to 400 μL of a methanol-acetonitrile mixed solution (volume ratio 1:1), and then vortexed to mix thoroughly, followed by centrifugation at 13000 rpm for 10 min at 4°C. The methanol-acetonitrile mixture effectively precipitated proteins in the serum and improved the extraction rate of the analyte. 450 μL of the supernatant was then collected, evaporated to dryness, and reconstituted with 50 μL of 80% methanol solution to obtain the sample for testing.
2.4. UHPLC-Q-Orbitrap HRMS Analysis
The ultra-performance liquid chromatography conditions were as follows. Serum sample analysis was performed on a Thermo Fisher Ultimate 3000 UPLC system, with separation carried out on a Waters ACQUITY UPLC HSS T3 C18 column (2.1 mm × 100 mm, 1.8 μm). The mobile phase consisted of 0.1% formic acid in water (A) and 0.1% formic acid in acetonitrile (B), with gradient elution as follows: 0-10 min, 100% B; 10-20 min, 100%-70% B; 10-25 min, 70%-60% B; 25-30 min, 60%-50% B; 30-40 min, 50%-30% B; 40-45 min, 30%-0% B; 45-60 min, 0% B; 60-60.1 min, 0%-100% B; 60.1-70 min, 100% B. The flow rate for sample separation was 0.2 mL/min, injection volume 10 μL, and column temperature 35°C.
The mass spectrometry data were acquired using a Thermo Fisher Scientific Q Exactive Orbitrap high-resolution mass spectrometer. The detection mode was Full MS-ddMS2, with separate scanning for positive and negative ion modes. The scan range was m/z 100-1500, with MS1 resolution set to 70,000 and MS2 resolution set to 17,500. The ion source voltage was 3.2 kV, capillary temperature 320°C, auxiliary gas heater temperature 350°C, sheath gas flow rate 40 L/min, and the auxiliary gas flow rate 15 L/min. The AGC target was set to 1e6, TopN was set to 5, while the collision energy for triggering MS2 scans used a stepped normalized collision energy (NCE) set at 30, 40, and 50.
2.5. Data Processing
Compound Discover 3.3 software was used for feature peak extraction of the raw mass spectrometry data, element matching of feature peaks, molecular formula prediction, and isotope distribution matching, with all the mass deviations set within 5 ppm. Feature peaks were identified using the online mzCloud database and a locally built mzVault traditional Chinese medicine natural product database (https://www.mzcloud.org/). 16 Identification criteria included a mass deviation of <5 ppm, compliance with isotope distribution, and an mzVault best match database score >70. Relevant literature and PubChem databases (https://pubchem.ncbi.nlm.nih.gov/) were consulted, 17 and the results then compared with the original components reported for ND (Figures S1–S10) and the original components of individual herbal medicines (Figure S11). The identification results were then manually cross-checked and confirmed.
2.6. Prototype Components and Disease Target Prediction
Due to the components absorbed into the blood being more likely to have pharmacological effects, the prototype components identified by UHPLC-Q-Orbitrap HRMS were used to construct the compound information database of ND to support network pharmacology research. The molecular files of the prototype compounds were downloaded from Chemical Book (https://www.chemicalbook.com/ProductIndex.aspx) and PubChem. 17 In addition, the relevant targets of these compounds were obtained from the SwissTargetPrediction database (https://www.swisstargetprediction.ch/). 18 To improve the reliability of the predicted structures, drug targets of the ND compounds with a score value > 0 were selected.
Screening of the depression-related targets from the GeneCards database (https://www.genecards.org/). 19 using a corresponding threshold of ≥16.518 identified the top 2,000 targets ranked by relevance score which were then selected as the research subjects. This threshold was chosen to ensure that the disease target set for subsequent network analysis was focused and of high confidence.
2.7. Network Diagram Construction
The HD component targets and disease targets were jointly imported into Venny 3.10.3 to draw a Venn diagram that showed their intersections. Excel was used to construct a network and type file of active components and intersecting targets that were then imported into Cytoscape 3.10.3 software for visualization of the interaction network and selection of core drug components based on their degree values.
The intersecting targets of ND and depression were then imported into the STRING database (https://stringdb.org/), 20 followed by selection of multiple proteins, setting the species to “Homo sapiens” and the interaction score threshold to 0.4 to generate the PPI network. 21 Cytoscape 3.10.3 software was used to visualize the protein interaction analysis data, with the CytoHubba plugin applied using the following seven algorithms to analyze the data: betweenness, degree, edge percolation component (EPC), maximal clique centrality (MCC), maximum neighborhood component (MNC), radiality, and stress. 22 The top 30 targets from each algorithm were obtained, with the intersection used to identify the core targets, which were then visualized using the MicroBioInfo online bioinformatics analysis and visualization platform (https://www.bioinformatics.com.cn/). 23
2.8. GO and KEGG Enrichment Analysis
The GO and KEGG enrichment analyses of the intersecting targets between ND components and depression were performed using the DAVID database (https://davidbioinformatics.nih.gov/), 24 with the species set to “Homo sapiens”. Fisher’s exact test was used and the Benjamini-Hochberg correction applied to account for multiple testing, witha FDR <0.05 used as the screening threshold.
2.9. Molecular Docking
Protein receptors with no mutations, a resolution of <2.5 Å, and a complete pocket structure were selected from the RCSB PDB database (https://www.rcsb.org/). 25 To validate the reliability of the molecular docking protocol, the co-crystallized ligands of the selected target proteins were re-docked into their respective binding pockets using Discovery Studio 2019 Client. The heavy atom root-mean-square deviation (RMSD) values between the re-docked ligand conformations and the original crystal structures were calculated, using a RMSD < 2.0 Å as the criterion for an effective docking protocol. The selected protein receptors were processed using the Protein Preparation Wizard module of Schrödinger software, which involved steps including import and process, refined hydrogen addition, structure optimization, and energy optimization. The molecular docking pocket was generated using the Receptor Grid Generation module. The 2D structures of the core components were obtained from the PubChem database 17 and imported into Schrödinger software for processing using the LigPrep module. Finally, molecular docking was performed using the Ligand Docking module of Schrödinger software, followed by the generation of molecular docking binding mode diagrams.26-28
2.10. Statistical Analysis
All data processing and analyses were performed using the bioinformatics tools described in the corresponding sections. Identification of blood components was based on mass spectrometry data, processed using Compound Discoverer 3.3 software, with quality control parameters set as described in Section 2.5. Screening of depression-related targets from the GeneCards database used a threshold of relevance score ≥16.518, which corresponded to the top 2000 targets. This threshold was chosen to ensure that the disease target set for subsequent network analysis was focused and of high confidence. For GO and KEGG pathway enrichment analysis, the DAVID database was used to perform a Fisher’s exact test. An adjusted P-value <0.05 (corrected for multiple testing using the Benjamini-Hochberg method) was adopted as the statistical significance threshold to identify significantly enriched terms and pathways. The STRING database was used for protein interaction network analysis, with the interaction score threshold set at 0.4 (medium confidence). Core targets were screened using the CytoHubba plugin of Cytoscape software and involved seven topological algorithms (betweenness, degree, edge percolation component [EPC], maximal clique centrality [MCC], maximum neighborhood component [MNC], radiality, and stress). The intersection of the top 30 targets from each algorithm were selected as the core targets. Schrödinger software was used to calculate the binding energies for molecular docking. The binding energy values were used directly to assess the potential binding affinity between compounds and targets, with lower values indicating stronger predicted binding.
3. Results
3.1. Analysis of Blood Components in Rats Treated With ND
Serum samples were collected at 0.5, 1, and 2 h after the last oral administration of ND to the rats. The total ion chromatograms of the serum samples collected at these times, compared to those of the control serum, are shown in Figure S12.
Basic Information of the 76 Compounds in the Blood Samples
Note. A, nutmeg; B, agarwood; C, Aucklandiae radix; D, Rhodiola rosea; E, Licorice; F, Terminalia chebula; G, Fructus chaerospondiatis; H, Corydalis impartiens; I, caraway; J, Schisandra chinensis.
3.2. Screening of the Chemical Components and Potential Targets
The compounds used in the network pharmacology analysis were selected from the 76 absorbed compounds. A total of 21 ingredients were excluded, including 17 prototype components detected in the rat serum (i.e., choline, betaine, L-tyrosine, 4-guanidinobutyric acid, and uric acid), three components whose compound targets had no overlap with the disease targets (L-aspartic acid, palmitoyl sphingomyelin, and trigonelline), and cyanuric acid for which relevant targets could not be identified. Finally, 55 blood-absorbed ingredients of the ND formula were included, corresponding to 747 unique drug-related targets. The structure of all the compounds used in the subsequent network pharmacology analysis is shown in Figure S13.
We then searched for targets related to depressive disorders using the GeneCards database, and obtained a total of 16,799 targets. The top 2,000 targets with a relevance score ≥ 16.518 were selected as the study subjects. Of these 2000 disorder targets, 212 overlapped with 747 drug targets. A Venn diagram of the overlapping targets is shown in Figure 1. Intersection targets of 55 absorption components of compound ND with depression
3.3. Compound-Intersection Target Network Construction
The 212 intersecting targets identified in the previous step, along with the 55 components in the blood that corresponded to these targets, were imported into Venny 3.10.3 software to construct a “ND chemical components-potential targets for depression” network. Statistical analysis by the software showed that this network consisted of 267 nodes and 736 edges. Network analysis revealed that individual ND compounds interacted with multiple targets, while several targets were influenced by more than one compound. As shown in Table S2, 16 core components (degree ≥ 20) were identified. The interaction network of the core components and intersecting targets is shown in Figure 2. Interaction network of core components and overlapping targets
3.4. PPI Network Construction and Core Target Screening
As shown in Figure 3, the PPI network contained a total of 212 nodes and 2,942 edges. By applying these intersections, 16 core targets were ultimately obtained, including ESR1, AKT serine/threonine kinase 1 (AKT1), axin1/beta-catenin (CTNNB1), glycogen synthase kinase-3 beta (GSK3B), and epidermal growth factor receptor erbB1 (EGFR) (Figure 4). PPI network diagram. Distribution of the core targets.
3.5. GO Function and KEGG Pathway Enrichment Analysis
GO enrichment analysis was performed on the 212 intersecting targets of depression and blood-absorbed components using FDR < 0.05 as the threshold. A total of 549 GO terms were identified, including 334 biological process (BP) terms, 78 cellular component (CC) terms, and 137 molecular function (MF) terms. The top 10 terms in each category were selected for annotation analysis. Representative BP terms included the insulin-like growth factor receptor signaling pathway (GO:0048009), epidermal growth factor receptor signaling pathway (GO:0007173), and positive regulation of the ERK1 and ERK2 cascade (GO:0070374). Representative CC terms included the plasma membrane (GO:0005886), postsynaptic membrane (GO:0045211), and dendrite (GO:0030425). Representative MF terms included identical protein binding (GO:0042802), protein kinase activity (GO:0004672), and protein tyrosine kinase activity (GO:0004713). The detailed results are shown in Figure 5. GO enrichment analysis
KEGG pathway enrichment analysis was performed on the 16 core targets derived from blood-absorbed components using the DAVID database. FDR < 0.05 was used as the screening threshold and 123 significantly enriched pathways were identified. As shown in Figure 6, the top 10 pathways with the smallest FDR values were selected for visualization. In the bubble plot, bubble size represents the number of enriched genes, while bubble color represents the FDR value (redder colors indicate higher significance), and the x-axis represents the fold enrichment. Smaller FDR values indicate stronger associations with depression, and larger bubble sizes indicate more target genes enriched in the pathway. KEGG pathway enrichment analysis of core targets
Based on the above results, the ND components in the blood samples, potential targets for depression, and the corresponding signaling pathways were imported into Cytoscape 3.10.3 software to construct a “compound-target-pathway” network (Figure 7). In this figure, green represents the chemical components of ND, red represents core targets, and orange represents signaling pathways. The connections indicate the relationships between signaling pathways, targets, and components. Network visualization analysis of the results identified 15 core protein targets, including GSK3B, IL6, interleukin-1 beta (IL1B), TNF, and MAPK3, corresponding to 28 components such as 18β-glycyrrhetintic acid, daidzein, liquiritigenin, naringenin chalcone, and sedanolide, and 10 pathways including pathways for cancer, inflammatory bowel disease, the C-type lectin receptor signaling pathway, breast cancer, and hepatitis C. Compound-key target-pathway network
3.6. Molecular Docking
Based on the PPI network analysis, four core targets (ESR1, IL1B, MAPK3, and TNF) were selected. The heavy atom root-mean-square deviation (RMSD) values for these targets were 1.55 Å, 0.83 Å, 1.73 Å, and 0.97 Å, respectively, all of which were below the commonly accepted threshold of 2.0 Å. This finding confirmed the reliability of the molecular docking method and the parameters used in this study. In addition, six core components (dehydrodiisoeugenol, liquiritigenin, naringenin chalcone, piperanine, piperine, and sedanolide) were selected. As shown in Figure 8 molecular docking was performed using Schrödinger software to obtain the binding energies between the ligands and receptors. It is generally believed that the lower the binding energy of molecular docking, the more stable the binding between the ligand and receptor, and also the greater binding ability.
29
The docking results of the six groups with the lowest binding energies are shown in Figure 9. Heatmap of binding energy between the components and targets Well-bound results in the molecular docking analysis.
4. Discussion
Depression is essentially a complex disease driven by multiple factors, with its pathogenesis involving dysregulation of multiple molecular controls that cause abnormalities in multiple cellular signaling pathways.
30
These changes lead ultimately to depression with the high heterogeneity and complex clinical phenotypes creating challenges for precise treatment (Figure 10). Schematic overview of the research workflow and antidepressant mechanism of Ninglung Deshee (ND)
In Tibetan medicine, depression falls under the category of “lung disease”, with its clinical symptoms being very similar to those of Ninglung disease. By drawing from clinical experience, Tibetan medicine often uses formulas to treat ND based on medicines that contain ingredients such as nutmeg, agarwood, Aucklandiae radix, and Fructus chaerospondiatis. Following the early stages of investigation, these treatments have gained considerable clinical recognition.
Our analysis identified MAPK3 as a key node in the ND-target network, consistent with the findings of previous reports. 31 Notably, the co-enrichment of the MAPK signaling pathway and the Prion disease signaling pathway suggests that the ND compound may exert its therapeutic effects by regulating multiple mechanisms related to neuronal survival and function. The MAPK pathway is involved mainly in cell proliferation, differentiation, and stress responses, while the Prion disease pathway is closely associated with neurodegenerative changes. The interplay between these two pathways provides a new perspective for further understanding the neurobiological mechanisms of depression and further confirms the significant role of the MAPK signaling pathway in the treatment of depression using ND.
There is substantial evidence of a strong association between ESR1, a key target of the compound ND, and the treatment of depression. Research has indicated that upregulation of miR-211 and miR-744 in the exosomes of postpartum depression patients may contribute to disease development by targeting ESR1. 32 Furthermore, ESR1 gene variants have been linked to reproductive disorders and the risk of depression, underscoring the importance of estrogen signaling in reproductive and mental health. 33
During the progression of major depressive disorder (MDD) and bipolar disorder (BD), the activation of pro-inflammatory immune responses and upregulation of the inflammatory setpoint are significantly associated with the onset of major depressive episodes (MDE) and severity of the disease. 34 Tumor necrosis factor α (TNF-α), as a key pro-inflammatory factor, not only synergistically regulates apoptosis and inflammatory responses through activation of TNF signaling and NF-κB pathways, 35 but can also be involved in the depressive process by modulating other inflammatory factors such as IL-6. These findings suggest ND may exert anti-inflammatory effects via modulation of TNF-α and IL-6, a possibility that aligns with the holistic approach of Tibetan medicine.
The results of the molecular docking analysis showed that the six core active components in the ND compound exhibited good binding characteristics with both TNF and ESR1 target proteins. Of these components, piperine showed the optimal binding characteristics with TNF, a finding consistent with those of previous studies on the anti-inflammatory effects of piperine. 36 When considering the central role of TNF-α in neuroinflammation associated with depression, it is likely that piperine may alleviate neuroinflammatory responses by directly acting on TNF and inhibiting the activation of downstream inflammatory pathways such as NF-κB.
Although this study provides a reference for the potential mechanisms of ND for treating depression, several limitations should be acknowledged. (1) Network pharmacology relies on the accuracy and completeness of public databases, with database bias possibly affecting the predicted results; (2) The bioactive compounds among the 76 blood-absorbed components that truly exert therapeutic effects remain to be confirmed. In this regard, the presence of a compound in the blood does not guarantee its activity at the target site; (3) Molecular docking is a computational prediction, and docking scores may not always correspond to actual in vivo biological efficacy; (4) The PPI and enrichment analyses were based on specific thresholds and different thresholds may have yielded different results.
5. Conclusions
This study used UPLC coupled with UHPLC-Q-Orbitrap HRMS to identify a total of 76 components of ND that entered the circulation. Combined with network pharmacology analysis, 16 core targets and 15 core components of ND were identified for treating depression. Pathway enrichment analysis showed that the MAPK signaling pathway and Prion disease signaling pathway may be key pathways through which ND exerts its effects. In addition, molecular docking results indicated that TNF and ESR1 play important roles in the antidepressant mechanism of ND. However, these computational predictions require further experimental validation to confirm the roles of these bioactive components and pathways in the treatment of depression.
Supplemental Material
Supplemental material - Study on the Mechanism of ND for Treating Depression Based on “Blood Components-Network Pharmacology”
Supplemental material for Study on the Mechanism of ND for Treating Depression Based on “Blood Components-Network Pharmacology” by Ma-jia Wan, Yecuo Lamao, Name Zhaxi, Jie-ben Jiao, Cairang Nanjia in Natural Product Communications
Footnotes
Acknowledgement
The authors thank Beijing Keweite Animal Technology Co., Ltd. for their support in the ethics review of the animal experimental procedures and animal quality certification (License No. SCXK (Jing) 2024-0001). We also appreciate the assistance provided by the animal facility personnel during the in vivo experiments. The authors gratefully acknowledge EditSprings (
) for their professional language editing services.
Ethical Considerations
This study was approved by the Laboratory Animal Welfare and Ethics Committee of Beijing Keweite Animal Technology Co., Ltd.
Consent to Participate
There are no human subjects in this article and informed consent is not applicable.
Author Contributions
Ma-Jia Wan; conceptualization, methodology, investigation, data curation, visualization, writing-original draft, and writing-review and editing. Lamao-Yecuo; data curation, validation, and writing-review and editing. Zhaxi-Name; data curation, validation, writing-review and editing. Jie-ben Jiao; data curation, and validation. Cairang-Nanjia; conceptualization, writing, review, editing, and supervision.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Social Science Fund project (Grant Nos: 20XMZ026) and the Key Laboratory for Tibet Plateau Phytochemistry of Qinghai Province.
Declaration of Conflicting Interests
This manuscript has been confirmed by all authors and has not been published or submitted for review in any domestic or international journals, nor is it currently under review by any other journal. All the authors have reviewed and agreed to the final content of the manuscript and the order of authorship, and declare that there are no conflicts of interest.
Data Availability Statement
Data will be made available on request.
Statement of Human and Animal Rights
The animal experimental procedures involved in this study were conducted in accordance with the guidelines for animal experiments of Beijing Keweite Animal Technology Co., Ltd. and approved by its Laboratory Animal Welfare and Ethics Committee.
Supplemental Material
Supplemental material for this article is available online.
Appendix
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
