Fundamentals of population pharmacokinetic modelling: validation methods

Population pharmacokinetics of factor IX in hemophilia B patients undergoing surgery

Essentials Factor IX (FIX) dosing using body weight frequently results in under and overdosing during surgery. We aimed to establish a population pharmacokinetic (PK) model describing the perioperative FIX levels. Population PK parameter values for clearance and V1 were 284 mL h-170 kg-1 and 5450 mL70 kg-1. Perioperative PK parameters differ from those during non-surgical prophylactic treatment. SUMMARY: Background Hemophilia B is a bleeding disorder characterized by a deficiency of coagulation factor IX (FIX). In the perioperative setting, patients receive FIX concentrates to ensure hemostasis. Although FIX is usually dosed according to bodyweight, under- and overdosing occurs frequently during surgery. Aim The objective was to quantify and explain the interpatient variability of perioperatively administered plasma-derived (pd) and recombinant (r) FIX concentrates. Methods Data were collected from 118 patients (median age, 40 years [range, 0.2-90]; weight, 79 kg [range, 5.3-132]) with moderate (28%) or severe hemophilia B (72%), undergoing 255 surgical procedures. Population pharmacokinetic (PK) parameters were estimated using nonlinear mixed-effect modeling in NONMEM. Results Measured perioperative FIX level vs. time profiles were adequately described using a three-compartment PK model. For a typical 34-year-old patient receiving rFIX, clearance (CL), intercompartmental clearance (Q2, Q3), distribution volume of the central compartment (V1) and peripheral compartments (V2, V3) plus interpatient variability (%CV) were: CL, 284 mL h-170 kg-1 (18%); V1, 5450 mL70 kg-1 (19%); Q2, 110 mL h-170 kg-1; V2, 4800 mL70 kg-1; Q3, 1610 mL h-170 kg-1; V3, 2040 mL70 kg-1. From 0.2 years, CL and V1 decreased 0.89% and 1.15% per year, respectively, until the age of 34 years. Patients receiving pdFIX exhibited a lower CL (11%) and V1 (17%) than patients receiving rFIX. Interpatient variability was successfully quantified and explained. Conclusions The estimated perioperative PK parameters of both pdFIX and rFIX are different from those reported for prophylactic treatment. The developed model may be used to apply PK-guided dosing of FIX concentrates during surgery.

Clinical applications of population pharmacokinetic models of antibiotics: Challenges and perspectives

Because of increasing antimicrobial resistance and the shortage of new antibiotics, there is a growing need to optimize the use of old and new antibiotics. Modelling of the pharmacokinetic/pharmacodynamic (PK/PD) characteristics of antibiotics can support the optimization of dosing regimens. Antimicrobial efficacy is determined by susceptibility of the drug to the microorganism and exposure to the drug, which relies on the PK and the dose. Population PK models describe relationships between patients characteristics and drug exposure. This article highlights three clinical applications of these models applied to antibiotics: 1) dosing evaluation of old antibiotics, 2) setting clinical breakpoints and 3) dosing individualization using therapeutic drug monitoring (TDM). For each clinical application, challenges regarding interpretation are discussed. An important challenge is to improve the understanding of the interpretation of modelling results for good implementation of the dosing recommendations, clinical breakpoints and TDM advices. Therefore, also background information on PK/PD principles and approaches to analyse PK/PD data are provided.

Population pharmacokinetics of Rilpivirine in HIV-1-infected patients treated with the single-tablet regimen rilpivirine/tenofovir/emtricitabine

PURPOSE

Rilpivirine, prescribed for the treatment of HIV infection, presents an important inter-individual pharmacokinetic variability. We aimed to determine population pharmacokinetic parameters of rilpivirine in adult HIV-infected patients and quantify their inter-individual variability.

METHODS

We conducted a multicenter, retrospective, and observational study in patients treated with the once-daily rilpivirine/tenofovir disoproxil fumarate/emtricitabine regimen. As part of routine therapeutic drug monitoring, rilpivirine concentrations were measured by UPLC-MS/MS. Population pharmacokinetic analysis was performed using NONMEM software. Once the compartmental and random effects models were selected, covariates were tested to explain the inter-individual variability in pharmacokinetic parameters. The final model qualification was performed by both statistical and graphical methods.

RESULTS

We included 379 patients, resulting in the analysis of 779 rilpivirine plasma concentrations. Of the observed trough individual plasma concentrations, 24.4% were below the 50 ng/ml minimal effective concentration. A one-compartment model with first-order absorption best described the data. The estimated fixed effect for plasma apparent clearance and distribution volume were 9 L/h and 321 L, respectively, resulting in a half-life of 25.2 h. The common inter-individual variability for both parameters was 34.1% at both the first and the second occasions. The inter-individual variability of clearance was 30.3%.

CONCLUSIONS

Our results showed a terminal half-life lower than reported and a high proportion of patients with suboptimal rilpivirine concentrations, which highlights the interest of using therapeutic drug monitoring in clinical practice. The population analysis performed with data from "real-life" conditions resulted in reliable post hoc estimates of pharmacokinetic parameters, suitable for individualization of dosing regimen.

Application of Population Pharmacokinetic Modeling for Individualized Infliximab Dosing Strategies in Crohn Disease

OBJECTIVES

The pharmacokinetics of infliximab (IFX) is highly variable in children with Crohn disease (CD), and a one-size-fits-all approach to dosing is inadequate. Model-based drug dosing can help individualize dosing strategies. We evaluated the predictive performance and clinical utility of a published population pharmacokinetic model of IFX in children with CD.

METHODS

Within a cohort of 34 children with CD who had IFX trough concentrations measured, the pharmacokinetics of each patient was estimated in NONMEM using a published population pharmacokinetic model. Infliximab concentrations were then predicted based on each patient's dosing history and compared with actual measured concentrations (n = 59). In addition, doses 5 to 10 mg/kg and dosing intervals every 4 to 8 weeks were simulated in each patient to examine dose-trough relationships.

RESULTS

Predicted concentrations were within ±1.0 μg/mL of actual measured concentrations for 88% of measurements. The median prediction error (ie, measure of bias) was -0.15 μg/mL (95% confidence interval -0.37 to -0.05 μg/mL) and absolute prediction error (ie, measure of precision) was 0.26 μg/mL (95% confidence interval 0.15 to 0.40 μg/mL). At standard maintenance dosing of 5 mg/kg every 8 weeks, a trough >3 μg/mL was predicted to be achieved in 32% of patients. To achieve a trough >3 μg/mL, a dosing interval ≤every 6 weeks was predicted to be required in 29% of patients.

CONCLUSIONS

A published IFX population pharmacokinetic model demonstrated accurate predictive performance in a pediatric CD population. Individualized IFX dosing strategies in children with CD will be critical to consistently achieve trough concentrations associated with optimal outcomes.

Revolutionizing Therapeutic Drug Monitoring with the Use of Interstitial Fluid and Microneedles Technology

While therapeutic drug monitoring (TDM) that uses blood as the biological matrix is the traditional gold standard, this practice may be impossible, impractical, or unethical for some patient populations (e.g., elderly, pediatric, anemic) and those with fragile veins. In the context of finding an alternative biological matrix for TDM, this manuscript will provide a qualitative review on: (1) the principles of TDM; (2) alternative matrices for TDM; (3) current evidence supporting the use of interstitial fluid (ISF) for TDM in clinical models; (4) the use of microneedle technologies, which is potentially minimally invasive and pain-free, for the collection of ISF; and (5) future directions. The current state of knowledge on the use of ISF for TDM in humans is still limited. A thorough literature review indicates that only a few drug classes have been investigated (i.e., anti-infectives, anticonvulsants, and miscellaneous other agents). Studies have successfully demonstrated techniques for ISF extraction from the skin but have failed to demonstrate commercial feasibility of ISF extraction followed by analysis of its content outside the ISF-collecting microneedle device. In contrast, microneedle-integrated biosensors built to extract ISF and perform the biomolecule analysis on-device, with a key feature of not needing to transfer ISF to a separate instrument, have yielded promising results that need to be validated in pre-clinical and clinical studies. The most promising applications for microneedle-integrated biosensors is continuous monitoring of biomolecules from the skin's ISF. Conducting TDM using ISF is at the stage where its clinical utility should be investigated. Based on the advancements described in the current review, the immediate future direction for this area of research is to establish the suitability of using ISF for TDM in human models for drugs that have been found suitable in pre-clinical experiments.

What is the role for population pharmacokinetics in hemophilia?

Prevention of bleeding in hemophilia requires that plasma levels of the deficient factor exceed the desired minimum target level. Large interindividual variability suggests that knowledge of individual pharmacokinetic (PK) would help to achieve this goal, simultaneously minimizing infusion frequency and the amount of concentrate used. Population PK (PopPK) allows for the incorporation of determinants of interpatient variability and eliminates the need for extensive postinfusion plasma sampling. Barriers to implementation of PopPK are the need for concentrate specific models, Bayesian calculation power, specific expertise for validation and appraisal of forecasted estimates. The Web Accessible Population Pharmacokinetic Service – Hemophilia ( www.wapps-hemo.org ), developed by an international research network of hemophilia centers will test if PK-guided dose individualization can improve patient important outcomes in hemophilia.

Evaluation of Amikacin Pharmacokinetics and Pharmacodynamics for Optimal Initial Dosing Regimen

Amikacin has been one of the important antimicrobial agents against Gram-negative pathogens. However, there is discrepancy regarding the amikacin initial dosage, with some reports recently recommending ≥25 mg/kg and others the conventional dosage (15-20 mg/kg). Hence, we evaluated the optimal initial dosing regimen of amikacin. Pharmacokinetic (PK) parameters were estimated using a population PK analysis. The pharmacodynamic (PD) target was a ratio of ≥8 between the concentration achieved 1 h after beginning the infusion (C peak) and the minimal inhibitory concentration (MIC) of the liable bacteria. Based on the population PK parameters, we simulated individual C peak for several dosing regimens by Monte Carlo method and analyzed the C peak/MIC ratio for MICs from 0.5 to 32 μg/mL. This study included 35 infected patients (25 males), with a median (range) age and body weight of 70 (15-95) years and 49.5 (32.5-78) kg, respectively. A two-compartment model was used, and total body clearance (CL) significantly correlated with creatinine clearance, and volume of distribution (V d) with body weight. Regarding the probability to achieve a C peak/MIC of ≥8, the 15 mg/kg regimen was sufficient to achieve the PK/PD target in ≥90% of patients for a MIC of 4 μg/mL or less. The cumulative fraction of response in Pseudomonas aeruginosa was that 76% of patients achieved a C peak/MIC of 8 with the amikacin dosage of 15 mg/kg/day. We suggest that the 15-mg/kg once-daily dosage of amikacin be recommended as the initial dosage. As its maintenance dosage, the 15 mg/kg/day amikacin dosage is needed for a MIC of ≤4 μg/mL, and amikacin monotherapy for a MIC of ≥8 μg/mL should be avoided.

Multi-scale Modeling in Clinical Oncology: Opportunities and Barriers to Success

Hierarchical processes spanning several orders of magnitude of both space and time underlie nearly all cancers. Multi-scale statistical, mathematical, and computational modeling methods are central to designing, implementing and assessing treatment strategies that account for these hierarchies. The basic science underlying these modeling efforts is maturing into a new discipline that is close to influencing and facilitating clinical successes. The purpose of this review is to capture the state-of-the-art as well as the key barriers to success for multi-scale modeling in clinical oncology. We begin with a summary of the long-envisioned promise of multi-scale modeling in clinical oncology, including the synthesis of disparate data types into models that reveal underlying mechanisms and allow for experimental testing of hypotheses. We then evaluate the mathematical techniques employed most widely and present several examples illustrating their application as well as the current gap between pre-clinical and clinical applications. We conclude with a discussion of what we view to be the key challenges and opportunities for multi-scale modeling in clinical oncology.

Pharmacotherapy during pediatric extracorporeal membrane oxygenation: a review

INTRODUCTION

Pediatric critical illness and associated alterations in organ function can change drug pharmacokinetics (PK). Extracorporeal membrane oxygenation (ECMO), a life-saving therapy for severe cardiac and/or respiratory failure, causes additional PK alterations that affect drug disposition.

AREAS COVERED

The purposes of this review are to discuss the PK changes that occur during ECMO, the associated therapeutic implications, and to review PK literature relevant to pediatric ECMO. We discuss various classes of drugs commonly used for pediatric patients on ECMO, including sedatives, analgesics, antimicrobials and cardiovascular drugs. Finally, we discuss future areas of research and recommend strategies for future pediatric ECMO pharmacologic investigations.

EXPERT OPINION

Clinicians caring for pediatric patients treated with ECMO must have an understanding of PK alterations that could lead to either therapeutic failures or increased drug toxicity during this life-saving therapy. Limited data currently exist for optimal drug dosing in pediatric populations who are treated with ECMO. While there are clear challenges to conducting and analyzing data associated with clinical pharmacokinetic-pharmacodynamic studies of children on ECMO, we present techniques to address these challenges. Improved understanding of the physiology and drug disposition during ECMO combined with PK-PD modeling will allow for more adaptable and individualized dosing schemes.

Assessing Predictive Performance of Published Population Pharmacokinetic Models of Intravenous Tobramycin in Pediatric Patients

Several population pharmacokinetic models describe the dose-exposure relationship of tobramycin in pediatric patients. Before the implementation of these models in clinical practice for dosage adjustment, their predictive performance should be externally evaluated. This study tested the predictive performance of all published population pharmacokinetic models of tobramycin developed for pediatric patients with an independent patient cohort. A literature search was conducted to identify suitable models for testing. Demographic and pharmacokinetic data were collected retrospectively from the medical records of pediatric patients who had received intravenous tobramycin. Tobramycin exposure was predicted from each model. Predictive performance was assessed by visual comparison of predictions to observations, by calculation of bias and imprecision, and through the use of simulation-based diagnostics. Eight population pharmacokinetic models were identified. A total of 269 concentration-time points from 41 pediatric patients with cystic fibrosis were collected for external evaluation. Three models consistently performed best in all evaluations and had mean errors ranging from -0.4 to 1.8 mg/liter, relative mean errors ranging from 4.9 to 29.4%, and root mean square errors ranging from 47.8 to 66.9%. Simulation-based diagnostics supported these findings. Models that allowed a two-compartment disposition generally had better predictive performance than those that used a one-compartment disposition model. Several published models of the pharmacokinetics of tobramycin showed reasonable low levels of bias, although all models seemed to have some problems with imprecision. This suggests that knowledge of typical pharmacokinetic behavior and patient covariate values alone without feedback concentration measurements from individual patients is not sufficient to make precise predictions.

Pharmacokinetic Variability of Daptomycin during Prolonged Therapy for Bone and Joint Infections

The interindividual and intraindividual variabilities in daptomycin pharmacokinetics were investigated in 23 patients (69 pharmacokinetic profiles) who were treated for several months for bone and joint infections. Population daptomycin clearance was significantly influenced by renal function and was significantly higher in male than in female patients. We observed significant intraindividual changes in daptomycin clearance, which were uncorrelated with changes in renal function, suggesting that therapeutic drug monitoring is important in patients receiving prolonged daptomycin therapy.

Pragmatic pharmacology: population pharmacokinetic analysis of fentanyl using remnant samples from children after cardiac surgery

AIMS

One barrier contributing to the lack of pharmacokinetic (PK) data in paediatric populations is the need for serial sampling. Analysis of clinically obtained specimens and data may overcome this barrier. To add evidence for the feasibility of this approach, we sought to determine PK parameters for fentanyl in children after cardiac surgery using specimens and data generated in the course of clinical care, without collecting additional blood samples.

METHODS

We measured fentanyl concentrations in plasma from leftover clinically-obtained specimens in 130 paediatric cardiac surgery patients and successfully generated a PK dataset using drug dosing data extracted from electronic medical records. Using a population PK approach, we estimated PK parameters for this population, assessed model goodness-of-fit and internal model validation, and performed subset data analyses. Through simulation studies, we compared predicted fentanyl concentrations using model-driven weight-adjusted per kg vs. fixed per kg fentanyl dosing.

RESULTS

Fentanyl clearance for a 6.4 kg child, the median weight in our cohort, is 5.7 l h(-1) (2.2-9.2 l h(-1) ), similar to values found in prior formal PK studies. Model assessment and subset analyses indicated the model adequately fit the data. Of the covariates studied, only weight significantly impacted fentanyl kinetics, but substantial inter-individual variability remained. In simulation studies, model-driven weight-adjusted per kg fentanyl dosing led to more consistent therapeutic fentanyl concentrations than fixed per kg dosing.

CONCLUSIONS

We show here that population PK modelling using sparse remnant samples and electronic medical records data provides a powerful tool for assessment of drug kinetics and generation of individualized dosing regimens.

External Validation of a Recently Developed Population Pharmacokinetic Model for Hydromorphone During Postoperative Pain Therapy

BACKGROUND AND OBJECTIVE

We recently developed a new population pharmacokinetic model for hydromorphone in patients including age and bodyweight as covariates. The aim of the present study was to evaluate prospectively the predictive performance of this new model during postoperative pain therapy.

METHODS

This was a prospective, single-blinded, randomized, single-center study with two parallel arms. Fifty patients aged 40-85 years undergoing cardiac surgery involving thoracotomy were enrolled. Hydromorphone was administered postoperatively on the intensive care unit as target controlled infusion (TCI) for patient controlled analgesia (TCI-PCA) using the new pharmacokinetic model, or as conventional patient controlled analgesia (PCA). Arterial blood samples were taken for measurement of the hydromorphone plasma concentration. The predictive performance of the pharmacokinetic model was assessed by the median performance error (MDPE), the median absolute performance error (MDAPE), wobble and divergence. For comparison, the performance indices were also determined for three older models from the literature.

RESULTS

903 plasma concentrations of 41 patients were analyzed. The mean values (95 % CI) of MDPE, MDAPE, wobble and divergence for the new pharmacokinetic model were 11.2 % (3.9 to 18.7 %), 28.5 % (23.9 to 33.0 %), 21.4 % (18.0 to 24.9 %) and -1.6 %/h (-2.3 to -0.8 %/h). When compared with older models from the literature, performance was better with less overshoot after bolus doses.

CONCLUSION

The new pharmacokinetic model of hydromorphone showed a good precision and a better performance than older models. It is therefore suitable for TCI with hydromorphone during postoperative pain therapy.

TRIAL REGISTRATION

EudraCT 2013-002875-16, Clinical Trials NCT02035709.

Model-based clinical pharmacology profiling of ipilimumab in patients with advanced melanoma

Aim

Ipilimumab is a fully human, monoclonal antibody that blocks cytotoxic T-lymphocyte antigen-4. The objective of the present study was to characterize the clinical pharmacology profile of ipilimumab using a population pharmacokinetic (PPK) approach.

Methods

The PPK model was developed using 2095 ipilimumab serum concentration values from 499 patients with unresectable stage III or IV melanoma from four phase II studies, with ipilimumab doses ranging from 0.3 to 10 mg kg⁻¹. The structural PK model was determined by developing a base PPK model. The effect of covariates on model parameters was assessed by a full covariate model, which incorporated all pre-specified covariate-parameter relationships into the base model. The final model was developed by backward elimination, followed by exclusion of covariates determined not to be of clinical relevance to ipilimumab, and was rigorously validated against both internal and external datasets.

Results

Ipilimumab PK was linear and time-invariant, with dose-proportional exposures over the available dose range, yielding a terminal half-life of approximately 15 days. Clearance of ipilimumab increased with increasing body weight and baseline serum lactate dehydrogenase concentrations, but was not affected by age, gender, concomitant budesonide, Eastern Cooperative Oncology Group performance status or prior systemic anticancer therapy. Furthermore, ipilimumab exposure was not affected by moderate renal impairment or mild hepatic impairment.

Conclusions

Ipilimumab concentration–time data were well described by a linear, two compartment, zero order i.v. infusion model. The model confirms that a body weight-normalized dosing regimen is appropriate for ipilimumab therapy in patients with advanced melanoma.

Establishing Good Practices for Exposure-Response Analysis of Clinical Endpoints in Drug Development

This tutorial aims at promoting good practices for exposure–response (E-R) analyses of clinical endpoints in drug development. The focus is on practical aspects of E-R analyses to assist modeling scientists with a process of performing such analyses in a consistent manner across individuals and projects and tailored to typical clinical drug development decisions. This includes general considerations for planning, conducting, and visualizing E-R analyses, and how these are linked to key questions.

Prediction of paraquat exposure and toxicity in clinically ill poisoned patients: a model based approach

AIMS

Paraquat poisoning is a medical problem in many parts of Asia and the Pacific. The mortality rate is extremely high as there is no effective treatment. We analyzed data collected during an ongoing cohort study on self-poisoning and from a randomized controlled trial assessing the efficacy of immunosuppressive therapy in hospitalized paraquat-intoxicated patients. The aim of this analysis was to characterize the toxicokinetics and toxicodynamics of paraquat in this population.

METHODS

A non-linear mixed effects approach was used to perform a toxicokinetic/toxicodynamic population analysis in a cohort of 78 patients.

RESULTS

The paraquat plasma concentrations were best fitted by a two compartment toxicokinetic structural model with first order absorption and first order elimination. Changes in renal function were used for the assessment of paraquat toxicodynamics. The estimates of toxicokinetic parameters for the apparent clearance, the apparent volume of distribution and elimination half-life were 1.17 l h(-1) , 2.4 l kg(-1) and 87 h, respectively. Renal function, namely creatinine clearance, was the most significant covariate to explain between patient variability in paraquat clearance.This model suggested that a reduction in paraquat clearance occurred within 24 to 48 h after poison ingestion, and afterwards the clearance was constant over time. The model estimated that a paraquat concentration of 429 μg l(-1) caused 50% of maximum renal toxicity. The immunosuppressive therapy tested during this study was associated with only 8% improvement of renal function.

CONCLUSION

The developed models may be useful as prognostic tools to predict patient outcome based on patient characteristics on admission and to assess drug effectiveness during antidote drug development.

Population pharmacokinetic modelling of recombinant factor IX Fc fusion protein (rFIXFc) in patients with haemophilia B

Background and Objectives

Recombinant factor IX Fc fusion protein (rFIXFc) is a clotting factor developed using monomeric Fc fusion technology to prolong the circulating half-life of factor IX. The objective of this analysis was to elucidate the pharmacokinetic characteristics of rFIXFc in patients with haemophilia B and identify covariates that affect rFIXFc disposition.

Methods

Population pharmacokinetic analysis using NONMEM^® was performed with clinical data from two completed trials in previously treated patients with severe to moderate haemophilia B. Twelve patients from a phase 1/2a study and 123 patients from a registrational phase 3 study were included in this population analysis.

Results

A three-compartment model was found to best describe the pharmacokinetics of rFIXFc. For a typical 73 kg patient, the clearance (CL), volume of the central compartment (V₁) and volume of distribution at steady state (V_ss) were 2.39 dL/h, 71.4 dL and 198 dL, respectively. Because of repeat pharmacokinetic profiles at week 26 for patients in a subgroup, inclusion of inter-occasion variability (IOV) on CL and V₁ were evaluated and significantly improved the model. The magnitude of IOV on CL and V₁ were both low to moderate (<20 %) and less than the corresponding inter-individual variability. Body weight (BW) was found to be the only significant covariate for rFIXFc disposition. However, the impact of BW was limited, as the BW power exponents on CL and V₁ were 0.436 and 0.396, respectively.

Conclusion

This is the first population pharmacokinetic analysis that systematically characterized the pharmacokinetics of long-lasting rFIXFc in patients with haemophilia B. The population pharmacokinetic model for rFIXFc can be utilized to evaluate and optimize dosing regimens for the treatment of patients with haemophilia B.

Electronic supplementary material

The online version of this article (doi:10.1007/s40262-013-0129-7) contains supplementary material, which is available to authorized users.

Absence of P-glycoprotein transport in the pharmacokinetics and toxicity of the herbicide paraquat

Genetic variation in the multidrug resistance gene ABCB1, which encodes the efflux transporter P-glycoprotein (P-gp), has been associated with Parkinson disease. Our goal was to investigate P-gp transport of paraquat, a Parkinson-associated neurotoxicant. We used in vitro transport models of ATPase activity, xenobiotic-induced cytotoxicity, transepithelial permeability, and rhodamine-123 inhibition. We also measured paraquat pharmacokinetics and brain distribution in Friend leukemia virus B-type (FVB) wild-type and P-gp-deficient (mdr1a(-/-)/mdr1b(-/-)) mice following 10, 25, 50, and 100 mg/kg oral doses. In vitro data showed that: 1) paraquat failed to stimulate ATPase activity; 2) resistance to paraquat-induced cytotoxicity was unchanged in P-gp-expressing cells in the absence or presence of P-gp inhibitors GF120918 [N-(4-[2-(1,2,3,4-tetrahydro-6,7-dimethoxy-2-isoquinolinyl)ethyl]-phenyl)-9,10-dihydro-5-methoxy-9-oxo-4-acridine carboxamide] and verapamil-37.0 [95% confidence interval (CI): 33.2-41.4], 46.2 (42.5-50.2), and 34.1 µM (31.2-37.2)-respectively; 3) transepithelial permeability ratios of paraquat were the same in P-gp-expressing and nonexpressing cells (1.55 ± 0.39 and 1.39 ± 0.43, respectively); and 4) paraquat did not inhibit rhodamine-123 transport. Population pharmacokinetic modeling revealed minor differences between FVB wild-type and mdr1a(-/-)/mdr1b(-/-) mice: clearances of 0.47 [95% confidence interval (CI): 0.42-0.52] and 0.78 l/h (0.58-0.98), respectively, and volume of distributions of 1.77 (95% CI: 1.50-2.04) and 3.36 liters (2.39-4.33), respectively; however, the change in clearance was in the opposite direction of what would be expected. It is noteworthy that paraquat brain-to-plasma partitioning ratios and total brain accumulation were the same across doses between FVB wild-type and mdr1a(-/-)/mdr1b(-/-) mice. These studies indicate that paraquat is not a P-gp substrate. Therefore, the association between ABCB1 pharmacogenomics and Parkinson disease is not attributed to alterations in paraquat transport.

Population pharmacokinetic and pharmacodynamic modeling of different formulations of ONO-5334, cathepsin K inhibitor, in Caucasian and Japanese postmenopausal females

ONO-5334, a selective inhibitor of cathepsin K, is a potential new treatment for osteoporosis. The objectives of this study were to (1) develop population pharmacokinetic-pharmacodynamic (PK-PD) models for ONO-5334 using dose-ascending data from healthy postmenopausal females, (2) examine comparability of PK and/or PD profile between Caucasian and Japanese, and (3) compare PK-PD profile between immediate release tablet (IRT) and sustained release tablet (SRT). The population PK-PD models were developed for each formulation for post-dose levels of bone resorption markers (serum CTX and NTX). The data were provided from 4 phase 1 studies with total of 201 Caucasian and 94 Japanese subjects. Plasma concentrations of ONO-5334 and bone resorption markers were thoroughly evaluated in those studies. An indirect response model described relationships between bone resorption markers and plasma concentrations of ONO-5334. There was no significant difference in PK and pharmacodynamic potency (IC50 ) between Caucasian and Japanese. Based on the developed model, serum CTX and NTX after administration of ONO-5334 IRT or SRT were simulated, and the results showed that ONO-5334 SRT would provide comparable PD effect on bone resorption markers with lower dose relative to IRT.

Establishing best practices and guidance in population modeling: an experience with an internal population pharmacokinetic analysis guidance

This tutorial describes the development of a population pharmacokinetic (Pop PK) analysis guidance within Pfizer, which strives for improved consistency and efficiency, and a more systematic approach to model building. General recommendations from the Pfizer internal guidance and a suggested workflow for Pop PK model building are discussed. A description is also provided for mechanisms by which conflicting opinions were captured and resolved across the organization to arrive at the final guidance.

CPT: Pharmacometrics & Systems Pharmacology (2013) 2, e51; doi:10.1038/psp.2013.26; advance online publication 3 July 2013