Dynamics of an HIV Model with Multiple Infection Stages and Treatment with Different Drug Classes

Highly active antiretroviral therapy can effectively control HIV replication in infected individuals. Some clinical and modeling studies suggested that viral decay dynamics may depend on the inhibited stages of the viral replication cycle. In this paper, we develop a general mathematical model incorporating multiple infection stages and various drug classes that can interfere with specific stages of the viral life cycle. We derive the basic reproductive number and obtain the global stability results of steady states. Using several simple cases of the general model, we study the effect of various drug classes on the dynamics of HIV decay. When drugs are assumed to be 100% effective, drugs acting later in the viral life cycle lead to a faster or more rapid decay in viremia. This is consistent with some patient and experimental data, and also agrees with previous modeling results. When drugs are not 100% effective, the viral decay dynamics are more complicated. Without a second population of long-lived infected cells, the viral load decline can have two phases if drugs act at an intermediate stage of the viral replication cycle. The slopes of viral load decline depend on the drug effectiveness, the death rate of infected cells at different stages, and the transition rate of infected cells from one to the next stage. With a second population of long-lived infected cells, the viral load decline can have three distinct phases, consistent with the observation in patients receiving antiretroviral therapy containing the integrase inhibitor raltegravir. We also fit modeling prediction to patient data under efavirenz (a nonnucleoside reverse-transcriptase inhibitor) and raltegravir treatment. The first-phase viral load decline under raltegravir therapy is longer than that under efavirenz, resulting in a lower viral load at initiation of the second-phase decline in patients taking raltegravir. This explains why patients taking a raltegravir-based therapy were faster to achieve viral suppression than those taking an efavirenz-based therapy. Taken together, this work provides a quantitative and systematic comparison of the effect of different drug classes on HIV decay dynamics and can explain the viral load decline in HIV patients treated with raltegravir-containing regimens.