Multiscale modeling of mucosal immune responses

The Future Exploring of Gut Microbiome-Immunity Interactions: From In Vivo/Vitro Models to In Silico Innovations

Investigating the complex interactions between microbiota and immunity is crucial for a fruitful understanding progress of human health and disease. This review assesses animal models, next-generation in vitro models, and in silico approaches that are used to decipher the microbiome-immunity axis, evaluating their strengths and limitations. While animal models provide a comprehensive biological context, they also raise ethical and practical concerns. Conversely, modern in vitro models reduce animal involvement but require specific costs and materials. When considering the environmental impact of these models, in silico approaches emerge as promising for resource reduction, but they require robust experimental validation and ongoing refinement. Their potential is significant, paving the way for a more sustainable and ethical future in microbiome-immunity research.

A multiscale mechanistic model of human dendritic cells for in-silico investigation of immune responses and novel therapeutics discovery

Dendritic cells (DCs) are professional antigen-presenting cells (APCs) with the unique ability to mediate inflammatory responses of the immune system. Given the critical role of DCs in shaping immunity, they present an attractive avenue as a therapeutic target to program the immune system and reverse immune disease disorders. To ensure appropriate immune response, DCs utilize intricate and complex molecular and cellular interactions that converge into a seamless phenotype. Computational models open novel frontiers in research by integrating large-scale interaction to interrogate the influence of complex biological behavior across scales. The ability to model large biological networks will likely pave the way to understanding any complex system in more approachable ways. We developed a logical and predictive model of DC function that integrates the heterogeneity of DCs population, APC function, and cell-cell interaction, spanning molecular to population levels. Our logical model consists of 281 components that connect environmental stimuli with various layers of the cell compartments, including the plasma membrane, cytoplasm, and nucleus to represent the dynamic processes within and outside the DC, such as signaling pathways and cell-cell interactions. We also provided three sample use cases to apply the model in the context of studying cell dynamics and disease environments. First, we characterized the DC response to Sars-CoV-2 and influenza co-infection by in-silico experiments and analyzed the activity level of 107 molecules that play a role in this co-infection. The second example presents simulations to predict the crosstalk between DCs and T cells in a cancer microenvironment. Finally, for the third example, we used the Kyoto Encyclopedia of Genes and Genomes enrichment analysis against the model's components to identify 45 diseases and 24 molecular pathways that the DC model can address. This study presents a resource to decode the complex dynamics underlying DC-derived APC communication and provides a platform for researchers to perform in-silico experiments on human DC for vaccine design, drug discovery, and immunotherapies.

Kinetics of antigen cross-presentation assessed experimentally and by a model of the complete endomembrane system

Dendritic cells (DCs) have a specialized endomembrane system capable of presenting exogenous antigens in the context of MHC class I (MHC-I) molecules. This process, named cross-presentation, is crucial to activate CD8+ T lymphocytes and initiate cytotoxic immune responses. In this report, we present an Agent-Based Model in combination with Ordinary Differential Equations with enough complexity to reproduce cross-presentation. The model embraces the secretory and endocytic pathways, in connection with the plasma membrane, the endoplasmic reticulum, and the cytosol. Key molecules required for cross-presentation were included as cargoes. In the simulations, the kinetics of MHC-I uptake and recycling, and cross-presentation accurately reproduced experimental values. The model proved to be a suitable tool to elaborate hypotheses and design experiments. In particular, the model predictions and the experimental results obtained indicate that the rate-limiting step in cross-presentation of soluble ovalbumin is MHC-I loading after proteasomal processing of the antigenic protein.

A roadmap for translational cancer glycoimmunology at single cell resolution

Cancer cells can evade immune responses by exploiting inhibitory immune checkpoints. Immune checkpoint inhibitor (ICI) therapies based on anti-CTLA-4 and anti-PD-1/PD-L1 antibodies have been extensively explored over the recent years to unleash otherwise compromised anti-cancer immune responses. However, it is also well established that immune suppression is a multifactorial process involving an intricate crosstalk between cancer cells and the immune systems. The cancer glycome is emerging as a relevant source of immune checkpoints governing immunosuppressive behaviour in immune cells, paving an avenue for novel immunotherapeutic options. This review addresses the current state-of-the-art concerning the role played by glycans controlling innate and adaptive immune responses, while shedding light on available experimental models for glycoimmunology. We also emphasize the tremendous progress observed in the development of humanized models for immunology, the paramount contribution of advances in high-throughput single-cell analysis in this context, and the importance of including predictive machine learning algorithms in translational research. This may constitute an important roadmap for glycoimmunology, supporting careful adoption of models foreseeing clinical translation of fundamental glycobiology knowledge towards next generation immunotherapies.

Nonlinear Reactor Design Optimization With Embedded Microkinetic Model Information

Despite the success of multiscale modeling in science and engineering, embedding molecular-level information into nonlinear reactor design and control optimization problems remains challenging. In this work, we propose a computationally tractable scale-bridging approach that incorporates information from multi-product microkinetic (MK) models with thousands of rates and chemical species into nonlinear reactor design optimization problems. We demonstrate reduced-order kinetic (ROK) modeling approaches for catalytic oligomerization in shale gas processing. We assemble a library of six candidate ROK models based on literature and MK model structure. We find that three metrics—quality of fit (e.g., mean squared logarithmic error), thermodynamic consistency (e.g., low conversion of exothermic reactions at high temperatures), and model identifiability—are all necessary to train and select ROK models. The ROK models that closely mimic the structure of the MK model offer the best compromise to emulate the product distribution. Using the four best ROK models, we optimize the temperature profiles in staged reactors to maximize conversions to heavier oligomerization products. The optimal temperature starts at 630–900K and monotonically decreases to approximately 560 K in the final stage, depending on the choice of ROK model. For all models, staging increases heavier olefin production by 2.5% and there is minimal benefit to more than four stages. The choice of ROK model, i.e., model-form uncertainty, results in a 22% difference in the objective function, which is twice the impact of parametric uncertainty; we demonstrate sequential eigendecomposition of the Fisher information matrix to identify and fix sloppy model parameters, which allows for more reliable estimation of the covariance of the identifiable calibrated model parameters. First-order uncertainty propagation determines this parametric uncertainty induces less than a 10% variability in the reactor optimization objective function. This result highlights the importance of quantifying model-form uncertainty, in addition to parametric uncertainty, in multi-scale reactor and process design and optimization. Moreover, the fast dynamic optimization solution times suggest the ROK strategy is suitable for incorporating molecular information in sequential modular or equation-oriented process simulation and optimization frameworks.

Computational models in inflammatory bowel disease

Inflammatory bowel disease (IBD) is a chronic and relapsing disease with multiple underlying influences and notable heterogeneity among its clinical and response-to-treatment phenotypes. There is no cure for IBD, and none of the currently available therapies have demonstrated clinical efficacies beyond 40%-60%. Data collected about its omics, pathogenesis, and treatment strategies have grown exponentially with time making IBD a prime candidate for artificial intelligence (AI) mediated discovery support. AI can be leveraged to further understand or identify IBD features to improve clinical outcomes. Various treatment candidates are currently under evaluation in clinical trials, offering further approaches and opportunities for increasing the efficacies of treatments. However, currently, therapeutic plans are largely determined using clinical features due to the lack of specific biomarkers, and it has become necessary to step into precision medicine to predict therapeutic responses to guarantee optimal treatment efficacy. This is accompanied by the application of AI and the development of multiscale hybrid models combining mechanistic approaches and machine learning. These models ultimately lead to the creation of digital twins of given patients delivering on the promise of precision dosing and tailored treatment. Interleukin-6 (IL-6) is a prominent cytokine in cell-to-cell communication in the inflammatory responses' regulation. Dysregulated IL-6-induced signaling leads to severe immunological or proliferative pathologies, such as IBD and colon cancer. This mini-review explores multiscale models with the aim of predicting the response to therapy in IBD. Modeling IL-6 biology and generating digital twins enhance the credibility of their prediction.

A Quantitative Systems Pharmacology Model for the Key Interleukins Involved in Crohn's Disease

Crohn's disease (CD) is a complex inflammatory bowel disease whose pathogenesis appears to involve several immunologic defects causing functional impairment of the gut. Its complexity and the reported loss of effectiveness over time of standard of care together with the increase in its worldwide incidence require the application of techniques aiming to find new therapeutic strategies. Currently, systems pharmacology modeling has been gaining importance as it integrates the available knowledge of the system into a single computational model. In this work, the following workflow for robust application of systems pharmacology modeling was followed: 1) scope definition; 2) species selection and circulating plasma levels based on a search in the literature; 3) representation of model topology and parametrization of the interactions, after literature data extraction and curation, and the implementation of ordinary differential equations in SimBiology (MATLAB version R2018b); and 4) model curation and evaluation by visual comparison of simulated interleukin (IL) concentrations with the reported levels in plasma, and sensitivity analysis performed to confirm model robustness and identify the most influential parameters. Finally, 5) exposure to two dose levels of recombinant human IL10 was evaluated by simulation and comparison with reported clinical study results. In summary, we present a quantitative systems pharmacology model for the main ILs involved in CD developed using a standardized methodology and supported by a comprehensive repository summarizing the most relevant literature in the field. However, it has to be taken into account that external validation is still pending as available clinical data were primarily used for model training. SIGNIFICANCE STATEMENT: Crohn's disease (CD) is a complex heterogeneous inflammatory bowel disorder. Systems pharmacology modeling offers a great opportunity for integration of the available knowledge on the disease using a computational framework. As a result of this work, a comprehensive repository along with a quantitative systems pharmacology model for the main interleukins involved in CD is provided. This model is useful for the in silico evaluation of biomarkers and potential therapeutic targets and can be adapted to address research gaps regarding CD.

High-resolution computational modeling of immune responses in the gut

Background

Helicobacter pylori causes gastric cancer in 1–2% of cases but is also beneficial for protection against allergies and gastroesophageal diseases. An estimated 85% of H. pylori–colonized individuals experience no detrimental effects. To study the mechanisms promoting host tolerance to the bacterium in the gastrointestinal mucosa and systemic regulatory effects, we investigated the dynamics of immunoregulatory mechanisms triggered by H. pylori using a high-performance computing–driven ENteric Immunity SImulator multiscale model. Immune responses were simulated by integrating an agent-based model, ordinary, and partial differential equations.

Results

The outputs were analyzed using 2 sequential stages: the first used a partial rank correlation coefficient regression–based and the second a metamodel-based global sensitivity analysis. The influential parameters screened from the first stage were selected to be varied for the second stage. The outputs from both stages were combined as a training dataset to build a spatiotemporal metamodel. The Sobol indices measured time-varying impact of input parameters during initiation, peak, and chronic phases of infection. The study identified epithelial cell proliferation and epithelial cell death as key parameters that control infection outcomes. In silico validation showed that colonization with H. pylori decreased with a decrease in epithelial cell proliferation, which was linked to regulatory macrophages and tolerogenic dendritic cells.

Conclusions

The hybrid model of H. pylori infection identified epithelial cell proliferation as a key factor for successful colonization of the gastric niche and highlighted the role of tolerogenic dendritic cells and regulatory macrophages in modulating the host responses and shaping infection outcomes.

Parallelisation strategies for agent based simulation of immune systems

BACKGROUND

In recent years, the study of immune response behaviour using bottom up approach, Agent Based Modeling (ABM), has attracted considerable efforts. The ABM approach is a very common technique in the biological domain due to high demand for a large scale analysis tools for the collection and interpretation of information to solve biological problems. Simulating massive multi-agent systems (i.e. simulations containing a large number of agents/entities) requires major computational effort which is only achievable through the use of parallel computing approaches.

RESULTS

This paper explores different approaches to parallelising the key component of biological and immune system models within an ABM model: pairwise interactions. The focus of this paper is on the performance and algorithmic design choices of cell interactions in continuous and discrete space where agents/entities are competing to interact with one another within a parallel environment.

CONCLUSIONS

Our performance results demonstrate the applicability of these methods to a broader class of biological systems exhibiting typical cell to cell interactions. The advantage and disadvantage of each implementation is discussed showing each can be used as the basis for developing complete immune system models on parallel hardware.

Agent-based models of inflammation in translational systems biology: A decade later

Agent-based modeling is a rule-based, discrete-event, and spatially explicit computational modeling method that employs computational objects that instantiate the rules and interactions among the individual components ("agents") of system. Agent-based modeling is well suited to translating into a computational model the knowledge generated from basic science research, particularly with respect to translating across scales the mechanisms of cellular behavior into aggregated cell population dynamics manifesting at the tissue and organ level. This capacity has made agent-based modeling an integral method in translational systems biology (TSB), an approach that uses multiscale dynamic computational modeling to explicitly represent disease processes in a clinically relevant fashion. The initial work in the early 2000s using agent-based models (ABMs) in TSB focused on examining acute inflammation and its intersection with wound healing; the decade since has seen vast growth in both the application of agent-based modeling to a wide array of disease processes as well as methodological advancements in the use and analysis of ABM. This report presents an update on an earlier review of ABMs in TSB and presents examples of exciting progress in the modeling of various organs and diseases that involve inflammation. This review also describes developments that integrate the use of ABMs with cutting-edge technologies such as high-performance computing, machine learning, and artificial intelligence, with a view toward the future integration of these methodologies. This article is categorized under: Translational, Genomic, and Systems Medicine > Translational Medicine Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models Models of Systems Properties and Processes > Organismal Models.

Review of the systems biology of the immune system using agent-based models

The immune system is an inherent protection system in vertebrate animals including human beings that exhibit properties such as self-organisation, self-adaptation, learning, and recognition. It interacts with the other allied systems such as the gut and lymph nodes. There is a need for immune system modelling to know about its complex internal mechanism, to understand how it maintains the homoeostasis, and how it interacts with the other systems. There are two types of modelling techniques used for the simulation of features of the immune system: equation-based modelling (EBM) and agent-based modelling. Owing to certain shortcomings of the EBM, agent-based modelling techniques are being widely used. This technique provides various predictions for disease causes and treatments; it also helps in hypothesis verification. This study presents a review of agent-based modelling of the immune system and its interactions with the gut and lymph nodes. The authors also review the modelling of immune system interactions during tuberculosis and cancer. In addition, they also outline the future research directions for the immune system simulation through agent-based techniques such as the effects of stress on the immune system, evolution of the immune system, and identification of the parameters for a healthy immune system.

A systems pharmacology model for inflammatory bowel disease

Motivation

The literature on complex diseases is abundant but not always quantitative. This is particularly so for Inflammatory Bowel Disease (IBD), where many molecular pathways are qualitatively well described but this information cannot be used in traditional quantitative mathematical models employed in drug development. We propose the elaboration and validation of a logic network for IBD able to capture the information available in the literature that will facilitate the identification/validation of therapeutic targets.

Results

In this article, we propose a logic model for Inflammatory Bowel Disease (IBD) which consists of 43 nodes and 298 qualitative interactions. The model presented is able to describe the pathogenic mechanisms of the disorder and qualitatively describes the characteristic chronic inflammation. A perturbation analysis performed on the IBD network indicates that the model is robust. Also, as described in clinical trials, a simulation of anti-TNFα, anti-IL2 and Granulocyte and Monocyte Apheresis showed a decrease in the Metalloproteinases node (MMPs), which means a decrease in tissue damage. In contrast, as clinical trials have demonstrated, a simulation of anti-IL17 and anti-IFNγ or IL10 overexpression therapy did not show any major change in MMPs expression, as corresponds to a failed therapy. The model proved to be a promising in silico tool for the evaluation of potential therapeutic targets, the identification of new IBD biomarkers, the integration of IBD polymorphisms to anticipate responders and non-responders and can be reduced and transformed in quantitative model/s.

Informatics for Nutritional Genetics and Genomics

While traditional nutrition science is focusing on nourishing population, modern nutrition is aiming at benefiting individual people. The goal of modern nutritional research is to promote health, prevent diseases, and improve performance. With the development of modern technologies like bioinformatics, metabolomics, and molecular genetics, this goal is becoming more attainable. In this chapter, we will discuss the new concepts and technologies especially in informatics and molecular genetics and genomics, and how they have been implemented to change the nutrition science and lead to the emergence of new branches like nutrigenomics, nutrigenetics, and nutritional metabolomics.

Model of Host-Pathogen Interaction Dynamics Links In Vivo Optical Imaging and Immune Responses

Tracking disease progression in vivo is essential for the development of treatments against bacterial infection. Optical imaging has become a central tool for in vivo tracking of bacterial population development and therapeutic response. For a precise understanding of in vivo imaging results in terms of disease mechanisms derived from detailed postmortem observations, however, a link between the two is needed. Here, we develop a model that provides that link for the investigation of Citrobacter rodentium infection, a mouse model for enteropathogenic Escherichia coli (EPEC). We connect in vivo disease progression of C57BL/6 mice infected with bioluminescent bacteria, imaged using optical tomography and X-ray computed tomography, to postmortem measurements of colonic immune cell infiltration. We use the model to explore changes to both the host immune response and the bacteria and to evaluate the response to antibiotic treatment. The developed model serves as a novel tool for the identification and development of new therapeutic interventions.

Solving Immunology?

Emergent responses of the immune system result from the integration of molecular and cellular networks over time and across multiple organs. High-content and high-throughput analysis technologies, concomitantly with data-driven and mechanistic modeling, hold promise for the systematic interrogation of these complex pathways. However, connecting genetic variation and molecular mechanisms to individual phenotypes and health outcomes has proven elusive. Gaps remain in data, and disagreements persist about the value of mechanistic modeling for immunology. Here, we present the perspectives that emerged from the National Institute of Allergy and Infectious Disease (NIAID) workshop 'Complex Systems Science, Modeling and Immunity' and subsequent discussions regarding the potential synergy of high-throughput data acquisition, data-driven modeling, and mechanistic modeling to define new mechanisms of immunological disease and to accelerate the translation of these insights into therapies.

Modeling the Role of Lanthionine Synthetase C-Like 2 (LANCL2) in the Modulation of Immune Responses to Helicobacter pylori Infection

Immune responses to Helicobacter pylori are orchestrated through complex balances of host-bacterial interactions, including inflammatory and regulatory immune responses across scales that can lead to the development of the gastric disease or the promotion of beneficial systemic effects. While inflammation in response to the bacterium has been reasonably characterized, the regulatory pathways that contribute to preventing inflammatory events during H. pylori infection are incompletely understood. To aid in this effort, we have generated a computational model incorporating recent developments in the understanding of H. pylori-host interactions. Sensitivity analysis of this model reveals that a regulatory macrophage population is critical in maintaining high H. pylori colonization without the generation of an inflammatory response. To address how this myeloid cell subset arises, we developed a second model describing an intracellular signaling network for the differentiation of macrophages. Modeling studies predicted that LANCL2 is a central regulator of inflammatory and effector pathways and its activation promotes regulatory responses characterized by IL-10 production while suppressing effector responses. The predicted impairment of regulatory macrophage differentiation by the loss of LANCL2 was simulated based on multiscale linkages between the tissue-level gastric mucosa and the intracellular models. The simulated deletion of LANCL2 resulted in a greater clearance of H. pylori, but also greater IFNγ responses and damage to the epithelium. The model predictions were validated within a mouse model of H. pylori colonization in wild-type (WT), LANCL2 whole body KO and myeloid-specific LANCL2-/- (LANCL2^Myeloid) mice, which displayed similar decreases in H. pylori burden, CX3CR1+ IL-10-producing macrophages, and type 1 regulatory (Tr1) T cells. This study shows the importance of LANCL2 in the induction of regulatory responses in macrophages and T cells during H. pylori infection.

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.

Mucosal Interactions between Genetics, Diet, and Microbiome in Inflammatory Bowel Disease

Numerous reviews have discussed gut microbiota composition changes during inflammatory bowel diseases (IBD), particularly Crohn’s disease (CD). However, most studies address the observed effects by focusing on studying the univariate connection between disease and dietary-induced alterations to gut microbiota composition. The possibility that these effects may reflect a number of other interconnected (i.e., pantropic) mechanisms, activated in parallel, particularly concerning various bacterial metabolites, is in the process of being elucidated. Progress seems, however, hampered by various difficult-to-study factors interacting at the mucosal level. Here, we highlight some of such factors that merit consideration, namely: (1) the contribution of host genetics and diet in altering gut microbiome, and in turn, the crosstalk among secondary metabolic pathways; (2) the interdependence between the amount of dietary fat, the fatty acid composition, the effects of timing and route of administration on gut microbiota community, and the impact of microbiota-derived fatty acids; (3) the effect of diet on bile acid composition, and the modulator role of bile acids on the gut microbiota; (4) the impact of endogenous and exogenous intestinal micronutrients and metabolites; and (5) the need to consider food associated toxins and chemicals, which can introduce confounding immune modulating elements (e.g., antioxidant and phytochemicals in oils and proteins). These concepts, which are not mutually exclusive, are herein illustrated paying special emphasis on physiologically inter-related processes.

Multi-class and multi-scale models of complex biological phenomena

Computational modeling has significantly impacted our ability to analyze vast (and exponentially increasing) quantities of experimental data for a variety of applications, such as drug discovery and disease forecasting. Single-scale, single-class models persist as the most common group of models, but biological complexity often demands more sophisticated approaches. This review surveys modeling approaches that are multi-class (incorporating multiple model types) and/or multi-scale (accounting for multiple spatial or temporal scales) and describes how these models, and combinations thereof, should be used within the context of the problem statement. We end by highlighting agent-based models as an intuitive, modular, and flexible framework within which multi-scale and multi-class models can be implemented.

Modeling-Enabled Systems Nutritional Immunology

This review highlights the fundamental role of nutrition in the maintenance of health, the immune response, and disease prevention. Emerging global mechanistic insights in the field of nutritional immunology cannot be gained through reductionist methods alone or by analyzing a single nutrient at a time. We propose to investigate nutritional immunology as a massively interacting system of interconnected multistage and multiscale networks that encompass hidden mechanisms by which nutrition, microbiome, metabolism, genetic predisposition, and the immune system interact to delineate health and disease. The review sets an unconventional path to apply complex science methodologies to nutritional immunology research, discovery, and development through “use cases” centered around the impact of nutrition on the gut microbiome and immune responses. Our systems nutritional immunology analyses, which include modeling and informatics methodologies in combination with pre-clinical and clinical studies, have the potential to discover emerging systems-wide properties at the interface of the immune system, nutrition, microbiome, and metabolism.

Sensitivity Analysis of an ENteric Immunity SImulator (ENISI)-Based Model of Immune Responses to Helicobacter pylori Infection

Agent-based models (ABM) are widely used to study immune systems, providing a procedural and interactive view of the underlying system. The interaction of components and the behavior of individual objects is described procedurally as a function of the internal states and the local interactions, which are often stochastic in nature. Such models typically have complex structures and consist of a large number of modeling parameters. Determining the key modeling parameters which govern the outcomes of the system is very challenging. Sensitivity analysis plays a vital role in quantifying the impact of modeling parameters in massively interacting systems, including large complex ABM. The high computational cost of executing simulations impedes running experiments with exhaustive parameter settings. Existing techniques of analyzing such a complex system typically focus on local sensitivity analysis, i.e. one parameter at a time, or a close “neighborhood” of particular parameter settings. However, such methods are not adequate to measure the uncertainty and sensitivity of parameters accurately because they overlook the global impacts of parameters on the system. In this article, we develop novel experimental design and analysis techniques to perform both global and local sensitivity analysis of large-scale ABMs. The proposed method can efficiently identify the most significant parameters and quantify their contributions to outcomes of the system. We demonstrate the proposed methodology for ENteric Immune SImulator (ENISI), a large-scale ABM environment, using a computational model of immune responses to Helicobacter pylori colonization of the gastric mucosa.

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

Clostridium difficile infections are associated with the use of broad-spectrum antibiotics and result in an exuberant inflammatory response, leading to nosocomial diarrhea, colitis and even death. To better understand the dynamics of mucosal immunity during C. difficile infection from initiation through expansion to resolution, we built a computational model of the mucosal immune response to the bacterium. The model was calibrated using data from a mouse model of C. difficile infection. The model demonstrates a crucial role of T helper 17 (Th17) effector responses in the colonic lamina propria and luminal commensal bacteria populations in the clearance of C. difficile and colonic pathology, whereas regulatory T (Treg) cells responses are associated with the recovery phase. In addition, the production of anti-microbial peptides by inflamed epithelial cells and activated neutrophils in response to C. difficile infection inhibit the re-growth of beneficial commensal bacterial species. Computational simulations suggest that the removal of neutrophil and epithelial cell derived anti-microbial inhibitions, separately and together, on commensal bacterial regrowth promote recovery and minimize colonic inflammatory pathology. Simulation results predict a decrease in colonic inflammatory markers, such as neutrophilic influx and Th17 cells in the colonic lamina propria, and length of infection with accelerated commensal bacteria re-growth through altered anti-microbial inhibition. Computational modeling provides novel insights on the therapeutic value of repopulating the colonic microbiome and inducing regulatory mucosal immune responses during C. difficile infection. Thus, modeling mucosal immunity-gut microbiota interactions has the potential to guide the development of targeted fecal transplantation therapies in the context of precision medicine interventions.