Pharmacokinetics in organs and the intact body: model validation and reduction

Putting square pegs in round holes: why traditional pharmacokinetic principles cannot universally be applied to iron-carbohydrate complexes

Intravenous iron-carbohydrate complexes are nanomedicines that are commonly used to treat iron deficiency and iron deficiency anemia of various etiologies. Many challenges remain regarding these complex drugs in the context of fully understanding their pharmacokinetic parameters. Firstly, the measurement of the intact iron nanoparticles versus endogenous iron concentration fundamentally limits the availability of data for computational modeling. Secondly, the models need to include several parameters to describe the iron metabolism which is not completely defined and those identified (e.g. ferritin) exhibit considerable interpatient variability. Additionally, modeling is further complicated by the lack of traditional receptor/enzyme interactions. The known parameters of bioavailability, distribution, metabolism, and excretion for iron-carbohydrate nanomedicines will be reviewed and future challenges that currently prevent the direct application of physiologically-based pharmacokinetic or other computational modeling techniques will be discussed.

Towards wearable and implantable continuous drug monitoring: A review

Continuous drug monitoring is a promising alternative to current therapeutic drug monitoring strategies and has a strong potential to reshape our understanding of pharmacokinetic variability and to improve individualised therapy. This review highlights recent advances in biosensing technologies that support continuous drug monitoring in real time. We focus primarily on aptamer-based biosensors, wearable and implantable devices. Emphasis is given to the approaches employed in constructing biosensors. We pay attention to sensors' biocompatibility, calibration performance, long-term characteristics stability and measurement quality. Last, we discuss the current challenges and issues to be addressed in continuous drug monitoring to make it a promising, future tool for individualised therapy. The ongoing efforts are expected to result in fully integrated implantable drug biosensing technology. Thus, we may anticipate an era of advanced healthcare in which wearable and implantable biochips will automatically adjust drug dosing in response to patient health conditions, thus enabling the management of diseases and enhancing individualised therapy.

Physiologically based structure of mean residence time

A mean residence time (MRT) is an important pharmacokinetic parameter. To the author's knowledge, however, a physiologically based structure of MRT (thereafter MRT structure) has not been published so far. Primarily this is because MRT structures cannot be identified by traditional pharmacokinetic methods used for the determination of MRT. Therefore, tools from the theory of linear dynamic systems were used for the structural identification of MRT in this study. The MRT structure identified is physiologically meaningful. Accordingly, it seems that the MRT structure identified may contribute to already established knowledge about MRT.

Drug structure-transport relationships

Malcolm Rowland has greatly facilitated an understanding of drug structure–pharmacokinetic relationships using a physiological perspective. His view points, covering a wide range of activities, have impacted on my own work and on my appreciation and understanding of our science. This overview summarises some of our parallel activities, beginning with Malcolm’s work on the pH control of amphetamine excretion, his work on the disposition of aspirin and on the application of clearance concepts in describing the disposition of lidocaine. Malcolm also spent a considerable amount of time developing principles that define solute structure and transport/pharmacokinetic relationships using in situ organ studies, which he then extended to involve the whole body. Together, we developed a physiological approach to studying hepatic clearance, introducing the convection–dispersion model in which there was a spread in blood transit times through the liver accompanied by permeation into hepatocytes and removal by metabolism or excretion into the bile. With a range of colleagues, we then further developed the model and applied it to various organs in the body. One of Malcolm’s special interests was in being able to apply this knowledge, together with an understanding of physiological differences in scaling up pharmacokinetics from animals to man. The description of his many other activities, such as the development of clearance concepts, application of pharmacokinetics to the clinical situation and using pharmacokinetics to develop new compounds and delivery systems, has been left to others.

Tissue-level modeling of xenobiotic metabolism in liver: An emerging tool for enabling clinical translational research

This review summarizes some of the recent developments and identifies critical challenges associated with in vitro and in silico representations of the liver and assesses the translational potential of these models in the quest of rationalizing the process of evaluating drug efficacy and toxicity. It discusses a wide range of research efforts that have produced, during recent years, quantitative descriptions and conceptual as well as computational models of hepatic processes such as biotransport and biotransformation, intra- and intercellular signal transduction, detoxification, etc. The above mentioned research efforts cover multiple scales of biological organization, from molecule-molecule interactions to reaction network and cellular and histological dynamics, and have resulted in a rapidly evolving knowledge base for a "systems biology of the liver." Virtual organ/organism formulations represent integrative implementations of particular elements of this knowledge base, usually oriented toward the study of specific biological endpoints, and provide frameworks for translating the systems biology concepts into computational tools for quantitative prediction of responses to stressors and hypothesis generation for experimental design.

Bioavailability, distribution and depletion of monensin in chickens

The pharmacokinetics of monensin including apparent volume of distribution, total body clearance, systemic bioavailability, partition coefficients and tissue residues were determined in chickens. The drug was given by intravenous injection in the left wing vein at the dose of 0.46 mg/kg and by intracrop administration at the dose of 4 mg/kg according to a destructive sampling. The pharmacokinetic variables were compared after noncompartmental, naïve averaged, naïve pooled and nonlinear mixed-effects modelling analyses. Partition coefficients and tissue residues were determined after a treatment with feed additives (125 mg/kg of feed) of 33 days. The clearance, volume of distribution and bioavailabilty were approximately 2.2 L/h/kg, approximately 9 L/kg and approximately 30% respectively except with nonlinear mixed effects models that presented values of 1.77 L/h/kg, 14.05 L/kg and 11.36% respectively. Tissue/plasma partition coefficients were estimated to 0.83, 3.39 and 0.51 for liver, fat and thigh muscle respectively. Monensin residues after treatment were not detected 6 h after withdrawal except for fat where monensin was still quantifiable 12 h after. Pharmacokinetic variables seem to be inaccurate when assessed with non linear mixed-effects modelling associated to destructive sampling in chickens. Values varied slightly with noncompartmental, naïve averaged and naïve pooled analyses. The absorption, elimination and partition parameters will be incorporated into a physiologically based pharmacokinetic model and the depletion study will be used to test the ability of this model to describe monensin residues in edible tissues under different dosage regimens.

Advanced pharmacokinetic models based on organ clearance, circulatory, and fractal concepts

Three advanced models of pharmacokinetics are described. In the first class are physiologically based pharmacokinetic models based on in vitro data on transport and metabolism. The information is translated as transporter and enzyme activities and their attendant heterogeneities into liver and intestine models. Second are circulatory models based on transit time distribution and plasma concentration time curves. The third are fractal models for nonhomogeneous systems and non-Fickian processes are presented. The usefulness of these pharmacokinetic models, with examples, is compared.

A minimal physiological model of thiopental distribution kinetics based on a multiple indicator approach

Currently available models of thiopental disposition kinetics using only plasma concentration-time data neglect the influence of intratissue diffusion and provide no direct information on tissue partitioning in individual subjects. Our approach was based on a lumped-organ recirculatory model that has recently been applied to unbound compounds. The goal was to find the simplest model that accounts for the heterogeneity in tissue partition coefficients and accurately describes initial distribution kinetics of thiopental in dogs. To ensure identifiability of the underlying axially distributed capillary-tissue exchange model, simultaneously measured disposition data of the vascular indicator, indocyanine green, and the marker of whole body water, antipyrine, were analyzed together with those of thiopental. A model obtained by grouping the systemic organs in two subsystems containing fat and nonfat tissues, successfully described all data and allowed an accurate estimation of model parameters. The estimated tissue partition coefficients were in accordance with those measured in rats. Because of the effect of tissue binding, the diffusional equilibration time characterizing intratissue distribution of thiopental is longer than that of antipyrine. The approach could potentially be used in clinical pharmacokinetics and could increase our understanding of the effect of obesity on the disposition kinetics of lipid-soluble drugs.

Mechanistic modeling of digoxin distribution kinetics incorporating slow tissue binding

This study aims to develop a mechanistic pharmacokinetic model that accounts for the kinetics of tissue binding in order to evaluate the effect of slow binding of digoxin to skeletal muscular Na(+)/K(+)-ATPase in humans. The approach is based on a minimal circulatory model with a systemic transit time density function that accounts for vascular mixing, transcapillary permeation and extravascular binding of the drug. The model parameters were estimated using previously published disposition data of digoxin in healthy volunteers and physiological distribution volumes taken from the literature. A time constant of the binding process of 34min was estimated indicating that receptor binding and not permeation clearance is the rate-limiting step of the distribution process. Model simulations suggest that up- or downregulation of sodium pumps, typically observed under physiological or pathophysiological conditions, could be detected with this method. The model allows a quantitative prediction of the effect of changes in skeletal muscular sodium pump activity on plasma levels of digoxin.

Whole body pharmacokinetic models

The aim of the current review is to summarise the present status of physiologically based pharmacokinetic (PBPK) modelling and its applications in drug research, and thus serve as a reference point to people interested in the methodology. The review is structured into three major sections. The first discusses the existing methodologies and techniques of PBPK model development. The second describes some of the most interesting PBPK model implementations published. The final section is devoted to a discussion of the current limitations and the possible future developments of the PBPK modelling approach. The current review is focused on papers dealing with the pharmacokinetics and/or toxicokinetics of medicinal compounds; references discussing PBPK models of environmental compounds are mentioned only if they represent considerable methodological developments or reveal interesting interpretations and/or applications.The major conclusion of the review is that, despite its significant potential, PBPK modelling has not seen the development and implementation it deserves, especially in the drug discovery, research and development processes. The main reason for this is that the successful development and implementation of a PBPK model is seen to require the investment of significant experience, effort, time and resources. Yet, a substantial body of PBPK-related research has been accumulated that can facilitate the PBPK modelling and implementation process. What is probably lagging behind is the expertise component, where the demand for appropriately qualified staff far outreaches availability.

Pharmacokinetics of idarubicin in the isolated perfused rat lung: effect of cinchonine and rutin

This study was designed to examine the effect of rutin and cinchonine on the uptake and metabolism of idarubicin (IDA) in the isolated perfused rat lung. IDA (2 mg) was infused for 2 min into the truncus pulmonalis in the presence of P-glycoprotein (P-gp) modulators cinchonine (1 microM) or rutin (6 microM). (Rutin is also known as an aldo-keto reductase inhibitor.) Venous outflow samples were collected up to 60 min, and the concentration of IDA and its primary metabolite idarubicinol (IDOL) were measured by high-performance liquid chromatography) with fluorescence detection. Thereafter, the tissue concentrations of IDA and IDOL were determined in the lung (n = 5 in each group). The estimated mean transit times for IDA in the treatment groups (MTT(cinchonine) = 21.8+/-3.5 min; MTT(rutin) = 20.1+/-5.0 min) were significantly higher than in the control group (11.6+/-2.1 min). Both cinchonine and rutin significantly enhanced the lung tissue concentrations of IDA (1.7- and 2.4-fold), as well as of IDOL (2.1- and 2.4-fold). Cinchonine and rutin also increased the outflow recovery of IDOL 2.6- and 2.7-fold, respectively. The results suggest that uptake kinetics of IDA into the rat lung is partly controlled by a P-gp efflux pump and its inhibition enhances the accumulation of IDA.

Prediction of the volume of distribution of a drug: which tissue-plasma partition coefficients are needed?

The aim of this study was to identify the tissue-plasma partition coefficients (Kp) needed for an initial prediction of the volume of distribution at steady state (Vd(ss)) of a drug in humans. Values of Kp were collected from the literature. Only Kp values plausibly representing true steady state distribution were accepted, and data had to be available for muscle, fat, skin and at least five other organs. The apparent volume of distribution of a drug in an organ/tissue (Vapp) was calculated as Kp multiplied by the volume of the organ/tissue, and the Vd(ss) as the sum of all available Vapp values. The percentage contribution of each Vapp to the Vd(ss) was estimated. In addition, linear regressions were calculated between Kp values of all drugs in a specific organ/tissue and Kp in muscle or fat. Finally, the Vd(ss) was re-calculated using (for basic drugs) the Kp in fat to calculate Vapp in fat and lungs and the Kp in muscle for the Vapp of all other organs/tissues. The two sets of estimates of Vd(ss) were compared by linear regression. The same calculations were performed for acidic drugs, except that muscle Kp was used also forthe lungs. Distribution to fat and muscle accounted for 84% (61-91%) (median and range) of the total estimated Vd(ss) of the basic drugs (n = 17). The regressions between Kp in organs/tissues and muscle Kp were statistically significant except in the case of liver. For acidic drugs (n = 18), distribution to fat and muscle accounted for 65% (42-92%) of Vd(ss), and the regressions of Kp were significant for all organs/tissues except kidney and bone. For both types of drugs, correlations between organ/tissue Kp values and Kp in fat were generally worse. There were excellent linear correlations between Vd(ss) calculated by means of only two Kp values and the originally calculated Vd(ss) (r2> or = 0.99 for both basic and acidic drugs; slopes were notsignificantly different from unity). Thus, initial estimation of the Vd(ss) of a new drug can normally be based on only two Kp values, those of muscle and fat. The muscle Kp can be used to represent all lean tissues, including the residual "carcass", with the exception that fat Kp can be used for distribution of basic drugs to lungs.

P-glycoprotein inhibitors enhance saturable uptake of idarubicin in rat heart: pharmacokinetic/pharmacodynamic modeling

Little is known about cardiac uptake kinetics of idarubicin, including a possible protective role of P-glycoprotein (Pgp)-mediated transport. This study therefore investigated uptake and negative inotropic action of idarubicin in the single-pass isolated perfused rat heart by using a pharmacokinetic/pharmacodynamic modeling approach. Idarubicin was administered as a 10-min constant infusion of 0.5 mg followed by a 70-min washout period in the absence and presence of the Pgp antagonists verapamil or amiodarone. Outflow concentration and left ventricular developed pressure were measured and the model parameters were estimated by simultaneous nonlinear regression. The results indicate the existence of a saturable, Michaelis-Menten type uptake process into the heart (K(m) = 3.06 microM, V(max) = 46.0 microM/min). Verapamil and amiodarone significantly enhanced the influx rate (V(max) increased 1.8-fold), suggesting that idarubicin is transported by Pgp directly out of the membrane before it gets into the cell. Verapamil and amiodarone attenuated the negative inotropic action of idarubicin, which was linked to the intracellular concentration of idarubicin.

Kinetics of propranolol uptake in muscle, skin, and fat of the isolated perfused rat hindlimb

The tissue distribution kinetics of a highly bound solute, propranolol, was investigated in a heterogeneous organ, the isolated perfused limb, using the impulse-response technique and destructive sampling. The propranolol concentration in muscle, skin, and fat as well as in outflow perfusate was measured up to 30 min after injection. The resulting data were analysed assuming (1) vascular, muscle, skin and fat compartments as well mixed (compartmental model) and (2) using a distributed-in-space model which accounts for the noninstantaneous intravascular mixing and tissue distribution processes but consists only of a vascular and extravascular phase (two-phase model). The compartmental model adequately described propranolol concentration-time data in the three tissue compartments and the outflow concentration-time curve (except of the early mixing phase). In contrast, the two-phase model better described the outflow concentration-time curve but is limited in accounting only for the distribution kinetics in the dominant tissue, the muscle. The two-phase model well described the time course of propranolol concentration in muscle tissue, with parameter estimates similar to those obtained with the compartmental model. The results suggest, first that the uptake kinetics of propranolol into skin and fat cannot be analysed on the basis of outflow data alone and, second that the assumption of well-mixed compartments is a valid approximation from a practical point of view (as, e.g., in physiological based pharmacokinetic modelling). The steady-state distribution volumes of skin and fat were only 16 and 4%, respectively, of that of muscle tissue (16.7 ml), with higher partition coefficient in fat (6.36) than in skin (2.64) and muscle (2.79).

19F-MRS studies of fluorinated drugs in humans

The use of 19F-NMR as a noninvasive probe to measure directly the pharmacokinetics of drugs at their target (effector) site(s) is illustrated in this article by human studies with 5-fluorouracil (5-FU). This drug, and several of its metabolites, have been measured in vivo in animals and in patients using standard clinical MRI systems. Using a pharmacokinetic imaging approach the parameter that can be measured most readily is the tumoral t(1/2) of 5-FU. Patients whose tumoral t(1/2) of 5-FU is equal to/greater than 20 min are designated as "trappers", and those whose tumoral t(1/2) of 5-FU is less are nontrappers. Trapping of 5-FU in tumors is a necessary, albeit not a sufficient condition, for response. Problems associated with the technical aspects of these measurements have been discussed, as well as how modulators and other agents will affect the tumoral t(1/2) of 5-FU. The rationale for the biological processes underlying the fate of 5-FU in humans has been illustrated with the use of a 12 compartment model, where several of the steps have been discussed and the consequences of their inhibition/stimulation related to the noninvasive studies that can be performed with modulators of the action of 5-FU. These 19F-NMR studies have now been extended to other fluoropyrimidines, some of which are prodrugs of 5-FU, and others where the fluorine atoms are on the ribose ring. These studies also reveal information that has both scientific and clinical significance. The studies presented here illustrate some of the potential and some of the usefulness of 19F-MRS in patient management and in drug development. It is a technique that has proven itself.

The anomalous pharmacokinetics of amiodarone explained by nonexponential tissue trapping

Conventional pharmacokinetic (PK) concepts fail to describe the long-term pharmacokinetics of the extremely cationic amphiphilic drug amiodarone. A nonclassical model based on the phenomenon of trapping at tissue binding sites with very long release times is presented, which implies that a volume of distribution and a steady-state level cannot be defined. In agreement with clinical PK data available in the literature, the model well describes not only single-dose disposition curves but also the persistently increasing plasma concentration-time curve during long-term treatment (up to 5 years) and the washout curve following cessation of therapy. The novel aspect is a long-tailed tissue residence time distribution which is incorporated into a recirculatory model leaving the initial distribution process and the clearance concept unchanged. The underlying theoretical approach, which is known as "strange or anomalous" kinetics in physical sciences, and the fractal scaling property of the model may enhance our understanding of the PK of extremely hydrophobic xenobiotics.

Cellular pharmacokinetics: effects of cytoplasmic diffusion and binding on organ transit time distribution

Distribution between well-stirred compartments is the classical paradigm in pharmacokinetics. Also in capillary-issue exchange modeling a barrier-limited approach is mostly adopted. As a consequence of tissue binding, however, drug distribution cannot be regarded as instantaneous even at the cellular level and the distribution process consists of at least two components: transmembrane exchange and cytoplasmic transport. Two concepts have been proposed for the cytoplasmic distribution process of hydrophobic or amphipathic molecules, (i) slowing of diffusion due to instantaneous binding to immobile cellular structures and (ii) slow binding after instantaneous distribution throughout the cytosol. The purpose of this study was to develop a general approach for comparing both models using a stochastic model of intra- and extravascular drug distribution. Criteria for model discrimination are developed using the first three central moments (mean, variance, and skewness) of the cellular residence time and organ transit time distribution, respectively. After matching the models for the relative dispersion the remaining differences in relative skewness are predicted, discussing the relative roles of membrane permeability, cellular binding and cytoplasmic transport. It is shown under which conditions the models are indistinguishable on the basis of venous organ outflow concentration-time curves. The relative dispersion of cellular residence times is introduced as a model-independent measure of cytoplasmic equilibration kinetics, which indicates whether diffusion through the cytoplasm is rate limiting. If differences in outflow curve shapes (their relative skewness) cannot be detected, independent information on binding and/or diffusion kinetics is necessary to avoid model misspecification. The method is applied to previously published hepatic outflow data of enalaprilat, triiodothyronine, and diclofenac. It provides a general framework for the modeling of cellular pharmacokinetics.