Translatability and transferability of in silico models: context of use switching to predict the effects of environmental chemicals on the immune system

Immunotoxicity hazard identification of chemicals aims to evaluate the potential for unintended effects of chemical exposure on the immune system. Perfluorinated alkylate substances (PFAS), such as perfluorooctane sulfonate (PFOS) and perfluorooctanoic acid (PFOA), are persistent, globally disseminated environmental contaminants known to be immunotoxic. Elevated PFAS exposure is associated with lower antibody responses to vaccinations in children and in adults. In addition, some studies have reported a correlation between PFAS levels in the body and lower resistance to disease, in other words an increased risk of infections or cancers.

In this context, modelling and simulation platforms could be used to simulate the human immune system with the aim to evaluate the adverse effects that immunotoxicants may have.

Here, we show the conditions under which a mathematical model developed for one purpose and application (e.g., in the pharmaceutical domain) can be successfully translated and transferred to another (e.g., in the chemicals domain) without undergoing significant adaptation. In particular, we demonstrate that the Universal Immune System Simulator was able to simulate the effects of PFAS on the immune system, introducing entities and new interactions that are biologically involved in the phenomenon. This also revealed a potentially exploitable pathway for assessing immunotoxicity through a computational model.

CaliPro: A Calibration Protocol That Utilizes Parameter Density Estimation to Explore Parameter Space and Calibrate Complex Biological Models

Introduction

Mathematical and computational modeling have a long history of uncovering mechanisms and making predictions for biological systems. However, to create a model that can provide relevant quantitative predictions, models must first be calibrated by recapitulating existing biological datasets from that system. Current calibration approaches may not be appropriate for complex biological models because: 1) many attempt to recapitulate only a single aspect of the experimental data (such as a median trend) or 2) Bayesian techniques require specification of parameter priors and likelihoods to experimental data that cannot always be confidently assigned. A new calibration protocol is needed to calibrate complex models when current approaches fall short.

Methods

Herein, we develop CaliPro, an iterative, model-agnostic calibration protocol that utilizes parameter density estimation to refine parameter space and calibrate to temporal biological datasets. An important aspect of CaliPro is the user-defined pass set definition, which specifies how the model might successfully recapitulate experimental data. We define the appropriate settings to use CaliPro.

Results

We illustrate the usefulness of CaliPro through four examples including predator-prey, infectious disease transmission, and immune response models. We show that CaliPro works well for both deterministic, continuous model structures as well as stochastic, discrete models and illustrate that CaliPro can work across diverse calibration goals.

Conclusions

We present CaliPro, a new method for calibrating complex biological models to a range of experimental outcomes. In addition to expediting calibration, CaliPro may be useful in already calibrated parameter spaces to target and isolate specific model behavior for further analysis.

Computer Modeling of Clonal Dominance: Memory-Anti-Naïve and Its Curbing by Attrition

Experimental and computational studies have revealed that T-cell cross-reactivity is a widespread phenomenon that can either be advantageous or detrimental to the host. In particular, detrimental effects can occur whenever the clonal dominance of memory cells is not justified by their infection-clearing capacity. Using an agent-based model of the immune system, we recently predicted the "memory anti-naïve" phenomenon, which occurs when the secondary challenge is similar but not identical to the primary stimulation. In this case, the pre-existing memory cells formed during the primary infection may be rapidly deployed in spite of their low affinity and can actually prevent a potentially higher affinity naïve response from emerging, resulting in impaired viral clearance. This finding allowed us to propose a mechanistic explanation for the concept of "antigenic sin" originally described in the context of the humoral response. However, the fact that antigenic sin is a relatively rare occurrence suggests the existence of evolutionary mechanisms that can mitigate the effect of the memory anti-naïve phenomenon. In this study we use computer modeling to further elucidate clonal dominance and the memory anti-naïve phenomenon, and to investigate a possible mitigating factor called attrition. Attrition has been described in the experimental and computational literature as a combination of competition for space and apoptosis of lymphocytes via type-I interferon in the early stages of a viral infection. This study systematically explores the relationship between clonal dominance and the mechanism of attrition. Our results suggest that attrition can indeed mitigate the memory anti-naïve effect by enabling the emergence of a diverse, higher affinity naïve response against the secondary challenge. In conclusion, modeling attrition allows us to shed light on the nature of clonal interaction and dominance.

Modeling genotypes in their microenvironment to predict single- and multi-cellular behavior

A cell's phenotype is the set of observable characteristics resulting from the interaction of the genotype with the surrounding environment, determining cell behavior. Deciphering genotype-phenotype relationships has been crucial to understanding normal and disease biology. Analysis of molecular pathways has provided an invaluable tool to such understanding; however, typically it does not consider the physical microenvironment, which is a key determinant of phenotype. In this study, we present a novel modeling framework that enables the study of the link between genotype, signaling networks, and cell behavior in a three-dimensional microenvironment. To achieve this, we bring together Agent-Based Modeling, a powerful computational modeling technique, and gene networks. This combination allows biological hypotheses to be tested in a controlled stepwise fashion, and it lends itself naturally to model a heterogeneous population of cells acting and evolving in a dynamic microenvironment, which is needed to predict the evolution of complex multi-cellular dynamics. Importantly, this enables modeling co-occurring intrinsic perturbations, such as mutations, and extrinsic perturbations, such as nutrient availability, and their interactions. Using cancer as a model system, we illustrate how this framework delivers a unique opportunity to identify determinants of single-cell behavior, while uncovering emerging properties of multi-cellular growth. This framework is freely available at http://www.microc.org.

Modeling the Effect of High Calorie Diet on the Interplay between Adipose Tissue, Inflammation, and Diabetes

Background

Type 2 diabetes (T2D) is a chronic metabolic disease potentially leading to serious widespread tissue damage. Human organism develops T2D when the glucose-insulin control is broken for reasons that are not fully understood but have been demonstrated to be linked to the emergence of a chronic inflammation. Indeed such low-level chronic inflammation affects the pancreatic production of insulin and triggers the development of insulin resistance, eventually leading to an impaired control of the blood glucose concentration. On the contrary, it is well-known that obesity and inflammation are strongly correlated.

Aim

In this study, we investigate in silico the effect of overfeeding on the adipose tissue and the consequent set up of an inflammatory state. We model the emergence of the inflammation as the result of adipose mass increase which, in turn, is a direct consequence of a prolonged excess of high calorie intake.

Results

The model reproduces the fat accumulation due to excessive caloric intake observed in two clinical studies. Moreover, while showing consistent weight gains over long periods of time, it reveals a drift of the macrophage population toward the proinflammatory phenotype, thus confirming its association with fatness.

Multiscale Agent-Based and Hybrid Modeling of the Tumor Immune Microenvironment

Multiscale systems biology and systems pharmacology are powerful methodologies that are playing increasingly important roles in understanding the fundamental mechanisms of biological phenomena and in clinical applications. In this review, we summarize the state of the art in the applications of agent-based models (ABM) and hybrid modeling to the tumor immune microenvironment and cancer immune response, including immunotherapy. Heterogeneity is a hallmark of cancer; tumor heterogeneity at the molecular, cellular, and tissue scales is a major determinant of metastasis, drug resistance, and low response rate to molecular targeted therapies and immunotherapies. Agent-based modeling is an effective methodology to obtain and understand quantitative characteristics of these processes and to propose clinical solutions aimed at overcoming the current obstacles in cancer treatment. We review models focusing on intra-tumor heterogeneity, particularly on interactions between cancer cells and stromal cells, including immune cells, the role of tumor-associated vasculature in the immune response, immune-related tumor mechanobiology, and cancer immunotherapy. We discuss the role of digital pathology in parameterizing and validating spatial computational models and potential applications to therapeutics.

Gene Regulatory Network Modeling of Macrophage Differentiation Corroborates the Continuum Hypothesis of Polarization States

Macrophages derived from monocyte precursors undergo specific polarization processes which are influenced by the local tissue environment: classically activated (M1) macrophages, with a pro-inflammatory activity and a role of effector cells in Th1 cellular immune responses, and alternatively activated (M2) macrophages, with anti-inflammatory functions and involved in immunosuppression and tissue repair. At least three different subsets of M2 macrophages, namely, M2a, M2b, and M2c, are characterized in the literature based on their eliciting signals. The activation and polarization of macrophages is achieved through many, often intertwined, signaling pathways. To describe the logical relationships among the genes involved in macrophage polarization, we used a computational modeling methodology, namely, logical (Boolean) modeling of gene regulation. We integrated experimental data and knowledge available in the literature to construct a logical network model for the gene regulation driving macrophage polarization to the M1, M2a, M2b, and M2c phenotypes. Using the software GINsim and BoolNet, we analyzed the network dynamics under different conditions and perturbations to understand how they affect cell polarization. Dynamic simulations of the network model, enacting the most relevant biological conditions, showed coherence with the observed behavior of in vivo macrophages. The model could correctly reproduce the polarization toward the four main phenotypes as well as to several hybrid phenotypes, which are known to be experimentally associated to physiological and pathological conditions. We surmise that shifts among different phenotypes in the model mimic the hypothetical continuum of macrophage polarization, with M1 and M2 being the extremes of an uninterrupted sequence of states. Furthermore, model simulations suggest that anti-inflammatory macrophages are resilient to shift back to the pro-inflammatory phenotype.

Facing the challenges of multiscale modelling of bacterial and fungal pathogen-host interactions

Recent and rapidly evolving progress on high-throughput measurement techniques and computational performance has led to the emergence of new disciplines, such as systems medicine and translational systems biology. At the core of these disciplines lies the desire to produce multiscale models: mathematical models that integrate multiple scales of biological organization, ranging from molecular, cellular and tissue models to organ, whole-organism and population scale models. Using such models, hypotheses can systematically be tested. In this review, we present state-of-the-art multiscale modelling of bacterial and fungal infections, considering both the pathogen and host as well as their interaction. Multiscale modelling of the interactions of bacteria, especially Mycobacterium tuberculosis, with the human host is quite advanced. In contrast, models for fungal infections are still in their infancy, in particular regarding infections with the most important human pathogenic fungi, Candida albicans and Aspergillus fumigatus. We reflect on the current availability of computational approaches for multiscale modelling of host-pathogen interactions and point out current challenges. Finally, we provide an outlook for future requirements of multiscale modelling.

IL-2 sensitivity and exogenous IL-2 concentration gradient tune the productive contact duration of CD8(+) T cell-APC: a multiscale modeling study

Background

The CD8⁺ T cell immune response fights acute infections by intracellular pathogens and, by generating an immune memory, enables immune responses against secondary infections. Activation of the CD8⁺ T cell immune response involves a succession of molecular events leading to modifications of CD8⁺ T cell population. To understand the endogenous and exogenous mechanisms controlling the activation of CD8⁺ T cells and to investigate the influence of early molecular events on the long-term cell population behavior, we developed a multiscale computational model. It integrates three levels of description: a Cellular Potts model describing the individual behavior of CD8⁺ T cells, a system of ordinary differential equations describing a decision-making molecular regulatory network at the intracellular level, and a partial differential equation describing the diffusion of IL-2 in the extracellular environment.

Results

We first calibrated the model parameters based on in vivo data and showed the model’s ability to reproduce early dynamics of CD8⁺ T cells in murine lymph nodes after influenza infection, both at the cell population and intracellular levels. We then showed the model’s ability to reproduce the proliferative responses of CD5^hi and CD5^lo CD8⁺ T cells to exogenous IL-2 under a weak TCR stimulation. This stressed the role of short-lasting molecular events and the relevance of explicitly describing both intracellular and cellular scale dynamics. Our results suggest that the productive contact duration of CD8⁺ T cell-APC is influenced by the sensitivity of individual CD8⁺ T cells to the activation signal and by the IL-2 concentration in the extracellular environment.

Conclusions

The multiscale nature of our model allows the reproduction and explanation of some acquired characteristics and functions of CD8⁺ T cells, and of their responses to multiple stimulation conditions, that would not be accessible in a classical description of cell population dynamics that would not consider intracellular dynamics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0323-y) contains supplementary material, which is available to authorized users.

Introduction of a Framework for Dynamic Knowledge Representation of the Control Structure of Transplant Immunology: Employing the Power of Abstraction with a Solid Organ Transplant Agent-Based Model

Agent-based modeling has been used to characterize the nested control loops and non-linear dynamics associated with inflammatory and immune responses, particularly as a means of visualizing putative mechanistic hypotheses. This process is termed dynamic knowledge representation and serves a critical role in facilitating the ability to test and potentially falsify hypotheses in the current data- and hypothesis-rich biomedical research environment. Importantly, dynamic computational modeling aids in identifying useful abstractions, a fundamental scientific principle that pervades the physical sciences. Recognizing the critical scientific role of abstraction provides an intellectual and methodological counterweight to the tendency in biology to emphasize comprehensive description as the primary manifestation of biological knowledge. Transplant immunology represents yet another example of the challenge of identifying sufficient understanding of the inflammatory/immune response in order to develop and refine clinically effective interventions. Advances in immunosuppressive therapies have greatly improved solid organ transplant (SOT) outcomes, most notably by reducing and treating acute rejection. The end goal of these transplant immune strategies is to facilitate effective control of the balance between regulatory T cells and the effector/cytotoxic T-cell populations in order to generate, and ideally maintain, a tolerant phenotype. Characterizing the dynamics of immune cell populations and the interactive feedback loops that lead to graft rejection or tolerance is extremely challenging, but is necessary if rational modulation to induce transplant tolerance is to be accomplished. Herein is presented the solid organ agent-based model (SOTABM) as an initial example of an agent-based model (ABM) that abstractly reproduces the cellular and molecular components of the immune response to SOT. Despite its abstract nature, the SOTABM is able to qualitatively reproduce acute rejection and the suppression of acute rejection by immunosuppression to generate transplant tolerance. The SOTABM is intended as an initial example of how ABMs can be used to dynamically represent mechanistic knowledge concerning transplant immunology in a scalable and expandable form and can thus potentially serve as useful adjuncts to the investigation and development of control strategies to induce transplant tolerance.

Multiscale modelling in immunology: a review

One of the greatest challenges in biomedicine is to get a unified view of observations made from the molecular up to the organism scale. Towards this goal, multiscale models have been highly instrumental in contexts such as the cardiovascular field, angiogenesis, neurosciences and tumour biology. More recently, such models are becoming an increasingly important resource to address immunological questions as well. Systematic mining of the literature in multiscale modelling led us to identify three main fields of immunological applications: host-virus interactions, inflammatory diseases and their treatment and development of multiscale simulation platforms for immunological research and for educational purposes. Here, we review the current developments in these directions, which illustrate that multiscale models can consistently integrate immunological data generated at several scales, and can be used to describe and optimize therapeutic treatments of complex immune diseases.

Strategies for efficient numerical implementation of hybrid multi-scale agent-based models to describe biological systems

Biologically related processes operate across multiple spatiotemporal scales. For computational modeling methodologies to mimic this biological complexity, individual scale models must be linked in ways that allow for dynamic exchange of information across scales. A powerful methodology is to combine a discrete modeling approach, agent-based models (ABMs), with continuum models to form hybrid models. Hybrid multi-scale ABMs have been used to simulate emergent responses of biological systems. Here, we review two aspects of hybrid multi-scale ABMs: linking individual scale models and efficiently solving the resulting model. We discuss the computational choices associated with aspects of linking individual scale models while simultaneously maintaining model tractability. We demonstrate implementations of existing numerical methods in the context of hybrid multi-scale ABMs. Using an example model describing Mycobacterium tuberculosis infection, we show relative computational speeds of various combinations of numerical methods. Efficient linking and solution of hybrid multi-scale ABMs is key to model portability, modularity, and their use in understanding biological phenomena at a systems level.

Agent-based model of fecal microbial transplant effect on bile acid metabolism on suppressing Clostridium difficile infection: an example of agent-based modeling of intestinal bacterial infection

Agent-based modeling is a computational modeling method that represents system-level behavior as arising from multiple interactions between the multiple components that make up a system. Biological systems are thus readily described using agent-based models (ABMs), as multi-cellular organisms can be viewed as populations of interacting cells, and microbial systems manifest as colonies of individual microbes. Intersections between these two domains underlie an increasing number of pathophysiological processes, and the intestinal tract represents one of the most significant locations for these inter-domain interactions, so much so that it can be considered an internal ecology of varying robustness and function. Intestinal infections represent significant disturbances of this internal ecology, and one of the most clinically relevant intestinal infections is Clostridium difficile infection (CDI). CDI is precipitated by the use of broad-spectrum antibiotics, involves the depletion of commensal microbiota, and alterations in bile acid composition in the intestinal lumen. We present an example ABM of CDI (the C. difficile Infection ABM, or CDIABM) to examine fundamental dynamics of the pathogenesis of CDI and its response to treatment with anti-CDI antibiotics and a newer treatment therapy, fecal microbial transplant. The CDIABM focuses on one specific mechanism of potential CDI suppression: commensal modulation of bile acid composition. Even given its abstraction, the CDIABM reproduces essential dynamics of CDI and its response to therapy, and identifies a paradoxical zone of behavior that provides insight into the role of intestinal nutritional status and the efficacy of anti-CDI therapies. It is hoped that this use case example of the CDIABM can demonstrate the usefulness of both agent-based modeling and the application of abstract functional representation as the biomedical community seeks to address the challenges of increasingly complex diseases with the goal of personalized medicine.

Computational modeling of heterogeneity and function of CD4+ T cells

The immune system is composed of many different cell types and hundreds of intersecting molecular pathways and signals. This large biological complexity requires coordination between distinct pro-inflammatory and regulatory cell subsets to respond to infection while maintaining tissue homeostasis. CD4+ T cells play a central role in orchestrating immune responses and in maintaining a balance between pro- and anti- inflammatory responses. This tight balance between regulatory and effector reactions depends on the ability of CD4+ T cells to modulate distinct pathways within large molecular networks, since dysregulated CD4+ T cell responses may result in chronic inflammatory and autoimmune diseases. The CD4+ T cell differentiation process comprises an intricate interplay between cytokines, their receptors, adaptor molecules, signaling cascades and transcription factors that help delineate cell fate and function. Computational modeling can help to describe, simulate, analyze, and predict some of the behaviors in this complicated differentiation network. This review provides a comprehensive overview of existing computational immunology methods as well as novel strategies used to model immune responses with a particular focus on CD4+ T cell differentiation.

Modeling biology spanning different scales: an open challenge

It is coming nowadays more clear that in order to obtain a unified description of the different mechanisms governing the behavior and causality relations among the various parts of a living system, the development of comprehensive computational and mathematical models at different space and time scales is required. This is one of the most formidable challenges of modern biology characterized by the availability of huge amount of high throughput measurements. In this paper we draw attention to the importance of multiscale modeling in the framework of studies of biological systems in general and of the immune system in particular.

Identification of core T cell network based on immunome interactome

Background

Data-driven studies on the dynamics of reconstructed protein-protein interaction (PPI) networks facilitate investigation and identification of proteins important for particular processes or diseases and reduces time and costs of experimental verification. Modeling the dynamics of very large PPI networks is computationally costly.

Results

To circumvent this problem, we created a link-weighted human immunome interactome and performed filtering. We reconstructed the immunome interactome and weighed the links using jackknife gene expression correlation of integrated, time course gene expression data. Statistical significance of the links was computed using the Global Statistical Significance (GloSS) filtering algorithm. P-values from GloSS were computed for the integrated, time course gene expression data. We filtered the immunome interactome to identify core components of the T cell PPI network (TPPIN). The interconnectedness of the major pathways for T cell survival and response, including the T cell receptor, MAPK and JAK-STAT pathways, are maintained in the TPPIN network. The obtained TPPIN network is supported both by Gene Ontology term enrichment analysis along with study of essential genes enrichment.

Conclusions

By integrating gene expression data to the immunome interactome and using a weighted network filtering method, we identified the T cell PPI immune response network. This network reveals the most central and crucial network in T cells. The approach is general and applicable to any dataset that contains sufficient information.

The role of low-grade inflammation and metabolic flexibility in aging and nutritional modulation thereof: a systems biology approach

Aging is a biological process characterized by the progressive functional decline of many interrelated physiological systems. In particular, aging is associated with the development of a systemic state of low-grade chronic inflammation (inflammaging), and with progressive deterioration of metabolic function. Systems biology has helped in identifying the mediators and pathways involved in these phenomena, mainly through the application of high-throughput screening methods, valued for their molecular comprehensiveness. Nevertheless, inflammation and metabolic regulation are dynamical processes whose behavior must be understood at multiple levels of biological organization (molecular, cellular, organ, and system levels) and on multiple time scales. Mathematical modeling of such behavior, with incorporation of mechanistic knowledge on interactions between inflammatory and metabolic mediators, may help in devising nutritional interventions capable of preventing, or ameliorating, the age-associated functional decline of the corresponding systems.

Competition for antigen between Th1 and Th2 responses determines the timing of the immune response switch during Mycobaterium avium subspecies paratuberulosis infection in ruminants

Johne's disease (JD), a persistent and slow progressing infection of ruminants such as cows and sheep, is caused by slow replicating bacilli Mycobacterium avium subspecies paratuberculosis (MAP) infecting macrophages in the gut. Infected animals initially mount a cell-mediated CD4 T cell response against MAP which is characterized by the production of interferon (Th1 response). Over time, Th1 response diminishes in most animals and antibody response to MAP antigens becomes dominant (Th2 response). The switch from Th1 to Th2 response occurs concomitantly with disease progression and shedding of the bacteria in feces. Mechanisms controlling this Th1/Th2 switch remain poorly understood. Because Th1 and Th2 responses are known to cross-inhibit each other, it is unclear why initially strong Th1 response is lost over time. Using a novel mathematical model of the immune response to MAP infection we show that the ability of extracellular bacteria to persist outside of macrophages naturally leads to switch of the cellular response to antibody production. Several additional mechanisms may also contribute to the timing of the Th1/Th2 switch including the rate of proliferation of Th1/Th2 responses at the site of infection, efficiency at which immune responses cross-inhibit each other, and the rate at which Th1 response becomes exhausted over time. Our basic model reasonably well explains four different kinetic patterns of the Th1/Th2 responses in MAP-infected sheep by variability in the initial bacterial dose and the efficiency of the MAP-specific T cell responses. Taken together, our novel mathematical model identifies factors of bacterial and host origin that drive kinetics of the immune response to MAP and provides the basis for testing the impact of vaccination or early treatment on the duration of infection.

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease, a chronic enteric disease of ruminants such as sheep and cows. Due to early culling and reduction in milk production of affected animals, MAP inflicts high economic cost to diary farms. MAP infection has a long incubation period of several years, and during the asymptomatic stage a strong cellular (T helper 1) immune response is thought to control MAP replication. Over time, Th1 response is lost and ineffective antibody response driven by Th2 cells becomes predominant. We develop the first mathematical model of helper T cell response to MAP infection to understand impact of various mechanisms on the dynamics of the switch from Th1 to Th2 response. Our results suggest that in contrast to the generally held belief, Th1/Th2 switch may be driven by the accumulation of long-lived extracellular bacteria, and therefore, may be the consequence of the disease progression of MAP-infected animals and not its cause. Our model highlights limitations of our current understanding of regulation of helper T cell responses during MAP infection and identifies areas for future experimental research.

Second order optimization for the inference of gene regulatory pathways

With the increasing availability of experimental data on gene interactions, modeling of gene regulatory pathways has gained special attention. Gradient descent algorithms have been widely used for regression and classification applications. Unfortunately, results obtained after training a model by gradient descent are often highly variable. In this paper, we present a new second order learning rule based on the Newton's method for inferring optimal gene regulatory pathways. Unlike the gradient descent method, the proposed optimization rule is independent of the learning parameter. The flow vectors are estimated based on biomass conservation. A set of constraints is formulated incorporating weighting coefficients. The method calculates the maximal expression of the target gene starting from a given initial gene through these weighting coefficients. Our algorithm has been benchmarked and validated on certain types of functions and on some gene regulatory networks, gathered from literature. The proposed method has been found to perform better than the gradient descent learning. Extensive performance comparison with the extreme pathway analysis method has underlined the effectiveness of our proposed methodology.

The onset of type 2 diabetes: proposal for a multi-scale model

Background

Type 2 diabetes mellitus (T2D) is a common age-related disease, and is a major health concern, particularly in developed countries where the population is aging, including Europe. The multi-scale immune system simulator for the onset of type 2 diabetes (MISSION-T2D) is a European Union-funded project that aims to develop and validate an integrated, multilevel, and patient-specific model, incorporating genetic, metabolic, and nutritional data for the simulation and prediction of metabolic and inflammatory processes in the onset and progression of T2D. The project will ultimately provide a tool for diagnosis and clinical decision making that can estimate the risk of developing T2D and predict its progression in response to possible therapies. Recent data showed that T2D and its complications, specifically in the heart, kidney, retina, and feet, should be considered a systemic disease that is sustained by a pervasive, metabolically-driven state of inflammation. Accordingly, there is an urgent need (1) to understand the complex mechanisms underpinning the onset of this disease, and (2) to identify early patient-specific diagnostic parameters and related inflammatory indicators.

Objective

We aim to accomplish this mission by setting up a multi-scale model to study the systemic interactions of the biological mechanisms involved in response to a variety of nutritional and metabolic stimuli and stressors.

Methods

Specifically, we will be studying the biological mechanisms of immunological/inflammatory processes, energy intake/expenditure ratio, and cell cycle rate. The overall architecture of the model will exploit an already established immune system simulator as well as several discrete and continuous mathematical methods for modeling of the processes critically involved in the onset and progression of T2D. We aim to validate the predictions of our models using actual biological and clinical data.

Results

This study was initiated in March 2013 and is expected to be completed by February 2016.

Conclusions

MISSION-T2D aims to pave the way for translating validated multilevel immune-metabolic models into the clinical setting of T2D. This approach will eventually generate predictive biomarkers for this disease from the integration of clinical data with metabolic, nutritional, immune/inflammatory, genetic, and gut microbiota profiles. Eventually, it should prove possible to translate these into cost-effective and mobile-based diagnostic tools.

Cellular and population plasticity of helper CD4(+) T cell responses

Vertebrates are constantly exposed to pathogens, and the adaptive immunity has most likely evolved to control and clear such infectious agents. CD4⁺ T cells are the major players in the adaptive immune response to pathogens. Following recognition of pathogen-derived antigens naïve CD4⁺ T cells differentiate into effectors which then control pathogen replication either directly by killing pathogen-infected cells or by assisting with generation of cytotoxic T lymphocytes (CTLs) or pathogen-specific antibodies. Pathogen-specific effector CD4⁺ T cells are highly heterogeneous in terms of cytokines they produce. Three major subtypes of effector CD4⁺ T cells have been identified: T-helper 1 (Th1) cells producing IFN-γ and TNF-α, Th2 cells producing IL-4 and IL-10, and Th17 cells producing IL-17. How this heterogeneity is maintained and what regulates changes in effector T cell composition during chronic infections remains poorly understood. In this review we discuss recent advances in our understanding of CD4⁺ T cell differentiation in response to microbial infections. We propose that a change in the phenotype of pathogen-specific effector CD4⁺ T cells during chronic infections, for example, from Th1 to Th2 response as observed in Mycobactrium avium ssp. paratuberculosis (MAP) infection of ruminants, can be achieved by conversion of T cells from one effector subset to another (cellular plasticity) or due to differences in kinetics (differentiation, proliferation, death) of different effector T cell subsets (population plasticity). We also shortly review mathematical models aimed at describing CD4⁺ T cell differentiation and outline areas for future experimental and theoretical research.

Interleukin 12 receptor deficiency in a child with recurrent bronchopneumonia and very high IgE levels

Interleukin-12 (IL-12) is involved in cellular immune responses against intracellular pathogens by promoting the generation of T naive in T helper 1 (Th1) cells and by increasing interferon-gamma (IFN-gamma) production from T and natural killer (NK) cells. A defective induction of a Th1 response may lead to a higher risk of infections, and, in particular, infections due to typical and atypical Mycobacteria. We report on the case of a girl with suffering from recurrent bronchopneumonia associated with very high serum IgE levels, who exhibited a profound impairment of the Th1 generation associated with a novel mutation in the exon 5 of the IL-12R β1 gene (R156H). Our data suggest that in children with severe and recurrent infections, even in the absence of a mycobacterial infection, functional and/or genetic alterations of the molecular mechanisms governing Th1/Th2 homeostasis might be responsible for an atypical immunodeficiency and, therefore, should be investigated in these patients.

How the interval between prime and boost injection affects the immune response in a computational model of the immune system

The immune system is able to respond more vigorously to the second contact with a given antigen than to the first contact. Vaccination protocols generally include at least two doses, in order to obtain high antibody titers. We want to analyze the relation between the time elapsed from the first dose (priming) and the second dose (boost) on the antibody titers. In this paper, we couple in vivo experiments with computer simulations to assess the effect of delaying the second injection. We observe that an interval of several weeks between the prime and the boost is necessary to obtain optimal antibody responses.

Intraindividual long term stability and response corridors of cytokines in healthy volunteers detected by a standardized whole-blood culture system for bed-side application

Background

The variation of immune cell activities over time is an immanent property of the human immune system, as can be measured by the stimulated secretion of cytokines in cell cultures. However, inter-individual variability is considerably higher. Especially the latter is the major reason why it has not been possible to establish international standard values for cytokines as was possible for other parameters, such as leukocyte sub-population numbers. In this trial, a highly standardized whole-blood culture model (TrueCulture®), developed to characterise drug effects on cells of the human immune system in clinical trials, was used to analyse cytokine patterns in the blood samples of 12 healthy subjects over a period of one month.

Methods

After an overnight fast, 12 healthy subjects donated blood three times a week on three consecutive days over a period of 4 weeks. TruCulture® blood collection and whole-blood culture systems were used to measure whole-blood leukocyte stimulation. The levels of IL-2, IL-5, IL-13, IL-6, IL-8, IL-10, IFNγ, and MCP-1 in the culture supernatants were quantified by sandwich ELISA.

Results

The pattern of cytokine concentrations in the supernatants of the stimulated whole-blood cultures was highly individual, but considerably stable over the whole observation period of 4 weeks.

Conclusions

By using TruCulture® it seems feasible to determine subject-specific cytokine reference patterns, for example under healthy conditions, or before starting an experimental treatment, e.g. during a clinical trial, against which changes in the behaviour of the immune system can be detected more accurately in future.

Development of immune-specific interaction potentials and their application in the multi-agent-system VaccImm

Peptide vaccination in cancer therapy is a promising alternative to conventional methods. However, the parameters for this personalized treatment are difficult to access experimentally. In this respect, in silico models can help to narrow down the parameter space or to explain certain phenomena at a systems level. Herein, we develop two empirical interaction potentials specific to B-cell and T-cell receptor complexes and validate their applicability in comparison to a more general potential. The interaction potentials are applied to the model VaccImm which simulates the immune response against solid tumors under peptide vaccination therapy. This multi-agent system is derived from another immune system simulator (C-ImmSim) and now includes a module that enables the amino acid sequence of immune receptors and their ligands to be taken into account. The multi-agent approach is combined with approved methods for prediction of major histocompatibility complex (MHC)-binding peptides and the newly developed interaction potentials. In the analysis, we critically assess the impact of the different modules on the simulation with VaccImm and how they influence each other. In addition, we explore the reasons for failures in inducing an immune response by examining the activation states of the immune cell populations in detail.In summary, the present work introduces immune-specific interaction potentials and their application to the agent-based model VaccImm which simulates peptide vaccination in cancer therapy.

Combining network modeling and gene expression microarray analysis to explore the dynamics of Th1 and Th2 cell regulation

Two T helper (Th) cell subsets, namely Th1 and Th2 cells, play an important role in inflammatory diseases. The two subsets are thought to counter-regulate each other, and alterations in their balance result in different diseases. This paradigm has been challenged by recent clinical and experimental data. Because of the large number of genes involved in regulating Th1 and Th2 cells, assessment of this paradigm by modeling or experiments is difficult. Novel algorithms based on formal methods now permit the analysis of large gene regulatory networks. By combining these algorithms with in silico knockouts and gene expression microarray data from human T cells, we examined if the results were compatible with a counter-regulatory role of Th1 and Th2 cells. We constructed a directed network model of genes regulating Th1 and Th2 cells through text mining and manual curation. We identified four attractors in the network, three of which included genes that corresponded to Th0, Th1 and Th2 cells. The fourth attractor contained a mixture of Th1 and Th2 genes. We found that neither in silico knockouts of the Th1 and Th2 attractor genes nor gene expression microarray data from patients with immunological disorders and healthy subjects supported a counter-regulatory role of Th1 and Th2 cells. By combining network modeling with transcriptomic data analysis and in silico knockouts, we have devised a practical way to help unravel complex regulatory network topology and to increase our understanding of how network actions may differ in health and disease.

Florfenicol inhibits allergic airway inflammation in mice by p38 MAPK-mediated phosphorylation of GATA 3

Florfenicol has been shown to possess anti-inflammatory activity. However, its possible use for asthma has not yet been studied. First we investigated the anti-inflammatory properties of florfenicol using mice asthma model. BALB/c mice were immunized and challenged by ovalbumin. Treatment with florfenicol caused a marked reduction in inflammatory cells and three Th2 type cytokines in the bronchoalveolar lavage fluids of mice. The levels of ovalbumin-specific IgE and airway hyperresponsiveness were significantly altered after treatment with florfenicol. Histological studies using H&E and AB-PAS staining demonstrate that florfenicol substantially inhibited ovalbumin-induced inflammatory cells infiltration in lung tissue and goblet cell hyperplasia in the airway. These results were similar to those obtained with dexamethasone treatment. We then investigated which signal transduction mechanisms could be implicated in florfenicol activity. Our results suggested that the protective effect of florfenicol was mediated by the inhibition of the p38 MAPK-mediated phosphorylation of GATA 3.

Gradient descent optimization in gene regulatory pathways

Background

Gene Regulatory Networks (GRNs) have become a major focus of interest in recent years. Elucidating the architecture and dynamics of large scale gene regulatory networks is an important goal in systems biology. The knowledge of the gene regulatory networks further gives insights about gene regulatory pathways. This information leads to many potential applications in medicine and molecular biology, examples of which are identification of metabolic pathways, complex genetic diseases, drug discovery and toxicology analysis. High-throughput technologies allow studying various aspects of gene regulatory networks on a genome-wide scale and we will discuss recent advances as well as limitations and future challenges for gene network modeling. Novel approaches are needed to both infer the causal genes and generate hypothesis on the underlying regulatory mechanisms.

Methodology

In the present article, we introduce a new method for identifying a set of optimal gene regulatory pathways by using structural equations as a tool for modeling gene regulatory networks. The method, first of all, generates data on reaction flows in a pathway. A set of constraints is formulated incorporating weighting coefficients. Finally the gene regulatory pathways are obtained through optimization of an objective function with respect to these weighting coefficients. The effectiveness of the present method is successfully tested on ten gene regulatory networks existing in the literature. A comparative study with the existing extreme pathway analysis also forms a part of this investigation. The results compare favorably with earlier experimental results. The validated pathways point to a combination of previously documented and novel findings.

Conclusions

We show that our method can correctly identify the causal genes and effectively output experimentally verified pathways. The present method has been successful in deriving the optimal regulatory pathways for all the regulatory networks considered. The biological significance and applicability of the optimal pathways has also been discussed. Finally the usefulness of the present method on genetic engineering is depicted with an example.

Integrative multicellular biological modeling: a case study of 3D epidermal development using GPU algorithms

BACKGROUND

Simulation of sophisticated biological models requires considerable computational power. These models typically integrate together numerous biological phenomena such as spatially-explicit heterogeneous cells, cell-cell interactions, cell-environment interactions and intracellular gene networks. The recent advent of programming for graphical processing units (GPU) opens up the possibility of developing more integrative, detailed and predictive biological models while at the same time decreasing the computational cost to simulate those models.

RESULTS

We construct a 3D model of epidermal development and provide a set of GPU algorithms that executes significantly faster than sequential central processing unit (CPU) code. We provide a parallel implementation of the subcellular element method for individual cells residing in a lattice-free spatial environment. Each cell in our epidermal model includes an internal gene network, which integrates cellular interaction of Notch signaling together with environmental interaction of basement membrane adhesion, to specify cellular state and behaviors such as growth and division. We take a pedagogical approach to describing how modeling methods are efficiently implemented on the GPU including memory layout of data structures and functional decomposition. We discuss various programmatic issues and provide a set of design guidelines for GPU programming that are instructive to avoid common pitfalls as well as to extract performance from the GPU architecture.

CONCLUSIONS

We demonstrate that GPU algorithms represent a significant technological advance for the simulation of complex biological models. We further demonstrate with our epidermal model that the integration of multiple complex modeling methods for heterogeneous multicellular biological processes is both feasible and computationally tractable using this new technology. We hope that the provided algorithms and source code will be a starting point for modelers to develop their own GPU implementations, and encourage others to implement their modeling methods on the GPU and to make that code available to the wider community.

Altered signaling through IL-12 receptor in children with very high serum IgE levels

An alteration of Th1/Th2 homeostasis may lead to diseases in humans. In this study, we investigated whether an impaired IL-12R signaling occurred in children with elevated serum IgE levels divided on the basis of the IgE levels (group A: >2000kU/l; group B: <2000kU/l). We evaluated the integrity of the IL-12R signaling through the analysis of phosphorylation/activation of STAT4, and mRNA expression and membrane assembly of the receptor chains. At a functional level, a proliferative defect of lymphocytes from group A patients was observed. In these patients, an abnormal IL-12R signaling was documented, and this finding was associated with abnormal expression of the IL-12Rbeta2 chain. Our data indicate that in patients with very high IgE levels the generation of Th1 response is impaired, and that this abnormality associates with abnormal IL-12R signaling.

Computational immunology meets bioinformatics: the use of prediction tools for molecular binding in the simulation of the immune system

We present a new approach to the study of the immune system that combines techniques of systems biology with information provided by data-driven prediction methods. To this end, we have extended an agent-based simulator of the immune response, C-ImmSim, such that it represents pathogens, as well as lymphocytes receptors, by means of their amino acid sequences and makes use of bioinformatics methods for T and B cell epitope prediction. This is a key step for the simulation of the immune response, because it determines immunogenicity. The binding of the epitope, which is the immunogenic part of an invading pathogen, together with activation and cooperation from T helper cells, is required to trigger an immune response in the affected host. To determine a pathogen's epitopes, we use existing prediction methods. In addition, we propose a novel method, which uses Miyazawa and Jernigan protein-protein potential measurements, for assessing molecular binding in the context of immune complexes. We benchmark the resulting model by simulating a classical immunization experiment that reproduces the development of immune memory. We also investigate the role of major histocompatibility complex (MHC) haplotype heterozygosity and homozygosity with respect to the influenza virus and show that there is an advantage to heterozygosity. Finally, we investigate the emergence of one or more dominating clones of lymphocytes in the situation of chronic exposure to the same immunogenic molecule and show that high affinity clones proliferate more than any other. These results show that the simulator produces dynamics that are stable and consistent with basic immunological knowledge. We believe that the combination of genomic information and simulation of the dynamics of the immune system, in one single tool, can offer new perspectives for a better understanding of the immune system.

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.

The development of locomotor kinematics in neonatal rats: an agent-based modeling analysis in group and individual contexts

An agent-based model of infant rat (pup) locomotion and aggregation was developed by modifying a previous model of pup aggregation [Schank, J.C., Alberts, J.R., 2000a. The developmental emergence of coupled activity as cooperative aggregation in rat pups. Proc. R. Soc. London B 267, 2307-2315]. The main difference between the earlier and current models is the incorporation of whole-body kinematics of directional locomotion. Data on locomotion and aggregation are presented for individuals and groups of 7- and 10-day-old pups and the data were used to evolve models (with a genetic algorithm) that fit these data. Aggregation between 7- and 10-day-old pups was considerably different and could be explained by agent-based models, in particular, models with directional-kinematic matrices specifying the probabilities of moving to adjacent cells. The directional kinematics of whole-body movement differed between the two age classes and differed between group and individual contexts for 10-day-old pups. This may indicate a developmental transition (by day 10) to more central control of behavior and the ability to change patterns of movement based on social context. The behavior analyzed with agent-based models may provide a precise way to measure motor and nervous system development in rats and other rodents.