Infections caused by multidrug resistant (MDR) pathogenic bacteria are a global health threat. Bacteriophages ("phage") are increasingly used as alternative or last-resort therapeutics to treat patients infected by MDR bacteria. However, the therapeutic outcomes of phage therapy may be limited by the emergence of phage resistance during treatment and/or by physical constraints that impede phage-bacteria interactions in vivo. In this work, we evaluate the role of lung spatial structure on the efficacy of phage therapy for Pseudomonas aeruginosa infections. To do so, we developed a spatially structured metapopulation network model based on the geometry of the bronchial tree, including host innate immune responses and the emergence of phage-resistant bacterial mutants. We model the ecological interactions between bacteria, phage, and the host innate immune system at the airway (node) level. The model predicts the synergistic elimination of a P. aeruginosa infection due to the combined effects of phage and neutrophils, given the sufficient innate immune activity and efficient phage-induced lysis. The metapopulation model simulations also predict that MDR bacteria are cleared faster at distal nodes of the bronchial tree. Notably, image analysis of lung tissue time series from wild-type and lymphocyte-depleted mice revealed a concordant, statistically significant pattern: infection intensity cleared in the bottom before the top of the lungs. Overall, the combined use of simulations and image analysis of in vivo experiments further supports the use of phage therapy for treating acute lung infections caused by P. aeruginosa, while highlighting potential limits to therapy in a spatially structured environment given impaired innate immune responses and/or inefficient phage-induced lysis.

Importance: Phage therapy is increasingly employed as a compassionate treatment for severe infections caused by multidrug-resistant (MDR) bacteria. However, the mixed outcomes observed in larger clinical studies highlight a gap in understanding when phage therapy succeeds or fails. Previous research from our team, using in vivo experiments and single-compartment mathematical models, demonstrated the synergistic clearance of acute P. aeruginosa pneumonia by phage and neutrophils despite the emergence of phage-resistant bacteria. In fact, the lung environment is highly structured, prompting the question of whether immunophage synergy explains the curative treatment of P. aeruginosa when incorporating realistic physical connectivity. To address this, we developed a metapopulation network model mimicking the lung branching structure to assess phage therapy efficacy for MDR P. aeruginosa pneumonia. The model predicts the synergistic elimination of P. aeruginosa by phage and neutrophils but emphasizes potential challenges in spatially structured environments, suggesting that higher innate immune levels may be required for successful bacterial clearance. Model simulations reveal a spatial pattern in pathogen clearance where P. aeruginosa are cleared faster at distal nodes of the bronchial tree than in primary nodes. Interestingly, image analysis of infected mice reveals a concordant and statistically significant pattern: infection intensity clears in the bottom before the top of the lungs. The combined use of modeling and image analysis supports the application of phage therapy for acute P. aeruginosa pneumonia while emphasizing potential challenges to curative success in spatially structured in vivo environments, including impaired innate immune responses and reduced phage efficacy.