Optimal control of drug therapy: Melding pharmacokinetics with viral dynamics

Pharmacokinetics were melded with a viral dynamical model to design an optimal drug administration regimen such that the basic reproductive number for the virus was minimized. One-compartmental models with two kinds of drug delivery routes, intravenous and extravascular with multiple dosages, and two drug elimination rates, first order and Michaelis-Menten rates, were considered. We defined explicitly the basic reproductive number for the viral dynamical model melded with pharmacokinetics. When the average plasma drug concentration was constant, intravenous administration of the drug with small dosages applied frequently minimized the basic reproductive number. For extravascular administration, the basic reproductive number initially decreases to a trough point and then increases as the drug dosage increases. When a therapeutic window is considered, numerical studies indicate that the wider the window, the smaller the basic reproductive number. Once the width of the therapeutic window is fixed, the basic reproductive number monotonously declines as the minimum therapeutic level increases. The findings suggest that the existence of drug dosage and drug administration interval that minimize the basic reproductive number could help design the optimal drug administration regimen.

[1]  Lynn McFadyen,et al.  Basic PK/PD principles of drug effects in circular/proliferative systems for disease modelling , 2010, Journal of Pharmacokinetics and Pharmacodynamics.

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

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

[4]  A. Ding,et al.  Assessing antiviral potency of anti-HIV therapies in vivo by comparing viral decay rates in viral dynamic models. , 2001, Biostatistics.

[5]  Gaston H. Gonnet,et al.  On the LambertW function , 1996, Adv. Comput. Math..

[6]  D. Bainov,et al.  Impulsive Differential Equations: Periodic Solutions and Applications , 1993 .

[7]  Robert J Smith Explicitly accounting for antiretroviral drug uptake in theoretical HIV models predicts long-term failure of protease-only therapy. , 2008, Journal of theoretical biology.

[8]  R. Anderson,et al.  The dynamics of drug action on the within-host population growth of infectious agents: melding pharmacokinetics with pathogen population dynamics. , 1998, Journal of theoretical biology.

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

[10]  D. Ho,et al.  Ordered accumulation of mutations in HIV protease confers resistance to ritonavir , 1996, Nature Medicine.

[11]  Jianhong Wu,et al.  Threshold virus dynamics with impulsive antiretroviral drug effects , 2011, Journal of Mathematical Biology.

[12]  J. Hale,et al.  Ordinary Differential Equations , 2019, Fundamentals of Numerical Mathematics for Physicists and Engineers.

[13]  Hulin Wu,et al.  Modeling HIV dynamics and antiviral response with consideration of time-varying drug exposures, adherence and phenotypic sensitivity. , 2003, Mathematical biosciences.

[14]  Yanni Xiao,et al.  Threshold dynamics for an HIV model in periodic environments , 2010 .

[15]  Nicolas Bacaër Approximation of the Basic Reproduction Number R0 for Vector-Borne Diseases with a Periodic Vector Population , 2007, Bulletin of mathematical biology.

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

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

[18]  A A Ding,et al.  Relationships between antiviral treatment effects and biphasic viral decay rates in modeling HIV dynamics. , 1999, Mathematical biosciences.

[19]  L M Wahl,et al.  Distinct effects of protease and reverse transcriptase inhibition in an immunological model of HIV-1 infection with impulsive drug effects , 2004, Bulletin of mathematical biology.

[20]  Michael A. Parniak,et al.  Effectiveness of 3TC in HIV clinical trials may be due in part to the M184V substitution in 3TC‐resistant HIV‐1 reverse transcriptase , 1996, AIDS.

[21]  Sanyi Tang,et al.  One-compartment model with Michaelis-Menten elimination kinetics and therapeutic window: an analytical approach , 2007, Journal of Pharmacokinetics and Pharmacodynamics.

[22]  Hulin Wu,et al.  Modeling Long-Term HIV Dynamics and Antiretroviral Response: Effects of Drug Potency, Pharmacokinetics, Adherence, and Drug Resistance , 2005, Journal of acquired immune deficiency syndromes.

[23]  J Elliott,et al.  Principles of pharmacodynamics and their applications in veterinary pharmacology. , 2004, Journal of veterinary pharmacology and therapeutics.

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

[25]  Yanni Xiao,et al.  Threshold dynamics for compartmental epidemic models with impulses , 2012 .

[26]  D. Kirschner,et al.  A methodology for performing global uncertainty and sensitivity analysis in systems biology. , 2008, Journal of theoretical biology.

[27]  S. Disney,et al.  On the Lambert W function: EOQ applications and pedagogical considerations , 2010 .

[28]  Yanni Xiao,et al.  The effects of population dispersal and pulse vaccination on disease control , 2010, Math. Comput. Model..

[29]  Xiao-Qiang Zhao,et al.  Threshold Dynamics for Compartmental Epidemic Models in Periodic Environments , 2008 .

[30]  Libin Rong,et al.  Mathematical modeling of viral kinetics under immune control during primary HIV-1 infection. , 2009, Journal of theoretical biology.