Combination of in vivo phage therapy data with in silico model highlights key parameters for treatment efficacy

The clinical (re)development of phage therapy to treat antibiotic resistant infections requires grasping specific biological properties of bacteriophages (phages) as antibacterial. However, identification of optimal dosing regimens is hampered by the poor understanding of phage-bacteria interactions in vivo. Here we developed a general strategy coupling in vitro and in vivo experiments with a mathematical model to characterize the interplay between phage and bacterial dynamics during pneumonia induced by a pathogenic strain of Escherichia coli. The model estimates some key parameters for phage therapeutic efficacy, in particular the impact of dose and route of administration on phage dynamics and the synergism of phage and the innate immune response on the bacterial clearance rate. Simulations predict a low impact of the intrinsic phage characteristics in agreement with the current semi-empirical choices of phages for compassionate treatments. Model-based approaches will foster the deployment of future phage therapy clinical trials.