Optimal Control of Combined Therapy in a Single Strain HIV-1 Model

Highly active antiretroviral therapy (HAART) is administered to symptomatic human immunodeficiency virus (HIV) infected individuals to improve their health. Various administration schemes are used to improve patients?lives and at the same time suppressing development of drug resistance, reduce evolution of new viral strains, minimize serious side effects, improve patient adherence and also reduce the costs of drugs. We deduce an optimal drug administration scheme useful in improving patients? health especially in poor resourced settings. In this paper we use the Pontryagin?s Maximum Principle to derive optimal drug dosages based on a mathematical dynamical model. We use methods of optimal control to determine optimal controls analytically, and then use the Runge-Kutta scheme of order four to numerically simulate different therapy effects. We simulate the different effects of a drug regimen composed of a protease inhibitor and a nucleoside reverse transcriptase inhibitor. Our results indicate that for highly toxic drugs, small dosage sizes and allowing drug holidays make a profound impact in both improving the quality of life and reducing economic costs of therapy. The results show that for drugs with less toxicity, continuous therapy is beneficial.

[1]  A. Oxenius,et al.  Cytotoxic T Lymphocyte Responses to Human Immunodeficiency Virus: Control and Escape , 2000, Stem cells.

[2]  A. Perelson,et al.  Dynamics of HIV infection of CD4+ T cells. , 1993, Mathematical biosciences.

[3]  A. Perelson Modelling viral and immune system dynamics , 2002, Nature Reviews Immunology.

[4]  S. Lenhart,et al.  OPTIMIZING CHEMOTHERAPY IN AN HIV MODEL , 1998 .

[5]  D. Kirschner,et al.  Optimal control of the chemotherapy of HIV , 1997, Journal of mathematical biology.

[6]  Denise E. Kirschner,et al.  Using Mathematics to Understand HIV Immune Dynamics , 1997 .

[7]  Pini Gurfil,et al.  Optimal control of HIV infection with a continuously-mutating viral population , 2002, Proceedings of the 2002 American Control Conference (IEEE Cat. No.CH37301).

[8]  D. Kirschner,et al.  Predicting differential responses to structured treatment interruptions during HAART , 2004, Bulletin of mathematical biology.

[9]  P. Klenerman,et al.  HIV: current opinion in escapology. , 2002, Current opinion in microbiology.

[10]  Mark Kot,et al.  Elements of Mathematical Ecology: Frontmatter , 2001 .

[11]  A. Perelson,et al.  Complex patterns of viral load decay under antiretroviral therapy: influence of pharmacokinetics and intracellular delay. , 2004, Journal of theoretical biology.

[12]  Hem Raj Joshi,et al.  Optimal control of an HIV immunology model , 2002 .

[13]  Chung Choo Chung,et al.  Optimal Scheduling of Drug Treatment for HIV Infection : Continuous Dose Control and Receding Horizon Control , 2003 .

[14]  M. Nowak,et al.  Immune Responses and Viral Phenotype: Do Replication Rate and Cytopathogenicity Influence Virus Load? , 2000 .

[15]  V. Robert-Hebmann,et al.  Apoptosis of uninfected cells induced by HIV envelope glycoproteins , 2004, Retrovirology.

[16]  Glenn F. Webb,et al.  Immunotherapy of HIV-1 Infection , 1998 .

[17]  D. Krakauer,et al.  T–cell induced pathogenesis in HIV: bystander effects and latent infection , 1999, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[18]  V. Jansen,et al.  Effector cytotoxic T lymphocyte numbers induced by vaccination should exceed levels in chronic infection for protection from HIV. , 2001, Vaccine.