Kinetic coevolutionary models predict the temporal emergence of HIV-1 resistance mutations under drug selection pressure

Drug resistance in HIV type 1 (HIV - 1) is a pervasive problem that affects the lives of millions of people worldwide. Although records of drug- resistant mutations (DRMs) have been extensively tabulated within public repositories, our understanding of the evolutionary kinetics of DRMs and how they evolve together remains limited. Epistasis, the interaction between a DRM and other residues in HIV - 1 protein sequences, is key to the temporal evolution of drug resistance. We use a Potts sequence- covariation statistical- energy model of HIV - 1 protein fitness under drug selection pressure, which captures epistatic interactions between all positions, combined with kinetic Monte - Carlo simulations of sequence evolutionary trajectories, to explore the acquisition of DRMs as they arise in an ensemble of drug-naive patient protein sequences. We follow the time course of 52 DRMs in the enzymes protease, RT, and integrase, the primary targets of antiretroviral therapy. The rates at which DRMs emerge are highly correlated with their observed acquisition rates reported in the literature when drug pressure is applied. This result highlights the central role of epistasis in determining the kinetics governing DRM emergence. Whereas rapidly acquired DRMs begin to accumulate as soon as drug pressure is applied, slowly acquired DRMs are contingent on accessory mutations that appear only after prolonged drug pressure. We provide a foundation for using computational methods to determine the temporal evolution of drug resistance using Potts statistical potentials, which can be used to gain mechanistic insights into drug resistance pathways in HIV - 1 and other infectious agents. Significance HIV - 1 affects the lives of millions of people worldwide; cases of panresistant HIV are emerging. We use a kinetic Monte - Carlo method to simulate the evolution of drug resistance based on HIV - 1 patient- derived sequence data available on public databases. Our simulations capture the reported time to acquire drug- resistance mutations (DRMs) across the major HIV - 1 drug- target enzymes: Protease, RT, and Integrase. The network of epistatic interactions with a primary DRM determines its acquisition rate, which is not explained by the overall fitness of the DRM or features of the genetic code, but instead by an "epistatic barrier". This work provides a framework for the development of computational methods that forecast the time- course over which resistance to antivirals develops in patients.