A pharmacokinetic-pharmacodynamic model for duodenal levodopa infusion

Objective Levodopa in presence of decarboxylase inhibitors is following two-compartment kinetics and its effect is typically modelled using sigmoid Emax models. Pharmacokinetic modelling of the absorption phase of oral distributions is problematic because of irregular gastric emptying. The purpose of this work was to identify and estimate a population pharmacokinetic- pharmacodynamic model for duodenal infusion of levodopa/carbidopa (Duodopa®) that can be used for in numero simulation of treatment strategies. Methods The modelling involved pooling data from two studies and fixing some parameters to values found in literature (Chan et al. J Pharmacokinet Pharmacodyn. 2005 Aug;32(3-4):307-31). The first study involved 12 patients on 3 occasions and is described in Nyholm et al. Clinical Neuropharmacology 2003:26:156-63. The second study, PEDAL, involved 3 patients on 2 occasions. A bolus dose (normal morning dose plus 50%) was given after a washout during night. Plasma samples and motor ratings (clinical assessment of motor function from video recordings on a treatment response scale between -3 and 3, where -3 represents severe parkinsonism and 3 represents severe dyskinesia.) were repeatedly collected until the clinical effect was back at baseline. At this point, the usual infusion rate was started and sampling continued for another two hours. Different structural absorption models and effect models were evaluated using the value of the objective function in the NONMEM package. Population mean parameter values, standard error of estimates (SE) and if possible, interindividual/interoccasion variability (IIV/IOV) were estimated. Results Our results indicate that Duodopa absorption can be modelled with an absorption compartment with an added bioavailability fraction and a lag time. The most successful effect model was of sigmoid Emax type with a steep Hill coefficient and an effect compartment delay. Estimated parameter values are presented in the table. Conclusions The absorption and effect models were reasonably successful in fitting observed data and can be used in simulation experiments.

Pharmacokinetic and safety evaluation of the use of ciprofloxacin on an isoniazid-rifampicin regimen in rabbits

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Pharmacokinetic Profile of A New Diclofenac Prodrug without Gastroulcerogenic Effect

Gastrotoxicity is a major problem for long-term therapy with non-steroidal anti-inflammatory drugs (NSAIDs). DICCIC (1-(2,6-dichlorophenyl)indolin-2-one) is a new diclofenac prodrug, which has proven anti-inflammatory activity without gastroulcerogenic effect. The aim of this work was to compare the pharmacokinetic profiles of diclofenac from DICCIC (7.6 mg/kg equivalent to 8.1 mg/kg diclofenac) and diclofenac (8.1 mg/kg) administration in Wistar rats weighing 250-300 g (n=20). The doses were calculated by interspecific allometric scaling based on the 2 mg/kg from diary human dose of diclofenac. Blood samples were collected in heparinized tubes via the femoral artery through the implanted catheter. The plasma was separated and quantitation was made in a HPLC system with a UV-Vis detector. The confidence limits of the bioanalytical method were appropriate for its application in a preclinical pharmacokinetic study. The AUC of diclofenac from DICCIC (53.7± 5.8 ug/mL.min) was significantly less (Mann Whitney test, p<0.05) than that of diclofenac from diclofenac administration (885.9 ± 124,8 ug/mL.min). Terminal half-life of diclofenac from DICCIC (50.1 ± 17.2 min) was significantly less (Mann Whitney test, p<0.05) than that of diclofenac from diclofenac administration (247.4 ± 100.9 min). Still the parameters clearance and distribution volume were calculated for diclofenac from diclofenac, whose results were 9.2 ±1.2 mL/min.kg and 3.3 ±1.2 L/kg, respectively. The results of DICCIC from DICCIC administration were 108.9 ± 19.6 mL/min.kg and 7.8 ± 2.4 L/kg for clearance and distribution volume, respectively. The pharmacokinetic profile demonstrated that there was an increase in diclofenac elimination and a lower exposure to diclofenac with administration of DICCIC compared to diclofenac. © 2013 Bentham Science Publishers.

The pharmacokinetic profile of fesoterodine: similarities and differences to tolterodine

BACKGROUND: Fesoterodine is a new antimuscarinic agent developed for the treatment of overactive bladder. Fesoterodine itself is inactive and is rapidly and extensively converted by ubiquitous esterases to its principal active moiety, 5-hydroxymethyl tolterodine (5-HMT). 5-HMT is formed via biotransformation of both fesoterodine and tolterodine, albeit by different metabolising enzymes, viz. esterases and CYP2D6 respectively. Tolterodine is a potent muscarinic receptor antagonist and has been used for the treatment of overactive bladder for over ten years. The objective of this study was to establish the pharmacokinetic profile of fesoterodine and to highlight ist potential pharmacokinetic advantages over tolterodine. DESIGN: Single-centre, open-label, randomised, 4-way crossover study in a total of 24 healthy male volunteers. Single oral doses of 4, 8, or 12 mg fesoterodine were administered after an overnight fast. In addition, the 8 mg dose was also administered after a standard high-fat and high-calorie breakfast. Blood and urine samples for the analysis of 5-HMT were collected before and multiple times after drug administration for pharmacokinetic analysis. RESULTS: The mean peak plasma concentration (Cmax) of 5-HMT and the mean area under the time versus concentration curve (AUC) increased proportionally with the fesoterodine dose. These two parameters were some 2-fold higher in CYP2D6 poor metabolisers, whereas the time to peak plasma concentration (tmax) and half life (t1/2) were not influenced by the dose or the CYP2D6 metaboliser status. If fesoterodine was taken following a high-fat breakfast, we observed small increases in Cmax and AUC. In spite of these modest genetic influences and food effects on the pharmacokinetics of fesoterodine, the overall interindividual variability in Cmax levels was relatively little compared to previously published reports using tolterodine. CONCLUSIONS: Due to the esterase-mediated cytochrome P450-independent formation of 5-HMT and involvement of multiple metabolic and renal excretion pathways in the elimination of 5-HMT, the effects of patient-intrinsic and -extrinsic factors on the pharmacokinetics of fesoterodine are only modest, with some 2-fold higher 5-HMT exposure. Therefore, in contrast to tolterodine, no reduction of fesoterodine dosage is required under conditions of reduced elimination. In most cases of drug interaction or renal/hepatic impairment, the fesoterodine dose may be increased to 8 mg/day based on individual patients' response, or patients may be required to remain at the initial recommended dose of 4 mg/day.

DCE@urLAB: a dynamic contrast-enhanced MRI pharmacokinetic analysis tool for preclinical data

Background DCE@urLAB is a software application for analysis of dynamic contrast-enhanced magnetic resonance imaging data (DCE-MRI). The tool incorporates a friendly graphical user interface (GUI) to interactively select and analyze a region of interest (ROI) within the image set, taking into account the tissue concentration of the contrast agent (CA) and its effect on pixel intensity. Results Pixel-wise model-based quantitative parameters are estimated by fitting DCE-MRI data to several pharmacokinetic models using the Levenberg-Marquardt algorithm (LMA). DCE@urLAB also includes the semi-quantitative parametric and heuristic analysis approaches commonly used in practice. This software application has been programmed in the Interactive Data Language (IDL) and tested both with publicly available simulated data and preclinical studies from tumor-bearing mouse brains. Conclusions A user-friendly solution for applying pharmacokinetic and non-quantitative analysis DCE-MRI in preclinical studies has been implemented and tested. The proposed tool has been specially designed for easy selection of multi-pixel ROIs. A public release of DCE@urLAB, together with the open source code and sample datasets, is available at http://www.die.upm.es/im/archives/DCEurLAB/ webcite.

What is the best size descriptor to use for pharmacokinetic studies in the obese?

The prevalence of obesity in the western world is dramatically rising, with many of these individuals requiring therapeutic intervention for a variety of disease states. Despite the growing prevalence of obesity there is a paucity of information describing how doses should be adjusted, or indeed whether they need to be adjusted, in the clinical setting. This review is aimed at identifying which descriptors of body size provide the most information about the relationship between dose and concentration in the obese. The size descriptors, weight, lean body weight, ideal body weight, body surface area, body mass index, fat-free mass, percent ideal body weight, adjusted body weight and predicted normal body weight were considered as potential size descriptors. We conducted an extensive review of the literature to identify studies that have assessed the quantitative relationship between the parameters clearance (CL) and volume of distribution (V) and these descriptors of body size. Surprisingly few studies have addressed the relationship between obesity and CL or V in a quantitative manner. Despite the lack of studies there were consistent findings: (i) most studies found total body weight to be the best descriptor of V. A further analysis of the studies that have addressed V found that total body weight or another descriptor that incorporated fat mass was the preferred descriptor for drugs that have high lipophilicity; (ii) in contrast, CL was best described by lean body mass and no apparent relationship between lipophilicity or clearance mechanism and preference for body size descriptor was found. In conclusion, no single descriptor described the influence of body size on both CL and V equally well. For drugs that are dosed chronically, and therefore CL is of primary concern, dosing for obese patients should not be based on their total weight. If a weight-based dose individualization is required then we would suggest that chronic drug dosing in the obese subject should be based on lean body weight, at least until a more robust size descriptor becomes available.

Some considerations on the design of population pharmacokinetic studies

The goal of this manuscript is to introduce a framework for consideration of designs for population pharmacokinetic orpharmacokinetic-pharmacodynamic studies. A standard one compartment pharmacokinetic model with first-order input and elimination is considered. A series of theoretical designs are considered that explore the influence of optimizing the allocation of sampling times, allocating patients to elementary designs, consideration of sparse sampling and unbalanced designs and also the influence of single vs. multiple dose designs. It was found that what appears to be relatively sparse sampling (less blood samples per patient than the number of fixed effects parameters to estimate) can also be highly informative. Overall, it is evident that exploring the population design space can yield many parsimonious designs that are efficient for parameter estimation and that may not otherwise have been considered without the aid of optimal design theory.

Bayesian population pharmacokinetic analysis of sirolimus

Objective: It is usual that data collected from routine clinical care is sparse and unable to support the more complex pharmacokinetic (PK) models that may have been reported in previous rich data studies. Informative priors may be a pre-requisite for model development. The aim of this study was to estimate the population PK parameters of sirolimus using a fully Bayesian approach with informative priors. Methods: Informative priors including prior mean and precision of the prior mean were elicited from previous published studies using a meta-analytic technique. Precision of between-subject variability was determined by simulations from a Wishart distribution using MATLAB (version 6.5). Concentration-time data of sirolimus retrospectively collected from kidney transplant patients were analysed using WinBUGS (version 1.3). The candidate models were either one- or two-compartment with first order absorption and first order elimination. Model discrimination was based on computation of the posterior odds supporting the model. Results: A total of 315 concentration-time points were obtained from 25 patients. Most data were clustered at trough concentrations with range of 1.6 to 77 hours post-dose. Using informative priors, either a one- or two-compartment model could be used to describe the data. When a one-compartment model was applied, information was gained from the data for the value of apparent clearance (CL/F = 18.5 L/h), and apparent volume of distribution (V/F = 1406 L) but no information was gained about the absorption rate constant (ka). When a two-compartment model was fitted to the data, the data were informative about CL/F, apparent inter-compartmental clearance, and apparent volume of distribution of the peripheral compartment (13.2 L/h, 20.8 L/h, and 579 L, respectively). The posterior distribution of the volume distribution of central compartment and ka were the same as priors. The posterior odds for the two-compartment model was 8.1, indicating the data supported the two-compartment model. Conclusion: The use of informative priors supported the choice of a more complex and informative model that would otherwise have not been supported by the sparse data.

A D-optimal designed population pharmacokinetic study of itraconazole capsules and solution in adults with cystic fibrosis

Background: Oral itraconazole (ITRA) is used for the treatment of allergic bronchopulmonary aspergillosis in patients with cystic fibrosis (CF) because of its antifungal activity against Aspergillus species. ITRA has an active hydroxy-metabolite (OH-ITRA) which has similar antifungal activity. ITRA is a highly lipophilic drug which is available in two different oral formulations, a capsule and an oral solution. It is reported that the oral solution has a 60% higher relative bioavailability. The influence of altered gastric physiology associated with CF on the pharmacokinetics (PK) of ITRA and its metabolite has not been previously evaluated. Objectives: 1) To estimate the population (pop) PK parameters for ITRA and its active metabolite OH-ITRA including relative bioavailability of the parent after administration of the parent by both capsule and solution and 2) to assess the performance of the optimal design. Methods: The study was a cross-over design in which 30 patients received the capsule on the first occasion and 3 days later the solution formulation. The design was constrained to have a maximum of 4 blood samples per occasion for estimation of the popPK of both ITRA and OH-ITRA. The sampling times for the population model were optimized previously using POPT v.2.0.[1] POPT is a series of applications that run under MATLAB and provide an evaluation of the information matrix for a nonlinear mixed effects model given a particular design. In addition it can be used to optimize the design based on evaluation of the determinant of the information matrix. The model details for the design were based on prior information obtained from the literature, which suggested that ITRA may have either linear or non-linear elimination. The optimal sampling times were evaluated to provide information for both competing models for the parent and metabolite and for both capsule and solution simultaneously. Blood samples were assayed by validated HPLC.[2] PopPK modelling was performed using FOCE with interaction under NONMEM, version 5 (level 1.1; GloboMax LLC, Hanover, MD, USA). The PK of ITRA and OH‑ITRA was modelled simultaneously using ADVAN 5. Subsequently three methods were assessed for modelling concentrations less than the LOD (limit of detection). These methods (corresponding to methods 5, 6 & 4 from Beal[3], respectively) were (a) where all values less than LOD were assigned to half of LOD, (b) where the closest missing value that is less than LOD was assigned to half the LOD and all previous (if during absorption) or subsequent (if during elimination) missing samples were deleted, and (c) where the contribution of the expectation of each missing concentration to the likelihood is estimated. The LOD was 0.04 mg/L. The final model evaluation was performed via bootstrap with re-sampling and a visual predictive check. The optimal design and the sampling windows of the study were evaluated for execution errors and for agreement between the observed and predicted standard errors. Dosing regimens were simulated for the capsules and the oral solution to assess their ability to achieve ITRA target trough concentration (Cmin,ss of 0.5-2 mg/L) or a combined Cmin,ss for ITRA and OH-ITRA above 1.5mg/L. Results and Discussion: A total of 241 blood samples were collected and analysed, 94% of them were taken within the defined optimal sampling windows, of which 31% where taken within 5 min of the exact optimal times. Forty six per cent of the ITRA values and 28% of the OH-ITRA values were below LOD. The entire profile after administration of the capsule for five patients was below LOD and therefore the data from this occasion was omitted from estimation. A 2-compartment model with 1st order absorption and elimination best described ITRA PK, with 1st order metabolism of the parent to OH-ITRA. For ITRA the clearance (ClItra/F) was 31.5 L/h; apparent volumes of central and peripheral compartments were 56.7 L and 2090 L, respectively. Absorption rate constants for capsule (kacap) and solution (kasol) were 0.0315 h-1 and 0.125 h-1, respectively. Comparative bioavailability of the capsule was 0.82. There was no evidence of nonlinearity in the popPK of ITRA. No screened covariate significantly improved the fit to the data. The results of the parameter estimates from the final model were comparable between the different methods for accounting for missing data, (M4,5,6)[3] and provided similar parameter estimates. The prospective application of an optimal design was found to be successful. Due to the sampling windows, most of the samples could be collected within the daily hospital routine, but still at times that were near optimal for estimating the popPK parameters. The final model was one of the potential competing models considered in the original design. The asymptotic standard errors provided by NONMEM for the final model and empirical values from bootstrap were similar in magnitude to those predicted from the Fisher Information matrix associated with the D-optimal design. Simulations from the final model showed that the current dosing regimen of 200 mg twice daily (bd) would provide a target Cmin,ss (0.5-2 mg/L) for only 35% of patients when administered as the solution and 31% when administered as capsules. The optimal dosing schedule was 500mg bd for both formulations. The target success for this dosing regimen was 87% for the solution with an NNT=4 compared to capsules. This means, for every 4 patients treated with the solution one additional patient will achieve a target success compared to capsule but at an additional cost of AUD $220 per day. The therapeutic target however is still doubtful and potential risks of these dosing schedules need to be assessed on an individual basis. Conclusion: A model was developed which described the popPK of ITRA and its main active metabolite OH-ITRA in adult CF after administration of both capsule and solution. The relative bioavailability of ITRA from the capsule was 82% that of the solution, but considerably more variable. To incorporate missing data, using the simple Beal method 5 (using half LOD for all samples below LOD) provided comparable results to the more complex but theoretically better Beal method 4 (integration method). The optimal sparse design performed well for estimation of model parameters and provided a good fit to the data.

Development of a physiologically-based pharmacokinetic model of the rat central nervous system

Central nervous system (CNS) drug disposition is dictated by a drug’s physicochemical properties and its ability to permeate physiological barriers. The blood–brain barrier (BBB), blood-cerebrospinal fluid barrier and centrally located drug transporter proteins influence drug disposition within the central nervous system. Attainment of adequate brain-to-plasma and cerebrospinal fluid-to-plasma partitioning is important in determining the efficacy of centrally acting therapeutics. We have developed a physiologically-based pharmacokinetic model of the rat CNS which incorporates brain interstitial fluid (ISF), choroidal epithelial and total cerebrospinal fluid (CSF) compartments and accurately predicts CNS pharmacokinetics. The model yielded reasonable predictions of unbound brain-to-plasma partition ratio (Kpuu,brain) and CSF:plasma ratio (CSF:Plasmau) using a series of in vitro permeability and unbound fraction parameters. When using in vitro permeability data obtained from L-mdr1a cells to estimate rat in vivo permeability, the model successfully predicted, to within 4-fold, Kpuu,brain and CSF:Plasmau for 81.5% of compounds simulated. The model presented allows for simultaneous simulation and analysis of both brain biophase and CSF to accurately predict CNS pharmacokinetics from preclinical drug parameters routinely available during discovery and development pathways.