Anti-Osteogenic Effect of Danshensu in Ankylosing Spondylitis: An in Vitro Study Based on Integrated Network Pharmacology

Ankylosing spondylitis (AS) is a chronic inflammatory disease characterized by abnormal bone metabolism, with few effective treatments available. Danshensu [3-(3,4-dihydroxy-phenyl) lactic acid) is a bioactive compound from traditional Chinese medicine with a variety of pharmacologic effects. In the present study, we investigated the pharmacologic effect and molecular mechanism of Danshensu in AS. Potential targets of Danshensu were identified in four drugs-genes databases; and potential pharmacologic target genes in AS were identified in three diseases-genes databases. Differentially expressed genes related to AS were obtained from the Gene Expression Omnibus database. Overlapping targets of Danshensu and AS were determined and a disease–active ingredient–target interaction network was constructed with Cytoscape software. Enrichment analyses of the common targets were performed using Bioconductor. To test the validity of the constructed network, an in vitro model was established by treating osteoblasts from newborn rats with low concentrations of tumor necrosis factor (TNF)-α. Then, the in vitro model and AS fibroblasts were treated with Danshensu (1–10 μM). Osteogenesis was evaluated by alkaline phosphatase staining and activity assay, alizarin red staining, quantitative PCR, and western blotting. We identified 2944 AS-related genes and 406 Danshensu targets, including 47 that were common to both datasets. The main signaling pathways associated with the targets were the c-Jun N-terminal kinase (JNK) and extracellular signal-regulated kinase (ERK) pathways. A low concentration of TNF-α (0.01 ng/ml) promoted the differentiation of osteoblasts; this was inhibited by Danshensu, which had the same effect on AS fibroblasts but had the opposite effect on normal osteoblasts. Danshensu also decreased the phosphorylation of JNK and ERK in AS fibroblasts. There results provide evidence that Danshensu exerts an anti-osteogenic effect via suppression of JNK and ERK signaling, highlighting its therapeutic potential for the treatment of AS.


INTRODUCTION
Ankylosing spondylitis (AS) is a chronic progressive inflammatory disease with insidious onset characterized by chronic inflammatory back pain, bone resorption, and new bone formation, with the latter leading to ankylosis and functional disability (Poddubnyy and Sieper, 2017;Lories, 2018). The pathogenesis of AS and therapeutic strategies targeting ectopic new bone formation have been the focus of many studies (Qin et al., 2018;van Tok et al., 2019). Heterotopic ossification involves endochondral and membranous bone formation and cartilage metaplasia, in which fibroblasts and bone cells such as mesenchymal stem cells, osteoblasts, and osteoclasts play an essential role (Van Mechelen et al., 2018;Liu et al., 2021). Inflammatory cytokines such as tumor necrosis factor (TNF)-α, interleukin (IL)-17A, and IL-6 are also implicated in new bone formation van Tok et al., 2019;Zou et al., 2020). TNF-α is known to inhibit osteoblast differentiation but at low concentrations, it not only enhanced proliferation but also promoted osteoblast differentiation in AS via Wnt signaling and subsequent bone formation via nuclear factor (NF)-κB (p65) and c-Jun N-terminal kinase (JNK)/activator protein (AP)-1 signaling (Frost et al., 1997;Li et al., 2018). Besides, low concentration of TNF-α enhanced the osteogenic differentiation of AS fibroblasts, which was accompanied by the expression of osteogenesis markers including Runt-related transcription factor (RUNX)2, osteopontin (OPN), and osteocalcin (OCN) (Zou et al., 2020).
As advanced ankylosis of the spine is the worst long-term outcome of AS, preventing ectopic new bone formation is a major focus of research. However, there are few drugs available that prevent ankylosis progression (Poddubnyy and Sieper, 2020). Biologic agents targeting TNF-α and IL-17A are important treatments that alleviate AS symptoms by suppressing inflammation and slowing spinal radiographic progression. However, not all patients respond to these treatments, which can increase the risk of infection; moreover, the high cost of these agents imposes a considerable economic burden on patients (Yin et al., 2020). The therapeutic effect of TNF-α blockers in new bone formation of AS still leave a question open. It has been shown to suppress inflammation but also enhance radiographic progression in some studies (Maksymowych et al., 2009;Pedersen et al., 2011). Therefore, there is an urgent need for alternative, effective medicines for the treatment of AS.
Danshensu [3-(3,4-dihydroxy-phenyl) lactic acid], is a bioactive compound from traditional Chinese medicine--Salvia miltiorrhiza (Danshen) with a variety of pharmacologic effects (Zeng et al., 2017;Wang et al., 2021;Xie et al., 2021). Danshensu has been integrated into the treatment of cardiovascular diseases (e.g., myocardial ischemia and reperfusion, atherosclerosis, and hypertension) and cerebrovascular diseases (e.g., cerebral embolism) owing to its capacity to modulate inflammation and exert antioxidative and anti-apoptotic effects . Danshensu was shown to suppress neuroinflammation by reducing the levels of proinflammatory cytokines such as TNF-α and IL-6 in the hippocampus , and exerted a protective effect against myocardial ischemia/reperfusion injury by blocking the phosphorylation and activation of JNK and nuclear translocation of nuclear factor (NF)-κB, which is involved in new bone formation in AS (Meng et al., 2016;Li et al., 2018). Danshensu plays an important role in bone physiology by regulating osteoblast differentiation and function. An H 2 S-releasing compound derived from Danshensu was found to have an antioxidant effect in an H 2 O 2 -induced cell injury model in MC3T3-E1 osteoblasts by attenuating the activation of the p38 mitogen-activated protein kinase (MAPK), extracellular signal-regulated kinase (ERK)1/2, and JNK-MAPK pathways (Yan et al., 2017). On the contrary, in our previous study we found that salvianolic acid B, another bioactive component of Danshen, promoted osteogenesis of normal human mesenchymal stem cells (hMSCs) by enhancing ERK1/2 signaling (Xu et al., 2014). However, little is known about the application or effect of Danshensu in AS.
Based on the above findings, we speculated that Danshensu can prevent ossification of osteoblasts and fibroblasts in AS. To test this hypothesis, in this study we investigated the effect and mechanism of action of Danshensu on osteoblast differentiation in AS. We first conducted a network pharmacology analysis to identify potential targets of Danshensu and analyzed the function and pathways related to the targets. We then constructed a disease-drug component-target-signaling pathway interaction network to identify important pathways, and evaluated the efficacy and mechanism of action of Danshensu using osteoblasts and AS fibroblasts (Figure 1).

Danshensu Target Identification
The Simplified Molecular Input Line Entry System (SMILES) and International Chemical Identifier key (InChiKey) for Danshensu were retrieved from the PubChem database. Danshensu was searched in BATMAN-TCP (prediction score of the target gene ≥20; adjusted p value of the prediction result <0.05), TCMSP, SwissTargetPrediction, and Pharmamapper databases based on the results of the search.

Identification of Pharmacologic Targets in AS
The keyword "ankylosing spondylitis" was used as input in the OMIM, GeneCards, and DisGeNET databases to identify genes related to AS. The original microarray data of GSE25101 and GSE73754 were analyzed with GEO2R to identify differentially expressed genes (DEGs) and generate a volcano plot; according to log 2 fold change (FC), the −log 10 (adjusted p values) were divided into upregulated genome, downregulated genome, no statistical difference, and nonexpressed difference genome (screening conditions: adjusted p value <0.01, |log 2 FC|≥1.0). Both AS-related gene and DEG names were normalized and merged and duplicate items were deleted, and the obtained gene sets were used as disease targets.

Standardization and Screening for Potential AS Targets of Danshensu
All targets were searched in the UniProt database with the species set to human, and the collected protein or gene information was corrected to the standard UniProt name (Symbol) of the target. The set of Danshensu targets was mapped to the set of disease targets and intersected to obtain the potential targets of Danshensu in AS.

Identification of Hub Genes Shared Between Danshensu Target and AS Genes
Topologic properties, degree, and interaction of the distribution network from the STRING database were analyzed using the cytoHubba plugin of Cytoscape.

Gene Ontology Functional Annotation and KEGG Pathway Enrichment Analysis
The drug component-disease target gene names were encoded by UniProt IDs. GO functional annotation of the targets and KEGG pathway enrichment analysis were performed using Rx64 v4.1.0 software with bioc-Lite.R, clusterProfliter.R, and pathview.R packages from the Bioconductor platform.

Construction of Disease-Drug Component-Target-Signaling Pathway Interaction Network
Cytoscape v3.8.2 software was used to construct the disease-drug component-target-signaling pathway interaction network, in which the nodes represented the disease, drug component, target, and signaling pathways and edges represented interactions between drug component and target, target and AS, or target and signaling pathway. The topologic properties of the network were analyzed using the Network Analysis plugin of Cytoscape.

Cell Culture
Hip capsular ligaments were collected from three AS patients at Guangdong Provincial People's Hospital. The study was approved by the ethics committee of the hospital. Fibroblasts were extracted from the hip capsular ligaments as previously described (Qin et al., 2018). Briefly, ligaments were washed twice with phosphate-buffered saline (PBS) and then cut into small pieces that were placed in 10-mm culture dishes. After 3-4 h, Dulbecco's Modified Eagle's Medium (DMEM; Gibco, Shanghai, China) containing 15% fetal bovine serum (FBS; Gibco) and 1% penicillin-streptomycin (Solarbio, Beijing, China) was added to the culture dishes. After cells had adhered, the ligament pieces were removed and the medium was replaced with fresh medium containing 10% FBS and 1% penicillin-streptomycin. The medium was changed every 3-4 days and cells from passage 3-5 were used for experiments.
Primary osteoblasts from the skull of 24-h newborn Sprague-Dawley rats were provided by Guangdong Medical University. Briefly, the calvaria was cut into small pieces and digested at 37°C with trypsin for 15 min, followed by 0.1% collagenase I for 20 min. The supernatant was centrifuged at 800 rpm for 10 min, and the cell precipitate was collected and inoculated in 25-cm 2 culture vials while the supernatant was transferred back to the centrifuge tube containing the bone fragments for additional digestion. The above steps were repeated 3 times. The cells were grown in DMEM containing 10% FBS and 1% penicillin-streptomycin at 37°C in a humidified atmosphere of 5% CO 2 . The culture medium was refreshed every 2-3 days. When the cells reach 80% confluence, they were digested with 0.25% EDTA-trypsin (Solarbio) and passaged in a 1:2 or 1:3 ratio. Cells from passage 3-5 were used for experiments.

Cell Proliferation Assay
The 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay was used to evaluate the effect of Danshensu on osteoblast and AS fibroblast proliferation. The cells were seeded in 96-well plates at a density of 3×10 3 (osteoblasts) and 5×10 3 /well (fibroblasts), respectively. After 24 h, the cells were divided into nine groups (5 wells per group) that were treated with Danshensu at a concentration of 0, 0.1, 0.3, 1, 3, 10, 30, 100, or 300 μM. After incubation for 7 days, 20 μl MTT (5 mg/ml) was added to each well. The mixed medium was discarded 4 h later and replaced with 150 μl dimethylsulfoxide (Solarbio), and the plate was agitated for 20 min at room temperature. A microplate reader (Bio-Rad Laboratories, Hercules, CA, United States) was used to measure the optical density (OD) of each well at a wavelength of 570 nm.

Alkaline Phosphatase Staining and Activity Assay
Cells were washed twice with PBS and fixed with 4% paraformaldehyde for 30 min. After three washes with PBS,

Alizarin Red Staining
Cells were washed twice with PBS and then fixed with 70% ethanol for 20 min at room temperature; they were then washed with water and stained with 0.2% alizarin red (Sigma, Tokyo, Japan) solution (pH 4.1) for 30 min. The solution was discarded and the cells were washed 3 times with distilled water.

RNA Isolation and Quantitative (q)PCR
After culturing for 14 days, total RNA was extracted from cells grown in 6-wells plates using TRIzol reagent (Ambion, Austin, TX, United States). The RNA was reverse transcribed using the PrimeScript RT reagent kit (Takara, Dalian, China) and qPCR was carried out to detect mRNA levels of Runt-related transcription factor (RUNX)2, osterix (SP7), osteocalcin (OCN), bone morphogenetic protein (BMP)2, and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) (internal control). Cycling conditions were as follows: 95°C for 5 min, followed by 40 cycles of 95°C for 15 s and 60°C for 60 s. Forward and reverse primers were synthesized by Sangon Biotech (Shanghai, China) and are shown in Table 1.

Western Blotting
After culturing for 14 days, cells were lysed with RIPA lysis buffer supplemented with phenylmethylsulfonyl fluoride (Beyotime), and 30 μg of protein per well were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis on a 10% acrylamide gel and electrotransferred to a polyvinylidene difluoride membrane that was incubated with antibodies against GAPDH (1:

Statistical Analysis
All data are presented as mean ± SD. SPSS v22.0 software was used to analyze the data. One-way analysis of variance was used to assess differences between groups, with a least significant difference post hoc test when groups had equal variance (p > 0.05) and Dunnett's post hoc test in the alternative case. Differences with a p value < 0.05 were considered statistically significant.

Identification of Danshensu Targets
We identified 21 putative targets of Danshensu ( Figure 2A) in BATMAN-TCM along with 22 in TCMSP, 100 in SwissTargetPrediction, and 273 in Pharmmapper ( Figure 2B). The identified genes were aggregated and duplicates and inconsistently named genes were removed, leaving 406 genes as potential targets.

Identification of AS-Related Genes
Using the keyword "ankylosing spondylitis," 108 AS-related genes were identified in OMIM along with 2,274 in GeneCards and 710 in DisGeNET. A total of 324 and 61 DEGs were obtained from GSE25101 and GSE73754, respectively ( Figure 2D). The identified genes were aggregated and duplicates and inconsistently named genes were removed, leaving 2,944 AS-related genes ( Figure 2C).

Identification and Verification of Hub Genes Among Danshensu Targets and AS Genes
Cytoscape was used to visualize key genes as a network from which overlapping genes were selected for further analysis by determining the node degree or interaction of sequenced candidate genes using the cytoHubba application. The top 10 hub genes in the network of Danshensu AS targets were IL6, ESR1, PTGS2, EGFR, IL10, TGFB1, TIMP1, MMP1, NR3C1, and RELA (Table 2 and Figure 3C).

Drug Component-Target GO Function Annotation and KEGG Pathway Enrichment Analysis
The target gene names (symbol) were transformed into gene IDs by searching the UniProt database, and the biocLite.R package of Bioconductor was used to perform GO functional annotation and KEGG pathway enrichment analysis of the target genes using Rx 64 v4.0.1 software, with p 0.05 and Q 0.05 as the threshold values. The top 40 GO terms are shown in Figure 4. Most were related to inflammatory response, oxidative stress, the MAPK cascade, etc. The KEGG pathway enrichment analysis yielded 13 signaling pathways among which ERK and JNK signaling pathways were highly enriched, suggesting that they are key pathways for the anti-osteogenic effect of Danshensu (Figure 4).

Construction of a Disease-Drug Component-Target-Signaling Pathway Interaction Network
The drug component-disease target, active ingredient molecular ID, and drug component-disease-signaling pathway interaction data were imported into Cytoscape, and a drug component-target interaction network was generated using the Merge function in the software (Figures 3D,E). These targets were presumed to be key targets of Danshensu in AS.

Danshensu has Little Effect on Osteoblast and AS Fibroblast Proliferation
To investigate the effect of Danshensu on osteoblast and AS fibroblast proliferation, the cells were exposed to increasing concentrations of Danshensu ranging from 0.1 to 300 μM. After culturing for 7 days, the MTT assay was performed. Danshensu did not affect the proliferation of osteoblasts and AS fibroblasts at concentrations between 0.1 and 30 μM, while at 100 μM osteoblast proliferation was significantly inhibited and at 300 μM, both osteoblast and AS fibroblasts were affected ( Figure 5). We therefore used 1, 3, and 10 μM Danshensu for the following experiment.

Low Concentrations of TNF-α Promote Osteoblast Differentiation
We investigated the effect of TNF-α on osteoblast differentiation by ALP staining and activity assay after 7 days of treatment. Both ALP staining and activity were increased in osteoblasts by treatment with 0.01 ng/ml TNF-α for 7 days ( Figure 6A). We also examined the effect of TNF-α on mineralization in osteoblasts by ARS. There were more calcium nodules in the group treated with 0.01 ng/ml TNF-α than in the other groups ( Figure 6B). TNF-α at concentrations 1 and 10 ng/ml was not toxic to osteoblasts, and there were no calcium nodules in the 10 ng/ml group.
To confirm the effect of low TNF-α concentrations on osteoblast differentiation, we evaluated the expression of BMP2 and SP7-2 genes related to this process-by qPCR. Application of 0.01 ng/ml TNF-α resulted in the upregulation of BMP2 and SP7 ( Figure 6C). We also examined the levels of RUNX2 and COL1A1 proteins by western blotting and found that they were significantly increased in the presence of 0.01 ng/ml TNF-α, whereas a concentration of 10 ng/ml reduced the protein level of RUNX2 ( Figure 6D). Based on these results, 0.01 ng/ml TNF-α was used in subsequent experiments to create a proinflammatory environment for inducing osteoblast differentiation.

Danshensu Inhibits TNF-α-Induced Osteogenic Differentiation
To investigate the anti-osteogenic effect of Danshensu, osteoblasts were treated with 0.01 ng/ml TNF-α and different concentrations of Danshensu, and ALP staining and activity were evaluated 7 days later. Danshensu inhibited ALP activity in TNFα-induced osteoblasts in a dose-dependent manner. ALP staining was also reduced in cells treated with Danshensu compared to those treated with 0.01 ng/ml TNF-α only ( Figure 7A). Similar results were obtained by ARS ( Figure 7B). Additionally, qPCR analysis showed that RUNX2, BMP2, and SP7 mRNA expression was reduced in cells treated with Danshensu compared to those treated with TNF-α only ( Figure 7C) and after 14 days, cells exposed to 0.01 ng/ml TNF-α and 1, 3, or 10 μM Danshensu had lower levels of COL1A1 and RUNX2 proteins compared to those treated with TNF-α only ( Figure 7D).

Danshensu Promotes Ossification of Normal Osteoblasts
To determine whether Danshensu affects normal osteoblasts, we assessed the osteogenic potential of osteoblasts treated with 1, 3, and 10 μM of Danshensu. The results of ALP staining and the activity assay showed that Danshensu increased osteogenic differentiation ( Figure 8A); this was confirmed by ARS ( Figure 8B). The qPCR analysis showed that RUNX2, BMP2, and SP7 mRNA expression was increased by Danshensu treatment, an effect that was concentration-dependent ( Figure 8C); and a similar effect was observed on COL1A1 and RUNX2 protein levels ( Figure 8D).

Danshensu Suppresses Ossification of AS Fibroblasts
As fibroblasts are known to contribute to heterotopic ossification, we extracted primary fibroblasts from hip ligaments of AS patients and treated these cells with Danshensu (1, 3, and 10 μM) to further assess whether Danshensu inhibits ectopic ossification in AS. The ALP staining and activity assay showed that Danshensu suppressed osteogenic differentiation of AS fibroblasts ( Figure 9A). This was accompanied by decreased expression of COL1 and RUNX2 proteins ( Figure 9B) and downregulation of OCN, SP7, and RUNX2 transcripts ( Figure 9C) as determined by western blotting and qPCR, respectively. ARS confirmed the results of ALP staining ( Figure 9D).

Danshensu Suppresses JNK and ERK Phosphorylation in AS Fibroblasts
Based on the results of the network pharmacology analysis, we examined the phosphorylation levels of JNK and ERK in AS fibroblasts treated with Danshensu for 14 days. The results TABLE 2 | Top 10 genes in the protein-protein interaction network ranked by the maximal clique centrality method.

Rank
Name Score showed that p-JNK and p-ERK levels were significantly downregulated in cells treated with 3 and 10 μM Danshensu ( Figure 10).

DISCUSSION
AS is a chronic inflammatory disease with complex pathophysiology and no effective treatments (Magrey and Khan, 2017). TNF-α and IL-17A play important roles in the development of AS . Although TNF-α blocker and IL-17A inhibitor are used for AS treatment, the side effects and high cost limit their widespread use. Traditional Chinese medicine is known for its multitarget, bidirectional, and systemic regulatory effects with fewer side effects but has not been applied to AS. Danshensu is a pure molecule with a clearly defined structure that is extracted from the widely used traditional Chinese medicine Danshen . This study provides the first evidence for the anti-osteogenic effect of Danshensu in AS based on in silico and in vitro data.
Network pharmacology allows visualization of the relationship between biological processes and drug action through construction of an interaction network . The mechanism of action of drugs can be explained by analysis of network components including drug targets, diseases, and molecular pathways (Li et al., 2020). It is particularly suitable for the study of multicomponent, multitarget, and multi-mechanism traditional Chinese medicines. In this study, we used a network pharmacologybased approach to identify targets, biological processes, and signaling pathways of Danshensu in AS. Danshensu target genes shared by AS were mostly related to the inflammatory response, cell proliferation and differentiation, fibrosis, and arthritis. The top gene, IL6, is known to promote AS fibroblast ossification via MAPK/ERK signaling . The results of the GO functional annotation and KEGG pathway enrichment analysis showed that the therapeutic effect of Danshensu may be related to the regulation of inflammatory response; protein transport; oxidative stress; cell apoptosis; and ERK, JNK, and MAPK signaling. The disease-drug component-target-signaling pathway interaction network revealed that the JNK and ERK pathways were closely related to the 10 hub genes, suggesting that they play a critical role in the antiosteogenic effect of Danshensu. The ERK and MAPK pathways were also related to most AS targets of Danshensu. The MAPK pathway including ERK, JNK and p38 pathways, is involved in inflammatory signaling, regulation of cell proliferation and differentiation, and activation of the immune system (Lai et al., 2020). The ERK1/2, JNK, and p38 pathways regulate osteoblast differentiation (Chen et al., 2016), and inhibition of MAPK signaling was shown to reduce the expression levels of ALP and OCN (Wang et al., 2019). JNK is a downstream effector of the proinflammatory cytokine-induced immune response in AS. Additionally, Wnt gene promoters contain p65 and c-Jun binding sites, and inhibiting the NF-κB and JNK/AP-1 pathways suppressed the TNF-α-induced upregulation of Wnt expression  and attenuated IL-17-induced inflammatory injury (Mao et al., 2019) in AS.
We confirmed the anti-osteogenic effect of Danshensu in AS using an in vitro model established by treating osteoblasts with low concentrations of TNF-α and AS fibroblasts. Fibroblasts and osteoblasts originate from mesenchymal stem cells and share a similar phenotype and differentiation pathway that makes their reciprocal transformation possible (Zou et al., 2016). AS fibroblasts have higher osteogenic potential than healthy control fibroblasts and may undergo ectopic ossification (Qin et al., 2018;Zhang et al., 2021). Consistent with the previous finding that mild inflammation promotes new bone formation of AS , we confirmed in the present study that treatment with a low concentration dose of TNF-α promoted osteoblast differentiation in normal osteoblasts. We therefore used 0.01 ng/ml TNF-α to mimic the AS inflammatory microenvironment in vitro. Our results showed that Danshensu prevented the ossification of TNF-α-treated osteoblasts in a dose-dependent manner. Notably, osteogenic activity markers were upregulated in osteoblasts exposed to Danshensu alone. These results indicate that Danshensu inhibits osteogenesis in osteoblasts exposed to low-intensity inflammatory stimulation while preserving their osteogenic potential. This effect of Danshensu is likely related to its regulation of inflammatory response as revealed in network pharmacology analysis. Importantly, Danshensu also suppressed osteogenic differentiation in AS fibroblasts, demonstrating that Danshensu can prevent abnormal new bone formation in AS.
The results of the GO and KEGG pathway enrichment analyses along with the target-signaling pathway interaction network indicated that JNK and ERK signaling pathways are involved in the anti-osteogenic effect of Danshensu in AS. We found that Danshensu treatment decreased the phosphorylation of JNK and ERK in AS fibroblasts, suggesting that Danshensu prevents ectopic ossification in AS via inhibition of these pathways. Interestingly, we previously found that salvianolic acid B, another bioactive component of Danshen, promoted osteogenesis of normal hMSCs by increasing ERK1/2 signaling (Xu et al., 2014). The opposite effects of Danshensu on normal vs abnormal cells may be related to the distinct microenvironments of the two cell types. The immune microenvironment in AS is characterized by activated inflammatory mediators that promote new bone formation (Schett et al., 2017). The proinflammatory cytokine TNF-α at a low concentration may exert this effect, which is countered by Danshensu via modulation of JNK and ERK signaling pathways. However, the detailed regulatory mechanism of Danshensu in the immune microenvironment in AS remains to be elucidated.
Besides JNK and ERK, other pathways that were found to be enriched in our analyses may be involved in the molecular mechanism of Danshensu in AS such as phosphatidylinositol 3-kinase (PI3K)-protein kinase B (AKT) and Janus kinase (JAK)-signal transducer and activator of transcription (STAT). PI3K-AKT signaling promotes the proliferation and differentiation of osteoblast precursors while p38 MAPK signaling regulates the proliferation and differentiation of Frontiers in Pharmacology | www.frontiersin.org November 2021 | Volume 12 | Article 772190 11 osteoclast progenitor cells (Han and Choi, 2017;Bai et al., 2018). The JAK-STAT pathway affects bone formation and the mechanical strength of bone (Sanpaolo et al., 2020). Forkhead box (Fox)O was shown to be associated with the regulation of PI3K-AKT signaling in mechanically stimulated osteoblasts (Li et al., 2010). Activation of hypoxia-inducible factor (HIF) along with an increased supply of oxygen and nutrients was proposed to underlie bone anabolism in mice lacking the oxygen sensor prolyl hydroxylase (PHD)2 (Stegen et al., 2018). Additionally, the enrichment of the apoptosis cascade and calcium signaling Frontiers in Pharmacology | www.frontiersin.org November 2021 | Volume 12 | Article 772190 13 pathway in our KEGG analysis suggests that Danshensu also exerts an anti-osteogenic effect in AS by inducing apoptosis in osteoblasts and promoting the uptake of external calcium ions by the body.
In conclusion, we identified potential targets of Danshensu in AS and showed that Danshensu suppresses the osteogenic differentiation of osteoblasts and AS fibroblasts via activation of JNK and ERK signaling pathways. Although the putative targets of Danshensu in AS need to be validated by in vivo studies, our results provide molecular-level evidence for the therapeutic potential of Danshensu in the treatment of AS.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the ethics committee of Guangdong Provincial