An integer programming approach for optimal drug dose computation

In this paper, we study the problem of determining the optimal drug administration strategy when only a finite number of different dosages are available, a lower bound is posed on the time intervals between two consecutive doses, and drug concentrations should not exceed the toxic concentration levels. The presence of only binary variables leads to the adoption of an integer programming (IP) scheme for the formulation and solution of the drug dose optimal control problem. The proposed method is extended to account for the stochastic formulation of the optimal control problem, so that it can be used in practical applications where large populations of patients are to be treated. A Finite Impulse Response (FIR) model derived from experimental pharmacokinetic data is employed to correlate the administered drug dose with the concentration-time profiles of the drug in the compartments (organs) of the body.

[1]  Malcolm Rowland,et al.  Hepatic clearance of drugs. I. Theoretical considerations of a “well-stirred” model and a “parallel tube” model. Influence of hepatic blood flow, plasma and blood cell binding, and the hepatocellular enzymatic activity on hepatic drug clearance , 1977, Journal of Pharmacokinetics and Biopharmaceutics.

[2]  Michael Nikolaou,et al.  FIR model identification: Parsimony through kernel compression with wavelets , 1998 .

[3]  Thomas J. McAvoy,et al.  Nonlinear FIR Modeling via a Neural Net PLS Approach , 1996 .

[4]  A Schumitzky,et al.  Design of dosage regimens: a multiple model stochastic control approach. , 1994, International journal of bio-medical computing.

[5]  W. Lovejoy,et al.  An efficient control strategy for dosage regimens , 1994, Journal of Pharmacokinetics and Biopharmaceutics.

[6]  E. Lee,et al.  Weighted target interval stochastic control methods with global optimization and their applications in individualizing therapy , 2007, Journal of Pharmacokinetics and Pharmacodynamics.

[7]  Panos M. Pardalos,et al.  Modeling and integer programming techniques applied to propositional calculus , 1990, Comput. Oper. Res..

[8]  Jan M. Maciejowski,et al.  Predictive control : with constraints , 2002 .

[9]  Younggil Kwon,et al.  Handbook of Essential Pharmacokinetics, Pharmacodynamics and Drug Metabolism for Industrial Scientists , 2001, Springer US.

[10]  Pantelis Sopasakis,et al.  A Model Predictive Control Approach for Optimal Drug Administration , 2009 .

[11]  Dimitri P. Bertsekas,et al.  Dynamic Programming and Optimal Control, Vol. II , 1976 .

[12]  R. Ocampo-Pérez,et al.  Adsorption of Fluoride from Water Solution on Bone Char , 2007 .

[13]  Alberto Bemporad,et al.  Control of systems integrating logic, dynamics, and constraints , 1999, Autom..

[14]  Selection of optimal dose. , 1985, Computers and biomedical research, an international journal.

[15]  J. Macgregor,et al.  Identification of finite impulse response models : Methods and robustness issues , 1996 .

[16]  Efstratios N. Pistikopoulos,et al.  Global Optimization Issues in Multiparametric Continuous and Mixed-Integer Optimization Problems , 2004, J. Glob. Optim..

[17]  C Kulcsár,et al.  Optimal experimental design and therapeutic drug monitoring. , 1994, International journal of bio-medical computing.

[18]  D Bayard,et al.  Achieving target goals most precisely using nonparametric compartmental models and "multiple model" design of dosage regimens. , 2000, Therapeutic drug monitoring.

[19]  J. Au,et al.  Physiologically Based Pharmacokinetic Models of 2′,3′-Dideoxyinosine , 1997, Pharmaceutical Research.

[20]  S S Hacisalihzade,et al.  Solution of the multiple dosing problem using linear programming. , 1985, International journal of bio-medical computing.

[21]  Paolo Zunino,et al.  Multidimensional Pharmacokinetic Models Applied to the Design of Drug-Eluting Stents , 2004 .

[22]  Dimitri P. Bertsekas,et al.  Dynamic Programming and Optimal Control, Two Volume Set , 1995 .

[23]  S S Hacisalihzade,et al.  A simple algorithm for the solution of the multiple dosing problem. , 1985, Computer methods and programs in biomedicine.

[24]  K. Black,et al.  Intravenous levodopa administration in humans based on a two-compartment kinetic model , 2007, Journal of Neuroscience Methods.

[25]  J. Hardman,et al.  Factors affecting drug absorption and distribution , 2008 .

[26]  Jacek M. Zurada,et al.  Individualization of pharmacological anemia management using reinforcement learning , 2005, Neural Networks.

[27]  Lennart Ljung,et al.  Linear approximations of nonlinear FIR systems for separable input processes , 2005, Autom..

[28]  K. Itoh,et al.  Factors Affecting the Pharmacokinetics of CPT-11: The Body Mass Index, Age and Sex are Independent Predictors of Pharmacokinetic Parameters of CPT-11 , 2004, Investigational New Drugs.

[29]  Richard Magin,et al.  Fractional kinetics in multi-compartmental systems , 2010, Journal of Pharmacokinetics and Pharmacodynamics.

[30]  J. Doull,et al.  Dose and time as variables of toxicity. , 2000, Toxicology.

[31]  M. Nomoto,et al.  Pharmacokinetics of l-dopa in plasma and extracellular fluid of striatum in common marmosets , 2003, Brain Research.

[32]  L Aarons,et al.  Population pharmacokinetics: theory and practice. , 1991, British journal of clinical pharmacology.

[33]  J. Rodman,et al.  Targeting the systemic exposure of teniposide in the population and the individual using a stochastic therapeutic objective , 1993, Journal of Pharmacokinetics and Biopharmaceutics.

[34]  P. Hansen,et al.  Essays and surveys in global optimization , 2005 .

[35]  V. A. Roshchin,et al.  Models and methods of solution of quadratic integer programming problems , 1987 .