- Open Access
The functional mechanisms of synchronizing royal jelly consumption and physical activity on rat with multiple sclerosis-like behaviors hallmarks based on bioinformatics analysis, and experimental survey
BMC Neuroscience volume 23, Article number: 34 (2022)
Natural nutrition and physical training have been defined as non-pharmacochemical complementary and alternative medicines to prevent and treat various pathogenesis. Royal jelly possesses various pharmacological properties and is an effective therapeutic supplement for halting neurodegeneration. Multiple sclerosis is a prevalent neurodegenerative disorder that manifests as a progressive neurological condition. Inflammation, hypoxia, and oxidative stress have been identified as significant hallmarks of multiple sclerosis pathology.
In the present study, based on artificial intelligence and bioinformatics algorithms, we marked hub genes, molecular signaling pathways, and molecular regulators such as non-coding RNAs involved in multiple sclerosis. Also, microRNAs as regulators can affect gene expression in many processes. Numerous pathomechanisms, including immunodeficiency, hypoxia, oxidative stress, neuroinflammation, and mitochondrial dysfunction, can play a significant role in the MSc pathogenesis that results in demyelination. Furthermore, we computed the binding affinity of bioactive compounds presented in Royal Jelly on macromolecules surfaces. Also, we predicted the alignment score of bioactive compounds over the pharmacophore model of candidate protein as a novel therapeutic approach. Based on the q-RT-PCR analysis, the expression of the Dnajb1/Dnajb1/Foxp1/Tnfsf14 and Hspa4 networks as well as miR-34a-5p and miR155-3p were regulated by the interaction of exercise training and 100 mg/kg Royal Jelly (ET-100RJ). Interestingly, characteristics, motor function, a proinflammatory cytokine, and demyelination were ameliorated by ET-100RJ.
Here, we indicated that interaction between exercise training and 100 mg/kg Royal jelly had a more effect on regulating the microRNA profiles and hub genes in rats with Multiple sclerosis.
Neurodegenerative diseases are progressive conditions that occur through impairment in the function or structure of neurons in the peripheral nervous system or brain . Moreover, This condition affects millions of individuals worldwide, and progressive degenerative diseases can influence body activities and functions such as breathing, moving, talking, and balance. For example, multiple sclerosis (MSc) is the prevalent neurodegenerative disorder caused by increasing neuroinflammation, which stimulates the immune cells and destroys neurons, axons, and myelin-producing oligodendrocytes .
Although several significant factors, including genetics, epigenetics, and environment, are associated with multiple sclerosis [3, 4], the pinpoint mechanism of MSc is not yet comprehended. However, several pathomechanisms such as immunodeficiency, hypoxia, oxidative stress, neuroinflammation, and mitochondrial dysfunction have been strongly involved in the MSc pathogenesis that leads to demyelination. Since increasing inflammation, hypoxia, and oxidative stress as dominant hallmarks are identified in multiple sclerosis pathology. Therefore, myelin Oligodendrocyte Glycoprotein (MOG) is widely applied as a model for demyelination. The dominant mechanism of MOG’s function is involved in inducing inflammation and oxidative stress, which leads to demyelination in the nervous system . On the other hand, the artificial intelligence survey and protein–protein interaction network of data’s multiple sclerosis highlighted many principal hub genes involved in multiple sclerosis incidence and marked candidate proteins for the design and discovery of the suitable drug.
Hence, based on previous evidence around multiple sclerosis pathogenesis, the Heat Shock Protein Families such as HSP40 and HSP70 are defined as potentially protein families involved in the demyelination process . Member B1 and B2 from the heat shock protein family with 40 kDa (HSP40) are termed DNAJB1and DNAJB2. Furthermore, HSPA4 from the heat shock protein family (70 kDa) is associated with demyelination in the peripheral nervous system and multiple sclerosis incidence . Moreover, large-scale genome-wide association studies (GWAS) showed forkhead box protein P1 loci (foxp1) as a susceptible genomic location in multiple sclerosis have found differential expression patterns . In addition, one of the crucial molecules involved in proliferation, differentiation, maturation, activation, and migration of the central nervous system is the tumor necrosis factor superfamily (TNFSF) . Also, the clinical study revealed a strong association between autoimmune disorders such as multiple sclerosis and TNFSF .
In addition to suscept genes in multiple sclerosis pathogenesis, other agents such as environmental factors, epigenetic factors, and regulatory factors are associated with multiple sclerosis pathogenesis [3, 4, 9]. The genome’s non-coding sequences are usually involved in the post-transcriptional regulation of genes and gene expression. The microRNAs usually bind to the 3´UTR of coding genes and lead to inhibition or suppression of target genes. On the other hand, several studies suggested that the binding of microRNAs in 5´UTR of genes could consider the positive regulation of gene expression, leading to the induction of transcription .
Notably, there has been no practical and effective treatment strategy to cure or halt the development of MSc recently, and management of MSc remains complex. Therefore, numerous attempts have been conducted to understand the molecular mechanisms responsible for MSc to progress new therapeutic approaches for this condition [11, 12].
As unique and exceptional food for bee queen larvae nutrition, Royal jelly is secreted from the hypopharyngeal glands of the worker bees [13, 14]. The Royal jelly with complex and variable bioactive compounds relevant to seasonal and geographical area conditions is influenced by growth, development, reproduction, and lifespan . Although the content of compounds is variable, Royal jelly includes proteins, lipids such as free fatty acids, sugar, various vitamins, essential amino acids, and water. Based on the data mining, we found that 10-hydroxy-2-decenoic acid is the principal compound and plays a crucial role in several biological processes such as inflammation and oxidative stress, estrogen-like function, antibacterial, and modulating the immune system . Hence, evidence proposed that Royal jelly with various pharmacological activities is identified as an efficient supplement in therapeutic medication.
Furthermore, based on clinical evidence, physical activity can ameliorate the physical fitness features, including mobility, fatigue, and central nervous system-related outcomes white matter integrity and cognition in people with MSc. Evidence has indicated that physical activity in rodents with experimental autoimmune encephalomyelitis (EAE) significantly decreased demyelination in the lumbar spinal cord and immune cell infiltration . In addition, exercise may consider a preventive therapy for developing the risk of the MSc . Moreover, exercise can delay disease onset and preserve the motor neurons in the lumbar spinal cord and axons, which reduces the severity of disability .
Based on literature review and artificial intelligence, the association of miR-155-3p and miR-34a-p with DNAJB1, DNAJB2, FOXP1, TNFSF14, and HSPA4 manifested. On the other hand, we found that miR-155-3p and miR-34a-p have been changed in multiple sclerosis compared to healthy individuals. Hence, we proposed that Royal jelly consumption and regular exercise might modulate these biological markers in multiple sclerosis and improve conditions through neuroinflammation reduction and the regeneration of oligodendrocytes that reprogram myelin.
Based on artificial intelligence surveys and analysis of microarray dataset in spinal-brain and blood samples of multiple sclerosis patients compared to healthy individuals, we detected 379 genes with significantly differential expression in blood samples which 293 genes are overexpressed and 86 genes are downregulated. We showed differential genes expression with P.value < 0.001 by heatmap diagram (Fig. 1A). after applying network visualize parameters, we obtained 70 hub genes. We have drawn a protein–protein interactions network construction of 70 hub genes with associated functions in different clusters (Fig. 1B). Moreover, analysis of spinal-brain samples of the multiple sclerosis dataset highlighted 380 genes with differential expression patterns compared to healthy individuals. Among 380 genes, we identified 218 genes with negative expression patterns and 162 genes with positive expression patterns. Outcomes of gene expression pattern analysis in multiple sclerosis samples were shown in the heatmap diagram considering P.value < 0.001 (Fig. 1C). The protein–protein interactions network and classification in clusters based on the function of hub genes have been drawn by Gephi software after applying the network’s visualized parameters, including betweenness centrality, degree, and closeness (Fig. 1D). Enrichment of hub genes marked several molecular signaling pathways in demyelination as pathogenesis hallmark in multiple sclerosis based on browsing enrichment databases. We found that inflammation elements and pathways, extracellular matrix organization, adipogenesis, TGF-B signaling and activation, JAK/STAT signaling, and TNF/stress oxidative pathway with principal roles are involved in multiple sclerosis development. We have shown highlighted results of enrichment in Fig. 2A–E. To determine common hub genes between the spinal-brain, blood, and hub genes obtained from these datasets, we designed a Venn graph and highlighted five common genes, including DNAJB1 and DNAJB2 FOXP1, TNFSF14, and HSPA4 in these genes lists (Fig. 3A). Browsing common hub genes in several databases such as KEGG, Panther, Reactome, WikiPathway, Biocarta, and NCI Nurture marked vital signaling pathways associated with hub genes involved in multiple sclerosis pathogenesis. Apoptosis and inflammation signaling pathway, hypoxia, oxidative stress, Wnt/b-catenin signaling pathway, P53 pathway, TGF-B signaling pathway, and IGF/insulin pathway are the most critical pathways in multiple sclerosis incidence (Fig. 3B–G).
We designed a microRNAs: mRNA interaction network based on selected genes in the present study. The microRNA-155-3p interacted with DNAJB2, FOXP1, and HSPA4. On the other hand, microRNA-34a-5p was associated with FOXP1, TNFSF14, HSPA4, and DNAJB2 genes (Fig. 4A).
Analysis of microRNAs array in multiple sclerosis revealed that miR-34a-5p and miR-155-3p were up-regulated in multiple sclerosis conditions. Moreover, the prediction of microRNAs that targeted selected genes highlighted the role of miR-34a-5p and miR155-3p in multiple sclerosis pathogenesis. Enrichment of selected microRNAs based on function in the miRNet algorithm manifested that these microRNAs were associated with cellular and molecular processes such as DNA damage repair, wound healing, vascular remodeling, dendritic cell differentiation, extracellular matrix remodeling, response to cytokines, inflammation, apoptosis, response to hypoxia, cell death, and aging processes.
We designed a heatmap diagram to illustrate microRNAs-molecular and cellular signaling pathway interactions based on miRPath V.3 servers (Fig. 4B).
Drug design results
Based on outcomes of computing binding affinity between bioactive compounds of Royal jelly and macromolecules involved in multiple sclerosis network with the highest degree and betweenness centrality by molecular docking practice, we concluded that all of the ten bioactive compounds of Royal jelly in this bioinformatics study bind to DNAJB1 and TNFSF14 with suitable binding affinity and acceptable RMSD threshold. Outcomes of molecular docking are presented in Table 1. Furthermore, we designed a pharmacophore model with the DDDRR template and survival score 5.665015 based on 18 active compounds documented of HSPA4 in the binding database, considering 80% of all of the library and 4 to 5 feature selection (Fig. 5A). The screening of Royal jelly’s bioactive compounds library over the pharmacophore model specified that Kaempferol (align score: 0.846716) and Vitamin B2 (align score: 0.713298) matched with pharmacophore models and could be proposed as the new active ligands for the management of disorders by targeting the HSP70 family (Fig. 5B, C).
The characteristic and behavior of the multiple sclerosis model of Rat
The gain weight of the multiple sclerosis model of rats (M.Sc group) significantly decreased compared with the control group (Fig. 6A). Moreover, the gain weight of the M.Sc group sharply reduced in the second week after induction of multiple sclerosis Vs the Control group (Fig. 6A). In addition, this study evaluated the balance and motor coordination with the rotarod and vertical pole tests (Fig. 6B, C). The data demonstrated that the motor coordination and balance decreased in the M.Sc group compared with the control group (Fig. 6B, C). In addition, the serum concentration of the proinflammatory cytokine Interferon-gamma has indicated in Fig. 6D. In addition, the results revealed that the level of concentration of IFNγ was enhanced in the M.Sc group compared with the control group (Fig. 6D). Furthermore, in the M.Sc group, the pathological analysis indicated that the number of damaged cells of the nucleus significantly increased compared with the control group, and also, the shape of the cell was changed in the MSc group (Fig. 6E). Based on these data, the mobility, gain of weight, macrophage activation, and demyelination of the multiple sclerosis model of rats were altered compared with the control group.
Physical activity and Royal jelly compound ameliorate the motor function, proinflammatory cytokine, and demyelination
Our data indicated that the gain weight significantly was improved in multiple sclerosis rats models were treated with 50 mg/kg (RJ-50 group) and 100 mg/kg Royal jelly (RJ-100 group) in comparison M.Sc group (Fig. 7A). Moreover, physical activity (ET group) ameliorated the gain weight compared with the M.Sc group (Fig. 7A). Notably, the interaction of the exercise training and consumption of Royal jelly (ET-RJ50, ET-RJ100 groups) predominantly improved the gain weight of the multiple sclerosis rats model compared with other groups (Fig. 7A). Furthermore, data demonstrated that balance and motor coordination multiple sclerosis rat models were significantly enhanced in both 50 mg/kg (RJ-50 group) and 100 mg/kg Royal jelly (RJ-100 group) compared to the Control group. Also, Multiple sclerosis rats were treated with exercise training and 50 mg/kg Royal jelly (ET-RJ50 group), and Multiple sclerosis rats were treated with exercise training, and 100 mg/kg Royal jelly (ET-RJ100 group) were increased balance and motor coordination (Fig. 7B, C). Moreover, balance and motor coordination of sclerosis rats models were significantly increased in the ET-RJ100 group compared with other groups (Fig. 7B, C).
The macrophage activation factor and immune interferon indicated that the Interferon-gamma (IFNγ) level was reduced by consumption of the Royal jelly and physical activity. In addition, the highest reduction in the serum IFNγ was accomplished when a combination of exercise training was conducted with 100 mg/kg Royal jelly (ET-RJ100 group). As depicted, our data indicated that the concentration of the IFNγ was decreased in the RJ-50 group, RJ-100 group, and ET group compared with the M.Sc group (Fig. 7D). Notably, the level of IFNγ was not significant between the ET and RJ-100 groups. Based on these data, combining the exercise training with 100 mg/kg Royal jelly (ET-RJ100 group) effectively improved multiple sclerosis rat models.
The effect of multiple sclerosis condition on the alternation of Dnajb1/Dnajb2/Foxp1/Tnfsf14 and Hspa4 network in the brain tissue
Our data revealed that the relative expression of the Dnajb1/Dnajb2/Foxp1/Tnfsf14 was increased in the M.Sc group compared with the Control group (Fig. 8A–D). Moreover, we found that the relative expression of the Hspa4 was significantly decreased in the Multiple sclerosis rats model compared with the Control group (Fig. 8E). Based on these results and artificial intelligence analyses, the Apoptosis and inflammation signaling pathway, hypoxia, oxidative stress, Wnt/b-catenin signaling pathway, P53 pathway, TGF-B signaling pathway, IGF/insulin pathway are dysregulated in multiple sclerosis.
Physical activity and Royal jelly compound altered the Dnajb1/Dnajb1/Foxp1/Tnfsf14 and Hspa4 network in the brain tissue
We selected a pool mRNA network based on artificial intelligence and in-silico analysis. The maximum reduction of Dnajb1/Dnajb2/Foxp1/Tnfsf14 network was observed when Multiple sclerosis rats were treated with exercise training and 100 mg/kg Royal jelly (ET-RJ100), but the relative expression of the Hspa4 significantly increased (Fig. 8A–E). Moreover, exercise training (ET group), 50 mg/kg Royal jelly (RJ-50), and 100 mg/kg Royal jelly (RJ-100) reduced the relative expression of the Dnajb1/Dnajb2/Foxp1/Tnfsf14 network and increased the expression level of the Hspa4 compared with the M.Sc group (Fig. 8A–E). Also, consumption of 100 mg/kg Royal jelly (RJ-100) significantly decreased the Dnajb1/Dnajb1/Foxp1/Tnfsf14 network and enhanced Hspa4 in comparison to RJ-50 and ET groups (Fig. 8A–E).
The effect of Multiple sclerosis condition on alternation microRNAs profiles in the brain and blood
This study identified the critical microRNAs that could regulate Dnajb1/Dnajb1/Foxp1/Tnfsf14 and Hspa4 networks based on the in-silico analysis. Moreover, the relative expression of miR-34a-5p and miR155-3p was detected in blood and brain rats. Furthermore, we indicated that the expression level of the miR-34a-5p and miR155-3p significantly enhanced in the M.Sc group Vs. Control group (Fig. 8F–I). Therefore, these microRNAs profiles might consider biomarkers for prognosis and diagnosis development of Multiple Sclerosis.
Physical activity and Royal jelly compound altered miR-34a-5p and miR155-3p in the brain and blood
Interestingly, the level expression of miR-34a-5p and miR155-3p was changed in the brain and blood by physical activity and consumption of the Royal jelly. Our data revealed that the expression level of the miR-34a-5p and miR155-3p significantly was decreased in Multiple sclerosis rats treated with exercise training and 100 mg/kg Royal jelly (ET-RJ100 group) compared with the other groups (Fig. 8F–I). Furthermore, Multiple sclerosis rats were treated with 50 mg/kg Royal jelly (RJ-50 group), 100 mg/kg Royal jelly (RJ-100 group), and exercise training (ET group) reduced the relative expression of the miR-34a-5p and miR155-3p compared with the control groups (Fig. 8F–I). Moreover, the interaction of the exercise training with 50 mg/kg Royal jelly and 100 mg/kg Royal jelly (ET-RJ50, ET-RJ100 groups) diminished the level of these microRNAs profiles, but the ET-RJ100 group compared with the ET-RJ50 group significantly decreased the relative expression of the miR-34a-5p and miR155-3p. Hence, we could conclude that the combination of the 100 mg/kg Royal jelly along with physical activity had most effective in comparison to other groups (Fig. 8F–I).
The effect of exercise along with the Royal jelly compound on the neurons
Based on the pathological analysis, we revealed that the Royal jelly compound could reduce the number of damaged cells of the nucleus compared with the M.Sc group, and also the shape of the cell was improved in the RJ-50 group, RJ-100 group compared with the M.Sc group (Fig. 9). Moreover, we demonstrated that physical activity (ET group) could decrease the number and shape of damaged cells of the nucleus. Notably, exercise and the Royal jelly compound significantly improved the number and shape of damaged cells and reduced demyelination (Fig. 9). Furthermore, based on these data, the ET-RJ100 group predominantly diminished the demyelination of the Multiple sclerosis rats (Fig. 9).
The biomarkers in multiple sclerosis are defined as candidates alterable in incidence, progression, and medication reaction [17, 18]. We highlighted genes and microRNAs with the potential role in prognosis and follow-up of multiple sclerosis status based on the artificial intelligence survey.
Analysis of the closest microarray datasets of the brain, spinal cord, and blood samples of multiple sclerosis conditions manifested essential hub genes involved in the demyelination and multiple sclerosis pathogenesis. Moreover, we highlighted molecular signaling pathways associated with these hub genes. Based on enrichment and data mining, we found that inflammation elements, extracellular matrix organization, adipogenesis, TGF-B signaling and activation, JAK/STAT signaling, and TNF/stress oxidative pathway with principal roles are involved in multiple sclerosis development.
The current study revealed that all ten bioactive compounds in Royal Jelly contain Genistein, Kaempferol, Formononetin, Vitamin B2, Vitamin A, and Naringenin Isosakuranetin, 24-methylenecholesterol, Hesperetin, and Vitamin E with suitability binding affinity and RMSD. Moreover, based on the artificial intelligence analysis, we found that these bioactive compounds could bind to DNAJB1 and TNFSF14. Molecular docking computing of bioactive compounds revealed that Royal Jelly’s consumption might change critical molecular signaling pathways to ameliorate life quality and increase muscle endurance in multiple sclerosis conditions. Based on the new result, we reported that binding affinity of molecular library docking, 24-methylenecholesterol with the lowest energy (− 7.3 kcal/mol) bound to TNFSF14 and DNAJB1 proteins surface. These results led to the change of TNF pathways via Royal Jelly consumption. Furthermore, Yujuan Sun et al. demonstrated that the 24-methylenecholesterol bioactive compound could reverse neuron injuries and central neuron system dysfunction in multiple sclerosis. Based on these data, Yujuan Sun et al. suggested that 24-methylenecholesterol had anti-aging activity, anti-oxidative stress function, and neuroprotective effects .
Moreover, the Naringenin compound has improved Alzheimer’s disease, Type 2 Diabetes, and Cancer in several clinical trials studies [20,21,22]. Experiments were carried out by incorporating 5% Naringenin compound into mice’s diets or injecting 20–80 mg/kg of Naringenin compound had immunomodulatory effects and could regulate T-regs and balance the Th1/Th2 ratio in autoimmune arthritis, resulting in less inflammation [23,24,25]. Also, growing evidence has been shown that the Hesperetin compound had neuroprotective properties in experimental models of neurodegenerative disorders . Moreover, the pharmacophore model design illustrated that the Kaempferol compound and vitamin B2 could align over the DDDRR models of the HSPA4 ligands library. Hence, vitamin B2 and Kaempferol compounds could be recommended as effective ligands for HSPA4.
This study revealed that the Royal Jell compound could improve mobility, a proinflammatory cytokine, and demyelination of Multiple sclerosis rats. Notably, among the different concentrations of Royal Jelly (50 mg/kg, 100 mg/kg), we found that 100 mg/mL had an effective and practical impact on proinflammatory cytokine and demyelination of Multiple sclerosis rats. Ample evidence has demonstrated that Royal Jelly suppressed the IL-1β, IL-18, and TNF-α levels. In addition, Royal Jelly has scavenging properties, which could reduce the antiradical and increase the anti-oxidative effect [27,28,29]. Interestingly, Hattori et al. has indicated that consumption of the Royal Jelly stimulated neuron differentiation and inhibited the production of astrocytes oligodendrocytes. In addition, Hattori found that the Royal Jelly compound had Neurotrophic effects [30, 31]. The other studies demonstrated that Royal Jelly attenuated the amyloid-beta deposition and cholesterol concentration as well as the cAMP/PKA/CREB/BDNF pathway significantly increased in mice treated with Royal Jelly. Furthermore, the autonomic nervous systems and blood–brain barrier were improved by consuming Royal Jelly [32, 33]. Moreover, we found that the relative expression of the Dnajb1/Dnajb1/Foxp1/Tnfsf14 and Hspa4 network was regulated by 100 mg/mL Royal compared with 50 mg/mL.
The characteristic and behavior finding showed that gain of weight significantly was improved in multiple sclerosis rats models were treated with 100 mg/kg Royal jelly (RJ-100 group) in comparison M.Sc group. This result indicated that consumption of 100 mg/kg of Royal jelly (RJ-100 group) is significantly more effective in gaining weight than 50 mg/kg (RJ-50 group). Notably, the interaction of the exercise training and consumption of Royal jelly (ET-RJ100 group) predominantly improved the gain weight of the multiple sclerosis rats model compared with other groups. Also, behavioral data displayed that balance and motor coordination in multiple sclerosis rat models were significantly enhanced in 100 mg/kg Royal jelly (RJ-100 group) compared to the Control group. Correspondingly, Multiple sclerosis rats were treated with exercise training, and 100 mg/kg Royal jelly (ET-RJ100 group) significantly increased balance and motor coordination compared with other groups. In addition, the highest reduction in the serum IFNγ was accomplished when a combination of exercise training was conducted with 100 mg/kg Royal jelly (ET-RJ100 group). Moreover, data demonstrated that 100 mg/kg of Royal jelly (ET-RJ100 group) is more efficient on expression levels of Dnajb1, Dnajb2, Foxp1, Tnfsf14, and Hspa4 in comparison with 50 mg/kg of Royal jelly (RJ-50 group). The maximum reduction of the Dnajb1/Dnajb2/Foxp1/Tnfsf14 network was observed when Multiple sclerosis rats were treated with exercise training and 100 mg/kg Royal jelly (ET-RJ100), but the relative expression of the Hspa4 significantly increased.
On the other hand, our data revealed that the expression level of the miR-34a-5p and miR155-3p significantly was decreased in Multiple sclerosis rats treated with exercise training and 100 mg/kg Royal jelly (ET-RJ100 group) compared with the other groups. Hence, we could conclude that the combination of the 100 mg/kg Royal jelly and physical activity had the most effective regulation in the microRNA profile compared to other groups. Based on these data, combining the exercise training with 100 mg/kg Royal jelly (ET-RJ100 group) effectively improved multiple sclerosis rat models.
Physical activity and mobility ameliorate modulating neurotransmitters in the central and peripheral nervous systems . Furthermore, based on the evidence, exercise and physical activity enhanced the muscle strength and motor capacity of MS patients [35,36,37]. Evidence has indicated that physical activity could modulate the brain’s structure and connectome [38, 39]. Moreover, based on the literature review, exercise training might ameliorate fatigue, depression, neurogenesis, neuroprotection, remyelination, and cognition [2, 40, 41]. Lourenco and colleagues demonstrated that aerobic exercise could increase the expression level of the Nrf2 as a vital transcription factor and antioxidant enzymes activity. Based on these data, they found that the oxidative stress status was significantly reduced in the spinal cord of rodents . In addition, physical activity-induced several signaling pathways in the different brain regions and neurotransmission . Based on the proteomic analysis, Lozinski et al. revealed that the alternation of genes is related to the metabolism pathway. They indicated that the expression protein of Pgam2 in glycolysis and Ndufb9, Ndufv3, and mt-Nd3 associated with mitochondria significantly changed by exercise training . In this regard, we found that physical activity improved the motor function, proinflammatory cytokine, and demyelination of the Multiple sclerosis rats. In addition, the relative expression of the Dnajb1/Dnajb1/Foxp1/Tnfsf14 and Hspa4 network was regulated by exercise training. Notably, we showed that the relative expression of the Dnajb1/Dnajb1/Foxp1/Tnfsf14 was down-regulated, and the level of the Hspa4 was up-regulated in ET compared with the M.Sc group. In addition, our data revealed that the exercise training and consumption of the Royal Jelly (ET-RJ50 and ET-RJ100) improved the Dnajb1/Dnajb1/Foxp1/Tnfsf14 and Hspa4 network. However, exercise training and 100 mg/kg Royal Jelly (ET-RJ100) significantly regulated the Dnajb1/Dnajb1/Foxp1/Tnfsf14 and Hspa4 network compared with ET-RJ50 (exercise training along with 50 mg/kg Royal Jelly). Hence, based on the characteristics and behavior (motor function, proinflammatory cytokine, and demyelination) and relative expression of the Dnajb1/Dnajb1/Foxp1/Tnfsf14 and Hspa4 network, we could conclude that exercise training along with 100 mg/kg Royal Jelly (ET-RJ100) had more effect on the Multiple sclerosis condition. Also, we have shown that Royal jelly with potent bioactive compounds might bind to hub genes, influence the improvement of multiple sclerosis hallmark and change the relative expression of hub genes involved in neurodegeneration to slight status. On the other hand, growing evidence revealed that physical activity could be rehabilitated processing speed in cognitive defense mechanisms, but there is no influential role in the cognitive domain [43, 44]. Our results found that synchronizing physical activity and Royal Jelly consumption could increase exercise's effect on ameliorating demyelination in multiple sclerosis hallmarks.
The microRNAs as regulators could affect gene expression in many processes. Thus, dysregulation of microRNAs and essential genes involved in the pathogenesis might develop diseases. Alteration of gene expression due to microRNAs binding to 3´UTR region of genes as a mechanism of microRNAs function is recognized that might degrade mRNA or block mRNA translation process. However, microRNAs could target the 5´UTR of coding sequences noncanonically. Here, microRNA might induce mRNA translation, and hence, upregulation of microRNA might lead to overexpressing the target gene [45,46,47]. In the present study, based on an in- silico survey, we found that miR-34a-5p and miR-155-3p binding to DNAJB2 in 5´UTR and CDS regions, respectively. Moreover, the FOXP1 targets miR-34a-5p and miR-155-3p in the CDS sequence. In addition, miR-34a-5p could target TNFSF14, and HSPA4 in CDS. On the other hand, HSPA4 is targeted by miR-155-3p in 3´UTR. Here, we revealed that the relative expression of the miR-34a-5p and miR-155-3p significantly increased in the blood and brain of Multiple sclerosis rats.
Moreover, in this study, we found the same pattern in the relative expression of these microRNA. In addition, the level expression of the miR-34a-5p and miR-155-3p were down-regulated by consumption of the Royal jelly (RJ-50, RJ-100) and exercise training (ET). Furthermore, consumption of the 100 mg/kg Royal jelly predominantly improved the level expression of the miR-34a-5p and miR-155-3p compared with the 50 mg/kg Royal jelly consumption. Marcin P. Mycko et al. demonstrated the upregulation of miR-155 in T-helper cells during multiple sclerosis in animal models with experimental autoimmune encephalomyelitis. Moreover, they showed that the miR-155 targeted the heat shock proteins family-40 containing Dnajb1 and Dnajb2, which could regulate Th17 cell differentiation and decrease experimental autoimmune encephalomyelitis and demyelination .
On the other hand, the literature review highlighted that miR-155 is a pivotal regulator in the inflammation pathway and modulated an autoimmune response. Furthermore, a recent study had shown that inhibition of miR-155 led to the prevention of pathophysiological mechanisms involved in multiple sclerosis incidence . Moreover, Forough Taheri reported that miR-34 has a potential role in Th17 cell differentiation induction, and up-regulation of miR-34 could lead to progressive multiple sclerosis status, and down-regulation of miR-34 might ameliorate multiple sclerosis condition .
Notably, we indicated that the relative expression of the miR-34a-5p and miR-155-3p were significantly reduced by combining the exercise training and consumption of the Royal jelly. Interestingly, our data found that interaction between exercise training and 100 mg/kg Royal jelly had more effect than the other groups. One of the limitations of this study was comparing different doses of Royal jelly. We selected 50 and 100 mg/kg Royal jelly based on the data mining and literature review and found that these doses were common in previous studies . In addition, after verifying our hypothesis, we suggested that consumption of the Royal jelly along with exercise training could conduct in human trials. Thus, these microRNAs could be affective biomarkers in various states of multiple sclerosis and, in addition, might be potential therapeutic candidates for the treatment of multiple sclerosis.
This study presented a list of bioactive compounds with suitable binding affinity over the DNAJB1 and TNFSF14 as the cutpoint proteins in the multiple sclerosis network based on a molecular docking survey. Nevertheless, we used total Royal Jelly as a crude drug in this study containing all of the bioactive compounds in the chemoinformatics section. Thus, we suggest that evaluating any compounds alone and/or together on multiple sclerosis should be examined in future studies to determine the mechanism and function of compounds in improving multiple sclerosis hallmarks.
The artificial intelligence survey results identified hub genes and microRNAs with a potential role in diagnosing and following multiple sclerosis status. This study aimed to demonstrate that the Royal Jell compound (100 mg/kg/day) could improve mobility, proinflammatory cytokine functions, and demyelination of rats exhibiting Multiple Sclerosis-like behaviors. In addition, this study demonstrated that the combination of exercise training and 100 mg/kg Royal jelly had a more effect on regulating microRNA profiles and hub genes in rats with multiple sclerosis.
Material and methods
We analyzed several microarray data sets obtained from the Gene Expression Omnibus database (46) in the central nervous system and blood samples in multiple sclerosis conditions to find principal hub genes involved in multiple sclerosis. Analysis of GSE52139 in spinal samples was applied by considering P.value < 0.05 and logarithm fold change (logFC) threshold + 1.5 highlighted 258 genes as significantly differential gene expression in multiple sclerosis. The genes with significant expression are shown in the heatmap diagram with P.value < 0.001. The protein–protein interaction in this network was first drawn in the STRING server with medium confidence and based on literature evidence . Then, estimating the visualizing parameters of the network such as degree = 5, betweenness centrality = 0.01, and closeness centrality = 0.10 marked 61 nodes as hub proteins in the multiple sclerosis network based on Cytoscape 3.6.0 algorithm .
On the other hand, we analyzed GSE21942 containing 12 blood samples of multiple sclerosis patients and 15 control subjects to identify differential gene expression in multiple sclerosis. Genes with significant expressions were collected based on P.value < 0.001. Moreover, gene expression patterns were determined based on logFC threshold + 1.5, and we designed a heatmap diagram based on differential gene expression with P.value < 0.001. The interaction network between 210 genes with significant differential gene expression was drawn in the STRING database . Considering visualizing parameters in the protein–protein interactions (PPIs) network, including betweenness centrality = 0.01, degree = 5, and closeness centrality = 0.10, we detected 70 hub nodes involved in multiple sclerosis pathogenesis. Enrichment of hub proteins in both samples of multiple sclerosis patients highlighted molecular signaling pathways involved in demyelination pathomechanisms. The Venn diagram determined the common genes between the spinal-brain dataset, blood dataset, hup genes obtained from spinal-brain samples, and hub genes collected from multiple sclerosis blood samples patients . In the next stage, we recovered the critical role of subscripted genes in signaling pathways by enrichment in several databases such as Panther, Reactome, Biocarta, WikiPathway, and KEGG [54, 55]. Finally, we selected five common hub genes with higher betweenness and the most degree among the multiple sclerosis network proteins. These genes contain HSPA4, TNFSF14, DNAJB1, DNAJB2, and FOXP1. Based on the literature review and strong evidence for diagnosis and follow-up of multiple sclerosis stages, the expression level of IFN-γ as a biomarker in multiple sclerosis was assessed in the serum of patients. Thus, in addition to five genes, we have selected measurements of IFN-γ for follow-up of multiple sclerosis patients.
Prediction of microRNAs profile with significant differential gene expression involved in multiple sclerosis highlighted 120 microRNAs with overexpressing pattern and 153 down-regulation in multiple sclerosis based on microarray dataset GSE17846 analysis (P.value < 0.05 and logFC + 0.2). Moreover, prediction based on target genes was conducted in miRWalk, miRNet, and miRDB servers . We selected miR-34a and miR-155-3p for measurement in our in-silico survey according to prediction outcomes. The gene-microRNAs network was designed in the miRNet algorithm. The miRNet algorithm (12) network indicated that miR-34a-5p and miR-155-3p were common between FOXP1 and DNAJB2 TNFSF14 and HSPA4. An analysis of the network was conducted in the miRPath V.3 database  to determine the roles of selected microRNAs. Furthermore, we drew molecular signaling pathways related to highlighted microRNAs based on the Kyoto Encyclopedia of Genes and Genomes analysis (KEGG) by considering P.value < 0.05, microT threshold 0.4, and fisher’s exact text enrichment method.
Computational Molecular Biology outcomes selected five genes, including DNAJB1, DNAJB2, FOXP1, TNFSF14, and HSPA4, based on the highest visualized parameters among proteins involved in the multiple sclerosis network. The three-dimensional structure of these proteins based on X-Ray diffraction was obtained with browsed in the protein data bank (PDB)  with Follow-up ID 3AGZ (DNAJB1), 4J6G (TNFSF14). After preparation and optimization (removal of solvent, ligands, non-complex compounds, and sub-chains) by dock prep tool in UCSF Chimera software 1.8.1 , these protein structures have been prepared for molecular docking technique. To predict the effects of Royal jelly bioactive compounds on macromolecules, we required computing binding affinity and estimating binding stability between ligands and target protein. Hence, we found ten bioactive compounds in Royal jelly based on the literature review [13, 14, 60, 61]. These bioactive molecules’ three-dimensional (3D) structure was browsed in the PubChem database and obtained in SDF format . We used Open Bable G.U.I. software to construct the compounds library . Estimating of molecular docking between 10 compounds and macromolecules (DNAJB1 and TNFSF14) were conducted in search space with dimensions (x = 79.3563, y = 95.4487, and z = 67.6549) and (x = 47.3635, y = 66.4363, and z = 49.6427) respectively for DNAJB1 and TNFSF14 by PyRx software . This molecular docking study considered the acceptable threshold for binding affinity < -5 and root-mean-square deviation of atomic positions < 2 .
On the other hand, we browsed reported ligands for the heat-shocked proteins family HSP70, in which HSPA4 is member 4 in this family in the Binding DataBase . Pharmacophore modeling was designed based on OPLS 2005, desalting, and generating tautomers in LigPrep, considering 80% of active ligands and 4 to 5 features in PhasePharma and LigPrep applications of maestro 11.5 . Based on the maximum survival score given to the pharmacophore model, we selected the suitable model and screened the alignment of 10 bioactive ligands obtained from Royal jelly on the pharmacophore model. The Royal jelly with complex and variable bioactive compounds relevant to seasonal and geographical area conditions is influenced by growth, development, reproduction, and lifespan. Although the content of compounds is variable, Royal jelly includes proteins, lipids such as free fatty acids, sugar, various vitamins, essential amino acids, and water. Based on the data mining, we found that 10-hydroxy-2-decenoic acid is the principal compound and could play a crucial role in several biological processes such as inflammation and oxidative stress, estrogen-like function, antibacterial, and modulating the immune system. Hence, evidence proposed that Royal jelly with various pharmacological activities is identified as an efficient supplement in therapeutic medication. The current study revealed that all ten bioactive compounds in Royal Jelly contain Genistein, Kaempferol, Formononetin, Vitamin B2, Vitamin A, and Naringenin Isosakuranetin, 24-methylenecholesterol, Hesperetin, and Vitamin E have an efficient supplement in therapeutic medication.
Moreover, in this chemoinformatic study, we described the binding affinity of 10 major bioactive compounds' Royal jelly over the cutpoint hub proteins, DNAJB1 and TNFSF14. We displayed that they all bind to cutpoint proteins with stable and suitable binding affinity. Since there are sufficient concentrations of any compounds in Royal Jelly, we used a cocktail of crude bioactive compounds in the experimental study.
On the other hand, it is worth noting that in this cheminformatics survey, we presented a list of bioactive compounds with suitable binding affinity over the cutpoint proteins. A list of bioactive compounds in Royal Jelly, including 17 compounds, was obtained by browsing the previous studies, that chemoinformatics computing highlighted ten bioactive compounds with the most effective molecular mechanism [68, 69].
All study procedures were approved and conducted by the Research Ethics Committees of Islamic Azad University Isfahan (Khorasgan) Branch (IR.IAU.KHUISF.REC.1400.027).
Animal grouping and treatments
Forty-two male-Sprague Dawley rats were provided from Islamic Azad University Marvdasht Branch and were kept under standard conditions (humidity of 65% ± 5, temperature 24 ± 3 °C, and under a 12 h light–dark cycle). Rats freely accessed regular diet in a whole experimental study. We implemented myelin oligodendrocyte glycoprotein (MOG) to induce MSc in rats. For induction autoimmune encephalomyelitis (EAE) as a multiple sclerosis model, 300 mg/kg myelin oligodendrocyte glycoprotein (MOG, Sigma–Aldrich) dissolved in 100 μl phosphate-buffered salines and then were subcutaneously injected. In this study, rats were randomly divided into seven groups (n = 6):
The rats without treatment are called the Control group (n = 6).
The Rat’s induced autoimmune encephalomyelitis (EAE) as multiple sclerosis model by Intraperitoneal injection (IP) myelin oligodendrocyte glycoprotein (MOG), which is called the M.Sc group (n = 6).
Multiple sclerosis rats were treated with 50 mg/kg Royal jelly called the RJ-50 group (n = 6).
Multiple sclerosis rats were treated with 100 mg/kg Royal jelly called the RJ-100 group (n = 6).
Multiple sclerosis rats were treated with exercise training called the ET group (n = 6).
Multiple sclerosis rats were treated with exercise training and 50 mg/kg Royal jelly, the ET-RJ50 group (n = 6).
Multiple sclerosis rats were treated with exercise training and 100 mg/kg Royal jelly, the ET-RJ100 group (n = 6).
The Royal jelly was dissolved in normal saline and administered by gavage in two doses for five weeks (50 mg/kg for RJ-50, ET-RJ50 groups, and 100 mg/kg for RJ-100, ET-RJ100 groups) . Moreover, During the experiment study, body weight was measured weekly. Finally, rats were euthanized using ketamine and xylazine. 4 ml of blood was extracted from each rat's right ventricle and placed in a sterile tube. To separate serum, blood-filled tubes were centrifuged at 1600 g for 15 min at 4 °C. Tissues and serum were collected, snap frozen, and transferred to − 80 °C for storing.
Exercise training protocol
Rats in exercise training groups conducted physical activity on a rat-running wheel for five weeks (5 days/week). The running speed and duration exercise training protocol was set for 30 min at 11 m/min. This device could record wheel revolutions and time. Following a one-week acclimatization period on the rat-running wheel, the rotations per minute (rpm) and duration of the exercise were gradually increased (30 min at a speed of 11 m/min) [70,71,72]. The entire workout session lasted 30 min. The detailed intervention protocol of exercise training is provided in Table 2.
Protocol of royal jelly
The Royal jelly was dissolved in normal saline and administered by gavage (500 μL) in two different doses (50 mg/kg and 100 mg/kg) for five weeks (50 mg/kg for RJ-50, ET-RJ50 groups, and 100 mg/kg for RJ-100, ET-RJ100 groups). Based on the evidence and literature review, we found that royal jelly's 50 mg/kg and 100 mg/kg doses affected the immune system, redox, oxidative stress, antioxidant balance, and free radicals conditions. Moreover, according to symptoms, Multiple sclerosis is an autoimmune, acute inflammatory, and oxidative stress-mediated demyelinating disease of the central nervous system (CNS). Hence, we selected 50 mg/kg and 100 mg/kg doses to evaluate the effect of Royal jelly on MS status [50, 73,74,75,76]
A Rotarod test was performed to evaluate balance and motor coordination. Rats were conducted three times per day for two days in this study. After rats were trained to acclimate to the rotarod device, the rotation speed gradually increased to 10 rpm for 30 min. The time score was calculated for each group as mean balance and motor coordination on the last day [77, 78].
Vertical Pole Test
According to the previous evidence, we also evaluated the Rat’s movement disorder by Vertical Pole Test. this device included a pole (3-cm diameter) and 150-cm height. Each Rat was placed on the pole and measured the time of rats who could remind on the pole .
According to the manufacturer’s protocol, the serum Interferon-gamma (IFNγ) was measured by commercially available enzyme-linked immunosorbent assay ELISA Kit (Abcam, ab239425) .
RNA isolation and quantitative real-time polymerase chain reaction
To confirm differential hub genes expression obtained from the bioinformatics survey, we sacrificed the multiple sclerosis Rattus models after the treatment stage with Royall jelly consumption and exercise; the brain tissues of animals were separated and snap-frozen in Liquid nitrogen and transferred to the genetics laboratory for storing in − 80 °C. In addition, the blood samples were collected, and the serum was separated for further experiments and tests. Trizol reagent (Thermo Scientific, USA) extracted the total RNA of brain tissues and blood samples based on the manufacturer’s protocol. Quantity and quality of total RNAs were verified in 260/280 nm ratio with NanoDrop spectrophotometry (Thermo Scientific, Waltham, MA, USA) and electrophoresis on the agarose gel. Treatment with DNase1 for extracted RNAs samples was performed to remove DNA contamination based on the manufacturer’s instructions (TaKaRa, Japan). For microRNA expression measurement, cDNA was synthesized by the exiqon manufacturer’s protocol without DNase1 treatment. The synthesized cDNAs were stored at − 20 °C until the next phase. Quantitative real-time PCR was performed using SYBR Green dye (TaKaRa, Kusatsu, Japan) on Rotor-Gene 6000 instrument (Corbett Life Science, Mortlake, Australia).
Sequences of primers were designed in BEACON designer and oligo seven software. The reference gene for normalization of microRNA expression is selected U6snRNA, and normalization of mRNAs expression was performed compared to the reference gene GAPDH 2−ΔΔCt methods analyzed real-time PCR data.
Immediately after the rats sacrificing, the fresh hippocamp region was separated from brain tissue. After that, the hippocamp was fixed in 10% buffered formalin for 24 h at room temperature, and then the hippocamp region was embedded in paraffin in this study. After that, fixed tissues were sliced into 5 m thick slices. Next, tissue was rehydrated and de-paraffinized with xylene, ethanol, and water based on the manuscript protocol. After that, they were stained with hematoxylin and eosin (H&E). Images were analyzed with an Olympus microscope.
We calculated the analysis of variance with GraphPad Prism Software (Version 9 Graph Pad Software Inc., La Jolla, CA). Moreover, the Shapiro–Wilk test was analyzed for normalizing distribution. Furthermore, Data were calculated by t-test and one-way analysis of variance (ANOVA) with Tukey’s post hoc test. Data are indicated as mean ± standard deviation (SD), and differences at P < 0.05 demonstrated the significance in all analyses.
Availability of data and materials
All of the raw data and the rest of the materials are remained in Islamic Azad University Isfahan (Khorasgan) Branch and are available upon request. The datasets analysed during the current study are available in the Gene Expression Omnibus database repository, https://www.ncbi.nlm.nih.gov/geo/, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE52139 PERSISTENT WEB, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE21942, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi, STRING server (https://string-db.org/), Venn diagram (https://bioinfogp.cnb.csic.es/tools/venny/), http://www.pantherdb.org/https://reactome.org/, https://maayanlab.cloud/Harmonizome/dataset/Biocarta+Pathways, https://www.wikipathways.org/index.php/WikiPathways, https://www.kegg.jp/, http://mirwalk.umm.uni-heidelberg.de/, https://www.mirnet.ca/, http://mirdb.org/, http://snf-515788.vm.okeanos.grnet.gr/, https://www.rcsb.org/, https://www.rcsb.org/structure/3AGZ, https://www.rcsb.org/structure/4J6G, https://pubchem.ncbi.nlm.nih.gov/compound/Genistein, https://pubchem.ncbi.nlm.nih.gov/compound/Kaempferol, https://pubchem.ncbi.nlm.nih.gov/compound/Formononetin, https://pubchem.ncbi.nlm.nih.gov/compound/Riboflavin, https://pubchem.ncbi.nlm.nih.gov/compound/Retinol, https://pubchem.ncbi.nlm.nih.gov/compound/Naringenin, https://pubchem.ncbi.nlm.nih.gov/compound/Isosakuranetin, https://pubchem.ncbi.nlm.nih.gov/compound/24Methylenecholesterol, https://pubchem.ncbi.nlm.nih.gov/compound/Hesperetin, https://pubchem.ncbi.nlm.nih.gov/compound/Vitamin-E
Quantitative real-time PCR
One-way analysis of variance
Myelin oligodendrocyte glycoprotein
Heat shock protein family with 40 kDa
Heat shock protein family with 70 kDa
Genome-wide association studies
forkhead box protein P1 loci
Tumor necrosis factor superfamily
logarithm fold change
Kyoto encyclopedia of genes and genomes analysis
Aldewachi H, Al-Zidan RN, Conner MT, Salman MM. High-throughput screening platforms in the discovery of novel drugs for neurodegenerative diseases. Bioengineering. 2021;8(2):30.
Lozinski BM, de Almeida LGN, Silva C, Dong Y, Brown D, Chopra S, Yong VW, Dufour A. Exercise rapidly alters proteomes in mice following spinal cord demyelination. Sci Rep. 2021;11(1):1–16.
Farrokhi M, Moeini P, Fazilati M, Nazem H, Faraji S, Saadatpour Z, Fadaei E, Saadatpour L, Rezaei A, Ansaripour S. Polymorphisms in CD14 gene may modify soluble CD14 Levels and represent risk factors for multiple sclerosis. Immunol Investig. 2016;2:1–8.
Mousavi SR, Tahmasebivand M, Khorrami M, Ayromlou H, Khalili SK, Khorvash F, Rikhtegar R, Khademi B, Bahmanpour Z, Emamalizadeh B. Connection of miR-185 and miR-320a expression levels with response to interferon-beta in multiple sclerosis patients. Mul Scler Relat Disord. 2020. https://doi.org/10.1016/j.msard.2020.102264.
Peschl P, Bradl M, Höftberger R, Berger T, Reindl M. Myelin oligodendrocyte glycoprotein: deciphering a target in inflammatory demyelinating diseases. Front Immunol. 2017;8:529.
Chiricosta L, Gugliandolo A, Bramanti P, Mazzon E. Could the heat shock proteins 70 family members exacerbate the immune response in multiple sclerosis? an in silico study. Genes. 2020;11(6):615.
Tian Z, Song Y, Yao Y, Guo J, Gong Z, Wang Z. Genetic etiology shared by multiple sclerosis and ischemic stroke. Front Genet. 2020;11:646.
Sonar S, Lal G. Role of tumor necrosis factor superfamily in neuroinflammation and autoimmunity. Front Immunol. 2015;6:364.
Mycko MP, Cichalewska M, Cwiklinska H, Selmaj KW. miR-155-3p drives the development of autoimmune demyelination by regulation of heat shock protein 40. J Neurosci. 2015;35(50):16504–15.
Annese T, Tamma R, De Giorgis M, Ribatti D. microRNAs biogenesis, functions and role in tumor angiogenesis. Frontiers Oncol. 2020. https://doi.org/10.3389/fonc.2020.581007.
Singh S, Singh TG, Rehni AK. An insight into molecular mechanisms and novel therapeutic approaches in epileptogenesis. CNS & Neurol Disord-Drug Targets. 2020;19(10):750–79 (Formerly Current Drug Targets-CNS & Neurological Disorders).
Rayatpour A, Farhangi S, Verdaguer E, Olloquequi J, Ureña J, Auladell C, Javan M. The cross talk between underlying mechanisms of multiple sclerosis and epilepsy may provide new insights for more efficient therapies. Pharmaceuticals. 2021;14(10):1031.
Gu L, Zeng H, Maeda K. 10-Hydroxy-2-decenoic acid in royal jelly extract induced both filaggrin and amino acid in a cultured human three-dimensional epidermis model. Cosmetics. 2017;4(4):48.
Kazemi V, Eskafi M, Saeedi M, Manayi A, Hadjiakhoondi A. Physicochemical properties of royal jelly and comparison of commercial with raw specimens. Jundishapur J Nat Pharm Prod. 2019. https://doi.org/10.5812/jjnpp.64920.
Wytrychowski M, Chenavas S, Daniele G, Casabianca H, Batteau M, Guibert S, Brion B. physicochemical characterisation of french royal jelly: comparison with commercial royal jellies and royal jellies produced through artificial bee-feeding. J Food Compos Anal. 2013;29(2):126–33.
Pryor WM, Freeman KG, Larson RD, Edwards GL, White LJ. Chronic exercise confers neuroprotection in experimental autoimmune encephalomyelitis. J Neurosci Res. 2015;93(5):697–706.
Ziemssen T, Akgün K, Brück W. Molecular biomarkers in multiple sclerosis. J Neuroinflammation. 2019;16(1):1–11.
Paap B, Hecker M, Koczan D, Zettl U. Molecular biomarkers in multiple sclerosis. J Clin Cell Immunol S. 2013;10(4):34–65.
Sun Y, Lin Y, Cao X, Xiang L, Qi J. Sterols from mytilidae show anti-aging and neuroprotective effects via anti-oxidative activity. Int J Mol Sci. 2014;15(12):21660–73.
Aroui S, Fetoui H, Kenani A. Natural dietary compound naringin inhibits glioblastoma cancer neoangiogenesis. BMC Pharmacol Toxicol. 2020;21(1):1–10.
Syed AA, Reza MI, Shafiq M, Kumariya S, Singh P, Husain A, Hanif K, Gayen JR. Naringin ameliorates type 2 diabetes mellitus-induced steatohepatitis by inhibiting RAGE/NF-κB mediated mitochondrial apoptosis. Life Sci. 2020;257: 118118.
Wu B-W, Guo J-D, Wu M-S, Liu Y, Lu M, Zhou Y-H, Han H-W. Osteoblast-derived lipocalin-2 regulated by miRNA-96-5p/Foxo1 advances the progression of Alzheimer’s disease. Epigenomics. 2020;12(17):1501–13.
Ahmad SF, Zoheir KM, Abdel-Hamied HE, Ashour AE, Bakheet SA, Attia SM, Abd-Allah AR. Amelioration of autoimmune arthritis by naringin through modulation of T regulatory cells and Th1/Th2 cytokines. Cell Immunol. 2014;287(2):112–20.
Wang J, Qi Y, Niu X, Tang H, Meydani SN, Wu D. Dietary naringenin supplementation attenuates experimental autoimmune encephalomyelitis by modulating autoimmune inflammatory responses in mice. J Nutr Biochem. 2018;54:130–9.
Guo Y-X, Zhang Y, Gao Y-H, Deng S-Y, Wang L-M, Li C-Q, Li X. Role of plant-derived natural compounds in experimental autoimmune encephalomyelitis: a review of the treatment potential and development strategy. Front Pharmacol. 2021;12:1626.
Khan A, Ikram M, Hahm JR, Kim MO. Antioxidant and anti-inflammatory effects of citrus flavonoid hesperetin: special focus on neurological disorders. Antioxidants. 2020;9(7):609.
Ahmad S, Campos MG, Fratini F, Altaye SZ, Li J. New insights into the biological and pharmaceutical properties of royal jelly. Int J Mol Sci. 2020;21(2):382.
Aslan Z, Aksoy L. Anti-inflammatory effects of royal jelly on ethylene glycol induced renal inflammation in rats. Int Braz J Urol. 2015;41:1008–13.
Kohno K, Okamoto I, Sano O, Arai N, Iwaki K, Ikeda M, Kurimoto M. Royal jelly inhibits the production of proinflammatory cytokines by activated macrophages. Biosci Biotechnol Biochem. 2004;68(1):138–45.
Hattori N, Nomoto H, Fukumitsu H, Mishima S, Furukawa S. Royal jelly and its unique fatty acid, 10-hydroxy-trans-2-decenoic acid, promote neurogenesis by neural stem/progenitor cells in vitro. Biomed Res. 2007;28(5):261–6.
Hattori N, Nomoto H, Fukumitsu H, Mishima S, Furukawa S. Royal jelly-induced neurite outgrowth from rat pheochromocytoma PC12 cells requires integrin signal independent of activation of extracellular signalregulated kinases. Biomed Res. 2007;28(3):139–46.
Zhang X, Yu Y, Sun P, Fan Z, Zhang W, Feng C. Royal jelly peptides: potential inhibitors of β-secretase in N2a/APP695swe cells. Sci Rep. 2019;9(1):1–11.
You M, Pan Y, Liu Y, Chen Y, Wu Y, Si J, Wang K, Hu F. Royal jelly alleviates cognitive deficits and β-amyloid accumulation in APP/PS1 mouse model via activation of the cAMP/PKA/CREB/BDNF pathway and inhibition of neuronal apoptosis. Frontiers Aging Neurosci. 2019;10:428.
Abedpoor N, Taghian F, Hajibabaie F. Physical activity ameliorates the function of organs via adipose tissue in metabolic diseases. Acta Histochem. 2022;124(2): 151844.
Sangelaji B, Kordi M, Banihashemi F, Nabavi SM, Khodadadeh S, Dastoorpoor M. A combined exercise model for improving muscle strength, balance, walking distance, and motor agility in multiple sclerosis patients: a randomized clinical trial. Iranian J neurol. 2016;15(3):111.
Halabchi F, Alizadeh Z, Sahraian MA, Abolhasani M. Exercise prescription for patients with multiple sclerosis; potential benefits and practical recommendations. BMC Neurol. 2017;17(1):1–11.
Sharif K, Watad A, Bragazzi NL, Lichtbroun M, Amital H, Shoenfeld Y. Physical activity and autoimmune diseases: get moving and manage the disease. Autoimmun Rev. 2018;17(1):53–72.
Prosperini L, Fanelli F, Petsas N, Sbardella E, Tona F, Raz E, Fortuna D, De Angelis F, Pozzilli C, Pantano P. Multiple sclerosis: changes in microarchitecture of white matter tracts after training with a video game balance board. Radiology. 2014;273(2):529–38.
Kjølhede T, Siemonsen S, Wenzel D, Stellmann J-P, Ringgaard S, Pedersen BG, Stenager E, Petersen T, Vissing K, Heesen C. Can resistance training impact MRI outcomes in relapsing-remitting multiple sclerosis? Mult Scler J. 2018;24(10):1356–65.
Sandroff BM, Baird JF, Silveira SL, Motl RW. Response heterogeneity in fitness, mobility and cognition with exercise-training in MS. Acta Neurol Scand. 2019;139(2):183–91.
Dalgas U, Langeskov-Christensen M, Stenager E, Riemenschneider M, Hvid LG. Exercise as medicine in multiple sclerosis—time for a paradigm shift: preventive, symptomatic, and disease-modifying aspects and perspectives. Curr Neurol Neurosci Rep. 2019;19(11):1–12.
Lourenco MV, Frozza RL, de Freitas GB, Zhang H, Kincheski GC, Ribeiro FC, Gonçalves RA, Clarke JR, Beckman D, Staniszewski A. Exercise-linked FNDC5/irisin rescues synaptic plasticity and memory defects in Alzheimer’s models. Nat Med. 2019;25(1):165–75.
Langeskov-Christensen M, Hvid LG, Jensen HB, Nielsen HH, Petersen T, Stenager E, Hämäläinen P, Dalgas U. Efficacy of high-intensity aerobic exercise on cognitive performance in people with multiple sclerosis: a randomized controlled trial. Mult Scler J. 2021;27(10):1585–96.
Gharakhanlou R, Wesselmann L, Rademacher A, Lampit A, Negaresh R, Kaviani M, Oberste M, Motl RW, Sandroff BM, Bansi J. Exercise training and cognitive performance in persons with multiple sclerosis: a systematic review and multilevel meta-analysis of clinical trials. Mult Scler J. 2021;27(13):1977–93.
Seok H, Ham J, Jang E-S, Chi SW. MicroRNA target recognition: insights from transcriptome-wide non-canonical interactions. Mol Cells. 2016;39(5):375.
Kim DH, Sætrom P, Snøve O, Rossi JJ. MicroRNA-directed transcriptional gene silencing in mammalian cells. Proc Natl Acad Sci. 2008;105(42):16230–5.
Xiao M, Li J, Li W, Wang Y, Wu F, Xi Y, Zhang L, Ding C, Luo H, Li Y. MicroRNAs activate gene transcription epigenetically as an enhancer trigger. RNA Biol. 2017;14(10):1326–34.
Maciak K, Dziedzic A, Miller E, Saluk-Bijak J. miR-155 as an important regulator of multiple sclerosis pathogenesis a review. Int J Mol Sci. 2021;22(9):4332.
Taheri F, Ebrahimi SO, Shareef S, Reiisi S. Regulatory and immunomodulatory role of miR-34a in T cell immunity. Life Sci. 2020;262: 118209.
Malekinejad H, Ahsan S, Delkhosh-Kasmaie F, Cheraghi H, Rezaei-Golmisheh A, Janbaz-Acyabar H. Cardioprotective effect of royal jelly on paclitaxel-induced cardio-toxicity in rats. Iran J Basic Med Sci. 2016;19(2):221.
Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, Doncheva NT, Legeay M, Fang T, Bork P. The STRING database in 2021: customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021;49(D1):D605–12.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Oliveros JC: VENNY. An interactive tool for comparing lists with Venn Diagrams. http://bioinfogp.cnb.csic.es/tools/venny/index.html. 2007.
Xie Z, Bailey A, Kuleshov MV, Clarke DJ, Evangelista JE, Jenkins SL, Lachmann A, Wojciechowicz ML, Kropiwnicki E, Jagodnik KM. Gene set knowledge discovery with enrichr. Current protocols. 2021;1(3): e90.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.
Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48(D1):D127–31.
Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T, Dalamagas T, Hatzigeorgiou AG. DIANA-miRPath v3. 0 deciphering microRNA function with experimental support. Nucleic Acids Res. 2015;43(1):460–6.
Burley SK, Bhikadiya C, Bi C, Bittrich S, Chen L, Crichlow GV, Christie CH, Dalenberg K, Di Costanzo L, Duarte JM. RCSB Protein Data Bank: powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences. Nucleic Acids Res. 2021;49(D1):D437–51.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12.
Oshvandi K, Aghamohammadi M, Kazemi F, Masoumi SZ, Mazdeh M, Vardanjani MM. 2020 Effect of royal jelly capsule on quality of life of patients with multiple sclerosis: A double-blind randomized controlled clinical trial. Iranian Red Crescent Medical Journal, 22(10).
Hegazi A, Al-Menabbawy K. Abd El Rahman E, Helal S: novel therapeutic modality employing apitherapy for controlling of multiple sclerosis. J Clin Cell Immunol. 2015;6(299):2.
Kim S, Chen J, Cheng T, Gindulyte A, He J, He S, Li Q, Shoemaker BA, Thiessen PA, Yu B. PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res. 2021;49(D1):D1388–95.
O'Boyle, N. M., Banck, M., James, C. A., Morley, C., Vandermeersch, T., & Hutchison, G. R. (2011). Open Babel: An open chemical toolbox. Journal of cheminformatics, 3(1), 1-14.
Dallakyan S, Olson AJ. Small-molecule library screening by docking with PyRx. In: Hempel JE, Williams CH, Hong CC, editors. Chemical biology. New York: Springer; 2015.
Mosaddeghi P, Eslami M, Farahmandnejad M, Akhavein M, Ranjbarfarrokhi R, Khorraminejad-Shirazi M, Shahabinezhad F, Taghipour M, Dorvash M, Sakhteman A. A systems pharmacology approach to identify the autophagy-inducing effects of Traditional Persian medicinal plants. Sci Rep. 2021;11(1):1–15.
Gilson MK, Liu T, Baitaluk M, Nicola G, Hwang L, Chong J. BindingDB in 2015: a public database for medicinal chemistry, computational chemistry and systems pharmacology. Nucleic Acids Res. 2016;44(D1):D1045–53.
Bhachoo J, Beuming T. Investigating protein–peptide interactions using the Schrödinger computational suite. In: Schueler-Furman O, London N, editors. Modeling peptide-protein interactions. New York: Springer; 2017.
Collazo N, Carpena M, Nuñez-Estevez B, Otero P, Simal-Gandara J, Prieto MA. Health promoting properties of bee royal jelly: food of the queens. Nutrients. 2021;13(2):543.
Ramadan MF, Al-Ghamdi A. Bioactive compounds and health-promoting properties of royal jelly: a review. J Funct Foods. 2012;4(1):39–52.
Tajiri N, Yasuhara T, Shingo T, Kondo A, Yuan W, Kadota T, Wang F, Baba T, Tayra JT, Morimoto T. Exercise exerts neuroprotective effects on Parkinson’s disease model of rats. Brain Res. 2010;1310:200–7.
Bernardes D, Oliveira ALRd. Regular exercise modifies histopathological outcomes of pharmacological treatment in experimental autoimmune encephalomyelitis. Frontiers Neurol. 2018. https://doi.org/10.3389/fneur.2018.00950.
Bernardes D, Oliveira ALR. Comprehensive catwalk gait analysis in a chronic model of multiple sclerosis subjected to treadmill exercise training. BMC Neurol. 2017;17(1):1–14.
Nejati V, Zahmatkesh E, Babaei M. Protective effects of royal jelly on oxymetholone-induced liver injury in mice. Iran Biomed J. 2016;20(4):229.
Almeer RS, AlBasher GI, Alarifi S, Alkahtani S, Ali D, Abdel Moneim AE. Royal jelly attenuates cadmium-induced nephrotoxicity in male mice. Sci Rep. 2019;9(1):1–12.
Kunugi H, Mohammed Ali A. Royal jelly and its components promote healthy aging and longevity: from animal models to humans. Int J Mol Sci. 2019;20(19):4662.
Abdelnour SA, Abd El-Hack ME, Alagawany M, Taha AE, Elnesr SS, Abd Elmonem OM, Swelum AA. Useful impacts of royal jelly on reproductive sides, fertility rate and sperm traits of animals. J Anim Physiol Anim Nutr. 2020;104(6):1798–808.
Abedpoor N, Taghian F, Ghaedi K, Niktab I, Safaeinejad Z, Rabiee F, Tanhaei S, Nasr-Esfahani MH. PPARγ/Pgc-1α-Fndc5 pathway up-regulation in gastrocnemius and heart muscle of exercised, branched chain amino acid diet fed mice. Nutr Metab. 2018;15(1):1–15.
Kim J-Y, Yi E-S, Lee H, Kim J-S, Jee Y-S, Kim S-E, Kim C-J, Ko I-G. Swimming exercise ameliorates symptoms of MOG-induced experimental autoimmune encephalomyelitis by inhibiting inflammation and demyelination in rats. Int Neurourol J. 2020;24(Suppl 1):S39.
Zhao Q, Yin L, Zhang L, Jiang D, Liu L, Ji H. Chitoheptaose promotes heart rehabilitation in a rat myocarditis model by improving antioxidant, anti-inflammatory, and antiapoptotic properties. Oxid Med Cell Longev. 2020. https://doi.org/10.1155/2020/2394704.
We thank our colleagues for their association and helpful discussions in this study.
There is no Funding.
Ethics approval and consent to participate
All study procedures were approved and conducted by the Research Ethics Committees of Islamic Azad University Isfahan (Khorasgan) Branch (IR.IAU.KHUISF.REC.1400.027). Moreover, All study procedures were confirmed that all methods were carried out in accordance with relevant guidelines and regulations. Also, all methods are reported in accordance with ARRIVE guidelines (https://arriveguidelines.org) for the reporting of animal experiments.
Consent for publication
There is no competing of interest to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lohrasbi, M., Taghian, F., Jalali Dehkordi, K. et al. The functional mechanisms of synchronizing royal jelly consumption and physical activity on rat with multiple sclerosis-like behaviors hallmarks based on bioinformatics analysis, and experimental survey. BMC Neurosci 23, 34 (2022). https://doi.org/10.1186/s12868-022-00720-0
- Multiple sclerosis
- Hub genes
- Royal Jelly
- Physical activity