AN AGE-STRUCTURED MODEL WITH IMMUNE RESPONSE OF HIV INFECTION

This paper develops and analyzes an age-structured model of HIV infection with various compartments, including target cells, infected cells, viral loads and immune effector cells, to provide a better understanding of the interaction between HIV and the immune system. We show that the proposed model has one uninfected steady state and several infected steady states. We conduct a local stability analysis of these steady states by using a generalized Jacobian matrix method in conjunction with the Laplace transform. In addition, we consider various techniques and ideas from optimal control theory to derive optimal therapy protocols by using two types of dynamic treatment methods representing reverse transcriptase inhibitors and protease inhibitors. We derive the necessary conditions (an optimality system) for optimal control functions by considering the first variations of the Lagrangian. Further, we obtain optimal therapy protocols by solving a large optimality system of equations through the use of a difference scheme based on the Runge-Kutta method. The results of numerical simulations indicate that the optimal therapy protocols can facilitate long-term control of HIV through a strong immune response after the discontinuation of the therapy.

[1]  A two-strain hiv-1 mathematical model to assess the effects of chemotherapy on disease parameters. , 2005, Mathematical biosciences and engineering : MBE.

[2]  M. L. Chambers The Mathematical Theory of Optimal Processes , 1965 .

[3]  L. D. Pillis,et al.  A Validated Mathematical Model of Cell-Mediated Immune Response to Tumor Growth , 2005 .

[4]  M A Nowak,et al.  Anti-viral drug treatment: dynamics of resistance in free virus and infected cell populations. , 1997, Journal of theoretical biology.

[5]  S. Nagata,et al.  Essential roles of the Fas-Fas ligand pathway in the development of pulmonary fibrosis. , 1999, The Journal of clinical investigation.

[6]  B. Adams,et al.  Dynamic multidrug therapies for hiv: optimal and sti control approaches. , 2004, Mathematical biosciences and engineering : MBE.

[7]  Alan S Perelson,et al.  HIV-1 infection and low steady state viral loads , 2002, Bulletin of mathematical biology.

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

[9]  J. C. López-Marcos,et al.  Runge-Katta methods for age-structured popualtion models , 1995 .

[10]  Alan S. Perelson,et al.  Mathematical Analysis of HIV-1 Dynamics in Vivo , 1999, SIAM Rev..

[11]  Angela R. McLean,et al.  Zidovudine and HIV: Mathematical models of within‐host population dynamics , 1995 .

[12]  Guanrong Chen,et al.  Feedback control of a biodynamical model of HIV-1 , 2001, IEEE Transactions on Biomedical Engineering.

[13]  J. Corbeil,et al.  Nevirapine resistance mutations of human immunodeficiency virus type 1 selected during therapy , 1994, Journal of virology.

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

[15]  Philip Laird,et al.  Dynamic Optimization , 1992, ML.

[16]  Denise Kirschner,et al.  A model for treatment strategy in the chemotherapy of AIDS , 1996 .

[17]  Helen Moore,et al.  A mathematical model for treatment-resistant mutations of hiv. , 2005, Mathematical biosciences and engineering : MBE.

[18]  M A Nowak,et al.  Viral dynamics in vivo: limitations on estimates of intracellular delay and virus decay. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Veerle Fack,et al.  A parallel iterated correction algorithm for the numerical solution of Volterra equations , 1995 .

[20]  Max Gunzburger,et al.  Perspectives in flow control and optimization , 1987 .

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

[22]  José Álvarez-Ramírez,et al.  Feedback Control of the chemotherapy of HIV , 2000, Int. J. Bifurc. Chaos.

[23]  Alan S. Perelson,et al.  Mathematical Analysis of Age-Structured HIV-1 Dynamics with Combination Antiretroviral Therapy , 2007, SIAM J. Appl. Math..

[24]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[25]  G. F. Webb,et al.  Understanding drug resistance for monotherapy treatment of HIV infection , 1997, Bulletin of mathematical biology.

[26]  A S Perelson,et al.  Modeling plasma virus concentration during primary HIV infection. , 2000, Journal of theoretical biology.

[27]  Jeehyun Lee,et al.  Optimal control of an age-structured model of HIV infection , 2012, Appl. Math. Comput..

[28]  Daniel Coombs,et al.  An age-structured model of hiv infection that allows for variations in the production rate of viral particles and the death rate of productively infected cells. , 2004, Mathematical biosciences and engineering : MBE.

[29]  Carlos Castillo-Chavez,et al.  How May Infection-Age-Dependent Infectivity Affect the Dynamics of HIV/AIDS? , 1993, SIAM J. Appl. Math..

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

[31]  S. Ruan,et al.  A delay-differential equation model of HIV infection of CD4(+) T-cells. , 2000, Mathematical biosciences.

[32]  G. Webb Theory of Nonlinear Age-Dependent Population Dynamics , 1985 .

[33]  A. Trkola,et al.  HIV-1-specific immune responses in subjects who temporarily contain virus replication after discontinuation of highly active antiretroviral therapy. , 1999, The Journal of clinical investigation.

[34]  B. Adams,et al.  HIV dynamics: Modeling, data analysis, and optimal treatment protocols , 2005 .

[35]  J. Lisziewicz,et al.  Control of HIV despite the discontinuation of antiretroviral therapy. , 1999, The New England journal of medicine.

[36]  Carlos Castillo-Chavez,et al.  Diseases with chronic stage in a population with varying size. , 2003, Mathematical biosciences.

[37]  W. Fleming,et al.  Deterministic and Stochastic Optimal Control , 1975 .

[38]  Harvey Thomas Banks,et al.  A state‐dependent Riccati equation‐based estimator approach for HIV feedback control , 2006 .

[39]  L. S. Pontryagin,et al.  Mathematical Theory of Optimal Processes , 1962 .

[40]  Z. Feng,et al.  A Two-Strain Tuberculosis Model with Age of Infection , 2002, SIAM J. Appl. Math..

[41]  Jeehyun Lee,et al.  Free Terminal Time Optimal Control Problem of an HIV Model Based on a Conjugate Gradient Method , 2011, Bulletin of mathematical biology.

[42]  Harvey Thomas Banks,et al.  Feedback control of HIV antiviral therapy with long measurement time , 2011 .