Bifurcations and multistability in a model of cytokine-mediated autoimmunity

This paper investigates the dynamics of immune response and autoimmunity with particular emphasis on the role of regulatory T cells (Tregs), T cells with different activation thresholds, and cytokines in mediating T cell activity. Analysis of the steady states yields parameter regions corresponding to regimes of normal clearance of viral infection, chronic infection, or autoimmune behavior, and the boundaries of stability and bifurcations of relevant steady states are found in terms of system parameters. Numerical simulations are performed to illustrate different dynamical scenarios, and to identify basins of attraction of different steady states and periodic solutions, highlighting the important role played by the initial conditions in determining the outcome of immune interactions.