A computational model of invasive aspergillosis in the lung and the role of iron
BMC Systems Biology volume 10, Article number: 34 (2016)
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.
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.
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.
Invasive aspergillosis represents a major and growing health problem in the U.S. and around the world. The growing population of immunocompromised patients, including those with haematologic malignancies, and stem cell- or solid organ-transplant recipients are at highest risk for this disease . In addition to conventionally immunosuppressed patients, other large populations are also at risk of this infection, including individuals with fibrocavitary tuberculosis in developing countries who develop chronic invasive aspergillosis as a secondary infection ; it is estimated that in 2007, 372,000 of the 7.7 million new cases of pulmonary tuberculosis world-wide also developed chronic pulmonary aspergillosis . The introduction of new antifungal drugs during the last decade, principally azole-based compounds capitalizing on new insights into the molecular structure of the fungal cell wall, has dramatically improved disease outcomes, but mortality rates remain approximately 30 % in recent surveys [1, 4]. In addition, increased resistance to these new drugs  raises the specter of a ‘perfect storm,’ as it has been called in , combining a rapidly growing patient population with a diminished repertoire of treatment options.
A substantial body of literature supports the critical role of iron homeostasis in Aspergillus biology. Aspergillus species adapt to iron-limited environments by activating a system of intracellular and secreted siderophores that scavenge iron from the environment and store it. In in vitro studies, Aspergillus siderophores remove iron from transferrin in human serum  and impair macrophage iron uptake ; conversely, neutrophil lactoferrin inhibits Aspergillus conidial growth by sequestering extracellular iron . In animal models, mutant Aspergillus species with defective siderophore systems are avirulent , and therapeutic iron chelation has an additive benefit to antifungal antibiotics . These mechanisms appear to be clinically important, since among immunocompromised stem cell transplant patients, clinically unsuspected iron overload is an independent risk factor with invasive aspergillosis [11, 12]. Taken together, these data suggest that the competition for iron is a key component of the pathogenesis of invasive aspergillosis.
The innate immune response to invasive aspergillosis is difficult to study. Interrogating dynamic cellular and molecular networks in a human host is, in most cases, impossible. In the study of the innate immune response to Aspergillus a number of in vitro and in vivo approaches have been used successfully. These include the in vitro interaction of Aspergillus with leukocytes and epithelial cells [13, 14]. In addition, animal models have been a valuable tool to investigate the complexities of cell-cell interactions and inflammatory pathways in a realistic system. These complementary approaches have led to recognition of neutrophils, macrophages, dendritic cells, and lung epithelial cells as key early players in host response to Aspergillus species [15, 16]. To date, the focus of the search for new therapeutics has been largely on fungal targets. But more recent promising efforts have looked to the host, in particular host immunity [17, 18]. However, a full exploration of the possibilities for anti-fungal therapeutics targeted at the host requires a better understanding of the innate host response. The complexity of the dynamic regulatory molecular networks and the multi-scale nature of the innate immune response strongly suggest taking a systems biology approach , as done in, e.g., [20, 21]. Here, we present a first step toward a multi-scale systems biology model of invasive aspergillosis in the lung, focused on the role of iron. In particular, we present the tissue level component of the model, validated with in vivo data from a mouse model of invasive aspergillosis.
Agent-based models (ABMs) are particularly well-suited for capturing the inherent heterogeneity of the immune system; an overview of such models is discussed in ; a more focused review of host-pathogen ABMs is presented in . In particular, the complexity of the lung physiology and its effect on dynamic interactions has been established , and specific interactions have been the focus of studies on intracellular regulatory networks  and the respiratory response to therapeutic interventions . A host-fungus interaction model incorporating data into clinically actionable therapeutic intervention in the case of invasive aspergillosis is presented in . The model presented here incorporates several parameter values and mechanistic behaviors from these models (see Additional files 1 and 2).
The literature indicates the importance of neutrophils and macrophages in the immune response to invasive aspergillosis. The critical importance of neutrophil involvement is detailed in , where chemotaxis is found to be the best strategy by which neutrophils find A. fumigatus conidia. This study also highlights the importance of pathogen distribution and spatial scale as critical factors, both of which have been incorporated into our model. Another recent ABM described the chemotactic recruitment of macrophages by epithelial cells [28, 29], focusing on the effectiveness of chemotaxis in the role of macrophage response to A. fumigatus. This model incorporates specific spatial structure as well as respiratory effects, and provides data on macrophage recruitment time and mechanisms of chemokine diffusion. Our model incorporates both macrophages and neutrophils, using chemokine diffusion and chemotaxis in a spatially heterogeneous domain to investigate the immune response to invasive aspergillosis over time. In addition, we introduce iron as a key factor in the survival of the fungus. ABMs have been used to study the lung in a variety of other contexts as well, including cancer , metastasis , fibrosis , and pneumococcal infection . Studies examining granuloma formation in tuberculosis infection have used ABMs of the lung to investigate chemokine diffusion  and pharmacokinetic/pharmacodynamic modeling of antibiotics .
There are many software packages for simulation of the immune system in various capacities. These include packages such as SIMMUNE , the Basic Immune Simulator , SIMISYS , and C-ImmSim [39–41]. These packages are general-purpose immune system simulators that, while effective and useful for their intended purposes, are ultimately not appropriate for our purposes.
Animals and in vivo procedures
Female wildtype C57Bl/6 mice were purchased from Jackson Laboratories (Bar Harbor, Maine) and maintained under pathogen-free conditions; experiments were performed in 6- to 10-week old animals. All animal experiments were approved by the Animal Care and Use Committee of University of Virginia.
We used a previously described animal model of invasive aspergillosis in neutropenic hosts [42–44]. Neutrophil depletion was achieved with a single i.p. injection of 80 μg of a monoclonal Ab (Gr-1, clone RB6-8C5) 1 day before an intratracheal challenge with A. fumigatus. We have previously reported that this protocol resulted in peripheral blood neutropenia (absolute circulating neutrophil count less than 50 cells/ μL) on days 1 and 3 after injection in both infected and uninfected mice, with a return of peripheral counts to pretreatment levels (>1000 cells/ μL) by day 5 [45, 46] and does not influence the number of other leukocyte subsets in the lung or spleen [44, 47].
Preparation and administration of A. fumigatus
A. fumigatus (strain 13073, American Type Culture Collection) conidia were collected in 0.1 % Tween in PBS from 7- to 14- day old cultures on Sabouraud’s dextrose agar plates, filtered through sterile gauze and counted under a hemacytometer. Fungal forms were administered intratracheally in inocula ranging from 2 to 5×106 in 30 μl saline per mouse.
Identification of leukocyte subsets
At pre-determined time points, animals were euthanized by CO2 asphyxiation, the pulmonary vasculature was perfused via the right ventricle with PBS containing 5mM EDTA, whole lungs were removed and single cell suspensions prepared as previously described [43–49]. The following antibodies were used to label cells for flow cytometry (from BD Biosciences, San Jose, CA, or eBiosciences, San Diego, CA): anti-CD11b-allophycocyanin-Cy7 (clone M1/70), anti-CD11c-PE-Cy7 (clone HL3), anti-CD45-peridinin chlorophyll protein (clone 30-F11), anti-Ly-6G-FITC (clone 1A8) and anti-Mac3-PE (clone M3/84). Samples were analyzed on a FACS Canto II instrument using Diva software (BD Biosciences). Neutrophils were identified as CD11b-hi Ly6G+ cells and recruited macrophages as CD11b+ CD11c- Mac3+ cells, as previously described [48, 49]. The absolute number of each leukocyte subset was determined as the product of the percentage of the cell type and the total number of cells in the sample, as determined using an automated cell counter (Countess, Invitrogen, Carlsbad, CA).
Since A. fumigatus grows as multicellular branching hyphae without forming distinct reproductive units in infected tissues, we used an assay for chitin, a carbohydrate component of hyphal wall that is absent from mammalian tissues and conidia, to quantify the burden of hyphae in infected lungs, as detailed previously . Organ chitin content has been validated as a readout of severity of infection in animal models of invasive aspergillosis by several groups [43, 51–54].
Measurement of lung iron
Lung iron content was measured as described in . Briefly, lungs were homogenized in 3 ml sterile water and 100 μL of each sample mixed with 100 μL of iron dissociation reagent (equal volumes 20 % trichloroacetic acid 2N HCl in distilled water and 4.5 % KMnO4, mixed immediately before use) in duplicate. Samples were incubated for 2 h at 60° C in a chemical fume hood and allowed to cool. Following incubation, 50 μL iron detection reagent (6.5 mM neucoproine, 6.5 mM ferrozine, 1 M L-ascorbic acid and 2.5 M ammonium acetate in distilled water) was added to each sample. After 30 min, samples were centrifuged at 10,000 g for ten minutes and plated on a 96-well plate with an iron standard for Atomic Absorption Spectrometry (Sigma-Aldrich) and optical density at 590 nm measured on a microplate reader.
Results and discussion
The simulation model
The agent-based model consists of a three-dimensional simulation of a 400 μ m×200 μ m×200 μ m section of lung tissue consisting of an alveolar duct, four adjacent capillaries, and surrounding lung parenchyma. The model was created in NetLogo , a popular platform for agent-based simulations. Figures 1, 2 and 3 provide still-frame snapshots from a three-dimensional dynamic simulation. A. fumigatus conidia enter at one end of the alveolar duct. As the simulation progresses, the conidia drift through the airway (see Fig. 1). The ABM simulates the epithelial clearance system via a parameter dictating the probability of conidia lodging in the epithelium. Left undisturbed, conidial spores enter a swelling phase prior to germination; subsequently, hyphal clusters begin to grow into the adjacent lung interstitium (see Fig. 2). In vivo, antimicrobial compounds in the airway surface fluid act as the first line of defense against conidia ; occasionally, conidia that are not cleared in this way can invade the interior of epithelial cells ; thus, the simulation also allows epithelial cells to internalize, damage, and kill conidia.
In the simulation model, epithelial cells act as the second line of defense against A. fumigatus. Epithelial cells recognize the presence of conidia and release inflammatory cytokines into the interstitial space in order to initiate an immune response. There are two immune cell types in the ABM: recruited monocytes/macrophages and neutrophils (see Fig. 3). We represent epithelial cells as releasing two types of chemotactic factors, for neutrophils and monocytes/macrophages, respectively. These chemotactic factors were tracked separately based on the literature [42, 59], and represent the aggregation of cytokines to which the respective immune cell types respond. The levels of these chemotactic factors are determined by the level of fungal burden, as measured by the number of conidia and hyphae. The chemotactic factors diffuse through the interstitial tissue, eventually reaching local capillaries. Once the local concentration for a particular immune cell rises above a threshold, they initiate recruitment of leukocytes adherent to the capillary vascular endothelium.
Recruited immune cells enter the interstitial space via the bloodstream. There, chemotaxis is simulated as the immune cells follow the gradient of the concentration of the chemotactic factor, a movement mechanism established in the literature [27, 28, 60]. At the source of the gradient, the immune cells may encounter fungal cells, whereupon such cells are attacked. Macrophages may internalize several fungal cells. Once internalized by a macrophage, cells are prohibited from escape and from germination. Over time, internalized fungal cells are damaged and ultimately destroyed. Damage inflicted by macrophages is limited by the number of internalized fungal cells, while neutrophil damage is limited by the number of available granules.
In the absence of sufficient chemotactic factors, immune cells move randomly throughout the tissue. Since the lifetime of neutrophils is between 24 and 48 h , neutrophils are represented as dying after this period. Macrophages leave the represented cross section when no conidial spores remain.
A key focus of the model is fungal acquisition of iron: the immune response induces hemorrhage, causing the tissue iron level to increase. The fungus is able to acquire iron both from the store of free iron and via a siderophore system. We model the level of available iron by having iron diffuse throughout the tissue. Once a certain level of iron is acquired by a fungal cell (determined by a parameter in the model), a new hyphal cell grows. Upon encountering hyphal cells, neutrophils sequester all of the iron in their immediate environment, preventing the fungus from acquiring it. If the immune system is not able to prevent the fungus from acquiring iron and consequently growing, invasive aspergillosis develops. The simulation continues in this manner until a pre-determined amount of time has passed. See Fig. 4 for a snapshot of the model interface, showing plots for various cell counts and model settings.
This description of the modeling of the immune response to A. fumigatus serves as an overview of the ABM and how the agents function. Many of the agent behaviors and parameter values are established in the literature; for a formal description of the model see Additional files 1 and 2. The ABM uses a combination of established behaviors, parameter settings, and the parameter analysis described below, building on previous work. Next we introduce results from parameter sensitivity analysis, which indicates the robustness of the model and offers insight into mechanisms of invasive aspergillosis in the lung.
We have chosen four quantities from the ABM to compare with laboratory time-course data: fungal burden, iron level throughout the tissue space, macrophage cell counts, and neutrophil cell counts. We simulate normal immune system conditions as well as neutropenia. Neutropenia is induced in the simulation model in a manner similar to the laboratory setting: presence of neutrophils in tissue is ablated for the first 48 h, and then resumes to normal levels over the remaining duration of the simulation. While certain aspects of the ABM simulation are fixed (e.g. number of initial conidia, lifespan of monocytes/neutrophils, killing rates of epithelial cells), there remains variation from run to run as many processes are pseudo-random. In order to obtain reliable results from simulation, we examined the quantities of interest over many runs, and used the average and the standard deviation to determine the typical behavior of the ABM with respect to these quantities. Data for all quantities except fungal growth were very reliable when averaged over 20 simulations (in both neutropenic and healthy patients); for fungal growth it was necessary to run 40 simulations in order to obtain reliable data. Figure 5 provides a summary of this data for both healthy and neutropenic patients. In this figure we see that the mean value from 40 simulations is very similar to the mean taken over 200 simulations, and the data indicate that there is no significant change expected beyond 200 runs. Thus all results were generated by using the mean over 40 runs at each setting, with error bars representing one standard deviation.
While many state variable values have been determined from the literature, there are nine parameters whose values had to be determined empirically. Parameter descriptions and baseline values are provided in Table 1. With one exception, each parameter was set to 10 %, 50 %, and 200 % of its baseline value (one at a time) in order to examine the robustness of model results with respect to the parameter value; simulations were run under both healthy and neutropenic conditions. Due to higher sensitivity, one parameter – the radius within which macrophages and neutrophils can detect fungal presence – was set to 12 μ m, 15 μ m, and 18 μ m. Full results for both normal and neutropenic conditions are provided in Additional file 3; an overview of the results follows.
Three of the nine parameters had no qualitative effect (and little quantitative effect) on model dynamics. The first of these, the cytokine production factor, is a multiplier that determines the unitless amount of macrophage- and neutrophil-specific cytokines that is produced by epithelial cells in response to fungal presence. The lack of effect is likely a result of the fact that macrophages and neutrophils follow cytokine gradients to determine movement; since chemotaxis is simulated based on relative levels of nearby cytokines, the actual multiplier is not critical. The second of these parameters is the proportion of nearby cytokines that are taken up by receptors on macrophages and neutrophils. This lack of a significant effect can be explained in a similar way: cytokine uptake does not have a sufficiently large effect on the relative nearby cytokine concentration to alter cell movement. The final parameter to have little effect is referred to as the recruitment threshold; this parameter determines both the minimal amount of cytokine that must be present in order for macrophages and neutrophils to move according to chemotaxis (if all nearby areas are below the threshold they move randomly) and the possibility for recruited macrophages and neutrophils to appear in the blood (new immune cells appear only if their respective cytokine levels are above the threshold). This parameter serves two roles. If there is no threshold for initiating chemotaxis, immune cells quickly become fixed in a single location and do not move at all. Additionally, there must be a minimal cytokine level in the bloodstream in order to initiate immune cell recruitment. This parameter’s lack of effect can possibly be due to it being relatively small with respect to the cytokine production factor: as long as cytokines are being produced, the minimal threshold requirement is more or less arbitrary, and serves only to keep the immune cells moving, and to ensure that they are not being continuously recruited without reference to what is happening in the system. Thus, in all three cases it is not surprising that the parameter scaling had little effect on model dynamics.
Four more parameters resulted in unsurprising quantitative (but not qualitative) differences. The first of these is the probability of conidial spores lodging in the epithelium – this can be thought of as the strength of the cilia in sweeping away the fungus. Figure 6 shows the effect of this parameter in an immunocompetent host simulation. As the probability of lodging increases, fungal growth increases. In turn, iron levels decrease as the fungus takes up more of the available iron. Higher fungus levels cause an increase in cytokine production, which leads to an increase in recruited macrophages. The second parameter that scales predictably is the proportion of available iron that is absorbed by fungal cells. Given a fixed amount of iron needed for growth, the lower this proportion is the less resultant fungal growth; at the same time, because there is less fungal growth there is an increase in systemic iron levels and a slightly lower macrophage count. However, when taking into account the variance of the data, the only significant difference is in fungal growth. The third parameter to scale quantitatively is the iron needed in order for fungus to grow. As this level is decreased, fungal growth increases; iron levels and macrophage counts remain approximately the same.
The fourth parameter in this category is the detection radius, which determines how close an immune cell needs to be in order to detect fungal presence. As noted above, model dynamics were more sensitive to this parameter than to any of the others. These results are highlighted in Fig. 7, which summarizes results in a neutropenic simulation. When the detection radius was set to be 12 μ m, fungal growth is rapid and hardly affected at all by the presence of immune cells. However, by changing the detection radius to just 15 μ m the dynamics were drastically altered: practically all of the fungus is removed by the immune cells. This is most likely an artifact of the way Netlogo handles distance – at 12 μ m, only the grid location of the immune cell is considered, whereas at 15 μ m each of the 26 neighbors (in three dimensions) are considered. Thus, in light of this, these results are not so surprising. From a detection radius of 15 μ m to 18 μ m, the results scale quantitatively in the expected manner: a higher detection radius results in less fungal growth and fewer macrophages, though iron levels remain fairly consistent.
The remaining two parameters subjected to sensitivity analysis are open to different interpretations, and offer interesting insights into the dynamics of the system. The first of these is the diffusion rate of iron and cytokines throughout the tissue. As the diffusion rate is lowered, the fungal growth rate is decreased, resulting in higher total iron levels in the system (as less iron is being consumed by the fungus); see Fig. 8 for results from neutropenic simulations. The diffusion rate does not have any significant effect on the macrophage count, as their presence is governed almost entirely by cytokine levels. Interestingly, at the highest diffusion rate (1.0), the total systemic iron is lower than when the diffusion rate is set to 0.25, but the fungal growth in each of these cases is the same. This effect is observed in both neutropenic and immunocompetent simulations. This may be because the macrophages are fighting off the fungus and thus controlling the fungal growth rate; in any event the result suggests that lowering iron diffusion might be an effective way to increase systemic iron levels without increasing fungal growth.
The last parameter subjected to sensitivity analysis is the maximum amount of iron that fungal cells can store; the data from a neutropenic simulation is presented in Fig. 9. As the amount of iron coming into the system is fixed (at a value of 1 unit per hour), the maximum iron level was set to 0.25, 1.25, and 5.0 (representing 10 %, 50 %, and 200 % of the fixed baseline value of 2.5 used in analysis of the other parameters). It is important to note that the amount of iron required for growth in each of these cases was fixed at the baseline value of 0.25. Thus when the maximum iron level is set to 0.25, the fungus is unable to grow; this allows systemic iron levels to rise, and since there is no fungus, fewer macrophages are recruited. These results suggest that even in neutropenic cases, fungal growth can be inhibited if the iron acquisition system is disabled.
Additional file 3 provides full data from the parameter analysis for both immunocompetent and neutropenic simulations; the results shown here summarize the most relevant findings. Data from neutropenic simulations are presented because in general, the effects are more pronounced in this setting. However, the same results hold for immunocompetent host simulations. It is worth noting that in all cases, neutrophil cell counts are not included as they are not significantly affected by different parameter values. In immunocompetent hosts, only the diffusion rate and the cytokine production factor affect neutrophil recruitment, and very slight effects are observed. In neutropenic simulations, no difference is observed. Due to the way in which neutropenia is simulated (as described above), it is not surprising that the parameters have little effect on neutrophil cell counts.
The parameter sensitivity analysis performed here is rudimentary and largely serves to examine the robustness of the system. Given that data are coming from stochastic simulations, a more traditional sensitivity analysis is computationally infeasible. Nevertheless, this analysis brings forth several interesting hypotheses. The dependence on iron for fungal growth cannot be considered a prediction of the model, as it is built into the code that dictates fungal behavior. However, the effect of iron and cytokine diffusion and iron storage on fungal growth indicate a possible means by which invasive aspergillosis can be avoided, even in the neutropenic case. It would be possible to investigate these predictions in future laboratory work. These results suggest that the model is robust with respect to the nine unvalidated parameters, indicating the strength and viability of using the simulation model to investigate invasive aspergillosis in the lung.
In vivo experiments were carried out for immunocompetent mice and those in whom neutropenia was induced, in order to study the effectiveness of the immune response under both conditions. Accordingly, simulations of the in silico model include data for both immunocompetent and neutropenic mice. In order to mimic in vivo conditions in neutropenic simulations, neutrophil levels are depleted gradually over the course of the first 24 h, n entirely for the following 48 h, and returned to full production over the course of the subsequent 24 h.
Chitin is a carbohydrate component of the hyphal wall that is absent from conidia and mammalian tissues. Chitin levels measured in vivo therefore indicate hyphal presence; as such, chitin levels are used as an indicator of total fungal burden. Figure 10 presents in vivo chitin levels and fungal cell counts from the agent-based model. Figures 11 and 12 present similar data (from the same simulations) for macrophage and neutrophil cell counts, respectively. In the data presented, a neutrophil-depleting antibody is given one day prior to the introduction of conidia; hence the time scales begin at day −1. Fungal spores are introduced in the ABM beginning at day 0. Figure 13 shows the qualitative similarity between experimental iron data and iron levels in the simulation. Parameter values for these figures were determined by a heuristic search method implemented directly in NetLogo; the values used to produce these figures are given in Table 1.
The immune response to fungal and other pathogens in the lung is multi-faceted and multi-scale. Thus, it lends itself well to a systems biology approach based on mathematical modeling. Our goal underlying the work presented here is to lay the groundwork for a multi-scale dynamic mathematical model of the innate immune response to A. fumigatus in the lung, focused on the “battle for iron” between host and pathogen. Iron acquisition has been shown to be an important virulence factor for the fungus, and several mechanisms of the immune response are designed to adversely affect fungal iron uptake. The model we constructed covers the tissue scale and incorporates recognition of the fungus by epithelial cells, production of chemotactic gradients and recruitment of macrophages and neutrophils, as well as changes in tissue-level iron concentrations. However, this is accomplished through rules that do not do full justice to the intracellular processes triggered by pathogen recognition. In a future multi-scale version of the model, these processes will be controlled through intracellular signaling and iron regulation mechanisms.
In any computational model, required approximations inevitably give rise to limitations. In the model presented here, one of these is the absence of volume exclusion; that is, that multiple cells may aggregate in a manner that is not physically feasible. While not indicated by preliminary simulation results, it is possible that cell crowding has an effect on the results of the model. The decision to not incorporate volume exclusion was motivated by the relatively low quantity of immune cells and the computational complexity required to model using volume exclusion. Incorporation of an increased variety of immune cell types would also increase the potential effectiveness of the model, as would incorporation of multiple cytokines rather than the cell-specific aggregates implemented here. A better understanding of how immune cells sequester iron, particularly in a quantitative sense, would also allow the model to be fine-tuned in a more realistic manner. The use of unit-less proportions to determine iron levels (and diffusion) act as another limitation, and one that could be addressed by biological experimentation. All of these limitations may be investigated in future versions of the model.
Parameter sensitivity analysis was conducted to ensure that the model is robust, and data was validated by fitting model parameters to data from a neutropenic mouse model. The computational model can capture the fungal burden qualitatively over time, using chitin levels as proxy, and the recruitment of macrophages, while reproducing leukocyte counts. Thus, this study can be considered as proof of the concept that a computational model of this kind can serve as a tool to study invasive aspergillosis. Some parameter values were determined from literature while others were determined experimentally; in this way, parameter values arose quite naturally from studies focusing on individual cells (for full model details see Additional files 1 and 2). In order to increase the realism of the model, we plan to add other cells (such as dendritic cells and natural killer cells); we then plan to endow all cells with intracellular networks that capture the signaling and metabolic processes involved. Such a multi-scale model can then be used to test hypotheses in silico and to explore therapeutic interventions.
Availability of data and materials
The data supporting the results of this article are included within the article.
Consent to publish
Animal studies were carried out in accordance with the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health and were approved by the Animal Studies Committee at University of Virginia School of Medicine.
Steinbach WJ, Marr KA, Anaissie EJ, Azie N, Quan SP, Meier-Kriesche HU, Apewokin S, Horn DL. Clinical epidemiology of 960 patients with invasive aspergillosis from the PATH Alliance registry. J Infect. 2012; 65(5):453–64.
Gangneux JP, Camus C, Philippe B. [Epidemiology of and risk factors for invasive aspergillosis in nonneutropenic patients]. Rev Mal Respir. 2008; 25(2):139–53.
Denning DW, Pleuvry A, Cole DC. Global burden of chronic pulmonary aspergillosis as a sequel to pulmonary tuberculosis. Bull World Health Organ. 2011; 89:864–72.
Herbrecht R, Denning DW, Patterson TF, Bennett JE, Greene RE, Oestmann JW, Kern WV, Marr KA, Ribaud P, Lortholary O, Sylvester R, Rubin RH, Wingard JR, Stark P, Durand C, Caillot D, Thiel E, Chandrasekar PH, Hodges MR, Schlamm HT, Troke PF, de Pauw B. Voriconazole versus amphotericin B for primary therapy of invasive aspergillosis. N Engl J Med. 2002; 347(6):408–15.
Hadrich I, Makni F, Neji S, Abbes S, Cheikhrouhou F, Trabelsi H, Sellami H, Ayadi A. Invasive aspergillosis: resistance to antifungal drugs. Mycopathologia. 2012; 174(2):131–41.
Hissen AH, Chow JM, Pinto LJ, Moore MM. Survival of Aspergillus fumigatus in serum involves removal of iron from transferrin: the role of siderophores. Infect Immun. 2004; 72(3):1402–8.
Seifert M, Nairz M, Schroll A, Schrettl M, Haas H, Weiss G. Effects of the Aspergillus fumigatus siderophore systems on the regulation of macrophage immune effector pathways and iron homeostasis. Immunobiology. 2008; 213(9-10):767–78.
Zarember KA, Sugui JA, Chang YC, Kwon-Chung KJ, Gallin JI. Human polymorphonuclear leukocytes inhibit Aspergillus fumigatus conidial growth by lactoferrin-mediated iron depletion. J Immunol. 2007; 178(10):6367–73.
Schrettl M, Bignell E, Kragl C, Joechl C, Rogers T, Arst HN, Haynes K, Haas H. Siderophore biosynthesis but not reductive iron assimilation is essential for Aspergillus fumigatus virulence. J Exp Med. 2004; 200(9):1213–9.
Ibrahim AS, Gebremariam T, French SW, Edwards JE, Spellberg B. The iron chelator deferasirox enhances liposomal amphotericin B efficacy in treating murine invasive pulmonary aspergillosis. J Antimicrob Chemother. 2010; 65(2):289–92.
Kontoyiannis DP, Chamilos G, Lewis RE, Giralt S, Cortes J, Raad II, Manning JT, Han X. Increased bone marrow iron stores is an independent risk factor for invasive aspergillosis in patients with high-risk hematologic malignancies and recipients of allogeneic hematopoietic stem cell transplantation. Cancer. 2007; 110(6):1303–6.
Altes A, Remacha AF, Sarda P, Sancho FJ, Sureda A, Martino R, Briones J, Brunet S, Canals C, Sierra J. Frequent severe liver iron overload after stem cell transplantation and its possible association with invasive aspergillosis. Bone Marrow Transplant. 2004; 34(6):505–9.
Hasenberg M, Behnsen J, Krappmann S, Brakhage A, Gunzer M. Phagocyte responses towards Aspergillus fumigatus. Int J Med Microbiol. 2011; 301(5):436–44.
Sun WK, Lu X, Li X, Sun QY, Su X, Song Y, Sun HM, Shi Y. Dectin-1 is inducible and plays a crucial role in Aspergillus-induced innate immune responses in human bronchial epithelial cells. Eur J Clin Microbiol Infect Dis. 2012; 31(10):2755–764.
Park SJ, Mehrad B. Innate immunity to Aspergillus species. Clin Microbiol Rev. 2009; 22(4):535–51.
Schrettl M, Ibrahim-Granet O, Droin S, Huerre M, Latge JP, Haas H. The crucial role of the Aspergillus fumigatus siderophore system in interaction with alveolar macrophages. Microbes Infect. 2010; 12(12-13):1035–41.
Carvalho A, Cunha C, Bistoni F, Romani L. Immunotherapy of aspergillosis. Clin Microbiol Infect. 2012; 18(2):120–5.
Carvalho A, Cunha C, Iannitti RG, Casagrande A, Bistoni F, Aversa F, Romani L. Host defense pathways against fungi: the basis for vaccines and immunotherapy. Front Microbiol. 2012; 3:176.
Horn F, Heinekamp T, Kniemeyer O, Pollmacher J, Valiante V, Brakhage AA. Systems biology of fungal infection. Front Microbiol. 2012; 3:108.
Albrecht D, Kniemeyer O, Mech F, Gunzer M, Brakhage A, Guthke R. On the way toward systems biology of Aspergillus fumigatus infection. Int J Med Microbiol. 2011; 301(5):453–9.
Linde J, Hortschansky P, Fazius E, Brakhage AA, Guthke R, Haas H. Regulatory interactions for iron homeostasis in Aspergillus fumigatus inferred by a Systems Biology approach. BMC Syst Biol. 2012; 6:6.
Chavali AK, Gianchandani EP, Tung KS, Lawrence MB, Peirce SM, Papin JA. Characterizing emergent properties of immunological systems with multi-cellular rule-based computational modeling. Trends Immunol. 2008; 29(12):589–99.
Bauer AL, Beauchemin CA, Perelson AS. Agent-based modeling of host-pathogen systems: The successes and challenges. Inf Sci (Ny). 2009; 179(10):1379–89.
Tawhai M, Clark A, Donovan G, Burrowes K. Computational modeling of airway and pulmonary vascular structure and function: development of a “lung physiome”. Crit Rev Biomed Eng. 2011; 39(4):319–36.
Lauzon AM, Bates JH, Donovan G, Tawhai M, Sneyd J, Sanderson MJ. A multi-scale approach to airway hyperresponsiveness: from molecule to organ. Front Physiol. 2012; 3:191.
Oliveira-Coelho A, Rodrigues F, Campos A, Lacerda JF, Carvalho A, Cunha C. Paving the way for predictive diagnostics and personalized treatment of invasive aspergillosis. Front Microbiol. 2015; 6:411.
Tokarski C, Hummert S, Mech F, Figge MT, Germerodt S, Schroeter A, Schuster S. Agent-based modeling approach of immune defense against spores of opportunistic human pathogenic fungi. Front Microbiol. 2012; 3:129.
Pollmächer J, Figge MT. Agent-based model of human alveoli predicts chemotactic signaling by epithelial cells during early Aspergillus fumigatus infection. PLoS ONE. 2014; 9(10):111630.
Pollmächer J, Figge MT. Deciphering chemokine properties by an agent-based model of Aspergillus fumigatus infection in human alveoli. Front Microbiol. 2015; 6(503). doi:http://dx.doi.org/10.3389/fmicb.2015.00503.
Wang Z, Zhang L, Sagotsky J, Deisboeck TS. Simulating non-small cell lung cancer with a multiscale agent-based model. Theor Biol Med Model. 2007; 4:50.
Pennisi M, Pappalardo F, Motta S. Agent based modeling of lung metastasis-immune system competition In: Andrews PS, Timmis J, Owens NDL, Aickelin U, Hart E, Hone A, Tyrrell AM, editors. Artificial Immune Systems: Lecture Notes in Computer Science. Berlin: Springer: 2009. p. 1–3.
Brown BN, Price IM, Toapanta FR, DeAlmeida DR, Wiley CA, Ross TM, Oury TD, Vodovotz Y. An agent-based model of inflammation and fibrosis following particulate exposure in the lung. Math Biosci. 2011; 231(2):186–96.
Smith AM, McCullers JA, Adler FR. Mathematical model of a three-stage innate immune response to a pneumococcal lung infection. J Theor Biol. 2011; 276(1):106–16.
Segovia-Juarez JL, Ganguli S, Kirschner D. Identifying control mechanisms of granuloma formation during m. tuberculosis infection using an agent-based model. J Theor Biol. 2004; 231(3):357–76. doi:http://dx.doi.org/10.1016/j.jtbi.2004.06.031.
Pienaar E, Cilfone NA, Lin PL, Dartois V, Mattila JT, Butler JR, Flynn JL, Kirschner DE, Linderman JJ. A computational tool integrating host immunity with antibiotic dynamics to study tuberculosis treatment. J Theor Biol. 2015; 367:166–79. doi:http://dx.doi.org/10.1016/j.jtbi.2014.11.021.
Meier-Schellersheim M, Xu X, Angermann B, Kunkel EJ, Jin T, Germain RN. (2006) Key Role of Local Regulation in Chemosensing Revealed by a New Molecular Interaction-Based Modeling Method. PLoS Comput Biol; 2(7):e82. doi:http://dx.doi.org/10.1371/journal.pcbi.0020082.
Folcik VA, An GC, Orosz CG. The Basic Immune Simulator: an agent-based model to study the interactions between innate and adaptive immunity. Theor Biol Med Model. 2007; 4:39.
Kalita JK, Chandrashekar K, Hans R, Selvam P. Computational modelling and simulation of the immune system. Int J Bioinform Res Appl. 2006; 2:63–88.
Castiglione F. C-ImmSim Simulator. Institute for Computing Applications, National Research Council (CNR) of Italy. 1995. http://www.iac.cnr.it/~filippo/c-immsim/. Accessed 18 Apr 2016.
Castiglione F, Bernaschi M, Succi S. Simulating the immune response on a distributed parallel computer. Int J Mod Phys C. 1997; 8(3):527–45.
Bernaschi M, Castiglione F. Design and implementation of an immune system simulator. Comput Biol Med. 2001; 31:303–31.
Mehrad B, Strieter RM, Moore TA, Tsai WC, Lira SA, Standiford TJ. CXC chemokine receptor-2 ligands are necessary components of neutrophil-mediated host defense in invasive pulmonary aspergillosis. J Immunol. 1999; 163(11):6086–94.
Park SJ, Hughes MA, Burdick M, Strieter RM, Mehrad B. Early NK cell-derived IFN-gamma is essential to host defense in neutropenic invasive aspergillosis. J Immunol. 2009; 182(7):4306–12.
Phadke AP, Akangire G, Park SJ, Lira SA, Mehrad B. The role of CC chemokine receptor 6 in host defense in a model of invasive pulmonary aspergillosis. Am J Respir Crit Care Med. 2007; 175(11):1165–72.
Park SJ, Wiekowski MT, Lira SA, Mehrad B. Neutrophils regulate airway responses in a model of fungal allergic airways disease. J Immunol. 2006; 176(4):2538–45.
Mehrad B, Park SJ, Akangire G, Standiford TJ, Wu T, Zhu J, Mohan C. The lupus-susceptibility locus, Sle3, mediates enhanced resistance to bacterial infections. J Immunol. 2006; 176(5):3233–9.
Park SJ, Burdick MD, Brix WK, Stoler MH, Askew DS, Strieter RM, Mehrad B. Neutropenia enhances lung dendritic cell recruitment in response to Aspergillus via a cytokine-to-chemokine amplification loop. J Immunol. 2010; 185(10):6190–7.
Chen L, Zhang Z, Barletta KE, Burdick MD, Mehrad B. Heterogeneity of lung mononuclear phagocytes during pneumonia: contribution of chemokine receptors. Am J Physiol Lung Cell Mol Physiol. 2013; 305(10):702–11.
Barletta KE, Cagnina RE, Wallace KL, Ramos SI, Mehrad B, Linden J. Leukocyte compartments in the mouse lung: distinguishing between marginated, interstitial, and alveolar cells in response to injury. J Immunol Methods. 2012; 375(1-2):100–10.
Lehmann PF, White LO. Chitin assay used to demonstrate renal localization and cortisone-enhanced growth of Aspergillus fumigatus mycelium in mice. Infect Immun. 1975; 12(5):987–92.
D’Angelo C, De Luca A, Zelante T, Bonifazi P, Moretti S, Giovannini G, Iannitti RG, Zagarella S, Bozza S, Campo S, Salvatori G, Romani L. Exogenous pentraxin 3 restores antifungal resistance and restrains inflammation in murine chronic granulomatous disease. J Immunol. 2009; 183(7):4609–18.
Hartigan AJ, Westwick J, Jarai G, Hogaboam CM. CCR7 deficiency on dendritic cells enhances fungal clearance in a murine model of pulmonary invasive aspergillosis. J Immunol. 2009; 183(8):5171–9.
Jaillon S, Peri G, Delneste Y, Fremaux I, Doni A, Moalli F, Garlanda C, Romani L, Gascan H, Bellocchio S, Bozza S, Cassatella MA, Jeannin P, Mantovani A. The humoral pattern recognition receptor PTX3 is stored in neutrophil granules and localizes in extracellular traps. J Exp Med. 2007; 204(4):793–804.
Stephens-Romero SD, Mednick AJ, Feldmesser M. The pathogenesis of fatal outcome in murine pulmonary aspergillosis depends on the neutrophil depletion strategy. Infect Immun. 2005; 73(1):114–25.
Riemer J, Hoepken HH, Czerwinska H, Robinson SR, Dringen R. Colorimetric ferrozine-based assay for the quantitation of iron in cultured cells. Anal Biochem. 2004; 331(2):370–5.
Wilensky U. NetLogo. Center for Connected Learning and Computer-Based Modeling. 2009. http://ccl.northwestern.edu/netlogo/.
Thompson AB, Robbins RA, Romberger DJ, Sisson JH, Spurzem JR, Teschler H, Rennard SI. Immunological functions of the pulmonary epithelium. Eur Respir J. 1995; 8(1):127–49.
Wasylnka JA, Moore MM. Uptake of Aspergillus fumigatus Conidia by phagocytic and nonphagocytic cells in vitro: quantitation using strains expressing green fluorescent protein. Infect Immun. 2002; 70(6):3156–63.
Espinosa V, Jhingran A, Dutta O, Kasahara S, Donnelly R, Du P, Rosenfeld J, Leiner I, Chen CC, Ron Y, Hohl TM, Rivera A. Inflammatory monocytes orchestrate innate antifungal immunity in the lung. PLoS Pathog. 2014; 10(2):1003940.
Dagenais TR, Keller NP. Pathogenesis of Aspergillus fumigatus in Invasive Aspergillosis. Clin Microbiol Rev. 2009; 22(3):447–65.
Stevens A, Lowe JS, Young B. Wheater’s Basic Histopathology: A Color Atlas and Text (Wheater’s Histology and Pathology), 4th edn. London: Churchill Livingstone; 2002.
The authors would like to acknowledge the help of Chris Gaskill and Eric Kernfeld, both of whom participated in a Research Experience for Undergraduates (REU) program in the summer of 2012 during which a preliminary version of the computer simulation was built. We also thank Shernita Lee and Madison Brandon for helpful discussions.
The authors declare that they have no competing interests.
MO contributed to the conception of the study, design of the agent-based model, implementation of the agent-based model, participated in analysis of agent-based model data and laboratory data, and helped to draft the manuscript. KRM and AMB contributed to the design and execution of the laboratory experimentation, and participated in its analysis. CL contributed to the conception and coordination of the study, design of the agent-based model, its statistical analysis, and analysis of biological relevance. BM contributed to the conception and coordination of the study, design and execution of the laboratory experimentation, supervised all laboratory procedures, participated in the design and statistical analysis of the agent-based model, and helped to draft the manuscript. RL contributed to the conception and coordination of the study, the design and analysis of the agent-based model, analysis of the laboratory data, and helped to draft the manuscript. All authors read and approved the final manuscript.
MO was supported in part by NSF Grant DMS-0931642 to the Mathematical Biosciences Institute. He was also supported partially by NSF Grant DMS-1062878 and US Department of Defense Grant Nr. W911NF0910538. KRM, AMB, and BM were supported in part by NIH Grants HL098526 and HL098329. RL was partially supported by NIH Grant Nr. 1R21AI101619-01. CL was supported by his affiliated institution. The funding institutions had no role in the design, collection, analysis, or interpretation of data; nor in the writing of the manuscript and the decision to submit the manuscript for publication.
Overview, Design Concepts, and Details (ODD) protocol for the agent-based model. A complete description of the agent-based model, including process ordering, state variable values, and algorithm pseudocode (in.pdf format). (PDF 145 kb)
State variables for the agent-based model. Descriptions and values for all state variable values, organized by cell type (with global variables indicated as such). N(μ,σ)indicates values taken from a normal distribution with mean μ and standard deviation σ (in.pdf format). (PDF 985 kb)
Experimental design and results of parameter sensitivity analysis. Parameters not validated by literature are altered to 10 %, 50 %, and 200 % of baseline values in order to determine the effect of each (one at a time). Relevant model dynamics are presented under both normal and neutropenic conditions (in.pdf format). (PDF 186 kb)
About this article
Cite this article
Oremland, M., Michels, K.R., Bettina, A.M. et al. A computational model of invasive aspergillosis in the lung and the role of iron. BMC Syst Biol 10, 34 (2016). https://doi.org/10.1186/s12918-016-0275-2
- Invasive aspergillosis
- Agent-based model