Complex patterns of viral load decay under antiretroviral therapy: influence of pharmacokinetics and intracellular delay

We present a model of HIV dynamics under antiretroviral therapy that combines drug pharmacokinetics and intracellular delay. A two compartment pharmacokinetic model is employed to determine the time evolution of the intracellular concentrations of the active forms of drugs, and thereby drug efficacy. The viral replication period is divided into pre- and post-drug action parts, allowing for the introduction of an intracellular delay in drug action. The standard model of viral dynamics is modified to account for the drug dependence of intracellular delay and continuously varying drug efficacy. Model calculations reveal that viral load decay in HIV infected patients under monotherapy can exhibit remarkably complex patterns depending on the relative magnitudes of the pharmacokinetic, intracellular, and intrinsic viral dynamic time-scales. The commonly assumed exponential decay is only a special case. However, uncertainties in measurement and the low sampling frequencies employed in present clinical studies preclude the identification of these patterns from existing clinical viral load data.