On the validity of the dispersion model of hepatic drug elimination when intravascular transit time densities are long-tailed

The dispersion model with mixed boundary conditions uses a single parameter, the dispersion number, to describe the hepatic elimination of xenobiotics and endogenous substances. An implicit a priori assumption of the model is that the transit time density of intravascular indicators is approximated by an inverse Gaussian distribution. This approximation is limited in that the model poorly describes the tail part of the hepatic outflow curves of vascular indicators. A sum of two inverse Gaussian functions is proposed as ail alternative, more flexible empirical model for transit time densities of vascular references. This model suggests that a more accurate description of the tail portion of vascular reference curves yields an elimination rate constant (or intrinsic clearance) which is 40% less than predicted by the dispersion model with mixed boundary conditions. The results emphasize the need to accurately describe outflow curves in using them as a basis for determining pharmacokinetic parameters using hepatic elimination models. (C) 1997 Society for Mathematical Biology.

Neonatal hepatic drug elimination

After the transition from in utero to newborn life, the neonate becomes solely reliant upon its own drug clearance processes to metabolise xenobiotics. Whilst most studies of neonatal hepatic drug elimination have focussed upon in vitro expression and activities of drug-metabolising enzymes, the rapid physiological changes in the early neonatal period of life also need to be considered. There are dramatic changes in neonatal liver blood how and hepatic oxygenation due to the loss of the umbilical blood supply, the increasing portal vein blood flow, and the gradual closure of the ductus venosus shunt during the first week of life. These changes which may well affect the capacity of neonatal hepatic drug metabolism. The hepatic expression of cytochromes P450 1A2, 2C, 2D6, 2E1 and 3A4 develop at different rates in the postnatal period, whilst 3A7 expression diminishes. Hepatic glucuronidation in the human neonate is relatively immature at birth, which contrasts with the considerably more mature neonatal hepatic sulfation activity. Limited in vivo studies show that the human neonate can significantly metabolise xenobiotics but clearance is considerably less compared with the older infant and adult. The neonatal population included in pharmacological studies is highly heterogeneous with respect to age, body weight, ductus venosus closure and disease processes, making it difficult to interpret data arising from human neonatal studies. Studies in the perfused foetal and neonatal sheep liver have demonstrated how the oxidative and conjugative hepatic elimination of drugs by the intact organ is significantly increased during the first week of life, highlighting that future studies will need to consider the profound physiological changes that may influence neonatal hepatic drug elimination shortly after birth.

Analytical validation for a series of marker compounds used to assess renal drug elimination processes

Renal drug elimination is determined by glomerular filtration, tubular secretion, and tubular reabsorption. Changes in the integrity of these processes influence renal drug clearance, and these changes may not be detected by conventional measures of renal function such as creatinine clearance. The aim of the current study was to examine the analytic issues needed to develop a cocktail of marker drugs (fluconazole, rac-pindolol, para-aminohippuric acid, sinistrin) to measure simultaneously the mechanisms contributing to renal clearance. High-performance liquid chromatographic methods of analysis for fluconazole, pindolol, para-aminohippuric acid, and creatinine and an enzymatic assay for sinistrin were developed or modified and then validated to allow determination of each of the compounds in both plasma and urine in the presence of all other marker drugs. A pilot clinical study in one volunteer was conducted to ensure that the assays were suitable for quantitating all the marker drugs to the sensitivity and specificity needed to allow accurate determination of individual renal clearances. The performance of all assays (plasma and urine) complied with published validation criteria. All standard curves displayed linearity over the concentration ranges required, with coefficients of correlation greater than 0.99. The precision of the interday and intraday variabilities of quality controls for each marker in plasma and urine were all less than 11.9% for each marker. Recoveries of markers (and internal standards) in plasma and urine were all at least 90%. All markers investigated were shown to be stable when plasma or urine was frozen and thawed. For all the assays developed, there were no interferences from other markers or endogenous substances. In a pilot clinical study, concentrations of all markers could be accurately and reproducibly determined for a sufficient duration of time after administration to calculate accurate renal clearance for each marker. This article presents details of the analytic techniques developed for measuring concentrations of marker drugs for different renal elimination processes administered as a single dose to define the processes contributing to renal drug elimination.

Simultaneous administration of a cocktail of markers to measure renal drug elimination pathways: absence of a pharmacokinetic interaction between fluconazole and sinistrin, p-aminohippuric acid and pindolol

Aims Previous studies suggest that estimated creatinine clearance, the conventional measure of renal function, does not adequately reflect charges in renal drug handling in some patients, including the immunosuppressed. The aim of this study was to develop and validate a cocktail of markers. to be given in a single administration, capable of detecting alterations in the renal elimination pathways of glomerular filtration, tubular secretion and tubular reabsorption. Methods Healthy male subjects (n = 12) received intravenously infused 2500 mg sinistrin (glomerular filtration) and 440 mg p-aminohippuric acid (PAH; anion secretion), and orally administered 100 mg fluconazole (reabsorption) and 15 mg rac-pindolol (cation secretion). The potential interaction between these markers was investigated in a pharmacokinetic study where markers (M) or fluconazole (F) were administered alone or together (M + F). Validated analytical methods were used to measure plasma and urine concentrations in order to quantify the renal handling of each marker. Plasma protein binding of fluconazole was measured by ultrafiltration. All subjects had an estimated creatinine clearance within the normal range. The renal clearance of each marker (Mean +/- s.d.) was calculated as the ratio of the amount excreted in urine and thearea-under-the-concentration-time curve. Statistical comparisons were made using a paired t-test and 95% confidence intervals were reported. Results The renal clearances of sinistrin (M: 119 +/- 31 ml min(-1); M + F: 130 +/- 40 ml min(-1); P = 0.32), PAH (M: 469 +/- 145 ml min(-1); M + F: 467 +/- 146 ml min(-1); P = 0.95), R-pindolol (M: 204 +/- 41 ml min(-1); M + F: 190 +/- 41 ml min(-1); P = 0.39; n = 11), S-pindolol (M: 225 +/- 55 ml min(-1); M + F: 209 +/- 60 ml min(-1); P = 0.27; n = 11) and fluconazole (F: 14.9 +/-3.8 ml min(-1); M + F: 13.6 +/- 3.4 ml min(-1); P = 0.16) were similar when the markers or fluconazole were administered alone (M or F) or as a cocktail (M + F). Conclusions This study found no interaction between markers and fluconazole in healthy male subjects, suggesting that a single administration of this cocktail of markers of different renal processes call be used to simultaneously investigate pathways of renal drug elimination.

Impulse-response studies on tracer doses of [C-14]lignocaine and its multiple metabolites in the perfused rat liver

The outflow-concentration-time profiles for lignocaine (lidocaine) and its metabolites have been measured after bolus impulse administration of [C-14]lignocaine into the perfused rat liver. Livers from female Sprague-Dawley rats were perfused in a once-through fashion with red-blood-cell-free Krebs-Henseleit buffer containing 0 or 2% bovine serum albumin. Perfusate flow rates of 20 and 30 mL min(-1) were used and both normal and retrograde flow directions were employed. Significant amounts of metabolite were detected in the effluent perfusate soon after lignocaine injection. The early appearance of metabolite contributed to bimodal outflow profiles observed for total C-14 radioactivity. The lignocaine outflow profiles were well characterized by the two-compartment dispersion model, with efflux rate << influx rate. The profiles for lignocaine metabolites were also characterized in terms of a simplified two-compartment dispersion model. Lignocaine was found to be extensively metabolized under the experimental conditions with the hepatic availability ranging between 0.09 and 0.18. Generally lignocaine and metabolite availability showed no significant change with alterations in perfusate flow rate from 20 to 30 mt min(-1) or protein content from 0 to 2%. A significant increase in lignocaine availability occurred when 1200 mu M unlabelled lignocaine was added to the perfusate. Solute mean transit times generally decreased with increasing flow rate and with increasing perfusate protein content. The results confirm that lignocaine pharmacokinetics in the liver closely follow the predictions of the well-stirred model. The increase in lignocaine availability when 1200 mu M unlabelled lignocaine was added to the perfusate is consistent with saturation of the hydroxylation metabolic pathways of lignocaine metabolism.

Metabolite mean transit times in the liver as predicted by various models of hepatic elimination

Predicted area under curve (AUC), mean transit time (MTT) and normalized variance (CV2) data have been compared for parent compound and generated metabolite following an impulse input into the liver, Models studied were the well-stirred (tank) model, tube model, a distributed tube model, dispersion model (Danckwerts and mixed boundary conditions) and tanks-in-series model. It is well known that discrimination between models for a parent solute is greatest when the parent solute is highly extracted by the liver. With the metabolite, greatest model differences for MTT and CV2 occur when parent solute is poorly extracted. In all cases the predictions of the distributed tube, dispersion, and tasks-in-series models are between the predictions of the rank and tube models. The dispersion model with mixed boundary conditions yields identical predictions to those for the distributed tube model (assuming an inverse gaussian distribution of tube transit times). The dispersion model with Danckwerts boundary conditions and the tanks-in series models give similar predictions to the dispersion (mixed boundary conditions) and the distributed tube. The normalized variance for parent compound is dependent upon hepatocyte permeability only within a distinct range of permeability values. This range is similar for each model but the order of magnitude predicted for normalized variance is model dependent. Only for a one-compartment system is the MIT for generated metabolite equal to the sum of MTTs for the parent compound and preformed metabolite administered as parent.

Modeling of hepatic elimination and organ distribution kinetics with the extended convection-dispersion model

The conventional convection-dispersion (also called axial dispersion) model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. An extended form of the convection-dispersion model has been developed to adequately describe the outflow concentration-time profiles for vascular markers at both short and long times after bolus injections into perfused livers. The model, based on flux concentration and a convolution of catheters and large vessels, assumes that solute elimination in hepatocytes follows either fast distribution into or radial diffusion in hepatocytes. The model includes a secondary vascular compartment, postulated to be interconnecting sinusoids. Analysis of the mean hepatic transit time (MTT) and normalized variance (CV2) of solutes with extraction showed that the discrepancy between the predictions of MTT and CV2 for the extended and conventional models are essentially identical irrespective of the magnitude of rate constants representing permeability, volume, and clearance parameters, providing that there is significant hepatic extraction. In conclusion, the application of a newly developed extended convection-dispersion model has shown that the unweighted conventional convection-dispersion model can be used to describe the disposition of extracted solutes and, in particular, to estimate hepatic availability and clearance in booth experimental and clinical situations.

Ion-trapping, microsomal binding, and unbound drug distribution in the hepatic retention of basic drugs

This study investigated the relative contribution of ion-trapping, microsomal binding, and distribution of unbound drug as determinants in the hepatic retention of basic drugs in the isolated perfused rat liver. The ionophore monensin was used to abolish the vesicular proton gradient and thus allow an estimation of ion-trapping by acidic hepatic vesicles of cationic drugs. In vitro microsomal studies were used to independently estimate microsomal binding and metabolism. Hepatic vesicular ion-trapping, intrinsic elimination clearance, permeability-surface area product, and intracellular binding were derived using a physiologically based pharmacokinetic model. Modeling showed that the ion-trapping was significantly lower after monensin treatment for atenolol and propranolol, but not for antipyrine. However, no changes induced by monensin treatment were observed in intrinsic clearance, permeability, or binding for the three model drugs. Monensin did not affect binding or metabolic activity in vitro for the drugs. The observed ion-trapping was similar to theoretical values estimated using the pHs and fractional volumes of the acidic vesicles and the pK(a) values of drugs. Lipophilicity and pK(a) determined hepatic drug retention: a drug with low pK(a) and low lipophilicity (e.g., antipyrine) distributes as unbound drug, a drug with high pK(a) and low lipophilicity (e.g., atenolol) by ion-trapping, and a drug with a high pK(a) and high lipophilicity (e.g., propranolol) is retained by ion-trapping and intracellular binding. In conclusion, monensin inhibits the ion-trapping of high pK(a) basic drugs, leading to a reduction in hepatic retention but with no effect on hepatic drug extraction.

A model to study intestinal and hepatic metabolism of propranolol in the dog

A model to investigate hepatic drug uptake and metabolism in the dog was developed for this study. Catheters were placed in the portal and hepatic veins during exploratory laparotomy to collect pre- and posthepatic blood samples at defined intervals. Drug concentrations in the portal vein were taken to reflect intestinal uptake and metabolism of an p.o. administered drug (propranolol), while differences in drug and metabolite concentrations between portal and hepatic veins reflected hepatic uptake and metabolism. A significant difference in propranolol concentration between hepatic and portal veins confirmed a high hepatic extraction of this therapeutic agent in the dog. This technically uncomplicated model may be used experimentally or clinically to determine hepatic function and metabolism of drugs that may be administered during anaesthesia and surgery.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

Principles and clinical application of assessing alterations in renal elimination pathways

Drugs and metabolites are eliminated from the body by metabolism and excretion. The kidney makes the major contribution to excretion of unchanged drug and also to excretion of metabolites. Net renal excretion is a combination of three processes - glomerular filtration, tubular secretion and tubular reabsorption. Renal function has traditionally been determined by measuring plasma creatinine and estimating creatinine clearance. However, estimated creatinine clearance measures only glomerular filtration with a small contribution from active secretion. There is accumulating evidence of poor correlation between estimated creatinine clearance and renal drug clearance in different clinical settings, challenging the 'intact nephron hypothesis' and suggesting that renal drug handling pathways may not decline in parallel. Furthermore, it is evident that renal drug handling is altered to a clinically significant extent in a number of disease states, necessitating dosage adjustment not just based on filtration. These observations suggest that a re-evaluation of markers of renal function is required. Methods that measure all renal handling pathways would allow informed dosage individualisation using an understanding of renal excretion pathways and patient characteristics. Methodologies have been described to determine individually each of the renal elimination pathways. However, their simultaneous assessment has only recently been investigated. A cocktail of markers to measure simultaneously the individual renal handling pathways have now been developed, and evaluated in healthy volunteers. This review outlines the different renal elimination pathways and the possible markers that can be used for their measurement. Diseases and other physiological conditions causing altered renal drug elimination are presented, and the potential application of a cocktail of markers for the simultaneous measurement of drug handling is evaluated. Further investigation of the effects of disease processes on renal drug handling should include people with HIV infection, transplant recipients (renal and liver) and people with rheumatoid arthritis. Furthermore, changes in renal function in the elderly, the effect of sex on renal function, assessment of living kidney donors prior to transplantation and the investigation of renal drug interactions would also be potential applications. Once renal drug handling pathways are characterised in a patient population, the implications for accurate dosage individualisation can be assessed. The simultaneous measurement of renal function elimination pathways of drugs and metabolites has the potential to assist in understanding how renal function changes with different disease states or physiological conditions. In addition, it will further our understanding of fundamental aspects of the renal elimination of drugs.

Physiologically based synthetic models of hepatic disposition

Current Physiologically based pharmacokinetic (PBPK) models are inductive. We present an additional, different approach that is based on the synthetic rather than the inductive approach to modeling and simulation. It relies on object-oriented programming A model of the referent system in its experimental context is synthesized by assembling objects that represent components such as molecules, cells, aspects of tissue architecture, catheters, etc. The single pass perfused rat liver has been well described in evaluating hepatic drug pharmacokinetics (PK) and is the system on which we focus. In silico experiments begin with administration of objects representing actual compounds. Data are collected in a manner analogous to that in the referent PK experiments. The synthetic modeling method allows for recognition and representation of discrete event and discrete time processes, as well as heterogeneity in organization, function, and spatial effects. An application is developed for sucrose and antipyrine, administered separately and together PBPK modeling has made extensive progress in characterizing abstracted PK properties but this has also been its limitation. Now, other important questions and possible extensions emerge. How are these PK properties and the observed behaviors generated? The inherent heuristic limitations of traditional models have hindered getting meaningful, detailed answers to such questions. Synthetic models of the type described here are specifically intended to help answer such questions. Analogous to wet-lab experimental models, they retain their applicability even when broken apart into sub-components. Having and applying this new class of models along with traditional PK modeling methods is expected to increase the productivity of pharmaceutical research at all levels that make use of modeling and simulation.

Hepatic pharmacokinetics of propranolol in rats with adjuvant-induced systemic inflammation

Systemic inflammation is known to affect drug disposition in the liver. This study sought to relate and quantitate changes in hepatic pharmacokinetics of propranolol with changes in hepatic architecture and physiology in adjuvant-treated rats. Transmission electron microscopy was used to assess morphological changes in mitochondria and lysosomes of adjuvant-treated rat livers. The disposition of propranolol was assessed in the perfused rat liver using the multiple indicator dilution technique. Hepatic extraction and mean transit time were determined from outflow-concentration profiles using a nonparametric method. Kinetic parameters were derived from a two-phase physiologically based organ pharmacokinetic model. Possible relationships were then explored between the changes in hepatic drug disposition and cytochrome P-450 activity and iron concentration. Adjuvant treatment induced the appearance of mitochondrial inclusions/tubularization and irregularly shaped lysosomes in rat livers. Livers from adjuvant-treated rats had (relative to normal) significantly higher alpha(1)-acid glycoprotein (orosomucoid) and iron tissue concentrations but lower cytochrome P-450 content. The hepatic extraction, metabolism, and ion trapping of propranolol were significantly impaired in adjuvant-treated rats and could be correlated with altered iron store and cytochrome P-450 activity. It is concluded that adjuvant-induced systemic inflammation alters hepatocellular morphology and biochemistry and consequently influences hepatic disposition of propranolol.

Pharmacokinetic software for the health sciences - Choosing the right package for teaching purposes

Computer assisted learning has an important role in the teaching of pharmacokinetics to health sciences students because it transfers the emphasis from the purely mathematical domain to an 'experiential' domain in which graphical and symbolic representations of actions and their consequences form the major focus for learning. Basic pharmacokinetic concepts can be taught by experimenting with the interplay between dose and dosage interval with drug absorption (e.g. absorption rate, bioavailability), drug distribution (e.g. volume of distribution, protein binding) and drug elimination (e.g. clearance) on drug concentrations using library ('canned') pharmacokinetic models. Such 'what if' approaches are found in calculator-simulators such as PharmaCalc, Practical Pharmacokinetics and PK Solutions. Others such as SAAM II, ModelMaker, and Stella represent the 'systems dynamics' genre, which requires the user to conceptualise a problem and formulate the model on-screen using symbols, icons, and directional arrows. The choice of software should be determined by the aims of the subject/course, the experience and background of the students in pharmacokinetics, and institutional factors including price and networking capabilities of the package(s). Enhanced learning may result if the computer teaching of pharmacokinetics is supported by tutorials, especially where the techniques are applied to solving problems in which the link with healthcare practices is clearly established.

Evolution of resistance to sulfadoxine-pyrimethamine in Plasmodium falciparum

The development of resistance to sulfadoxine-pyrimethamine by Plasmodium parasites is a major problem for the effective treatment of malaria, especially P. falciparum malaria. Although the molecular basis for parasite resistance is known, the factors promoting the development and transmission of these resistant parasites are less clear. This paper reports the results of a quantitative comparison of factors previously hypothesized as important for the development of drug resistance, drug dosage, time of treatment, and drug elimination half-life, with an in-host dynamics model of P. falciparum malaria in a malaria-naive host. The results indicate that the development of drug resistance can be categorized into three stages. The first is the selection of existing parasites with genetic mutations in the dihydrofolate reductase or dihydropteroate synthetase gene. This selection is driven by the long half-life of the sulfadoxine-pyrimethamine combination. The second stage involves the selection of parasites with allelic types of higher resistance within the host during an infection. The timing of treatment relative to initiation of a specific anti-P. falciparum EMP1 immune response is an important factor during this stage, as is the treatment dosage. During the third stage, clinical treatment failure becomes prevalent as the parasites develop sufficient resistance mutations to survive therapeutic doses of the drug combination. Therefore, the model output reaffirms the importance of correct treatment of confirmed malaria cases in slowing the development of parasite resistance to sulfadoxine-pyrimethamine.