Fluctuation theory of immune response: A statistical mechanical approach to understand pathogen induced T-cell...

In this period of intense interest in human immunity, we attempt here to quantify the immune response against pathogen invasion through T-cell population dynamics. Borrowing concepts from equilibrium statistical mechanics, we introduce a new description of the immune response function (IMRF) in terms of fluctuations in the population number of relevant biological cells (effector and regulatory T-cells). We use a coarse-grained chemical reaction network model (CG-CRNM) to calculate the number fluctuations and show that the response function derived as such can, indeed, capture the crossover observed in a T-cell driven immune response. We employ the network model to learn the effect of vitamin-D as an immunomodulator. We solve our CG-CRNM using a stochastic Gillespie algorithm. Depending on the effector T-cell concentration, we can classify immune regulation regimes into three categories: weak, strong, and moderate. The IMRF is found to behave differently in these three regimes. A damped cross-regulatory behavior found in the dynamics of effector and regulatory T-cell concentration in the diseased states correlates well with the same found in a cohort of patients with specific malignancies and autoimmune diseases. Importantly, the crossover from the weakly regulated steady state to the other (the strongly regulated) is accompanied by a divergence-like growth in the fluctuation of both the effector and the regulatory T-cell concentration, characteristic of a dynamic phase transition. We believe such steady-state IMRF analyses could help not only to phase-separate different immune stages but also aid in the valuable connection between autoimmunity, optimal vitamin-D, and consequences of immunosuppressive stress and malignancy.