PKQuest: a general physiologically based pharmacokinetic model. Introduction and application to propranolol

BACKGROUND: A "physiologically based pharmacokinetic" (PBPK) approach uses a realistic model of the animal to describe the pharmacokinetics. Previous PBPKs have been designed for specific solutes, required specification of a large number of parameters and have not been designed for general use. METHODS: This new PBPK program (PKQuest) includes a "Standardhuman" and "Standardrat" data set so that the user input is minimized. It has a simple user interface, graphical output and many new features: 1) An option that uses the measured plasma concentrations to solve for the time course of the gastrointestinal, intramuscular, intraperotineal or skin absorption and systemic availability of a drug – for a general non-linear system. 2) Capillary permeability limitation defined in terms of the permeability-surface area products. 4) Saturable plasma and tissue protein binding. 5) A lung model that includes perfusion-ventilation mismatch. 6) A general optimization routine using either a global (simulated annealing) or local (Powell) minimization applicable to all model parameters. RESULTS: PKQuest was applied to measurements of human propranolol pharmacokinetics and intestinal absorption. A meal has two effects: 1) increases portal blood flow by 50%; and 2) decreases liver metabolism by 20%. There is a significant delay in the oval propranolol absorption in fasting subjects that is absent in fed subjects. The oral absorption of the long acting form of propranolol continues for a period of more than 24 hours. CONCLUSIONS: PKQuest provides a new general purpose, easy to use, freely distributed and physiologically rigorous PBPK software routine

Physiologically based pharmacokinetic modeling of arterial – antecubital vein concentration difference

BACKGROUND: Modeling of pharmacokinetic parameters and pharmacodynamic actions requires knowledge of the arterial blood concentration. In most cases, experimental measurements are only available for a peripheral vein (usually antecubital) whose concentration may differ significantly from both arterial and central vein concentration. METHODS: A physiologically based pharmacokinetic (PBPK) model for the tissues drained by the antecubital vein (referred to as "arm") is developed. It is assumed that the "arm" is composed of tissues with identical properties (partition coefficient, blood flow/gm) as the whole body tissues plus a new "tissue" representing skin arteriovenous shunts. The antecubital vein concentration depends on the following parameters: the fraction of "arm" blood flow contributed by muscle, skin, adipose, connective tissue and arteriovenous shunts, and the flow per gram of the arteriovenous shunt. The value of these parameters was investigated using simultaneous experimental measurements of arterial and antecubital concentrations for eight solutes: ethanol, thiopental, (99)Tc(m)-diethylene triamine pentaacetate (DTPA), ketamine, D(2)O, acetone, methylene chloride and toluene. A new procedure is described that can be used to determine the arterial concentration for an arbitrary solute by deconvolution of the antecubital concentration. These procedures are implemented in PKQuest, a general PBPK program that is freely distributed . RESULTS: One set of "standard arm" parameters provides an adequate description of the arterial/antecubital vein concentration for ethanol, DTPA, thiopental and ketamine. A significantly different set of "arm" parameters was required to describe the data for D(2)O, acetone, methylene chloride and toluene – probably because the "arm" is in a different physiological state. CONCLUSIONS: Using the set of "standard arm" parameters, the antecubital vein concentration can be used to determine the whole body PBPK model parameters for an arbitrary solute without any additional adjustable parameters. Also, the antecubital vein concentration can be used to estimate the arterial concentration for an arbitrary input for solutes for which no arterial concentration data is available

Tikhonov adaptively regularized gamma variate fitting to assess plasma clearance of inert renal markers

The Tk-GV model fits Gamma Variates (GV) to data by Tikhonov regularization (Tk) with shrinkage constant, λ, chosen to minimize the relative error in plasma clearance, CL (ml/min). Using 169Yb-DTPA and 99mTc-DTPA (n = 46, 8–9 samples, 5–240 min) bolus-dilution curves, results were obtained for fit methods: (1) Ordinary Least Squares (OLS) one and two exponential term (E1 and E2), (2) OLS-GV and (3) Tk-GV. Four tests examined the fit results for: (1) physicality of ranges of model parameters, (2) effects on parameter values when different data subsets are fit, (3) characterization of residuals, and (4) extrapolative error and agreement with published correction factors. Test 1 showed physical Tk-GV results, where OLS-GV fits sometimes-produced nonphysical CL. Test 2 showed the Tk-GV model produced good results with 4 or more samples drawn between 10 and 240 min. Test 3 showed that E1 and E2 failed goodness-of-fit testing whereas GV fits for t > 20 min were acceptably good. Test 4 showed CLTk-GV clearance values agreed with published CL corrections with the general result that CLE1 > CLE2 > CLTk-GV and finally that CLTk-GV were considerably more robust, precise and accurate than CLE2, and should replace the use of CLE2 for these renal markers

The pharmacokinetics of the interstitial space in humans

BACKGROUND: The pharmacokinetics of extracellular solutes is determined by the blood-tissue exchange kinetics and the volume of distribution in the interstitial space in the different organs. This information can be used to develop a general physiologically based pharmacokinetic (PBPK) model applicable to most extracellular solutes. METHODS: The human pharmacokinetic literature was surveyed to tabulate the steady state and equilibrium volume of distribution of the solutes mannitol, EDTA, morphine-6-glucuronide, morphine-3-glucuronide, inulin and β-lactam antibiotics with a range of protein binding (amoxicillin, piperacillin, cefatrizine, ceforanide, flucloxacillin, dicloxacillin). A PBPK data set was developed for extracellular solutes based on the literature for interstitial organ volumes. The program PKQuest was used to generate the PBPK model predictions. The pharmacokinetics of the protein (albumin) bound β-lactam antibiotics were characterized by two parameters: 1) the free fraction of the solute in plasma; 2) the interstitial albumin concentration. A new approach to estimating the capillary permeability is described, based on the pharmacokinetics of the highly protein bound antibiotics. RESULTS: About 42% of the total body water is extracellular. There is a large variation in the organ distribution of this water – varying from about 13% of total tissue water for skeletal muscle, up to 70% for skin and connective tissue. The weakly bound antibiotics have flow limited capillary-tissue exchange kinetics. The highly protein bound antibiotics have a significant capillary permeability limitation. The experimental pharmacokinetics of the 11 solutes is well described using the new PBPK data set and PKQuest. CONCLUSIONS: Only one adjustable parameter (systemic clearance) is required to completely characterize the PBPK for these extracellular solutes. Knowledge of just this systemic clearance allows one to predict the complete time course of the absolute drug concentrations in the major organs. PKQuest is freely available

The comparative biology of New Zealand oystercatchers

Oystercatchers comprise a distinctive group of mollusc-eating shorebirds. They form an extremely uniform monogeneric family which has not undergone any major adaptive radiations into a diversity of ecological niches, but rather has dispersed from original centres of distribution to occupy identical niches in new geographical localities. The uniformity of structure and habit displayed within the group has been attributed by Larson (1957) to a high ecobiotic specialisation with centripetal selection involved. Throughout their range, oystercatchers exploit identical ecological niches which require specialised habits for successful utilisation. The specialised feeding habits of oystercatchers are well documented (Murphy, 1925; Dewar, 1940; Larson, 1957; Tinbergen and Norton-Griffiths, 1964; Dare, 1966), and a natural consequence of this specialisation is that it is restrictive to adaptive radiation