Translational potential of systems-based models of inflammation

Computational modeling of the synergistic role of GCN2 and the HPA axis in regulating the integrated stress response in the central circadian timing system

The circadian timing system and integrated stress response (ISR) systems are fundamental regulatory mechanisms that maintain body homeostasis. The central circadian pacemaker in the suprachiasmatic nucleus (SCN) governs daily rhythms through interactions with peripheral oscillators via the hypothalamus-pituitary-adrenal (HPA) axis. On the other hand, ISR signaling is pivotal for preserving cellular homeostasis in response to physiological changes. Notably, disrupted circadian rhythms are observed in cases of impaired ISR signaling. In this work, we examine the potential interplay between the central circadian system and the ISR, mainly through the SCN and HPA axis. We introduce a semimechanistic mathematical model to delineate SCN's capacity for indirectly perceiving physiological stress through glucocorticoid-mediated feedback from the HPA axis and orchestrating a cellular response via the ISR mechanism. Key components of our investigation include evaluating general control nonderepressible 2 (GCN2) expression in the SCN, the effect of physiological stress stimuli on the HPA axis, and the interconnected feedback between the HPA and SCN. Simulation revealed a critical role for GCN2 in linking ISR with circadian rhythms. Experimental findings have demonstrated that a Gcn2 deletion in mice leads to rapid re-entrainment of the circadian clock following jetlag as well as to an elongation of the circadian period. These phenomena are well replicated by our model, which suggests that both the swift re-entrainment and prolonged period can be ascribed to a reduced robustness in neuronal oscillators. Our model also offers insights into phase shifts induced by acute physiological stress and the alignment/misalignment of physiological stress with external light-dark cues. Such understanding aids in strategizing responses to stressful events, such as nutritional status changes and jetlag.NEW & NOTEWORTHY This study is the first theoretical work to investigate the complex interaction between integrated stress response (ISR) sensing and central circadian rhythm regulation, encompassing the suprachiasmatic nucleus (SCN) and hypothalamus-pituitary-adrenal (HPA) axis. The findings carry implications for the development of dietary or pharmacological interventions aimed at facilitating recovery from stressful events, such as jetlag. Moreover, they provide promising prospects for potential therapeutic interventions that target circadian rhythm disruption and various stress-related disorders.

The Physiological and Pharmacological Significance of the Circadian Timing of the HPA Axis: A Mathematical Modeling Approach

As a potent endogenous regulator of homeostasis, the circadian time-keeping system synchronizes internal physiology to periodic changes in the external environment to enhance survival. Adapting endogenous rhythms to the external time is accomplished hierarchically with the central pacemaker located in the suprachiasmatic nucleus (SCN) signaling the hypothalamus-pituitary-adrenal (HPA) axis to release hormones, notably cortisol, which help maintain the body's circadian rhythm. Given the essential role of HPA-releasing hormones in regulating physiological functions, including immune response, cell cycle, and energy metabolism, their daily variation is critical for the proper function of the circadian timing system. In this review, we focus on cortisol and key fundamental properties of the HPA axis and highlight their importance in controlling circadian dynamics. We demonstrate how systems-driven, mathematical modeling of the HPA axis complements experimental findings, enhances our understanding of complex physiological systems, helps predict potential mechanisms of action, and elucidates the consequences of circadian disruption. Finally, we outline the implications of circadian regulation in the context of personalized chronotherapy. Focusing on the chrono-pharmacology of synthetic glucocorticoids, we review the challenges and opportunities associated with moving toward personalized therapies that capitalize on circadian rhythms.

Light-induced synchronization of the SCN coupled oscillators and implications for entraining the HPA axis

The suprachiasmatic nucleus (SCN) synchronizes the physiological rhythms to the external light-dark cycle and tunes the dynamics of circadian rhythms to photoperiod fluctuations. Changes in the neuronal network topologies are suggested to cause adaptation of the SCN in different photoperiods, resulting in the broader phase distribution of neuron activities in long photoperiods (LP) compared to short photoperiods (SP). Regulated by the SCN output, the level of glucocorticoids is elevated in short photoperiod, which is associated with peak disease incidence. The underlying coupling mechanisms of the SCN and the interplay between the SCN and the HPA axis have yet to be fully elucidated. In this work, we propose a mathematical model including a multiple-cellular SCN compartment and the HPA axis to investigate the properties of the circadian timing system under photoperiod changes. Our model predicts that the probability-dependent network is more energy-efficient than the distance-dependent network. Coupling the SCN network by intra-subpopulation and inter-subpopulation forces, we identified the negative correlation between robustness and plasticity of the oscillatory network. The HPA rhythms were predicted to be strongly entrained to the SCN rhythms with a pro-inflammatory high-amplitude glucocorticoid profile under SP. The fast temporal topology switch of the SCN network was predicted to enhance synchronization when the synchronization is not complete. These synchronization and circadian dynamics alterations might govern the seasonal variation of disease incidence and its symptom severity.

Teaching computational systems biology with an eye on quantitative systems pharmacology at the undergraduate level: Why do it, who would take it, and what should we teach?

Computational systems biology (CSB) is a field that emerged primarily as the product of research activities. As such, it grew in several directions in a distributed and uncoordinated manner making the area appealing and fascinating. The idea of not having to follow a specific path but instead creating one fueled innovation. As the field matured, several interdisciplinary graduate programs emerged attempting to educate future generations of computational systems biologists. These educational initiatives coordinated the dissemination of information across student populations that had already decided to specialize in this field. However, we are now entering an era where CSB, having established itself as a valuable research discipline, is attempting the next major step: Entering undergraduate curricula. As interesting as this endeavor may sound, it has several difficulties, mainly because the field is not uniformly defined. In this manuscript, we argue that this diversity is a significant advantage and that several incarnations of an undergraduate-level CSB biology course could, and should, be developed tailored to programmatic needs. In this manuscript, we share our experiences creating a course as part of a Biomedical Engineering program.

On the analysis of complex biological supply chains: From Process Systems Engineering to Quantitative Systems Pharmacology

The use of models in biology has become particularly relevant as it enables investigators to develop a mechanistic framework for understanding the operating principles of living systems as well as in quantitatively predicting their response to both pathological perturbations and pharmacological interventions. This application has resulted in a synergistic convergence of systems biology and pharmacokinetic-pharmacodynamic modeling techniques that has led to the emergence of quantitative systems pharmacology (QSP). In this review, we discuss how the foundational principles of chemical process systems engineering inform the progressive development of more physiologically-based systems biology models.

Determining the role of inflammation in the selection of JAK2 mutant cells in myeloproliferative neoplasms

Myeloproliferative neoplasm (MPN) is a hematologic malignancy characterized by the clonal outgrowth of hematopoietic cells with a somatically acquired mutation most commonly in JAK2 (JAK2^V617F). This mutation endows upon myeloid progenitors cytokine independent growth and consequently leads to excessive production of myeloid lineage cells. It has been previously suggested that inflammation may play a role in the clonal evolution of JAK2^V617F mutants. In particular, it is possible that one or more cellular kinetic parameters of hematopoietic stem cells (HSCs) are affected by inflammation, such as division or death rates of cells, and the probability of HSC differentiation. This suggests a mechanism that can steer the outcome of the cellular competition in favor of the mutants, initiating the disease. In this paper we create a number of mathematical evolutionary models, from very abstract to more concrete, that describe cellular competition in the context of inflammation. It is possible to build a model axiomatically, where only very general assumptions are imposed on the modeling components and no arbitrary (and generally unknown) functional forms are used, and still generate a set of testable predictions. In particular, we show that, if HSC death is negligible, the evolutionary advantage of mutant cells can only be conferred by an increase in differentiation probability of HSCs in the presence of inflammation, and if death plays a significant role in the dynamics, an additional mechanism may be an increase of HSC's division-to-death ratio in the presence of inflammation. Further, we show that in the presence of inflammation, the wild type cell population is predicted to shrink under inflammation (even in the absence of mutants). Finally, it turns out that if only the differentiation probability is affected by the inflammation, then the resulting steady state population of wild type cells will contain a relatively smaller percentage of HSCs under inflammation. If the division-to-death rate is also affected, then the percentage of HSCs under inflammation can either decrease or increase, depending on other parameters.

Understanding Physiology in the Continuum: Integration of Information from Multiple -Omics Levels

In this paper, we discuss approaches for integrating biological information reflecting diverse physiologic levels. In particular, we explore statistical and model-based methods for integrating transcriptomic, proteomic and metabolomics data. Our case studies reflect responses to a systemic inflammatory stimulus and in response to an anti-inflammatory treatment. Our paper serves partly as a review of existing methods and partly as a means to demonstrate, using case studies related to human endotoxemia and response to methylprednisolone (MPL) treatment, how specific questions may require specific methods, thus emphasizing the non-uniqueness of the approaches. Finally, we explore novel ways for integrating -omics information with PKPD models, toward the development of more integrated pharmacology models.

Augmenting Surgery via Multi-scale Modeling and Translational Systems Biology in the Era of Precision Medicine: A Multidisciplinary Perspective

In this era of tremendous technological capabilities and increased focus on improving clinical outcomes, decreasing costs, and increasing precision, there is a need for a more quantitative approach to the field of surgery. Multiscale computational modeling has the potential to bridge the gap to the emerging paradigms of Precision Medicine and Translational Systems Biology, in which quantitative metrics and data guide patient care through improved stratification, diagnosis, and therapy. Achievements by multiple groups have demonstrated the potential for (1) multiscale computational modeling, at a biological level, of diseases treated with surgery and the surgical procedure process at the level of the individual and the population; along with (2) patient-specific, computationally-enabled surgical planning, delivery, and guidance and robotically-augmented manipulation. In this perspective article, we discuss these concepts, and cite emerging examples from the fields of trauma, wound healing, and cardiac surgery.

A Chemical Engineer's Perspective on Health and Disease

Chemical process systems engineering considers complex supply chains which are coupled networks of dynamically interacting systems. The quest to optimize the supply chain while meeting robustness and flexibility constraints in the face of ever changing environments necessitated the development of theoretical and computational tools for the analysis, synthesis and design of such complex engineered architectures. However, it was realized early on that optimality is a complex characteristic required to achieve proper balance between multiple, often competing, objectives. As we begin to unravel life's intricate complexities, we realize that that living systems share similar structural and dynamic characteristics; hence much can be learned about biological complexity from engineered systems. In this article, we draw analogies between concepts in process systems engineering and conceptual models of health and disease; establish connections between these concepts and physiologic modeling; and describe how these mirror onto the physiological counterparts of engineered systems.

Computational Modeling of Inflammation and Wound Healing

OBJECTIVE

Inflammation is both central to proper wound healing and a key driver of chronic tissue injury via a positive-feedback loop incited by incidental cell damage. We seek to derive actionable insights into the role of inflammation in wound healing in order to improve outcomes for individual patients.

APPROACH

To date, dynamic computational models have been used to study the time evolution of inflammation in wound healing. Emerging clinical data on histo-pathological and macroscopic images of evolving wounds, as well as noninvasive measures of blood flow, suggested the need for tissue-realistic, agent-based, and hybrid mechanistic computational simulations of inflammation and wound healing.

INNOVATION

We developed a computational modeling system, Simple Platform for Agent-based Representation of Knowledge, to facilitate the construction of tissue-realistic models.

RESULTS

A hybrid equation-agent-based model (ABM) of pressure ulcer formation in both spinal cord-injured and -uninjured patients was used to identify control points that reduce stress caused by tissue ischemia/reperfusion. An ABM of arterial restenosis revealed new dynamics of cell migration during neointimal hyperplasia that match histological features, but contradict the currently prevailing mechanistic hypothesis. ABMs of vocal fold inflammation were used to predict inflammatory trajectories in individuals, possibly allowing for personalized treatment.

CONCLUSIONS

The intertwined inflammatory and wound healing responses can be modeled computationally to make predictions in individuals, simulate therapies, and gain mechanistic insights.

A Systems Engineering Perspective on Homeostasis and Disease

Engineered systems are coupled networks of interacting sub-systems, whose dynamics are constrained to requirements of robustness and flexibility. They have evolved by design to optimize function in a changing environment and maintain responses within ranges. Analysis, synthesis, and design of complex supply chains aim to identify and explore the laws governing optimally integrated systems. Optimality expresses balance between conflicting objectives while resiliency results from dynamic interactions among elements. Our increasing understanding of life’s multi-scale architecture suggests that living systems share similar characteristics with much to be learned about biological complexity from engineered systems. If health reflects a dynamically stable integration of molecules, cell, tissues, and organs; disease indicates displacement compensated for and corrected by activation and combination of feedback mechanisms through interconnected networks. In this article, we draw analogies between concepts in systems engineering and conceptual models of health and disease; establish connections between these concepts and physiologic modeling; and describe how these mirror onto the physiological counterparts of engineered systems.

Modeling physiologic variability in human endotoxemia

The control and management of inflammation is a key aspect of clinical care for critical illnesses such as sepsis. In an ideal reaction to injury, the inflammatory response provokes a strong enough response to heal the injury and then restores homeostasis. When inflammation becomes dysregulated, a persistent inflammatory state can lead to significant deleterious effects and clinical challenges. Thus, gaining a better biological understanding of the mechanisms driving the inflammatory response is of the utmost importance. In this review, we discuss our work with the late Stephen F. Lowry to investigate systemic inflammation through systems biology of human endotoxemia. We present our efforts in modeling the human endotoxemia response with a particular focus on physiologic variability. Through modeling, with a focus ultimately on translational applications, we obtain more fundamental understanding of relevant physiological processes. And by taking advantage of the information embedded in biological rhythms, ranging in time scale from high-frequency autonomic oscillations reflected in heart rate variability to circadian rhythms in inflammatory mediators, we gain insight into the underlying physiology.

Linking Inflammation, Cardiorespiratory Variability, and Neural Control in Acute Inflammation via Computational Modeling

Acute inflammation leads to organ failure by engaging catastrophic feedback loops in which stressed tissue evokes an inflammatory response and, in turn, inflammation damages tissue. Manifestations of this maladaptive inflammatory response include cardio-respiratory dysfunction that may be reflected in reduced heart rate and ventilatory pattern variabilities. We have developed signal-processing algorithms that quantify non-linear deterministic characteristics of variability in biologic signals. Now, coalescing under the aegis of the NIH Computational Biology Program and the Society for Complexity in Acute Illness, two research teams performed iterative experiments and computational modeling on inflammation and cardio-pulmonary dysfunction in sepsis as well as on neural control of respiration and ventilatory pattern variability. These teams, with additional collaborators, have recently formed a multi-institutional, interdisciplinary consortium, whose goal is to delineate the fundamental interrelationship between the inflammatory response and physiologic variability. Multi-scale mathematical modeling and complementary physiological experiments will provide insight into autonomic neural mechanisms that may modulate the inflammatory response to sepsis and simultaneously reduce heart rate and ventilatory pattern variabilities associated with sepsis. This approach integrates computational models of neural control of breathing and cardio-respiratory coupling with models that combine inflammation, cardiovascular function, and heart rate variability. The resulting integrated model will provide mechanistic explanations for the phenomena of respiratory sinus-arrhythmia and cardio-ventilatory coupling observed under normal conditions, and the loss of these properties during sepsis. This approach holds the potential of modeling cross-scale physiological interactions to improve both basic knowledge and clinical management of acute inflammatory diseases such as sepsis and trauma.

Sepsis: from pattern to mechanism and back

Sepsis is a clinical entity in which complex inflammatory and physiological processes are mobilized, not only across a range of cellular and molecular interactions, but also in clinically relevant physiological signals accessible at the bedside. There is a need for a mechanistic understanding that links the clinical phenomenon of physiologic variability with the underlying patterns of the biology of inflammation, and we assert that this can be facilitated through the use of dynamic mathematical and computational modeling. An iterative approach of laboratory experimentation and mathematical/computational modeling has the potential to integrate cellular biology, physiology, control theory, and systems engineering across biological scales, yielding insights into the control structures that govern mechanisms by which phenomena, detected as biological patterns, are produced. This approach can represent hypotheses in the formal language of mathematics and computation, and link behaviors that cross scales and domains, thereby offering the opportunity to better explain, diagnose, and intervene in the care of the septic patient.

Multimarker panels in sepsis

Biomarkers differentiate between 2 or more biologic states. The complexity of diseases like sepsis makes it unlikely that any single marker will allow for precise disease specification. Combining several biomarkers into a single classification rule should help to improve their accuracy and, therefore, their usefulness. This article reviews several studies using multimarker panels, and highlights the potential of more sophisticated diagnostic and prognostic techniques in future multimarker panels. More complex algorithms should accelerate the adoption of multimarker panels into the routine management of patients with sepsis, provided that clinicians understand the multimarker approach.

Modeling autonomic regulation of cardiac function and heart rate variability in human endotoxemia

Heart rate variability (HRV), the quantification of beat-to-beat variability, has been studied as a potential prognostic marker in inflammatory diseases such as sepsis. HRV normally reflects significant levels of variability in homeostasis, which can be lost under stress. Much effort has been placed in interpreting HRV from the perspective of quantitatively understanding how stressors alter HRV dynamics, but the molecular and cellular mechanisms that give rise to both homeostatic HRV and changes in HRV have received less focus. Here, we develop a mathematical model of human endotoxemia that incorporates the oscillatory signals giving rise to HRV and their signal transduction to the heart. Connections between processes at the cellular, molecular, and neural levels are quantitatively linked to HRV. Rhythmic signals representing autonomic oscillations and circadian rhythms converge to modulate the pattern of heartbeats, and the effects of these oscillators are diminished in the acute endotoxemia response. Based on the semimechanistic model developed herein, homeostatic and acute stress responses of HRV are studied in terms of these oscillatory signals. Understanding the loss of HRV in endotoxemia serves as a step toward understanding changes in HRV observed clinically through translational applications of systems biology based on the relationship between biological processes and clinical outcomes.

Multiscale model for the assessment of autonomic dysfunction in human endotoxemia

Severe injury and infection are associated with autonomic dysfunction. The realization that a dysregulation in autonomic function may predispose a host to excessive inflammatory processes has renewed interest in understanding the role of central nervous system (CNS) in modulating systemic inflammatory processes. Assessment of heart rate variability (HRV) has been used to evaluate systemic abnormalities and as a predictor of the severity of illness. Dissecting the relevance of neuroimmunomodulation in controlling inflammatory processes requires an understanding of the multiscale interplay between CNS and the immune response. A vital enabler in that respect is the development of a systems-based approach that integrates data across multiple scales, and models the emerging host response as the outcome of interactions of critical modules. Thus, a multiscale model of human endotoxemia, as a prototype model of systemic inflammation in humans, is proposed that integrates processes across the host from the cellular to the systemic host response level. At the cellular level interacting components are associated with elementary signaling pathways that propagate extracellular signals to the transcriptional response level. Further, essential modules associated with the neuroendocrine immune crosstalk are considered. Finally, at the systemic level, phenotypic expressions such as HRV are incorporated to assess systemic decomplexification indicative of the severity of the host response. Thus, the proposed work intends to associate acquired endocrine dysfunction with diminished HRV as a critical enabler for clarifying how cellular inflammatory processes and neural-based pathways mediate the links between patterns of autonomic control (HRV) and clinical outcomes.

Agent-based modeling of endotoxin-induced acute inflammatory response in human blood leukocytes

BACKGROUND

Inflammation is a highly complex biological response evoked by many stimuli. A persistent challenge in modeling this dynamic process has been the (nonlinear) nature of the response that precludes the single-variable assumption. Systems-based approaches offer a promising possibility for understanding inflammation in its homeostatic context. In order to study the underlying complexity of the acute inflammatory response, an agent-based framework is developed that models the emerging host response as the outcome of orchestrated interactions associated with intricate signaling cascades and intercellular immune system interactions.

METHODOLOGY/PRINCIPAL FINDINGS

An agent-based modeling (ABM) framework is proposed to study the nonlinear dynamics of acute human inflammation. The model is implemented using NetLogo software. Interacting agents involve either inflammation-specific molecules or cells essential for the propagation of the inflammatory reaction across the system. Spatial orientation of molecule interactions involved in signaling cascades coupled with the cellular heterogeneity are further taken into account. The proposed in silico model is evaluated through its ability to successfully reproduce a self-limited inflammatory response as well as a series of scenarios indicative of the nonlinear dynamics of the response. Such scenarios involve either a persistent (non)infectious response or innate immune tolerance and potentiation effects followed by perturbations in intracellular signaling molecules and cascades.

CONCLUSIONS/SIGNIFICANCE

The ABM framework developed in this study provides insight on the stochastic interactions of the mediators involved in the propagation of endotoxin signaling at the cellular response level. The simulation results are in accordance with our prior research effort associated with the development of deterministic human inflammation models that include transcriptional dynamics, signaling, and physiological components. The hypothetical scenarios explored in this study would potentially improve our understanding of how manipulating the behavior of the molecular species could manifest into emergent behavior of the overall system.

NETWORKS, BIOLOGY AND SYSTEMS ENGINEERING: A CASE STUDY IN INFLAMMATION

Biological systems can be modeled as networks of interacting components across multiple scales. A central problem in computational systems biology is to identify those critical components and the rules that define their interactions and give rise to the emergent behavior of a host response. In this paper we will discuss two fundamental problems related to the construction of transcription factor networks and the identification of networks of functional modules describing disease progression. We focus on inflammation as a key physiological response of clinical and translational importance.

In silico simulation of corticosteroids effect on an NFkB- dependent physicochemical model of systemic inflammation

BACKGROUND

During the onset of an inflammatory response signaling pathways are activated for "translating" extracellular signals into intracellular responses converging to the activation of nuclear factor (NF)-kB, a central transcription factor in driving the inflammatory response. An inadequate control of its transcriptional activity is associated with the culmination of a hyper-inflammatory response making it a desired therapeutic target. Predicated upon the nature of the response, a systems level analysis might provide rational leads for the development of strategies that promote the resolution of the response.

METHODOLOGY AND FINDINGS

A physicochemical host response model is proposed to integrate biological information in the form of kinetic rules and signaling cascades with pharmacokinetic models of drug action for the modulation of the response. The unifying hypothesis is that the response is triggered by the activation of the NFkB signaling module and corticosteroids serve as a template for assessing anti-inflammatory strategies. The proposed in silico model is evaluated through its ability to predict and modulate uncontrolled responses. The pre-exposure of the system to hypercortisolemia, i.e. 6 hr before or simultaneously with the infectious challenge "reprograms" the dynamics of the host towards a balanced inflammatory response. However, if such an intervention occurs long before the inflammatory insult a symptomatic effect is observed instead of a protective relief while a steroid infusion after inducing inflammation requires much higher drug doses.

CONCLUSIONS AND SIGNIFICANCE

We propose a reversed engineered inflammation model that seeks to describe how the system responds to a multitude of external signals. Timing of intervention and dosage regimes appears to be key determinants for the protective or symptomatic effect of exogenous corticosteroids. Such results lie in qualitative agreement with in vivo human studies exposed both to LPS and corticosteroids under various time intervals thus improving our understanding of how interacting modules generate a behavior.