MyD88 dimerization inhibitor ST2825 targets the aggressiveness of synovial fibroblasts in rheumatoid arthritis patients
Arthritis Research & Therapy volume 25, Article number: 180 (2023)
Dimerization of the myeloid differentiation primary response 88 protein (MyD88) plays a pivotal role in the exacerbated response to innate immunity-dependent signaling in rheumatoid arthritis (RA). ST2825 is a highly specific inhibitor of MyD88 dimerization, previously shown to inhibit the pro-inflammatory gene expression in peripheral blood mononuclear cells from RA patients (RA PBMC). In this study, we elucidated the effect of disrupting MyD88 dimerization by ST2825 on the pathological properties of synovial fibroblasts from RA patients (RA SFs).
RA SFs were treated with varying concentrations of ST2825 in the presence or absence of bacterial lipopolysaccharides (LPS) to activate innate immunity-dependent TLR signaling. The DNA content of the RA SFs was quantified by imaging cytometry to investigate the effect of ST2825 on different phases of the cell cycle and apoptosis. RNA-seq was used to assess the global response of the RA SF toward ST2825. The invasiveness of RA SFs in Matrigel matrices was measured in organoid cultures. SFs from osteoarthritis (OA SFs) patients and healthy dermal fibroblasts were used as controls.
ST2825 reduced the proliferation of SFs by arresting the cells in the G0/G1 phase of the cell cycle. In support of this finding, transcriptomic analysis by RNA-seq showed that ST2825 may have induced cell cycle arrest by primarily inhibiting the expression of critical cell cycle regulators Cyclin E2 and members of the E2F family transcription factors. Concurrently, ST2825 also downregulated the genes encoding for pain, inflammation, and joint catabolism mediators while upregulating the genes required for the translocation of nuclear proteins into the mitochondria and members of the mitochondrial respiratory complex 1. Finally, we demonstrated that ST2825 inhibited the invasiveness of RA SFs, by showing decreased migration of LPS-treated RA SFs in spheroid cultures.
The pathological properties of the RA SFs, in terms of their aberrant proliferation, increased invasiveness, upregulation of pain and inflammation mediators, and disruption of mitochondrial homeostasis, were attenuated by ST2825 treatment. Taken together with the previously reported anti-inflammatory effects of ST2825 in RA PBMC, this study strongly suggests that targeting MyD88 dimerization could mitigate both systemic and synovial pathologies in a variety of inflammatory arthritic diseases.
Rheumatoid arthritis (RA) is a chronic autoimmune form of joint disease mediated by systemic and local inflammatory mechanisms which ultimately promote pathogenic and destructive processes in the synovium of RA patients . These processes are mediated by the interplay between immune cells and synovial fibroblasts (SFs) in the arthritic synovium [1, 2]. SFs turn into imprinted aggressors that promote and maintain inflammatory and degenerative processes in RA [3, 4]. The aggressiveness of RA SFs includes increased proliferation, aberrant cell cycle progression, reduced apoptosis, and exacerbated production of pro-inflammatory and joint catabolism mediators [5,6,7,8]. Together, these cellular and molecular changes in the RA SFs promote joint damage in RA and other chronic autoimmune joint diseases [9,10,11]. The increased invasive capacity of the RA SFs is another major characteristic of their aggressiveness, which promotes the formation of the hyperplasic and granulated synovial tissue called pannus and drives cartilage and bone degeneration [12, 13]. These invasive properties of the RA SFs are regulated by the action of several inflammatory mediators. Collectively, these aggressive pathogenic features that resemble the transformed properties of cancer-associated fibroblasts contribute to the destruction of the joint architecture, altered function, and pain sensitization in RA patients.
Innate immune mechanisms that act locally at the level of joints, and systemically at the level of immune cells, sustain the chronic inflammatory milieu in the RA joints . Toll-like receptors (TLRs) belong to a well-characterized group of innate immune receptors that sense distinct damage-associated molecular patterns (DAMPs), pathogen-associated molecular patterns (PAMPs), and pro-inflammatory mediators in the joint and systemic microenvironments . Various autoantibody immune complexes have been defined as the activators of TLR2 and TLR4 pathways [16, 17], strongly supporting their role in systemic and local inflammation in RA. As a result, several efforts have been made to target upstream and downstream molecules involved in the TLR pathways for RA therapy [14, 18, 19]. Preclinical animal models have evidenced the therapeutic potential of targeting TLR2 and TLR4 [20,21,22]. However, despite the attempts to successfully achieve amelioration of RA remission, clinical trials targeting TLR4 have only shown moderate effects , while others remain in the early phase . The lack of significant results in clinical trials suggests that successful targeting of innate immune pathways may require extensive multipronged inhibitory approaches. Therefore, the inhibition of common adaptor proteins that participate in TLR signaling, rather than the receptors themselves, may be a more suitable approach to observe significant clinical effects.
In this regard, myeloid differentiation primary response 88 (MyD88) is a critical adaptor protein for TLR signaling. Intracellular MyD88 dimerization mediates the recruitment of interleukin receptor-associated kinases 1 and 4 (IRAK1/4) and subsequently led to the activation of intracellular signaling cascades, perpetuating inflammatory signaling in the RA synovium . MyD88 activation downstream of TLR2 and TLR4 signaling has been proposed as a key mechanism promoting RA SF proliferation, invasion, and increased production of inflammatory mediators [26,27,28]. Importantly, activating mutations in MyD88 was shown to cause severe inflammatory arthritis in humans . Therefore, targeting MyD88 could potentially represent a promising approach to control SF aggressiveness in RA patients.
Our previous study revealed that the MyD88 dimerization inhibitor, ST2825, can suppress systemic inflammation in RA and suggested that it may also likely modulate gene expression in the RA synovium . The ST2825 structure has been widely described, and its mechanism-of-action for specifically targeting MyD88 is mediated by interfering with homo-oligomerization of the BB loop at the MyD88 TIR domain, thereby affecting its dimerization [31, 32]. Thus, in the present study, we performed an unbiased transcriptomic evaluation of the effect of ST2825 in SFs from RA patients, supported by biological function studies. Our results demonstrated that the blockade of MyD88 dimerization by ST2825 causes cell cycle arrest, reduces invasiveness, downregulates the expression of pain and inflammatory mediators, and upregulates the expression of genes that support mitochondrial function in SFs from RA patients. Taken together, these findings strongly suggest that ST2825 could potentially mitigate both local and systemic inflammation in RA patients.
We aimed to elucidate the effect of disrupting MyD88 dimerization by ST2825 on the pathological properties of SFs from RA patients in an unbiased manner using transcriptomic analysis and biological assays. Our working hypothesis is that the synthetic small molecule ST2825 mitigates the aggressive behavior of RA SFs mediated by their proliferative and antiapoptotic capacity, altered expression of inflammatory mediators, and invasive properties.
LPS from Escherichia coli (CAT-L-2880, SIGMA®) and ST2825 inhibitor of MyD88 dimerization (Cat. No. A3840, APExBIO) were used for hDF and SFs stimulation. All reagents were reconstituted according to the manufacturer’s instructions. DAPI (CS1-0127-2 mL, Nexcelom Bioscience) was used for cell cycle assays.
Experiments on human SFs were performed according to the guidelines and approval of the Emory University Institutional Review Board. OA SFs were prepared from the knee synovium collected from OA patients undergoing total joint replacement as described in our previous publications [33, 34]. RA SFs isolated from the knee synovium of RA patients were collected from cadaveric donors within 48 h of death and were purchased from Articular Engineering, LLC. Dermal fibroblasts from healthy human donors (hDF) were purchased from Lonza. Cells were cultured in DMEM with 10% fetal bovine serum (Corning) and 1% penicillin/streptomycin and utilized between passages 3 and 6. Age and sex of the patients from which the SFs were derived are as described in Fig. S2. Cells were incubated at 37 °C in a humidified 5% CO2 atmosphere and 95% humidity.
Cell cycle and apoptosis analysis
The hDF, OA SFs, and RA SFs were cultured in a 12-well plate after the density adjustment at 1 × 105 cells/mL. The percentages of cells in the G0/G1, S, G2/M phases and apoptotic cells were determined at 24, 48, and 72 h after the incubation with 0, 5, and 10 μM of ST2825. Briefly, cells were permeabilized with 70% ethanol, and their nuclei were stained with DAPI (4 nM) and examined by measuring the amount of DNA per single cell in a multifluorescent channel analysis (blue 377/477) on the Celigo Imaging Cytometer, according to Nexcelom Bioscience protocols (Assay ID: Celigo_02_0001). The percentage of apoptotic cells was determined by gating for intact cells containing less than 2 × amount of DNA calculated at the interface of DAPI integrated intensity on the Celigo Imaging Cytometer (Additional file 1: Fig. S5).
Cell invasion assays were carried out on RA SFs 3D spheroids according to a previously published protocol using the Celigo Imaging Cytometer . RA SFs were plated in Nexcelom3D Ultra-low Attachment Round Bottom 96-well Plates at a density of 4000 cells/well in 200 μL of medium and cultured for a period of 4 days, following which the formation of spheroids was visually confirmed using an Olympus CKX53 Microscope (× 10 magnification). On day 4, the Corning® Matrigel® matrix was thawed on ice according to the manufacturer’s protocols. LPS alone at a concentration of 30 ng/mL or LPS containing 10 μM ST2825 were added to the Matrigel and distributed into the 96-well plates containing RA SF spheroids. For unstimulated RA SFs, pure Matrigel was added. Matrigel containing 10 μM ST2825 was added as a control. Automated image acquisition and quantification of the invaded area in μm2 were performed on the Celigo Imaging Cytometer at 0, 48, 72, and 96 h after the addition of Matrigel, according to the Nexcelom Bioscience protocols (Assay ID: Celigo_03_0002) .
Bulk RNA sequencing
RA SFs were cultured in 6-well plates after the density of 2.5 × 105 cells/well followed by stimulation with 30 ng/mL of LPS, with or without 10 μM of ST2825 for 24 h. ST2825 was added 45 min prior to LPS stimulation. Unstimulated RA SFs were taken as the control group. Additionally, unstimulated OA SFs were cultured and collected after 24 h and were considered as control cells for unstimulated RA SFs. The total RNA was extracted and purified using Direct-zol RNA MicroPrep (Zymo Research) following the manufacturer’s protocol. RNA quality and quantity were assessed using a 2100 Bioanalyzer (Agilent Technologies). Only samples with an RNA integrity number (RIN) > 8.6 were used for sequencing. Three biological replicates were used per condition. Libraries were generated from 250 ng RNA using TruSeq Stranded Total RNA Sample Prep Kit (Illumina). Sequencing was carried out using the NovaSeq 6000 system (Novogene UC Davis Sequencing Center, Novogene Corporation Inc.). FASTQ files from these samples were uploaded to the DNASTAR Lasergene (version 22.214.171.124) and ArrayStar software for analysis. Paired-end reads were mapped to the GRCh37 human genome assembly. RNA levels were normalized using Log2 RPKM (reads per kilobase of exon model per million mapped sequence) by performing Student’s t-test for genes at 95% confidence.
Differential expression and pathway analysis
Differentially expressed genes (DEGs) were identified using the BioJupies web tool . DEGs with p-value < 0.05 and fold change (FC) ≥ 1.5 were used to perform volcano plots, heatmaps, and subsequent pathway analysis. QIAGEN Ingenuity Pathway Analysis (IPA) was used to identify distinct up- and downregulated canonical pathways between different conditions in our study. Data was visualized using GraphPad Prism version 9.4.1, VolcaNoseR, and BioJupies web tools [36, 37]. Protein interaction networks of relevant DEGs were identified and analyzed in STRING version 11.5 , the interaction score was set at high confidence = 0.700, and a false discovery rate (FDR) cutoff of 0.05 was used to determine up- or downregulated pathways in the network.
Gene expression analysis by real-time PCR
The cDNA was synthesized from 1 μg of RNA using qScript™ cDNA SuperMix 5X (QuantaBio). Quantitative real-time polymerase chain reaction (qRT-PCR) was performed using PerfeCTa SYBR® Green FastMix 2X (QuantaBio) on the QuantStudio™ 6 Flex Real-Time PCR System (Applied Biosystems™). For this analysis, at least two biological replicates and 3 technical replicates were used. Fold change (FC) was determined through the 2−ΔΔCt method, and TBP mRNA level was used as the reference gene. Primer sequences are detailed in Additional file 1: Table S1.
Cells were harvested and lysed using Pierce™ IP lysis buffer (Thermo Fisher Scientific) containing 1% Halt™ protease and phosphatase inhibitors (Thermo Fisher Scientific). Equal amounts of protein were electrophoresed on 10% SDS–polyacrylamide gels. The primary antibodies used for the western blot are listed in Additional file 1: Table S2. Bands were visualized using ChemiDoc MP Imaging System (Bio-Rad). For luciferase reporter assays plasmids encoding NF-κB (pNL3.2.NF-kB-RE (Promega) and IL-15 , RA SFs were transiently transfected with Viromer Red reagent (Origene). The pSV2bgal plasmid was used as a control for transfection efficiency. Cells were treated with LPS (35 ng/mL), ST2825 (10 μM), or both for 24 h. Luciferase and β-galactosidase activities were assayed using the GloMax Discover System (Promega). Reporter activities were normalized for transfection efficiency and reported as fold change in luciferase activity. Assays were performed as triplicates.
Normality tests were performed to determine the data distribution and appropriate descriptive and inferential statistics. The results were analyzed using the software GraphPad Prism version 9.4.1. The specific statistical analysis is described in the corresponding figure caption. One-way ANOVA, two-way ANOVA and Dunnett’s multiple comparisons, Student’s t-test, and one-way ANOVA and Tukey’s multiple comparisons test were used. A p-value less than 0.05 was considered statistically significant.
Effect of inhibiting MyD88 dimerization on hDF and OA SFs
To identify the effects of inhibiting MyD88 dimerization, the viability of hDF incubated with 0, 5, 10, and 30 μM of ST2825 for 24 h was determined (Fig. 1A). The results indicated that 5 and 10 μM of ST2825 did not compromise cell viability; therefore, those two concentrations were used for subsequent experiments. To evaluate the long-term effect of ST2825 on proliferation, hDF and OA SFs were monitored for 24, 48, and 72 h after incubation with 5 and 10 μM of ST2825 (Fig. 1B, D). The analysis reported no statistically significant differences among conditions, suggesting that cell viability was not compromised. Although not significant, we observed a drop in cell number, particularly in hDF treated with 10 μM of ST2825. Thus, the DNA content of hDF and OA SFs was quantified to investigate the effect of ST2825 on different phases of the cell cycle (Fig. 1C, E and Additional file 1: Fig. S1). Our analysis revealed that 10 μM of ST2825 induced a significant decrease in the percentage of phase S cells in hDF and OA SFs, with a corresponding increase in the percentage of G0/G1 phase OA SFs at 48 h. Interestingly, these changes were accentuated in OA SFs treated with both 5 and 10 μM of ST2825 at 72 h, suggesting that ST2825 may influence cell proliferation through induction of G0/G1 cell cycle arrest.
Transcriptomic analysis to confirm the inflammatory status of the RA SFs used in this study
Since our goal was to determine whether the MyD88 dimerization inhibition mediated by ST2825 modulates pathogenic processes in RA SFs, we confirmed the inflammatory status of the RA SFs samples used in this study by comparing them with OA SFs isolated from OA synovium of patients undergoing total knee replacement (Additional file 1: Fig. S2A). We identified a total of 1879 differentially expressed genes (DEGs) by 1.5-fold change (p < 0.05) between RA and OA SFs. The analysis revealed 826 downregulated and 1053 upregulated genes (Additional file 1: Fig. S2B). Hierarchical cluster analysis of normalized gene expression from RA and OA SFs shows 2129 genes at 95% confidence that clustered together in the appropriate group (Additional file 1: Fig. S2C). Upregulated and downregulated genes were further analyzed to identify canonical pathways underlying pathogenic features in RA SFs (Additional file 1: Fig. S2D). The overexpression of important inflammatory genes, along with the increased expression of MYD88 observed in RA SFs, established the suitability of the RA SFs samples for this study (Additional file 1: Fig. S2E). Importantly, it supports our hypothesis for targeting this protein in SF as a potential pharmacological approach in the arthritic joints.
ST2825 arrests cell cycle progression in RA SFs
To determine the biological effects of ST2825 on RA SFs, we monitored cell proliferation for 24, 48, and 72 h after incubation with 5 and 10 μM of ST2825 (Fig. 2A). Although cell numbers remained unchanged, the treatment with 10 μM of ST2825 caused a significant increase in the percentage of G0/G1 phase and significantly decreased the percentage of S and G2/M phase RA SFs (Fig. 2B). The same analysis was performed on RA SFs incubated with ST2825 for 48 and 72 h, even though a similar trend was observed, the changes in cell cycle phases were not significant suggesting that the ST2825 acted within 24 h of the treatment (Additional file 1: Figs. S3 and S4). We, therefore, selected the 24 h time point for transcriptomic analysis by RNA-seq. Genes were selected based on the reactome pathway enrichment (Additional file 1: Fig. S6 and Additional file 2: Supplementary Data 1), and the main downregulated pathways were related to the cell cycle, which matches with the results obtained for the analysis of the cell cycle states. Those downregulated genes with 1.5 FC and p < 0.05 were further analyzed through QIAGEN Ingenuity Pathway Analysis (IPA), and common genes across downregulated pathways were selected for heatmaps. The transcriptomic analysis revealed cell cycle regulation among the top five IPA canonical pathways downregulated by ST2825 (Fig. 2C). We next performed a paired differential gene expression analysis in the 3 independent RA SFs samples; heatmaps show that ST2825 treatment successfully downregulated cell cycle-related genes (Fig. 2D). The predicted protein interaction network showed the main gene cluster involved in the regulation of G1 cell cycle arrest (green), G0 and early G1 (red), and Cyclin D-associated events in G1 (blue) at high interaction confidence = 0.700 and PPI enrichment p-value = 1.43E − 07 (Fig. 2E). The percentage of apoptotic cells was evaluated in all time points among conditions; significant differences in apoptosis were only observed at 72 h with 10 μM of ST2825 (Additional file 1: Fig. S5). Taken together, these observations suggest that MyD88 dimerization inhibition arrests cell cycle progression by inducing G0/G1 cell cycle arrest, which may potentially modulate the exacerbated proliferation and apoptosis resistance observed in RA SFs.
ST2825 inhibits the pathological gene expression signature of LPS-stimulated RA SFs
We next investigated the effect of ST2825 under inflammatory conditions by stimulating RA SFs with LPS, a well-known pro-inflammatory inductor known to activate TLR signaling upstream of MyD88 dimerization. RNA sequencing was performed from LPS-stimulated RA SFs after 24 h (Additional file 1: Fig. S7 and Additional file 3: Supplementary Data 2). In order to investigate the possible role of ST2825 in modulating the expression of critical genes involved with RA SF aggressiveness, we manually selected some genes associated with catabolic and inflammatory pathways along with pain mediators (Fig. 3A). The analysis revealed upregulation of critical genes involved in catabolic processes such as CXCL10, TNFSF13B, ADAMTS4 (ADAM metallopeptidase with thrombospondin type 1 motif 4), and IL1B (Interleukin 1 beta). Similarly, the upregulation of several genes related to inflammatory pathways such as TLR1 (Toll-like receptor 1), TLR3 (Toll-like receptor 3), CASP1, and MYD88 was observed. It is important to mention that some of these genes were also upregulated in RA SF vs. OA SF in Additional file 1: Fig. S2, supporting our hypothesis that ST2825 can target distinct aspects of RA SF aggressiveness. One of the most interesting findings was the LPS-mediated upregulation of interleukin 33 (IL33), C–C motif chemokine ligand 2 (CCL2), and prostaglandin-endoperoxide synthase 2 (PTGS2), genes that have been described as pain mediators (Fig. 3A). Differential gene expression analysis also showed that the most downregulated genes by the effect ST2825 on LPS-stimulated RA SFs were those related to cell cycle regulation and DNA mismatch repair (Fig. 3B, E). We successfully validated the ST2825-induced downregulation of the pro-inflammatory genes IL1B and MYD88 and the cell cycle regulators CCNE2 and MYBL2 (Myb-related protein B) genes by qRT-PCR (Fig. 3C, D). The predicted STRING protein interaction network showed that the main gene cluster involved in the regulation of G1 cell cycle arrest (yellow), G0 and early G1 (green), and DNA repair (blue) at high interaction confidence = 0.700 and PPI enrichment p-value = 0.00431 (Fig. 3F). Upregulated genes by the effect ST2825 on LPS-stimulated RA SFs were related to Sirtuin signaling, mitochondrial metabolism, pluripotency, apoptosis, and chondroitin sulfate biosynthesis (Fig. 3G). The predicted STRING protein interaction network also showed that upregulated genes participate in distinct mitochondrial processes; the analysis was performed at high interaction confidence = 0.700, and the PPI enrichment p-value was 4.25E − 08 (Fig. 3H).
ST2825 suppresses the invasive features of LPS-stimulated RA SFs in 3D spheroid cultures
Invasiveness represents a key pathogenic feature of RA SFs due to their contribution to pannus formation and joint destruction. In this study, LPS was used to potentiate the invasive properties of RA SFs. In order to evaluate long-term changes, spheroids were performed and monitored for 48, 72, and 96 h after incubation with LPS or LPS plus ST2825 (Fig. 4A). The area of the invasion reached statistical significance at 96 h, an effect that was significantly suppressed by the MyD88 inhibitor ST2825 (Fig. 4B, C). Taken together, the data obtained from this analysis demonstrates that local treatment with ST2825 may suppress the invasiveness of RA SFs contributing to ameliorating joint damage.
ST2825 targets NF-κB-dependent mechanisms and IKK-related kinases in RA SFs
In order to elucidate the effect of ST2825 on MyD88 and NF-κB signaling pathways in RA SFs, we treated cells with LPS and ST2825 and performed western blot and luciferase reporter assays. Our results show that the levels of MyD88, IRAK4, and TRAF6 were not significantly changed by ST2825 (Fig. 5A and Additional file 1: Fig. S8). However, this result was expected as ST2825 inhibits MyD88 dimerization and subsequent signaling events by disrupting the formation of the Myddosome complex but does not directly affect the gene expression or protein level of MyD88 or its interacting partners. Nonetheless, we were able to demonstrate that the level of p65, which is the downstream effector of MyD88, was downregulated by ST2825 (Fig. 5A). In addition, our results showed that targeting MyD88 by ST2825 significantly inhibits p65-dependent NF-κB binding site reporter luciferase activity (Fig. 5B) and may have an additional effect by inactivating the IKK-related kinases TAK1 and TBK1, and NF-κB p65 in RA SFs treated with LPS (Fig. 5C and Additional file 1: Fig. S8).
ST2825 is a halogenated heptapeptide that competitively binds to the Toll interleukin-1 receptor (TIR) domain of MyD88 and blocks its dimerization, a critical event in the TLR/IL-1 signaling pathway, which culminates in the activation of NF-κB signaling . We previously showed that ST2825 inhibited the overexpression and release of chemokines and cytokines on PBMC [30, 40]. In this study, we show that MyD88 inhibition by ST2825 effectively reduced the aggressiveness of RA SFs via arresting cell cycle progression, inhibiting the expression of pain and inflammatory mediators, decreasing their invasiveness, and promoting mitochondrial function. Based on these findings from our previous and current studies performed in RA PBMC and RA SFs, respectively, we hypothesize that ST2825 could potentially mitigate the chronic activation of both systemic and local inflammatory processes in RA patients (Fig. 6). Our data also show that ST2825 inhibited cell cycle progression in OA SFs and hDF, suggesting that it could be developed as a therapeutic for OA and psoriatic arthritis.
One of the major findings of this study is that ST2825 decreased cell proliferation of RA SFs by arresting the cells in the G0/G1 phase of the cell cycle as early as 24 h after treatment. This arrest in the cell cycle may have contributed to the increased number of apoptotic nuclei at 72 h after treatment . Moreover, pro-apoptotic genes BCL2 like 1 (BCL2L1, also known as BCLXL) and BCL2-binding component 3 (BBC3, also known as PUMA) were upregulated by ST2825 in LPS-stimulated RA SFs, suggesting that ST2825 may relieve the apoptotic resistance of RA SFs. Several studies have reported that MyD88 is linked to cell proliferation. MYD88 activating mutations were shown to result in the aberrant proliferation of multiple types of cancer cells [42,43,44,45]. In particular, MyD88 was reported to play a role in the occurrence and development of breast cancer , epithelial ovarian carcinoma , human hepatocellular carcinoma , osteosarcoma , adenocarcinoma of the colon , and hematopoietic cancer-like lymphoma and leukemia . ST2825 was reported to inhibit the growth and apoptosis resistance of lymphoid malignancies, but the underlying molecular mechanism is incompletely understood . Our findings show that ST2825 may play an inhibitory role in cell cycle progression by primarily downregulating the expression of critical cell cycle checkpoint regulators such as Cyclin E2 (CCNE2) and the members of the E2F family transcription factors (E2F1, E2F2, E2F7, and E2F8). CCNE2 is an important component in cell proliferation  and cancer , which may explain the anti-proliferative activity of ST2825. Regarding the E2F transcription factor family, previous reports have revealed its crucial regulatory role in the timely expression of genes needed for cell cycle progression and proliferation . RA synovial tissue exhibits much greater levels of E2F2 expression. E2F2 has also been linked to RA pathology by potentiating the aggressive characteristics of RA SFs, including proliferation, invasion, and in vitro cytokine production . Moreover, E2F2 binds to the MYD88 promoter and upregulates its expression, which subsequently activates the PI3K/AKT/NF-κB pathway .
The current knowledge suggests that abrogating MyD88 functions will most likely mitigate cell proliferation. According to the increased percentage of cells in phase G0/G1 and decreased percentage in S and G2/M phases found in our study, along with downregulated genes associated with cell cycle progression; the treatment with 10 μM of ST2825 causes a cell cycle arrest in G1. Therefore, cell proliferation can be prevented. These results match with several reports which have demonstrated that TLR4/MyD88 blockade can suppress the cell proliferation of gastric cancer cell lines , hepatocarcinoma HEPG2 cells , human hepatoma cell lines , and head kidney lymphocytes in an infectious model . In RA, abrogating the TLR4/MyD88/NF-κB axis has demonstrated beneficial effects in SFs [58,59,60]. Due to the tumor-like phenotype displayed by RA SFs, we speculate that a similar mechanism is displayed in the arthritic joint.
In the present study, invasiveness was successfully targeted by ST2825 on LPS-stimulated RA SF spheroids. In this regard, previous studies have reported suppression of LPS-induced migration and invasion of breast cancer cells after targeting MyD88-dependent signaling pathways, for instance, MyD88/GSK-3β/Snail and MyD88/NF-κB/Snail . Osteosarcoma cell migration was also abolished after treatment with ST2825 suggesting that MyD88 plays a major role in increasing the proliferation rate and reducing apoptosis . Increased expression of E2F1 in RA SFs has been strongly associated with pro-proliferative and pro-invasive effects through activation of the TGFβ-activated MAPK signaling. E2F1 seems to participate in modulating the expression of proangiogenic factors, vascular endothelial growth factor (VEGF) and hypoxia-inducible factor 1 subunit alpha (HIF1A), and cell cycle progression genes (CCNE2 and MYBL2), which are critical for synovial hyperplasia in RA . ST2825 has been previously reported as a repressor of MYBL2 expression, which is a central regulator of cell proliferation, cell survival, and differentiation in tumorigenesis . Our transcriptomic analysis has revealed that ST2825 downregulates those critical genes involved in migration, which suggests that the downstream effect of ST2825 will most likely result in suppressing the expression of cell cycle and proliferation-dependent genes on LPS-stimulated RA SFs.
Expression of several pro-inflammatory genes in LPS-stimulated RA SFs was downregulated after ST2825 treatment. The overexpression of IL1B, a key inflammatory mediator involved in RA pathogenesis, was significantly decreased by ST2825. Beyond the MyD88 dimerization inhibition, ST2825 significantly decreased the expression of MYD88. This effect might be attributed to the decreased level of its transcriptional activator E2F2, suggesting an additional regulatory mechanism for inflammation in RA orchestrated by ST2825 . The anti-inflammatory effects of ST2825 described in neuroinflammation suggest that this inhibitor is capable of suppressing pro-inflammatory mediators and oxidative stress by inhibiting the ROS/NLRP3/Caspase-1 signaling pathway . Taken together, these and our previous findings observed in RA PBMC, we speculated that ST2825 anti-inflammatory effects are able to prevent both local and systemic inflammation in RA patients.
Pain is a clinically relevant consequence of chronic inflammation in RA patients. In this regard, preclinical studies performed in animal models of RA and OA have shown that the blockade of IL-33 has the potential to reduce both pain and joint damage [64, 65]. At the mechanistic level, IL-33 may increase pain sensitivity via upregulating cyclooxygenase 2 (COX2) expression [66, 67]. CCL2-CCR2 interaction and downstream signaling were reported to excite sensory neurons and contribute to knee hyperalgesia in an experimental OA mouse model [68,69,70]. Our data show that ST2825 downregulates IL33, PTGS2 (encodes for COX2), and CCL2 expression in RA SFs. Thus, ST2825 may have an added therapeutic advantage as a pain medication for RA patients.
Studies in the BV2 neuronal cell line showed that ST2825 downregulated COX2 expression and decreased the production of reactive oxygen species . Our data suggest that ST2825 not only acts by downregulating COX2 expression but also exerts its mito-protective effects by upregulating genes encoding for mitochondrial proteins like TOMM7, TIMM8B, and NADH:ubiquinone oxidoreductase subunits, NDUFA1 and NDUFA2. TOMM7 is a mitochondrial translocase that plays an essential role in the translocation of nuclear-encoded mitochondrial proteins into mitochondria and regulates the coupling between oxygen consumption and ATP synthesis [71, 72]. Mutations in TOMM7 have resulted in syndromic short stature in humans [71, 73]. NDUFA1 and NDUFA2 are components of the mitochondrial respiratory complex I, which utilizes the NADH produced by Krebs Cycle to play a central role in energy metabolism [71, 74]. Interestingly, an expression quantitative trait locus (eQTL) study identified phosphorylation-associated single nucleotide polymorphisms (SNPs) linked with the NDUFA2 gene as significant RA risk alleles . Although the role of these mitochondrial proteins in RA SF pathology is not well understood, based on reports from other cell and tissue types, we speculate that their upregulation by ST2825 may counteract the mitochondrial dysfunction and restore energy homeostasis in RA SFs.
This study also elucidated the potential effect of ST2825 on targeting NF-κB-dependent mechanisms in RA SFs. Our results indicated that disrupting the Myddosome formation decreased the level of p65 and its activation. In addition, the activation of the IKK-related kinases, TAK1 and TBK1, was decreased by the effect of ST2825. Supporting our results, a previous study performed in bone marrow-derived macrophages (BMDMs) reported that LPS stimulation was able to induce the activation of TBK1 via the MyD88-dependent pathway, while in MyD88−/− BMDMs the activation was reduced . Additionally, TAK1 can initiate the activation of the IKK canonical complex by TLRs in a MyD88-dependent pathway in distinct cell types [76,77,78]. TAK1 inhibition has also been associated with decreased cell proliferation and cytokine secretion in malignant B cells . Taken together, these findings suggest that ST2825 mechanism-of-action could be dependent on inhibition of the IKK-related kinases, TBK1 and TAK1, and subsequent suppression of IKK canonical kinases (IKKα/β) converging on NF-κB p65 inactivation.
Our study did not consider the disease stage or disease activity score, which are important limitations. Additionally, the inclusion of treatments used for RA and OA patients may have revealed additional findings in SFs. Experiments with MYD88 siRNA could have provided proof for our speculation regarding the vital role of MyD88 driving the aggressive behavior of RA SFs; however, those experiments were not technically feasible due to the low transfection efficiency of MYD88 siRNA in RA SFs.
In summary, our results show that ST2825 inhibits multiple aspects of the RA SF aggressive features by arresting their cell cycle progression, inhibiting the production of inflammation and pain mediators, decreasing invasion, promoting apoptosis, and likely restoring mitochondrial homeostasis. Although this study is limited to in vitro experiments in a small sample size, our results provide a robust and unbiased justification to undertake future studies to evaluate the effects of ST2825 in a greater number of RA patient-derived SFs. Furthermore, our comprehensive RNA-seq analysis will serve as a resource to understand the importance of MyD88 dimerization in other inflammatory diseases. In the future, preclinical studies in animal models must be coupled with single-cell RNA sequencing to decipher the safety and efficacy of ST2825 on various subtypes of RA SFs and immune cell populations.
Availability of data and materials
All data underlying the findings of this article are fully available without restriction. RNA-seq data is deposited at the international public repository of Gene Expression Omnibus (GEO/NCBI) database as the study number GSE217012.
Peripheral blood mononuclear cells
Myeloid differentiation primary response 88
E2F family transcription factors
Pathogen-associated molecular patterns
Damage-associated molecular patterns
Tumor necrosis factor
- CXCL10 :
C-X-C motif chemokine ligand 10
- TLR3 :
Toll-like receptor 3
- TNFSF13B :
TNF superfamily member 13b
- CASP1 :
- IL1B :
- IL33 :
- ADAMTS4 :
ADAM metallopeptidase with thrombospondin type 1 motif 4
- IL1B :
Interleukin 1 beta
- TLR1 :
Toll-like receptor 1
- TLR3 :
Toll-like receptor 3
- CCL2 :
C-C motif chemokine ligand 2
- PTGS2 :
Prostaglandin-endoperoxide synthase 2
Interleukin receptor-associated kinases 1 and 4
Nuclear factor kappa light chain enhancer of activated B cells
B cell lymphoma 2
- BCL2L1 :
BCL2 like 1
- BBC3 :
BCL2-binding component 3
Nicotinamide adenine dinucleotide (reduced form)
Expression quantitative trait locus
Single nucleotide polymorphisms
Toll interleukin-1 receptor
Human dermal fibroblasts
RNA integrity number
Differentially expressed genes
Coordinated lysosomal expression and regulation
- VEGF :
Vascular endothelial growth factors
- HIF1A :
Hypoxia-inducible factor 1 subunit alpha
- MYBL2 :
Myb-related protein B
Inhibitor of nuclear factor kappa-B (IκB) kinase complex
TANK-binding kinase 1
Transforming growth factor-β (TGF-β)-activated kinase 1
Komatsu N, Takayanagi H. Mechanisms of joint destruction in rheumatoid arthritis - immune cell-fibroblast-bone interactions. Nat Rev Rheumatol. 2022;18(7):415–29.
Tu J, Huang W, Zhang W, Mei J, Zhu C. Two main cellular components in rheumatoid arthritis: communication between T cells and fibroblast-like synoviocytes in the joint synovium. Front Immunol. 2022;13:922111.
Bottini N, Firestein GS. Duality of fibroblast-like synoviocytes in RA: passive responders and imprinted aggressors. Nat Rev Rheumatol. 2013;9(1):24–33.
Nygaard G, Firestein GS. Restoring synovial homeostasis in rheumatoid arthritis by targeting fibroblast-like synoviocytes. Nat Rev Rheumatol. 2020;16(6):316–33.
Woods JM, Klosowska K, Spoden DJ, Stumbo NG, Paige DJ, Scatizzi JC, et al. A cell-cycle independent role for p21 in regulating synovial fibroblast migration in rheumatoid arthritis. Arthritis Res Ther. 2006;8(4):R113.
Xu N, Wang Y, Li D, Chen G, Sun R, Zhu R, et al. MDM4 overexpression contributes to synoviocyte proliferation in patients with rheumatoid arthritis. Biochem Biophys Res Commun. 2010;401(3):417–21.
Taranto E, Xue JR, Lacey D, Hutchinson P, Smith M, Morand EF, et al. Detection of the p53 regulator murine double-minute protein 2 in rheumatoid arthritis. J Rheumatol. 2005;32(3):424–9.
Hu Z, Chen Y, Zhu S, Feng X, Zhang B, Huang J. Sonic Hedgehog promotes proliferation and migration of fibroblast-like synoviocytes in rheumatoid arthritis via Rho/ROCK signaling. J Immunol Res. 2022;2022:3423692.
Mizoguchi F, Slowikowski K, Wei K, Marshall JL, Rao DA, Chang SK, et al. Functionally distinct disease-associated fibroblast subsets in rheumatoid arthritis. Nat Commun. 2018;9(1):789.
Zhang F, Wei K, Slowikowski K, Fonseka CY, Rao DA, Kelly S, et al. Defining inflammatory cell states in rheumatoid arthritis joint synovial tissues by integrating single-cell transcriptomics and mass cytometry. Nat Immunol. 2019;20(7):928–42.
Korsunsky I, Wei K, Pohin M, Kim EY, Barone F, Major T, et al. Cross-tissue, single-cell stromal atlas identifies shared pathological fibroblast phenotypes in four chronic inflammatory diseases. Med (N Y). 2022;3(7):481–518.e14.
Tolboom TC, van der Helm-Van Mil AH, Nelissen RG, Breedveld FC, Toes RE, Huizinga TW. Invasiveness of fibroblast-like synoviocytes is an individual patient characteristic associated with the rate of joint destruction in patients with rheumatoid arthritis. Arthritis Rheum. 2005;52(7):1999–2002.
Laragione T, Brenner M, Sherry B, Gulko PS. CXCL10 and its receptor CXCR3 regulate synovial fibroblast invasion in rheumatoid arthritis. Arthritis Rheum. 2011;63(11):3274–83.
Joosten LA, Abdollahi-Roodsaz S, Dinarello CA, O’Neill L, Netea MG. Toll-like receptors and chronic inflammation in rheumatic diseases: new developments. Nat Rev Rheumatol. 2016;12(6):344–57.
Elshabrawy HA, Essani AE, Szekanecz Z, Fox DA, Shahrara S. TLRs, future potential therapeutic targets for RA. Autoimmun Rev. 2017;16(2):103–13.
Sohn DH, Rhodes C, Onuma K, Zhao X, Sharpe O, Gazitt T, et al. Local Joint inflammation and histone citrullination in a murine model of the transition from preclinical autoimmunity to inflammatory arthritis. Arthritis Rheumatol. 2015;67(11):2877–87.
Manivel VA, Sohrabian A, Rönnelid J. Granulocyte-augmented chemokine production induced by type II collagen containing immune complexes is mediated via TLR4 in rheumatoid arthritis patients. Eur J Immunol. 2016;46(12):2822–34.
Lee JH, Cho ML, Kim JI, Moon YM, Oh HJ, Kim GT, et al. Interleukin 17 (IL-17) increases the expression of Toll-like receptor-2, 4, and 9 by increasing IL-1beta and IL-6 production in autoimmune arthritis. J Rheumatol. 2009;36(4):684–92.
Achek A, Yesudhas D, Choi S. Toll-like receptors: promising therapeutic targets for inflammatory diseases. Arch Pharm Res. 2016;39(8):1032–49.
Joosten LA, Koenders MI, Smeets RL, Heuvelmans-Jacobs M, Helsen MM, Takeda K, et al. Toll-like receptor 2 pathway drives streptococcal cell wall-induced joint inflammation: critical role of myeloid differentiation factor 88. J Immunol. 2003;171(11):6145–53.
Abdollahi-Roodsaz S, Joosten LA, Roelofs MF, Radstake TR, Matera G, Popa C, et al. Inhibition of Toll-like receptor 4 breaks the inflammatory loop in autoimmune destructive arthritis. Arthritis Rheum. 2007;56(9):2957–67.
Achek A, Shah M, Seo JY, Kwon HK, Gui X, Shin HJ, et al. Linear and rationally designed stapled peptides abrogate TLR4 pathway and relieve inflammatory symptoms in rheumatoid arthritis rat model. J Med Chem. 2019;62(14):6495–511.
Monnet E, Choy EH, McInnes I, Kobakhidze T, de Graaf K, Jacqmin P, et al. Efficacy and safety of NI-0101, an anti-toll-like receptor 4 monoclonal antibody, in patients with rheumatoid arthritis after inadequate response to methotrexate: a phase II study. Ann Rheum Dis. 2020;79(3):316–23.
Hernández-Jiménez M, Martín-Vílchez S, Ochoa D, Mejía-Abril G, Román M, Camargo-Mamani P, et al. First-in-human phase I clinical trial of a TLR4-binding DNA aptamer, ApTOLL: safety and pharmacokinetics in healthy volunteers. Mol Ther Nucleic Acids. 2022;28:124–35.
Sacre SM, Andreakos E, Kiriakidis S, Amjadi P, Lundberg A, Giddins G, et al. The Toll-like receptor adaptor proteins MyD88 and Mal/TIRAP contribute to the inflammatory and destructive processes in a human model of rheumatoid arthritis. Am J Pathol. 2007;170(2):518–25.
Zhan H, Chen H, Tang Z, Liu S, Xie K, Wang H. SIX1 attenuates inflammation and rheumatoid arthritis by silencing MyD88-dependent TLR1/2 signaling. Int Immunopharmacol. 2022;106:108613.
Wang S, Wang L, Wu C, Sun S, Pan JH. E2F2 directly regulates the STAT1 and PI3K/AKT/NF-κB pathways to exacerbate the inflammatory phenotype in rheumatoid arthritis synovial fibroblasts and mouse embryonic fibroblasts. Arthritis Res Ther. 2018;20(1):225.
Cho ML, Ju JH, Kim HR, Oh HJ, Kang CM, Jhun JY, et al. Toll-like receptor 2 ligand mediates the upregulation of angiogenic factor, vascular endothelial growth factor and interleukin-8/CXCL8 in human rheumatoid synovial fibroblasts. Immunol Lett. 2007;108(2):121–8.
Sikora KA, Bennett JR, Vyncke L, Deng Z, Tsai WL, Pauwels E, et al. Germline gain-of-function myeloid differentiation primary response gene-88 (MYD88) mutation in a child with severe arthritis. J Allergy Clin Immunol. 2018;141(5):1943–7.e9.
Ramirez-Perez S, Oregon-Romero E, Reyes-Perez IV, Bhattaram P. Targeting MyD88 downregulates inflammatory mediators and pathogenic processes in PBMC from DMARDs-naïve rheumatoid arthritis patients. Front Pharmacol. 2021;12:800220.
Loiarro M, Capolunghi F, Fantò N, Gallo G, Campo S, Arseni B, et al. Pivotal advance: inhibition of MyD88 dimerization and recruitment of IRAK1 and IRAK4 by a novel peptidomimetic compound. J Leukoc Biol. 2007;82(4):801–10.
Loiarro M, Ruggiero V, Sette C. Targeting the Toll-like receptor/interleukin 1 receptor pathway in human diseases: rational design of MyD88 inhibitors. Clin Lymphoma Myeloma Leuk. 2013;13(2):222–6.
Jones K, Angelozzi M, Gangishetti U, Haseeb A, de Charleroy C, Lefebvre V, et al. Human adult fibroblast-like synoviocytes and articular chondrocytes exhibit prominent overlap in their transcriptomic signatures. ACR Open Rheumatol. 2021;3(6):359–70.
Bhattaram P, Muschler G, Wixler V, Lefebvre V. Inflammatory cytokines stabilize SOXC transcription factors to mediate the transformation of fibroblast-like synoviocytes in arthritic disease. Arthritis Rheumatol. 2018;70(3):371–82.
Vinci M, Box C, Eccles SA. Three-dimensional (3D) tumor spheroid invasion assay. J Vis Exp. 2015;99:e52686.
Torre D, Lachmann A, Ma’ayan A. BioJupies: automated generation of interactive notebooks for RNA-Seq data analysis in the cloud. Cell Syst. 2018;7(5):556–61.e3.
Goedhart J, Luijsterburg MS. VolcaNoseR is a web app for creating, exploring, labeling and sharing volcano plots. Sci Rep. 2020;10(1):20560.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Jones K, Ramirez-Perez S, Niu S, Gangishetti U, Drissi H, Bhattaram P. SOX4 and RELA function as transcriptional partners to regulate the expression of TNF-responsive genes in fibroblast-like synoviocytes. Front Immunol. 2022;13:789349.
Ramírez-Pérez S, Hernández-Palma LA, Oregon-Romero E, Anaya-Macías BU, García-Arellano S, González-Estevez G, et al. Downregulation of inflammatory cytokine release from IL-1β and LPS-stimulated PBMC orchestrated by ST2825, a MyD88 dimerisation inhibitor. Molecules. 2020;25(18):4322.
Pucci B, Kasten M, Giordano A. Cell cycle and apoptosis. Neoplasia. 2000;2(4):291–9.
Zhu G, Cheng Z, Huang Y, Zheng W, Yang S, Lin C, et al. MyD88 mediates colorectal cancer cell proliferation, migration and invasion via NF-κB/AP-1 signaling pathway. Int J Mol Med. 2020;45(1):131–40.
Xu X, Yin Y, Tang J, Xie Y, Han Z, Zhang X, et al. Long non-coding RNA Myd88 promotes growth and metastasis in hepatocellular carcinoma via regulating Myd88 expression through H3K27 modification. Cell Death Dis. 2017;8(10):e3124.
Liu JH, Chen C, Li ZY, Zou ZM, Gao DC, Zhang X, et al. The MyD88 inhibitor TJ-M2010-2 suppresses proliferation, migration and invasion of breast cancer cells by regulating MyD88/GSK-3β and MyD88/NF-κB signalling pathways. Exp Cell Res. 2020;394(2):112157.
Shiratori E, Itoh M, Tohda S. MYD88 inhibitor ST2825 suppresses the growth of lymphoma and leukaemia cells. Anticancer Res. 2017;37(11):6203–9.
Xiang F, Ni Z, Zhan Y, Kong Q, Xu J, Jiang J, et al. Increased expression of MyD88 and association with paclitaxel resistance in breast cancer. Tumour Biol. 2016;37(5):6017–25.
Zhu Y, Huang JM, Zhang GN, Zha X, Deng BF. Prognostic significance of MyD88 expression by human epithelial ovarian carcinoma cells. J Transl Med. 2012;10:77.
Liang B, Chen R, Wang T, Cao L, Liu Y, Yin F, et al. Myeloid differentiation factor 88 promotes growth and metastasis of human hepatocellular carcinoma. Clin Cancer Res. 2013;19(11):2905–16.
Chen J, He J, Yang Y, Jiang J. An analysis of the expression and function of myeloid differentiation factor 88 in human osteosarcoma. Oncol Lett. 2018;16(4):4929–36.
Salcedo R, Worschech A, Cardone M, Jones Y, Gyulai Z, Dai RM, et al. MyD88-mediated signaling prevents development of adenocarcinomas of the colon: role of interleukin 18. J Exp Med. 2010;207(8):1625–36.
Fukuda K, Miura Y, Maeda T, Hayashi S, Kuroda R. Expression profiling of genes in rheumatoid fibroblast-like synoviocytes regulated by tumor necrosis factor-like ligand 1A using cDNA microarray analysis. Biomed Rep. 2019;1(1):1–5.
Lee C, Fernandez KJ, Alexandrou S, Sergio CM, Deng N, Rogers S, et al. Cyclin E2 promotes whole genome doubling in breast cancer. Cancers (Basel). 2020;12(8):2268.
Ren B, Cam H, Takahashi Y, Volkert T, Terragni J, Young RA, et al. E2F integrates cell cycle progression with DNA repair, replication, and G(2)/M checkpoints. Genes Dev. 2002;16(2):245–56.
Zhang R, Wang L, Pan JH, Han J. A critical role of E2F transcription factor 2 in proinflammatory cytokines-dependent proliferation and invasiveness of fibroblast-like synoviocytes in rheumatoid Arthritis. Sci Rep. 2018;8(1):2623.
Yue Y, Zhou T, Gao Y, Zhang Z, Li L, Liu L, et al. High mobility group box 1/toll-like receptor 4/myeloid differentiation factor 88 signaling promotes progression of gastric cancer. Tumour Biol. 2017;39(3):1010428317694312.
Liu Y, Li T, Xu Y, Xu E, Zhou M, Wang B, et al. Effects of TLR4 gene silencing on the proliferation and apotosis of hepatocarcinoma HEPG2 cells. Oncol Lett. 2016;11(5):3054–60.
Zhou Y, Chen X, Cao Z, Li J, Long H, Wu Y, et al. R848 is involved in the antibacterial immune response of golden pompano (Trachinotus ovatus) through TLR7/8-MyD88-NF-κB-signaling pathway. Front Immunol. 2020;11:617522.
Li Y, Xu JZ, Gu CX, Liu GL, Tian K. Carvacrol suppresses inflammatory responses in rheumatoid arthritis fibroblast-like synoviocytes. J Cell Biochem. 2019;120(5):8169–76.
Yao RB, Zhao ZM, Zhao LJ, Cai H. Sinomenine inhibits the inflammatory responses of human fibroblast-like synoviocytes via the TLR4/MyD88/NF-κB signaling pathway in rheumatoid arthritis. Pharmazie. 2017;72(6):355–60.
Wu R, Long L, Chen Q, Wu X, Zhu J, Zhou B, et al. Effects of Tim-3 silencing on the viability of fibroblast-like synoviocytes and lipopolysaccharide-induced inflammatory reactions. Exp Ther Med. 2017;14(3):2721–7.
Ainsworth RI, Hammaker D, Nygaard G, Ansalone C, Machado C, Zhang K, et al. Systems-biology analysis of rheumatoid arthritis fibroblast-like synoviocytes implicates cell line-specific transcription factor function. Nat Commun. 2022;13(1):6221.
Musa J, Aynaud MM, Mirabeau O, Delattre O, Grünewald TG. MYBL2 (B-Myb): a central regulator of cell proliferation, cell survival and differentiation involved in tumorigenesis. Cell Death Dis. 2017;8(6):e2895.
Zhang SS, Liu M, Liu DN, Shang YF, Wang YH, Du GH. ST2825, a small molecule inhibitor of MyD88, suppresses NF-κB activation and the ROS/NLRP3/cleaved caspase-1 signaling pathway to attenuate lipopolysaccharide-stimulated neuroinflammation. Molecules. 2022;27(9):2990.
He Z, Song Y, Yi Y, Qiu F, Wang J, Li J, et al. Blockade of IL-33 signalling attenuates osteoarthritis. Clin Transl Immunology. 2020;9(10):e1185.
Li Y, Fu Y, Chen H, Liu X, Li M. Blocking interleukin-33 alleviates the joint inflammation and inhibits the development of collagen-induced arthritis in mice. J Immunol Res. 2020;2020:4297354.
Li Y, Shi J, Qi S, Zhang J, Peng D, Chen Z, et al. IL-33 facilitates proliferation of colorectal cancer dependent on COX2/PGE(2). J Exp Clin Cancer Res. 2018;37(1):196.
Chen L, Song Z, Cao X, Fan M, Zhou Y, Zhang G. Interleukin-33 regulates the endoplasmic reticulum stress of human myometrium via an influx of calcium during initiation of labor. Elife. 2022;11:e75072.
Ishihara S, Obeidat AM, Wokosin DL, Ren D, Miller RJ, Malfait AM, et al. The role of intra-articular neuronal CCR2 receptors in knee joint pain associated with experimental osteoarthritis in mice. Arthritis Res Ther. 2021;23(1):103.
Miller RE, Tran PB, Das R, Ghoreishi-Haack N, Ren D, Miller RJ, et al. CCR2 chemokine receptor signaling mediates pain in experimental osteoarthritis. Proc Natl Acad Sci U S A. 2012;109(50):20602–7.
Miotla Zarebska J, Chanalaris A, Driscoll C, Burleigh A, Miller RE, Malfait AM, et al. CCL2 and CCR2 regulate pain-related behaviour and early gene expression in post-traumatic murine osteoarthritis but contribute little to chondropathy. Osteoarthritis Cartilage. 2017;25(3):406–12.
Young C, Batkovskyte D, Kitamura M, Shvedova M, Mihara Y, Akiba J, et al. A hypomorphic variant in the translocase of the outer mitochondrial membrane complex subunit TOMM7 causes short stature and developmental delay. HGG Adv. 2023;4(1):100148.
Hasson SA, Kane LA, Yamano K, Huang CH, Sliter DA, Buehler E, et al. High-content genome-wide RNAi screens identify regulators of parkin upstream of mitophagy. Nature. 2013;504(7479):291–5.
Garg A, Keng WT, Chen Z, Sathe AA, Xing C, Kailasam PD, et al. Autosomal recessive progeroid syndrome due to homozygosity for a TOMM7 variant. J Clin Invest. 2022;132:e156864.
Huttula S, Väyrynen H, Helisalmi S, Kytövuori L, Luukkainen L, Hiltunen M, et al. NDUFA1 p.Gly32Arg variant in early-onset dementia. Neurobiol Aging. 2022;114:113–6.
He P, Jiang F, Guo W, Guo YF, Lei SF, Deng FY. PhosSNPs-regulated gene network and pathway significant for rheumatoid arthritis. Hum Hered. 2021;86(1–4):10–20.
Clark K, Takeuchi O, Akira S, Cohen P. The TRAF-associated protein TANK facilitates cross-talk within the IkappaB kinase family during Toll-like receptor signaling. Proc Natl Acad Sci U S A. 2011;108(41):17093–8.
He A, Ji R, Shao J, He C, Jin M, Xu Y. TLR4-MyD88-TRAF6-TAK1 complex-mediated NF-κB activation contribute to the anti-inflammatory effect of V8 in LPS-induced human cervical cancer SiHa cells. Inflammation. 2016;39(1):172–81.
Ansell SM, Hodge LS, Secreto FJ, Manske M, Braggio E, Price-Troska T, et al. Activation of TAK1 by MYD88 L265P drives malignant B-cell growth in non-Hodgkin lymphoma. Blood Cancer J. 2014;4(2):e183.
We thank the support from the National Institutes of Health/National Institute of Arthritis, Musculoskeletal and Skin Disease grant (R01 AR070736), and Startup funds from the Department of Orthopaedics, Emory University School of Medicine assigned to PB. We also thank the support from the National Council of Science and Technology (CONACYT) grant (800977) assigned to SRP (CVU: 660472).
Ethics approval and consent to participate
The study was conducted according to the guidelines and approval of the Emory Institution Review Board.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Cell cycle analysis on hDF and OA SFs treated with ST2825. Representative histograms show cell cycle differences among different ST2825 concentrations and time points. G0/G1-phase cells, S-phase cells, and G2/M-phase cells are shown in the figure for hDF (A) and OA SFs (B). The y-axes represent the total count of cells, and the x-axes represent integrated intensity of DAPI measured by whole-well image cytometry. Fig. S2. Transcriptomic analysis revealed the upregulation of critical inflammatory mediators in RA SFs. (A) Description of OA and RA SFs used in RNA-seq experiments. (B) Volcano plot of up and downregulated genes in RA SFs compared OA SFs by 1.5 FC and p-value < 0.05 (C) Hierarchical clustering of the genes differentially expressed genes by Euclidean distance and centroid linkage method (D). Ingenuity Pathway Analysis (Qiagen) showing the top 5 canonical pathways predicted to be associated with the genes up and downregulated in RA SFs. (E) Heatmaps Log2 RPKM values of genes associated with canonical pathways identified in panel D. Fig. S3. Cell cycle analysis of RA-FLS treated with ST2825 at 48 and 72 h. Quantification of the percentage of RA SFs in various phases of the cell cycle upon ST2825 treatment at 48 (A) and 72 (B) h by imaging cytometry. One-way ANOVA and Dunnett’s multiple comparisons test were performed to determine statistical significance. *p<0.05, ****p<0.0001. Fig. S4. Cell cycle analysis on RA SFs treated with ST2825. Representative histograms show cell cycle differences among different ST2825 concentrations and time points. Apoptotic cells (yellow), G0/G1-phase cells (green), S-phase cells (blue), and G2/Mphase cells (purple) are shown in the figure. The y axes represent the total count of cells, and the x axes represents integrated intensity of DAPI measured by whole-well image cytometry. Fig. S5. Apoptosis analysis of RA-FLS treated with ST2825. Apoptosis was determined on RA SFs treated with 0, 5, and 10 μM of ST2825 at 0, 24, 48, and 72 h. No statistically significant differences were observed after 24 and 48 h of incubation with ST2825. The percentage of apoptotic cells significantly increased after 72 h of incubation with 10 μM of ST2825 (p=0.0476). Two-way ANOVA and Dunnett’s multiple comparisons test were used to determine statistical significane. *p<0.05. Fig. S6. Transcriptomic analysis identified DEGs and canonical pathways in RA SFs treated with ST2825. (A) Volcano plot of up and downregulated genes in RA SFs compared ST2825-treated RA SFs by 1.5 FC and p-value < 0.05. Enriched pathway analysis showing the top 10 canonical downregulated (B) and upregulated (C) pathways predicted to be associated with the genes up and downregulated in RA SFs treated with ST2825. Fig. S7. Transcriptomic analysis identified the effect of ST2825 on DEGs and canonical pathways in LPS-treated RA SFs. (A) Volcano plot of up and downregulated genes in LPS-stimulated RA SFs after treatment with ST2825 by 1.5 FC and p-value < 0.05. Enriched pathway analysis showing the top 10 canonical downregulated (B) and upregulated (C) pathways predicted to be associated with the genes up and downregulated in RA SFs. Fig S8. Raw data of western blot images. (A) Fig 5A and (B) Fig 5C. Table S1. Primer sequences for qRT-PCR. Table S2. Antibodies used for western blot.
Reactome Pathway Enrichment.
Reactome Pathway Enrichment LPS-stimulated RA.
About this article
Cite this article
Ramirez-Perez, S., Vekariya, R., Gautam, S. et al. MyD88 dimerization inhibitor ST2825 targets the aggressiveness of synovial fibroblasts in rheumatoid arthritis patients. Arthritis Res Ther 25, 180 (2023). https://doi.org/10.1186/s13075-023-03145-0