Systems Pharmacology and Molecular Docking Reveals the Mechanisms of Nux Vomica for the Prevention of Myasthenia Gravis

Qiu, Chao; Chen, Qiang; Hou, Qun; Qi, Guanshu

doi:https://doi.org/10.1155/2022/9043822

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Natural Products and Herbal Medicine: Safety, Efficacy, and Bioactivity

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 9043822 | https://doi.org/10.1155/2022/9043822

Systems Pharmacology and Molecular Docking Reveals the Mechanisms of Nux Vomica for the Prevention of Myasthenia Gravis

Chao Qiu,¹Qiang Chen,¹Qun Hou,¹and Guanshu Qi¹

Academic Editor: Shoib Baba

Received10 May 2022

Accepted20 May 2022

Published26 Jun 2022

Abstract

Background. Myasthenia gravis (MG) is a rare autoimmune disease with clinical symptoms of fluctuating muscle weakness. Due to the side effects of current therapies, there is an urgent need for a new medication for MG treatment. Nux vomica is a traditional Chinese medicine used in various diseases. However, the mechanism of action of Nux vomica against MG remains unclear. Methods. Network pharmacology was used to explore the underlying mechanisms of Nux vomica in MG treatment, which was validated using molecular docking and in vivo experiments in mice. Results. Twelve bioactive compounds and 72 targets in Nux vomica were screened. Seventy-nine myasthenia-related targets were obtained from the GENECARD and DisGeNET databases. PPI networks of Nux vomica- and myasthenia-related targets were constructed using Bisogenet, and these two networks were subsequently merged into an intersection to establish a core-target PPI network that consisted of 204 nodes and 4,668 edges. KEGG enrichment analysis indicated that 132 pathways were enriched in 204 core targets. In addition, we obtained 50 docking pairs via molecular docking. In vivo experiments revealed that Nux vomica can improve the symptoms of MG. Conclusion. Nux vomica is involved in the pathogenesis of MG through the “multicomponent-target-pathway” mechanism.

1. Introduction

Myasthenia gravis (MG) is an autoantibody-mediated neuromuscular disease that causes fluctuating muscle weakness predominantly in the ocular, bulbar, and proximal limb muscles [1]. The prevalence rate is approximately one in 5,000 worldwide, and the incidence is higher in younger women and older men [2]. Currently, patients with MG are mainly treated with corticosteroids, long-term immunosuppressive drugs, plasma exchange, or thymectomy [3]. However, their clinical application is limited due to their high costs and side effects, such as cardiac rhythm disorders and hypotension [4]. Hence, alternative therapies with higher efficacy and fewer side effects are needed for MG treatment.

Traditional Chinese medicines (TCMs) have received increasing attention because of their high efficacy, safety, and few adverse effects. According to the TCM theory, MG is classified as “Weizheng” and diagnosed as “flaccidity syndrome” and “Qi-deficiency pattern” [5]. Thus, Qi-supplementation therapy is important for MG treatment. An increasing number of studies have shown that TCMs as an alternative or adjuvant of Western medicines could significantly improve the clinical efficacy of MG treatment [6–8]. Jianpiyiqi granules have been shown to remarkably alleviate the clinical symptoms in children with ocular MG accompanied by a decrease in acetylcholine receptor autoantibodies [6]. In addition, the Qiangji Jianli decoction has been shown to alleviate oxidative stress injury and facilitate mitochondrial biogenesis via the AMPK/PGC-1α-dependent signaling pathway, which restored muscle energy supply and ameliorated gastrocnemius muscle contractility in experimental autoimmune MG rats [7]. Another well-known TCM treatment for MG, Bu Zhong Yi Qi decoction, has been reported to decrease simulated weightlessness-induced muscle atrophy and promote myogenesis by reducing nuclear receptor corepressor 1 (NCoR1)-associated gene expression [8]. Nux vomica, also called Semen Strychni or Ma qianzi (Chinese: 马钱子), is the seed of Strychnos Nux-vomica L, belonging to the genus Strychnos of the family Loganiaceae. Studies have illustrated that Nux vomica has numerous biological activities, such as analgesic, anti-inflammatory, antitumor, and immunomodulatory effects [9]. A previous study found that Nux vomica could alleviate inflammatory response and thereby ameliorates MG in experimental autoimmune MG rats [10]. However, the molecular mechanisms of Nux vomica against MG remain unclear. Therefore, it is meaningful to explore the targets and further clarify the underlying mechanisms of action of Nux vomica on MG treatment based on network pharmacology.

Network pharmacology, first introduced by British pharmacologist Hopkins, has proven to be an effective strategy for revealing the complicated relationship between TCMs and diseases [11, 12]. TCMs exert their therapeutic functions through a multicompound-target-pathway network. Thus, the molecular mechanisms of TCMs can be clarified using network pharmacology systematically and integrally [13]. In this study, we screened the active components of Nux vomica and explored its underlying molecular mechanisms on MG based on network pharmacology. Our study provides a theoretical basis for the future clinical application of Nux vomica for treating MG.

2. Materials and Methods

2.1. Collection of Therapeutic Targets for MG

Therapeutic targets for MG were obtained from GeneCards [14] (https://www.genecards.org/) and DisGeNET [15] (https://www.disgenet.org/) databases. Candidate targets were collected after the deletion of duplicates.

2.2. Collection of Bioactive Ingredients and Therapeutic Targets for Nux Vomica

TCMSP (https://tcmspw.com/tcmsp.php) is an online platform for integrating pharmacokinetics, medicinal chemistry, therapeutic targets, and drug-target-disease networks [16]. Based on the TCMSP platform, bioactive ingredients and targets of Nux vomica were obtained using the screening criteria of oral bioavailability (OB) > 30% and drug-likeness (DL) > 0.18 [17]. In addition, strychnine was considered as a candidate bioactive ingredient based on a previous study. The pharmacokinetics and DL properties of strychnine were further assessed using the online SwissADME database (https://www.swissadme.ch/index.php). The therapeutic targets of strychnine were then predicted using the SwissTargetPrediction online tool (https://swisstargetprediction.ch/).

2.3. PPI Network Construction

A PPI network was constructed and visually analyzed using the Cytoscape plugin Bisogenet, which contains the following six PPI databases: the IntAct Molecular Interaction Database (IntAct), Biomolecular Interaction Network Database (BIND), Biological General Repository for Interaction Datasets (BioGRID), Database of Interacting Proteins (DIP), Molecular INTeraction Database (MINT), and Human Protein Reference Database (HPRD) [18]. PPI networks of the Nux vomica-related targets and MG-related targets were constructed using Cytoscape.

2.4. Central Network Evaluation

Central network evaluation is considered a primary method for core protein screening in PPI networks. First, the PPI networks of the Nux vomica-related targets and MG-related targets were merged into an intersection. Next, the Cytoscape plugin CytoNCA was applied to assess the intersection from the following six centrality measures: betweenness centrality (BC), closeness centrality (CC), degree centrality (DC), eigenvector centrality (EC), network centrality (NC), and local average connectivity-based method (LAC) [19]. In addition, CytoHubba, a Cytoscape plugin, was used to select the hub genes [20].

2.5. GO and KEGG Pathway Enrichment Analyses

Gene ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed based on the public database GO (https://geneontology.org/) and KEGG (https://www.genome.jp/kegg/).

A hypergeometric distribution model was used to evaluate whether the target gene set was significantly associated with a specific gene ontology or biological pathway using the following formula:

Here, N represents the total number of genes from the reference terms, M represents the number of genes annotated to specific GO terms or pathways, n represents the number of identified Nux vomica-target genes, and k represents the number of common genes between Nux vomica-target genes and the reference set. The clusterProfiler package in R (4.0) (https://www.bioconductor.org/packages/release/bioc/html/clusterProfiler.html) was adopted to perform the GO and KEGG pathway enrichment analyses with P values <0.01, adjusted by the Bonferroni correction.

2.6. Molecular Docking Analysis

The zinc database (https://zinc15.docking.org) was used to download two-dimensional (2D) structures of bioactive compounds [21]. AutoDockTools-1.5.6 was used to add nonpolar hydrogen, calculate charges, and set rotatable keys for the structure and save it as a “pdbqt” file [22]. The 3D structures of the hub targets were downloaded from the Protein Data Bank (PDB) database (https://www.rcsb.org/) [23]. The hub targets and the original ligands were separated and the water molecules were removed using PyMol 2.3.0 [24]. Next, the ligand and receptor were imported into AutoDockTools-1.5.6 to add nonpolar hydrogen, calculate charges, and save as a “pdbqt” file. The protein proligand was set as the center of the docking box; if no proligand was present in the target protein, the docking region was identified based on key amino acid residues near the active pocket of the target, as reported in previous studies, or the fpocket [25] online pocket parameter prediction website for analysis. The prepared protein files, ligand files, compound files, and scripts were placed in the appropriate folders, and molecular docking was performed using AutoDock Vina 1.1.2 [26]. Finally, the binding affinity was calculated using AutoDock Vina 1.1.2, and the docking results were visualized using the open-source version PyMol 2.3.0.

2.7. Animal Experiments

Female C57BL/6J mice (aged 7–8 weeks) were purchased from the Animal Resources Center (Murdoch, WA, Australia). All mice were housed in a controlled environment under a 12-h light/dark cycle with 21°C and 45–65% humidity and had free access to food and water.

After anaesthesia, mice in the model group were emulsified with 40 μg MuSK in 100 μL phosphate-buffered saline (PBS) and 100 μL complete Freund’s adjuvant (CFA; Becton, Dickinson and Company, NJ, USA), which were intraperitoneally injected into the hind foot. The control group was injected with 100 μL of PBS and 100 μL of CFA. The mice in both groups were injected with cyclophosphamide monohydrate (300 mg/kg, dissolved in 0.9% saline to a final concentration of 10 mg/mL; Sigma-Aldrich, MO, USA) within 24 h to suppress immune resistance. The aforementioned operation was repeated after 30 days. The body weight of the mice was measured weekly, and a clinical examination was performed every 2 weeks to determine the model-building cycle. After the MG mouse model was successfully constructed, the mice were randomly divided into five groups (n = 8). The mice in the low-, mid-, and high-dose Nux vomica groups were administered low (4 mg/kg), mid (8 mg/kg), and high (16 mg/kg) doses of Nux vomica. The mice in the control and model groups were intragastrically administered equivalent doses of saline. The body weight of the mice was measured weekly, and a clinical examination was performed every 2 weeks. After 30 days, mice were sacrificed by cervical dislocation after injecting sodium pentobarbital (50 mg/kg), and tendon tissues were removed for subsequent experiments.

2.8. Clinical Evaluation

For clinical examination, mice were left for 3 min on a flat platform and observed for signs of MG. The clinical muscle weakness was graded as follows: Grade 0.5, the mouse exhibited slightly lower activity than normal after the standard exercise state, and the hanging time was >60 s; Grade 1, the mouse exhibited normal movement, lying prone, staying still, and/or hanging time <60 s; Grade 2, before and after exercise, the mouse was lying prone, staying still, and/or hanging test <30 s; Grade 3, severe weakness, dehydration and paralysis, weight loss >15%, and/or hanging time 0 s; Grade 4, the mouse died. The grip of the mice was observed every 2 weeks. First, the tail of each mouse was gently grabbed and the mouse was suspended above a metal grid. When the mouse grasped the metal grid, the mouse was lifted and the tension was measured 20 times. One end was fixed on the tension device, and the front end of the force measurement device was made into a grid shape. The tail of the mouse was grasped, the front end of the force measuring device was grasped, and the tail of the mouse was pulled forcefully. When the mouse was not strong enough to grasp the front end of the force measuring device, the force measuring device showed the grip force at this time. In addition, clinical progress was assessed using an inverted screening test every 2 weeks. In this test, the mouse was placed in the center of the grid and immediately rotated to the upside-down position, and stably maintained 50 cm above the soft cushion. The time when the mouse fell was recorded and the observation time was limited to 300 s.

2.9. Quantitative Real-Time PCR (qRT-PCR) Assay

Total RNA was extracted from the tendon tissues of the mice using TRIzol® Reagent (Ambion™) (Applied Biosystems, MA, USA), and 500 ng of RNA was reverse-transcribed into cDNA using TaqMan® Reverse Transcription Reagents (Thermo Fisher Scientific, MA, USA). The expression of target genes was examined using ABI Prism SDS 7000 (Thermo Fisher Scientific) at 95°C for 1 min, followed by 40 cycles of 95°C for 15 s, and 60°C for 1 min, and then calculated using the 2^−ΔΔCt method. Primer sequences for qRT-PCR are listed in Table S1. GAPDH was used as an internal control.

2.10. Western Blot Analysis

The predicted target protein expression was detected by western blotting. Briefly, total proteins were extracted from tendon tissue using RIPA lysis buffer (Applied Biosystems, MA, USA). Then, the proteins were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and transferred to a polyvinylidene fluoride membrane (Millipore, MA, USA). The membranes were blocked with 5% skimmed milk and incubated with primary antibodies anti-EGFR, anti-TP53, anti-AKT1, anti-MYC, and anti-GAPDH (1 : 1,000, Abcam, UK) at 4°C overnight. After washing with PBS, the membranes were incubated with a secondary antibody (1 : 1,000, Abcam). Enhanced chemiluminescent reagents (Millipore, MA, USA) were used to detect the protein bands. GAPDH was used as an internal control.

2.11. Statistical Analysis

Data are shown as the mean ± standard deviation (SD) and analyzed by one-way ANOVA followed by Tukey’s test using SPSS 27.0 software (IBM, IL, USA). Differences were considered significant at .

3. Results

3.1. Nux Vomica Bioactive Compounds-Targets Network

The TCMSP database was used to search for the bioactive compounds of Nux vomica with the criteria of OB ≥ 30% and DL ≥ 0.18. The TCMSP database was also used to search for the targets of Nux vomica. As a result, 11 active compounds with the targets were found in Nux vomica. Previous studies have shown that strychnine is one of the major bioactive and principal toxic compounds of Nux vomica [9]. Therefore, strychnine was considered as a candidate in this study. The DL and pharmacokinetics of strychnine were subsequently evaluated using the online tool SwissADME. The results indicated that strychnine exhibited good gastrointestinal absorption and high DL properties, confirming that the strain was an active compound of Nux vomica. The targets of strychnine were predicted by SwissTargetPrediction, and the top 25 were viewed as key targets of strychnine. Collectively, we found 12 bioactive compounds and 72 targets in Nux vomica, and the details are shown in Tables 1 and 2. Subsequently, we constructed a Nux vomica bioactive compound-target network using Cytoscape (Figure 1).

3.2. GO and KEGG Pathway Enrichment Analysis of Nux Vomica

To further investigate the 72 targets of the Nux vomica bioactive compounds, GO enrichment analysis was performed using biological processes (BP), molecular function (MF), and cellular component (CC) terms. A total of 103 MFs, 454 BPs, and 48 CCs were found to be enriched (). The GO analysis with the top 15 markedly enriched MF, BP, and CC terms is shown in Figures 2(a)–2(c). KEGG pathway enrichment analysis demonstrated that 72 targets of Nux vomica bioactive compounds were enriched in 17 associated pathways (). The top 15 KEGG pathways in which the targets of Nux vomica bioactive compounds were enriched are presented in Figure 2(d) and Table 3.

(a)

(b)

(c)

(d)

3.3. Target Genes of MG

Thirty-seven MG-related targets were obtained from the GENECARD database, and 42 were from the DisGeNET database. Seventy-nine disease target genes were retrieved after combining the results of both databases, and the duplicates were deleted. The details are listed in Table 4.

3.4. Nux Vomica-MG PPI Network

The PPI networks of Nux vomica- and MG-related targets were constructed using Bisogenet. The PPI network of Nux vomica-related genes was constructed with 2,275 nodes and 92,252 edges (Figure 3(a)). In addition, the MG-related target network was constructed with 3,480 nodes and 121,919 edges (Figure 3(b)). Both networks were subsequently merged into an intersection with 1,517 nodes and 50,040 edges (Figure 3(c)). CytoNCA was used to assess the intersection of the PPI network by topological analysis. A Nux vomica on the myasthenia PPI network was first screened based on the criteria of “BC ≥ 902.385, CC ≥ 0.471, DC ≥ 52” (Figure 3(d)). A core-target PPI network was further screened using the criteria of “BC ≥ 1985.705, CC ≥ 0.500, DC ≥ 90,” which consists of 204 nodes and 4,668 edges (Figure 3(e)).

Figure 3

Identification of core targets for Nux vomica against MG. (a) The PPI network of Nux vomica-related targets (2,275 nodes and 92,252 edges). (b) The PPI network of MG-related targets (3,480 nodes and 121,919 edges). (c) Intersection of PPI networks (1,517 nodes and 5,0040 edges). (d) PPI network by the screening criteria of BC ≥ 902.385, CC ≥ 0.471, DC ≥ 52 (567 nodes and 19,098 edges). (e) Core-target PPI network by the screening criteria of BC ≥ 1985.705, CC ≥ 0.500, DC ≥ 90 (204 nodes and 4,668 edges). BC, betweenness centrality; CC, closeness centrality; DC, degree centrality.

3.5. Enrichment Analysis of the Core-Target PPI Network

In the GO enrichment analysis, a total of 144 molecular functions, 1,570 biological processes, and 151 cell components were identified in 204 core targets of Nux vomica following MG treatment (). The top 15 significant GO terms, including MF, BP, and CC terms, are shown in Figures 4(a)–4(c). KEGG pathway enrichment analysis demonstrated that 132 pathways were enriched in 204 core targets (). The top 15 KEGG pathways enriched for the targets of bioactive compounds from Nux vomica are shown in Figure 4(d).

(a)

(b)

(c)

(d)

3.6. Enrichment Analysis of the Hub Genes

Hub genes were selected using CytoHubba based on the core-target PPI network. The results showed that TP53, EGFR, UBC, MYC, HSP90AA1, EP300, HSPA8, RPS27A, AKT1, and GAPDH were the major hub nodes in the core-target PPI network, and the details are shown in Table 5. The PPI network of the hub genes is shown in Figure 5.

The results of the most representative GO terms included regulation of protein modification by small protein conjugation or removal, cellular response to drugs, activation of the innate immune response, intrinsic apoptotic signaling pathway, regulation of protein ubiquitination, regulation of the innate immune response, positive regulation of cellular protein localization, regulation of protein stability, positive regulation of the innate immune response, regulation of the mRNA metabolic process, positive regulation of establishment of protein localization, regulation of binding, regulation of sequence-specific DNA binding transcription factor activity, regulation of apoptotic signaling pathway, and rhythmic process (Figure 6).

In addition, the top 15 pathways in which the core targets were enriched were identified (Figure 7), mainly including the neurotrophin signaling pathway, PI3K-Akt signaling pathway, cell cycle, apoptosis, viral carcinogenesis, prostate cancer, hepatitis C, pancreatic cancer, Epstein-Barr (EB) virus infection, human cytomegalovirus infection, proteoglycans in cancer, hepatitis B, chronic myeloid leukemia, Kaposi sarcoma-associated herpes, and microRNAs (miRNAs) in cancer.

3.7. Molecular Docking Analysis

The five bioactive compounds (strychnine, stigmasterol, (S)-stylopine, brucine-N-oxide, and isostrychnineN-oxide (II)) are shown in Table 6. The top ten hug targets (TP53, EGFR, UBC, MYC, HSP90AA1, EP300, HSPA8, RPS27A, AKT1, and GAPDH) were selected as the core compounds and core targets for molecular docking. AutoDock Vina assesses the binding of small molecules mainly by binding energy (affinity) to assess the binding strength of a small molecule to a protein. If the affinity is <0, the ligand can spontaneously bind to the receptor, and the smaller the value, the easier it is for the active ingredient to bind to the receptor. Therefore, the interaction between the active compound and its corresponding target was analyzed by the binding energy, and a threshold affinity <–5 kcal/mol was set in this study. A total of 50 docking pairs were obtained, as shown in Table 7. Figure 8 shows the binding patterns of the top ten binding core targets and their corresponding compounds.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 8

Molecular docking of hub genes with active compounds. (a) The binding poses of strychnine complexed with EGFR. (b) The binding poses of Stigmasterol complexed with EP300. (c) The binding poses of (S)-Stylopine complexed with EP300. (d) The binding poses of (S)-Stylopine complexed with EGFR. (e) The binding poses of Stigmasterol complexed with EGFR. (f) The binding poses of brucine-N-oxide complexed with EGFR. (g) The binding poses of strychnine complexed with EP300. (h) The binding poses of (S)-Stylopine complexed with TP53. (i) The binding poses of IsostrychnineN-oxide (II) complexed with EGFR. (j) The binding poses of Stigmasterol complexed with GAPDH.

3.8. MG Mouse Model Establishment

As shown in Figure 9, the body weight of the mice increased every week, and after 42 days, the body weight of the mice in the model group was significantly lower than that of the mice in the control group (). Clinical examinations were performed every 2 weeks. After 42 days, compared to the control group, the clinical scores were increased, the grip strengths were weakened, and the inverted screen times were shortened in the model group ().

(a)

(b)

(c)

(d)

3.9. Nux Vomica Improves Clinical Examination Indexes of MG Model Mice

After treatment with Nux vomica, the clinical examination indicators of the MG model (weight, clinical score, grip strength, and inverted screen time) significantly improved. After 28 days of treatment, the weight of the MG mice increased significantly, clinical scores reduced, grip strength increased, and the inverted screen times were increased with the increased dose of Nux vomica (Figure 10).

(a)

(b)

(c)

(d)

3.10. Expression of Hub Targets

To further verify the effect of Nux vomica against MG, we detected the expression of hub targets by qRT-PCR and western blotting. The results showed that compared to the model group, Nux vomica treatment significantly reduced EGFR and TP53 expression and increased MYC expression in a dose-dependent manner. Moreover, HSPA8 expression was not significantly changed, while AKT1 expression was significantly reduced, and the effect was particularly pronounced in the high-dose group. The qRT-PCR results also showed that the mRNA expression of RPS27, UBC, HSP90AA1, and EP300 was not significantly different between the groups (Figure 11).

(a)

(b)

4. Discussion

MG is an autoimmune disease of the neuromuscular junction, which is mainly mediated by acetylcholine receptor (AChR) antibody and muscle-specific kinase (MuSK) [27, 28]. Presently, traditional therapies for long-term MG treatment have several adverse effects in some populations. Therefore, it is vital to explore new methods for treating MG. Owing to their high efficiency and few side effects, TCMs have been widely used in MG treatment as an alternative or adjuvant therapy. In addition, a previous meta-analysis indicated that compared to using Western medication alone, TCMs could function as adjuvant therapies to improve the efficacy and decrease the side effects and relapse rate in treating MG [29, 30]. As a TCM, Nux vomica is warm in nature, with a slightly bitter flavor due to its various pharmacological activities. However, the effects and mechanism of action of Nux vomica on MG treatment remain unclear. Therefore, our study systematically explored the underlying mechanism of Nux vomica treatment for MG based on network pharmacology and further validated the therapeutic effect of Nux vomica on MG in combination with in vivo experiments.

Based on the public databases, 12 bioactive compounds of Nux vomica were collected by setting OB ≥ 30% and DL ≥ 0.18 (Table 1). Studies have isolated many chemical compounds from Nux vomica, including alkaloids, glycosides, steroids, and organic acids [9]. Among them, alkaloids, including strychnine, brucine-N-oxide, isobrucine, vomicine, and brucine N-oxide, are the primary chemical components of Nux vomica and have been reported to have anticancer effects. Additionally, strychnine has been shown to exhibit significant anti-inflammatory activity [31], while brucine-N-oxide possesses a strong antinociceptive potential [32]. Moreover, stigmasterol can inhibit the activity of acetylcholinesterase (AChE) and has antioxidant, anti-inflammatory, and neuroprotective effects [33–35]. Thus, many compounds from Nux vomica possess favorable biological activity, including inhibiting AChE and immunomodulation and may exert a mitigative effect on MG.

MG occurrence is closely related to continuous viral infections, genetic factors, and abnormal immune responses. Growing evidence has demonstrated that EBV infection is correlated with the occurrence of MG [36, 37]. For example, EBV infection could affect TLR7/9-dependent innate immune responses, which enhance the inflammatory response, B-cell dysfunction, and autoimmunity in the thymus of patients with MG [36, 37]. Moreover, a recent pilot study showed that hepatitis E infection was relevant to MG, at least in China, and may contribute to MG by inducing replication in the thymus, bystander activation, and stimulation of antigen-specific T cells [38]. Target identification is a critical step for the mechanism investigation on Nux vomica treating MG. According to the PPI network of active compound targets, Nux vomica has the characteristics of multiple components and targets in the treatment of MG. Top 10 potential targets (EP300, EGFR, TP53, RPS27A, MYC, HSP90AA1, UBC, AKT1, GAPDH, HSPA8) of active compounds from Nux vomica were predicted to be associated with MG treatment. These 10 hub targets were validated in vivo, and results showed that Nux vomica treatment downregulated the expression of EGFR, TP53, and AKT1, while upregulated the MYC level in MG model mice. These findings indicated that EGFR, TP53, AKT1, and MYC may be the key targets for Nux vomica treating MG. In addition, the pathways of Nux vomica against MG were related to EBV infection, hepatitis B, hepatitis C, and other viral infections based on KEGG analysis of the core genes. We also identified the miRNA pathway in cancer. An increasing number of studies have suggested that miRNAs are involved in the occurrence and development of MG [39]. The miRNA profile of peripheral blood from MG patients was sequenced using a next-generation sequencing approach, and 41 miRNAs were found to be differentially expressed in nonresponders versus responders to immunosuppressive therapies. More importantly, miR-181d-5p, -323b-3p, -340-3p, -409-3p, and -485-3p can modulate the immune response and drug metabolism processes in patients with MG [40]. In another study, miR-612, miR-3651, and miR-3653 were found to be highly expressed in the peripheral blood mononuclear cells of patients with AChR-positive early-onset MG (AChR-EOMG), suggesting that these miRNAs may be involved in the pathogenesis of AChR-EOMG [41]. It was predicted that the active ingredients of Nux vomica may regulate the immune response by regulating the miRNA signaling pathway, thereby playing a role in the treatment of MG.

5. Conclusions

Our study systematically analyzed the potential effects of Nux vomica in the treatment of MG based on network pharmacology and validated its effects in combination with in vivo experiments. The findings showed that the bioactive compounds of Nux vomica exerted therapeutic effects on MG by acting on multiple targets and multiple pathways. Based on the network pharmacology, this study confirmed the therapeutic effects of Nux vomica on MG probably via downregulating the expression of EGFR, TP53, and AKT1 and upregulating MYC level. Our data may provide a valuable theoretical basis for future MG diagnosis and treatment and also serve as a reference for further mechanistic research on Nux vomica treating MG. However, we also found that the mechanisms of Nux vomica against are involved in the PI3K-Akt signaling pathway and miRNAs, so further experiments are needed to prove these predictions. Besides, the therapeutic efficacy of Nux vomica on MG is needed to be verified in a clinical trial.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Medical and Health Research Projects of Zhejiang Province (number 2019RC058).

Supplementary Materials

Table S1. The primer sequences for qRT-PCR. (Supplementary Materials)

References

M. K. Hehir and N. J. Silvestri, “Generalized myasthenia gravis,” Neurologic Clinics, vol. 36, no. 2, pp. 253–260, 2018.
View at: Publisher Site | Google Scholar
E. Ciafaloni, “Myasthenia gravis and congenital myasthenic syndromes,” Continuum: Lifelong Learning in Neurology, vol. 25, no. 6, pp. 1767–1784, 2019.
View at: Publisher Site | Google Scholar
R. Mantegazza and C. Antozzi, “From traditional to targeted immunotherapy in myasthenia gravis: prospects for research,” Frontiers in Neurology, vol. 11, p. 981, 2020.
View at: Publisher Site | Google Scholar
V. Voinov, Plasmapheresis in Treatment of Myasthenia Gravis, IntecOpen, London, UK, 2018.
X. Jiang, G. Chen, J. Huang et al., “Modified Buzhong Yiqi decoction for myasthenia gravis,” Medicine (Baltimore), vol. 97, no. 50, Article ID e13677, 2018.
View at: Publisher Site | Google Scholar
G. Qi, P. Liu, S. Gu, H. Yang, H. Dong, and Y. Xue, “Effect of Jianpiyiqi granule on ocular myasthenia gravis in children prepared with a formula from traditional Chinese medicine,” Journal of Traditional Chinese Medicine, vol. 37, pp. 546–553, 2017.
View at: Google Scholar
W. Jiao, F. Hu, J. Li et al., “Qiangji Jianli decoction promotes mitochondrial biogenesis in skeletal muscle of myasthenia gravis rats via AMPK/PGC-1α signaling pathway,” Biomedicine & Pharmacotherapy, vol. 129, Article ID 110482, 2020.
View at: Publisher Site | Google Scholar
M. Zhu, Z. Liu, M. Gao et al., “The effect of Bu Zhong Yi Qi decoction on simulated weightlessness-induced muscle atrophy and its mechanisms,” Molecular Medicine Reports, vol. 16, no. 4, pp. 5165–5174, 2017.
View at: Publisher Site | Google Scholar
R. Guo, T. Wang, G. Zhou et al., “Botany, phytochemistry, pharmacology and toxicity of strychnos nux-vomica L.: a review,” The American Journal of Chinese Medicine, vol. 46, pp. 1–23, 2018.
View at: Publisher Site | Google Scholar
X. H. Jiang, Y. Chen, Y. Y. Ding, H. Qiu, D. Y. Zhou, and C. L. Qiu, “Effect of grilled nux vomica on differential RNA expression profile of gastrocnemius muscle and toll-like receptor 4 (tlr-4)/nuclear factor kappa b (nf-κb) signaling in experimental autoimmune myasthenia gravis rats,” Medical Science Monitor, vol. 26, Article ID e919150, 2020.
View at: Publisher Site | Google Scholar
A. L. Hopkins, “Network pharmacology: the next paradigm in drug discovery,” Nature Chemical Biology, vol. 4, no. 11, pp. 682–690, 2008.
View at: Publisher Site | Google Scholar
W. Y. Lee, C. Y. Lee, Y. S. Kim, and C. E. Kim, “The methodological trends of traditional herbal medicine employing network pharmacology,” Biomolecules, vol. 9, 2019.
View at: Publisher Site | Google Scholar
B. Liu, Z. Song, J. Yu, P. Li, Y. Tang, and J. Ge, “The atherosclerosis-ameliorating effects and molecular mechanisms of BuYangHuanWu decoction,” Biomedicine & Pharmacotherapy, vol. 123, Article ID 109664, 2020.
View at: Publisher Site | Google Scholar
N. Rappaport, M. Twik, I. Plaschkes et al., “Malacards: an amalgamated human disease compendium with diverse clinical and genetic annotation and structured search,” Nucleic Acids Research, vol. 45, pp. D877–D887, 2017.
View at: Publisher Site | Google Scholar
J. Piñero, À. Bravo, N. Queralt-Rosinach et al., “Disgenet: a comprehensive platform integrating information on human disease-associated genes and variants,” Nucleic Acids Research, vol. 45, pp. D833–D839, 2017.
View at: Publisher Site | Google Scholar
J. Zhang, Q. Huang, R. Zhao, and Z. Ma, “A network pharmacology study on the Tripteryguim wilfordii hook for treatment of Crohn’s disease,” BMC Complementary Medicine and Therapies, vol. 20, no. 1, p. 95, 2020.
View at: Publisher Site | Google Scholar
J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, p. 13, 2014.
View at: Publisher Site | Google Scholar
A. Martin, M. E. Ochagavia, L. C. Rabasa, J. Miranda, J. Fernandez-de-Cossio, and R. Bringas, “BisoGenet: a new tool for gene network building, visualization and analysis,” BMC Bioinformatics, vol. 11, no. 1, p. 91, 2010.
View at: Publisher Site | Google Scholar
Y. Tang, M. Li, J. Wang, Y. Pan, and F.-X. Wu, “CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks,” Biosystems, vol. 127, pp. 67–72, 2015.
View at: Publisher Site | Google Scholar
C. H. Chin, S. H. Chen, H. H. Wu, C. W. Ho, M. T. Ko, and C. Y. Lin, “Cytohubba: identifying hub objects and sub-networks from complex interactome,” BMC Systems Biology, vol. 8, no. 4, p. S11, 2014.
View at: Publisher Site | Google Scholar
T. Sterling and J. J. Irwin, “ZINC 15-ligand discovery for everyone,” Journal of Chemical Information and Modeling, vol. 55, no. 11, pp. 2324–2337, 2015.
View at: Publisher Site | Google Scholar
G. M. Morris, R. Huey, W. Lindstrom et al., “AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility,” Journal of Computational Chemistry, vol. 30, no. 16, pp. 2785–2791, 2009.
View at: Publisher Site | Google Scholar
S. K. Burley, H. M. Berman, G. J. Kleywegt, J. L. Markley, H. Nakamura, and S. Velankar, “Protein data bank (PDB): the single global macromolecular structure archive,” Methods in Molecular Biology, vol. 1607, pp. 627–641, 2017.
View at: Publisher Site | Google Scholar
D. Seeliger and B. L. de Groot, “Ligand docking and binding site analysis with PyMOL and Autodock/Vina,” Journal of Computer-Aided Molecular Design, vol. 24, no. 5, pp. 417–422, 2010.
View at: Publisher Site | Google Scholar
P. Schmidtke, V. Le Guilloux, J. Maupetit, and P. Tuffery, “fpocket: online tools for protein ensemble pocket detection and tracking,” Nucleic Acids Research, vol. 38, pp. W582–W589, 2010.
View at: Publisher Site | Google Scholar
O. Trott and A. J. Olson, “AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading,” Journal of Computational Chemistry, vol. 31, pp. 455–461, 2010.
View at: Publisher Site | Google Scholar
N. E. Gilhus and J. J. Verschuuren, “Myasthenia gravis: subgroup classification and therapeutic strategies,” The Lancet Neurology, vol. 14, no. 10, pp. 1023–1036, 2015.
View at: Publisher Site | Google Scholar
D. B. Sanders, G. I. Wolfe, M. Benatar et al., “International consensus guidance for management of myasthenia gravis,” Neurology, vol. 87, no. 4, pp. 419–425, 2016.
View at: Publisher Site | Google Scholar
S. Chen, M.-B. Xu, X.-L. Zhou, P.-Q. Rong, T.-Y. Jin, and G.-Q. Zheng, “Chinese herbal medicine for myasthenia gravis: a systematic review and meta-analysis,” Frontiers in Pharmacology, vol. 9, p. 969, 2018.
View at: Publisher Site | Google Scholar
R. Xie, L. Liu, R. Wang, and C. Huang, “Traditional Chinese medicine for myasthenia gravis,” Medicine (Baltimore), vol. 99, no. 29, Article ID e21294, 2020.
View at: Publisher Site | Google Scholar
J. Chen, X. Wang, Y.-G. Qu et al., “Analgesic and anti-inflammatory activity and pharmacokinetics of alkaloids from seeds of Strychnos nux-vomica after transdermal administration: effect of changes in alkaloid composition,” Journal of Ethnopharmacology, vol. 139, no. 1, pp. 181–188, 2012.
View at: Publisher Site | Google Scholar
Y. Wu, T.-S. Wang, F.-Z. Yin, and B.-C. Cai, “Analgesic and anti-inflammatory properties of brucine and brucine n-oxide extracted from seeds of Strychnos nux-vomica,” Journal of Ethnopharmacology, vol. 88, pp. 205–214, 2003.
View at: Google Scholar
S. Panda, M. Jafri, A. Kar, and B. K. Meheta, “Thyroid inhibitory, antiperoxidative and hypoglycemic effects of stigmasterol isolated from butea monosperma,” Fitoterapia, vol. 80, no. 2, pp. 123–126, 2009.
View at: Publisher Site | Google Scholar
C. I. B. Walker, S. M. Oliveira, R. Tonello et al., “Anti-nociceptive effect of stigmasterol in mouse models of acute and chronic pain,” Naunyn-Schmiedeberg’s Archives of Pharmacology, vol. 390, no. 11, pp. 1163–1172, 2017.
View at: Publisher Site | Google Scholar
J. Sun, X. Li, J. Liu, X. Pan, and Q. Zhao, “Stigmasterol exerts neuro-protective effect against ischemic/reperfusion injury through reduction of oxidative stress and inactivation of autophagy,” Neuropsychiatric Disease and Treatment, vol. 15, pp. 2991–3001, 2019.
View at: Publisher Site | Google Scholar
P. Cavalcante, B. Galbardi, S. Franzi et al., “Increased expression of Toll-like receptors 7 and 9 in myasthenia gravis thymus characterized by active Epstein–Barr virus infection,” Immunobiology, vol. 221, no. 4, pp. 516–527, 2016.
View at: Publisher Site | Google Scholar
P. Cavalcante, C. Barzago, F. Baggi et al., “Toll-like receptors 7 and 9 in myasthenia gravis thymus: amplifiers of autoimmunity?” Annals of the New York Academy of Sciences, vol. 1413, no. 1, pp. 11–24, 2018.
View at: Publisher Site | Google Scholar
L. Wang, F. Gao, G. Lin et al., “Association of hepatitis e virus infection and myasthenia gravis: a pilot study,” Journal of Hepatology, vol. 68, no. 6, pp. 1318–1320, 2018.
View at: Publisher Site | Google Scholar
L. Wang and L. Zhang, “Emerging roles of dysregulated micrornas in myasthenia gravis,” Frontiers in Neuroscience, vol. 14, p. 507, 2020.
View at: Publisher Site | Google Scholar
P. Cavalcante, T. Mizrachi, C. Barzago et al., “Microrna signature associated with treatment response in myasthenia gravis: a2q further step towards precision medicine,” Pharmacological Research, vol. 148, Article ID 104388, p. 104388, 2019.
View at: Publisher Site | Google Scholar
C. Barzago, J. Lum, P. Cavalcante et al., “A novel infection- and inflammation-associated molecular signature in peripheral blood of myasthenia gravis patients,” Immunobiology, vol. 221, no. 11, pp. 1227–1236, 2016.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Chao Qiu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

303

Downloads

451

Citations