Bioinformatics Advance Access originally published online on October 31, 2008
Bioinformatics 2009 25(8):1040-1047; doi:10.1093/bioinformatics/btn568
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Stochastic modelling of genotypic drug-resistance for human immunodeficiency virus towards long-term combination therapy optimization
1Department of Computer Science and Automation, University of Roma TRE, 2Informa Contract Research Organisation, 3Infectious Disease Clinic, Catholic University of Sacred Heart, Rome, Italy and 4Department of Molecular Biology, University of Siena, Siena, Italy
*To whom correspondence should be addressed.
| Abstract |
|---|
Motivation: Several mathematical models have been investigated for the description of viral dynamics in the human body: HIV-1 infection is a particular and interesting scenario, because the virus attacks cells of the immune system that have a role in the antibody production and its high mutation rate permits to escape both the immune response and, in some cases, the drug pressure. The viral genetic evolution is intrinsically a stochastic process, eventually driven by the drug pressure, dependent on the drug combinations and concentration: in this article the viral genotypic drug resistance onset is the main focus addressed. The theoretical basis is the modelling of HIV-1 population dynamics as a predator–prey system of differential equations with a time-dependent therapy efficacy term, while the viral genome mutation evolution follows a Poisson distribution. The instant probabilities of drug resistance are estimated by means of functions trained from in vitro phenotypes, with a roulette-wheel-based mechanisms of resistant selection. Simulations have been designed for treatments made of one and two drugs as well as for combination antiretroviral therapies. The effect of limited adherence to therapy was also analyzed. Sequential treatment change episodes were also exploited with the aim to evaluate optimal synoptic treatment scenarios.
Results: The stochastic predator–prey modelling usefully predicted long-term virologic outcomes of evolved HIV-1 strains for selected antiretroviral therapy combinations. For a set of widely used combination therapies, results were consistent with findings reported in literature and with estimates coming from analysis on a large retrospective data base (EuResist).
Contact: prosperi{at}dia.uniroma3.it
Supplementary information: Supplementary data are available at Bioinformatics online.
Associate Editor: Martin Bishop
Received on July 22, 2008; revised on October 29, 2008; accepted on October 29, 2008