Bioinformatics Advance Access first published online on October 31, 2008
This version published online on November 3, 2008
Bioinformatics, doi:10.1093/bioinformatics/btn568
Stochastic Modelling of Genotypic Drug-Resistance for Human Immunodeficiency Virus towards Long-Term Combination Therapy Optimisation
1; Department of Computer Science and Automation, University of Roma TRE, Rome, Italy; 2 Informa Contract Research Organisation, Rome, Italy; 3 Infectious Disease Clinic, Catholic University of Sacred Heart, Rome, Italy; 4 Department of Molecular Biology, University of Siena, Italy
*To whom correspondence should be addressed. Dr. Mattia C.F. Prosperi, E-mail: prosperi{at}dia.uniroma3.it
| 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 paper 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
Associate Editor: Prof. Martin Bishop
Received on July 22, 2008; revised on October 29, 2008; accepted on October 29, 2008