Optimal switching strategies in multi-drug therapies for chronic diseases

Abstract

Antimicrobial resistance is a threat to public health with millions of deaths linked to drug resistant infections every year. To mitigate resistance, common strategies that are used are combination therapies and therapy switching. However, the stochastic nature of pathogenic mutation makes the optimization of these strategies challenging. Here, we propose a two-scale stochastic model that considers the effective evolution of therapies in a multidimensional efficacy space, where each dimension represents the efficacy of a specific drug in the therapy. The diffusion of therapies within this space is subject to stochastic resets, representing therapy switches. The boundaries of the space, inferred from coarser pathogen-host dynamics, can be either reflecting or absorbing. Reflecting boundaries impede full recovery of the host, while absorbing boundaries represent the development of antimicrobial resistance, leading to therapy failure. We derive analytical expressions for the average absorption times, accounting for both continuous and discrete genomic changes using the frameworks of Langevin and Master equations, respectively. These expressions allow us to evaluate the relevance of times between drug-switches and the number of simultaneous drugs in relation to typical timescales for drug resistance development. We also explore realistic scenarios where therapy constraints are imposed to the number of administered therapies and/or their costs, finding non-trivial optimal drug-switching protocols that maximize the time before antimicrobial resistance develops while reducing therapy costs.

Abstract: arXiv descriptive metadata (CC0). Thank you to arXiv for use of its open access interoperability. This service was not reviewed or approved by, nor does it necessarily express or reflect the policies or opinions of, arXiv.