A reduced mathematical model of the acute inflammatory response: I. Derivation of model and analysis of anti-inflammation

Modelling the continuum of macrophage phenotypes and their role in inflammation

Macrophages are a type of white blood cell that play a significant role in determining the inflammatory response associated with a wide range of medical conditions. They are highly plastic, having the capacity to adopt numerous polarisation states or 'phenotypes' with disparate pro- or anti-inflammatory roles. Many previous studies divide macrophages into two categorisations: M1 macrophages are largely pro-inflammatory in nature, while M2 macrophages are largely restorative. However, there is a growing body of evidence that the M1 and M2 classifications represent the extremes of a much broader spectrum of phenotypes, and that intermediate phenotypes can play important roles in the progression or treatment of many medical conditions. In this article, we present a model of macrophage dynamics that includes a continuous description of phenotype, and hence incorporates intermediate phenotype configurations. We describe macrophage phenotype switching via nonlinear convective flux terms that scale with background levels of generic pro- and anti-inflammatory mediators. Through numerical simulation and bifurcation analysis, we unravel the model's resulting dynamics, paying close attention to the system's multistability and the extent to which key macrophage-mediator interactions provide bifurcations that act as switches between chronic states and restoration of health. We show that interactions that promote M1-like phenotypes generally result in a greater array of stable chronic states, while interactions that promote M2-like phenotypes can promote restoration of health. Additionally, our model admits oscillatory solutions reminiscent of relapsing-remitting conditions, with macrophages being largely polarised toward anti-inflammatory activity during remission, but with intermediate phenotypes playing a role in inflammatory flare-ups. We conclude by reflecting on our observations in the context of the ongoing pursuance of novel therapeutic interventions.

A mathematical model for wound healing in the reef-building coral Pocillopora damicornis

Coral reefs, among the most diverse ecosystems on Earth, currently face major threats from pollution, unsustainable fishing practices , and perturbations in environmental parameters brought on by climate change. Corals also sustain regular wounding from other sea life and human activity. Recent reef restoration practices have even involved intentional wounding by systematically breaking coral fragments and relocating them to revitalize damaged reefs, a practice known as microfragmentation. Despite its importance, very little research has explored the inner mechanisms of wound healing in corals. Some reef-building corals have been observed to initiate an immunological response to wounding similar to that observed in mammalian species. Utilizing prior models of wound healing in mammalian species as the mathematical basis, we formulated a mechanistic model of wound healing, including observations of the immune response and tissue repair in scleractinian corals for the species Pocillopora damicornis. The model consists of four differential equations which track changes in remaining wound debris, number of cells involved in inflammation, number of cells involved in proliferation, and amount of wound closure through re-epithelialization. The model is fit to experimental wound size data from linear and circular shaped wounds on a live coral fragment. Mathematical methods, including numerical simulations and local sensitivity analysis, were used to analyze the resulting model. The parameter space was also explored to investigate drivers of other possible wound outcomes. This model serves as a first step in generating mathematical models for wound healing in corals that will not only aid in the understanding of wound healing as a whole, but also help optimize reef restoration practices and predict recovery behavior after major wounding events.

Characterization of differences in immune responses during bolus and continuous infusion endotoxin challenges using mathematical modelling

Endotoxin administration is commonly used to study the inflammatory response, and though traditionally given as a bolus injection, it can be administered as a continuous infusion over multiple hours. Several studies hypothesize that the latter better represents the prolonged and pronounced inflammation observed in conditions like sepsis. Yet very few experimental studies have administered endotoxin using both strategies, leaving significant gaps in determining the underlying mechanisms responsible for their differing immune responses. We used mathematical modelling to analyse cytokine data from two studies administering a 2 ng kg-1 dose of endotoxin, one as a bolus and the other as a continuous infusion over 4 h. Using our model, we simulated the dynamics of mean and subject-specific cytokine responses as well as the response to long-term endotoxin administration. Cytokine measurements revealed that the bolus injection led to significantly higher peaks for interleukin (IL)-8, while IL-10 reaches higher peaks during continuous administration. Moreover, the peak timing of all measured cytokines occurred later with continuous infusion. We identified three model parameters that significantly differed between the two administration methods. Monocyte activation of IL-10 was greater during the continuous infusion, while tumour necrosis factor α $ {\alpha} $ and IL-8 recovery rates were faster for the bolus injection. This suggests that a continuous infusion elicits a stronger, longer-lasting systemic reaction through increased stimulation of monocyte anti-inflammatory mediator production and decreased recovery of pro-inflammatory catalysts. Furthermore, the continuous infusion model exhibited prolonged inflammation with recurrent peaks resolving within 2 days during long-term (20-32 h) endotoxin administration.

Agent-based vs. equation-based multi-scale modeling for macrophage polarization

Macrophages show high plasticity and result in heterogenic subpopulations or polarized states identified by specific cellular markers. These immune cells are typically characterized as pro-inflammatory, or classically activated M1, and anti-inflammatory, or alternatively activated M2. However, a more precise definition places them along a spectrum of activation where they may exhibit a number of pro- or anti-inflammatory roles. To understand M1-M2 dynamics in the context of a localized response and explore the results of different mathematical modeling approaches based on the same biology, we utilized two different modeling techniques, ordinary differential equation (ODE) modeling and agent-based modeling (ABM), to simulate the spectrum of macrophage activation to general pro- and anti-inflammatory stimuli on an individual and multi-cell level. The ODE model includes two hallmark pro- and anti-inflammatory signaling pathways and the ABM incorporates similar M1-M2 dynamics but in a spatio-temporal platform. Both models link molecular signaling with cellular-level dynamics. We then performed simulations with various initial conditions to replicate different experimental setups. Similar results were observed in both models after tuning to a common calibrating experiment. Comparing the two models' results sheds light on the important features of each modeling approach. When more data is available these features can be considered when choosing techniques to best fit the needs of the modeler and application.

Mathematical models of cystic fibrosis as a systemic disease

Cystic fibrosis (CF) is widely known as a disease of the lung, even though it is in truth a systemic disease, whose symptoms typically manifest in gastrointestinal dysfunction first. CF ultimately impairs not only the pancreas and intestine but also the lungs, gonads, liver, kidneys, bones, and the cardiovascular system. It is caused by one of several mutations in the gene of the epithelial ion channel protein CFTR. Intense research and improved antimicrobial treatments during the past eight decades have steadily increased the predicted life expectancy of a person with CF (pwCF) from a few weeks to over 50 years. Moreover, several drugs ameliorating the sequelae of the disease have become available in recent years, and notable treatments of the root cause of the disease have recently generated substantial improvements in health for some but not all pwCF. Yet, numerous fundamental questions remain unanswered. Complicating CF, for instance in the lung, is the fact that the associated insufficient chloride secretion typically perturbs the electrochemical balance across epithelia and, in the airways, leads to the accumulation of thick, viscous mucus and mucus plaques that cannot be cleared effectively and provide a rich breeding ground for a spectrum of bacterial and fungal communities. The subsequent infections often become chronic and respond poorly to antibiotic treatments, with outcomes sometimes only weakly correlated with the drug susceptibility of the target pathogen. Furthermore, in contrast to rapidly resolved acute infections with a single target pathogen, chronic infections commonly involve multi-species bacterial communities, called "infection microbiomes," that develop their own ecological and evolutionary dynamics. It is presently impossible to devise mathematical models of CF in its entirety, but it is feasible to design models for many of the distinct drivers of the disease. Building upon these growing yet isolated modeling efforts, we discuss in the following the feasibility of a multi-scale modeling framework, known as template-and-anchor modeling, that allows the gradual integration of refined sub-models with different granularity. The article first reviews the most important biomedical aspects of CF and subsequently describes mathematical modeling approaches that already exist or have the potential to deepen our understanding of the multitude aspects of the disease and their interrelationships. The conceptual ideas behind the approaches proposed here do not only pertain to CF but are translatable to other systemic diseases. This article is categorized under: Congenital Diseases > Computational Models.

A mathematical model of SARS-CoV-2 immunity predicts paxlovid rebound

Nirmatrelvir/ritonavir (Paxlovid), an oral antiviral medication targeting SARS-CoV-2, remains an important treatment for COVID-19. Initial studies of nirmatrelvir/ritonavir were performed in SARS-CoV-2 unvaccinated patients without prior confirmed SARS-CoV-2 infection; however, most individuals have now either been vaccinated and/or have experienced SARS-CoV-2 infection. After nirmatrelvir/ritonavir became widely available, reports surfaced of "Paxlovid rebound," a phenomenon in which symptoms (and SARS-CoV-2 test positivity) would initially resolve, but after finishing treatment, symptoms and test positivity would return. We used a previously described parsimonious mathematical model of immunity to SARS-CoV-2 infection to model the effect of nirmatrelvir/ritonavir treatment in unvaccinated and vaccinated patients. Model simulations show that viral rebound after treatment occurs only in vaccinated patients, while unvaccinated (SARS-COV-2 naïve) patients treated with nirmatrelvir/ritonavir do not experience any rebound in viral load. This work suggests that an approach combining parsimonious models of the immune system could be used to gain important insights in the context of emerging pathogens.

Modeling of adaptive immunity uncovers disease tolerance mechanisms

When an organism is challenged with a pathogen a cascade of events unfolds. The innate immune system rapidly mounts a preliminary nonspecific defense, while the acquired immune system slowly develops microbe-killing specialists. These responses cause inflammation, and along with the pathogen cause direct and indirect tissue damage, which anti-inflammatory mediators seek to temper. This interplay of systems is credited for maintaining homeostasis but may produce unexpected results such as disease tolerance. Tolerance is characterized by the persistence of pathogen and damage mitigation, where the relevant mechanisms are poorly understood. In this work we develop an ordinary differential equations model of the immune response to infection in order to identify key components in tolerance. Bifurcation analysis uncovers health, immune- and pathogen-mediated death clinical outcomes dependent on pathogen growth rate. We demonstrate that decreasing the inflammatory response to damage and increasing the strength of the immune system gives rise to a region in which limit cycles, or periodic solutions, are the only biological trajectories. We then describe areas of parameter space corresponding to disease tolerance by varying immune cell decay, pathogen removal, and lymphocyte proliferation rates.

Propagation of an idealized infection in an airway tree, consequences of the inflammation on the oxygen transfer to blood

A mathematical model of infection, inflammation and immune response in an idealized bronchial tree is developed. This work is based on a model from the literature that is extended to account for the propagation dynamics of an infection between the airways. The inflammation affects the size of the airways, the air flows distribution in the airway tree, and, consequently, the oxygen transfers to blood. We test different infections outcomes and propagation speed. In the hypotheses of our model, the inflammation can reduce notably and sometimes drastically the oxygen flow to blood. Our model predicts how the air flows and oxygen exchanges reorganize in the tree during an infection. Our results highlight the links between the localization of the infection and the amplitude of the loss of oxygen flow to blood. We show that a compensation phenomena due to the reorganization of the flow exists, but that it remains marginal unless the power produced the ventilation muscles is increased. Our model forms a first step towards a better understanding of the dynamics of bronchial infections.

Mathematical Model of Triple-Negative Breast Cancer in Response to Combination Chemotherapies

Triple-negative breast cancer (TNBC) is a heterogenous disease that is defined by its lack of targetable receptors, thus limiting treatment options and resulting in higher rates of metastasis and recurrence. Combination chemotherapy treatments, which inhibit tumor cell proliferation and regeneration, are a major component of standard-of-care treatment of TNBC. In this manuscript, we build a coupled ordinary differential equation model of TNBC with compartments that represent tumor proliferation, necrosis, apoptosis, and immune response to computationally describe the biological tumor affect to a combination of chemotherapies, doxorubicin (DRB) and paclitaxel (PTX). This model is parameterized using longitudinal [18F]-fluorothymidine positron emission tomography (FLT-PET) imaging data which allows for a noninvasive molecular imaging approach to quantify the tumor proliferation and tumor volume measurements for two murine models of TNBC. Animal models include a human cell line xenograft model, MDA-MB-231, and a syngeneic 4T1 mammary carcinoma model. The mathematical models are parameterized and the percent necrosis at the end time point is predicted and validated using histological hematoxylin and eosin (H&E) data. Global Sobol' sensitivity analysis is conducted to further understand the role each parameter plays in the model's goodness of fit to the data. In both the MDA-MB-231 and the 4T1 tumor models, the designed mathematical model can accurately describe both tumor volume changes and final necrosis volume. This can give insight into the ordering, dosing, and timing of DRB and PTX treatment. More importantly, this model can also give insight into future novel combinations of therapies and how the immune system plays a role in therapeutic response to TNBC, due to its calibration to two types of TNBC murine models.

MATHEMATICAL INSIGHTS INTO THE DYNAMICS OF INNATE IMMUNE RESPONSE DURING INFLAMMATION

Innate immune system cells activate in response to infection and trigger an acute inflammatory reaction to restore tissue homeostasis and promote subsequent tissue repair. Their activation and functions must be very well regulated to avoid tissue damage, organ dysfunction, or even death. In this work, a new set of mathematical models is presented to examine the dynamics of the innate immune system response to tissue damage and provide further understanding of the role of the innate immune system during the early stages of an inflammatory response. Different damaged cells production functions are proposed to represent the effect of secondary tissue damage by the innate immune system. The stability and bifurcation analyses of the model reveal that there is an important threshold parameter that can be controlled in order to avoid sustained chronic inflammation and secure a successful healing outcome. A set of numerical simulations is also performed to support the presented theoretical results and demonstrate the medical applicability of the new mathematical model.

Chemical compositions and pharmacological activities of natural musk (Moschus) and artificial musk: A review

ETHNOPHARMACOLOGICAL RELEVANCE

Natural musk (Moschus), derived naturally from male musk deer (Moschus berezovskii Flerov, Moschus sifanicus Przewalski, or Moschus moschiferus Linnaeus), has long been an important component of traditional Chinese medicine (TCM), and was used as resuscitation, blood circulation, and collateral drainage. detumescence and pain relief. Artificial musk was researched and applied into TCM as natural musk being as unsustainable resources.

AIM OF THE STUDY

We mainly summarized chemical compositions, pharmacological activities and mechanism of action of natural and artificial musk, and designed to serve as a foundation for further research into musk chemical compositions and pharmacological effect.

MATERIALS AND METHODS

Those mainstream scientific databases including Google Scholar, ScienceDirect, SpringerLink, CNKI, Wiley Online Library, web of science, were used for searching with below "Keywords", as well as literature-tracking. Literatures spanned 1962 to 2021, and involved into Chinese, English, Janpanese, Korean.

RESULTS

Natural musk contains some very desirable but scarce compounds, as well as their biological features, which led to the development of artificial musk. The chemical ingredients, pharmacological activities, and mechanisms of action of natural and artificial musk are summarized and compared in this paper. Polypeptide and protein, muscone, musclide, steroids, muscopyridine, and other chemical constituents of musk demonstrated important therapeutic properties against inflammation, immune system disorders, neurological disorders, cardiovascular system disorders, and so on. The mechanism of action contributed to effect on mediators, acceptors and relative signal pathways.

CONCLUSIONS

Natural and artificial musk were revealed having some activated compounds, and showed excellent pharmacological effect. Meantime, above two sides of natural and artificial musk ought to get further research.

A data-driven model of the role of energy in sepsis

Exposure to pathogens elicits a complex immune response involving multiple interdependent pathways. This response may mitigate detrimental effects and restore health but, if imbalanced, can lead to negative outcomes including sepsis. This complexity and need for balance pose a challenge for clinicians and have attracted attention from modelers seeking to apply computational tools to guide therapeutic approaches. In this work, we address a shortcoming of such past efforts by incorporating the dynamics of energy production and consumption into a computational model of the acute immune response. With this addition, we performed fits of model dynamics to data obtained from non-human primates exposed to Escherichia coli. Our analysis identifies parameters that may be crucial in determining survival outcomes and also highlights energy-related factors that modulate the immune response across baseline and altered glucose conditions.

Divergent COVID-19 Disease Trajectories Predicted by a DAMP-Centered Immune Network Model

COVID-19 presentations range from mild to moderate through severe disease but also manifest with persistent illness or viral recrudescence. We hypothesized that the spectrum of COVID-19 disease manifestations was a consequence of SARS-CoV-2-mediated delay in the pathogen-associated molecular pattern (PAMP) response, including dampened type I interferon signaling, thereby shifting the balance of the immune response to be dominated by damage-associated molecular pattern (DAMP) signaling. To test the hypothesis, we constructed a parsimonious mechanistic mathematical model. After calibration of the model for initial viral load and then by varying a few key parameters, we show that the core model generates four distinct viral load, immune response and associated disease trajectories termed “patient archetypes”, whose temporal dynamics are reflected in clinical data from hospitalized COVID-19 patients. The model also accounts for responses to corticosteroid therapy and predicts that vaccine-induced neutralizing antibodies and cellular memory will be protective, including from severe COVID-19 disease. This generalizable modeling framework could be used to analyze protective and pathogenic immune responses to diverse viral infections.

An Agent-Based Model of Combination Oncolytic Viral Therapy and Anti-PD-1 Immunotherapy Reveals the Importance of Spatial Location When Treating Glioblastoma

Simple Summary

A combination of oncolytic viral therapy and immunotherapy provides an alternative option to the standard of care for treating the lethal brain tumor glioblastoma (GBM). Although this combination therapy shows promise, there are many unknown questions regarding how the tumor landscape and spatial dosing strategies impact the effectiveness of the treatment. Our study aims to shed light on these questions using a novel spatially explicit computational model of GBM response to treatment. Our results suggest that oncolytic viral dosing in the location of highest tumor cell density leads to substantial tumor size reduction over viral dosing in the center of the tumor. These results can help to inform future clinical trials and more effective treatment strategies for oncolytic viral therapy in GBM patients.

Abstract

Oncolytic viral therapies and immunotherapies are of growing clinical interest due to their selectivity for tumor cells over healthy cells and their immunostimulatory properties. These treatment modalities provide promising alternatives to the standard of care, particularly for cancers with poor prognoses, such as the lethal brain tumor glioblastoma (GBM). However, uncertainty remains regarding optimal dosing strategies, including how the spatial location of viral doses impacts therapeutic efficacy and tumor landscape characteristics that are most conducive to producing an effective immune response. We develop a three-dimensional agent-based model (ABM) of GBM undergoing treatment with a combination of an oncolytic Herpes Simplex Virus and an anti-PD-1 immunotherapy. We use a mechanistic approach to model the interactions between distinct populations of immune cells, incorporating both innate and adaptive immune responses to oncolytic viral therapy and including a mechanism of adaptive immune suppression via the PD-1/PD-L1 checkpoint pathway. We utilize the spatially explicit nature of the ABM to determine optimal viral dosing in both the temporal and spatial contexts. After proposing an adaptive viral dosing strategy that chooses to dose sites at the location of highest tumor cell density, we find that, in most cases, this adaptive strategy produces a more effective treatment outcome than repeatedly dosing in the center of the tumor.

Towards in silico Models of the Inflammatory Response in Bone Fracture Healing

In silico modeling is a powerful strategy to investigate the biological events occurring at tissue, cellular and subcellular level during bone fracture healing. However, most current models do not consider the impact of the inflammatory response on the later stages of bone repair. Indeed, as initiator of the healing process, this early phase can alter the regenerative outcome: if the inflammatory response is too strongly down- or upregulated, the fracture can result in a non-union. This review covers the fundamental information on fracture healing, in silico modeling and experimental validation. It starts with a description of the biology of fracture healing, paying particular attention to the inflammatory phase and its cellular and subcellular components. We then discuss the current state-of-the-art regarding in silico models of the immune response in different tissues as well as the bone regeneration process at the later stages of fracture healing. Combining the aforementioned biological and computational state-of-the-art, continuous, discrete and hybrid modeling technologies are discussed in light of their suitability to capture adequately the multiscale course of the inflammatory phase and its overall role in the healing outcome. Both in the establishment of models as in their validation step, experimental data is required. Hence, this review provides an overview of the different in vitro and in vivo set-ups that can be used to quantify cell- and tissue-scale properties and provide necessary input for model credibility assessment. In conclusion, this review aims to provide hands-on guidance for scientists interested in building in silico models as an additional tool to investigate the critical role of the inflammatory phase in bone regeneration.

Modeling Neuroimmune Interactions in Human Subjects and Animal Models to Predict Subtype-Specific Multidrug Treatments for Gulf War Illness

Gulf War Illness (GWI) is a persistent chronic neuroinflammatory illness exacerbated by external stressors and characterized by fatigue, musculoskeletal pain, cognitive, and neurological problems linked to underlying immunological dysfunction for which there is no known treatment. As the immune system and the brain communicate through several signaling pathways, including the hypothalamic–pituitary–adrenal (HPA) axis, it underlies many of the behavioral and physiological responses to stressors via blood-borne mediators, such as cytokines, chemokines, and hormones. Signaling by these molecules is mediated by the semipermeable blood–brain barrier (BBB) made up of a monocellular layer forming an integral part of the neuroimmune axis. BBB permeability can be altered and even diminished by both external factors (e.g., chemical agents) and internal conditions (e.g., acute or chronic stress, or cross-signaling from the hypothalamic–pituitary–gonadal (HPG) axis). Such a complex network of regulatory interactions that possess feed-forward and feedback connections can have multiple response dynamics that may include several stable homeostatic states beyond normal health. Here we compare immune and hormone measures in the blood of human clinical samples and mouse models of Gulf War Illness (GWI) subtyped by exposure to traumatic stress for subtyping this complex illness. We do this via constructing a detailed logic model of HPA–HPG–Immune regulatory behavior that also considers signaling pathways across the BBB to neuronal–glial interactions within the brain. We apply conditional interactions to model the effects of changes in BBB permeability. Several stable states are identified in the system beyond typical health. Following alignment of the human and mouse blood profiles in the context of the model, mouse brain sample measures were used to infer the neuroinflammatory state in human GWI and perform treatment simulations using a genetic algorithm to optimize the Monte Carlo simulations of the putative treatment strategies aimed at returning the ill system back to health. We identify several ideal multi-intervention strategies and potential drug candidates that may be used to treat chronic neuroinflammation in GWI.

Dynamically linking influenza virus infection kinetics, lung injury, inflammation, and disease severity

Influenza viruses cause a significant amount of morbidity and mortality. Understanding host immune control efficacy and how different factors influence lung injury and disease severity are critical. We established and validated dynamical connections between viral loads, infected cells, CD8⁺ T cells, lung injury, inflammation, and disease severity using an integrative mathematical model-experiment exchange. Our results showed that the dynamics of inflammation and virus-inflicted lung injury are distinct and nonlinearly related to disease severity, and that these two pathologic measurements can be independently predicted using the model-derived infected cell dynamics. Our findings further indicated that the relative CD8⁺ T cell dynamics paralleled the percent of the lung that had resolved with the rate of CD8⁺ T cell-mediated clearance rapidly accelerating by over 48,000 times in 2 days. This complimented our analyses showing a negative correlation between the efficacy of innate and adaptive immune-mediated infected cell clearance, and that infection duration was driven by CD8⁺ T cell magnitude rather than efficacy and could be significantly prolonged if the ratio of CD8⁺ T cells to infected cells was sufficiently low. These links between important pathogen kinetics and host pathology enhance our ability to forecast disease progression, potential complications, and therapeutic efficacy.

A Turing mechanism in order to explain the patchy nature of Crohn's disease

Crohn's disease is an inflammatory bowel disease (IBD) that is not well understood. In particular, unlike other IBDs, the inflamed parts of the intestine compromise deep layers of the tissue and are not continuous but separated and distributed through the whole gastrointestinal tract, displaying a patchy inflammatory pattern. In the present paper, we introduce a toy-model which might explain the appearance of such patterns. We consider a reaction-diffusion system involving bacteria and phagocyte and prove that, under certain conditions, this system might reproduce an activator-inhibitor dynamic leading to the occurrence of Turing-type instabilities. In other words, we prove the existence of stable stationary solutions that are spatially periodic and do not vanish in time. We also propose a set of parameters for which the system exhibits such phenomena and compare it with realistic parameters found in the literature. This is the first time, as far as we know, that a Turing pattern is investigated in inflammatory models.

Compartmental Model Suggests Importance of Innate Immune Response to COVID-19 Infection in Rhesus Macaques

The pandemic outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has quickly spread worldwide, creating a serious health crisis. The virus is primarily associated with flu-like symptoms but can also lead to severe pathologies and death. We here present an ordinary differential equation model of the intrahost immune response to SARS-CoV-2 infection, fitted to experimental data gleaned from rhesus macaques. The model is calibrated to data from a nonlethal infection, but the model can replicate behavior from various lethal scenarios as well. We evaluate the sensitivity of the model to biologically relevant parameters governing the strength and efficacy of the immune response. We also simulate the effect of both anti-inflammatory and antiviral drugs on the host immune response and demonstrate the ability of the model to lessen the severity of a formerly lethal infection with the addition of the appropriately calibrated drug. Our model emphasizes the importance of tight control of the innate immune response for host survival and viral clearance.

Dynamical Simulation of Effective Stem Cell Transplantation for Modulation of Microglia Responses in Stroke Treatment

Stem cell transplantation therapy may inhibit inflammation during stroke and increase the presence of healthy cells in the brain. The novelty of this work, is to introduce a new mathematical model of stem cells transplanted to treat stroke. This manuscript studies the stability of the mathematical model by using the current biological information on stem cell therapy as a possible treatment for inflammation from microglia during stroke. The model is proposed to represent the dynamics of various immune brain cells (resting microglia, pro-inflammation microglia, and anti-inflammation microglia), brain tissue damage and stem cells transplanted. This model is based on a set of five ordinary differential equations and explores the beneficial effects of stem cells transplanted at early stages of inflammation during stroke. The Runge–Kutta method is used to discuss the model analytically and solve it numerically. The results of our simulations are qualitatively consistent with those observed in experiments in vivo, suggesting that the transplanted stem cells could contribute to the increase in the rate of ant-inflammatory microglia and decrease the damage from pro-inflammatory microglia. It is found from the analysis and simulation results that stem cell transplantation can help stroke patients by modulation of the immune response during a stroke and decrease the damage on the brain. In conclusion, this approach may increase the contributions of stem cells transplanted during inflammation therapy in stroke and help to study various therapeutic strategies for stem cells to reduce stroke damage at the early stages.

Bistable Mathematical Model of Neutrophil Migratory Patterns After LPS-Induced Epigenetic Reprogramming

The highly controlled migration of neutrophils toward the site of an infection can be altered when they are trained with lipopolysaccharides (LPS), with high dose LPS enhancing neutrophil migratory pattern toward the bacterial derived source signal and super-low dose LPS inducing either migration toward an intermediary signal or dysregulation and oscillatory movement. Empirical studies that use microfluidic chemotaxis-chip devices with two opposing chemoattractants showed differential neutrophil migration after challenge with different LPS doses. The epigenetic alterations responsible for changes in neutrophil migratory behavior are unknown. We developed two mathematical models that evaluate the mechanistic interactions responsible for neutrophil migratory decision-making when exposed to competing chemoattractants and challenged with LPS. The first model, which considers the interactions between the receptor densities of two competing chemoattractants, their kinases, and LPS, displayed bistability between high and low ratios of primary to intermediary chemoattractant receptor densities. In particular, at equilibrium, we observe equal receptor densities for low LPS (< 15ng/mL); and dominance of receptors for the primary chemoattractant for high LPS (> 15ng/mL). The second model, which included additional interactions with an extracellular signal-regulated kinase in both phosphorylated and non-phosphorylated forms, has an additional dynamic outcome, oscillatory dynamics for both receptors, as seen in the data. In particular, it found equal receptor densities in the absence of oscillation for super-low and high LPS challenge (< 0.4 and 1.1 376 ng/mL). Predicting the mechanisms and the type of external LPS challenge responsible for neutrophils migration toward pro-inflammatory chemoattractants, migration toward pro-tolerant chemoattractants, or oscillatory movement is necessary knowledge in designing interventions against immune diseases, such as sepsis.

Aging-induced fragility of the immune system

The adaptive and innate branches of the vertebrate immune system work in close collaboration to protect organisms from harmful pathogens. As an organism ages its immune system undergoes immunosenescence, characterized by declined performance or malfunction in either immune branch, which can lead to disease and death. In this study we develop a mathematical framework of coupled innate and adaptive immune responses, namely the integrated immune branch (IIB) model. This model describes dynamics of immune components in both branches, uses a shape-space representation to encode pathogen-specific immune memory, and exhibits three steady states - health, septic death, and chronic inflammation - qualitatively similar to clinically-observed immune outcomes. In this model, the immune system (initialized in the health state) is subjected to a sequence of pathogen encounters, and we use the number of prior pathogen encounters as a proxy for the "age" of the immune system. We find that repeated pathogen encounters may trigger a fragility in which any encounter with a novel pathogen will cause the system to irreversibly switch from health to chronic inflammation. This transition is consistent with the onset of "inflammaging", a condition observed in aged individuals who experience chronic low-grade inflammation even in the absence of pathogens. The IIB model predicts that the onset of chronic inflammation strongly depends on the history of encountered pathogens; the timing of onset differs drastically when the same set of infections occurs in a different order. Lastly, the coupling between the innate and adaptive immune branches generates a trade-off between rapid pathogen clearance and a delayed onset of immunosenescence. Overall, by considering the complex feedback between immune compartments, our work suggests potential mechanisms for immunosenescence and provides a theoretical framework at the system level and on the scale of an organism's lifetime to account for clinical observations.

Dynamic analysis of the role of innate immunity in SEIS epidemic model

Consideration of every important aspect while modeling a disease makes the model more precise and the disease eradication strategy more powerful. In the present paper, we analyze the importance of innate immunity on SEIS modeling. We propose an SEIS model with Holling type II and type III functions representing innate immunity. We find the existence and stability conditions for the equilibria. When innate immunity is in the form of Holling type II function, the disease-free equilibrium exists for reproduction number less than unity and is locally asymptotically stable, and supercritical transcritical (forward) as well as subcritical transcritical (backward) bifurcation may occur where the contact rate β = β ∗ acts as the bifurcation parameter. Hence, disease-free equilibrium need not be globally stable. For reproduction number greater than unity unique endemic equilibrium exists which is locally asymptotically stable. The global stability conditions for the same are deduced with the help of Lozinski i ˘ measure. When innate immunity is considered a Holling type III function, the disease-free equilibrium point exists for reproduction number less than unity and is locally as well as globally stable. The existence of either unique or multiple endemic equilibria is found when reproduction number is greater than unity, and there exists at least one locally asymptotically stable equilibrium point and bistability can also be encountered. The conditions for the existence of Andronov-Hopf bifurcation are deduced for both cases. Moreover, we observe that ignoring innate immunity annihilates the possibility of Andronov-Hopf bifurcation. Numerical simulation is performed to validate the mathematical findings. Comparing the obtained results to the case when innate immunity is ignored, it is deduced that ignoring it ends the possibility of backward bifurcation, Andronov-Hopf bifurcation as well as the existence of multiple equilibria, and it also leads to the prediction of higher infection than the actual which may deflect the accuracy of the model to a high extent. This would further lead to false predictions and inefficient disease control strategies which in turn would make disease eradication a difficult and more expensive task.

Modeling and simulations of CoViD-19 molecular mechanism induced by cytokines storm during SARS-CoV2 infection

It is highly desired to explore the interventions of COVID-19 for early treatment strategies. Such interventions are still under consideration. A model is benchmarked research and comprises target cells, virus infected cells, immune cells, pro-inflammatory cytokines, and, anti-inflammatory cytokine. The interaction of the drug with the inflammatory sub-system is analyzed with the aid of kinetic modeling. The impact of drug therapy on the immune cells is modelled and the computational framework is verified with the aid of numerical simulations. The work includes a significant hypothesis that quantifies the complex dynamics of the infection, by relating it to the effect of the inflammatory syndrome generated by IL-6. In this paper we use the cancer immunoediting process: a dynamic process initiated by cancer cells in response to immune surveillance of the immune system that it can be conceptualized by an alternating movement that balances immune protection with immune evasion. The mechanisms of resistance to immunotherapy seem to broadly overlap with those used by cancers as they undergo immunoediting to evade detection by the immune system. In this process the immune system can both constrain and promote tumour development, which proceeds through three phases termed: (i) Elimination, (ii) Equilibrium, and, (iii) Escape [1]. We can also apply these concepts to viral infection, which, although it is not exactly “immunoediting”, has many points in common and helps to understand how it expands into an “untreated” host and can help in understanding the SARS-CoV2 virus infection and treatment model.

Spatial considerations in the resolution of inflammation: Elucidating leukocyte interactions via an experimentally-calibrated agent-based model

Many common medical conditions (such as cancer, arthritis, chronic obstructive pulmonary disease (COPD), and others) are associated with inflammation, and even more so when combined with the effects of ageing and multimorbidity. While the inflammatory response varies in different tissue types, under disease and in response to therapeutic interventions, it has common interactions that occur between immune cells and inflammatory mediators. Understanding these underlying inflammatory mechanisms is key in progressing treatments and therapies for numerous inflammatory conditions. It is now considered that constituent mechanisms of the inflammatory response can be actively manipulated in order to drive resolution of inflammatory damage; particularly, those mechanisms related to the pro-inflammatory role of neutrophils and the anti-inflammatory role of macrophages. In this article, we describe the assembly of a hybrid mathematical model in which the spatial spread of inflammatory mediators is described through partial differential equations, and immune cells (neutrophils and macrophages) are described individually via an agent-based modelling approach. We pay close attention to how immune cells chemotax toward pro-inflammatory mediators, presenting a model for cell chemotaxis that is calibrated against experimentally observed cell trajectories in healthy and COPD-affected scenarios. We illustrate how variations in key model parameters can drive the switch from resolution of inflammation to chronic outcomes, and show that aberrant neutrophil chemotaxis can move an otherwise healthy outcome to one of chronicity. Finally, we reflect on our results in the context of the on-going hunt for new therapeutic interventions.

Inflammation is the body’s primary defence to harmful stimuli such as infections, toxins and tissue strain but also underlies a much broader range of conditions, including asthma, arthritis and cancer. The inflammatory response is key in resolving injury to facilitate recovery, and involves a range of interactions between immune cells (leukocytes, neutrophils and macrophages in particular) and inflammatory mediators. Immune cells are recruited from the blood stream in response to injury. Once in tissue, neutrophils release toxins to kill invading agents and resolve damage; however, if not carefully managed by other immune cells (mainly macrophages), their responses can increase inflammation instead of helping to resolve it. We model these interactions in response to damage using a spatial model, examining how a healthy response can prevent localised inflammation from spreading. We pay close attention to how cells migrate toward the damaged area, as many inflammatory conditions are associated with impairment of this process. We calibrate our model against experimentally-observed cell trajectories from healthy patients and patients with chronic obstructive pulmonary disease. We illustrate that a healthy outcome depends strongly upon efficient cell migration and a delicate balance between the pro- and anti-inflammatory effects of neutrophils and macrophages.

A mathematical model of the multiple sclerosis plaque

Multiple sclerosis is an autoimmune disease that affects white matter in the central nervous system. It is one of the primary causes of neurological disability among young people. Its characteristic pathological lesion is called a plaque, a zone of inflammatory activity and tissue destruction that expands radially outward by destroying the myelin and oligodendrocytes of white matter. The present paper develops a mathematical model of the multiple sclerosis plaques. Although these plaques do not provide reliable information of the clinical disability in MS, they are nevertheless useful as a primary outcome measure of Phase II trials. The model consists of a system of partial differential equations in a simplified geometry of the lesion, consisting of three domains: perivascular space, demyelinated plaque, and white matter. The model describes the activity of various pro- and anti-inflammatory cells and cytokines in the plaque, and quantifies their effect on plaque growth. We show that volume growth of plaques are in qualitative agreement with reported clinical studies of several currently used drugs. We then use the model to explore treatments with combinations of such drugs, and with experimental drugs. We finally consider the benefits of early vs. delayed treatment.

Use of Angiotensin-Converting Enzyme Inhibitors and Angiotensin II Receptor Blockers During the COVID-19 Pandemic: A Modeling Analysis

Angiotensin-converting enzyme inhibitors (ACEi) and angiotensin II receptor blockers (ARB) are frequently prescribed for a range of diseases including hypertension, proteinuric chronic kidney disease, and heart failure. There is evidence indicating that these drugs upregulate ACE2, a key component of the renin-angiotensin system (RAS) and is found on the cells of a number of tissues, including the epithelial cells in the lungs. While ACE2 has a beneficial role in many diseases such as hypertension, diabetes, and cardiovascular disease, it also serves as a receptor for both SARS-CoV and SARS-CoV-2 via binding with the spike protein of the virus, thereby allowing it entry into host cells. Thus, it has been suggested that these therapies can theoretically increase the risk of SARS- CoV-2 infection and cause more severe COVID-19. Given the success of ACEi and ARBs in cardiovascular diseases, we seek to gain insights into the implications of these medications in the pathogenesis of COVID-19. To that end, we have developed a mathematical model that represents the RAS, binding of ACE2 with SARS-CoV-2 and the subsequent cell entry, and the host's acute inflammatory response. The model can simulate different levels of SARS-CoV-2 exposure, and represent the effect of commonly prescribed anti-hypertensive medications, ACEi and ARB, and predict tissue damage. Model simulations indicate that whether the extent of tissue damage may be exacerbated by ACEi or ARB treatment depends on a number of factors, including the level of existing inflammation, dosage, and the effect of the drugs on ACE2 protein abundance. The findings of this study can serve as the first step in the development of appropriate and more comprehensive guidelines for the prescription of ACEi and ARB in the current and future coronavirus pandemics.

Spectral clustering of risk score trajectories stratifies sepsis patients by clinical outcome and interventions received

Sepsis is not a monolithic disease, but a loose collection of symptoms with diverse outcomes. Thus, stratification and subtyping of sepsis patients is of great importance. We examine the temporal evolution of patient state using our previously-published method for computing risk of transition from sepsis into septic shock. Risk trajectories diverge into four clusters following early prediction of septic shock, stratifying by outcome: the highest-risk and lowest-risk groups have a 76.5% and 10.4% prevalence of septic shock, and 43% and 18% mortality, respectively. These clusters differ also in treatments received and median time to shock onset. Analyses reveal the existence of a rapid (30–60 min) transition in risk at the time of threshold crossing. We hypothesize that this transition occurs as a result of the failure of compensatory biological systems to cope with infection, resulting in a bifurcation of low to high risk. Such a collapse, we believe, represents the true onset of septic shock. Thus, this rapid elevation in risk represents a potential new data-driven definition of septic shock.

Mathematical modelling of trastuzumab-induced immune response in an in vivo murine model of HER2+ breast cancer

The goal of this study is to develop an integrated, mathematical-experimental approach for understanding the interactions between the immune system and the effects of trastuzumab on breast cancer that overexpresses the human epidermal growth factor receptor 2 (HER2+). A system of coupled, ordinary differential equations was constructed to describe the temporal changes in tumour growth, along with intratumoural changes in the immune response, vascularity, necrosis and hypoxia. The mathematical model is calibrated with serially acquired experimental data of tumour volume, vascularity, necrosis and hypoxia obtained from either imaging or histology from a murine model of HER2+ breast cancer. Sensitivity analysis shows that model components are sensitive for 12 of 13 parameters, but accounting for uncertainty in the parameter values, model simulations still agree with the experimental data. Given theinitial conditions, the mathematical model predicts an increase in the immune infiltrates over time in the treated animals. Immunofluorescent staining results are presented that validate this prediction by showing an increased co-staining of CD11c and F4/80 (proteins expressed by dendritic cells and/or macrophages) in the total tissue for the treated tumours compared to the controls ($p < 0.03$). We posit that the proposed mathematical-experimental approach can be used to elucidate driving interactions between the trastuzumab-induced responses in the tumour and the immune system that drive the stabilization of vasculature while simultaneously decreasing tumour growth-conclusions revealed by the mathematical model that were not deducible from the experimental data alone.

The ups and downs of S. aureus nasal carriage

Staphylococcus aureus infections are a growing concern worldwide due to the increasing number of strains that exhibit antibiotic resistance. Recent studies have indicated that some percentage of people carry the bacteria in the nasal cavity and therefore are at a higher risk of subsequent, and more serious, infections in other parts of the body. However, individuals carrying the infection can be classified as only intermittent carriers versus persistent carriers, being able to eliminate the bacteria and later colonized again. Using a model of bacterial colonization of the anterior nares, we investigate oscillatory patterns related to intermittent carriage of S. aureus. Following several studies using global sensitivity analysis techniques, various insights into the model's behaviour were made including interacting effects of the bacteria's growth rate and movement in the mucus, suggesting parameter connections associated with biofilm-like behaviour. Here the bacterial growth rate and bacterial movement are explicitly connected, leading to expanded oscillatory behaviour in the model. We suggest possible implications that this oscillatory behaviour can have on the definition of intermittent carriage and discuss differences in the bacterial virulence dependent upon individual host health. Furthermore, we show that connecting the bacterial growth and movement also expands the region of the parameter space for which the bacteria are able to survive and persist.

Dynamical modeling of pro- and anti-inflammatory cytokines in the early stage of septic shock

A dynamical model of the pathophysiological behaviors of IL18 and IL10 cytokines with their receptors is tested against data for the case of early sepsis. The proposed approach considers the surroundings (organs and bone marrow) and the different subsystems (cells and cyctokines). The interactions between blood cells, cytokines and the surroundings are described via mass balances. Cytokines are adsorbed onto associated receptors at the cell surface. The adsorption is described by the Langmuir model and gives rise to the production of more cytokines and associated receptors inside the cell. The quantities of pro and anti-inflammatory cytokines present in the body are combined to give global information via an inflammation level function which describes the patient’s state. Data for parameter estimation comes from the Sepsis 48 H database. Comparisons between patient data and simulations are presented and are in good agreement. For the IL18/IL10 cytokine pair, 5 key parameters have been found. They are linked to pro-inflammatory IL18 cytokine and show that the early sepsis is driven by components of inflammatory character.

Mechanisms and Points of Control in the Spread of Inflammation: A Mathematical Investigation

Understanding the mechanisms that control the body’s response to inflammation is of key importance, due to its involvement in myriad medical conditions, including cancer, arthritis, Alzheimer’s disease and asthma. While resolving inflammation has historically been considered a passive process, since the turn of the century the hunt for novel therapeutic interventions has begun to focus upon active manipulation of constituent mechanisms, particularly involving the roles of apoptosing neutrophils, phagocytosing macrophages and anti-inflammatory mediators. Moreover, there is growing interest in how inflammatory damage can spread spatially due to the motility of inflammatory mediators and immune cells. For example, impaired neutrophil chemotaxis is implicated in causing chronic inflammation under trauma and in ageing, while neutrophil migration is an attractive therapeutic target in ailments such as chronic obstructive pulmonary disease. We extend an existing homogeneous model that captures interactions between inflammatory mediators, neutrophils and macrophages to incorporate spatial behaviour. Through bifurcation analysis and numerical simulation, we show that spatially inhomogeneous outcomes can present close to the switch from bistability to guaranteed resolution in the corresponding homogeneous model. Finally, we show how aberrant spatial mechanisms can play a role in the failure of inflammation to resolve and discuss our results within the broader context of seeking novel inflammatory treatments.

Complex Stability and an Irrevertible Transition Reverted by Peptide and Fibroblasts in a Dynamic Model of Innate Immunity

We here apply a control analysis and various types of stability analysis to an in silico model of innate immunity that addresses the management of inflammation by a therapeutic peptide. Motivation is the observation, both in silico and in experiments, that this therapy is not robust. Our modeling results demonstrate how (1) the biological phenomena of acute and chronic modes of inflammation may reflect an inherently complex bistability with an irrevertible flip between the two modes, (2) the chronic mode of the model has stable, sometimes unique, steady states, while its acute-mode steady states are stable but not unique, (3) as witnessed by TNF levels, acute inflammation is controlled by multiple processes, whereas its chronic-mode inflammation is only controlled by TNF synthesis and washout, (4) only when the antigen load is close to the acute mode's flipping point, many processes impact very strongly on cells and cytokines, (5) there is no antigen exposure level below which reduction of the antigen load alone initiates a flip back to the acute mode, and (6) adding healthy fibroblasts makes the transition from acute to chronic inflammation revertible, although (7) there is a window of antigen load where such a therapy cannot be effective. This suggests that triple therapies may be essential to overcome chronic inflammation. These may comprise (1) anti-immunoglobulin light chain peptides, (2) a temporarily reduced antigen load, and (3a) fibroblast repopulation or (3b) stem cell strategies.

Modeling Oncolytic Viral Therapy, Immune Checkpoint Inhibition, and the Complex Dynamics of Innate and Adaptive Immunity in Glioblastoma Treatment

Oncolytic viruses are of growing interest to cancer researchers and clinicians, due to their selectivity for tumor cells over healthy cells and their immunostimulatory properties. The immune response to an oncolytic virus plays a critical role in treatment efficacy. However, uncertainty remains regarding the circumstances under which the immune system either assists in eliminating tumor cells or inhibits treatment via rapid viral clearance, leading to the cessation of the immune response. In this work, we develop an ordinary differential equation model of treatment for a lethal brain tumor, glioblastoma, using an oncolytic Herpes Simplex Virus. We use a mechanistic approach to model the interactions between distinct populations of immune cells, incorporating both innate and adaptive immune responses to oncolytic viral therapy (OVT), and including a mechanism of adaptive immune suppression via the PD-1/PD-L1 checkpoint pathway. We focus on the tradeoff between viral clearance by innate immune cells and the innate immune cell-mediated recruitment of antiviral and antitumor adaptive immune cells. Our model suggests that when a tumor is treated with OVT alone, the innate immune cells' ability to clear the virus quickly after administration has a much larger impact on the treatment outcome than the adaptive immune cells' antitumor activity. Even in a highly antigenic tumor with a strong innate immune response, the faster recruitment of antitumor adaptive immune cells is not sufficient to offset the rapid viral clearance. This motivates our subsequent incorporation of an immunotherapy that inhibits the PD-1/PD-L1 checkpoint pathway by blocking PD-1, which we combine with OVT within the model. The combination therapy is most effective for a highly antigenic tumor or for intermediate levels of innate immune localization. Extreme levels of innate immune cell activity either clear the virus too quickly or fail to activate a sufficiently strong adaptive response, yielding ineffective combination therapy of GBM. Hence, we show that the innate and adaptive immune interactions significantly influence treatment response and that combining OVT with an immune checkpoint inhibitor expands the range of immune conditions that allow for tumor size reduction or clearance.

Dynamic Modelling of Interactions between Microglia and Endogenous Neural Stem Cells in the Brain during a Stroke

In this paper, we study the interactions between microglia and neural stem cells and the impact of these interactions on the brain cells during a stroke. Microglia cells, neural stem cells, the damage on brain cells from the stroke and the impacts these interactions have on living brain cells are considered in the design of mathematical models. The models consist of ordinary differential equations describing the effects of microglia on brain cells and the interactions between microglia and neural stem cells in the case of a stroke. Variables considered include: resident microglia, classically activated microglia, alternatively activated microglia, neural stem cells, tissue damage on cells in the brain, and the impacts these interactions have on living brain cells. The first model describes what happens in the brain at the stroke onset during the first three days without the generation of any neural stem cells. The second model studies the dynamic effect of microglia and neural stem cells on the brain cells following the generation of neural stem cells and potential recovery after this stage. We look at the stability and the instability of the models which are both studied analytically. The results show that the immune cells can help the brain by cleaning dead cells and stimulating the generation of neural stem cells; however, excessive activation may cause damage and affect the injured region. Microglia have beneficial and harmful functions after ischemic stroke. The microglia stimulate neural stem cells to generate new cells that substitute dead cells during the recovery stage but sometimes the endogenous neural stem cells are highly sensitive to inflammatory in the brain.

A Whole-Body Mathematical Model of Sepsis Progression and Treatment Designed in the BioGears Physiology Engine

Sepsis is a debilitating condition associated with a high mortality rate that greatly strains hospital resources. Though advances have been made in improving sepsis diagnosis and treatment, our understanding of the disease is far from complete. Mathematical modeling of sepsis has the potential to explore underlying biological mechanisms and patient phenotypes that contribute to variability in septic patient outcomes. We developed a comprehensive, whole-body mathematical model of sepsis pathophysiology using the BioGears Engine, a robust open-source virtual human modeling project. We describe the development of a sepsis model and the physiologic response within the BioGears framework. We then define and simulate scenarios that compare sepsis treatment regimens. As such, we demonstrate the utility of this model as a tool to augment sepsis research and as a training platform to educate medical staff.

Responses of pro- and anti-inflammatory cytokines in zebrafish liver exposed to sublethal doses of Aphanizomenon flosaquae DC-1 aphantoxins

Blooms of the dominant cyanobacterium Aphanizomenon flosaquae are frequently encountered in natural waters, and their secretion of neurotoxic paralytic shellfish toxins called aphantoxins threatens environmental safety and human health worldwide. The liver is the primary detoxification organ in animals, and its pro- and anti-inflammatory responses are important functions in the detoxification of toxins. Therefore, we investigated the response of these inflammatory factors to aphantoxins in the liver of zebrafish (Danio rerio). A. flosaquae DC-1 was sampled during blooms in Dianchi Lake, China and cultured, and the toxin was extracted and analyzed using high performance liquid chromatography. The primary constituents were gonyautoxins 1 (34.04%) and 5 (21.28%) and neosaxitoxin (12.77%). Zebrafish were injected intraperitoneally with 5.3 μg (low dose) or 7.61 μg (high dose) of saxitoxin equivalents [equivalents (eq.)]/kg body weight of A. flosaquae DC-1 aphantoxins. Hyperemia, the hepatosomatic index (HSI), and physiological and molecular responses of pro- and anti-inflammatory cytokines in the zebrafish liver were investigated at different time points 1-24 h post-exposure. Aphantoxins significantly enhanced hepatic hyperemia and altered the HSI 3-24 h post-exposure, suggesting that inflammation caused morphological changes. Subsequent investigations using the enzyme-linked immunosorbent assay showed that the pro-inflammatory cytokines tumor necrosis factor-α, interleukin-1β (IL-1β), IL-6, and IL-8 and anti-inflammatory cytokines IL-10 and transforming growth factor β were higher in the liver of zebrafish exposed to aphantoxins, which indicated physiological inflammatory responses. Further analysis by real-time fluorescence quantitative polymerase chain reaction demonstrated upregulated mRNA expression of these cytokines, suggesting molecular inflammatory responses in the zebrafish liver. These changes showed dose- and time-dependent patterns. These results indicated that aphantoxins induced hyperemia and altered the HSI, and subsequently increased the levels of proinflammatory cytokines TNF-α, IL-1β, IL-6 and IL-8 to induce physiological inflammatory responses. These changes activated the anti-inflammatory cytokines IL-10 and TGF-β to suppress inflammatory damage. The induced changes were the result of upregulated mRNA expression of these inflammatory cytokines caused by aphantoxins. Aphantoxins resulted in hepatic immunotoxicity and response by inducing pro-inflammatory cytokines. Zebrafish liver in turn suppressed the inflammatory damage by upregulating the activities of anti-inflammatory cytokines. In the future, these pro- and anti-inflammatory cytokines in the zebrafish liver may be prove to be useful biomarkers of aphantoxins and blooms in nature.

Models of cytokine dynamics in the inflammatory response of viral zoonotic infectious diseases

Inflammatory responses to an infection from a zoonotic pathogen, such as avian influenza viruses, hantaviruses and some coronaviruses, are distinctly different in their natural reservoir versus human host. While not as well studied in the natural reservoirs, the pro-inflammatory response and viral replication appear controlled and show no obvious pathology. In contrast, infection in humans results in an initial high viral load marked by an aggressive pro-inflammatory response known as a cytokine storm. The key difference in the course of the infection between the reservoir and human host is the inflammatory response. In this investigation, we apply a simple two-component differential equation model for pro-inflammatory and anti-inflammatory responses and a detailed mathematical analysis to identify specific regions in parameter space for single stable endemic equilibrium, bistability or periodic solutions. The extensions of the deterministic model to two stochastic models account for variability in responses seen at the cell (local) or tissue (global) levels. Numerical solutions of the stochastic models exhibit outcomes that are typical of a chronic infection in the natural reservoir or a cytokine storm in human infection. In the chronic infection, occasional flare-ups between high and low responses occur when model parameters are in a region of bistability or periodic solutions. The cytokine storm with a vigorous pro-inflammatory response and less vigorous anti-inflammatory response occurs in the parameter region for a single stable endemic equilibrium with a strong pro-inflammatory response. The results of the model analyses and the simulations are interpreted in terms of the functional role of the cytokines and the inflammatory responses seen in infection of the natural reservoir or of the human host.

Feature Engineering for ICU Mortality Prediction Based on Hourly to Bi-Hourly Measurements

Mortality prediction for intensive care unit (ICU) patients is a challenging problem that requires extracting discriminative and informative features. This study presents a proof of concept for exploring features that can provide clinical insight. Through a feature engineering approach, it is attempted to improve ICU mortality prediction in field conditions with low frequently measured data (i.e., hourly to bi-hourly). Features are explored by investigating the vital signs measurements of ICU patients, labelled with mortality or survival at discharge. The vital signs of interest in this study are heart and respiration rate, oxygen saturation and blood pressure. The latter comprises systolic, diastolic and mean arterial pressure. In the feature exploration process, it is aimed to extract simple and interpretable features that can provide clinical insight. For this purpose, a classifier is required that maximises the margin between the two classes (i.e., survival and mortality) with minimum tolerance to misclassification errors. Moreover, it preferably has to provide a linear decision surface in the original feature space without mapping to an unlimited dimensionality feature space. Therefore, a linear hard margin support vector machine (SVM) classifier is suggested. The extracted features are grouped in three categories: statistical, dynamic and physiological. Each category plays an important role in enhancing classification error performance. After extracting several features within the three categories, a manual feature fine-tuning is applied to consider only the most efficient features. The final classification, considering mortality as the positive class, resulted in an accuracy of 91.56 % , sensitivity of 90.59 % , precision of 86.52 % and F 1 -score of 88.50 % . The obtained results show that the proposed feature engineering approach and the extracted features are valid to be considered and further enhanced for the mortality prediction purpose. Moreover, the proposed feature engineering approach moved the modelling methodology from black-box modelling to grey-box modelling in combination with the powerful classifier of SVMs.

Identifying important parameters in the inflammatory process with a mathematical model of immune cell influx and macrophage polarization

In an inflammatory setting, macrophages can be polarized to an inflammatory M1 phenotype or to an anti-inflammatory M2 phenotype, as well as existing on a spectrum between these two extremes. Dysfunction of this phenotypic switch can result in a population imbalance that leads to chronic wounds or disease due to unresolved inflammation. Therapeutic interventions that target macrophages have therefore been proposed and implemented in diseases that feature chronic inflammation such as diabetes mellitus and atherosclerosis. We have developed a model for the sequential influx of immune cells in the peritoneal cavity in response to a bacterial stimulus that includes macrophage polarization, with the simplifying assumption that macrophages can be classified as M1 or M2. With this model, we were able to reproduce the expected timing of sequential influx of immune cells and mediators in a general inflammatory setting. We then fit this model to in vivo experimental data obtained from a mouse peritonitis model of inflammation, which is widely used to evaluate endogenous processes in response to an inflammatory stimulus. Model robustness is explored with local structural and practical identifiability of the proposed model a posteriori. Additionally, we perform sensitivity analysis that identifies the population of apoptotic neutrophils as a key driver of the inflammatory process. Finally, we simulate a selection of proposed therapies including points of intervention in the case of delayed neutrophil apoptosis, which our model predicts will result in a sustained inflammatory response. Our model can therefore provide hypothesis testing for therapeutic interventions that target macrophage phenotype and predict outcomes to be validated by subsequent experimentation.

Using experimental data and mathematical analysis, we develop a model for the inflammatory response that includes macrophage polarization between M1 and M2 phenotypes. Dysfunction of this phenotypic switch can disrupt the timely influx and egress of immune cells during the healing process and lead to chronic wounds or disease. The modulation of macrophage population has been suggested as a strategy to dampen inflammation in diseases that feature chronic inflammation, such as diabetes and atherosclerosis. It is therefore important that we learn more about which components of the system drive the population level switch in phenotype. Our model is able to reproduce the expected timing of sequential influx of neutrophils and macrophages in response to an inflammatory stimulus. Model parameters were estimated with weighted least squares fitting to in vivo experimental data from a mouse model of peritonitis while considering identifiability of parameter sets. We perform sensitivity analysis that identifies primary drivers of the system, and predict the effects of variations in these key parameters on immune cell populations.

A mathematical model of the use of supplemental oxygen to combat surgical site infection

Infections are a common complication of any surgery, often requiring a recovery period in hospital. Supplemental oxygen therapy administered during and immediately after surgery is thought to enhance the immune response to bacterial contamination. However, aerobic bacteria thrive in oxygen-rich environments, and so it is unclear whether oxygen has a net positive effect on recovery. Here, we develop a mathematical model of post-surgery infection to investigate the efficacy of supplemental oxygen therapy on surgical-site infections. A 4-species, coupled, set of non-linear partial differential equations that describes the space-time dependence of neutrophils, bacteria, chemoattractant and oxygen is developed and analysed to determine its underlying properties. Through numerical solutions, we quantify the efficacy of different supplemental oxygen regimes on the treatment of surgical site infections in wounds of different initial bacterial load. A sensitivity analysis is performed to investigate the robustness of the predictions to changes in the model parameters. The numerical results are in good agreement with analyses of the associated well-mixed model. Our model findings provide insight into how the nature of the contaminant and its initial density influence bacterial infection dynamics in the surgical wound.

Mathematical modeling of septic shock based on clinical data

BACKGROUND

Mathematical models of diseases may provide a unified approach for establishing effective treatment strategies based on fundamental pathophysiology. However, models that are useful for clinical practice must overcome the massive complexity of human physiology and the diversity of patients' environmental conditions. With the aim of modeling a complex disease, we choose sepsis, which is highly complex, life-threatening systemic disease with high mortality. In particular, we focused on septic shock, a subset of sepsis in which underlying circulatory and cellular/metabolic abnormalities are profound enough to substantially increase mortality. Our model includes cardiovascular, immune, nervous system models and a pharmacological model as submodels and integrates them to create a sepsis model based on pathological facts.

RESULTS

Model validation was done in two steps. First, we established a model for a standard patient in order to confirm the validity of our approach in general aspects. For this, we checked the correspondence between the severity of infection defined in terms of pathogen growth rate and the ease of recovery defined in terms of the intensity of treatment required for recovery. The simulations for a standard patient showed good correspondence. We then applied the same simulations to a patient with heart failure as an underlying disease. The model showed that spontaneous recovery would not occur without treatment, even for a very mild infection. This is consistent with clinical experience. We next validated the model using clinical data of three sepsis patients. The model parameters were tuned for these patients based on the model for the standard patient used in the first part of the validation. In these cases, the simulations agreed well with clinical data. In fact, only a handful parameters need to be tuned for the simulations to match with the data.

CONCLUSIONS

We have constructed a model of septic shock and have shown that it can reproduce well the time courses of treatment and disease progression. Tuning of model parameters for each patient could be easily done. This study demonstrates the feasibility of disease models, suggesting the possibility of clinical use in the prediction of disease progression, decisions on the timing of drug dosages, and the estimation of time of infection.

Mathematical modeling of energy consumption in the acute inflammatory response

When a pathogen invades the body, an acute inflammatory response is activated to eliminate the intruder. In some patients, runaway activation of the immune system may lead to collateral tissue damage and, in the extreme, organ failure and death. Experimental studies have found an association between severe infections and depletion in levels of adenosine triphosphate (ATP), increase in nitric oxide production, and accumulation of lactate, suggesting that tissue energetics is compromised. In this work we present a differential equations model that incorporates the dynamics of ATP, nitric oxide, and lactate accompanying an acute inflammatory response and employ this model to explore their roles in shaping this response. The bifurcation diagram of the model system with respect to the pathogen growth rate reveals three equilibrium states characterizing the health, aseptic and septic conditions. We explore the domains of attraction of these states to inform the instantiation of heterogeneous virtual patient populations utilized in a survival analysis. We then apply the model to study alterations in the inflammatory response and survival outcomes in metabolically altered conditions such as hypoglycemia, hyperglycemia, and hypoxia.

A mathematical model of neuroinflammation in severe clinical traumatic brain injury

Background

Understanding the interdependencies among inflammatory mediators of tissue damage following traumatic brain injury (TBI) is essential in providing effective, patient-specific care. Activated microglia and elevated concentrations of inflammatory signaling molecules reflect the complex cascades associated with acute neuroinflammation and are predictive of recovery after TBI. However, clinical TBI studies to date have not focused on modeling the dynamic temporal patterns of simultaneously evolving inflammatory mediators, which has potential in guiding the design of future immunomodulation intervention studies.

Methods

We derived a mathematical model consisting of ordinary differential equations (ODE) to represent interactions between pro- and anti-inflammatory cytokines, M1- and M2-like microglia, and central nervous system (CNS) tissue damage. We incorporated variables for several cytokines, interleukin (IL)-1β, IL-4, IL-10, and IL-12, known to have roles in microglial activation and phenotype differentiation. The model was fit to cerebrospinal fluid (CSF) cytokine data, collected during the first 5 days post-injury in n = 89 adults with severe TBI. Ensembles of model fits were produced for three patient subgroups: (1) a favorable outcome group (GOS = 4,5) and (2) an unfavorable outcome group (GOS = 1,2,3) both with lower pro-inflammatory load, and (3) an unfavorable outcome group (GOS = 1,2,3) with higher pro-inflammatory load. Differences in parameter distributions between subgroups were ranked using Bhattacharyya metrics to identify mechanistic differences underlying the neuroinflammatory patterns of patient groups with different TBI outcomes.

Results

Optimal model fits to data showed different microglial and damage responses by patient subgroup. Upon comparison of model parameter distributions, unfavorable outcome groups were characterized by either a prolonged, pathophysiological or a transient, sub-physiological course of neuroinflammation.

Conclusion

By developing a mathematical characterization of inflammatory processes informed by clinical data, we have created a system for exploring links between acute neuroinflammatory components and patient outcome in severe TBI.

Electronic supplementary material

The online version of this article (10.1186/s12974-018-1384-1) contains supplementary material, which is available to authorized users.

Inflammation and Disease: Modelling and Modulation of the Inflammatory Response to Alleviate Critical Illness

Critical illness, a constellation of interrelated inflammatory and physiological derangements occurring subsequent to severe infection or injury, affects a large number of individuals in both developed and developing countries. The prototypical complex system embodied in critical illness has largely defied therapy beyond supportive care. We have focused on the utility of data-driven and mechanistic computational modelling to help address the complexity of critical illness and provide pathways towards discovering potential therapeutic options and combinations. Herein, we review recent progress in this field, with a focus on both animal and computational models of critical illness. We suggest that therapy for critical illness can be posed as a model-based dynamic control problem, and discuss novel theoretical and experimental approaches involving biohybrid devices aimed at reprogramming inflammation dynamically. Together, these advances offer the potential for Model-based Precision Medicine for critical illness.

Modeling the Bistable Dynamics of the Innate Immune System

The size of primary challenge with lipopolysaccharide induces changes in the innate immune cells phenotype between pro-inflammatory and pro-tolerant states when facing a secondary lipopolysaccharide challenge. To determine the molecular mechanisms governing this differential response, we propose a mathematical model for the interaction between three proteins involved in the immune cell decision making: IRAK-1, PI3K, and RelB. The mutual inhibition of IRAK-1 and PI3K in the model leads to bistable dynamics. By using the levels of RelB as indicative of strength of the immune responses, we connect the size of different primary lipopolysaccharide doses to the differential phenotypical outcomes following a secondary challenge. We further predict under what circumstances the primary LPS dose does not influence the response to a secondary challenge. Our results can be used to guide treatments for patients with either autoimmune disease or compromised immune system.

Supplemented Alkaline Phosphatase Supports the Immune Response in Patients Undergoing Cardiac Surgery: Clinical and Computational Evidence

Alkaline phosphatase (AP) is an enzyme that exhibits anti-inflammatory effects by dephosphorylating inflammation triggering moieties (ITMs) like bacterial lipopolysaccharides and extracellular nucleotides. AP administration aims to prevent and treat peri- and post-surgical ischemia reperfusion injury in cardiothoracic surgery patients. Recent studies reported that intravenous bolus administration and continuous infusion of AP in patients undergoing coronary artery bypass grafting with cardiac valve surgery induce an increased release of liver-type “tissue non-specific alkaline phosphatase” (TNAP) into the bloodstream. The release of liver-type TNAP into circulation could be the body's way of strengthening its defense against a massive ischemic insult. However, the underlying mechanism behind the induction of TNAP is still unclear. To obtain a deeper insight into the role of AP during surgery, we developed a mathematical model of systemic inflammation that clarifies the relation between supplemented AP and TNAP and describes a plausible induction mechanism of TNAP in patients undergoing cardiothoracic surgery. The model was validated against clinical data from patients treated with bovine Intestinal AP (bIAP treatment) or without AP (placebo treatment), in addition to standard care procedures. We performed additional in-silico experiments adding a secondary source of ITMs after surgery, as observed in some patients with complications, and predicted the response to different AP treatment regimens. Our results show a strong protective effect of supplemented AP for patients with complications. The model provides evidence of the existence of an induction mechanism of liver-type tissue non-specific alkaline phosphatase, triggered by the supplementation of AP in patients undergoing cardiac surgery. To the best of our knowledge this is the first time that a quantitative and validated numerical model of systemic inflammation under clinical treatment conditions is presented.

A review of inflammatory mechanism in airway diseases

BACKGROUND

Inflammation in the lung is the body's natural response to injury. It acts to remove harmful stimuli such as pathogens, irritants, and damaged cells and initiate the healing process. Acute and chronic pulmonary inflammation are seen in different respiratory diseases such as; acute respiratory distress syndrome, chronic obstructive pulmonary disease (COPD), asthma, and cystic fibrosis (CF).

FINDINGS

In this review, we found that inflammatory response in COPD is determined by the activation of epithelial cells and macrophages in the respiratory tract. Epithelial cells and macrophages discharge transforming growth factor-β (TGF-β), which trigger fibroblast proliferation and tissue remodeling. Asthma leads to airway hyper-responsiveness, obstruction, mucus hyper-production, and airway-wall remodeling. Cytokines, allergens, chemokines, and infectious agents are the main stimuli that activate signaling pathways in epithelial cells in asthma. Mutation of the CF transmembrane conductance regulator (CFTR) gene results in CF. Mutations in CFTR influence the lung epithelial innate immune function that leads to exaggerated and ineffective airway inflammation that fails to abolish pulmonary pathogens. We present mechanistic computational models (based on ordinary differential equations, partial differential equations and agent-based models) that have been applied in studying the complex physiological and pathological mechanisms of chronic inflammation in different airway diseases.

CONCLUSION

The scope of the present review is to explore the inflammatory mechanism in airway diseases and highlight the influence of aging on airways' inflammation mechanism. The main goal of this review is to encourage research collaborations between experimentalist and modelers to promote our understanding of the physiological and pathological mechanisms that control inflammation in different airway diseases.