A mathematical modelling of rheumatoid arthritis
Abstract
We developed a mathematical model to describe the inflammatory dynamics during arthritis development and resolution, focusing on the interaction between three key immune cell populations within the joint tissue, namely neutrophils (N), sublining synovial tissue-resident macrophages (SL-TRM, T), and circulating monocytes that are recruited into the inflamed joint and differentiate into monocyte-derived macrophages (MoMac, M). An inflammatory signal (S), representing the response to an external arthritogenic stimulus, drives the recruitment and activation of neutrophils and MoMac, while SL-TRM exert a regulatory and protective role. We propose a minimal ordinary differential equation model of acute inflammatory response into a joint that couples an inflammatory signal S to the three cell populations (N, T, M) and incorporates recruitment, stimulation, nonlinear self-limitation and mutual regulation. We show that the equilibrium structure depends on three adimensional parameters A, B, C and derive a cubic equation characterizing nontrivial steady states. We perform an analytical study of equilibria and their stability, and identify a bifurcation at A = 1 separating a disease-free regime from a sustained inflammatory state. Numerical simulations illustrate typical transient dynamics under various scenarios and show the ability of the model to capture the transition from acute to chronic inflammation, as well as the resolution of inflammation under certain conditions. In conclusion, our model provides a framework to explore mechanisms underlying acute versus sustained joint inflammation and offers a basis for parameter estimation from experimental data, as well as future extensions toward evaluation of therapeutic strategies.
Related
- Methodological quality of systematic reviews and meta-analyses of rheumatoid arthritis pharmacotherapy: A protocol for a meta-epidemiological study using AMSTAR 2
- Fitness advantages, prediction of resistances and protocol design: exploiting the synergy between mathematical modelling and statistics
- Mathematical modelling of a multicluster open system: transfer-operator recovery from aggregate data under a least-action prior
- Mathematical folklore
- Digital Learning Platforms and Mathematical Problem-Solving Ability: The Moderating Role of Self-Regulated Learning
BibTeX
@misc{azerad2026mathematical,
title = {{A mathematical modelling of rheumatoid arthritis}},
author = {Pascal Azerad and Oana Iosifescu and Andrea Parmeggiani and Olivier Bortolotti and Florence Apparailly and Gabriel Courties},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.15543},
eprint = {2609.15543},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.15543},
}
RIS
TY - UNPB
TI - A mathematical modelling of rheumatoid arthritis
AU - Pascal Azerad
AU - Oana Iosifescu
AU - Andrea Parmeggiani
AU - Olivier Bortolotti
AU - Florence Apparailly
AU - Gabriel Courties
PY - 2026
DA - 2026/09/14
PB - arXiv
DO - 10.48550/arxiv.2609.15543
UR - https://arxiv.org/abs/2609.15543
AB - We developed a mathematical model to describe the inflammatory dynamics during arthritis development and resolution, focusing on the interaction between three key immune cell populations within the joint tissue, namely neutrophils (N), sublining synovial tissue-resident macrophages (SL-TRM, T), and circulating monocytes that are recruited into the inflamed joint and differentiate into monocyte-derived macrophages (MoMac, M). An inflammatory signal (S), representing the response to an external arthritogenic stimulus, drives the recruitment and activation of neutrophils and MoMac, while SL-TRM exert a regulatory and protective role. We propose a minimal ordinary differential equation model of acute inflammatory response into a joint that couples an inflammatory signal S to the three cell populations (N, T, M) and incorporates recruitment, stimulation, nonlinear self-limitation and mutual regulation. We show that the equilibrium structure depends on three adimensional parameters A, B, C and derive a cubic equation characterizing nontrivial steady states. We perform an analytical study of equilibria and their stability, and identify a bifurcation at A = 1 separating a disease-free regime from a sustained inflammatory state. Numerical simulations illustrate typical transient dynamics under various scenarios and show the ability of the model to capture the transition from acute to chronic inflammation, as well as the resolution of inflammation under certain conditions. In conclusion, our model provides a framework to explore mechanisms underlying acute versus sustained joint inflammation and offers a basis for parameter estimation from experimental data, as well as future extensions toward evaluation of therapeutic strategies.
ER -