John Bartels

Translational Systems Biology of Inflammation

Inflammation is a complex, multi-scale biologic response to stress that is also required for repair and regeneration after injury. Despite the repository of detailed data about the cellular and molecular processes involved in inflammation, including some understanding of its pathophysiology, little progress has been made in treating the severe inflammatory syndrome of sepsis. To address the gap between basic science knowledge and therapy for sepsis, a community of biologists and physicians is using systems biology approaches in hopes of yielding basic insights into the biology of inflammation. “Systems biology” is a discipline that combines experimental discovery with mathematical modeling to aid in the understanding of the dynamic global organization and function of a biologic system (cell to organ to organism). We propose the term translational systems biology for the application of similar tools and engineering principles to biologic systems with the primary goal of optimizing clinical practice. We describe the efforts to use translational systems biology to develop an integrated framework to gain insight into the problem of acute inflammation. Progress in understanding inflammation using translational systems biology tools highlights the promise of this multidisciplinary field. Future advances in understanding complex medical problems are highly dependent on methodological advances and integration of the computational systems biology community with biologists and clinicians.

THE ACUTE INFLAMMATORY RESPONSE IN DIVERSE SHOCK STATES

A poorly controlled acute inflammatory response can lead to organ dysfunction and death. Severe systemic inflammation can be induced and perpetuated by diverse insults such as the administration of toxic bacterial products (e.g., endotoxin), traumatic injury, and hemorrhage. Here, we probe whether these varied shock states can be explained by a universal inflammatory system that is initiated through different means and, once initiated, follows a course specified by the cellular and molecular mechanisms of the immune and endocrine systems. To examine this question, we developed a mathematical model incorporating major elements of the acute inflammatory response in C57Bl/6 mice, using input from experimental data. We found that a single model with different initiators including the autonomic system could describe the response to various insults. This model was able to predict a dose range of endotoxin at which mice would die despite having been calibrated only in nonlethal inflammatory paradigms. These results show that the complex biology of inflammation can be modeled and supports the hypothesis that shock states induced by a range of physiologic challenges could arise from a universal response that is differently initiated and modulated.

IN SILICO MODELS OF ACUTE INFLAMMATION IN ANIMALS

ABSTRACT Trauma and hemorrhagic shock elicit an acute inflammatory response, predisposing patients to sepsis, organ dysfunction, and death. Few approved therapies exist for these acute inflammatory states, mainly due to the complex interplay of interacting inflammatory and physiological elements working at multiple levels. Various animal models have been used to simulate these phenomena, but these models often do not replicate the clinical setting of multiple overlapping insults. Mathematical modeling of complex systems is an approach for understanding the interplay among biological interactions. We constructed a mathematical model using ordinary differential equations that encompass the dynamics of cells and cytokines of the acute inflammatory response, as well as global tissue dysfunction. The model was calibrated in C57Bl/6 mice subjected to (1) various doses of lipopolysaccharide (LPS) alone, (2) surgical trauma, and (3) surgery + hemorrhagic shock. We tested the models predictive ability in scenarios on which it had not been trained, namely, (1) surgery ± hemorrhagic shock + LPS given at times after the beginning of surgical instrumentation, and (2) surgery + hemorrhagic shock + bilateral femoral fracture. Software was created that facilitated fitting of the mathematical model to experimental data, as well as for simulation of experiments with various inflammatory challenges and associated variations (gene knockouts, inhibition of specific cytokines, etc.). Using this software, the C57Bl/6-specific model was recalibrated for inflammatory analyte data in CD14−/− mice and was used to elucidate altered features of inflammation in these animals. In other experiments, rats were subjected to surgical trauma ± LPS or to bacterial infection via fibrin clots impregnated with various inocula of Escherichia coli. Mathematical modeling may provide insights into the complex dynamics of acute inflammation in a manner that can be tested in vivo using many fewer animals than has been possible previously.

The role of initial trauma in the host's response to injury and hemorrhage: insights from a correlation of mathematical simulations and hepatic transcriptomic analysis.

ABSTRACT Trauma and hemorrhagic shock (HS) elicit severe physiological disturbances that predispose the victims to subsequent organ dysfunction and death. The general lack of effective therapeutic options for these patients is mainly due to the complex interplay of interacting inflammatory and physiological elements working at multiple levels. Systems biology has emerged as a new paradigm that allows the study of large portions of physiological networks simultaneously. Seeking a better understanding of the interplay among known inflammatory pathways, we constructed a mathematical model encompassing the dynamics of the acute inflammatory response that incorporates the intertwined effects of inflammation and global tissue damage. The model was calibrated using data from C57Bl/6 mice subjected to endotoxemia, sham operation (i.e., surgical trauma induced by cannulation [ST]) or ST + HS+ resuscitation (ST-HS-R). An in silico simulation, made at whole-organism level, suggested that similar pathways of different magnitudes were operant as the degree of total body damage increased. We sought to validate this hypothesis by subjecting mice to HS and comparing the models predictions to circulating markers of inflammation and tissue injury as well as the global transcriptomic response of the liver. C57Bl/6 mice were subjected to ST or ST-HS (without resuscitation). Liver gene expression was assessed using an Affymetrix DNA microarray (GeneChip Mouse Expression Set 430A, Affymetrix, Santa Clara, CA), which contains 22,621 probe sets and effectively interrogates 12,341 mouse genes. The microarray data sets were subjected to hierarchical clustering and pathway analysis. In agreement with model predictions, circulating levels of inflammation/tissue injury markers and the microarray analysis both demonstrated that ST alone accounts for a substantial proportion of the observed phenotypic and genetic/molecular changes versus untreated animals. The addition of HS further increased the magnitude of gene expression, but relatively few additional genes were recruited. Mathematical simulations and DNA microarrays, both systems biology tools, may provide valuable insight into the complex global physiological interactions that occur in response to trauma and hemorrhagic shock.ABBREVIATIONS-ST-surgical trauma, HS-hemorrhagic shock, ST-HS-surgical trauma + hemorrhagic shock, R-resuscitation, MAP-mean arterial blood pressure, IL-interleukin, TNF-&agr;-tumor necrosis factor alpha, AST-aspartate aminotransferase, iNOS-inducible nitric oxide synthase, eNOS-endothelial nitric oxide synthase, NO-nitric oxide, NO2−-nitrate, NO3−-nitrite, LPS-bacterial lipopolysaccharide (endotoxin), IP-ingenuity pathway, IPKB-ingenuity pathway knowledge base

Translational systems approaches to the biology of inflammation and healing.

Inflammation is a complex, non-linear process central to many of the diseases that affect both developed and emerging nations. A systems-based understanding of inflammation, coupled to translational applications, is therefore necessary for efficient development of drugs and devices, for streamlining analyses at the level of populations, and for the implementation of personalized medicine. We have carried out an iterative and ongoing program of literature analysis, generation of prospective data, data analysis, and computational modeling in various experimental and clinical inflammatory disease settings. These simulations have been used to gain basic insights into the inflammatory response under baseline, gene-knockout, and drug-treated experimental animals for in silico studies associated with the clinical settings of sepsis, trauma, acute liver failure, and wound healing to create patient-specific simulations in polytrauma, traumatic brain injury, and vocal fold inflammation; and to gain insight into host-pathogen interactions in malaria, necrotizing enterocolitis, and sepsis. These simulations have converged with other systems biology approaches (e.g., functional genomics) to aid in the design of new drugs or devices geared towards modulating inflammation. Since they include both circulating and tissue-level inflammatory mediators, these simulations transcend typical cytokine networks by associating inflammatory processes with tissue/organ impacts via tissue damage/dysfunction. This framework has now allowed us to suggest how to modulate acute inflammation in a rational, individually optimized fashion. This plethora of computational and intertwined experimental/engineering approaches is the cornerstone of Translational Systems Biology approaches for inflammatory diseases.

In Silico Augmentation of the Drug Development Pipeline: Examples from the study of Acute Inflammation

The clinical translation of promising basic biomedical findings, whether derived from reductionist studies in academic laboratories or as the product of extensive high‐throughput and high‐content screens in the biotechnology and pharmaceutical industries, has reached a period of stagnation in which ever higher research and development costs are yielding ever fewer new drugs. Systems biology and computational modeling have been touted as potential avenues by which to break through this logjam. However, few mechanistic computational approaches are used in a manner that is fully cognizant of the inherent clinical realities in which the drugs developed through this ostensibly rational process will be ultimately used. This article presents a translational systems biology approach to inflammation. This approach is based on the use of mechanistic computational modeling centered on inherent clinical applicability, namely that a unified suite of models can be applied to generate in silico clinical trials, individualized computational models as tools for personalized medicine, and rational drug and device design based on disease mechanism. Drug Dev Res 72: 187–200, 2011.   © 2010 Wiley‐Liss, Inc.

Mathematical modeling of posthemorrhage inflammation in mice: studies using a novel, computer-controlled, closed-loop hemorrhage apparatus.

Hemorrhagic shock (HS) elicits a global acute inflammatory response, organ dysfunction, and death. We have used mathematical modeling of inflammation and tissue damage/dysfunction to gain insight into this complex response in mice. We sought to increase the fidelity of our mathematical model and to establish a platform for testing predictions of this model. Accordingly, we constructed a computerized, closed-loop system for mouse HS. The intensity, duration, and time to achieve target MAP could all be controlled using a software. Fifty-four male C57/black mice either were untreated or underwent surgical cannulation. The cannulated mice were divided into 8 groups: (a) 1, 2, 3, or 4 h of surgical cannulation alone and b) 1, 2, 3, or 4 h of cannulation + HS (25 mmHg). MAP was sustained by the computer-controlled reinfusion and withdrawal of shed blood within ±2 mmHg. Plasma was assayed for the cytokines TNF, IL-6, and IL-10 as well as the NO reaction products NO2−/NO3−. The cytokine and NO2−/NO3− data were compared with predictions from a mathematical model of post-hemorrhage inflammation, which was calibrated on different data. To varying degrees, the levels of TNF, IL-6, IL-10, and NO2−/NO3− predicted by the mathematical model matched these data closely. In conclusion, we have established a hardware/software platform that allows for highly accurate, reproducible, and mathematically predictable HS in mice.ABBREVIATIONS-HS-hemorrhagic shock; HR-heart rate

Combined In Silico, In Vivo, and In Vitro Studies Shed Insights into the Acute Inflammatory Response in Middle-Aged Mice

We combined in silico, in vivo, and in vitro studies to gain insights into age-dependent changes in acute inflammation in response to bacterial endotoxin (LPS). Time-course cytokine, chemokine, and NO2 −/NO3 − data from “middle-aged” (6–8 months old) C57BL/6 mice were used to re-parameterize a mechanistic mathematical model of acute inflammation originally calibrated for “young” (2–3 months old) mice. These studies suggested that macrophages from middle-aged mice are more susceptible to cell death, as well as producing higher levels of pro-inflammatory cytokines, vs. macrophages from young mice. In support of the in silico-derived hypotheses, resident peritoneal cells from endotoxemic middle-aged mice exhibited reduced viability and produced elevated levels of TNF-α, IL-6, IL-10, and KC/CXCL1 as compared to cells from young mice. Our studies demonstrate the utility of a combined in silico, in vivo, and in vitro approach to the study of acute inflammation in shock states, and suggest hypotheses with regard to the changes in the cytokine milieu that accompany aging.

Mathematical and Statistical Modeling of Acute Inflammation

A mathematical model involving a system of ordinary differential equations has been developed with the goal of assisting the design of therapies directed against the inflammatory consequences of infection and trauma. Though the aim is to build a model of greater complexity, which would eventually overcome some existing limitations (such as a reduced subset of inflammatory interactions, the use of mass action kinetics, and calibration to circulating but not local levels of cytokines), the model can at this time simulate certain disease scenarios qualitatively as well as predicting the time course of cytokine levels in distinct paradigms of inflammation in mice. A parameter search algorithm is developed that aids in the identification of different regimes of behaviour of the model and helps with its calibration to data. Extending this mathematical model, with validation in humans, may lead to the in silico development of novel therapeutic approaches and real-time diagnostics.

A linear code parameter search algorithm with applications to immunology

Abstract The immune system is modeled by way of a system of ordinary differential equations involving a large number of parameters, such as growth rates and initial conditions. Key to successful implementation of the model is the estimation of such parameters from available data. A parameter search algorithm based on linear codes is developed having as aim the identification of different regimes of behaviour of the model, the estimation of parameters in a high dimensional space, and the model calibration to data.