Aspergillus fumigatus is one of the most important human fungal pathogens, causing life-threatening diseases. Since humans inhale hundreds to thousands of fungal conidia every day, the lower respiratory tract is the primary site of infection. Current interaction networks of the innate immune response attribute fungal recognition and detection to alveolar macrophages, which are thought to be the first cells to get in contact with the fungus. At present, these networks are derived from in vitro or in situ assays, as the peculiar physiology of the human lung makes in vivo experiments, including imaging on the cell-level, hard to realize. We implemented a spatio-temporal agent-based model of a human alveolus in order to perform in silico experiments of a virtual infection scenario, for an alveolus infected with A. fumigatus under physiological conditions. The virtual analog captures the three-dimensional alveolar morphology consisting of the two major alveolar epithelial cell types and the pores of Kohn as well as the dynamic process of respiration. To the best of our knowledge this is the first agent-based model of a dynamic human alveolus in the presence of respiration. A key readout of our simulations is the first-passage-time of alveolar macrophages, which is the period of time that elapses until the first physical macrophage-conidium contact is established. We tested for random and chemotactic migration modes of alveolar macrophages and varied their corresponding parameter sets. The resulting first-passage-time distributions imply that randomly migrating macrophages fail to find the conidium before the start of germination, whereas guidance by chemotactic signals derived from the alveolar epithelial cell associated with the fungus enables a secure and successful discovery of the pathogen in time.

Elucidation of pathogenicity mechanisms of the most important human-pathogenic fungi, Aspergillus fumigatus and Candida albicans, has gained great interest in the light of the steadily increasing number of cases of invasive fungal infections. A key feature of these infections is the interaction of the different fungal morphotypes with epithelial and immune effector cells in the human host. Because of the high level of complexity, it is necessary to describe and understand invasive fungal infection by taking a systems biological approach, i.e., by a comprehensive quantitative analysis of the non-linear and selective interactions of a large number of functionally diverse, and frequently multifunctional, sets of elements, e.g., genes, proteins, metabolites, which produce coherent and emergent behaviors in time and space. The recent advances in systems biology will now make it possible to uncover the structure and dynamics of molecular and cellular cause-effect relationships within these pathogenic interactions. We review current efforts to integrate omics and image-based data of host-pathogen interactions into network and spatio-temporal models. The modeling will help to elucidate pathogenicity mechanisms and to identify diagnostic biomarkers and potential drug targets for therapy and could thus pave the way for novel intervention strategies based on novel antifungal drugs and cell therapy.

Opportunistic fungal pathogens can cause bloodstream infection and severe sepsis upon entering the blood stream of the host. The early immune response in human blood comprises the elimination of pathogens by antimicrobial peptides and innate immune cells, such as neutrophils or monocytes. Mathematical modeling is a predictive method to examine these complex processes and to quantify the dynamics of pathogen-host interactions. Since model parameters are often not directly accessible from experiment, their estimation is required by calibrating model predictions with experimental data. Depending on the complexity of the mathematical model, parameter estimation can be associated with excessively high computational costs in terms of run time and memory. We apply a strategy for reliable parameter estimation where different modeling approaches with increasing complexity are used that build on one another. This bottom-up modeling approach is applied to an experimental human whole-blood infection assay for Candida albicans. Aiming for the quantification of the relative impact of different routes of the immune response against this human-pathogenic fungus, we start from a non-spatial state-based model (SBM), because this level of model complexity allows estimating a priori unknown transition rates between various system states by the global optimization method simulated annealing. Building on the non-spatial SBM, an agent-based model (ABM) is implemented that incorporates the migration of interacting cells in three-dimensional space. The ABM takes advantage of estimated parameters from the non-spatial SBM, leading to a decreased dimensionality of the parameter space. This space can be scanned using a local optimization approach, i.e., least-squares error estimation based on an adaptive regular grid search, to predict cell migration parameters that are not accessible in experiment. In the future, spatio-temporal simulations of whole-blood samples may enable timely stratification of sepsis patients by distinguishing hyper-inflammatory from paralytic phases in immune dysregulation.