Tacrolimus dose individualization in "de novo" patients after 10 years of experience in liver transplantation: pharmacokinetic considerations and patient pathophysiology.

AIM: To determine how changes in tacrolimus (TAC) immunosuppression clinical practice, in the first 15 days post liver transplantation (LT) and across a decade, impact a clinical covariate - pharmacokinetic (PK) model, developed in data from 1998, thus testing its utility in dose individualization across time. Patient cohorts from 1998 (Reference: R-1998) and 2007 (EVALUATION: E-2007) were compared. METHODS: Analysis of monitoring observations (Cmin and Cmin/dose) and the biochemical variables aspartate aminotransferase (AST), hematocrit (HCT), albumin (ALB) and serum creatinine (SCr) was done for 0 - 3 and 4 - 15 days post transplantation (PT). The population PK model developed for R-1998 [1] was re-evaluated for the two cohorts. RESULTS: Significant differences in R-1998 vs. E-2007 existed in Cmin and Cmin/dose and in covariates AST (as hepatic function marker) and SCr (as toxicity marker). E-2007 had lower levels of Cmin and Cmin/dose (1/CL), lower AST with faster recovery and lower variability in Cmin/dose for similar dose. AST was a covariate on CL/F in the 0 - 3 day PT period. In 4 - 15 days PT for E-2007, low levels of HCT and ALB as CL/F predictors confirmed a subgroup with higher CL/F (23.8 l/h vs. 19.3 l/h). The R-1998 model's original structure was confirmed. CONCLUSIONS: Ten years of use of TAC shows gain in therapeutic targeting efficiency, due to improvement in LT methods, knowledge of the drug and consideration of PK steady state. The remaining uncertainty with TAC monitoring in LT can be resolved with application of PK principles combined with patients' diosyncrasies in the model developed for TAC dose individualization in R-1998. The applicability of the model as nucleus in Bayes individualization remains intact across a decade.

A predictive pharmacokinetic/pharmacodynamic model of fentanyl for analgesia/sedation in neonates based on a semi-physiologic approach.

BACKGROUND AND OBJECTIVES: Fentanyl is a synthetic opioid commonly used as an anesthetic and also increasingly popular as a sedative agent in neonates. Initial dosage regimens in this population are often empirically derived from adults on a body weight basis. However, ontogenic maturation processes related to drug disposition are not necessarily always body weight correlates. We developed a predictive pharmacokinetic/pharmacodynamic model that includes growth and maturation physiologic changes for fentanyl in neonatal care. METHODS: Key pharmacokinetic variables and principles (protein binding, clearance, distribution) as related to fentanyl pharmacokinetic/pharmacodynamic behavior in adults (tricompartmental model) and to neonatal physiologic data (organ weights and blood flows, body composition, renal and hepatic function, etc.) were used to guide the building of a semi-physiologic ontogenic model. The model applies to a normal-term neonate without any other intervention. Then, extension to a pharmacokinetic/pharmacodynamic link model for fentanyl was made. The final model was evaluated by predicting the time course of plasma concentrations and the effect of a standard regimen of 10.5 µg/kg over a 1-h period followed by 1.5 µg/kg/h for 48 h. RESULTS: Hepatic clearance was linked to ontogeny of unbound fraction and of α1-acid glycoprotein. All parameters were reduced in the neonate compared to adults but in a differing proportion due to qualitative changes in physiology that are analyzed and accounted for. Systemic clearance (CLS), volume of the central compartment (V1) and steady-state volume of distribution predicted by the model were 0.028 L/min, 1.26 L, and 22.04 L, respectively. Weight-corrected parameters generally decreased in adults compared with neonates, but differentially, e.g., CLS = 0.0093 versus 0.0088 L/min/kg, while V1 = 0.42 versus 0.18 L/kg (neonates vs. adults). Under such complexity a pharmacokinetic/pharmacodynamic model is the appropriate method for rational efficacy targeting. Fentanyl pharmacodynamics in neonates were considered to be similar to those in adults except for the equilibrium rate constant, which was also scaled on an ontogenic basis. The model adequately predicted the reported mean expected concentration-time profiles for the standard regimen. CONCLUSIONS: Integrated pharmacokinetic/pharmacodynamic modeling showed that the usually prescribed dosage regimens of fentanyl in neonates may not always provide the optimum degree of sedation. The model could be used in optimal design of clinical trials for this vulnerable population. Prospective clinical testing is the reasonable next step.

Tacrolimus pharmacokinetics in the early post-liver transplantation period and clinical applicability via Bayesian prediction.

PURPOSE: To define and validate a pharmacokinetic (PK) model for tacrolimus (TAC) that includes patient pathophysiology and has clinical applicability in the first 2 weeks post-liver transplantation (PLT). METHODS: Routine monitoring records [dose, trough levels (C(min)), demographics, biochemistry] from 75 patients treated with TAC (Prograf®) PLT were used to develop a population PK model (employing NONMEM®) testing for predictors of oral clearance (CL/F) according to bedside evidence and primarily with aspartate aminotransferase (AST), albumin (ALB), and hematocrit (HCT). Patients were catergorized into subgroups with above and below "normal" thresholds for AST (500 U/L), ALB (2.5 g/dL), and HCT (28 %), respectively. The model was validated with ten patients from the same period and 15 more recent patients. An empirical Bayes method was developed and applied to the prediction of individual profiles serving as a dose adjustment tool. RESULTS: The number of days PLT (Days PLT) was a key variable during the first 2 weeks, with a dichotomy in the mono-compartmental parameters for 0-3 Days PLT and 4-15 Days PLT. During 0-3 Days PLT, AST levels, indicative of allograft functionality (and TAC metabolism), were crucial predictors of elimination. Three groups were identified with the following clearances: CL/F0₋3 = 8.93 L/h for AST ≥ 500 U/L and CL/F0₋3 = 11.0 L/h for AST <500 U/L. During 4-15 Day PLT, low values of ALB (<2.5 g/dL) and HCT (<28 %) combined were determinant of a patient subgroup with a tendency to underexposure and complexity in empirical dose adjustment. The CL/F4₋15 = 25.1 L/h for this subgroup compared to CL/F4₋15 = 17.1 L/h for the others in that period. The elimination half-life for individual patients varied over tenfold so that a large number of subjects were not at steady state, making the use of a PK model necessary to achieve rapidly and safely the target concentration for TAC in LT. Validation of the model demonstrated that both bias and precision were within acceptable limits. CONCLUSION: For TAC therapy, covariate models using mixed effects methods are most useful when combined with patient-specific biochemical assays as well as clinical evidence. In such cases, the observed C(min) and Bayes methods can provide the most likely individual PK parameters, hence the optimal next dose to reach individualized target levels for each patient.

Pathophysiological idiosyncrasies and pharmacokinetic realities may interfere with tacrolimus dose titration in liver transplantation.

PURPOSE: To explore the main factors that make it difficult to empirically monitor tacrolimus (TAC) in the early period post-liver transplantation (LTx), with a specific focus on those aspects related to patient idiosyncrasy and clinical status as well as to the pharmacokinetic (PK) assumptions on which drug individualization in clinical practice is based. METHODS: Retrospective monitoring data from 75 de novo liver transplant patients treated with twice daily with TAC and followed for up to 15 days were analyzed. An extensive battery of laboratory measurements were available. Dose adjustment was performed empirically using trough levels (C(min)). The population was separated into two major background groups according to low or high values of aspartate aminotransferase (AST) (Group 1 and 2, respectively) based on AST measurements made during the first 4 days post-LTx. Each of these two major groups was then further subdivided into two subgroups based on elevated (Groups 1A, 2A) or reduced (Groups 1B, 2B) combined albumin (cut-off 2.5 g/dl) and hematocrit (cut-off 28%). RESULTS: The C(min)/Dose ratio [inversely proportional to systemic clearance (CL)] had a variability [coefficient of variation (CV) >80%) that was incongruently higher for the ratio than for C(min) and Dose separately. This was attributed to most patients not being at steady state or physiologically stable in the early post-LTx period. Group 1 patients were more predictable than Group 2 patients, who were responsible for the variability in the ratio. C(min) was lower in the reduced ALB and HCT patient groups when AST conditions were similar (1A vs. 1B and 2A vs. 2B), likely due to increased TAC metabolic clearance (reduced C(min)/Dose). This situation existed for two periods: 0-15 days post-LTx and 4-15 days post-LTx observations. Group 2A patients were the main source of the paradoxical variability in C(min)/Dose (higher ratio of 2.7; CV = 100%), suggesting a lower clearance and difficulty in the recovery of stability. In contrast, Group 2B patients had the lower ratio (1.4; 47%) but required the highest number of dose adjustments as the variability was hard to identify clinically. Group 1A patients were the most predictable empirically. When observations from 15 new patients who entered the clinic in 2007 and 2008 were used for the analysis, the same sub-groups existed in the same proportions in both years. CONCLUSION: The difficulty in empirical dose adjustment of TAC is associated to the inevitable non-fulfillment of PK assumptions early post-LTx and also to the inherent complexity of the clinical condition, leading to increased uncertainty for the clinician regarding dose selection. Identifying these sub-categories provides a rational means of classifying patients akin to a phenotype. The complexity of the kinetics in LTx and TAC treatment does not invalidate C(min) as a biomarker, but a Bayes algorithm including a full PK structure and these covariates would be optimal.

Simulation of sirolimus exposures and population variability immediately post renal transplantation: importance of the patient's CYP3A5 genotype in tailoring treatment.

The rapid achievement of efficacious exposure to sirolimus (SRL) after renal transplantation is crucial. However, there is high unpredictability in the dose to exposure relationship. Part of the variation is related to patients originating from subpopulations of fast or slow metabolizers via the CYP3A5*1/*3 genotype. The probability of achieving therapeutic SRL blood concentrations for each subpopulation under two equal-intensity increasing-frequency protocols after the start of treatment was explored with Monte Carlo simulation. The population pharmacokinetic model and inter-patient variability distributions of Djebli et al. (DRH2006) were sampled. They developed a base and final model with a genotype covariate for CL/F in patients receiving calcineurin inhibitor (CNI)-free therapy with SRL, mycophenolate mofetil and corticosteroids. Fast metabolizers (expressers) had a CL/F of 28.3 l/h whilst slow metabolizers (non-expressers) had a CL/F of 14.1 l/h. Here, in simulation, a standard 10 mg QD SRL was contrasted with a higher frequency of 5 mg BID SRL as related to the proportion of next dosed patients being within the 15-30 ng/ml trough levels on day 7 after transplantation. Near 0% of expressers on either regimen reached or exceeded the 30 ng/ml trough on day 7. Expressers showed protocol dependence for the chance of being within the 15-30 ng/ml range with the 5 mg BID protocol doubling those chances. Non-expressers appeared less protocol dependent for the probability of being above or below the 15-30 ng/ml range. The ability to determine the genotype early on may help to rationalize the initial titration of individual patients receiving CNI-free renal transplantation treatment with SRL.

Modeling methadone pharmacokinetics in rats in presence of P-glycoprotein inhibitor valspodar.

PURPOSE: To quantify the in vivo role of P-glycoprotein (P-gp) in the pharmacokinetics of methadone after intravenous and oral administration, using valspodar as a P-gp inhibitor. MATERIALS AND METHODS: Methadone plasma concentrations after intravenous (0.35 mg/kg) and oral (6 mg/kg) administration were analyzed, in absence and presence of valspodar, using nonlinear mixed effects modeling (NONMEM V). Non-parametric bootstrap analysis and posterior predictive check were employed as model evaluation techniques. RESULTS: The pharmacokinetics of methadone in the rat was successfully modeled using a two-compartmental model with a linear elimination from the central compartment and a first-order absorption process with lag time. Valspodar increased methadone F by 122% (95%CI: 34-269%) and decreased the V ( c ) and V ( p ) by 35% (95%CI: 16-49%) and 81% (95%CI: 63-93%), respectively. No effect of valspodar on other pharmacokinetic parameters was discernible. The non-parametric bootstrap analysis confirmed the absence of bias on the parameter estimates, and visual predictive check evidence the adequacy of the model to reproduce the observed time course of methadone plasma concentrations. CONCLUSION: Valspodar increased methadone's bioavailability as consequence of P-gp inhibition, which resulted in an increased analgesic effect of methadone.

The role of unbound drug in pharmacokinetics/pharmacodynamics and in therapy.

The evolution of research on drug protein binding is discussed with the unbound concentration (Cu) and the unbound fraction (fu) as protagonists. Particular attention is paid to the mechanisms via which alterations in binding affect the pharmacokinetics (PK) and the effect, or independently the pharmacodynamics (PD). Apart from albumin, the important alpha-acid glycoprotein (AGP), as well as specific drug classes and applications in the clinic and development (routine monitoring, cancer and HIV therapy, allometry) are addressed. The flaws with the classical method of indirectly calculating the Cu or the unbound PK/PD parameters, based on the fu in vitro, are related to the intrinsic complexity and variability in the outcomes. Increased focus is urged on directly estimating the unbound PK/PD and also on using population statistical methods.

Pharmacokinetic-pharmacodynamic modeling of the hydroxy lerisetron metabolite L6-OH in rats: an integrated parent-metabolite model.

PURPOSE: The twofold aim of this study was to characterize in vivo in rats the pharmacokinetics (PK) and pharmacodynamics (PD) of L6-OH, a metabolite of lerisetron with in vitro pharmacological activity, and evaluate the extent to which L6-OH contributes to the overall effect. METHODS: The PK of L6-OH was determined directly postmetabolite i.v. dose (PK-1), and also simultaneously for L (lerisetron concentration) and for generated L6-OH after lerisetron dose (200 microg kg(-1), i.v.), using Nonlinear Mixed Effects Modeling with an integrated parent-metabolite PK model (PK-2). Surrogate effect was measured by inhibition of serotonin-induced bradycardia. Protein binding was assayed via ultrafiltration and all quantification was performed via liquid chromatography-electrospray ionization-mass spectrometry. RESULTS: L6-OH showed elevated plasma and renal clearances, and volume of distribution (PK-1). The in vivo potency (PD) of L6-OH was high (EC(50) = 0.098 ng mL(-1) and EC(50unbound) = 0.040 ng mL(-1)). Total clearance for L (PK-2) in the presence of generated L6-OH (CL(L) = CL(-->L6-OH) + CL(n)) was 0.0139 L min(-1). Most of this clearance was L6-OH formation (F(c) = 99.6%), but only an 8.6% fraction of L6-OH was released into the bloodstream. The remainder undergoes biliar and fecal elimination. The parameters estimated from PK-2 were used to predict concentrations of L6-OH (Cp(L6)) generated after a lerisetron therapeutic dose (10 microg kg(-1)) in the rat. These concentrations are needed for the PD model and are below the quantification limit. Cp(L6max) was less than the EC(50) of L6-OH. CONCLUSIONS: We conclude that after lerisetron administration, L6-OH is extensively formed in the rat but it is quickly eliminated; therefore, besides being equipotent with the parent drug, the L6-OH metabolite does not influence the effect of lerisetron.

Bicompartmental kinetics of tobramycin analysed with a wide range of covariates.

The pharmacokinetics of tobramycin was studied in adult patients (N = 151) admitted either for initial suspicion of Gram-negative infection or for prophylaxis. In addition to age, weight, height and creatinine clearance (CrCL), a range of other covariates were also analysed, including type of pathology, co-medication, fever, sex and ethnicity (Basque or not). All patients received 100mg tobramycin every 8 h and samples were collected at three time points after the first dose and at two time points after the fourth dose and assayed with a fluorescence polarisation immunoassay. The population mixed effects bicompartmental parameters were obtained from 725 concentration measurements using NONMEM, FOCE method, and were: systemic clearance, CL = 6.03 L/h (between-subject coefficient of variation (CV) %, 29.4%); volume of distribution, V = 15.04 L (7.3%); and intercompartmental constants, k(12) = 0.192 h(-1) (56%) and k(21) = 0.55 h(-1) (no CV% determined). Covariate modelling was performed within NONMEM. Two alternative significant covariate models (Models 1 and 2) are proposed, with functions of CrCL and/or sex (Model 2). However, for clinical purposes, differentiation by sex is insignificant. Model 1 is for CL = 3.1 + 0.05.CrCLL/h (17.3%); V = 14.6 L (12%); k(12) = 0.224 h(-1) (63%) and k(21) = 0.468 h(-1). Stochastic simulation was used to predict the expected concentration 95th percentiles after the recommended 7 mg/kg dose and for minimum inhibitory concentration (MIC) = 1 mg/L, as well as alternative once-daily dosing regimens for MIC = 2 mg/L. It is seen that once-daily high-dose tobramycin is an appropriate strategy with respect to pharmacodynamic indices, C(peak)/MIC or AUC/MIC (where C(peak) is the peak plasma concentration and AUC is the area under the concentration-time curve).

Alpha-1-acid glycoprotein directly affects the pharmacokinetics and the analgesic effect of methadone in the rat beyond protein binding.

Methadone is a basic drug highly bound to alpha1-acid glycoprotein (AGP), a plasma protein that increases in several pathological situations. Our aims were to evaluate the processes (pharmacokinetics-PK and/or pharmacodynamics-PD) associated with changes of methadone analgesia under conditions of increased AGP, and whether these changes are dependent on binding, secondary to a pathology, or directly attributable to AGP. AGP levels, in rats, were increased by two different methods: (a) experimental inflammation with turpentine oil (TP), and (b) by directly infusing the protein (exo-AGP). Both had a corresponding control group. Tail-flick analgesia and PK were evaluated after methadone dose (0.35 mg/kg i.v.). Bicompartmental PK parameters as well as interanimal and assay variabilities were estimated using NONMEM. The relationship between Cp and analgesic effect (PD) was analyzed with WINNONLIN. AGP levels in both pretreated groups (TP and exo-AGP) were significantly increased, and the unbound fraction (fu) was decreased, compared to controls. All PK parameters were lower in the pretreated groups, but in exo-AGP the difference was maintained even when corrected by fu. Paradoxically, also in exo-AGP the analgesic effect was practically nonexistent, although the unbound Cp remained high, possibly associated to a change in the PD. AGP appears responsible for alterations in both PK and PD, beyond protein binding and inflammatory processes.

Effect of P-glycoprotein inhibition on methadone analgesia and brain distribution in the rat.

Methadone is an opiate drug that has been identified as an in-vitro substrate of the efflux pump P-glycoprotein (P-gp), active in the intestinal epithelium and in the blood-brain barrier (BBB), among other sites. The objective of this study was to test in vivo, in the rat model, the role of P-gp modulation on the analgesic effect and brain uptake of methadone, as well as identify the most relevant site via dual oral and intravenous (i.v.) experiments. The P-gp specific inhibitor (valspodar or PSC833) was preadministered (10 mg kg(-1) i.v.) to test groups. Analgesia was measured using the tailflick test. The ED50 for oral methadone (2, 3, 6 and 8 mg kg(-1)) decreased three-fold in valspodar groups compared with controls (2.23 +/- 0.002 mg kg(-1) and 6.07 +/- 0.07 mg kg(-1); P < 0.0001). The overall analgesic effect (% antinociception) was elevated 3.1 times in pretreated compared with control rats (90.65% +/- 0.22 vs 29.23% +/- 14.0; P < 0.01) after 6 mg kg(-1) oral methadone and 2.8 times after i.v. (0.35 mg kg(-1)) administration (91.75% +/- 4.27 vs 32.45% +/- 9.0; P < 0.01). The brain:plasma distribution ratio was higher in pretreated animals and AUCbrain (overall brain concentration) was 6 times higher after oral methadone and 4 times higher after i.v. compared with controls, disproportionally increased relative to plasma, implying an active process at the BBB. P-gp, and hence substrate comedication, plays a critical role in the evolution of the methadone analgesic effect and in its brain uptake, independent of the administration route.

Age-related changes in pharmacokinetics and pharmacodynamics of lerisetron in the rat: a population pharmacokinetic model.

BACKGROUND: The importance of studying the effects of age on the pharmacokinetics and pharmacodynamics of lerisetron - a new 5-hydroxytryptamine-3 (serotonin) receptor antagonist - comes from the facts that lerisetron will be administered to patients that are being treated with cytotoxic drugs and that the elderly frequently suffer from neoplastic diseases. OBJECTIVE: The present study was designed to explore the effects of age on the pharmacokinetics and pharmacodynamics of lerisetron by using an aged rat model. A mixed-effects population study was carried out in order to analyze the sparse data and to create covariate models which could be used to derive dosage recommendations. METHODS: Fischer 344 rats (n = 44) were divided into three groups, depending on their age: 5, 13, and 25 months. Blood samples were collected before administration of 200 micro g/kg of lerisetron for measurements of albumin, alpha(1)-acid glycoprotein, and unbound fraction of lerisetron. The lerisetron plasma concentrations were measured by high-performance liquid chromatography. A two-compartment model was fitted to the data using the nonlinear mixed-effects computer program WinNonMix. The population analysis was performed with the complete set of the collected data, and the potential sources of variability in the population parameters were investigated. Additionally, a pharmacodynamic study was performed. The effect of lerisetron (inhibition of the von Bezold-Jarisch reflex) was evaluated in young, adult, and senescent Fischer 344 rats. RESULTS: The mean values of the individual Bayes estimates of the parameters showed a decrease in total clearance and distribution volume of the central compartment in old rats. The lerisetron free (unbound) fraction remained unchanged among the groups, and there were no significant differences in alpha(1)-acid glycoprotein levels. The concentration-effect relationship was best described by a sigmoid E(max) model. Since the drug concentration in plasma at half-maximal effect (EC(50)) decreased in old rats, an increased sensitivity to the effect of lerisetron in old animals could be expected. CONCLUSION: Both pharmacokinetic changes (decreased volume of distribution and clearance and increased elimination half-life) and pharmacodynamic alterations (decrease in total and unbound EC(50)) may be responsible for the different responses to lerisetron observed in old rats.

Population pharmacokinetics of netilmicin in short-term prophylactic treatment.

AIMS: To characterize the population pharmacokinetics of netilmicin, an aminoglycoside antibiotic, in adult urology patients and to develop a covariate model for improved dose titration. METHODS: Data from 62 adult patients (55 male, seven female), undergoing urological surgery and treated with netilmicin for short-term prophylaxis, were evaluated retrospectively. The group had (median, range) ages 68, 31-92 years, weights 72, 43-106 kg and heights 167, 148-182 cm. No patient showed renal impairment before netilmicin treatment (serum creatinine

'In vitro' binding of propofol to serum lipoproteins in thyroid dysfunction.

OBJECTIVE: To determine a regression relationship between the unbound fraction (fu percent) of propofol, a highly lipophilic intravenous anaesthetic agent, and demographic and biochemical variables in a thyroid dysfunction population. METHODS: Serum samples from patients with hypo- (n=33) and hyperthyroidism (n=33) and also from healthy volunteers (control group; n=9) were spiked with propofol to a total (bound + unbound propofol) concentration of 10 microg/ml. The unbound concentration was determined using ultrafiltration followed by high-performance liquid chromatography with fluorimetric detection and the unbound fraction percent (fu) and the binding ratio [bound/free concentration (B/F)] were calculated. Albumin, alpha(1)-acid glycoprotein, cholesterol, triglycerides, HDL, LDL, VLDL lipoproteins, and T3 and T4 hormone concentrations were measured. B/F has a linear relationship with lipoprotein concentration--unlike that for fu which is curvilinear - so, we developed a linear regression model of B/F for propofol with the above biochemical variables and with gender. RESULTS: The fu in hypothyroidism was significantly lower than in the healthy control (mean +/- standard deviation 0.74 +/- 0.13% vs 0.87 +/- 0.12%, P < 0.01) but the fu was no different in hyperthyroid subjects than controls (0.94 +/- 0.16% vs 0.87 +/- 0.12%, P > 0.05). HDL, LDL and VLDL lipoproteins were significant predictors of B/F (P < 0.05) and were capable of explaining 54% of the variability in the 'in vitro' binding of propofol in thyroid disease. CONCLUSION: These results could help explain the interindividual variability in the protein binding of propofol in thyroid dysfunction patients and suggest the inclusion of lipoproteins in covariate model development for the pharmacokinetic/pharmacodynamic parameters of propofol.

Prediction of unbound propofol concentrations in a diabetic population.

Propofol is a short-acting general intravenous anesthetic characterized by a wide interindividual variability in the response after the same dose. Its binding to serum proteins exceeds 98%, so small changes in protein concentrations can be amplified in the unbound fraction of the drug and hence possibly in the effect. It is then likely that part of the variability in the response could be attributed to differences in protein levels among individuals and particularly among those with pathologies such as diabetes. The aim of this study was to establish predictive regression models in a diabetes mellitus (DM) population between unbound:bound propofol ratios and demographic and biochemical indices. Unbound:bound propofol ratios can be routinely obtained in the clinic as opposed to the free fraction of the drug. In DM patients (30 women and 37 men aged between 17 and 78 y) with mellitus type 1 (n = 37) and type 2 (n = 30) diabetes, the authors measured the lipoproteins (HDL, LDL, and VLDL), cholesterol, triglycerides, albumin, alpha1-acid glycoprotein (AAG), free fatty acids (FFA), glycosylated hemoglobin, and the unbound fraction (Fu) and the bound/free ratio (B/F) of propofol. A linearized regression model between the above variables--as well as age, sex, and type of diabetes--and Fu was then developed. Patients had blood drawn and sera separated by centrifugation and spiked with propofol to a concentration of 10 microg/mL. The Fu was determined via ultrafiltration. Multiple linear regression analysis was used to identify significant predictor variables of Fu in this population and two models were originated: one with lipoprotein serum concentrations as explanatory variables (Model A) and another that depended on cholesterol and triglycerides (Model B). Both models presented high correlation coefficients (r2 = 0.71 and 0.68, respectively; P < 0.0001), and each was used to predict Fu in an independent group of 15 DM patients of similar characteristics and biochemical indices as the model development group. Bias and precision were for Model A, 0.9% and 7.8%, and for Model B, 3.0% and 8.7%, respectively. Both models were compared with each other and to a naive predictor (the mean) and each was better than the naive model in predicting the unbound fraction of propofol. Model A and model B could be used in estimating Fu of propofol in DM patients based on the more routine clinical measures of lipoprotein serum concentrations or cholesterol and triglyceride levels.

Sex specificity in methadone analgesia in the rat: a population pharmacokinetic and pharmacodynamic approach.

PURPOSE: To quantify the extent to which a sex-specific dichotomy in the temporal evolution of the analgesic effect, after intravenous (i.v.) methadone injection in the rat, relates to the pharmacokinetics (PK) and pharmacodynamics (PD) that mediate the dose-to-effect pathway. METHODS: Tail-flick analgesia was measured after i.v. methadone injection (0.35 mg/kg) in female (n = 16) and male (n = 16) Sprague-Dawley rats. The PK were evaluated in separate female (n = 56) and male (n = 56) rats after they had received the same dose of methadone i.v. (0.35 mg/kg). A bicompartmental model described the kinetics and a sigmoid Emax model-related drug effect vs. simulated concentrations (pharmacodynamics) at the times of effect measurement. All model parameters as well as interanimal and assay variabilities were estimated with a mixed-effects population method using the program NONMEM. RESULTS: The area under the effect-time curve (AUCE0-120) was (mean +/- interanimal SD) 1859+/-346 min in the females, which was significantly lower than the 4871+/-393 min in the males (P < 0.0001). On the contrary, the profiles of concentration vs. time were higher in females and, therefore, corresponded inversely to the effect vs. time-relative magnitudes. The central volume of distribution, V1, was 1.94+/-0.37 l/kg for female rats and 3.01+/-0.33 l/kg for male rats. Also, the central clearance was 0.077+/-0.006 l/min/kg and 0.102+/-0.005 l/min/ kg, respectively, for female and male rats. Both parameters differed significantly between sexes (P < 0.0001). The pharmacodynamic maximum observed effect parameter (Emax) was 37%+/-29% in female rats and 85%+/-16% in male rats, and these values were significantly different (P < 0.0001). The parameter for the concentration eliciting half of Emax (EC50) was 24.1+/-7.5 microg/l in female rats and 20.3+/-2.9 microg/l in male rats, and the Hill-related exponent, gamma, was 6.3+/-3.9 in female rats and 5.5+/-4.1 in male rats. These parameters did not differ significantly (at the P < 0.05 level). CONCLUSIONS: A sex-specific dichotomy in the methadone antinociceptive effect, in the rat, was not proportionally related to plasma concentrations. Each sex corresponded to a distinct subpopulation of the PK parameters and one of the pharmacodynamic parameters (Emax). When the course of a drug involves PK or PD subpopulations, PK/PD modeling can afford the safest prediction of the effect-time evolution for a particular dose.

Effect of omeprazole on oral and intravenous RS-methadone pharmacokinetics and pharmacodynamics in the rat.

The effect of omeprazole on oral and intravenous (iv) RS-methadone pharmacokinetics and pharmacodynamics was studied in awake, freely moving rats, which were divided in four groups: oral RS-methadone (3 mg/kg) was given (a) to a control group (CO(oral)) (n = 65) and (b) to an omeprazole pretreated group (OP(oral)) (n = 77), and iv RS-methadone (0.35 mg/kg) was administered (c) to a control group (CO(iv)) (n = 86) and (d) to an omeprazole pretreated group (OP(iv)) (n = 86). Omeprazole (2 mg/kg) was given iv 2 h before RS-methadone. Plasma concentrations of RS-methadone (Cp) were determined by high-performance liquid chromatography and analgesic response by tail flick for 0-180 min (oral) and 0-120 min (iv). RS-Methadone rate of absorption (mean +/- SE) was faster in OP(oral) (k(01) = 0.31 +/- 0.04 min(-1)) than in CO(oral) (k(01) = 0.05 +/- 0.007 min(-1)), consequently plasma peak concentrations (C(max)) were greater (197.41 +/- 33.70 ng/mL versus 83.54 +/- 7.97 ng/mL) and the time to reach C(max) (t(max)) was shorter (11.23 +/- 1.32 min versus 39.18 +/- 1.74 min). Mean area under the Cp-time curve (AUC(0-infinity)) and hence bioavailability of oral RS-methadone were increased by omeprazole without significant changes in the elimination. Omeprazole did not affect the pharmacokinetics of iv RS-methadone. The changes of the analgesic effect of RS-methadone as a function of time were similar in all four groups. In the CO(oral) group, Cp and analgesic effect were defined by the E(max) model. The relationship between Cp and drug effect in the OP(oral) group showed a counterclockwise hysteresis (k(e0) of 0.018 +/- 0.006 min(-1)). For the iv groups (CO(iv) and OP(iv)), the Cp-analgesic effect relationship was described by an E(max) sigmoid model and omeprazole did not affect the pharmacodynamic parameters. It is concluded that omeprazole causes an increase in the bioavailability of oral RS-methadone without modifying the analgesic response but affecting the Cp-effect relationship.

Pharmacokinetics and pharmacological effect of lerisetron, a new 5-HT3 antagonist, in rats.

The pharmacokinetics of lerisetron, a novel 5HT(3) antagonist, are studied together with its efficacy in inhibiting the serotonin (5-HT)-evoked transient bradycardia reflex (von Bezold-Jarisch reflex) in Sprague Dawley rats. [(14)C]Lerisetron (50, 100, and 200 microg/kg) was given to rats by intravenous (iv) injection and plasma levels of unchanged (UL) and total (unchanged + changed, TL) drug were measured by liquid chromatography with radioactivity monitoring and scintillation counting, respectively. Linearity of UL and TL pharmacokinetics over the dose range was established by noncompartmental analysis. Protein binding of lerisetron was measured in vitro by ultrafiltration. The unbound fraction was 14.4 +/- 1.4%. A nonlinear mixed effects ("population") bicompartmental pharmacokinetic analysis showed that volume of distribution and clearance (CL) were high for both forms of the drug, but CL was significantly smaller for TL [(mean +/- SEM) 0.014 +/- 0.03 L/min for UL versus 0.006 +/- 0.03 L/min for TL, p < 0.05)]. Large interindividual variabilities were observed for both forms. The response to lerisetron administration (inhibition of bradycardia) was evaluated at different doses (2, 3, 5, 6, and 10 microg/kg, iv) at times 2-180 min after dose administration and related to simulated concentrations. Inhibition was 100% 5 min after the 10-microg/kg dose and, 3 h later, it was still > 10%. Response to lerisetron was dose related in the range studied. Pharmacodynamic parameters were estimated by a sigmoid E(max) naive-pooled model. The parameters were also different between the two forms: EC(50) was 0.44 ng/mL (CV = 5.9%) for UL and 0.88 ng/mL (CV = 4.9%) for TL. We conclude that UL and TL pharmacokinetics were linear and that the differences in the kinetics and dynamics between the two forms suggest the presence of at least one metabolite.