Fitness landscapes emerging from pharmacodynamic functions in the evolution of multidrug resistance

Adaptive evolution often involves beneficial mutations at more than one locus. In this case, the trajectory and rate of adaptation is determined by the underlying fitness landscape, that is, the fitness values and mutational connectivity of all genotypes under consideration. Drug resistance, especially resistance to multiple drugs simultaneously, is also often conferred by mutations at several loci so that the concept of fitness landscapes becomes important. However, fitness landscapes underlying drug resistance are not static but dependent on drug concentrations, which means they are influenced by the pharmacodynamics of the drugs administered. Here, I present a mathematical framework for fitness landscapes of multidrug resistance based on Hill functions describing how drug concentrations affect fitness. I demonstrate that these ‘pharmacodynamic fitness landscapes’ are characterized by pervasive epistasis that arises through (i) fitness costs of resistance (even when these costs are additive), (ii) nonspecificity of resistance mutations to drugs, in particular cross‐resistance, and (iii) drug interactions (both synergistic and antagonistic). In the latter case, reciprocal drug suppression may even lead to reciprocal sign epistasis, so that the doubly resistant genotype occupies a local fitness peak that may be difficult to access by evolution. Simulations exploring the evolutionary dynamics on some pharmacodynamic fitness landscapes with both constant and changing drug concentrations confirm the crucial role of epistasis in determining the rate of multidrug resistance evolution.