Control of human immune response function by T-cell population fluctuation and relaxation dynamics

Clinical studies have indicated that in malignant surveillances fluctuations in the population of certain effector T-cell repertoire become suppressed. Motivated by such observations and in an attempt to quantify adaptive human response to pathogens, we define an immune response function (IMRF) in terms of mean square fluctuations of T-cell concentrations. We employ a recently developed kinetic model of T-cell regulation that contains the essential immunosuppressive effects of vitamin-D. We employ Gillespie algorithm to make the first study of fluctuations along the stochastic trajectories. This fluctuation-based IMRF can differentiate responses of different individuals after pathogenic incursion both under healthy and disease conditions. We find that relative fluctuations in T-cells (and hence IMRF) are different in strongly regulated (malignant prone) and weakly regulated (autoimmune prone) regions. The cross-over from one steady state (weakly regulated) to the other (strongly regulated) is accompanied by a divergence-like growth in the fluctuation of both the effector and regulatory T-cell concentration over a wide range of pathogenic stimulation, displaying a dynamical phase transition like behavior. The growth in fluctuation in this desired immune response regime is found to arise from an intermittent fluctuation between regulatory and effector T-cells that results in a bimodal distribution of population of each, indicating bistability. The signature of intermittent behavior is further confirmed by calculating the power spectrum of the corresponding fluctuation of time correlation function. The calculated time correlation functions of fluctuations show that the slow fluctuation causes the bistabilty in healthy state. Thus, in diseases diagnosis process, such steady state response parameters can provide immense information which might become helpful to define an immune status.

💡 Research Summary

The paper introduces a novel quantitative metric, the Immune Response Function (IMRF), to capture the dynamic variability of T‑cell populations during immune surveillance. Traditional immunological assessments often rely on average cell counts or cytokine levels, which overlook the stochastic fluctuations that can be critical for distinguishing healthy from diseased states. By defining IMRF as the mean‑square fluctuation of effector (Teff) and regulatory (Treg) T‑cell concentrations normalized by their mean numbers, the authors provide a measure that directly reflects the system’s intrinsic noise and its susceptibility to external perturbations.

To explore the behavior of IMRF, the authors employ a recently developed kinetic model of T‑cell regulation that explicitly incorporates the immunosuppressive role of active vitamin‑D. The model consists of five interacting species: inactive T‑cells, activated Teff, activated Treg, active vitamin‑D, and a pathogenic antigen. Vitamin‑D promotes Treg activation while simultaneously inhibiting Teff proliferation, creating an asymmetric feedback loop that can generate two distinct steady‑states: a “strongly regulated” state (characteristic of cancer‑prone individuals) and a “weakly regulated” state (characteristic of autoimmune‑prone individuals).

Because the system operates at low cell numbers where stochastic effects dominate, the authors simulate its dynamics using the Gillespie stochastic simulation algorithm. They run tens of thousands of independent trajectories to obtain reliable statistics on both mean values and fluctuations. The simulations reveal that, as the pathogen stimulation strength (k) and vitamin‑D concentration (VD) are varied, the IMRF exhibits a sharp, divergence‑like increase at a specific region of parameter space. This region corresponds to the crossover between the weakly and strongly regulated steady‑states. Near this crossover, both Teff and Treg populations display dramatically amplified fluctuations, indicating a dynamical phase‑transition‑like behavior.

Further analysis shows that the probability distributions of Teff and Treg counts become bimodal in the crossover region, confirming the presence of bistability: the system can reside in either the high‑Teff/low‑Treg (weak regulation) or low‑Teff/high‑Treg (strong regulation) mode, and stochastic switching between these modes occurs intermittently. Time‑correlation functions of the fluctuations decay slowly, exhibiting long‑range correlations, and their Fourier transforms (power spectra) display a pronounced 1/f‑type low‑frequency component. These signatures are hallmarks of intermittent dynamics and support the interpretation that the system’s slow fluctuations underlie the observed bistability.

Clinically, the authors argue that IMRF can serve as a powerful biomarker. In strongly regulated (cancer‑prone) individuals, the IMRF remains low because fluctuations are suppressed, reflecting an immune system that is less responsive to emerging threats. Conversely, in weakly regulated (autoimmune‑prone) individuals, the IMRF is elevated, mirroring heightened variability and a propensity for over‑active immune responses. Thus, measuring IMRF could help differentiate between these immunological phenotypes, guide early diagnosis, and monitor therapeutic interventions such as vitamin‑D supplementation.

The study acknowledges several limitations. The current model captures only the vitamin‑D mediated feedback and does not include the full cytokine network, spatial heterogeneity, or cell migration, all of which can influence immune dynamics. Moreover, the translation of simulated fluctuations to experimentally measurable quantities (e.g., flow cytometry‑derived cell count variance) remains to be validated. Future work is proposed to collect longitudinal patient data, calibrate the model against real‑world measurements, and integrate IMRF into machine‑learning frameworks for personalized immune profiling.

In summary, by focusing on stochastic fluctuations rather than static averages, the authors demonstrate that the immune system exhibits a dynamical phase transition between distinct regulatory regimes. The IMRF, grounded in mean‑square T‑cell fluctuations, provides a novel, information‑rich descriptor of immune status that could augment existing diagnostic tools and deepen our understanding of immune regulation in health and disease.