Quantification of plasma cell dynamics using mathematical modelling

Plasma cells (PCs) are the main antibody-producing cells in humans. They are long-lived so that specific antibodies against either pathogens or vaccines are produced for decades. PC longevity is attributed to specific areas within the bone marrow micro-environment, the so-called ‘niche’, providing the cells with required growth and survival factors. With antigen encounters, e.g. infection or vaccination, new PCs are generated and home to the bone marrow where they compete with resident PCs for the niche. We propose a parametrized mathematical model describing healthy PC dynamics in the bone marrow. The model accounts for competition for the niche between newly produced PCs owing to vaccination and resident PCs. Mathematical analysis and numerical simulations of the model allow explanation of the recovery of PC homoeostasis after a vaccine-induced perturbation, and the fraction of vaccine-specific PCs inside the niche. The model enables quantification of the niche-related dynamics of PCs, i.e. the duration of PC transition into the niche and the impact of different rates for PC transitions into and out of the niche on the observed cell dynamics. Ultimately, it provides a potential basis for further investigations in health and disease.