Recurrent aphthous stomatitis is a complex inflammatory disease caused by multiple factors, which severely influences the life of patients. The etiology of RAS lesions is unknown, but some local trauma, systemic, immune, genetic, allergic, nutritional, and microbiological factors are considered to be causative factors(Chiang et al., 2019). However, there are no effective drugs with low side effects against RAS.
HB has several biological activities against many Gram-positive and some Gram-negative bacteria(García-de-la-Mària et al., 2018; Higaki, Hasegawa, Morohashi, & Takayoshi, 1995). Besides its antibacterial activity, HB also acts as a well-established anti-inflammatory agent. Recent studies have shown that HB has both in vitro and in vivo anti-inflammatory effects(Mao et al., 2010; Park, Chung, Kim, Kwon, & Joh, 2007). Previous studies have shown the effectiveness of HB against RAS(Chen CH, 2009; L, 2021; Yang JM, 2013). Nevertheless, the mechanism is unclear. Presently, a comprehensive network pharmacology method was established to study the therapeutic targets and possible mechanisms of HB for treating RAS.
The chemical composition of HB mainly includes alkaloids, flavonoids, sterols, lactones, phenolic acids, terpenes, phenylpropanoids, and other chemical components using liquid chromatography-tandem mass spectrometry(Sun et al., 2016; H. Wang et al., 2013; X. Xian et al., 2014). Alkaloids are the main active ingredients in HB. About 42 alkaloids have been reported in the current studies. In the present study, quercetin and rutaecarpine were among the important active compounds of HB. Quercetin is a common flavonoid in Chinese herbs with therapeutic properties, including antioxidant and anti-inflammatory activities(Harwood et al., 2007). Quercetin has a long history of consumption as part of the normal human diet(Harwood et al., 2007). In a prospective active control trial, researchers confirmed the effectiveness of quercetin in reducing ulcer size and pain in patients with RAS(Pandya, Kalappanavar, Annigeri, & Rao, 2017). Quercetin was shown to be a safe and acceptable highly effective novel adjuvant therapy in two randomized controlled trials(Arafa, Ghalwash, El-Kersh, & Elmazar, 2018; Hamdy & Ibrahem, 2010). In addition, 90 percent of patients responded that they appreciated the ease of application when using the topical quercetin, and they did not object to its consistency or taste(Hamdy & Ibrahem, 2010). It may be promising to formulate a quercetin-based drug treatment for RAS in the future. Rutaecarpine is a major indolequinazoline alkaloid isolated from Evodiae fructus the dried fruit of Evodia rutaecarpa (Juss.) Benth(Nguyen et al., 2013). An increasing number of studies have revealed that rutaecarpine has extensive pharmacological actions, ranging from anti-inflammatory, anti-oxidant effects to vasodilatory effects(Cho, Shim, Lee, Mar, & Kim, 2005; Hibino, Yuzurihara, Kase, & Takeda, 2009; Ko et al., 2007). Mechanistic studies have shown evodiamine inhibits hypoxia induced inflammatory responses (COX-2 and iNOS expression) via suppressing the activation of hypoxia inducible factor-alpha (HIF-α) in macrophages(Y. N. Liu et al., 2009). Notably, many studies have shown that rutaecarpine can reduce pro-inflammatory mediators such as TNF-α, IL-1β, IL-6, and NF-κB to improve the symptoms of inflammation and oxidative stress(Han, Hu, & Chen, 2019; Jin et al., 2017; Y. Li et al., 2019). Taken together, we believe that these compounds play an important role in the healing process of RAS.
In the present study, integrated bioinformatics methods assisted in an analysis of how critical genes change in their expression to uncover potential RAS pathways based on GEO datasets “GSE37265”, and we identified 578 DEGs, including 535 up-regulated DEGs and 43 down-regulated DEG. Further, using the network pharmacology approach, we identified 27 intersection genes with HB treatment against RAS. The DGE-analysis identified 23 up-regulated DEGs and 4 down-regulated DEGs. PPI network analysis demonstrated that one active ingredient matches multiple targets and, in turn, one target matches multiple ingredients; reflecting the complex mechanism of action of multi-ingredient and multi-target against diseases in Chinese medicine. Gene pathway network analysis revealed that TNF, IL6, IL1B, PTGS2, CCL2, and IL10 were among the core targets of the HB in the treatment of RAS. Previous studies suggested that TNF-α plays a vital role in the pathological process of RAS, Oliveira(Oliveira et al., 2016) et al. found in their experiment that reducing TNF-α expression optimized the healing of traumatic oral ulcers in diabetic rats. TNF-α, as a pro-inflammatory factor, is involved in the development of RAS and can be used as an indicator of the severity of the disease. PTGS2 encoding cyclooxygenase-2(COX-2) has been documented to play a role in the development of RAS, Experiments in rats showed a decrease in PTGS2 activity and a decrease in COX-2 protein expression during the healing of oral ulcers(Slomiany & Slomiany, 2002).
GO enrichment analysis suggested that the HB regulates a variety of BPs and affects various CCs and MFs. GO functional enrichment analysis were mainly on biological processes, including response to lipopolysaccharide, response to molecule of bacterial origin, cellular response to lipopolysaccharide, and cellular response to molecule of bacterial origin. Cellular response to molecule of bacterial origin and reactive oxygen species metabolic process are important BPs involved in the development of RAS. wound healing and aging are key contributors that cause RAS. Reactive oxygen species metabolic process may trigger the stress response by regulating gene expression to shift redox homeostasis to a more oxidative state. This process facilitates the onset of RAS(Apel & Hirt, 2004; Liochev, 2013). Positive regulation of protein phosphorylation is a central signaling pathway that regulates cell survival, playing significant roles in the inhibition of RAS(Balasuriya et al., 2020; W. Yang, Zhang, & Wang, 2017).
The mechanisms of HB for treating RAS are fulfilled at least two function modules, including the inhibition of inflammation and the regulation of immunological response. The implementation of the functions relies on a smooth run of the complex multi-pathways network, particularly TNF signaling pathway and HIF-1 signaling pathway. CCL2, CXCL10, IL6, MMP3, TNF, IRF1, SELE, VCAM1, IL1B, PTGS2, ICAM1and FOS are involved in TNF-1 signaling pathway. Previous studies have been reported that excessive production of TNF-α, IL-6, and IL-1β was associated with an increased risk of RAS development(Bazrafshani, Hajeer, Ollier, & Thornhill, 2002; Guimarães et al., 2007; Scully & Porter, 2008). It has also been reported that matrix metalloproteinases play a negative role in mucosal pathology of the oral cavity(Karasneh, Bani-Hani, Alkhateeb, Hassan, & Thornhill, 2014; Skulason, Holbrook, & Kristmundsdottir, 2009). In summary, the anti-inflammatory effect is an important factor for HB contributing to the treatment of oral ulcers. Three targets are involved in HIF-1 signaling pathway which is an important regulator of the immune system, including IL6, NOS2, and IFNG. NOS2 is the isoforms of NOS, which produce high-level sustained NO synthesis and strongly affect adaptive immune responses(Bogdan, 2015). NO regulates the responses of many immune and inflammatory cell types like macrophages which are associated with some of the most important immune pathologies(Niedbala et al., 2014). Increased serum NO level is thought to take part in the pathogenesis of RAS(Gurel, Altinyazar, Unalacak, Armutcu, & Koca, 2007). Thus, NOS2 plays an important role in the HIF-1 signaling pathway for the regulation of anti-inflammation and immune protection functions of HB. Immunologic disturbances play a crucial role in the etiopathogenesis of oral ulcers. Hence, it is worthwhile to explore the impact of HB on immune function. According to the pathway enriched results, many target genes are enriched in the pathways related to immune inflammatory diseases. For example, nine targets are enriched in Toll-like receptor signaling pathway (P-value = 1.45 × 10 − 10). Six targets are enriched in NF-kappa B signaling pathway (P-value = 4.94× 10 − 6) and four targets in T cell receptor signaling pathway (P-value = 1.16 × 10 − 6). Therefore, it can be presumed that one pharmacological function of HB is related to the regulation of immune response for RAS treatment.
Molecular docking displayed components have good binding to core target genes. All selected core proteins had good affinity for the ligand (<-5 kcal/mol) and all effective docking fractions were > -7 kcal/mol, indicating that the compound has a strong binding affinity for the docked protein(Li, Fu, & Zhang, 2019). In silico validation of key targets for bioactive components offers an alternative path for the exploration of ligand-target interactions and action mechanisms. We could conclude that the different components display diverse interactions with these targets in the predicted pathways. The results demonstrated a synergistic effect mode of the HB in its oral ulcer treatment effects. PTGS2, which has the highest docking score, is a biomarker of iron death and inhibits the expression of inflammatory factors and apoptosis(H. Zhou et al., 2019; Y. Zhou et al., 2021). Future research directions can start from iron death mechanism pathway.
The limitation of this study is the lack of clinical or animal experiments. Further studies will combine molecular biology and pathophysiology to validate the predicted potential key targets and pathways, and explore the mechanism of action of the active ingredients.