An agent-based model of inflammation and fibrosis following particulate exposure in the lung

Elucidating the interaction between stretch and stiffness using an agent-based spring network model of progressive pulmonary fibrosis

Pulmonary fibrosis is a deadly disease that involves the dysregulation of fibroblasts and myofibroblasts, which are mechanosensitive. Previous computational models have succeeded in modeling stiffness-mediated fibroblasts behaviors; however, these models have neglected to consider stretch-mediated behaviors, especially stretch-sensitive channels and the stretch-mediated release of latent TGF-β. Here, we develop and explore an agent-based model and spring network model hybrid that is capable of recapitulating both stiffness and stretch. Using the model, we evaluate the role of mechanical signaling in homeostasis and disease progression during self-healing and fibrosis, respectively. We develop the model such that there is a fibrotic threshold near which the network tends towards instability and fibrosis or below which the network tends to heal. The healing response is due to the stretch signal, whereas the fibrotic response occurs when the stiffness signal overpowers the stretch signal, creating a positive feedback loop. We also find that by changing the proportional weights of the stretch and stiffness signals, we observe heterogeneity in pathological network structure similar to that seen in human IPF tissue. The system also shows emergent behavior and bifurcations: whether the network will heal or turn fibrotic depends on the initial network organization of the damage, clearly demonstrating structure's pivotal role in healing or fibrosis of the overall network. In summary, these results strongly suggest that the mechanical signaling present in the lungs combined with network effects contribute to both homeostasis and disease progression.

An agent-based modeling approach for lung fibrosis in response to COVID-19

The severity of the COVID-19 pandemic has created an emerging need to investigate the long-term effects of infection on patients. Many individuals are at risk of suffering pulmonary fibrosis due to the pathogenesis of lung injury and impairment in the healing mechanism. Fibroblasts are the central mediators of extracellular matrix (ECM) deposition during tissue regeneration, regulated by anti-inflammatory cytokines including transforming growth factor beta (TGF-β). The TGF-β-dependent accumulation of fibroblasts at the damaged site and excess fibrillar collagen deposition lead to fibrosis. We developed an open-source, multiscale tissue simulator to investigate the role of TGF-β sources in the progression of lung fibrosis after SARS-CoV-2 exposure, intracellular viral replication, infection of epithelial cells, and host immune response. Using the model, we predicted the dynamics of fibroblasts, TGF-β, and collagen deposition for 15 days post-infection in virtual lung tissue. Our results showed variation in collagen area fractions between 2% and 40% depending on the spatial behavior of the sources (stationary or mobile), the rate of activation of TGF-β, and the duration of TGF-β sources. We identified M2 macrophages as primary contributors to higher collagen area fraction. Our simulation results also predicted fibrotic outcomes even with lower collagen area fraction when spatially-localized latent TGF-β sources were active for longer times. We validated our model by comparing simulated dynamics for TGF-β, collagen area fraction, and macrophage cell population with independent experimental data from mouse models. Our results showed that partial removal of TGF-β sources changed the fibrotic patterns; in the presence of persistent TGF-β sources, partial removal of TGF-β from the ECM significantly increased collagen area fraction due to maintenance of chemotactic gradients driving fibroblast movement. The computational findings are consistent with independent experimental and clinical observations of collagen area fractions and cell population dynamics not used in developing the model. These critical insights into the activity of TGF-β sources may find applications in the current clinical trials targeting TGF-β for the resolution of lung fibrosis.

An agent-based modeling approach for lung fibrosis in response to COVID-19

The severity of the COVID-19 pandemic has created an emerging need to investigate the long-term effects of infection on patients. Many individuals are at risk of suffering pulmonary fibrosis due to the pathogenesis of lung injury and impairment in the healing mechanism. Fibroblasts are the central mediators of extracellular matrix (ECM) deposition during tissue regeneration, regulated by anti-inflammatory cytokines including transforming growth factor beta (TGF-β). The TGF-β-dependent accumulation of fibroblasts at the damaged site and excess fibrillar collagen deposition lead to fibrosis. We developed an open-source, multiscale tissue simulator to investigate the role of TGF-β sources in the progression of lung fibrosis after SARS-CoV-2 exposure, intracellular viral replication, infection of epithelial cells, and host immune response. Using the model, we predicted the dynamics of fibroblasts, TGF-β, and collagen deposition for 15 days post-infection in virtual lung tissue. Our results showed variation in collagen area fractions between 2% and 40% depending on the spatial behavior of the sources (stationary or mobile), the rate of activation of TGF-β, and the duration of TGF-β sources. We identified M2 macrophages as primary contributors to higher collagen area fraction. Our simulation results also predicted fibrotic outcomes even with lower collagen area fraction when spatially-localized latent TGF-β sources were active for longer times. We validated our model by comparing simulated dynamics for TGF-β, collagen area fraction, and macrophage cell population with independent experimental data from mouse models. Our results showed that partial removal of TGF-β sources changed the fibrotic patterns; in the presence of persistent TGF-β sources, partial removal of TGF-β from the ECM significantly increased collagen area fraction due to maintenance of chemotactic gradients driving fibroblast movement. The computational findings are consistent with independent experimental and clinical observations of collagen area fractions and cell population dynamics not used in developing the model. These critical insights into the activity of TGF-β sources may find applications in the current clinical trials targeting TGF-β for the resolution of lung fibrosis.

Author summary

COVID-19 survivors are at risk of lung fibrosis as a long-term effect. Lung fibrosis is the excess deposition of tissue materials in the lung that hinder gas exchange and can collapse the whole organ. We identified TGF-β as a critical regulator of fibrosis. We built a model to investigate the mechanisms of TGF-β sources in the process of fibrosis. Our results showed spatial behavior of sources (stationary or mobile) and their activity (activation rate of TGF-β, longer activation of sources) could lead to lung fibrosis. Current clinical trials for fibrosis that target TGF-β need to consider TGF-β sources' spatial properties and activity to develop better treatment strategies.

A 3D Agent-Based Model of Lung Fibrosis

Understanding the pathophysiology of lung fibrosis is of paramount importance to elaborate targeted and effective therapies. As it onsets, the randomly accumulating extracellular matrix (ECM) breaks the symmetry of the branching lung structure. Interestingly, similar pathways have been reported for both idiopathic pulmonary fibrosis and radiation-induced lung fibrosis (RILF). Individuals suffering from the disease, the worldwide incidence of which is growing, have poor prognosis and a short mean survival time. In this context, mathematical and computational models have the potential to shed light on key underlying pathological mechanisms, shorten the time needed for clinical trials, parallelize hypotheses testing, and improve personalized drug development. Agent-based modeling (ABM) has proven to be a reliable and versatile simulation tool, whose features make it a good candidate for recapitulating emergent behaviors in heterogeneous systems, such as those found at multiple scales in the human body. In this paper, we detail the implementation of a 3D agent-based model of lung fibrosis using a novel simulation platform, namely, BioDynaMo, and prove that it can qualitatively and quantitatively reproduce published results. Furthermore, we provide additional insights on late-fibrosis patterns through ECM density distribution histograms. The model recapitulates key intercellular mechanisms, while cell numbers and types are embodied by alveolar segments that act as agents and are spatially arranged by a custom algorithm. Finally, our model may hold potential for future applications in the context of lung disorders, ranging from RILF (by implementing radiation-induced cell damage mechanisms) to COVID-19 and inflammatory diseases (such as asthma or chronic obstructive pulmonary disease).

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.

Mathematical modeling of ventilator-induced lung inflammation

Despite the benefits of mechanical ventilators, prolonged or misuse of ventilators may lead to ventilation-associated/ventilation-induced lung injury (VILI). Lung insults, such as respiratory infections and lung injuries, can damage the pulmonary epithelium, with the most severe cases needing mechanical ventilation for effective breathing and survival. Damaged epithelial cells within the alveoli trigger a local immune response. A key immune cell is the macrophage, which can differentiate into a spectrum of phenotypes ranging from pro- to anti-inflammatory. To gain a greater understanding of the mechanisms of the immune response to VILI and post-ventilation outcomes, we developed a mathematical model of interactions between the immune system and site of damage while accounting for macrophage phenotype. Through Latin hypercube sampling we generated a collection of parameter sets that are associated with a numerical steady state. We then simulated ventilation-induced damage using these steady state values as the initial conditions in order to evaluate how baseline immune state and lung health affect outcomes. We used a variety of methods to analyze the resulting parameter sets, transients, and outcomes, including a random forest decision tree algorithm and parameter sensitivity with eFAST. Analysis shows that parameters and properties of transients related to epithelial repair and M1 activation are important factors. Using the results of this analysis, we hypothesized interventions and used these treatment strategies to modulate the response to ventilation for particular parameters sets.

A physiological model of the inflammatory-thermal-pain-cardiovascular interactions during an endotoxin challenge

KEY POINTS

Inflammation in response to bacterial endotoxin challenge impacts physiological functions, including cardiovascular, thermal and pain dynamics, although the mechanisms are poorly understood. We develop an innovative mathematical model incorporating interaction pathways between inflammation and physiological processes observed in response to an endotoxin challenge. We calibrate the model to individual data from 20 subjects in an experimental study of the human inflammatory and physiological responses to endotoxin, and we validate the model against human data from an independent study. Using the model to simulate patient responses to different treatment modalities reveals that a multimodal treatment combining several therapeutic strategies gives the best recovery outcome.

ABSTRACT

Uncontrolled, excessive production of pro-inflammatory mediators from immune cells and traumatized tissues can cause systemic inflammatory conditions such as sepsis, one of the ten leading causes of death in the USA, and one of the three leading causes of death in the intensive care unit. Understanding how inflammation affects physiological processes, including cardiovascular, thermal and pain dynamics, can improve a patient's chance of recovery after an inflammatory event caused by surgery or a severe infection. Although the effects of the autonomic response on the inflammatory system are well-known, knowledge about the reverse interaction is lacking. The present study develops a mathematical model analyzing the inflammatory system's interactions with thermal, pain and cardiovascular dynamics in response to a bacterial endotoxin challenge. We calibrate the model with individual data from an experimental study of the inflammatory and physiological responses to a one-time administration of endotoxin in 20 healthy young men and validate it against data from an independent endotoxin study. We use simulation to explore how various treatments help patients exposed to a sustained pathological input. The treatments explored include bacterial endotoxin adsorption, antipyretics and vasopressors, as well as combinations of these. Our findings suggest that the most favourable recovery outcome is achieved by a multimodal strategy, combining all three interventions to simultaneously remove endotoxin from the body and alleviate symptoms caused by the immune system as it fights the infection.

Mathematical modeling of ventilator-induced lung inflammation

Respiratory infections, such as the novel coronavirus (SARS-COV-2) and other lung injuries, damage the pulmonary epithelium. In the most severe cases this leads to acute respiratory distress syndrome (ARDS). Due to respiratory failure associated with ARDS, the clinical intervention is the use of mechanical ventilation. Despite the benefits of mechanical ventilators, prolonged or misuse of these ventilators may lead to ventilation-associated/ventilation-induced lung injury (VILI). Damage caused to epithelial cells within the alveoli can lead to various types of complications and increased mortality rates. A key component of the immune response is recruitment of macrophages, immune cells that differentiate into phenotypes with unique pro- and/or anti-inflammatory roles based on the surrounding environment. An imbalance in pro- and anti-inflammatory responses can have deleterious effects on the individual's health. To gain a greater understanding of the mechanisms of the immune response to VILI and post-ventilation outcomes, we develop a mathematical model of interactions between the immune system and site of damage while accounting for macrophage polarization. Through Latin hypercube sampling we generate a virtual cohort of patients with biologically feasible dynamics. We use a variety of methods to analyze the results, including a random forest decision tree algorithm and parameter sensitivity with eFAST. Analysis shows that parameters and properties of transients related to epithelial repair and M1 activation and de-activation best predicted outcome. Using this new information, we hypothesize inter-ventions and use these treatment strategies to modulate damage in select virtual cases.

Multi-scale models of lung fibrosis

The architectural complexity of the lung is crucial to its ability to function as an organ of gas exchange; the branching tree structure of the airways transforms the tracheal cross-section of only a few square centimeters to a blood-gas barrier with a surface area of tens of square meters and a thickness on the order of a micron or less. Connective tissue comprised largely of collagen and elastic fibers provides structural integrity for this intricate and delicate system. Homeostatic maintenance of this connective tissue, via a balance between catabolic and anabolic enzyme-driven processes, is crucial to life. Accordingly, when homeostasis is disrupted by the excessive production of connective tissue, lung function deteriorates rapidly with grave consequences leading to chronic lung conditions such as pulmonary fibrosis. Understanding how pulmonary fibrosis develops and alters the link between lung structure and function is crucial for diagnosis, prognosis, and therapy. Further information gained could help elaborate how the healing process breaks down leading to chronic disease. Our understanding of fibrotic disease is greatly aided by the intersection of wet lab studies and mathematical and computational modeling. In the present review we will discuss how multi-scale modeling has facilitated our understanding of pulmonary fibrotic disease as well as identified opportunities that remain open and have produced techniques that can be incorporated into this field by borrowing approaches from multi-scale models of fibrosis beyond the lung.

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.

Influence of Pathogens and Mechanical Stimuli in Inflammation

Inflammation is a process driven by underlying cell-cell communication and many other factors. In this study, a model of cell-cell communications was proposed to study factors driving the inflammation time course. Analyses of inflammations that are driven by the combined effects of strain (mechanical stimuli) and/or pathogens are considered in this paper. An agent-based model was employed to simulate inflammation where macrophages and fibroblasts influence each other through cell signaling cytokines that diffuse and spread in the tissue space. The communication network of macrophages and fibroblasts was then inferred and its network model (termed TE network) was generated and analyzed. The results suggest that factors driving inflammation time course can be discriminated by the characteristics of TE networks. Inflammation driven only by pathogens has certain TE network characteristics indicating slower and lower information exchange among cells. Multiple stimuli can help to maintain sufficient information exchange among cells, which is beneficial for inflammation resolution. The TE network captures the unfolding of the innate immune system over time, and the history of pathogens invasion. The resulting network leads to an improved understanding of the resilience of the system to future pathogen invasion.

Gynura Rhizoma containing pyrrolizidine alkaloids induces the hepatic sinusoidal obstruction syndrome in mice via upregulating fibrosis-related factors

Recently, hepatic sinusoidal obstruction syndrome (HSOS) caused by herbal preparations containing pyrrolizidine alkaloids (PAs), such as Gynura Rhizoma (Tusanqi), has gained global attention. However, the lack of a reliable and reproducible animal model has greatly hampered mechanistic studies. Therefore, we aimed to establish a reproducible HSOS mouse model and investigate the hepatotoxic mechanism. The model was established by intragastrical administration of Gynura Rhizoma extract, i.e., 1.0 g extract/kg per day (equal to 16.7 g crude drug/kg per day based on extraction rate and 49.1 mg PA/kg per day based on the total PA content in the extract determined) for 40 successive days. Then, the mice were sacrificed, and their blood samples and livers were collected for analyses. Using hematoxylin-eosin (HE) and Masson staining, scanning electron microscopy imaging, clinical biomarkers, and other assays, we showed that the HSOS was successfully induced in our mouse model. Furthermore, we detected the key factors involved in liver fibrosis in the mice, revealing significantly increased hydroxyproline concentration; elevated expression of α-smooth muscle actin (α-SMA) and fibrosis-related genes such as Collagen-1, Collagen-3, Mmp2, Mmp13, Timp1, Timp3, and Activin, upregulated Smad3 phosphorylation, and increased serum TGF-β levels. Moreover, pro-inflammatory cytokines, including Tnf-α, Il-1β, and Il-6, were also increased in the model. All these results demonstrate the key roles of the TGF-β-Smad3 and inflammatory signaling pathways in this Gynura Rhizoma-induced HSOS mouse model, suggesting that blockade of fibrosis and/or inflammation should be an effective treatment for HSOS.

Integrated lung tissue mechanics one piece at a time: Computational modeling across the scales of biology

The lung is a delicately balanced and highly integrated mechanical system. Lung tissue is continuously exposed to the environment via the air we breathe, making it susceptible to damage. As a consequence, respiratory diseases present a huge burden on society and their prevalence continues to rise. Emergent function is produced not only by the sum of the function of its individual components but also by the complex feedback and interactions occurring across the biological scales - from genes to proteins, cells, tissue and whole organ - and back again. Computational modeling provides the necessary framework for pulling apart and putting back together the pieces of the body and organ systems so that we can fully understand how they function in both health and disease. In this review, we discuss models of lung tissue mechanics spanning from the protein level (the extracellular matrix) through to the level of cells, tissue and whole organ, many of which have been developed in isolation. This is a vital step in the process but to understand the emergent behavior of the lung, we must work towards integrating these component parts and accounting for feedback across the scales, such as mechanotransduction. These interactions will be key to unlocking the mechanisms occurring in disease and in seeking new pharmacological targets and improving personalized healthcare.

Overproduction of TNF-α and lung structural remodelling due to chronic exposure to volcanogenic air pollution

Volcanogenic air pollution studies and their effects on the respiratory system are still outnumbered by studies regarding the effects of anthropogenic air pollution, representing an unknown risk to human population inhabiting volcanic areas worldwide (either eruptive or non-eruptive areas). This study was carried in the archipelago of the Azores- Portugal, in two areas with active volcanism (Village of Furnas and Village of Ribeira Quente) and a reference site (Rabo de Peixe). The hydrothermal volcanism of Furnas volcanic complex is responsible for the release of 1000 t d^-1 of CO₂, H₂S, the radioactive gas - radon, among others. Besides the gaseous emissions, particulate matter and metals (Hg, Cd, Zn, Al, Ni, etc.) are also released into the environment. We tested a hypothesis whether chronic exposure to volcanogenic air pollution causes lung structural remodelling, in the house mouse, Mus musculus, as a bioindicator species. Histopathological evaluations were performed to assess the amount of macrophages, mononuclear leukocyte infiltrate, pulmonary emphysema, and the production of pro-inflammatory cytokine TNF-α. Also, the percentage of collagen and elastin fibers was calculated. Mice chronically exposed to volcanogenic air pollution presented an increased score in the histopathological evaluations for the amount of macrophages, mononuclear leukocyte infiltrate, pulmonary emphysema and production of TNF-α; and also increased percentages of collagen and elastin. For the first time, we demonstrate that non-eruptive active volcanism has a high potential to cause lung structural remodelling. This study also highlights the Mus musculus as a useful bioindicator for future biomonitoring programs in these type of volcanic environments.

Risk Analysis Implications of Dose-Response Thresholds for NLRP3 Inflammasome-Mediated Diseases: Respirable Crystalline Silica and Lung Cancer as an Example

Chronic inflammation mediates an extraordinarily wide range of diseases. Recent progress in understanding intracellular inflammasome assembly, priming, activation, cytokine signaling, and interactions with mitochondrial reactive oxygen species, lysosome disruption, cell death, and prion-like polymerization and spread of inflammasomes among cells, has potentially profound implications for dose-response modeling. This article discusses mechanisms of exposure concentration and duration thresholds for NOD-like receptor protein 3 (NLRP3)-mediated inflammatory responses and develops a simple biomathematical model of the onset of exposure-related tissue-level chronic inflammation and resulting disease risks, focusing on respirable crystalline silica (RCS) and lung cancer risk as an example. An inflammation-mediated 2-stage clonal expansion model of RCS-induced lung cancer is proposed that explains why relatively low estimated concentrations of RCS (eg, <1 mg/m³) do not increase lung cancer risk and why even high occupational concentrations increase risk only modestly (typically relative risk <2). The model of chronic inflammation implies a dose-response threshold for excess cancer risk, in contrast to traditional linear-no-threshold assumptions. If this implication is correct, then concentrations of crystalline silica (or amphibole asbestos fibers, or other environmental challenges that act via the NLRP3 inflammasome) below the threshold do not cause chronic inflammation and resulting elevated risks of inflammation-mediated diseases.

A review of inflammatory mechanism in airway diseases

BACKGROUND

Inflammation in the lung is the body's natural response to injury. It acts to remove harmful stimuli such as pathogens, irritants, and damaged cells and initiate the healing process. Acute and chronic pulmonary inflammation are seen in different respiratory diseases such as; acute respiratory distress syndrome, chronic obstructive pulmonary disease (COPD), asthma, and cystic fibrosis (CF).

FINDINGS

In this review, we found that inflammatory response in COPD is determined by the activation of epithelial cells and macrophages in the respiratory tract. Epithelial cells and macrophages discharge transforming growth factor-β (TGF-β), which trigger fibroblast proliferation and tissue remodeling. Asthma leads to airway hyper-responsiveness, obstruction, mucus hyper-production, and airway-wall remodeling. Cytokines, allergens, chemokines, and infectious agents are the main stimuli that activate signaling pathways in epithelial cells in asthma. Mutation of the CF transmembrane conductance regulator (CFTR) gene results in CF. Mutations in CFTR influence the lung epithelial innate immune function that leads to exaggerated and ineffective airway inflammation that fails to abolish pulmonary pathogens. We present mechanistic computational models (based on ordinary differential equations, partial differential equations and agent-based models) that have been applied in studying the complex physiological and pathological mechanisms of chronic inflammation in different airway diseases.

CONCLUSION

The scope of the present review is to explore the inflammatory mechanism in airway diseases and highlight the influence of aging on airways' inflammation mechanism. The main goal of this review is to encourage research collaborations between experimentalist and modelers to promote our understanding of the physiological and pathological mechanisms that control inflammation in different airway diseases.

Multi-scale immunological and biomechanical model of emphysema progression

This work presents a multi-scale agent-based model of emphysema progression that includes both the slow action of the immune system and the fast action of force redistribution and fracture propagation of the biological tissue. The two scales are coupled because the immune response causes inflammation and adaptation, which affects the biomechanical parameters of the tissue such as his elasticity. During repeated inflammation and breathing cycles, the tissue weakens and breaks down. We found that macrophages lifespan and cytokynes diffusion ratio are the parameters that influence the outcome of the model the most.

Coupled Immunological and Biomechanical Model of Emphysema Progression

Chronic Obstructive Pulmonary Disease (COPD) is a disabling respiratory pathology, with a high prevalence and a significant economic and social cost. It is characterized by different clinical phenotypes with different risk profiles. Detecting the correct phenotype, especially for the emphysema subtype, and predicting the risk of major exacerbations are key elements in order to deliver more effective treatments. However, emphysema onset and progression are influenced by a complex interaction between the immune system and the mechanical properties of biological tissue. The former causes chronic inflammation and tissue remodeling. The latter influences the effective resistance or appropriate mechanical response of the lung tissue to repeated breathing cycles. In this work we present a multi-scale model of both aspects, coupling Finite Element (FE) and Agent Based (AB) techniques that we would like to use to predict the onset and progression of emphysema in patients. The AB part is based on existing biological models of inflammation and immunological response as a set of coupled non-linear differential equations. The FE part simulates the biomechanical effects of repeated strain on the biological tissue. We devise a strategy to couple the discrete biological model at the molecular /cellular level and the biomechanical finite element simulations at the tissue level. We tested our implementation on a public emphysema image database and found that it can indeed simulate the evolution of clinical image biomarkers during disease progression.

Evaluation of Ventilation-Induced Lung Inflammation Through Multi-Scale Simulations

Ventilation-induced lung injury is a common problem faced by patients with respiratory problems who require mechanical ventilation (MV). This injury may lead to a greater chance of developing or exacerbating the acute respiratory distress syndrome which further complicates the therapeutic use of MV. The chain of events begins with the MV initiating an immune response that leads to inflammation induced tissue material alteration (stiffening) and eventually the loss of lung resistance. It is clear from this sequence of events that the phenomenon of ventilation induced injury is multi-scale by nature and, hence, requires holistic analysis involving simulations and informatics. An effective approach to this problem is to break it down into several major physical models. Each physical model is developed separately and can be seen as a component in a larger system that comprises the scale of the problem being investigated. In this paper, a multi-scale system consisting of breathing mechanics, tissue deformation, and cellular mechanics models is developed to assess the immune response. To demonstrate the potential of the model, a fluid–solid model is employed for breathing mechanics, a plane-strain elasticity model is applied to assess tissue deformation, and a cellular automata (CA) model is developed to account for immune response. A case study of three lower airways is presented. The CA model shows that this increased the immune response by five times, which correlates with alteration in the tissue microstructure. This alteration in turn is reflected in the material constant value obtained in the tissue mechanics model. However, the changes in strain rates in the airways after inflammation (and hence, lung compliance) were not as significant as the rates of change in immune response. Finally, results from the fluid–solid model demonstrate its potential for airflow characterization caused by tissue deformation that could lead to disease identification.

Multiplicity of Mathematical Modeling Strategies to Search for Molecular and Cellular Insights into Bacteria Lung Infection

Even today two bacterial lung infections, namely pneumonia and tuberculosis, are among the 10 most frequent causes of death worldwide. These infections still lack effective treatments in many developing countries and in immunocompromised populations like infants, elderly people and transplanted patients. The interaction between bacteria and the host is a complex system of interlinked intercellular and the intracellular processes, enriched in regulatory structures like positive and negative feedback loops. Severe pathological condition can emerge when the immune system of the host fails to neutralize the infection. This failure can result in systemic spreading of pathogens or overwhelming immune response followed by a systemic inflammatory response. Mathematical modeling is a promising tool to dissect the complexity underlying pathogenesis of bacterial lung infection at the molecular, cellular and tissue levels, and also at the interfaces among levels. In this article, we introduce mathematical and computational modeling frameworks that can be used for investigating molecular and cellular mechanisms underlying bacterial lung infection. Then, we compile and discuss published results on the modeling of regulatory pathways and cell populations relevant for lung infection and inflammation. Finally, we discuss how to make use of this multiplicity of modeling approaches to open new avenues in the search of the molecular and cellular mechanisms underlying bacterial infection in the lung.

Modeling the dynamics of chromosomal alteration progression in cervical cancer: A computational model

Computational modeling has been applied to simulate the heterogeneity of cancer behavior. The development of Cervical Cancer (CC) is a process in which the cell acquires dynamic behavior from non-deleterious and deleterious mutations, exhibiting chromosomal alterations as a manifestation of this dynamic. To further determine the progression of chromosomal alterations in precursor lesions and CC, we introduce a computational model to study the dynamics of deleterious and non-deleterious mutations as an outcome of tumor progression. The analysis of chromosomal alterations mediated by our model reveals that multiple deleterious mutations are more frequent in precursor lesions than in CC. Cells with lethal deleterious mutations would be eliminated, which would mitigate cancer progression; on the other hand, cells with non-deleterious mutations would become dominant, which could predispose them to cancer progression. The study of somatic alterations through computer simulations of cancer progression provides a feasible pathway for insights into the transformation of cell mechanisms in humans. During cancer progression, tumors may acquire new phenotype traits, such as the ability to invade and metastasize or to become clinically important when they develop drug resistance. Non-deleterious chromosomal alterations contribute to this progression.

In Situ Visualization for 3D Agent-Based Vocal Fold Inflammation and Repair Simulation

A fast and insightful visualization is essential in modeling biological system behaviors and understanding underlying inter-cellular mechanisms. High fidelity models produce billions of data points per time step, making in situ visualization techniques extremely desirable as they mitigate I/O bottlenecks and provide computational steering capability. In this work, we present a novel high-performance scheme to couple in situ visualization with the simulation of the vocal fold inflammation and repair using little to no extra cost in execution time or computing resources. The visualization component is first optimized with an adaptive sampling scheme to accelerate the rendering process while maintaining the precision of the displayed visual results. Our software employs VirtualGL to perform visualization in situ. The scheme overlaps visualization and simulation, resulting in the optimal utilization of computing resources. This results in an in situ system biology simulation suite capable of remote simulation of 17 million biological cells and 1.2 billion chemical data points, remote visualization of the results, and delivery of visualized frames with aggregated statistics to remote clients in real-time.

Fibroblasts: Diverse Cells Critical to Biomaterials Integration

Fibroblasts are key participants in wound healing and inflammation, and are capable of driving the progression of tissue repair to fully functional tissue or pathologic scar, or fibrosis, depending on the specific mechanical and biochemical cues with which they are presented. Thus, understanding and modulating the fibroblastic response to implanted materials is paramount to achieving desirable outcomes, such as long-term implant function or tissue regeneration. However, fibroblasts are remarkably heterogeneous and can differ vastly in their contributions to regeneration and fibrosis. This heterogeneity exists between tissues and within tissues, down to the level of individual cells. This review will discuss the role of fibroblasts, the pitfalls of describing them as a collective, the specifics of their function, and potential future directions to better understand and organize their highly variable biology.

Simulation of Healing Threshold in Strain-Induced Inflammation Through a Discrete Informatics Model

Respiratory diseases such as asthma and acute respiratory distress syndrome as well as acute lung injury involve inflammation at the cellular level. The inflammation process is very complex and is characterized by the emergence of cytokines along with other changes in cellular processes. Due to the complexity of the various constituents that makes up the inflammation dynamics, it is necessary to develop models that can complement experiments to fully understand inflammatory diseases. In this study, we developed a discrete informatics model based on cellular automata (CA) approach to investigate the influence of elastic field (stretch/strain) on the dynamics of inflammation and account for probabilistic adaptation based on statistical interpretation of existing experimental data. Our simulation model investigated the effects of low, medium, and high strain conditions on inflammation dynamics. Results suggest that the model is able to indicate the threshold of innate healing of tissue as a response to strain experienced by the tissue. When strain is under the threshold, the tissue is still capable of adapting its structure to heal the damaged part. However, there exists a strain threshold where healing capability breaks down. The results obtained demonstrate that the developed discrete informatics based CA model is capable of modeling and giving insights into inflammation dynamics parameters under various mechanical strain/stretch environments.

Strain-induced inflammation in pulmonary alveolar tissue due to mechanical ventilation

Inflammation is the body's attempt at self-protection to remove harmful stimuli, including damaged cells, irritants, or pathogens and begin the healing process. In this study, strain-induced inflammation in pulmonary alveolar tissue under high tidal volume is investigated through a combination of an inflammation model and fluid structure interaction (FSI) analysis. A realistic three-dimensional organ model for alveolar sacs is built, and FSI is employed to evaluate strain distribution in alveolar tissue for different tidal volume (TV) values under the mechanical ventilation (MV) condition. The alveolar tissue is treated as a hyperelastic solid and provides the environment for the tissue constituents. The influence of different strain distributions resulting from different tidal volumes is investigated. It is observed that strain is highly distributed in the inlet area. In addition, strain versus time curves in different locations through the alveolar model reveals that middle layers in the alveolar region would undergo higher levels of strain during breathing under the MV condition. Three different types of strain distributions in the alveolar region from the FSI simulation are transferred to the CA model to study population dynamics of cell constituents under MV for different TVs; 200, 500 and 1000 mL, respectively. The CA model results suggests that strain distribution plays a significant role in population dynamics. An interplay between strain magnitude and distribution appears to influence healing capability. Results suggest that increasing TV leads to an exponential rise in tissue damage by inflammation.

Systems Medicine for Lung Diseases: Phenotypes and Precision Medicine in Cancer, Infection, and Allergy

Lung diseases cause an enormous socioeconomic burden. Four of them are among the ten most important causes of deaths worldwide: Pneumonia has the highest death toll of all infectious diseases, lung cancer kills the most people of all malignant proliferative disorders, chronic obstructive pulmonary disease (COPD) ranks third in mortality among the chronic noncommunicable diseases, and tuberculosis is still one of the most important chronic infectious diseases. Despite all efforts, for example, by the World Health Organization and clinical and experimental researchers, these diseases are still highly prevalent and harmful. This is in part due to the specific organization of tissue homeostasis, architecture, and immunity of the lung. Recently, several consortia have formed and aim to bring together clinical and molecular data from big cohorts of patients with lung diseases with novel experimental setups, biostatistics, bioinformatics, and mathematical modeling. This "systems medicine" concept will help to match the different disease modalities with adequate therapeutic and possibly preventive strategies for individual patients in the sense of precision medicine.

A Computational Model of Cellular Engraftment on Lung Scaffolds

The possibility that stem cells might be used to regenerate tissue is now being investigated for a variety of organs, but these investigations are still essentially exploratory and have few predictive tools available to guide experimentation. We propose, in this study, that the field of lung tissue regeneration might be better served by predictive tools that treat stem cells as agents that obey certain rules of behavior governed by both their phenotype and their environment. Sufficient knowledge of these rules of behavior would then, in principle, allow lung tissue development to be simulated computationally. Toward this end, we developed a simple agent-based computational model to simulate geographic patterns of cells seeded onto a lung scaffold. Comparison of the simulated patterns to those observed experimentally supports the hypothesis that mesenchymal stem cells proliferate preferentially toward the scaffold boundary, whereas alveolar epithelial cells do not. This demonstrates that a computational model of this type has the potential to assist in the discovery of rules of cellular behavior.

A computational model of invasive aspergillosis in the lung and the role of iron

Background

Invasive aspergillosis is a severe infection of immunocompromised hosts, caused by the inhalation of the spores of the ubiquitous environmental molds of the Aspergillus genus. The innate immune response in this infection entails a series of complex and inter-related interactions between multiple recruited and resident cell populations with each other and with the fungal cell; in particular, iron is critical for fungal growth.

Results

A computational model of invasive aspergillosis is presented here; the model can be used as a rational hypothesis-generating tool to investigate host responses to this infection. Using a combination of laboratory data and published literature, an in silico model of a section of lung tissue was generated that includes an alveolar duct, adjacent capillaries, and surrounding lung parenchyma. The three-dimensional agent-based model integrates temporal events in fungal cells, epithelial cells, monocytes, and neutrophils after inhalation of spores with cellular dynamics at the tissue level, comprising part of the innate immune response. Iron levels in the blood and tissue play a key role in the fungus’ ability to grow, and the model includes iron recruitment and consumption by the different types of cells included. Parameter sensitivity analysis suggests the model is robust with respect to unvalidated parameters, and thus is a viable tool for an in silico investigation of invasive aspergillosis.

Conclusions

Using laboratory data from a mouse model of invasive aspergillosis in the context of transient neutropenia as validation, the model predicted qualitatively similar time course changes in fungal burden, monocyte and neutrophil populations, and tissue iron levels. This model lays the groundwork for a multi-scale dynamic mathematical model of the immune response to Aspergillus species.

Electronic supplementary material

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

Computational modeling of cardiac fibroblasts and fibrosis

Altered fibroblast behavior can lead to pathologic changes in the heart such as arrhythmia, diastolic dysfunction, and systolic dysfunction. Computational models are increasingly used as a tool to identify potential mechanisms driving a phenotype or potential therapeutic targets against an unwanted phenotype. Here we review how computational models incorporating cardiac fibroblasts have clarified the role for these cells in electrical conduction and tissue remodeling in the heart. Models of fibroblast signaling networks have primarily focused on fibroblast cell lines or fibroblasts from other tissues rather than cardiac fibroblasts, specifically, but they are useful for understanding how fundamental signaling pathways control fibroblast phenotype. In the future, modeling cardiac fibroblast signaling, incorporating -omics and drug-interaction data into signaling network models, and utilizing multi-scale models will improve the ability of in silico studies to predict potential therapeutic targets against adverse cardiac fibroblast activity.

A computational model of unresolved allergic inflammation in chronic asthma

We have previously developed an agent-based computational model to demonstrate the feasibility of a novel hypothesis we term the inflammatory twitch. This hypothesis potentially explains the dynamics of the normal response to allergic inflammation in the lung (Pothen JJ, Poynter ME, Bates JH. J Immunol 190: 3510-3516, 2013) on the basis that antigenic stimulation sets in motion both the onset of inflammation and its subsequent resolution. The result is a self-limited inflammatory event that is similar in a formal sense to a skeletal muscle twitch. We hypothesize here that the chronic airway inflammation characteristic of asthma may represent the failure of the inflammatory twitch to resolve back to baseline. Our model provides a platform with which to perform virtual experiments aimed at investigating possible mechanisms leading to accentuation and/or prolongation of the inflammatory twitch. We used our model to determine how the inflammatory twitch is modified by knocking out certain cell types, interfering with cell activity, and altering cell lifetimes. Increasing the duration of activation of proinflammatory cells (considered to be chiefly neutrophils and eosinophils) markedly accentuated and prolonged the inflammatory twitch. This aberrant twitch behavior was largely abrogated by knocking out T-helper cells (simulating the effect of corticosteroids). The aberrant inflammatory twitch was also normalized by reducing the lifetime of the proinflammatory cells, suggesting that increasing apoptosis of these cells may be a therapeutic target in asthma.

A Multiscale Agent-Based in silico Model of Liver Fibrosis Progression

Chronic hepatic inflammation involves a complex interplay of inflammatory and mechanical influences, ultimately manifesting in a characteristic histopathology of liver fibrosis. We created an agent-based model (ABM) of liver tissue in order to computationally examine the consequence of liver inflammation. Our liver fibrosis ABM (LFABM) is comprised of literature-derived rules describing molecular and histopathological aspects of inflammation and fibrosis in a section of chemically injured liver. Hepatocytes are modeled as agents within hexagonal lobules. Injury triggers an inflammatory reaction, which leads to activation of local Kupffer cells and recruitment of monocytes from circulation. Portal fibroblasts and hepatic stellate cells are activated locally by the products of inflammation. The various agents in the simulation are regulated by above-threshold concentrations of pro- and anti-inflammatory cytokines and damage-associated molecular pattern molecules. The simulation progresses from chronic inflammation to collagen deposition, exhibiting periportal fibrosis followed by bridging fibrosis, and culminating in disruption of the regular lobular structure. The ABM exhibited key histopathological features observed in liver sections from rats treated with carbon tetrachloride (CCl₄). An in silico “tension test” for the hepatic lobules predicted an overall increase in tissue stiffness, in line with clinical elastography literature and published studies in CCl₄-treated rats. Therapy simulations suggested differential anti-fibrotic effects of neutralizing tumor necrosis factor alpha vs. enhancing M2 Kupffer cells. We conclude that a computational model of liver inflammation on a structural skeleton of physical forces can recapitulate key histopathological and macroscopic properties of CCl₄-injured liver. This multiscale approach linking molecular and chemomechanical stimuli enables a model that could be used to gain translationally relevant insights into liver fibrosis.

Agent-based modeling of the immune system: NetLogo, a promising framework

Several components that interact with each other to evolve a complex, and, in some cases, unexpected behavior, represents one of the main and fascinating features of the mammalian immune system. Agent-based modeling and cellular automata belong to a class of discrete mathematical approaches in which entities (agents) sense local information and undertake actions over time according to predefined rules. The strength of this approach is characterized by the appearance of a global behavior that emerges from interactions among agents. This behavior is unpredictable, as it does not follow linear rules. There are a lot of works that investigates the immune system with agent-based modeling and cellular automata. They have shown the ability to see clearly and intuitively into the nature of immunological processes. NetLogo is a multiagent programming language and modeling environment for simulating complex phenomena. It is designed for both research and education and is used across a wide range of disciplines and education levels. In this paper, we summarize NetLogo applications to immunology and, particularly, how this framework can help in the development and formulation of hypotheses that might drive further experimental investigations of disease mechanisms.

An agent-based model of the response to angioplasty and bare-metal stent deployment in an atherosclerotic blood vessel

PURPOSE

While animal models are widely used to investigate the development of restenosis in blood vessels following an intervention, computational models offer another means for investigating this phenomenon. A computational model of the response of a treated vessel would allow investigators to assess the effects of altering certain vessel- and stent-related variables. The authors aimed to develop a novel computational model of restenosis development following an angioplasty and bare-metal stent implantation in an atherosclerotic vessel using agent-based modeling techniques. The presented model is intended to demonstrate the body's response to the intervention and to explore how different vessel geometries or stent arrangements may affect restenosis development.

METHODS

The model was created on a two-dimensional grid space. It utilizes the post-procedural vessel lumen diameter and stent information as its input parameters. The simulation starting point of the model is an atherosclerotic vessel after an angioplasty and stent implantation procedure. The model subsequently generates the final lumen diameter, percent change in lumen cross-sectional area, time to lumen diameter stabilization, and local concentrations of inflammatory cytokines upon simulation completion. Simulation results were directly compared with the results from serial imaging studies and cytokine levels studies in atherosclerotic patients from the relevant literature.

RESULTS

The final lumen diameter results were all within one standard deviation of the mean lumen diameters reported in the comparison studies. The overlapping-stent simulations yielded results that matched published trends. The cytokine levels remained within the range of physiological levels throughout the simulations.

CONCLUSION

We developed a novel computational model that successfully simulated the development of restenosis in a blood vessel following an angioplasty and bare-metal stent deployment based on the characteristics of the vessel cross-section and stent. A further development of this model could ultimately be used as a predictive tool to depict patient outcomes and inform treatment options.

Mathematical modeling of apoptosis

Apoptosis is a form of programmed cell death, which is fundamental to all multicellular organisms. Deregulation of apoptosis leads to a number of severe diseases including cancer. Apoptosis is initiated either by extrinsic signals via stimulation of receptors at the cellular surface or intrinsic signals, such as DNA damage or growth factor withdrawal. Apoptosis has been extensively studied using systems biology which substantially contributed to the understanding of this death signaling network. This review gives an overview of mathematical models of apoptosis and the potential of systems biology to contribute to the development of novel therapies for cancer or other apoptosis-related diseases.

Chronic alcohol induces M2 polarization enhancing pulmonary disease caused by exposure to particulate air pollution

BACKGROUND

Chronic alcohol consumption causes persistent oxidative stress in the lung, leading to impaired alveolar macrophage (AM) function and impaired immune responses. AMs play a critical role in protecting the lung from particulate matter (PM) inhalation by removing particulates from the airway and secreting factors which mediate airway repair. We hypothesized AM dysfunction caused by chronic alcohol consumption increases the severity of injury caused by PM inhalation.

METHODS

Age- and sex-matched C57BL/6 mice were fed the Lieber-DeCarli liquid diet containing either alcohol or an isocaloric substitution (control diet) for 8 weeks. Mice from both diet groups were exposed to combustion-derived PM (CDPM) for the final 2 weeks. AM number, maturation, and polarization status were assessed by flow cytometry. Noninvasive and invasive strategies were used to assess pulmonary function and correlated with histomorphological assessments of airway structure and matrix deposition.

RESULTS

Co-exposure to alcohol and CDPM decreased AM number and maturation status (CD11c expression), while increasing markers of M2 activation (interleukin [IL]-4Rα, Ym1, Fizz1 expression, and IL-10 and transforming growth factor [TGF]-β production). Changes in AM function were accompanied by decreased airway compliance and increased elastance. Altered lung function was attributable to elevated collagen content localized to the small airways and loss of alveolar integrity. Intranasal administration of neutralizing antibody to TGF-β during the CDPM exposure period improved changes in airway compliance and elastance, while reducing collagen content caused by co-exposure.

CONCLUSIONS

Combustion-derived PM inhalation causes enhanced disease severity in the alcoholic lung by stimulating the release of latent TGF-β stores in AMs. The combinatorial effect of elevated TGF-β, M2 polarization of AMs, and increased oxidative stress impairs pulmonary function by increasing airway collagen content and compromising alveolar integrity.

Multiscale computational models of complex biological systems

Integration of data across spatial, temporal, and functional scales is a primary focus of biomedical engineering efforts. The advent of powerful computing platforms, coupled with quantitative data from high-throughput experimental methodologies, has allowed multiscale modeling to expand as a means to more comprehensively investigate biological phenomena in experimentally relevant ways. This review aims to highlight recently published multiscale models of biological systems, using their successes to propose the best practices for future model development. We demonstrate that coupling continuous and discrete systems best captures biological information across spatial scales by selecting modeling techniques that are suited to the task. Further, we suggest how to leverage these multiscale models to gain insight into biological systems using quantitative biomedical engineering methods to analyze data in nonintuitive ways. These topics are discussed with a focus on the future of the field, current challenges encountered, and opportunities yet to be realized.

The inflammatory twitch as a general strategy for controlling the host response

Allergic inflammation is a general host-defense mechanism for dealing with perceived foreign invaders. Although most effort has been directed toward understanding how this response gets turned on, how it gets turned off again when no longer needed is just as important to an organism's survival. We postulate that the control of the allergic inflammatory response is achieved via frequency modulation whereby a sequence of self-resolving events is repetitively invoked only so long as Ag is present. This leads to the notion of a unitary inflammatory event that we argue has formal similarity to the skeletal muscle twitch, albeit manifest over a much longer time scale. To test the plausibility of this hypothesis, we created an agent-based computational model of the allergic inflammatory response in the lungs. Continual stimulation of the model results in cycles of tissue damage and repair interspersed with periods of nonresponsiveness indicative of a refractory period. These findings are consistent with the inflammatory twitch hypothesis and the notion that the allergic inflammatory response is controlled via frequency modulation. We speculate that chronic inflammatory diseases may represent a failure of the inflammatory twitch to resolve toward baseline.

Functional complexity of the Leishmania granuloma and the potential of in silico modeling

In human and canine visceral leishmaniasis and in various experimental models of this disease, host resistance is strongly linked to efficient granuloma development. However, it is unknown exactly how the granuloma microenvironment executes an effective antileishmanial response. Recent studies, including using advanced imaging techniques, have improved our understanding of granuloma biology at the cellular level, highlighting heterogeneity in granuloma development and function, and hinting at complex cellular, temporal, and spatial dynamics. In this mini-review, we discuss the factors involved in the formation and function of Leishmania donovani-induced hepatic granulomas, as well as their importance in protecting against inflammation-associated tissue damage and the generation of immunity to rechallenge. Finally, we discuss the role that computational, agent-based models may play in answering outstanding questions within the field.

Translational applications of evaluating physiologic variability in human endotoxemia

Dysregulation of the inflammatory response is a critical component of many clinically challenging disorders such as sepsis. Inflammation is a biological process designed to lead to healing and recovery, ultimately restoring homeostasis; however, the failure to fully achieve those beneficial results can leave a patient in a dangerous persistent inflammatory state. One of the primary challenges in developing novel therapies in this area is that inflammation is comprised of a complex network of interacting pathways. Here, we discuss our approaches towards addressing this problem through computational systems biology, with a particular focus on how the presence of biological rhythms and the disruption of these rhythms in inflammation may be applied in a translational context. By leveraging the information content embedded in physiologic variability, ranging in scale from oscillations in autonomic activity driving short-term heart rate variability to circadian rhythms in immunomodulatory hormones, there is significant potential to gain insight into the underlying physiology.