Global sensitivity analysis of a mathematical model of acute inflammation identifies nonlinear dependence of cumulative tissue damage on host interleukin-6 responses

Abstract

The precise inflammatory role of the cytokine interleukin (IL)-6 and its utility as a biomarker or therapeutic target have been the source of much debate, presumably due to the complex pro- and anti-inflammatory effects of this cytokine. We previously developed a nonlinear ordinary differential equation (ODE) model to explain the dynamics of endotoxin (lipopolysaccharide; LPS)-induced acute inflammation and associated whole-animal damage/dysfunction (a proxy for the health of the organism), along with the inflammatory mediators tumor necrosis factor (TNF)-α, IL-6, IL-10, and nitric oxide (NO). The model was partially calibrated using data from endotoxemic C57Bl/6 mice. Herein, we investigated the sensitivity of the area under the damage curve (AUCD) to the 51 rate parameters of the ODE model for different levels of simulated LPS challenges using a global sensitivity approach called Random Sampling High Dimensional Model Representation (RS-HDMR). We explored sufficient parametric Monte Carlo samples to generate the variance-based Sobol' global sensitivity indices, and found that inflammatory damage was highly sensitive to the parameters affecting the activity of IL-6 during the different stages of acute inflammation. The AUCIL6 showed a bimodal distribution, with the lower peak representing healthy response and the higher peak representing sustained inflammation. Damage was minimal at low AUCIL6, giving rise to a healthy response. In contrast, intermediate levels of AUCIL6 resulted in high damage, and this was due to the insufficiency of damage recovery driven by anti-inflammatory responses and the activation of positive feedback sustained by IL-6. At high AUCIL6, damage recovery was interestingly restored in some population of simulated animals due to the NO-mediated anti-inflammatory responses. These observations suggest that the host's health status during acute inflammation depends in a nonlinear fashion on the magnitude of the inflammatory stimulus, on the host's propensity to produce IL-6, and on NO-mediated downstream responses.

Systems Modeling of Interactions between Mucosal Immunity and the Gut Microbiome during Clostridium difficile Infection

“…Consequently, the discovered marker could indicate the stage of infection or the need for or extension of treatment. Methods have been described based on the correlational evaluation of model quantities in virtual populations and global sensitivities in the analysis and discovery of potential biomarkers in the context of lipoprotein metabolism and acute inflammation, respectively [ 53 , 54 ]. The sensitivity analysis allows for the determination model parameters that are most impactful on immune response and infection outcome.…”

Neuro-Inflammation in Pediatric Traumatic Brain Injury—from Mechanisms to Inflammatory Networks

“…Most relationships or correlations established between cytokines in a network do not reflect this as they are done using correlation coefficients that assume a linear, monotonic relationship. This is demonstrably false, with a notable example being IL-6 displaying a bimodal distribution of expression following exposure to LPS in mice [174]. Instead of using data transformation to fit our parametric, linear modeling techniques like a Procrustean bed, we suggest it is worth exploring non-linear patterns of correlations using non-parametric coefficients that are less susceptible to outliers, such as Kendall's tau or Hoeffding's D [175].…”

Systems biology of IL-6, IL-12 family cytokines

“…A positive feedback driven by a damage-mediated increase in IL-6 maintains the inflammation. The model includes a representative but limited subset of mediators whereas other inflammatory triggers such as hemorrhagic shock might modulate the role of IL-6 [38].…”