Search

The physiological changes following Roux-en-Y gastric bypass (RYGB) surgery may impact drug release from mechanistically different controlled-release tablets, making generic substitution inappropriate. This study aimed to characterise the pharmacokinetic-pharmacodynamic relationships of oxycodone from a lipid-based and water-swellable controlled-release tablet in RYGB patients. Twenty RYGB patients received 10-mg oral solution oxycodone or 20-mg controlled-release (water-swellable or lipid-based) oxycodone in a three-way, randomised, semiblinded and cross-over study. Blood sampling and pupillary recordings were conducted over a 24-h period. A previously established pharmacokinetic-pharmacodynamic model of these three formulations in healthy volunteers was used in the analysis as a reference model. No differences in absorption kinetics were seen between controlled-release formulations in patients. However, the absorption lag time was 11.5 min in patients vs 14 min in healthy volunteers for controlled-release tablets (P

Physiologically based pharmacokinetic (PBPK) models allow integrated representations of the time course of drug concentrations in the important organs and sites of drug action, toxicity, absorption, metabolism, and excretion. PBPK models have multiple compartments representing defined organs or physiological spaces, with parameters that are potentially directly measurable. They lie on a continuum of pharmacokinetic model types, which includes Compartmental (least complex), Semi-PBPK, PBPK, and Systems Biology (most complex) models. Compartmental models are typically parameterized using a “top-down” approach, where the model is fitted to a specific data set to estimate parameter values. PBPK models are typically parameterized using a “bottom-up” approach, where model parameters are derived from literature, in vitro data, or scaled data from another species. Although there is increasing availability and use of commercial software for PBPK modeling, it is feasible to construct bespoke PBPK for specific projects using general-purpose software platforms. The challenge here is often collecting, collating, and justifying the data used to parameterize the model. In this chapter, the fundamental equations for a simple PBPK models are presented with respect to key “submodels” of an example whole-body model. The current rise in the rate of publication and application of PBPK models is likely to be sustained in the foreseeable future. Contributing factors will likely be the wider use of commercial PBPK models, the advent of broader collective efforts to advance and coordinate PBPK and Systems Biology modeling, and the increasing ease with which bespoke PBPK models can be coded and shared

Oral controlled-release formulations are playing an ever-increasing role in opioid therapy; however, little is known about their influence on the relationship between pharmacokinetics and pharmacodynamics. The study aim was to characterize the pharmacokinetic-pharmacodynamics of two controlled-release tablet formulations and a liquid formulation of oxycodone in healthy, opioid-naïve volunteers, which can serve as a reference for future patient studies. A semi-double-blinded, three-way crossover study was conducted, with fifteen healthy volunteers receiving two differently designed 20 mg monophasic controlled-release oxycodone tablets and 10 mg oral solution oxycodone in a randomized order. Venous plasma concentrations and pupil diameter were determined pre-dose and 0.25, 0.5, 0.75, 1, 1.5, 2, 2.33, 2.66, 3, 3.33, 3.66, 4, 5, 6, 8, 12 and 24 hour post-dose. Oxycodone pharmacokinetics was best described by a two-compartment model with first-order absorption. The controlled-release formulations had an absorption lag of 0.23 hour and a slower absorption rate constant (k aCR = 0.19 hour -1 ) compared to the oral solution (k aSOL = 0.94 hour -1 ). Effects on pupil diameter were delayed relative to plasma (14 minutes half-life) for all formulations and were best described by a proportional E max model. The plasma concentration of oxycodone at half-maximum effect was lower in males (31.1 μg/L) compared to females (52.8 μg/L; P < .001). The absorption profile of controlled-release oxycodone formulations provided a prolonged onset and offset of action compared to oral solution oxycodone. The controlled-release formulations showed no differences in pharmacokinetic and pharmacodynamic parameters suggesting that both may be used interchangeably in human beings with normal gastrointestinal function.

The aim of this study was to develop a population in vitro–in vivo pharmacokinetic model that simultaneously describe the absorption and accumulation kinetics of itraconazole (ICZ) and hydroxy-itraconazole (HICZ) in healthy subjects. The model integrated meta-models of gastrointestinal pH and gastrointestinal transit time and in vitro dissolution models of ICZ with the absorption and disposition kinetics of ICZ and HICZ. Mean concentration intravenous data, and single- and multi-dose oral data were used for model development. Model development was conducted in NONMEM in a stepwise manner. First, a model of intravenous data (systemic kinetics) was established and then extended to include the oral data. The latter was then extended to establish the in vitro–in vivo pharmacokinetic model. The systemic disposition of ICZ was best described by a 3-compartment model with oral absorption described by 4-transit compartments and HICZ distribution by a 1-compartment model. ICZ clearance was best described using a mixed inhibition model that allowed HICZ concentrations to inhibit the clearance of parent drug. HICZ clearance was described by Michaelis–Menten elimination kinetics. An in vitro–in vivo model was successfully established for both formulations. The presented model was able to describe ICZ and HICZ plasma concentrations over a wide range of oral and intravenous doses and allowed the exploration of complexities associated with the non-linear ICZ and HICZ kinetics. The model may provide insight into the variability in exposure of ICZ with respect to relating in vivo dissolution characteristics with in vivo disposition kinetics. Refereed/Peer-reviewed

Exposure-response (ER) modeling for fixed-dose combinations (FDC) has previously been found to have an inflated false positive rate (FP), i.e., observing a significant effect of FDC components when no true effect exists. Longitudinal exposure-response (LER)analysis utilizes the time course of the data and is valid for several clinical endpoints for FDCs. The aim of the study was to investigate if LER is applicable for the validation of FDCs by demonstrating the contribution of each component to the overall effect without inflation of FP rates. FP and FN rates associated with ER and LER analysis were investigated using stochastic simulation and estimation. Four hundred thirty-two scenarios with varying numbers of patients, duration, sampling frequency, dose distribution, design, and drug activity were analyzed using a range of linear, log-linear, and non-linear models to asses FP and FN rates. Lastly, the impact of the clinical trial parameters was investigated. LER analyses provided well-controlled FP rates of the expected 5% or less; however, in low information clinical trials consisting of 30 patients, 4 samples, and 20 days, LER analyses lead to inflated FN rates. Parameter investigation showed that when the clinical trial includes sufficient patients, duration, samples, and an appropriate trial design, the FN rates are in general below the expected 5% for LER analysis. Based on the results, LER analysis can be used for the validation of FDCs and fixed ratio drug combinations. The method constitutes anew avenue for providing evidence that demonstrates the contribution of each component to the overall clinical effect. Refereed/Peer-reviewed

Purpose: Lenalidomide is used widely in B-cell malignancies for its immunomodulatory activity. It is primarily eliminated via the kidneys, with a significant proportion of renal elimination attributed to active processes. Lenalidomide is a weak substrate of P-glycoprotein (P-gp), though it is unclear whether P-gp is solely responsible for lenalidomide transport. This study aimed to determine whether the current knowledge of lenalidomide was sufficient to describe the pharmacokinetics of lenalidomide in multiple tissues. Methods: A physiologically based pharmacokinetic model was developed using the Open Systems Pharmacology Suite to explore the pharmacokinetics of lenalidomide in a variety of tissues. Data were available for mice dosed intravenously at 0.5, 1.5, 5, and 10 mg/kg, with concentrations measured in plasma, brain, heart, kidney, liver, lung, muscle, and spleen. P-gp expression and activity were sourced from the literature. Results: The model predictions in plasma, liver, and lung were representative of the observed data (median prediction error 13%, − 10%, and 30%, respectively, with 90% confidence intervals including zero), while other tissue predictions showed sufficient similarity to the observed data. Contrary to the data, model predictions for the brain showed no drug reaching brain tissue when P-gp was expressed at the blood-brain barrier. The data were better described by basolateral transporters at the intracellular wall. Local sensitivity analysis showed that transporter activity was the most sensitive parameter in these models for exposure. Conclusion: As P-gp transport at the blood-brain barrier did not explain the observed brain concentrations alone, there may be other transporters involved in lenalidomide disposition. Refereed/Peer-reviewed

Introduction To investigate the pharmacokinetics and pharmacodynamics of oral pimobendan in conscious, healthy cats. Animals Eight healthy adult cats. Materials and methods A randomised, single-blinded, crossover design was used. Two oral doses of pimobendan (0.625-mg [LD], 1.25-mg [HD]) and a control substance (3-mL water) were administered to each cat. Blood collection, echocardiography, and oscillometric blood pressure measurements were performed repeatedly for 12 h following each dose. Plasma concentrations of pimobendan and the active metabolite, O-desmethylpimobendan (ODMP), were quantified using ultra-high-performance liquid chromatography tandem mass spectrometry. Cardiovascular parameters were evaluated for between- and within-treatment effects over time using linear mixed modelling. Results Pimobendan was rapidly absorbed and converted to ODMP with the pimobendan AUC0–∞ greater than ODMP AUC0–∞ (ODMP:pimobendan AUC0–∞ ratio 0.6 [LD] and 0.5 [HD]) despite a longer elimination half-life of ODMP (pimobendan t1/2 0.8 h vs. ODMP t1/2 1.6 h [LD]; pimobendan t1/2 0.7 h vs. ODMP t1/2 1.3 h [HD]). Averaged across all time points, pimobendan increased several measures of systolic function; however, its effect could not be further characterised. Although treatment was well-tolerated, two cats vomited following HD and another had a ventricular premature beat recorded following LD. Conclusions The lower ODMP:pimobendan AUC0–∞ ratio compared to that observed previously in dogs suggests reduced metabolism in cats. Treatment effects were observed in measures of systolic function; however, the duration of action and differences in effects between the two pimobendan doses could not be characterised. Further studies are required to evaluate pimobendan in feline cardiovascular medicine.

Morphine is commonly used for pain management in preterm neonates. The aims of this study were to compare published models of neonatal pharmacokinetics of morphine and its metabolites with a new dataset, and to combine the characteristics of the best predictive models to design a meta-model for morphine and its metabolites in preterm neonates. Moreover, the concentration-analgesia relationship for morphine in this clinical setting was also investigated. A population of 30 preterm neonates (gestational age: 23-32weeks) received a loading dose of morphine (50-100μg/kg), followed by a continuous infusion (5-10μg/kg/h) until analgesia was no longer required. Pain was assessed using the Premature Infant Pain Profile. Five published population models were compared using numerical and graphical tests of goodness-of-fit and predictive performance. Population modelling was conducted using NONMEM® and the $PRIOR subroutine to describe the time-course of plasma concentrations of morphine, morphine-3-glucuronide, and morphine-6-glucuronide, and the concentration-analgesia relationship for morphine. No published model adequately described morphine concentrations in this new dataset. Previously published population pharmacokinetic models of morphine, morphine-3-glucuronide, and morphine-6-glucuronide were combined into a meta-model. The meta-model provided an adequate description of the time-course of morphine and the concentrations of its metabolites in preterm neonates. Allometric weight scaling was applied to all clearance and volume terms. Maturation of morphine clearance was described as a function of postmenstrual age, while maturation of metabolite elimination was described as a function of postnatal age. A clear relationship between morphine concentrations and pain score was not established.

AIMS Composite indices for quantifying rheumatoid arthritis (RA) disease activity such as the 28-joint disease activity score (DAS28) are comprised of single parameters ('metrics') in various combinations. Population modelling methods were used to evaluate single metrics for their ability to reflect changes in disease activity with a view to understanding and improving composite indices. METHODS A total of 11 single metrics of RA disease activity (tender and swollen joint counts, acute phase reactants and global health, pain and physical function assessments) were obtained from 203 patients with recent onset RA. Participants received combination disease-modifying anti-rheumatic drugs (DMARDs) according to a treat-to-target approach with a pre-defined protocol for treatment intensification. Models describing each metric's magnitude and variability of change from baseline to a single 'treated' state in the population were developed using nonmem(®) . Measures that displayed uniformly large changes between states across the population were ranked higher in terms of discriminatory capacity. RESULTS Joint counts demonstrated a greater ability to discriminate changes in RA disease activity than others. Correlations between metrics demonstrated that erythrocyte sedimentation rate (ESR) had limited relationships with others for baseline scores and changes in RA disease activity (r generally

Objectives The selection of sample times for a pharmacokinetic study is important when trapezoidal integration (e.g. non-compartmental analysis) is used to determine the area under the concentration–time curve (AUC). The aim of this study was to develop an algorithm that determines optimal times that provide the most accurate AUC by minimising trapezoidal integration error. Methods The algorithm required initial single individual or mean pooled concentration data but did not specifically require a prior pharmacokinetic model. Optimal sample intervals were determined by minimising trapezoidal error using a genetic algorithm followed by a quasi-Newton method. The method was evaluated against simulated and clinical datasets to determine the method's ability to estimate the AUC. Key findings The sample times produced by the algorithm were able to accurately estimate the AUC of pharmacokinetic profiles, with the relative AUC having 90% confidence intervals of 0.919–1.05 for profiles with two-compartment kinetics. When comparing the algorithm with rich sampling (e.g. phase I trial), the algorithm provided equivalent or superior sample times with fewer observations. Conclusions The creation of the algorithm and its companion web application allows users with limited pharmacometric or programming training can obtain optimal sampling times for pharmacokinetic studies.

In February 2019, the Population Approach Group of Australia and New Zealand (PAGANZ) celebrated 20 years of providing education, mentoring, and support to the Australasian pharmacometrics community. This article provides insight into the history of this group and its operations to the present date. Refereed/Peer-reviewed

Background and Objective: SUBA-itraconazole and Sporanox are two oral formulations of itraconazole. Drug–drug interactions with omeprazole have been previously reported; however, mechanistic understanding of the pharmacological and physiological interactions of omeprazole with orally administered itraconazole within a population modeling paradigm is lacking. The objective of this analysis was to mechanistically describe and quantify the effect of omeprazole on the pharmacokinetics of itraconazole and its major metabolite, hydroxyitraconazole from the SUBA itraconazole and Sporanox formulations. Methods: An in vitro–in vivo (IVIV) pharmacokinetic model of itraconazole and hydroxyitraconazole was developed including data from an omeprazole interaction study with SUBA itraconazole. Meta-models of gastric pH for healthy subjects and subjects receiving omeprazole were integrated into the IVIV model to capture omeprazole-mediated gastric pH changes on itraconazole dissolution and absorption. Results: Omeprazole influenced the kinetics of itraconazole through altering the dissolution and absorption due to the pH-dependent solubility of itraconazole, inhibition of efflux transporters, and inhibiting the metabolism of itraconazole and hydroxyitraconazole. The model-predicted population effects of omeprazole on itraconazole from SUBA-itraconazole were to increase the area under the concentration–time curve (AUC0–24) and maximum concentration (Cmax) by 35 and 31%, respectively, and to decrease AUC0–24 and Cmax from Sporanox by 68 and 76%, respectively. Conclusion: Unlike SUBA itraconazole, which requires basic pH for itraconazole release, the omeprazole-induced pH-mediated reduction in Sporanox dissolution overrides any increased exposure from the drug–drug interaction at hepatic metabolizing enzymes or efflux transporters. The model presented here is the most complete quantitative description of the pharmacokinetics of itraconazole and hydroxyitraconazole currently available. Refereed/Peer-reviewed

Aims: Lenalidomide is an immunomodulatory imide drug used broadly in the treatment of multiple myeloma and lymphoma. It continues to be evaluated in chroniclymphocytic leukaemia (CLL) at lower doses due to dose‐related toxicities including tumour flare and tumour lysis syndrome. This study aimed to develop a population pharmacokinetic model for lenalidomide in multiple cancers, including CLL, to identify any disease‐related differences in disposition. Methods: Lenalidomide concentrations from 4 clinical trials were collated (1999 samples, 125 subjects), covering 4 cancers (multiple myeloma, CLL, acute myeloid leukaemia and acute lymphoblastic leukaemia) and a large dose range (2.5–75 mg).A population pharmacokinetic model was developed with NONMEM and patient demographics were tested as covariates. Results: The data were best fitted by a 1‐compartment kinetic model with absorption described by 7 transit compartments. Clearance and volume of distribution were allometrically scaled for fat‐free mass. The population parameter estimates for apparent clearance, apparent volume of distribution and transit rate constant were 12 L/h(10.8–13.6), 68.8 L (61.8–76.3), and 13.5 h−1 (11.9–36.8) respectively. Patients with impaired renal function (creatinine clearance

Aim Leflunomide, via its active metabolite teriflunomide, is used in rheumatoid arthritis (RA) treatment, yet approximately 20 to 40% of patients cease due to toxicity. The aim was to develop a time-to-event model describing leflunomide cessation due to toxicity within a clinical cohort and to investigate potential predictors of cessation such as total and free teriflunomide exposure and pharmacogenetic influences. Methods This study included individuals enrolled in the Early Arthritis inception cohort at the Royal Adelaide Hospital between 2000 and 2013 who received leflunomide. A time-to-event model in nonmem was used to describe the time until leflunomide cessation and the influence of teriflunomide exposure and pharmacogenetic variants. Random censoring of individuals was simultaneously described. The clinical relevance of significant covariates was visualized via simulation. Results Data from 105 patients were analyzed, with 34 ceasing due to toxicity. The baseline dropout hazard and baseline random censoring hazard were best described by step functions changing over discrete time intervals. No statistically significant associations with teriflunomide exposure metrics were identified. Of the screened covariates, carriers of the C allele of CYP1A2 rs762551 had a 2.29 fold increase in cessation hazard compared with non-carriers (95% CI 2.24, 2.34, P = 0.016). Conclusions A time-to-event model described the time between leflunomide initiation and cessation due to side effects. The C allele of CYP1A2 rs762551 was linked to increased leflunomide toxicity, while no association with teriflunomide exposure was identified. Future research should continue to investigate exposure–toxicity relationships, as well as potentially toxic metabolites.

Itraconazole is an orally active antifungal agent that has complex and highly variable absorption kinetics that is highly affected by food. This study aimed to develop a population pharmacokinetic model for itraconazole and the active metabolite hydroxyitraconazole, in particular, quantifying the effects of food and formulation on oral absorption. Plasma pharmacokinetic data were collected from seven phase I crossover trials comparing the SUBA-itraconazole and Sporanox formulations of itraconazole. First, a model of single-dose itraconazole data was developed, which was then extended to the multidose data. Covariate effects on itraconazole were then examined before extending the model to describe hydroxyitraconazole. The final itraconazole model was a 2-compartment model with oral absorption described by 4-transit compartments. Multidose kinetics was described by total effective daily dose- and time-dependent changes in clearance and bioavailability. Hydroxyitraconazole was best described by a 1-compartment model with mixed first-order and Michaelis-Menten elimination for the single-dose data and a time-dependent clearance for the multidose data. The relative bioavailability of SUBA-itraconazole compared to that of Sporanox was 173% and was 21% less variable between subjects. Food resulted in a 27% reduction in bioavailability and 58% reduction in the transit absorption rate constant compared to that with the fasted state, irrespective of the formulation. This analysis presents the most extensive population pharmacokinetic model of itraconazole and hydroxyitraconazole in the literature performed in healthy subjects. The presented model can be used for simulating food effects on itraconazole exposure and for performing prestudy power analysis and sample size estimation, which are important aspects of clinical trial design of bioequivalence studies.

Pimobendan is an inodilator used in the treatment of canine congestive heart failure (CHF). The aim of this study was to investigate the pharmacokinetics and cardiovascular effects of a nonaqueous oral solution of pimobendan using a single-dose, operator-blinded, parallel-dose study design. Eight healthy dogs were divided into two treatment groups consisting of water (negative control) and pimobendan solution. Plasma samples and noninvasive measures of cardiovascular function were obtained over a 24-h period following dosing. Pimobendan and its active metabolite were quantified using an ultra-high-performance liquid chromatography–mass spectrometer (UHPLC-MS) assay. The oral pimobendan solution was rapidly absorbed [time taken to reach maximum concentration (Tmax) 1.1 h] and readily converted to the active metabolite (metabolite Tmax 1.3 h). The elimination half-life was short for both pimobendan and its active metabolite (0.9 and 1.6 h, respectively). Maximal cardiovascular effects occurred at 2–4 h after a single oral dose, with measurable effects occurring primarily in echocardiographic indices of systolic function. Significant effects persisted for

Aims To develop a population model describing the disease activity (DAS28) time course in patients with early rheumatoid arthritis (RA) treated with triple disease-modifying anti-rheumatic drug (DMARD) therapy (methotrexate, sulfasalazine and hydroxychloroquine). Methods DAS28 was obtained in 263 patients with early RA from initiation of therapy until 60 weeks. Using nonmem®, base models (DAS28 vs. time) and covariate influences were investigated for the population. Results The best model was an exponential model of DAS28 vs. time that was additive to baseline DAS28, with covariance between parameters, and a combined residual error model. Age and patient smoking history were covariates significantly affecting response to therapy. Population estimates were baseline DAS28 (5.7), extent of change in DAS28 (−2.8) and the half-life of disease activity (6.2 weeks; time to steady disease state achieved within approximately 30 weeks). Older individuals exhibited more severe baseline DAS28, described by a power function centred around 57 years (baseline DAS28 for 40- and 70-year-old patients were 5.4 vs. 5.8, respectively) and current smokers took longer to achieve a steady disease state (approximately 50 weeks). There was considerable within-patient random variability in DAS28 over time (empirical 90% CI for DAS28 in a population typical patient at 60 weeks: 1.8, 4.2 with median value of 2.8). Conclusions This is the first report of a disease activity model for early RA treated with triple DMARD therapy. Smoking and age were identified as covariates.

Poor profiles of potential drug candidates, including pharmacokinetic properties, have been acknowledged as a significant hindrance to the development of modern therapeutics. Contemporary drug discovery and development would be incomplete without the aid of molecular modeling (in-silico) techniques, allowing the prediction of pharmacokinetic properties such as clearance, unbound fraction, volume of distribution and bioavailability. As with all models, in-silico approaches are subject to their interpretability, a trait that must be balanced with accuracy when considering the development of new methods. The best models will always require reliable data to inform them, presenting significant challenges, particularly when appropriate in-vitro or in-vivo data may be difficult or time-consuming to obtain. This article seeks to review some of the key in-silico techniques used to predict key pharmacokinetic properties and give commentary on the current and future directions of the field.

Infliximab is an anti-tumour necrosis factor alpha monoclonal antibody used to treat inflammatory diseases. Many patients fail during induction and others respond initially but relapse during maintenance therapy. Although anti-drug antibodies (ADA) are associated with some clinical failures, there is evidence that some failures may be due to subtherapeutic exposure. Adapting doses based on clinical outcomes and trough concentrations can improve response and reduce the proportion that develop ADA, but identification of appropriate doses in the presence of time-varying patient factors is complicated. Several adaptive dosing strategies (label recommendations versus therapeutic drug monitoring with an established stepwise algorithm or proportional dose adjustments or Bayesian population pharmacokinetic model-based dosing) were simulated on a virtual population (constructed with time-varying covariates and random effects on individual pharmacokinetic parameters) using R to assess their relative performance. Strategies were evaluated on their ability to maintain trough infliximab concentrations above an established target, 3 mg/L, during maintenance phase. Model-based dosing was superior in maintaining target trough concentrations, showing individuals in maintenance achieving concentrations above the target faster and a lower proportion of individuals who developed ADA. Model-based dosing results were consistent across a range of baseline covariate groups. This in silico assessment of adaptive dosing strategies demonstrated that, when challenged with dynamic covariate and random effect changes occurring in individual pharmacokinetic parameters, model-based approaches were superior to other strategies. Model-based dosing has not been tested clinically; however, the potential benefits of model-based dosing for infliximab suggest that it should be investigated to reduce subtherapeutic exposure. Refereed/Peer-reviewed

Leflunomide's active metabolite teriflunomide inhibits dihydro-oroate dehydrogenase, an enzyme essential to proliferation of T lymphocytes. As teriflunomide must reach the target site to have this effect, this study assessed the distribution of teriflunomide into T lymphocytes, as intracellular concentrations may be a superior response biomarker to plasma concentrations. CD3 MicroBeads (Miltenyi Biotec, Bergisch Gladbach, Germany) were used to extract CD3(+) T cells from the peripheral blood of patients with rheumatoid arthritis who were taking a stable dose of leflunomide. Unbound plasma and intra-CD3(+) T cell teriflunomide concentrations were quantified using liquid chromatography-mass spectrometry. Concentration (log transformed) and partition differences were assessed through paired Student t tests. Sixteen patients provided plasma steady-state teriflunomide samples, and eight provided a sample 6-12 weeks later. At time-point one, the geometric mean teriflunomide concentration (range) in CD3(+) T cells was 18.12mg/L (6.15-42.26mg/L) compared with 69.75 mg/L (32.89-263.1mg/L) unbound in plasma (P < 0.001). The mean partition coefficient (range) for unbound plasma teriflunomide into CD3(+) T cells was 0.295 (0.092-0.632), which was significantly different from unity (P < 0.001). The median (range) change in teriflunomide concentration between the two time points was 14% (210% to 40%) in unbound plasma and 229% (269 to 138%) for CD3(+) T cells. Because teriflunomide concentrations in CD3(+) T cells were lower and displayed a higher intraindividual variability than the unbound plasma concentrations, its applicability as a therapeutic drug-monitoring marker may be limited. Refereed/Peer-reviewed

This article reviews the major physiological and physicochemical principles of the effect of food and gastrointestinal (GI) pH on the absorption and bioavailability of oral drugs, and the various absorption models that are used to describe/predict oral drug absorption. The rate and extent of oral drug absorption is determined by a complex interaction between a drug's physicochemical properties, GI physiologic factors, and the nature of the formulation administered. GI pH is an important factor that can markedly affect oral drug absorption and bioavailability as it may have significant influence on drug dissolution & solubility, drug release, drug stability, and intestinal permeability. Different regions of the GI tract have different drug absorptive properties. Thus, the transit time in each GI region and its variability between subjects may contribute to the variability in the rate and/or extent of drug absorption. Food-drug interactions can result in delayed, decreased, increased, and sometimes un-altered drug absorption. Food effects on oral absorption can be achieved by direct and indirect mechanisms. Various models have been proposed to describe oral absorption ranging from empirical models to the more sophisticated “mechanism-based” models. Through understanding of the physicochemical and physiological rate-limiting factors affecting oral absorption, modellers can implement simplified population-based modelling approaches that are less complex than whole-body physiologically-based models but still capture the essential elements in a physiological way and hence will be more suited for population modelling of large clinical data sets. It will also help formulation scientists to better predict formulation performance and to develop formulations that maximize oral bioavailability. Refereed/Peer-reviewed

Cats are the most common mammal to develop invasive fungal rhinosinusitis caused by cryptic species in Aspergillus section Fumigati that are resistant to azoles but susceptible to caspofungin. In this study nonlinear mixed-effects pharmacokinetic modeling and simulation was used to investigate caspofungin pharmacokinetics and explore dosing regimens in cats using caspofungin minimum effective concentrations (MECs). Plasma concentrations in healthy cats were determined using HPLC-MS/MS after administration of a single and seven consecutive daily intravenous doses of 1 mg/kg caspofungin. In the final pharmacokinetic model an optimum maximum concentration (Cmax): MEC ratio of 10-20 was used to guide caspofungin efficacy. Simulations were performed for dosing regimens (doses 0.25-2 mg/kg and 6-72 h dosing intervals) with and without inclusion of a loading dose. Using a 1 mg/kg dose Cmax first dose was 14.8 μg/mL, Cmax at steady state was 19.8 μg/mL, Cmin was 5 μg/mL and Cmax: MEC was >20 in 42.6% of cats after multiple doses. An optimal Cmax: MEC ratio was achieved in caspofungin simulations using 0.75 mg/kg q 24 h or 1 mg/kg q 72h. However, at 1 mg/kg q 72h, Cmin was < MEC ( 2.5 μg/mL for 98% of the population. Based on the modeling data this dosing regimen is likely to achieve target therapeutic concentrations, meet the proposed Cmax: MEC window and provide consistent exposure between doses. Refereed/Peer-reviewed

Non-compartmental analysis (NCA) is regarded as the standard for establishing bioequivalence, despite its limitations and the existence of alternative methods such as non-linear mixed effects modelling (NLMEM). Comparisons of NCA and NLMEM in bioequivalence testing have been limited to drugs with one-compartment kinetics and have included a large number of different approaches. A simulation tool was developed with the ability to rapidly compare NCA and NLMEM methods in determining bioequivalence using both R and NONMEM and applied to a drug with two-compartment pharmacokinetics. Concentration–time profiles were simulated where relative bioavailability, random unexplained variability (RUV) at the lower limit of quantification (LLOQ) differed between simulations. NLMEM analyses employed either the M1 or M3 methods for dealing with values below the LLOQ. It was used to elucidate the impact of changes in (i) RUV at the LLOQ, (ii) the extent of censoring data below the LLOQ and (iii) the concentration sampling times. The simulations showed NLMEM having a consistent 20–40% higher accuracy and sensitivity in identifying bioequivalent studies when compared to NCA, while NCA was found to have a 1–10% higher specificity than NLMEM. Increasing data censoring by increasing the LLOQ resulted in decreases of ~10% to the accuracy and sensitivity of NCA, with minimal effects on NLMEM. The tool provides a platform for comparing NCA and NLMEM methods and its use can be extended beyond the scenarios reported here. In the situations examined it is seen that NLMEM is more accurate than NCA and may offer some advantages in the determination of bioequivalence. Refereed/Peer-reviewed

The literature presently lacks a population pharmacokinetic analysis of doxycycline. This study aimed to develop a population pharmacokinetic model of doxycycline plasma concentrations that could be used to assess the power of bioequivalence between Doryx delayed-release tablets and Doryx MPC. Doxycycline pharmacokinetic data were available from eight phase 1 clinical trials following single/multiple doses of conventional-release doxycycline capsules, Doryx delayed-release tablets, and Doryx MPC under fed and fasted conditions. A population pharmacokinetic model was developed in a stepwise manner using NONMEM, version 7.3. The final covariate model was developed according to a forward inclusion ( P < 0.01) and then backward deletion ( P < 0.001) procedure. The final model was a two-compartment model with two-transit absorption compartments. Structural covariates in the base model included formulation effects on relative bioavailability ( F ), absorption lag (ALAG), and the transit absorption rate (KTR) under the fed status. An absorption delay (lag) for the fed status (FTLAG2 = 0.203 h) was also included in the model as a structural covariate. The fed status was observed to decrease F by 10.5%, and the effect of female sex was a 14.4% increase in clearance. The manuscript presents the first population pharmacokinetic model of doxycycline plasma concentrations following oral doxycycline administration. The model was used to assess the power of bioequivalence between Doryx delayed-release tablets and Doryx MPC, and it could potentially be used to critically examine and optimize doxycycline dose regimens.

This study aimed to conduct a quantitative meta-analysis for the values of, and variability in, gastrointestinal (GI) pH in the different GI segments; characterize the effect of food on the values and variability in these parameters; and present quantitative meta-models of distributions of GI pH to help inform models of oral drug absorption. The literature was systemically reviewed for the values of, and the variability in, GI pH under fed and fasted conditions. The GI tract was categorized into the following 10 distinct regions: stomach (proximal, mid-distal), duodenum (proximal, mid-distal), jejunum and ileum (proximal, mid, and distal small intestine), and colon (ascending, transverse, and descending colon). Meta-analysis used the “metafor” package of the R language. The time course of postprandial stomach pH was modeled using NONMEM. Food significantly influenced the estimated meta-mean stomach and duodenal pH but had no significant influence on small intestinal and colonic pH. The time course of postprandial pH was described using an exponential model. Increased meal caloric content increased the extent and duration of postprandial gastric pH buffering. The different parts of the small intestine had significantly different pH. Colonic pH was significantly different for descending but not for ascending and transverse colon. Knowledge of GI pH is important for the formulation design of the pH-dependent dosage forms and in understanding the dissolution and absorption of orally administered drugs. The meta-models of GI pH may also be used as part of semi-physiological pharmacokinetic models to characterize the effect of GI pH on the in vivo drug release and pharmacokinetics. Refereed/Peer-reviewed

This paper describes the development of a sensitive high performance liquid chromatography (HPLC) method for quantification of rosiglitazone in sheep plasma and amniotic fluid. Samples were prepared by liquid-liquid extraction using tert-butyl methyl ether, and rosiglitazone was quantitated by HPLC using a C18 column and fluorescence detector with an excitation wavelength of 247 nm and emission wavelength of 367 nm. The mobile phase consisted of ammonium acetate (10 mM, pH 5.2) and acetonitrile (56.5:43.5, v/v) with a flow rate of 1 ml/min. Ketoconazole was used as the internal standard (IS). The plasma calibration curve was linear over the range of 2.5-250 ng/ml (mean r2 = 0.9940 ± 0.0024; n = 6) with accuracy of 99.4-102.8% over the calibration range. The intra-day and inter-day coefficient of variation (%CV, percent coefficient of variation) were in the range of 0.01-8.68% in sheep plasma. Similar performance was achieved for amniotic fluid. The described method was successfully applied to quantitate rosiglitazone concentrations in the pregnant ewe and her fetus. Refereed/Peer-reviewed

Pharmacokinetic/pharmacodynamic (PK/PD) modelling can be used to characterize the relationship between dose regimen of opioids, plasma concentration and effect of opioids, which in turn can lead to more rational treatment regimens of pain. The aim of this study was to investigate the concentration-effect relationship for transdermal buprenorphine and fentanyl in experimentally induced pain. Twenty-two healthy volunteers were randomized to receive transdermal patches with fentanyl (25 μg/hr, 72 hr), buprenorphine (20 μg/hr, 144 hr) or placebo. The experimental pain tests were pressure at the tibial bone, cutaneous thermal stimulation, cold pressor test (conditioning stimulus (3 ± 0.3°C cold water), nerve growth factor-induced muscle soreness and intradermal capsaicin-induced hyperalgesia and allodynia. Experiments were carried out at baseline, 24, 48, 72 and 144 hr after application of patches. Time-course of placebo was described first and was afterwards added to the description of the time-courses of buprenorphine and fentanyl. This was either described by zero (no drug effect), linear or E(max) model concentration-effect relationships. Time-dependent changes in pain measures in the placebo arm were described by linear or quadratic functions. The time-course of fentanyl and buprenorphine plasma concentrations was complex but could be represented by cubic spline interpolation in the models. Buprenorphine significantly attenuated bone-associated pain, heat pain, nerve growth factor-induced soreness and cold pressor pain. Fentanyl significantly attenuated cold pressor pain for the administered dose regimens. Although the PK/PD relationship for both drugs could be described with similar models, tissue-differentiated analgesic effects between buprenorphine and fentanyl was shown.

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT Both methadone and buprenorphine cross the placenta during pregnancy which can result in neonatal abstinence syndrome (NAS). Reduced NAS in buprenorphine-exposed infants may be as a result of relatively less buprenorphine reaching the foetal circulation compared with methadone. WHAT THIS STUDY ADDS Under chronic dosing conditions in humans, the transfer of methadone and buprenorphine appears greater than previously found in human placental tissue in vitro models. Infants born to buprenorphine maintained women may be exposed to relatively less of the maternal dose compared with infants born to women maintained on methadone during pregnancy. There is stereoselectivity in the transfer of the individual enantiomers of methadone across the placenta to the foetal circulation. AIMS The aim of this study was to compare the transfer of buprenorphine and methadone between maternal and cord blood in women under chronic dosing conditions and to determine if differences exist in the transfer of the two methadone enantiomers. METHODS Maternal and cord blood samples were collected at delivery from women maintained on methadone (35, 25–140 mg day−1) (median; range) or buprenorphine (6.00, 2–20 mg day−1) during pregnancy. Plasma concentration ratios are presented as an indicator of foetal exposure relative to the mother. RESULTS Methadone was quantified in all samples, with cord : maternal plasma methadone concentration ratios (n= 15 mother-infant pairs) being significantly higher (P < 0.0001; mean difference (MD) 0.07; 95% confidence interval (CI) 0.048, 0.092) for the active (R)-methadone enantiomer (0.41; 0.19, 0.56) (median; range) compared with (S)-methadone (0.36; 0.15, 0.53). (R)- : (S)-methadone concentration ratios were also significantly higher (P < 0.0001; MD 0.24 95% CI 0.300, 0.180) for cord (1.40; 0.95, 1.67) compared with maternal plasma (1.16; 0.81, 1.38). Half the infant buprenorphine samples were below the assay lower limit of quantification (LLOQ) (0.125 ng ml−1). The latter was four-fold lower than the LLOQ for methadone (0.50 ng ml−1). The cord : maternal plasma buprenorphine concentration ratio (n= 9 mother-infant pairs) was 0.35; 0.14, 0.47 and for norbuprenorphine 0.49; 0.24, 0.91. CONCLUSIONS The transfer of the individual methadone enantiomers to the foetal circulation is stereoselective. Infants born to buprenorphine maintained women are not exposed to a greater proportion of the maternal dose compared with methadone and may be exposed to relatively less of the maternal dose compared with infants born to women maintained on methadone during pregnancy.

The aim of this study was to compare the pharmacokinetic profile, as well as the efficacy and tolerability, of i.n. and i.v. administration of fentanyl in acute, episodic pain in patients undergoing third-molar extraction.In this randomized, double-blind, double-dummy, 2-way, crossover study, patients were randomized to receive 1 of 4 doses (75, 100, 150, or 200 microg) by both the i.n. and i.v. routes in random order, after each of 2 separate molar extractions (interval,or=1 week). Venous blood samples were obtained for quantification of plasma fentanyl concentrations before and at 1, 3, 5, 7, 9, 12, 15, 25, 40, 60, 90, 120, and 180 minutes after administration. Pain scores (on an 11-point numeric rating scale) were recorded before and at 15, 30, 45, 60, 75, 90, 105, 120, 150, 180, 210, and 240 minutes. Patients indicated the times at which they perceived meaningful pain relief (onset of action) and at which analgesia ended (duration of effect), after which they were able to use rescue medication (time to rescue medication use).A total of 24 patients were enrolled (in all, 47 extractions) (46% male; mean age, 24.1 years; 94% white, 6% Asian). Mean T(max) values were 12.8 and 6.0 minutes (P0.001), times to onset of analgesia were 7 and 2 minutes (P0.001), and durations of effect were 56 and 59 minutes after i.n. and i.v. administration (P=NS), respectively. Differences in the onsets and durations of analgesia after i.n. and i.v. administration of single doses were not significantly different, and neither was the difference in overall analgesia, with pain scores returning to near-predose values at statistically similar times after dosing. Duration of effect was directly related to i.n. fentanyl dose, with significantly less use of rescue medication after i.n. than after i.v. administration (P0.005). The i.n. and i.v. formulations were both well tolerated, with similar numbers of nasally related adverse events recorded for both routes of administration.Onsets and durations of analgesia were not significantly different between single doses of i.n. and i.v. fentanyl in these adults undergoing third-molar extraction. Both i.n. and i.v. administration were generally well tolerated.

BACKGROUND: Polymorphisms in the opioid receptor genes may affect the pharmacodynamics (PD) of oxycodone and be part of the reason behind the diversity in clinical response. The aim of the analysis was to model the exposure-response profile of oxycodone for three different pain variables and search for genetic covariates. Model simulations were used to predict how population and effect-size impact the power to detect clinical significant SNPs.METHOD: The population pharmacokinetic-pharmacodynamic (PKPD) model of oral single-dosed oxycodone was based on pooled data from three published studies in healthy volunteers. Pain tolerance data from muscle pressure (n=36), visceral pressure (n=54) and skin pinch (n=34) were included. Genetic associations with 18 opioid-receptor SNPs were explored using a stepwise covariate approach. Model simulations were performed using the estimated model parameters.RESULTS: None of the selected SNPs were associated with analgesic response of oxycodone at P50% for similar population sizes) or large populations (n>200 for a 20% response difference) are necessary to identify clinical significant SNP effects due to high population variability.CONCLUSION: A population PKPD model has been developed for oral oxycodone using three different pain variables to explore impact of genetic covariates and study design. None of the selected polymorphisms were significantly associated with analgesic response of oxycodone, but an association of baseline response in muscle cuff pressure with two OPRK1 SNPs was identified.

Objective To investigate the influence of CYP2D6 genotype on the oral clearance of (R)-, (S)- and rac-methadone. Methods In this retrospective study, CYP2D6 genotypes were identified in 56 methadone maintained subjects. Plasma concentrations of (R)-, (S)- and rac-methadone were determined by stereoselective HPLC and sufficient data were available to estimate the apparent oral clearances of (R)-, (S)- and rac-methadone using a population kinetic model in 37 of the genotyped subjects. Results The CYP2D6 allele frequencies were similar to those previously reported in Caucasians, the most common being: CYP2D6*1 (35.2%), CYP2D6*2 (12.0%) and CYP2D6*4 (22.2%). Three unknown SNPs were found in four subjects: 1811G > A (n = 1), 1834C > T (n = 1) and 2720G > C (n = 2). The oral clearances of (R)-, (S)- and rac-methadone varied 5.4-, 6.8- and 6.1-fold, respectively. No significant differences in methadone oral clearance were found between CYP2D6 genotypic PM, IM and EM (p = 0.57, 0.40 and 0.43 for (R)-, (S)- and rac-methadone, respectively). Only 1 subject had duplication of functional CYP2D6 alleles and the oral clearance of the three analytes was not markedly altered. Conclusions CYP2D6 poor, intermediate and extensive metabolizer genotypes did not appear to impact on the oral clearance of (R)-, (S)- or rac-methadone. In addition, methadone dosage requirements were not influenced by CYP2D6 genotypes in these subjects. However, the impact of duplication of functional CYP2D6 alleles on oral clearance and dosage requirements requires further investigation.

Background and purpose: At present there are few data regarding the rate and extent of brain–blood partitioning of the opioid active metabolite of morphine, morphine-6-glucuronide (M6G). In this study the cerebral kinetics of M6G were determined, after a short-term intravenous infusion, in chronically instrumented conscious sheep. Experimental approach: Five sheep received an intravenous infusion of M6G 2.2 mg kg-1 over a four-minute period. Non-linear mixed-effects analysis, with hybrid physiologically based kinetic models, was used to estimate cerebral kinetics from the arterio-sagittal sinus concentration gradients and cerebral blood flow measurements. Key results: A membrane limited model was selected as the final model. The blood-brain equilibration of M6G was relatively slow (time to reach 50% equilibration of the deep compartment 5.8 min), with low membrane permeability (PS, population mean, 2.5 ml min-1) from the initial compartment (V1, 13.7 ml) to a small deep distribution volume (V2) of 18.4 ml. There was some between-animal variability (%CV) in the initial distribution volume (29%), but this was not identified for PS or V2. Conclusion and Implications: Pharmacokinetic modelling of M6G showed a delayed equilibration between brain and blood of a nature that is primarily limited by permeability across the blood-brain-barrier, in accordance with its physico-chemical properties. British Journal of Pharmacology (2006) 149, 754–760. doi:10.1038/sj.bjp.0706916

Summary Background The sources of variability in the blood concentrations of morphine after subcutaneous administration are unclear. Methods A population pharmacokinetic model (using NONMEM) was used to examine data on morphine disposition following its administration (5 mg) via an indwelling s.c. cannula to elderly post-operative patients. Results A two-compartment model with first order absorption was preferred. The between subject variability (BSV) was between 45 and 95% for the systemic parameters, and 71% for the absorption rate constant. The rate of subcutaneous absorption of morphine was not inherently more variable than other factors affecting morphine disposition. A sensitivity analysis supported these findings. In eight patients who received a second dose of morphine approximately 300 min after the first dose, the between occasion variability (BOV) of the parameters was up to 81% of the total variability. Conclusions A significant component of the variability in morphine disposition was due to factors that could alter within a patient between doses. Taken together, this reinforces the need to titrate the dose of morphine to clinical endpoints, and questions the concept of a fixed dose for a given patient over time, especially if given at fixed dose intervals (i.e. on a time-contingent basis).

The delayed onset and long duration of action of the opioid agonist levo-α-acetyl-methadol (LAAM) has been attributed to the formation of active metabolites. However, at present, little is known about the time course of blood–brain equilibration of LAAM itself. The cerebral kinetics of LAAM were quantified using physiologically based kinetic models and a conscious chronically instrumented sheep preparation. Seven sheep were administered 4 min intravenous infusions of 30 mg LAAM. Concentrations of LAAM and N-demethylated metabolites (nor-LAAM and di-nor-LAAM) in whole blood (0–75 min) were measured using a validated HPLC assay. LAAM did not alter cerebral blood flow, mean arterial pressure or cause significant respiratory depression. Cardiac output was similar to baseline at 4 min, but decreased by 30% at 10 min and remained at this level for the duration of the 75 min study period. Cerebral kinetics were best described by a membrane-limited model, with a relatively slow blood–brain tissue equilibration half-life of 22 min due to intermediate permeability (56 ml min−1) and a large cerebral distribution volume (724 ml). In conclusion, pharmacokinetic–pharmacodynamic modelling of LAAM should account for the large equilibration delay between brain and blood caused by slow equilibration kinetics. This may account for some of the delay in onset of effect previously attributed to the delayed appearance of active metabolites in blood. British Journal of Pharmacology (2006) 147, 209–217. doi:10.1038/sj.bjp.0706470

To construct a population pharmacokinetic model for methadone enantiomers in the setting of methadone maintenance treatment for opioid dependence.A population pharmacokinetic model was developed using P-Pharm software for rac-, (R)- and (S)-methadone using data (8-13 plasma samples per subject) obtained from 59 methadone maintenance patients during one interdosing interval at steady state. The patients were randomly assigned to either a development (n = 38) or a validation dataset (n = 21). The model was refined by inclusion of all subjects to construct a final basic model, which was used to construct a covariate model.A population-based two-compartment open model with first-order absorption and lag time was developed and validated for all analytes. The population geometric mean (coefficient of variation) of maximum a posteriori probability Bayesian estimated values for clearance, terminal half-life and volume of distribution at steady-state of the active (R)-enantiomer were 8.7 (42%) l h(-1), 51 (45%) h and 597 (45%) l, respectively. For all analytes, the volume of the central compartment was decreased with increasing plasma alpha(1)-acid glycoprotein concentration and was lower in females, while the delay in absorption was longer at higher doses. No covariates were identified for apparent oral clearance. The apparent oral clearance of (R)-methadone (geometric mean ratio; 95% confidence interval) was 105% (99, 110), that of (S)-methadone (P = 0.19), while (R)-methadone V(c)/F (154%; 151, 157), V(dss) /F (173%; 164, 183), t(1/2beta) (162%; 153, 172) and mean residence time (166%; 156, 176) were significantly greater (P0.0001) than for (S)-methadone. The population pharmacokinetic models were able to predict accurately oral clearance values from limited (one or two samples) blood sampling protocols.The substantial stereoselectivity in methadone disposition reinforces the potential for misinterpretation of racemic methadone disposition data. The marked interindividual variability in (R)-methadone clearance, with no covariates identified, highlights the need for alternative methods to determine an individual's metabolic clearance. The ability to predict (R)-methadone clearance from one to two blood samples at steady state may prove clinically useful if a drug-drug interaction or poor adherence are suspected and guide the prescriber in deciding if a client's request for a dose increase is warranted or whether an alternative opioid would be more appropriate.

Aim: To determine the Michaelis-Menten kinetics of hydrocodone metabolism to its O- and N-demethylated products, hydromorphone and norhydrocodone, to determine the individual cytochrome p450 enzymes involved, and to predict the in vivo hepatic intrinsic clearance of hydrocodone via these pathways. Methods: Liver microsomes from six CYP2D6 extensive metabolizers (EM) and one CYP2D6 poor metabolizer (PM) were used to determine the kinetics of hydromorphone and norhydrocodone formation. Chemical and antibody inhibitors were used to identify the cytochrome p450 isoforms catalyzing these pathways. Expressed recombinant cytochrome p450 enzymes were used to characterize further the metabolism of hydrocodone. Results: Hydromorphone formation in liver microsomes from CYP2D6 EMs was dependent on a high affinity enzyme (Km = 26 microm) contributing 95%, and to a lesser degree a low affinity enzyme (Km = 3.4 mm). In contrast, only a low affinity enzyme (Km = 8.5 mm) formed this metabolite in the liver from the CYP2D6 PM, with significantly decreased hydromorphone formation compared with the livers from the EMs. Norhydrocodone was formed by a single low affinity enzyme (Km = 5.1 mm) in livers from both CYP2D6 EM and PM. Recombinant CYP2D6 and CYP3A4 formed only hydromorphone and only norhydrocodone, respectively. Hydromorphone formation was inhibited by quinidine (a selective inhibitor of CYP2D6 activity), and monoclonal antibodies specific to CYP2D6. Troleandomycin, ketoconazole (both CYP3A4 inhibitors) and monoclonal antibodies specific for CYP3A4 inhibited norhydrocodone formation. Extrapolation of in vitro to in vivo data resulted in a predicted total hepatic clearance of 227 ml x h-1 x kg-1 and 124 ml x h-1 x kg-1 for CYP2D6 EM and PM, respectively. Conclusions: The O-demethylation of hydrocodone is predominantly catalyzed by CYP2D6 and to a lesser extent by an unknown low affinity cytochrome p450 enzyme. Norhydrocodone formation was attributed to CYP3A4. Comparison of recalculated published clearance data for hydrocodone, with those predicted in the present work, indicate that about 40% of the clearance of hydrocodone is via non-CYP pathways. Our data also suggest that the genetic polymorphisms of CYP2D6 may influence hydrocodone metabolism and its therapeutic efficacy.

The aim of this study was to develop population pharmacokinetic-pharmacodynamic models for morphine in experimental pain induced by skin heat and muscle pressure, and to evaluate the experimental pain models with regard to assessment of morphine pharmacodynamics. In a randomised, double-blind, placebo-controlled, crossover study, 39 healthy volunteers received an oral dose of 30 mg morphine hydrochloride or placebo. Non-linear mixed effects modelling was used to describe the plasma concentrations of morphine and metabolites, and the analgesic effect of morphine on experimental pain in skin and muscle. Baseline pain metrics varied between individuals and occasions, and were described with interindividual and interoccasion variability. Placebo-response did not change with time. For both pain metrics, morphine effect was proportional to baseline pain and was described with a linear model with interindividual variability on drug effect slope and linked to an effect compartment for muscle pressure. The models indicate that a steady-state morphine concentration of 21 ng/ml causes 33% and 0.84% increases in stimulus intensity from baseline for muscle pressure and skin heat, respectively. The population pharmacokinetic-pharmacodynamic models developed in this study indicate that mechanical stimulation of muscle is a more clinically relevant pain stimulus for the assessment of morphine pharmacodynamics than thermal stimulation of skin. Refereed/Peer-reviewed

Introduction The analytical solutions to compartmental pharmacokinetic models are well known, but have not been presented in a form that easily allows for complex dosing regimen and changes in covariate/parameter values that may occur at discrete times within and/or between dosing intervals. Methods Laplace transforms were used to derive ADVAN-style analytical solutions for 1, 2, and 3 compartment pharmacokinetic linear models of intravenous and first-order absorption drug administration. The equations calculate the change in drug amounts in each compartment of the model over a time interval (t; t = t2 − t1) accounting for any dose or covariate events acting in the time interval. The equations were coded in the R language and used to simulate the time-course of drug amounts in each compartment of the systems. The equations were validated against commercial software [NONMEM (Beal, Sheiner, Boeckmann, & Bauer, 2009)] output to assess their capability to handle both complex dosage regimens and the effect of changes in covariate/parameter values that may occur at discrete times within or between dosing intervals. Results For all tested pharmacokinetic models, the time-course of drug amounts using the ADVAN-style analytical solutions were identical to NONMEM outputs to at least four significant figures, confirming the validity of the presented equations. Discussion To our knowledge, this paper presents the ADVAN-style equations for common pharmacokinetic models in the literature for the first time. The presented ADVAN-style equations overcome obstacles to implementing the classical analytical solutions in software, and have speed advantages over solutions using differential equation solvers. The equations presented in this paper fill a gap in the pharmacokinetic literature, and it is expected that these equations will facilitate the investigation of useful open-source software for modelling pharmacokinetic data.

Using a modelling approach, this study aimed to (i) examine whether the pharmacodynamics of the analgesic and antihyperalgesic effects of morphine differ; (ii) investigate the influence of demographic, pain sensitivity and genetic (OPRM1) variables on between-subject variability of morphine pharmacokinetics and pharmacodynamics in human experimental pain models. The study was a randomized, double-blind, 5-arm, cross-over, placebo-controlled study. The psychophysical cutaneous pain tests, electrical pain tolerance (EPTo) and secondary hyperalgesia areas (2HA) were studied in 28 healthy individuals (15 males). The subjects were chosen based on a previous trial where 100 subjects rated (VAS) their pain during a heat injury (47°C, 7 min., 12.5 cm2). The 33% lowest- and highest pain-sensitive subjects were offered participation in the present study. A two-compartment linear model with allometric scaling for weight provided the best description of the plasma concentration-time profile of morphine. Changes in the EPTo and 2HA responses with time during the placebo treatment were best described by a linear model and a quadratic model, respectively. The model discrimination process showed clear evidence for adding between-occasion variability (BOV) on baseline and the placebo slope for EPTo and 2HA, respectively. The sensitivity covariate was significant on baseline EPTo values and genetics as a covariate on the placebo slope for 2HA. The analgesic and antihyperalgesic effects of morphine were pharmacologically distinct as the models had different effect site equilibration half-lives and different covariate effects. Morphine had negligible effect on 2HA, but significant effect on EPTo. Refereed/Peer-reviewed

Aims To investigate the steady-state pharmacokinetics of (R)- and (S)-methadone in a methadone maintenance population. Methods Eighteen patients recruited from a public methadone maintenance program underwent an interdosing interval pharmacokinetic study. Plasma and urine samples were collected and analysed for methadone and its major metabolite (EDDP) using stereoselective h.p.l.c. Methadone plasma protein binding was examined using ultrafiltration, and plasma α 1 -acid glycoprotein concentrations were quantified by radial immunoassay. Results (R)-methadone had a significantly (P< 0.05) greater unbound fraction (mean 173%) and total renal clearance (182%) compared with (S)-methadone, while maximum measured plasma concentrations (83%) and apparent partial clearance of methadone to EDDP (76%) were significantly (P< 0.001) lower. When protein binding was considered (R)-methadone plasma clearance of the unbound fraction (59%) and apparent partial intrinsic clearance to EDDP (44%) were significantly (P

A stereoselective HPLC assay was developed for the quantification of the enantiomers of methadone and its major oxidative metabolite, 2-ethylidene-1,5-dimethyl-3,3-diphenylpyrrolidine (EDDP) in human urine. The compounds were quantified in a single assay following liquid-liquid extraction and stereoselective HPLC with UV detection. Calibration curve concentrations ranged from 0.125 to 12.5 microM for each enantiomer. Assay performance was assessed using quality control samples, and the inter- and intra-assay bias (10%) and precision (15%) were acceptable for all compounds. The assay was successfully used to quantitate the enantiomers of methadone and EDDP in urine samples obtained from subjects receiving methadone maintenance therapy.

Aims To investigate the kinetics of CYP-mediated N-demethylation of methadone in human liver microsomes, and examine the role of stereoselectivity and CYP isoforms involved.

The purpose of the study was to investigate placebo and buprenorphine effects on event-related potentials (ERPs) in experimental pain and the potential benefit of population pharmacodynamic modelling in data analysis. Nineteen healthy volunteers received transdermal placebo and buprenorphine in a cross-over study. Drug plasma concentrations and ERPs after electrical stimulation at the median nerve with intensity adjusted to pain detection threshold were recorded until 144hrs after administration. Placebo and concentration-effect models were fitted to data using non-linear mixed-effects modelling implemented in NONMEM (V7.2.0.). Pharmacodynamic models were developed to adequately describe both placebo and buprenorphine ERP data. Models predicted significant placebo effects, but did not predict significant effects related to buprenorphine concentration. Models revealed that ERPs varied both between subjects and between study occasions. ERPs were found to be reproducible within subjects and occasions as population variance was found to be eight times higher than the unexplained variances. Between-subject variance accounted for more than 75% of the population variance. In conclusion, pharmacodynamic modelling was successfully implemented to allow for placebo and variability correction in ERP of experimental pain. Improved outcome of ERP studies can be expected if variation between subjects and study occasions can be identified and described. Refereed/Peer-reviewed

Summary What is known and objective A common polymorphism (C1858T) in the gene that encodes the protein tyrosine phosphatase non-receptor type 22 (PTPN22) is associated with altered T-cell responses and increased susceptibility to rheumatoid arthritis (RA) and other autoimmune diseases. Teriflunomide, the active metabolite of leflunomide, reduces T-cell responses through inhibition of tyrosine kinase p56LCK. We examined a potential association between PTPN22 genotype and response or toxicity to leflunomide in Caucasian RA patients taking leflunomide in combination with other disease-modifying antirheumatic drugs (DMARDs). Methods Patients enrolled in the Royal Adelaide Hospital RA inception cohort and taking leflunomide were eligible for inclusion. Participants were followed for 12 months after leflunomide initiation or until either another DMARD was added or leflunomide was ceased. Clinical response according to change in 28-joint Disease Activity Score (DAS28) and cessation due to toxicity were assessed. Results and discussion A total of 94 participants were included in the study, 75 of whom carried the CC genotype, 18 the CT, whereas one individual carried the TT genotype. Over the first 12 months of leflunomide treatment, there was no statistically significant relationship between carrying the T allele and change in DAS28 (−0·84 vs. −1·15, P = 0·446) nor with cessation of leflunomide treatment due to side effects (P = 0·433). These results indicate that PTPN22 C1858T genotype has no effect on response or toxicity outcomes in leflunomide-treated RA patients. What is new and conclusion This is the first study to evaluate the biologically plausible hypothesis that PTPN22 genotype might be a predictor of response/toxicity to leflunomide therapy. Despite this, PTPN22 genotype was not associated with leflunomide response or toxicity in patients with RA.

Irinotecan eluting embolization beads (DEBIRI) are currently being evaluated in the clinic for the treatment of colorectal cancer metastases to the liver. The aim of this study was to determine the safety and pharmacokinetics associated with two cycles of hepatic embolization using DEBIRI followed by intravenous administration of irinotecan. Pigs were embolized with DEBIRI (100-300 μm, 100 mg dose, n = 6) and blood samples taken over 24 h to determine plasma levels of irinotecan and SN-38 metabolite and for haematology and biochemistry. At 24 h an IV infusion of 250 mg/m(2) of irinotecan was administered and the plasma levels taken again. This cycle was repeated 3 weeks later. A single animal was subjected to a more aggressive regimen of embolization with 200 mg bead dose and IV of 350 mg/m(2) for two cycles. Three animals were sacrificed at 6 weeks and the remaining four (n = 3 standard dose, n = 1 high dose) animals at 12 weeks and detailed histopathology performed. All animals tolerated the treatments well, with only minor changes in haematological and biochemical parameters. There was no overlap in drug plasma levels observed from the bead and IV treatments when given 24 h apart and no difference between the pharmacokinetic profiles of the two cycles separated by 3 weeks. Irinotecan plasma AUC values were similar in both the embolization and IV arms of the study. C(max) values obtained during the IV arms of the study are approximately double that of the embolization arms whilst T(max) times are shorter in the IV arms, supporting extended release of drug from the beads. Bioavailability for bead-based delivery was double that for IV administration, which was attributed to reduced clearance of the drug when delivered by this route. No additive toxicity was observed as a consequence of the combined treatments. The combination of irinotecan delivery via drug eluting bead and IV was well-tolerated with no significant clinical effects. Pharmacokinetic analyses suggest the bioavailability from bead-based delivery of drug is double that of IV infusion, attributable to reduced drug clearance for the former.

Pharmacokinetic (PK) and pharmacodynamic (PD) data were available from a study of a nasal delivery system for the opioid analgesic fentanyl, together with data on the kinetics of fentanyl in arterial blood in man, and in the lung and brain of sheep. Our aim was to reconcile these data using a physiologically-based population recirculatory PK-PD model, with emphasis on achieving a meta-model that could simultaneously account for the arterial and venous (arm) concentrations of fentanyl, could relate PD effects (pain scores) to the CNS concentrations of fentanyl, and could account for the effect of body size and age on fentanyl kinetics. Data on the concentration gradients of fentanyl across brain, lung and muscle were used to develop sub-models of fentanyl kinetics in these organs. The sub-models were incorporated into a "whole body" recirculatory model by adding additional sub-models for a venous mixing compartment, the liver and gut, the kidney and the "rest of the body" with blood flows and organ volumes based on values for a Standard Man. Inter-individual variability was achieved by allometric scaling of organ size and blood flows, evidence-based assumptions about the effect of weight and age on cardiac output, and inter-individual variability in the free fraction in plasma and hepatic extraction of fentanyl. Post-operative pain scores were found to be temporally related to the predicted brain concentrations of fentanyl.We conclude that a physiologically-based metamodelling approach was able to describe clinical PK-PD studies of fentanyl whilst providing a mechanistic interpretation of key aspects of its disposition Refereed/Peer-reviewed

The aim of this study was to compare the transfer of buprenorphine and methadone between maternal and cord blood in women under chronic dosing conditions and to determine if differences exist in the transfer of the two methadone enantiomers.Maternal and cord blood samples were collected at delivery from women maintained on methadone (35, 25-140 mg day⁻¹) (median; range) or buprenorphine (6.00, 2-20 mg day⁻¹) during pregnancy. Plasma concentration ratios are presented as an indicator of foetal exposure relative to the mother.Methadone was quantified in all samples, with cord : maternal plasma methadone concentration ratios (n= 15 mother-infant pairs) being significantly higher (P0.0001; mean difference (MD) 0.07; 95% confidence interval (CI) 0.048, 0.092) for the active (R)-methadone enantiomer (0.41; 0.19, 0.56) (median; range) compared with (S)-methadone (0.36; 0.15, 0.53). (R)- : (S)-methadone concentration ratios were also significantly higher (P0.0001; MD 0.24 95% CI 0.300, 0.180) for cord (1.40; 0.95, 1.67) compared with maternal plasma (1.16; 0.81, 1.38). Half the infant buprenorphine samples were below the assay lower limit of quantification (LLOQ) (0.125 ng ml⁻¹). The latter was four-fold lower than the LLOQ for methadone (0.50 ng ml⁻¹). The cord : maternal plasma buprenorphine concentration ratio (n= 9 mother-infant pairs) was 0.35; 0.14, 0.47 and for norbuprenorphine 0.49; 0.24, 0.91.The transfer of the individual methadone enantiomers to the foetal circulation is stereoselective. Infants born to buprenorphine maintained women are not exposed to a greater proportion of the maternal dose compared with methadone and may be exposed to relatively less of the maternal dose compared with infants born to women maintained on methadone during pregnancy.