Optimization of sequential therapies to maximize extinction of resistant bacteria through collateral sensitivity
Pith reviewed 2026-05-18 11:05 UTC · model grok-4.3
The pith
Switching between two antibiotics at tuned intervals drives resistant bacteria to extinction by exploiting collateral sensitivity.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In a four-genotype stochastic birth-death model with two bacteriostatic antibiotics exhibiting strong reciprocal collateral sensitivity, the extinction probability of the bacterial population under subinhibitory concentrations varies nonlinearly with the switching period. This dependence manifests as stepwise increases in extinction probability that align with discrete switch events. Fast sequential therapies prove suboptimal because they fail to allow for the evolution of resistance, which serves as a crucial component in achieving eradication through collateral sensitivity. A geometric distribution framework provides accurate predictions for cumulative extinction probabilities, with the 0
What carries the argument
Four-genotype stochastic birth-death model that tracks wild-type, single-resistant, and double-resistant bacteria under alternating subinhibitory antibiotic concentrations linked by reciprocal collateral sensitivity.
Load-bearing premise
The two antibiotics exhibit strong reciprocal collateral sensitivity so that resistance to one substantially increases sensitivity to the other.
What would settle it
Laboratory experiments that measure extinction rates for bacterial populations with documented collateral sensitivity would falsify the claim if lengthening the switching period to match single-resistant mutant fixation times produced no stepwise rise in extinction probability.
Figures
read the original abstract
Antimicrobial resistance (AMR) threatens global health. A promising and underexplored strategy to tackle this problem is sequential therapies exploiting collateral sensitivity (CS), whereby resistance to one drug increases sensitivity to another. Here, we develop a four-genotype stochastic birth-death model with two bacteriostatic antibiotics to identify switching periods that maximize bacterial extinction under subinhibitory concentrations. We show that extinction probability depends nonlinearly on switching period, with stepwise increases aligned to discrete switch events: fast sequential therapies are suboptimal as they do not allow for the evolution of resistance, a key ingredient in these therapies. A geometric distribution framework accurately predicts cumulative extinction probabilities, where the per-switch extinction probability rises with switching period. We further derive a heuristic approximation for the extinction probability based on times to fixation of single-resistant mutants. Sensitivity analyses reveal that strong reciprocal CS is required for this strategy to work, and we explore how increasing antibiotic doses and higher mutation rates modulate extinction in a nonmonotonic manner. Finally, we discuss how longer therapies maximize extinction but also cause higher resistance, leading to a Pareto front of optimal switching periods. Our results provide quantitative design principles for in vitro and clinical sequential antibiotic therapies, underscoring the potential of CS-guided regimens to suppress resistance evolution and eradicate infections.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a four-genotype stochastic birth-death model for bacterial populations under sequential therapies with two bacteriostatic antibiotics exhibiting collateral sensitivity. It claims that extinction probability depends nonlinearly on switching period, with stepwise increases aligned to discrete switch events; fast sequential therapies are suboptimal because they do not allow resistance evolution. A geometric distribution framework predicts cumulative extinction probabilities (with per-switch extinction rising with period), and a heuristic approximation is derived from times to fixation of single-resistant mutants. Sensitivity analyses show that strong reciprocal collateral sensitivity is required for high extinction probabilities, while increasing doses and mutation rates modulate extinction nonmonotonically; longer therapies maximize extinction but increase resistance, yielding a Pareto front of optimal periods.
Significance. If the central claims hold, the work supplies quantitative design principles for collateral-sensitivity-guided sequential therapies that could inform in vitro and clinical protocols. The largely parameter-free qualitative structure, the geometric and heuristic approximations, and the explicit Pareto-front trade-off between extinction and resistance are strengths. The stochastic modeling and sensitivity analyses add rigor, though the dependence on strong reciprocal CS narrows the immediate translational scope.
major comments (1)
- [Sensitivity analyses] Sensitivity analyses section: the headline result (nonlinear rise in extinction probability with switching period and suboptimality of fast switching) is shown to collapse when reciprocal collateral sensitivity is weakened. Because this assumption is load-bearing for the design principle and the manuscript provides no discussion or citations on the empirical prevalence of strong reciprocal CS in clinical or environmental isolates, the applicability of the proposed switching-period optima remains uncertain.
minor comments (2)
- [Abstract] The abstract states that the geometric framework 'accurately predicts' cumulative extinction probabilities; a brief quantitative comparison (e.g., mean absolute error or R^{2} against simulations) in the main text would strengthen this claim.
- [Methods] Notation for the four genotypes and the fitness parameters under each drug could be introduced with a single table or diagram early in the Methods to improve readability.
Simulated Author's Rebuttal
We thank the referee for their constructive and insightful review of our manuscript. We address the major comment below in a point-by-point manner.
read point-by-point responses
-
Referee: [Sensitivity analyses] Sensitivity analyses section: the headline result (nonlinear rise in extinction probability with switching period and suboptimality of fast switching) is shown to collapse when reciprocal collateral sensitivity is weakened. Because this assumption is load-bearing for the design principle and the manuscript provides no discussion or citations on the empirical prevalence of strong reciprocal CS in clinical or environmental isolates, the applicability of the proposed switching-period optima remains uncertain.
Authors: We appreciate the referee drawing attention to this aspect of our sensitivity analyses. The manuscript already states that strong reciprocal collateral sensitivity is required for the strategy to produce high extinction probabilities and that the headline nonlinear dependence on switching period collapses under weaker CS. This is presented as a central result rather than an afterthought. We agree that the current version lacks explicit discussion of how commonly strong reciprocal CS occurs in clinical or environmental isolates, which limits immediate claims about applicability. In the revised manuscript we will add a dedicated paragraph in the Discussion that reviews published empirical examples of strong reciprocal CS (e.g., in Pseudomonas aeruginosa and other opportunistic pathogens under beta-lactam/aminoglycoside or fluoroquinolone pairings), cites the relevant experimental literature, and explicitly qualifies the translational scope of the predicted switching-period optima pending further strain-specific validation. revision: yes
Circularity Check
No significant circularity; derivation self-contained from model dynamics
full rationale
The paper constructs a four-genotype stochastic birth-death process with explicit fitness parameters encoding collateral sensitivity, then computes extinction probabilities via direct simulation and standard approximations (geometric distribution over per-switch events, heuristic based on fixation times of single-resistant mutants). These steps follow from solving or approximating the underlying Markov process rather than redefining inputs as outputs. No load-bearing self-citations, no parameters fitted to the target extinction probability and then relabeled as predictions, and no ansatz or uniqueness claim imported from prior author work. Sensitivity analyses independently vary the CS strength, confirming the result depends on an explicit modeling assumption rather than embedding it tautologically. The chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Strong reciprocal collateral sensitivity exists between the two bacteriostatic antibiotics.
- domain assumption Bacterial population dynamics can be captured by a four-genotype stochastic birth-death process.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
four-genotype stochastic birth-death model with two bacteriostatic antibiotics... birth rates β1,B = k_CS k_B β
-
IndisputableMonolith/Foundation/ArithmeticFromLogic.leanLogicNat induction and orbit embedding unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
extinction probability depends nonlinearly on switching period, with stepwise increases aligned to discrete switch events
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
In this scenario, rings, rather than stars, are the dominant structure
fit naturally within the same framework. CS rela- tionships may vary during treatment [45–47]; such time- dependence can be incorporated by allowing parameters to evolve between switches and calibrating them from data. Empirical evaluation can proceedin vitrowith- out altering the theoretical structure. A suitable test system is the ciprofloxacin–tobramyc...
-
[2]
P. Catal´ an, E. Wood, J. M. Blair, I. Gudelj, J. R. Iredell, and R. E. Beardmore, Seeking patterns of antibiotic re- sistance in atlas, an open, raw mic database with patient metadata, Nature Communications13, 2917 (2022)
work page 2022
-
[3]
M. Naghavi, S. E. Vollset, K. S. Ikuta, L. R. Swetschin- ski, A. P. Gray, E. E. Wool, G. R. Aguilar, T. Mestrovic, G. Smith, C. Han,et al., Global burden of bacterial an- timicrobial resistance 1990–2021: a systematic analysis with forecasts to 2050, The Lancet404, 1199 (2024)
work page 1990
-
[4]
M. Miethke, M. Pieroni, T. Weber, M. Br¨ onstrup, P. Hammann, L. Halby, P. B. Arimondo, P. Glaser, B. Aigle, H. B. Bode,et al., Towards the sustainable dis- covery and development of new antibiotics, Nature Re- views Chemistry5, 726 (2021)
work page 2021
-
[5]
G. Liu, D. B. Catacutan, K. Rathod, K. Swanson, W. Jin, J. C. Mohammed, A. Chiappino-Pepe, S. A. Syed, M. Fragis, K. Rachwalski,et al., Deep learning- guided discovery of an antibiotic targeting acinetobacter baumannii, Nature Chemical Biology19, 1342 (2023)
work page 2023
-
[6]
A. Martins, F. Jud´ ak, Z. Farkas, P. Szili, G. Gr´ ezal, B. Cs¨ org˝ o, M. S. Czikkely, E. Maharramov, L. Daruka, R. Spohn,et al., Antibiotic candidates for gram-positive bacterial infections induce multidrug resistance, Science Translational Medicine17, eadl2103 (2025)
work page 2025
- [7]
-
[8]
M. Baym, L. K. Stone, and R. Kishony, Multidrug evolu- tionary strategies to reverse antibiotic resistance, Science 351, aad3292 (2016)
work page 2016
-
[9]
M. Tyers and G. D. Wright, Drug combinations: a strat- egy to extend the life of antibiotics in the 21st century, Nature Reviews Microbiology17, 141 (2019)
work page 2019
-
[10]
R. Roemhild, T. Bollenbach, and D. I. Andersson, The physiology and genetics of bacterial responses to antibi- otic combinations, Nature Reviews Microbiology20, 478 (2022)
work page 2022
-
[11]
C. E. Rosenkilde, C. Munck, A. Porse, M. Linkevicius, D. I. Andersson, and M. O. Sommer, Collateral sensi- tivity constrains resistance evolution of the ctx-m-15β- lactamase, Nature Communications10, 618 (2019)
work page 2019
-
[12]
B. Kavˇ ciˇ c, G. Tkaˇ cik, and T. Bollenbach, Mechanisms of drug interactions between translation-inhibiting antibi- otics, Nature Communications11, 4013 (2020)
work page 2020
-
[13]
P. D. Tamma, S. E. Cosgrove, and L. L. Maragakis, Com- bination therapy for treatment of infections with gram- negative bacteria, Clinical Microbiology Reviews25, 450 (2012)
work page 2012
-
[14]
R. Pena-Miller, D. Laehnemann, G. Jansen, A. Fuentes- Hernandez, P. Rosenstiel, H. Schulenburg, and R. Beard- more, When the most potent combination of antibiotics selects for the greatest bacterial load: the smile-frown transition, PLoS Biology11, e1001540 (2013)
work page 2013
-
[15]
M. Vestergaard, W. Paulander, R. L. Marvig, J. Clasen, N. Jochumsen, S. Molin, L. Jelsbak, H. Ingmer, and A. Folkesson, Antibiotic combination therapy can se- lect for broad-spectrum multidrug resistance in pseu- domonas aeruginosa, International Journal of Antimicro- bial Agents47, 48 (2016)
work page 2016
-
[16]
R. Roemhild and H. Schulenburg, Evolutionary ecol- ogy meets the antibiotic crisis: Can we control pathogen adaptation through sequential therapy?, Evo- lution, Medicine, and Public Health2019, 37 (2019)
work page 2019
-
[17]
C. Barbosa, V. Trebosc, C. Kemmer, P. Rosenstiel, R. Beardmore, H. Schulenburg, and G. Jansen, Alterna- tive evolutionary paths to bacterial antibiotic resistance cause distinct collateral effects, Molecular Biology and Evolution34, 2229 (2017)
work page 2017
-
[18]
L. Imamovic, M. M. H. Ellabaan, A. M. D. Machado, L. Citterio, T. Wulff, S. Molin, H. K. Johansen, and M. O. A. Sommer, Drug-driven phenotypic convergence supports rational treatment strategies of chronic infec- tions, Cell172, 121 (2018)
work page 2018
-
[19]
N. L. Podnecky, E. G. Fredheim, J. Kloos, V. Sørum, R. Primicerio, A. P. Roberts, D. E. Rozen, Ø. Samuelsen, and P. J. Johnsen, Conserved collateral antibiotic suscep- tibility networks in diverse clinical strains of escherichia coli, Nature Communications9, 3673 (2018)
work page 2018
-
[20]
S. Hernando-Amado, F. Sanz-Garc´ ıa, and J. L. Mart´ ınez, Rapid and robust evolution of collateral sensitivity in pseudomonas aeruginosa antibiotic-resistant mutants, Science Advances6, eaba5493 (2020)
work page 2020
-
[21]
A. Fuentes-Hernandez, J. Plucain, F. Gori, R. Pena- Miller, C. Reding, G. Jansen, H. Schulenburg, I. Gudelj, and R. Beardmore, Using a sequential regimen to elimi- nate bacteria at sublethal antibiotic dosages, PLoS Biol- ogy13, e1002104 (2015)
work page 2015
- [22]
-
[23]
R. E. Beardmore, R. Pe˜ na-Miller, F. Gori, and J. Iredell, Antibiotic cycling and antibiotic mixing: which one best mitigates antibiotic resistance?, Molecular Biology and Evolution34, 802 (2017)
work page 2017
-
[24]
J. Maltas and K. B. Wood, Pervasive and diverse collat- eral sensitivity profiles inform optimal strategies to limit antibiotic resistance, PLoS Biology17, e3000515 (2019)
work page 2019
-
[25]
B. Morsky and D. C. Vural, Suppressing evolution of an- tibiotic resistance through environmental switching, The- oretical Ecology15, 115 (2022)
work page 2022
-
[26]
G. Katriel, Optimizing antimicrobial treatment sched- ules: some fundamental analytical results, Bulletin of Mathematical Biology86, 1 (2024)
work page 2024
-
[27]
D. T. Weaver, E. S. King, J. Maltas, and J. G. Scott, Reinforcement learning informs optimal treatment strategies to limit antibiotic resistance, Proceedings of the National Academy of Sciences of the U.S.A.121, e2303165121 (2024)
work page 2024
- [28]
-
[29]
L. B. Aulin, A. Liakopoulos, P. H. van der Graaf, D. E. Rozen, and J. C. van Hasselt, Design principles of collat- eral sensitivity-based dosing strategies, Nature Commu- nications12, 5691 (2021)
work page 2021
-
[30]
D. C. Angst, B. Tepekule, L. Sun, B. Bogos, and S. Bon- hoeffer, Comparing treatment strategies to reduce an- tibiotic resistance in an in vitro epidemiological setting, 15 Proceedings of the National Academy of Sciences of the U.S.A.118, e2023467118 (2021)
work page 2021
-
[31]
T. Gil-Gil, B. A. Berryhill, J. A. Manuel, A. P. Smith, I. C. McCall, F. Baquero, and B. R. Levin, The evo- lution of heteroresistance via small colony variants in Escherichia colifollowing long-term exposure to bacte- riostatic antibiotics, Nature Communications15, 7936 (2024)
work page 2024
-
[32]
D. T. Gillespie, Stochastic simulation of chemical kinet- ics, Annual Review of Physical Chemistry58, 35 (2007)
work page 2007
-
[33]
D. T. Gillespie, A general method for numerically sim- ulating the stochastic time evolution of coupled chemi- cal reactions, Journal of Computational Physics22, 403 (1976)
work page 1976
-
[34]
P. Ankomah and B. R. Levin, Exploring the collabora- tion between antibiotics and the immune response in the treatment of acute, self-limiting infections, Proceedings of the National Academy of Sciences of the U.S.A.111, 8331 (2014)
work page 2014
-
[35]
D. I. Andersson and D. Hughes, Evolution of antibiotic resistance at non-lethal drug concentrations, Drug Resis- tance Updates15, 162 (2012)
work page 2012
- [36]
-
[37]
F. Sanz-Garc´ ıa, T. Gil-Gil, P. Laborda, P. Blanco, L.- E. Ochoa-S´ anchez, F. Baquero, J. L. Mart´ ınez, and S. Hernando-Amado, Translating eco-evolutionary biol- ogy into therapy to tackle antibiotic resistance, Nature Reviews Microbiology21, 671 (2023)
work page 2023
-
[38]
D. Nichol, P. Jeavons, A. G. Fletcher, R. A. Bonomo, P. K. Maini, J. L. Paul, R. A. Gatenby, A. R. Ander- son, and J. G. Scott, Steering evolution with sequential therapy to prevent the emergence of bacterial antibiotic resistance, PLoS Computational Biology11, e1004493 (2015)
work page 2015
-
[39]
R. Roemhild, C. S. Gokhale, P. Dirksen, C. Blake, P. Rosenstiel, A. Traulsen, D. I. Andersson, and H. Schu- lenburg, Cellular hysteresis as a principle to maximize the efficacy of antibiotic therapy, Proceedings of the National Academy of Sciences of the U.S.A.115, 9767 (2018)
work page 2018
-
[40]
R. E. Beardmore and R. Pe˜ na-Miller, Rotating antibi- otics selects optimally against antibiotic resistance, in theory, Mathematical Biosciences and Engineering7, 527 (2010)
work page 2010
-
[41]
R. Beardmore and R. Pena-Miller, Antibiotic cycling ver- sus mixing: the difficulty of using mathematical models to definitively quantify their relative merits., Mathemat- ical Biosciences and Engineering7, 923 (2010)
work page 2010
-
[42]
C. Nyhoegen and H. Uecker, Sequential antibiotic ther- apy in the laboratory and in the patient, Journal of the Royal Society Interface20, 20220793 (2023)
work page 2023
-
[43]
F. Girardini, L. Antozzi, and R. Raiteri, Impact of antibi- otic changes in empirical therapy on antimicrobial resis- tance in intensive care unit-acquired infections, Journal of Hospital Infection52, 136 (2002)
work page 2002
-
[44]
M. M¨ uller, A. dela Pena, and H. Derendorf, Issues in pharmacokinetics and pharmacodynamics of anti- infective agents: distribution in tissue, Antimicrobial Agents and Chemotherapy48, 1441 (2004)
work page 2004
-
[45]
P. Czuppon, T. Day, F. D´ ebarre, and F. Blanquart, A stochastic analysis of the interplay between antibiotic dose, mode of action, and bacterial competition in the evolution of antibiotic resistance, PLoS Computational Biology19, e1011364 (2023)
work page 2023
-
[46]
M. Yoshida, S. G. Reyes, S. Tsuda, T. Horinouchi, C. Fu- rusawa, and L. Cronin, Time-programmable drug dosing allows the manipulation, suppression and reversal of an- tibiotic drug resistance in vitro, Nature communications 8, 15589 (2017)
work page 2017
-
[47]
C. Barbosa, R. Roemhild, P. Rosenstiel, and H. Schu- lenburg, Evolutionary stability of collateral sensitivity to antibiotics in the model pathogen pseudomonas aerugi- nosa, Elife8, e51481 (2019)
work page 2019
- [48]
-
[49]
S. Hernando-Amado, C. L´ opez-Causap´ e, P. Laborda, F. Sanz-Garc´ ıa, A. Oliver, and J. L. Mart´ ınez, Rapid phe- notypic convergence towards collateral sensitivity in clin- ical isolates of pseudomonas aeruginosa presenting dif- ferent genomic backgrounds, Microbiology Spectrum11, e02276 (2023)
work page 2023
-
[50]
S. Hernando-Amado, P. Laborda, R. La Rosa, S. Molin, H. K. Johansen, and J. L. Mart´ ınez, Ciprofloxacin re- sistance rapidly declines in nfxb defective clinical strains of pseudomonas aeruginosa, Nature Communications16, 4992 (2025)
work page 2025
-
[51]
J. Molina-Hern´ andez, Optimization of sequential thera- pies to maximize extinction of resistant bacteria through collateral sensitivity, GitHub 10.5281/zenodo.17485976 (2025)
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.