Large-scale computational modelling of the M1 and M2 synovial macrophages in rheumatoid arthritis

Naouel Zerrouk,Rachel Alcraft,Benjamin A. Hall,Franck Augé,Anna Niarakis

DOI: https://doi.org/10.1038/s41540-024-00337-5

2024-01-26

npj Systems Biology and Applications

Abstract:Abstract Macrophages play an essential role in rheumatoid arthritis. Depending on their phenotype (M1 or M2), they can play a role in the initiation or resolution of inflammation. The M1/M2 ratio in rheumatoid arthritis is higher than in healthy controls. Despite this, no treatment targeting specifically macrophages is currently used in clinics. Thus, devising strategies to selectively deplete proinflammatory macrophages and promote anti-inflammatory macrophages could be a promising therapeutic approach. State-of-the-art molecular interaction maps of M1 and M2 macrophages in rheumatoid arthritis are available and represent a dense source of knowledge; however, these maps remain limited by their static nature. Discrete dynamic modelling can be employed to study the emergent behaviours of these systems. Nevertheless, handling such large-scale models is challenging. Due to their massive size, it is computationally demanding to identify biologically relevant states in a cell- and disease-specific context. In this work, we developed an efficient computational framework that converts molecular interaction maps into Boolean models using the CaSQ tool. Next, we used a newly developed version of the BMA tool deployed to a high-performance computing cluster to identify the models’ steady states. The identified attractors are then validated using gene expression data sets and prior knowledge. We successfully applied our framework to generate and calibrate the M1 and M2 macrophage Boolean models for rheumatoid arthritis. Using KO simulations, we identified NFkB, JAK1/JAK2, and ERK1/Notch1 as potential targets that could selectively suppress proinflammatory macrophages and GSK3B as a promising target that could promote anti-inflammatory macrophages in rheumatoid arthritis.

mathematical & computational biology

What problem does this paper attempt to address?

The paper attempts to address the issue of functional imbalance between M1 and M2 synovial macrophages in Rheumatoid Arthritis (RA). Specifically, the paper aims to develop an efficient computational framework to convert molecular interaction maps into Boolean models and use these models to identify and validate potential therapeutic targets to selectively inhibit pro-inflammatory M1 macrophages and promote anti-inflammatory M2 macrophages.

### Background
1. **Rheumatoid Arthritis**: RA is a complex inflammatory autoimmune disease characterized by persistent inflammatory responses in the synovial tissue, leading to joint pain, swelling, stiffness, and subsequent extra-articular manifestations such as gastrointestinal and cardiovascular diseases.
2. **Role of Macrophages**: Macrophages play a key role in the pathology of RA. Depending on their phenotype (M1 or M2), they can be involved in either the initiation or resolution of inflammation.
- **M1 Macrophages**: Responsible for the excessive production of inflammatory cytokines and matrix-degrading enzymes, leading to cartilage destruction, and attracting and activating pro-inflammatory T cells.
- **M2 Macrophages**: Alleviate inflammation by producing anti-inflammatory cytokines (such as IL-10 and TGF-β), maintaining tissue homeostasis and repair, and activating regulatory T cell functions.
3. **Imbalance of M1/M2 Ratio**: In RA patients, the proportion of M1 macrophages is higher than that of M2 macrophages, which may lead to persistent inflammation and disease progression.
4. **Limitations of Existing Therapies**: Currently, there are no therapies specifically targeting macrophages in clinical practice. Existing treatments mainly focus on reducing inflammation and alleviating symptoms but have various adverse effects.

### Research Objectives
1. **Develop Computational Framework**: Construct an efficient computational framework to convert molecular interaction maps into Boolean models and use high-performance computing clusters to identify the stable states (attractors) of the models.
2. **Validate Model Behavior**: Use gene expression datasets and prior knowledge to validate the identified attractors.
3. **Identify Therapeutic Targets**: Identify potential therapeutic targets that can selectively inhibit M1 macrophages and promote M2 macrophages through simulated knockout (KO) experiments.

### Main Findings
1. **Model Generation and Calibration**: Successfully generated and calibrated Boolean models for RA M1 and M2 macrophages.
2. **Potential Therapeutic Targets**: Identified NFkB, JAK1/JAK2, and ERK1/Notch1 as potential therapeutic targets through KO simulations, which can selectively inhibit pro-inflammatory M1 macrophages. Additionally, GSK3B was identified as a potential target that can promote anti-inflammatory M2 macrophages.

### Conclusion
The study successfully generated and validated Boolean models for RA M1 and M2 macrophages by developing an efficient computational framework and identified potential therapeutic targets, providing new insights and methods for the treatment of RA.