Development of a population simulation model for HIV monotherapy virological outcomes using lamivudine

Current highly active antiretroviral therapy (HAART) requires the use of combinations of three drugs to minimize the early emergence of drug-resistant HIV strains. Therefore, long-term monotherapy data with new agents are unavailable. However, the development of computer models for Monte-Carlo-type simulations of antiviral monotherapy, which incorporate HIV infection dynamic distributions from previously studied populations, together with pharmacokinetics and pharmacodynamic parameters of the new agent, could serve as an important tool. The nucleoside lamivudine (3TC) was used as a representative drug to standardize an improved pharmacodynamic and infection dynamic monotherapy model. 3TC plasma concentration versus time profiles was used to drive the cellular accumulation of 3TC-triphosphate (TP) in primary human lymphocytes in the model, over a 16 week period. The fraction of HIV reverse transcription inhibited was calculated using the median inhibitory concentration and intracellular 3TC-TP levels. Virus loads and activated CD4+ T-cell counts were generated for 2,200 theoretical individuals and compared with the outcomes of an actual 3TC monotherapy trial at the same dose. Pharmacokinetic variance alone did not account for the interindividual HIV-load variability. However, selection of appropriate distributions of the various pharmacokinetic and infection dynamics parameters produced a similar range of virus load reductions to actual observations. Therefore, once parameter and variance distributions are standardized, this modelling approach could be helpful in planning clinical trials and predicting the antiviral contribution of each agent in a HAART modality.