{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,22]],"date-time":"2026-09-22T21:48:53Z","timestamp":1790113733837,"version":"4.0.1"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1014408","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2026,6,25]],"date-time":"2026-06-25T00:00:00Z","timestamp":1782345600000}}],"reference-count":41,"publisher":"Public Library of Science (PLoS)","issue":"6","license":[{"start":{"date-parts":[[2026,6,22]],"date-time":"2026-06-22T00:00:00Z","timestamp":1782086400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/publicdomain\/zero\/1.0\/"}],"funder":[{"DOI":"10.13039\/100008902","name":"Los Alamos National Laboratory","doi-asserted-by":"publisher","award":["LDRD grant 20200695E"],"award-info":[{"award-number":["LDRD grant 20200695E"]}],"id":[{"id":"10.13039\/100008902","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>\n                    Patients infected with life-threatening multi-drug resistant (MDR) bacteria have been treated with cocktails of bacteriophages. This is a complicated form of personalized medicine as the phages given to a patient have to be selected beforehand on the basis of their lytic capacity of the infecting bacteria. Because bacteria rapidly become resistant, the evolution of resistance to a diverse cocktail of phages is a complicated dynamical process, during which competing bacterial strains replace one another by accumulating several resistance mechanisms, each of which may involve a fitness cost. As a consequence, it is typically not known why a particular phage therapy succeeded or failed, and how one can optimize the composition of the cocktails to maximize the rate of success. To improve upon this, we extend an existing\n                    <jats:italic>in vivo<\/jats:italic>\n                    -calibrated mouse model into a novel mathematical model for the human situation, and include multiple phages infecting multiple bacterial strains, differing in their resistance to each of the phages. We adjust several parameter estimates of the bacterial model to the human situation, and use the model to describe a successful case of phage therapy involving several cocktails, each containing several phages. In the model, treatment success crucially depended on pretreatment resistance levels, and on the diversity and the timing of the cocktails. Once an appropriate cocktail is found, it is less important to further optimize the infection rates of the phages. Resistant bacterial strains expand rapidly when sensitive strains decline, and the higher the infectivity of the phages, the faster resistant strains expand. Because resistance evolves rapidly, it is best to provide a diverse set of phages right from the start of therapy, i.e., to hit hard and early, and create a high genetic barrier to bacterial resistance.\n                  <\/jats:p>","DOI":"10.1371\/journal.pcbi.1014408","type":"journal-article","created":{"date-parts":[[2026,6,22]],"date-time":"2026-06-22T17:52:17Z","timestamp":1782150737000},"page":"e1014408","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":1,"title":["Towards modeling phage therapy"],"prefix":"10.1371","volume":"22","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-2130-691X","authenticated-orcid":true,"given":"Rob J.","family":"de Boer","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Robert","family":"Schooley","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Alan S.","family":"Perelson","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"340","published-online":{"date-parts":[[2026,6,22]]},"reference":[{"key":"pcbi.1014408.ref001","doi-asserted-by":"crossref","first-page":"197","DOI":"10.1146\/annurev-med-080219-122208","article-title":"Phage Therapy for Antibiotic-Resistant Bacterial Infections","volume":"73","author":"GF Hatfull","year":"2022","journal-title":"Annu Rev Med"},{"issue":"6","key":"pcbi.1014408.ref002","doi-asserted-by":"crossref","first-page":"1434","DOI":"10.1038\/s41564-024-01705-x","article-title":"Personalized bacteriophage therapy outcomes for 100 consecutive cases: a multicentre, multinational, retrospective observational study","volume":"9","author":"J-P Pirnay","year":"2024","journal-title":"Nat Microbiol"},{"issue":"3","key":"pcbi.1014408.ref003","doi-asserted-by":"crossref","first-page":"391","DOI":"10.1038\/s41564-019-0666-4","article-title":"Treat phage like living antibiotics","volume":"5","author":"RT Schooley","year":"2020","journal-title":"Nat Microbiol"},{"issue":"10","key":"pcbi.1014408.ref004","doi-asserted-by":"crossref","DOI":"10.1128\/AAC.00954-17","article-title":"Development and Use of Personalized Bacteriophage-Based Therapeutic Cocktails To Treat a Patient with a Disseminated Resistant Acinetobacter baumannii Infection","volume":"61","author":"RT Schooley","year":"2017","journal-title":"Antimicrob Agents Chemother"},{"issue":"1","key":"pcbi.1014408.ref005","doi-asserted-by":"crossref","first-page":"60","DOI":"10.1093\/emph\/eoy005","article-title":"Phage treatment of an aortic graft infected with Pseudomonas aeruginosa","volume":"2018","author":"BK Chan","year":"2018","journal-title":"Evol Med Public Health"},{"issue":"5","key":"pcbi.1014408.ref006","doi-asserted-by":"crossref","first-page":"730","DOI":"10.1038\/s41591-019-0437-z","article-title":"Engineered bacteriophages for treatment of a patient with a disseminated drug-resistant Mycobacterium abscessus","volume":"25","author":"RM Dedrick","year":"2019","journal-title":"Nat Med"},{"issue":"1","key":"pcbi.1014408.ref007","doi-asserted-by":"crossref","first-page":"3776","DOI":"10.1038\/s41467-022-31455-5","article-title":"Comparative genomics of Acinetobacter baumannii and therapeutic bacteriophages from a patient undergoing phage therapy","volume":"13","author":"M Liu","year":"2022","journal-title":"Nat Commun"},{"issue":"2","key":"pcbi.1014408.ref008","doi-asserted-by":"crossref","first-page":"153","DOI":"10.1111\/j.1558-5646.1961.tb03139.x","article-title":"Conditions for the Existence of Bacteriophage","volume":"15","author":"A Campbell","year":"1961","journal-title":"Evolution"},{"issue":"977","key":"pcbi.1014408.ref009","doi-asserted-by":"crossref","first-page":"3","DOI":"10.1086\/283134","article-title":"Resource-Limited Growth, Competition, and Predation: A Model and Experimental Studies with Bacteria and Bacteriophage","volume":"111","author":"BR Levin","year":"1977","journal-title":"The American Naturalist"},{"issue":"1","key":"pcbi.1014408.ref010","doi-asserted-by":"crossref","DOI":"10.1016\/j.chom.2017.06.018","article-title":"Synergy between the Host Immune System and Bacteriophage Is Essential for Successful Phage Therapy against an Acute Respiratory Pathogen","volume":"22","author":"DR Roach","year":"2017","journal-title":"Cell Host Microbe"},{"issue":"6","key":"pcbi.1014408.ref011","doi-asserted-by":"crossref","first-page":"581","DOI":"10.1038\/s41587-019-0133-z","article-title":"Phage therapy\u2019s latest makeover","volume":"37","author":"C Schmidt","year":"2019","journal-title":"Nat Biotechnol"},{"key":"pcbi.1014408.ref012","doi-asserted-by":"crossref","first-page":"293","DOI":"10.4039\/Ent91293-5","article-title":"The components of predation as revealed by a study of small mammal predation of the European pine sawfly","volume":"91","author":"CS Holling","year":"1959","journal-title":"Canadian Entomologist"},{"issue":"5","key":"pcbi.1014408.ref013","doi-asserted-by":"crossref","first-page":"869","DOI":"10.1006\/bulm.2000.0181","article-title":"Modeling immune responses with handling time","volume":"62","author":"SS Pilyugin","year":"2000","journal-title":"Bull Math Biol"},{"issue":"5","key":"pcbi.1014408.ref014","doi-asserted-by":"crossref","DOI":"10.1152\/ajplung.00477.2004","article-title":"Sequential recruitment of neutrophils into lung and bronchoalveolar lavage fluid in LPS-induced acute lung injury","volume":"289","author":"J Reutershan","year":"2005","journal-title":"Am J Physiol Lung Cell Mol Physiol"},{"issue":"6","key":"pcbi.1014408.ref015","doi-asserted-by":"crossref","first-page":"989","DOI":"10.1172\/JCI104338","article-title":"Leukokinetic studies. IV. The total blood, circulating and marginal granulocyte pools and the granulocyte turnover rate in normal subjects","volume":"40","author":"JW Athens","year":"1961","journal-title":"J Clin Invest"},{"issue":"26","key":"pcbi.1014408.ref016","doi-asserted-by":"crossref","first-page":"3476","DOI":"10.1182\/blood-2016-07-727669","article-title":"Human CD62Ldim neutrophils identified as a separate subset by proteome profiling and in vivo pulse-chase labeling","volume":"129","author":"T Tak","year":"2017","journal-title":"Blood"},{"key":"pcbi.1014408.ref017","volume-title":"The Perfect Predator: A Scientist\u2019s Race to Save Her Husband from a Deadly Superbug","author":"S Strathdee","year":"2019"},{"issue":"7","key":"pcbi.1014408.ref018","doi-asserted-by":"crossref","first-page":"450","DOI":"10.1056\/NEJM199508173330710","article-title":"Time to hit HIV, early and hard","volume":"333","author":"DD Ho","year":"1995","journal-title":"N Engl J Med"},{"issue":"11","key":"pcbi.1014408.ref019","doi-asserted-by":"crossref","first-page":"1953","DOI":"10.1086\/430005","article-title":"Estimating HIV evolutionary pathways and the genetic barrier to drug resistance","volume":"191","author":"N Beerenwinkel","year":"2005","journal-title":"J Infect Dis"},{"issue":"12","key":"pcbi.1014408.ref020","doi-asserted-by":"crossref","first-page":"8289","DOI":"10.1073\/pnas.122244799","article-title":"A critical concentration of neutrophils is required for effective bacterial killing in suspension","volume":"99","author":"Y Li","year":"2002","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"5","key":"pcbi.1014408.ref021","doi-asserted-by":"crossref","first-page":"613","DOI":"10.1084\/jem.20040725","article-title":"Determination of the critical concentration of neutrophils required to block bacterial growth in tissues","volume":"200","author":"Y Li","year":"2004","journal-title":"J Exp Med"},{"issue":"8","key":"pcbi.1014408.ref022","doi-asserted-by":"crossref","first-page":"2776","DOI":"10.1172\/JCI63939","article-title":"How many neutrophils are enough (redux, redux)?","volume":"122","author":"SC Silverstein","year":"2012","journal-title":"J Clin Invest"},{"issue":"8","key":"pcbi.1014408.ref023","doi-asserted-by":"crossref","first-page":"3002","DOI":"10.1172\/JCI59832","article-title":"Evidence for bistable bacteria-neutrophil interaction and its clinical implications","volume":"122","author":"R Malka","year":"2012","journal-title":"J Clin Invest"},{"issue":"6","key":"pcbi.1014408.ref024","doi-asserted-by":"crossref","first-page":"2693","DOI":"10.1128\/AAC.01687-10","article-title":"Saturability of granulocyte kill of Pseudomonas aeruginosa in a murine model of pneumonia","volume":"55","author":"GL Drusano","year":"2011","journal-title":"Antimicrob Agents Chemother"},{"issue":"10","key":"pcbi.1014408.ref025","doi-asserted-by":"crossref","first-page":"4368","DOI":"10.1128\/AAC.00133-10","article-title":"Impact of burden on granulocyte clearance of bacteria in a mouse thigh infection model","volume":"54","author":"GL Drusano","year":"2010","journal-title":"Antimicrob Agents Chemother"},{"key":"pcbi.1014408.ref026","doi-asserted-by":"crossref","first-page":"241","DOI":"10.1016\/j.jtbi.2017.06.037","article-title":"Modeling the synergistic elimination of bacteria by phage and the innate immune system","volume":"429","author":"CYJ Leung","year":"2017","journal-title":"J Theor Biol"},{"issue":"3","key":"pcbi.1014408.ref027","doi-asserted-by":"crossref","first-page":"430","DOI":"10.1128\/mr.56.3.430-481.1992","article-title":"Bacteriophage lysis: mechanism and regulation","volume":"56","author":"R Young","year":"1992","journal-title":"Microbiol Rev"},{"issue":"1","key":"pcbi.1014408.ref028","doi-asserted-by":"crossref","first-page":"86","DOI":"10.1128\/MMBR.67.1.86-156.2003","article-title":"Bacteriophage T4 genome","volume":"67","author":"ES Miller","year":"2003","journal-title":"Microbiol Mol Biol Rev"},{"issue":"1","key":"pcbi.1014408.ref029","doi-asserted-by":"crossref","first-page":"43","DOI":"10.1007\/BF02458281","article-title":"Extending the quasi-steady state approximation by changing variables","volume":"58","author":"JA Borghans","year":"1996","journal-title":"Bull Math Biol"},{"issue":"4","key":"pcbi.1014408.ref030","doi-asserted-by":"crossref","first-page":"567","DOI":"10.1006\/jtbi.1995.0165","article-title":"Towards a general function describing T cell proliferation","volume":"175","author":"RJ De Boer","year":"1995","journal-title":"J Theor Biol"},{"issue":"2","key":"pcbi.1014408.ref031","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1012793","article-title":"Multi-strain phage induced clearance of bacterial infections","volume":"21","author":"J Marchi","year":"2025","journal-title":"PLoS Comput Biol"},{"issue":"12","key":"pcbi.1014408.ref032","first-page":"001110","article-title":"Functional diversity increases the efficacy of phage combinations","volume":"167","author":"RCT Wright","year":"2021","journal-title":"Microbiology (Reading)"},{"issue":"2","key":"pcbi.1014408.ref033","doi-asserted-by":"crossref","first-page":"113","DOI":"10.1038\/s41579-019-0278-2","article-title":"The pan-immune system of bacteria: antiviral defence as a community resource","volume":"18","author":"A Bernheim","year":"2020","journal-title":"Nat Rev Microbiol"},{"issue":"10","key":"pcbi.1014408.ref034","doi-asserted-by":"crossref","first-page":"104","DOI":"10.1007\/s11538-021-00936-x","article-title":"Atto-Foxes and Other Minutiae","volume":"83","author":"AC Fowler","year":"2021","journal-title":"Bull Math Biol"},{"key":"pcbi.1014408.ref035","article-title":"Stable coexistence and transport of lytic phage infections with migrating bacterial hosts","author":"J Marchi","year":"2025","journal-title":"bioRxiv"},{"issue":"16","key":"pcbi.1014408.ref036","doi-asserted-by":"crossref","first-page":"5646","DOI":"10.1128\/AEM.00757-12","article-title":"Effects of sequential and simultaneous applications of bacteriophages on populations of Pseudomonas aeruginosa in vitro and in wax moth larvae","volume":"78","author":"AR Hall","year":"2012","journal-title":"Appl Environ Microbiol"},{"issue":"1","key":"pcbi.1014408.ref037","doi-asserted-by":"crossref","first-page":"1057","DOI":"10.1186\/s13063-022-07047-5","article-title":"Safety and microbiological activity of phage therapy in persons with cystic fibrosis colonized with Pseudomonas aeruginosa: study protocol for a phase 1b\/2, multicenter, randomized, double-blind, placebo-controlled trial","volume":"23","author":"PD Tamma","year":"2022","journal-title":"Trials"},{"issue":"1","key":"pcbi.1014408.ref038","doi-asserted-by":"crossref","first-page":"5725","DOI":"10.1038\/s41467-025-61268-1","article-title":"Macrophage-induced reduction of bacteriophage density limits the efficacy of in vivo pulmonary phage therapy","volume":"16","author":"S Zborowsky","year":"2025","journal-title":"Nat Commun"},{"key":"pcbi.1014408.ref039","doi-asserted-by":"crossref","DOI":"10.1093\/cid\/ciad516","article-title":"Pharmacokinetic and Pharmacodynamic Obstacles for Phage Therapy From the Perspective of Clinical Practice","volume":"77","author":"R Mi\u0119dzybrodzki","year":"2023","journal-title":"Clin Infect Dis"},{"issue":"16","key":"pcbi.1014408.ref040","doi-asserted-by":"crossref","first-page":"7066","DOI":"10.1128\/JVI.00306-16","article-title":"Notwithstanding circumstantial alibis, cytotoxic T cells can be major killers of HIV-1-infected cells","volume":"90","author":"S Gadhamsetty","year":"2016","journal-title":"J Virol"},{"issue":"1","key":"pcbi.1014408.ref041","first-page":"1","article-title":"Inverse Modelling, Sensitivity and Monte Carlo Analysis in R Using Package FME","volume":"33","author":"K Soetaert","year":"2010","journal-title":"Journal of Statistical Software"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1014408","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2026,6,25]],"date-time":"2026-06-25T00:00:00Z","timestamp":1782345600000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1014408","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,6,25]],"date-time":"2026-06-25T17:51:17Z","timestamp":1782409877000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1014408"}},"subtitle":[],"editor":[{"given":"Roland R","family":"Regoes","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2026,6,22]]},"references-count":41,"journal-issue":{"issue":"6","published-online":{"date-parts":[[2026,6,22]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1014408","relation":{},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,6,22]]}}}