Overall biological trends upon LMWF5A treatment
RNA sequencing results in an extensive, complicated dataset that can be difficult to interpret holistically, highlighting the importance of in silico analysis tools. To understand the overall biological trends present in a dataset of interest, IPA software employs an algorithm to construct a graphical summary that visualizes the most significant canonical pathways, upstream regulators, and biological functions from its core analysis and depicts their relationships to each other.
The RNA sequencing datasets from this study provide a snapshot of all genes that are regulated by LMWF5A in PBMC compared to saline vehicle control under three immunostimulatory conditions after 24 hours of treatment. The number of differentially expressed transcripts (adjusted p < 0.05) were 55, 64, and 139 for the three following comparisons, respectively: LMWF5A + LPS vs. saline + LPS, LMWF5A + LPS/IFNγ vs. saline + LPS/IFNγ, and LMWF5A + IL-4/IL-13 vs. saline + IL-4/IL-13. The graphical summaries of these datasets present the main biological targets and downstream activities of LMWF5A in this immune cell model. The graphical summary of the LPS/IFNγ data is shown as a representative network (Fig. 1); summaries for LPS and IL-4/IL-13 are shown in Supplemental Fig. 1a and b.
Several nodes are common between the three stimulation conditions. All nodes corresponding to biological targets represent molecules that are important in the immune response, such as cytokines, pattern recognition receptors, and inflammation-related transcription factors. Each graphical summary predicts inhibition of multiple interferons (IFNγ, IFNα2), interleukins (IL-1β, IL-6), TNFα, the hormone and cytokine prolactin (PRL), toll-like receptors (TLRs; TLR7), and transcription factors (interferon regulatory factor [IRF]7, STAT1). In addition, many of the implicated biological functions are associated with a decrease in pro-inflammatory activities. These include the general immune response of cells, cytotoxicity, monocyte/macrophage phagocytosis, integrin-mediated signaling, as well as the development, trafficking, and adhesion of various immune cell types. Furthermore, one pathway demonstrates predicted activation with LMWF5A treatment in IL-4/IL-13-stimulated PBMC, the macrophage-stimulating protein (MSP)- recepteur d'origine nantais (RON) signaling pathway, which has been described as an important negative regulator of inflammation due to its inhibition of pro-inflammatory cytokine production by macrophages [28].
Thus, the overall biological trend observed in PBMC stimulated with LPS, LPS/IFNγ, or IL-4/IL-13 in the presence of LMWF5A compared to saline control is an overarching suppression of inflammation, as numerous key pro-inflammatory regulators and functions are predicted to be inhibited and some anti-inflammatory regulators are predicted to be activated in this in silico analysis. This drug activity supports the use of LMWF5A in a range of clinical conditions that result from a hyperactive or chronic immune response.
Correlation of upstream regulators with LMWF5A
The identification of molecules that act in a similar or opposite manner to LMWF5A is valuable to further understand how LMWF5A acts on immune cells. The differential gene expression due to LMWF5A treatment in the RNA sequencing datasets was correlated with the activities of known substances using IPA upstream regulator analysis.
Upstream regulator analysis in IPA was utilized to connect LMWF5A-regulated genes with chemicals and proteins that share a subset of targets with analogous biological regulation (Fig. 2a). The significance of the correlation between the molecule in the IPA knowledge base and the uploaded dataset was measured by z-score, with a positive score indicating similar activity to the dataset and a negative score indicating an anti-match in activity. A comparison analysis of upstream regulators was performed to visualize regulators that are relevant to all stimulation conditions simultaneously. Figure 2b lists the top upstream regulators that correlate with LMWF5A-regulated transcripts, while Fig. 2c lists the top regulators that are inversely correlated.
Interestingly, the upstream regulator with the highest total z-score was the corticosteroid dexamethasone (Fig. 2b), which acts as a glucocorticoid receptor agonist and affects a multitude of indirect downstream targets [29]. The primary action of dexamethasone is to suppress inflammatory cell activity and inflammatory mediator levels, which mirrors the biological trends predicted by IPA for LMWF5A in this study as well as previously published work on LMWF5A. Dexamethasone, as a synthetic steroid, has multiple indications, including inflammatory conditions (such as rheumatoid and psoriatic arthritis, systemic lupus erythematosus, and Crohn’s disease), multiple sclerosis, cerebral edema, shock, and allergies, amongst others [30]. Importantly, dexamethasone has also been administered as an intra-articular injection for joint inflammation [30] and has been recently studied as a potent anti-inflammatory for the treatment of COVID-19 [22]; both conditions are currently being investigated in LMWF5A clinical trials.
The other top five most positively correlated upstream regulators with LMWF5A (Fig. 2b) include sirtuin 1 (SIRT1), prostaglandin E receptor 4 (PTGER4), filgrastim, and SB203580. Like dexamethasone, each of these compounds has demonstrated anti-inflammatory activity.
SIRT1 is a NAD-dependent protein deacetylase that negatively regulates inflammation by altering cytokine levels and immune cell recruitment and activation by deacetylating and suppressing transcription factors, including NFκB [31]. Similarly, LMWF5A has also been proven to modulate transcription factor activity to decrease pro-inflammatory cytokine levels. Its ability to reduce the release of both TNFα and IL-1β has been linked to its effects on the NFκB-repressing PPARγ and AhR pathways in LPS-stimulated PBMC, and it has also been shown to prevent NFκB reporter activity [8].
PTGER4 (or EP4) is a transmembrane, G-coupled protein receptor that becomes activated when bound to the cyclooxygenase (COX) pathway product PGE2, which is induced during inflammation. Although the role of PGE2 is pleotropic, it exerts its anti-inflammatory effects via PTGER4. PTGER4 binds to EP4 receptor-associated protein, which in turn, reduces the phosphorylation and increases the stability of p105, an inhibitor of NFκB and mitogen-activated protein/extracellular signal-regulated kinase (MEK) [32]. Overall, PTGER4 downregulates inflammation by modulating macrophage cytokine and chemokine secretion as well as T cell proliferation, differentiation, and cytokine production [33]. LMWF5A has been associated with an upregulation of the COX pathway and its products, including the PTGER4 ligand PGE2, in PBMC and primary human osteoarthritic cells [7, 9]. With respect to the COX pathway, the activity of LMWF5A may be unique to the anti-inflammatory drug class, as most inhibit both cytokine and prostaglandin release while LMWF5A inhibits cytokine release but promotes prostaglandin release. Inhibition of the COX/prostaglandin pathway can result in harmful side effects [34], and the stimulation of this pathway with concomitant inhibition of cytokine production by LMWF5A may, conversely, promote anti-inflammation, resolution, and healing.
Filgrastim, recombinant human granulocyte colony stimulating factor, is commonly used clinically as a complement to chemotherapy due to its ability to stimulate granulocyte production, thus preventing low white blood cell counts [35]. In addition, this protein has confirmed anti-inflammatory properties in vivo with respect to the cytokine response and has been suggested as a potential treatment for chronic inflammatory conditions; for example, LPS-induced cytokine release has been shown to be attenuated in healthy human volunteers treated with filgrastim [36, 37].
Finally, SB203580 is a potent p38 mitogen-activated protein kinase (MAPK) inhibitor, with strong effects on cytokine production [38]. The p38 MAPK signaling pathway is critical to the regulation of many cellular processes, including inflammation, as it is activated in response to inflammatory mediators and other stress-related molecules and acts as a major regulator of cytokine production [39]. Overall, the top upstream regulators that significantly match LMWF5A activity have been empirically proven to downregulate components critical to inflammation, including pro-inflammatory cytokines and transcription factors.
In contrast, the top five upstream regulators that do not match the actions of LMWF5A and have the lowest total z-scores are lipopolysaccharide (LPS), IFNγ, STAT1, poly rI:rC-RNA, and IFNα2 (Fig. 2c). This list consists of established pro-inflammatory pathogen- and damage-associated molecular patterns (P/DAMP), transcription factors, and cytokines. During the inflammatory cascade, recognition of PAMPs and DAMPs by pattern recognition receptors, like TLRs, results in signal transduction to turn on multiple transcription factors (NF-κB, MAPK, STAT) that increase the expression of key cytokines that promote inflammation [1].
LPS and poly rI:rC-RNA are both categorized as PAMPs and function as signals of infection to initiate immune signaling via TLRs. LPS is a major component of Gram-negative bacterial cell walls. It acts as a PAMP upon bacterial infection in vivo and is widely utilized in vitro to stimulate the TLR4-induced immune response, which can occur upon recognition of either pathogens or endogenous molecules that are released upon tissue damage [40]. LPS was used in two of the three conditions in this study to stimulate an immune response in PBMC, and its ranking as the most negatively correlated upstream regulator to LMWF5A highlights the fact that LMWF5A strongly counteracts the inflammatory outcomes of TLR4 signaling. Moreover, the anti-inflammatory effects of LMWF5A on cytokine release and transcription factor activity have been extensively studied using cells treated with LPS as a TLR4 stimulant [4, 8, 41]. poly rI:rC-RNA is a dsRNA mimic that can simulate infection with a dsRNA virus and activate TLR3 [42], suggesting that LMWF5A may also offset the actions of other TLR-driven pathways; the relationships between LMWF5A and other TLRs are currently under investigation.
As additional opposing upstream regulators, the cytokines IFNγ and IFNα2 represent both classes of interferons, type II and I, respectively. They are secreted upon viral infection to limit viral replication and regulate the subsequent immune response [43]. The identification of these IFNs as highly ranked, negatively correlated upstream regulators to LMWF5A emphasizes the potential benefits of LMWF5A for the treatment of viral infections that involve a hyperimmune response or cytokine storm, including COVID-19 [44]. The Janus kinase (JAK)-STAT signaling pathway is the most studied IFN-related transcription factor pathway, but IFNs also activate other signaling cascades, including p38 MAPK and phosphatidylinositol 3'-kinase, to exert their anti-viral and pro-inflammatory effects [45].
Because both type I and II IFNs activate STAT complexes, it is not surprising that STAT1 is also part of this list of opposing upstream regulators. Upon recognition of IFNs, a variety of interleukins, or other cytokines, STAT1 is phosphorylated, mainly by JAK kinases, and activated to drive a pro-inflammatory cascade [46]. LMWF5A was also previously demonstrated to inhibit the ability of STAT1 to bind its cognate DNA sequence in LPS-stimulated PBMC, implicating regulation of this transcription factor as part of the LMWF5A mode of action [8]. Inhibition of the JAK-STAT pathway has been suggested as a therapeutic strategy for the treatment of inflammatory conditions, including rheumatoid arthritis (RA), psoriasis, inflammatory bowel disease [47] as well as COVID-19 [48]. In summary, the analysis of inversely correlated upstream regulators of LMWF5A represent pro-inflammatory factors, emphasizing the anti-inflammatory activity of LMWF5A.
Comparison of dexamethasone and LMWF5A targets
In the three different immunostimulatory conditions used to evaluate the effect of LMWF5A on PBMC, dexamethasone was determined to be the most significant positively correlated upstream regulator to LWMF5A. Due to this finding and the longstanding and widespread utility of dexamethasone in the clinic, a detailed comparison of its targets to the LMWF5A datasets was performed.
The IPA Ingenuity Knowledge Base is an extensive database that is comprised of empirically proven interactions and relationships between genes, proteins, drugs, etc. By querying the Knowledge Base for factors both linked to dexamethasone and differentially expressed upon LMWF5A treatment, it was determined that LMWF5A regulates 29, 26, and 42 established dexamethasone targets for the LPS/IFNγ, LPS, and IL-4/IL-13 stimulation conditions, respectively. These common targets represent a large proportion of each LMWF5A dataset, with 53%, 41%, and 30% of the transcripts significantly regulated by LMWF5A are also targeted by dexamethasone. However, based on IPA computational modeling, directional regulation by LMWF5A of 7, 10, and 17 of the common targets were inconsistent with that of dexamethasone for the LPS/IFNγ, LPS, and IL-4/IL-13 stimulation conditions, respectively. Nevertheless, overall IPA modeling predicted significantly similar regulation of downstream molecules by dexamethasone and LMWF5A based on directionally correlated transcripts. The common targets between dexamethasone and LMWF5A are shown for the LPS/IFNγ stimulation conditions in a representative figure (Fig. 3). The common targets for the LPS and IL-4/IL-13 conditions are presented in Supplemental Fig. 2a and b.
In Fig. 3, 29 transcripts are presented that were regulated by LMWF5A in LPS/IFNγ-stimulated PBMC at 24 h that are also known targets of dexamethasone. Of these transcripts, 26 were decreased by LMWF5A (green) and 3 were increased (red). Further, of the 29 transcripts, 19 were directionally consistent with the action of dexamethasone (orange or blue lines) and 7 were directionally inconsistent (yellow lines); 3 dexamethasone targets have unpredicted directional regulation (gray lines).
An example of a target that is regulated by LMWF5A in a direction consistent with dexamethasone is C–C motif chemokine ligand 2 (CCL2), also called monocyte chemoattractant-1. Upon recognition of inflammatory stimuli, CCL2 expression is induced, and this chemokine drives migration of immune cells, particularly monocytes, to the site of infection or tissue injury [49]. Dysregulated, increased CCL2 expression is linked to the pathology of many diseases, including heart failure, RA, and diabetes, due its overpromotion of immune cell infiltration and downstream pro-inflammatory effects [50]. Serum CCL2 has been demonstrated to be increased in OA patients versus healthy controls, suggesting its importance to OA pathogenesis [51]. Some studies have reported increased CCL2 levels in OA synovial fluid as well, and CCL2 has been linked to OA-associated pain in addition to factors influencing cartilage catabolism [52]. With respect to COVID-19, CCL2, along with many other inflammatory cytokines, contributes to the cytokine storm that occurs during the body’s dysregulated response to SARS-CoV2 infection and its level has been correlated with increased disease severity [48]. Thus, the similar activities of both dexamethasone and LMWF5A on CCL2 should provide benefit to patients with CCL2-related diseases, including OA and COVID-19.
Conversely, an example of a target that is regulated by LMWF5A in a direction inconsistent with dexamethasone is cathepsin B (CTSB). CTSB is member of the cathepsin family of cysteine proteases, which are localized in the lysosome [53]. It is well studied in the context of cancer [54] but has also been implicated in cartilage degradation and OA pathogenesis due to its proteolytic activity of extracellular matrix components and its ability to promote the activity of other proteases [55]. In fact, the enzymatic activity of CTSB was demonstrated to be enhanced in vitro in stimulated primary chondrocytes and in vivo in cartilage, serum, and synovial fluid from OA patients, in which CTSB activity levels were associated with OA disease severity [56]. In the IPA upstream regulator analysis, the Knowledge Database records CTSB as being upregulated by dexamethasone, supporting the cartilage-degrading effects observed with corticosteroids, while the differential expression analysis showed that CTSB was downregulated by LMWF5A. This and other differences may be important distinctions between these two anti-inflammatory treatments. While it is possible to explore each specific directionally consistent or inconsistent target, the present in silico analysis allows for a broad understanding of overall LMWF5A activity and its comparison to dexamethasone as well as the interplay between common targets and their downstream pathways.
Further comparison of dexamethasone and LMWF5A revealed 71 total common targets (Fig. 4). Of the 71 common targets, only five transcripts were found to be affected by LMWF5A under all three stimulation conditions: β-actin (ACTB), γ-actin (ACTG1), myristoylated alanine-rich C-kinase substrate (MARCKS), and tryptophanyl-tRNA synthetase 1 (WARS1), all associated with the cytoskeleton, and cytochrome b-245 b chain (CYBB; also called NADPH oxidase 2), which functions to produce reactive oxygen species to eliminate pathogens during infection and is linked to the development of acute respiratory distress syndrome [57, 58]. Ten of the targets were shared between the datasets containing LPS, and three were shared between the IL-4/IL-13 and each of the LPS and LPS/IFNγ datasets. Moreover, 50 of the 71 common targets were unique to one stimulation condition, highlighting the wide-ranging effects of LMWF5A as well as its different activities depending on the inflammatory milieu. This result suggests that LMWF5A treatment may be beneficial in multiple inflammatory environments.
The differences in common targets between the datasets points to a limitation of the current study. This analysis solely examined the effects of LMWF5A in one cell model at one time point, albeit with three immunostimulation conditions. Additional timepoints or cell types would provide even further insight into the mode of action of LMWF5A and may identify more overlapping targets. For example, due to the known effects of LMWF5A on cytokine release from PBMC, an earlier timepoint may have captured more differentially expressed cytokine RNAs as well as any temporal differences in cytokine RNA levels due to LMWF5A treatment. However, extensive information has been gained with the large datasets in this investigation.