Modeling intercellular interactions in early Mycobacterium infection

In silico agent-based modeling approach to characterize multiple in vitro tuberculosis infection models

In vitro models of Mycobacterium tuberculosis (Mtb) infection are a valuable tool for examining host-pathogen interactions and screening drugs. With the development of more complex in vitro models, there is a need for tools to help analyze and integrate data from these models. To this end, we introduce an agent-based model (ABM) representation of the interactions between immune cells and bacteria in an in vitro setting. This in silico model was used to simulate both traditional and spheroid cell culture models by changing the movement rules and initial spatial layout of the cells in accordance with the respective in vitro models. The traditional and spheroid simulations were calibrated to published experimental data in a paired manner, by using the same parameters in both simulations. Within the calibrated simulations, heterogeneous outputs are seen for bacterial count and T cell infiltration into the macrophage core of the spheroid. The simulations also predict that equivalent numbers of activated macrophages do not necessarily result in similar bacterial reductions; that host immune responses can control bacterial growth in both spheroid structure dependent and independent manners; that STAT1 activation is the limiting step in macrophage activation in spheroids; and that drug screening and macrophage activation studies could have different outcomes depending on the in vitro culture used. Future model iterations will be guided by the limitations of the current model, specifically which parts of the output space were harder to reach. This ABM can be used to represent more in vitro Mtb infection models due to its flexible structure, thereby accelerating in vitro discoveries.

Prediction of Tuberculosis From Lung Tissue Images of Diversity Outbred Mice Using Jump Knowledge Based Cell Graph Neural Network

Tuberculosis (TB), primarily affecting the lungs, is caused by the bacterium Mycobacterium tuberculosis and poses a significant health risk. Detecting acid-fast bacilli (AFB) in stained samples is critical for TB diagnosis. Whole Slide (WS) Imaging allows for digitally examining these stained samples. However, current deep-learning approaches to analyzing large-sized whole slide images (WSIs) often employ patch-wise analysis, potentially missing the complex spatial patterns observed in the granuloma essential for accurate TB classification. To address this limitation, we propose an approach that models cell characteristics and interactions as a graph, capturing both cell-level information and the overall tissue micro-architecture. This method differs from the strategies in related cell graph-based works that rely on edge thresholds based on sparsity/density in cell graph construction, emphasizing a biologically informed threshold determination instead. We introduce a cell graph-based jumping knowledge neural network (CG-JKNN) that operates on the cell graphs where the edge thresholds are selected based on the length of the mycobacteria's cords and the activated macrophage nucleus's size to reflect the actual biological interactions observed in the tissue. The primary process involves training a Convolutional Neural Network (CNN) to segment AFBs and macrophage nuclei, followed by converting large (42831*41159 pixels) lung histology images into cell graphs where an activated macrophage nucleus/AFB represents each node within the graph and their interactions are denoted as edges. To enhance the interpretability of our model, we employ Integrated Gradients and Shapely Additive Explanations (SHAP). Our analysis incorporated a combination of 33 graph metrics and 20 cell morphology features. In terms of traditional machine learning models, Extreme Gradient Boosting (XGBoost) was the best performer, achieving an F1 score of 0.9813 and an Area under the Precision-Recall Curve (AUPRC) of 0.9848 on the test set. Among graph-based models, our CG-JKNN was the top performer, attaining an F1 score of 0.9549 and an AUPRC of 0.9846 on the held-out test set. The integration of graph-based and morphological features proved highly effective, with CG-JKNN and XGBoost showing promising results in classifying instances into AFB and activated macrophage nucleus. The features identified as significant by our models closely align with the criteria used by pathologists in practice, highlighting the clinical applicability of our approach. Future work will explore knowledge distillation techniques and graph-level classification into distinct TB progression categories.

Aerosol Transport Modeling: The Key Link Between Lung Infections of Individuals and Populations

The recent COVID-19 pandemic has propelled the field of aerosol science to the forefront, particularly the central role of virus-laden respiratory droplets and aerosols. The pandemic has also highlighted the critical need, and value for, an information bridge between epidemiological models (that inform policymakers to develop public health responses) and within-host models (that inform the public and health care providers how individuals develop respiratory infections). Here, we review existing data and models of generation of respiratory droplets and aerosols, their exhalation and inhalation, and the fate of infectious droplet transport and deposition throughout the respiratory tract. We then articulate how aerosol transport modeling can serve as a bridge between and guide calibration of within-host and epidemiological models, forming a comprehensive tool to formulate and test hypotheses about respiratory tract exposure and infection within and between individuals.

Modeling insights into SARS-CoV-2 respiratory tract infections prior to immune protection

Mechanistic insights into human respiratory tract (RT) infections from SARS-CoV-2 can inform public awareness as well as guide medical prevention and treatment for COVID-19 disease. Yet the complexity of the RT and the inability to access diverse regions pose fundamental roadblocks to evaluation of potential mechanisms for the onset and progression of infection (and transmission). We present a model that incorporates detailed RT anatomy and physiology, including airway geometry, physical dimensions, thicknesses of airway surface liquids (ASLs), and mucus layer transport by cilia. The model further incorporates SARS-CoV-2 diffusivity in ASLs and best-known data for epithelial cell infection probabilities, and, once infected, duration of eclipse and replication phases, and replication rate of infectious virions. We apply this baseline model in the absence of immune protection to explore immediate, short-term outcomes from novel SARS-CoV-2 depositions onto the air-ASL interface. For each RT location, we compute probability to clear versus infect; per infected cell, we compute dynamics of viral load and cell infection. Results reveal that nasal infections are highly likely within 1-2 days from minimal exposure, and alveolar pneumonia occurs only if infectious virions are deposited directly into alveolar ducts and sacs, not via retrograde propagation to the deep lung. Furthermore, to infect just 1% of the 140 m2 of alveolar surface area within 1 week, either 103 boluses each with 106 infectious virions or 106 aerosols with one infectious virion, all physically separated, must be directly deposited. These results strongly suggest that COVID-19 disease occurs in stages: a nasal/upper RT infection, followed by self-transmission of infection to the deep lung. Two mechanisms of self-transmission are persistent aspiration of infected nasal boluses that drain to the deep lung and repeated rupture of nasal aerosols from infected mucosal membranes by speaking, singing, or cheering that are partially inhaled, exhaled, and re-inhaled, to the deep lung.

Spatially distributed infection increases viral load in a computational model of SARS-CoV-2 lung infection

A key question in SARS-CoV-2 infection is why viral loads and patient outcomes vary dramatically across individuals. Because spatial-temporal dynamics of viral spread and immune response are challenging to study in vivo, we developed Spatial Immune Model of Coronavirus (SIMCoV), a scalable computational model that simulates hundreds of millions of lung cells, including respiratory epithelial cells and T cells. SIMCoV replicates viral growth dynamics observed in patients and shows how spatially dispersed infections can lead to increased viral loads. The model also shows how the timing and strength of the T cell response can affect viral persistence, oscillations, and control. By incorporating spatial interactions, SIMCoV provides a parsimonious explanation for the dramatically different viral load trajectories among patients by varying only the number of initial sites of infection and the magnitude and timing of the T cell immune response. When the branching airway structure of the lung is explicitly represented, we find that virus spreads faster than in a 2D layer of epithelial cells, but much more slowly than in an undifferentiated 3D grid or in a well-mixed differential equation model. These results illustrate how realistic, spatially explicit computational models can improve understanding of within-host dynamics of SARS-CoV-2 infection.

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

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

An agent-based simulator for the gastrointestinal pathway of Listeria monocytogenes

We developed an agent-based gastric simulator for a human host to illustrate the within host survival mechanisms of Listeria monocytogenes. The simulator incorporates the gastric physiology and digestion processes that are critical for pathogen survival in the stomach. Mathematical formulations for the pH dynamics, stomach emptying time, and survival probability in the presence of gastric acid are integrated in the simulator to evaluate the portion of ingested bacteria that survives in the stomach and reaches the small intestine. The parameters are estimated using in vitro data relevant to the human stomach and L. monocytogenes. The simulator predicts that 5%-29% of ingested bacteria can survive a human stomach and reach the small intestine. In the absence of extensive scientific experiments, which are not feasible on the grounds of ethical and safety concerns, this simulator may provide a supplementary tool to evaluate pathogen survival and subsequent infection, especially with regards to the ingestion of small doses.

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

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

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

A mathematical model shows macrophages delay Staphylococcus aureus replication, but limitations in microbicidal capacity restrict bacterial clearance

S. aureus is a leading cause of bacterial infection. Macrophages, the first line of defence in the human immune response, phagocytose and kill S. aureus but the pathogen can evade these responses. Therefore, the exact role of macrophages is incompletely defined. We develop a mathematical model of macrophage - S. aureus dynamics, built on recent experimental data. We demonstrate that, while macrophages may not clear infection, they significantly delay its growth and potentially buy time for recruitment of further cells. We find that macrophage killing is a major obstacle to controlling infection and ingestion capacity also limits the response. We find bistability such that the infection can be limited at low doses. Our combination of experimental data, mathematical analysis and model fitting provide important insights in to the early stages of S. aureus infections, showing macrophages play an important role limiting bacterial replication but can be overwhelmed with large inocula.

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

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

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.

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.

Strategies for calibrating models of biology

Computational and mathematical modelling has become a valuable tool for investigating biological systems. Modelling enables prediction of how biological components interact to deliver system-level properties and extrapolation of biological system performance to contexts and experimental conditions where this is unknown. A model's value hinges on knowing that it faithfully represents the biology under the contexts of use, or clearly ascertaining otherwise and thus motivating further model refinement. These qualities are evaluated through calibration, typically formulated as identifying model parameter values that align model and biological behaviours as measured through a metric applied to both. Calibration is critical to modelling but is often underappreciated. A failure to appropriately calibrate risks unrepresentative models that generate erroneous insights. Here, we review a suite of strategies to more rigorously challenge a model's representation of a biological system. All are motivated by features of biological systems, and illustrative examples are drawn from the modelling literature. We examine the calibration of a model against distributions of biological behaviours or outcomes, not only average values. We argue for calibration even where model parameter values are experimentally ascertained. We explore how single metrics can be non-distinguishing for complex systems, with multiple-component dynamic and interaction configurations giving rise to the same metric output. Under these conditions, calibration is insufficiently constraining and the model non-identifiable: multiple solutions to the calibration problem exist. We draw an analogy to curve fitting and argue that calibrating a biological model against a single experiment or context is akin to curve fitting against a single data point. Though useful for communicating model results, we explore how metrics that quantify heavily emergent properties may not be suitable for use in calibration. Lastly, we consider the role of sensitivity and uncertainty analysis in calibration and the interpretation of model results. Our goal in this manuscript is to encourage a deeper consideration of calibration, and how to increase its capacity to either deliver faithful models or demonstrate them otherwise.

Understanding the pathophysiology of the human TB lung granuloma using in vitro granuloma models

Tuberculosis remains a major human health threat that infects one in three individuals worldwide. Infection with Mycobacterium tuberculosis is a standoff between host and bacteria in the formation of a granuloma. This review will introduce a variety of bacterial and host factors that impact individual granuloma fates. The authors describe advances in the development of in vitro granuloma models, current evidence surrounding infection and granuloma development, and the applicability of existing in vitro models in the study of human disease. In vitro models of infection help improve our understanding of pathophysiology and allow for the discovery of other potential models of study.

A spatial model of the efficiency of T cell search in the influenza-infected lung

Emerging strains of influenza, such as avian H5N1 and 2009 pandemic H1N1, are more virulent than seasonal H1N1 influenza, yet the underlying mechanisms for these differences are not well understood. Subtle differences in how a given strain interacts with the immune system are likely a key factor in determining virulence. One aspect of the interaction is the ability of T cells to locate the foci of the infection in time to prevent uncontrolled expansion. Here, we develop an agent based spatial model to focus on T cell migration from lymph nodes through the vascular system to sites of infection. We use our model to investigate whether different strains of influenza modulate this process. We calibrate the model using viral and chemokine secretion rates we measure in vitro together with values taken from literature. The spatial nature of the model reveals unique challenges for T cell recruitment that are not apparent in standard differential equation models. In this model comparing three influenza viruses, plaque expansion is governed primarily by the replication rate of the virus strain, and the efficiency of the T cell search-and-kill is limited by the density of infected epithelial cells in each plaque. Thus for each virus there is a different threshold of T cell search time above which recruited T cells are unable to control further expansion. Future models could use this relationship to more accurately predict control of the infection.

Oxygen Modulates the Effectiveness of Granuloma Mediated Host Response to Mycobacterium tuberculosis: A Multiscale Computational Biology Approach

Mycobacterium tuberculosis associated granuloma formation can be viewed as a structural immune response that can contain and halt the spread of the pathogen. In several mammalian hosts, including non-human primates, Mtb granulomas are often hypoxic, although this has not been observed in wild type murine infection models. While a presumed consequence, the structural contribution of the granuloma to oxygen limitation and the concomitant impact on Mtb metabolic viability and persistence remains to be fully explored. We develop a multiscale computational model to test to what extent in vivo Mtb granulomas become hypoxic, and investigate the effects of hypoxia on host immune response efficacy and mycobacterial persistence. Our study integrates a physiological model of oxygen dynamics in the extracellular space of alveolar tissue, an agent-based model of cellular immune response, and a systems biology-based model of Mtb metabolic dynamics. Our theoretical studies suggest that the dynamics of granuloma organization mediates oxygen availability and illustrates the immunological contribution of this structural host response to infection outcome. Furthermore, our integrated model demonstrates the link between structural immune response and mechanistic drivers influencing Mtbs adaptation to its changing microenvironment and the qualitative infection outcome scenarios of clearance, containment, dissemination, and a newly observed theoretical outcome of transient containment. We observed hypoxic regions in the containment granuloma similar in size to granulomas found in mammalian in vivo models of Mtb infection. In the case of the containment outcome, our model uniquely demonstrates that immune response mediated hypoxic conditions help foster the shift down of bacteria through two stages of adaptation similar to thein vitro non-replicating persistence (NRP) observed in the Wayne model of Mtb dormancy. The adaptation in part contributes to the ability of Mtb to remain dormant for years after initial infection.

The importance of liver microcirculation in promoting autoimmune hepatitis via maintaining an inflammatory cytokine milieu--a mathematical model study

In autoimmune diseases, inflammatory cytokine concentrations are important for initiating and maintaining the status of autoimmunity. Autoimmune hepatitis (AIH) is an inflammatory liver disease characterized by a loss of immune tolerance against specific antigens located in hepatocytes. During the progression of the disease, antigen-presenting cells and different classes of T-helper cells secrete specific cytokines important for maintaining the disease. As these cytokines are secreted into the local liver environment, the blood flow in liver sinusoids might influence the local cytokine concentration. Considering the liver tissue as a porous medium, based on Darcy׳s law, the microcirculation within a liver lobule was modelled. Using realistic physiological pressure differences and tissue permeabilities, the blood velocity inside the sinusoids could be calculated and validated with blood velocity data obtained via Orthogonal Polarization Spectral Imaging (OPSI). Furthermore, oxygen consumption is modelled to obtain Rappaport׳s acinus model. Finally, steady state spatial distributions of secreted cytokines within the liver lobule could be estimated for specified realistic production rates of T-helper cells. It could be demonstrated that the characteristics of the liver microcirculation have an important impact on establishing inflammatory cytokine levels within the portal fields and the vascular septa promoting the occurrence of interface hepatitis.

An agent-based model of cellular dynamics and circadian variability in human endotoxemia

As cellular variability and circadian rhythmicity play critical roles in immune and inflammatory responses, we present in this study an agent-based model of human endotoxemia to examine the interplay between circadian controls, cellular variability and stochastic dynamics of inflammatory cytokines. The model is qualitatively validated by its ability to reproduce circadian dynamics of inflammatory mediators and critical inflammatory responses after endotoxin administration in vivo. Novel computational concepts are proposed to characterize the cellular variability and synchronization of inflammatory cytokines in a population of heterogeneous leukocytes. Our results suggest that there is a decrease in cell-to-cell variability of inflammatory cytokines while their synchronization is increased after endotoxin challenge. Model parameters that are responsible for IκB production stimulated by NFκB activation and for the production of anti-inflammatory cytokines have large impacts on system behaviors. Additionally, examining time-dependent systemic responses revealed that the system is least vulnerable to endotoxin in the early morning and most vulnerable around midnight. Although much remains to be explored, proposed computational concepts and the model we have pioneered will provide important insights for future investigations and extensions, especially for single-cell studies to discover how cellular variability contributes to clinical implications.

Agent-based dynamic knowledge representation of Pseudomonas aeruginosa virulence activation in the stressed gut: Towards characterizing host-pathogen interactions in gut-derived sepsis

Background

There is a growing realization that alterations in host-pathogen interactions (HPI) can generate disease phenotypes without pathogen invasion. The gut represents a prime region where such HPI can arise and manifest. Under normal conditions intestinal microbial communities maintain a stable, mutually beneficial ecosystem. However, host stress can lead to changes in environmental conditions that shift the nature of the host-microbe dialogue, resulting in escalation of virulence expression, immune activation and ultimately systemic disease. Effective modulation of these dynamics requires the ability to characterize the complexity of the HPI, and dynamic computational modeling can aid in this task. Agent-based modeling is a computational method that is suited to representing spatially diverse, dynamical systems. We propose that dynamic knowledge representation of gut HPI with agent-based modeling will aid in the investigation of the pathogenesis of gut-derived sepsis.

Methodology/Principal Findings

An agent-based model (ABM) of virulence regulation in Pseudomonas aeruginosa was developed by translating bacterial and host cell sense-and-response mechanisms into behavioral rules for computational agents and integrated into a virtual environment representing the host-microbe interface in the gut. The resulting gut milieu ABM (GMABM) was used to: 1) investigate a potential clinically relevant laboratory experimental condition not yet developed - i.e. non-lethal transient segmental intestinal ischemia, 2) examine the sufficiency of existing hypotheses to explain experimental data - i.e. lethality in a model of major surgical insult and stress, and 3) produce behavior to potentially guide future experimental design - i.e. suggested sample points for a potential laboratory model of non-lethal transient intestinal ischemia. Furthermore, hypotheses were generated to explain certain discrepancies between the behaviors of the GMABM and biological experiments, and new investigatory avenues proposed to test those hypotheses.

Conclusions/Significance

Agent-based modeling can account for the spatio-temporal dynamics of an HPI, and, even when carried out with a relatively high degree of abstraction, can be useful in the investigation of system-level consequences of putative mechanisms operating at the individual agent level. We suggest that an integrated and iterative heuristic relationship between computational modeling and more traditional laboratory and clinical investigations, with a focus on identifying useful and sufficient degrees of abstraction, will enhance the efficiency and translational productivity of biomedical research.

A multifaceted approach to modeling the immune response in tuberculosis

Tuberculosis (TB) is a deadly infectious disease caused by Mycobacterium tuberculosis (Mtb). No available vaccine is reliable and, although treatment exists, approximately 2 million people still die each year. The hallmark of TB infection is the granuloma, a self-organizing structure of immune cells forming in the lung and lymph nodes in response to bacterial invasion. Protective immune mechanisms play a role in granuloma formation and maintenance; these act over different time/length scales (e.g., molecular, cellular, and tissue scales). The significance of specific immune factors in determining disease outcome is still poorly understood, despite incredible efforts to establish several animal systems to track infection progression and granuloma formation. Mathematical and computational modeling approaches have recently been applied to address open questions regarding host-pathogen interaction dynamics, including the immune response to Mtb infection and TB granuloma formation. This provides a unique opportunity to identify factors that are crucial to a successful outcome of infection in humans. These modeling tools not only offer an additional avenue for exploring immune dynamics at multiple biological scales but also complement and extend knowledge gained via experimental tools. We review recent modeling efforts in capturing the immune response to Mtb, emphasizing the importance of a multiorgan and multiscale approach that has tuneable resolution. Together with experimentation, systems biology has begun to unravel key factors driving granuloma formation and protective immune response in TB. WIREs Syst Biol Med 2011 3 479-489 DOI: 10.1002/wsbm.131

Multiscale computational modeling reveals a critical role for TNF-α receptor 1 dynamics in tuberculosis granuloma formation

Multiple immune factors control host responses to Mycobacterium tuberculosis infection, including the formation of granulomas, which are aggregates of immune cells whose function may reflect success or failure of the host to contain infection. One such factor is TNF-α. TNF-α has been experimentally characterized to have the following activities in M. tuberculosis infection: macrophage activation, apoptosis, and chemokine and cytokine production. Availability of TNF-α within a granuloma has been proposed to play a critical role in immunity to M. tuberculosis. However, in vivo measurement of a TNF-α concentration gradient and activities within a granuloma are not experimentally feasible. Further, processes that control TNF-α concentration and activities in a granuloma remain unknown. We developed a multiscale computational model that includes molecular, cellular, and tissue scale events that occur during granuloma formation and maintenance in lung. We use our model to identify processes that regulate TNF-α concentration and cellular behaviors and thus influence the outcome of infection within a granuloma. Our model predicts that TNF-αR1 internalization kinetics play a critical role in infection control within a granuloma, controlling whether there is clearance of bacteria, excessive inflammation, containment of bacteria within a stable granuloma, or uncontrolled growth of bacteria. Our results suggest that there is an interplay between TNF-α and bacterial levels in a granuloma that is controlled by the combined effects of both molecular and cellular scale processes. Finally, our model elucidates processes involved in immunity to M. tuberculosis that may be new targets for therapy.

Higher level of replication efficiency of 2009 (H1N1) pandemic influenza virus than those of seasonal and avian strains: kinetics from epithelial cell culture and computational modeling

The pathogenicity and transmission of influenza A viruses are likely determined in part by replication efficiency in human cells, which is the net effect of complex virus-host interactions. H5N1 avian, H1N1 seasonal, and H1N1 2009 pandemic influenza virus strains were compared by infecting human differentiated bronchial epithelial cells in air-liquid interface cultures at relatively low virus particle/cell ratios. Differential equation and computational models were used to characterize the in vitro kinetic behaviors of the three strains. The models were calibrated by fitting experimental data in order to estimate difficult-to-measure parameters. Both models found marked differences in the relative values of p, the virion production rate per cell, and R(0), an index of the spread of infection through the monolayer, with the values for the strains in the following rank order (from greatest to least): pandemic strain, followed by seasonal strain, followed by avian strain, as expected. In the differential equation model, which treats virus and cell populations as well mixed, R(0) and p varied proportionately for all 3 strains, consistent with a primary role for productivity. In the spatially explicit computational model, R(0) and p also varied proportionately except that R(0) derived for the pandemic strain was reduced, consistent with constrained viral spread imposed by multiple host defenses, including mucus and paracrine antiviral effects. This synergistic experimental-computational strategy provides relevant parameters for identifying and phenotyping potential pandemic strains.

Tuberculosis: global approaches to a global disease

Mycobacterium tuberculosis is a remarkably successful human pathogen. The interaction with the human host is complex and much remains unknown. Recent advances in systems biology have allowed the integration of data from humans and animal models into computational approaches. For example, mathematical models provide a platform for in silico manipulation of host-pathogen interactions to gain insight into this infection across temporal and biologic scales. Here, we review recent studies on global approaches toward identifying comprehensive responses of both host and bacillus during infection, and the potential for incorporation of these data into many types of useful computational systems. Systems biology approaches provide a unique opportunity to study interventions that may improve therapy and vaccines against this major killer.

Agent-based modeling of host-pathogen systems: The successes and challenges

Agent-based models have been employed to describe numerous processes in immunology. Simulations based on these types of models have been used to enhance our understanding of immunology and disease pathology. We review various agent-based models relevant to host-pathogen systems and discuss their contributions to our understanding of biological processes. We then point out some limitations and challenges of agent-based models and encourage efforts towards reproducibility and model validation.

Synergy between individual TNF-dependent functions determines granuloma performance for controlling Mycobacterium tuberculosis infection

Mycobacterium tuberculosis is one of the world's most deadly human pathogens; an integrated understanding of how it successfully survives in its host is crucial to developing new treatment strategies. One notable characteristic of infection with M. tuberculosis is the formation of granulomas, aggregates of immune cells whose structure and function may reflect success or failure of the host to contain infection. One central regulator of host responses to infection, including granuloma formation, is the pleiotropic cytokine TNF-alpha. Experimental work has characterized roles for TNF in macrophage activation; regulation of apoptosis; chemokine and cytokine production; and regulation of cellular recruitment via transendothelial migration. Separating the effects of these functions is presently difficult or impossible in vivo. To this end, we applied a computational model to understand specific roles of TNF in control of tuberculosis in a single granuloma. In the model, cells are represented as discrete entities on a spatial grid responding to environmental stimuli by following programmed rules determined from published experimental studies. Simulated granulomas emerge as a result of these rules. After confirming the importance of TNF in this model, we assessed the effects of individual TNF functions. The model predicts that multiple TNF activities contribute to control of infection within the granuloma, with macrophage activation as a key effector mechanism for controlling bacterial growth. Results suggest that bacterial numbers are a strong contributing factor to granuloma structure with TNF. Finally, TNF-dependent apoptosis may reduce inflammation at the cost of impairing mycobacterial clearance.

Functional group/guild modelling of inter-specific pathogen interactions: a potential tool for predicting the consequences of co-infection

Although co-infection is the norm in most human and animal populations, clinicians currently have no practical tool to assist them in choosing the best treatment strategy for such patients. Given the vast range of potential pathogens which may co-infect the host, obtaining such a practical tool may seem an intractable problem. In ecology the joint concepts of functional groups and guilds have been used to conceptually simplify complex ecosystems, in order to understand how their component parts interact and may be manipulated. Here we propose a mechanism by which to apply these concepts to pathogen co-infection systems. Further, we describe how these groups could be incorporated into a mathematical modelling framework which, after validation, could be used as a clinical tool to predict the outcome of any particular combination of pathogens co-infecting a host.

Computer immunology

This review describes a body of work on computational immune systems that behave analogously to the natural immune system. These artificial immune systems (AIS) simulate the behavior of the natural immune system and in some cases have been used to solve practical engineering problems such as computer security. AIS have several strengths that can complement wet lab immunology. It is easier to conduct simulation experiments and to vary experimental conditions, for example, to rule out hypotheses; it is easier to isolate a single mechanism to test hypotheses about how it functions; agent-based models of the immune system can integrate data from several different experiments into a single in silico experimental system.