Training the next generation of pharmacometric modelers: a multisector perspective.

The current demand for pharmacometricians outmatches the supply provided by academic institutions and considerable investments are made to develop the competencies of these scientists on-the-job. Even with the observed increase in academic programs related to pharmacometrics, this need is unlikely to change in the foreseeable future, as the demand and scope of pharmacometrics applications keep expanding. Further, the field of pharmacometrics is changing. The field largely started when Lewis Sheiner and Stuart Beal published their seminal papers on population pharmacokinetics in the late 1970's and early 1980's and has continued to grow in impact and use since its inception. Physiological-based pharmacokinetics and systems pharmacology have grown rapidly in scope and impact in the last decade and machine learning is just on the horizon. While all these methodologies are categorized as pharmacometrics, no one person can be an expert in everything. So how do you train future pharmacometricians? Leading experts in academia, industry, contract research organizations, clinical medicine, and regulatory gave their opinions on how to best train future pharmacometricians. Their opinions were collected and synthesized to create some general recommendations.

A Physiological-Based Pharmacokinetic Model Embedded with a Target-Mediated Drug Disposition Mechanism Can Characterize Single-Dose Warfarin Pharmacokinetic Profiles in Subjects with Various CYP2C9 Genotypes under Different Cotreatments.

Warfarin, a commonly prescribed oral anticoagulant medication, is highly effective in treating deep vein thrombosis and pulmonary embolism. However, the clinical dosing of warfarin is complicated by high interindividual variability in drug exposure and response and its narrow therapeutic index. CYP2C9 genetic polymorphism and drug-drug interactions (DDIs) are substantial contributors to this high variability of warfarin pharmacokinetics (PK), among numerous factors. Building a physiology-based pharmacokinetic (PBPK) model for warfarin is not only critical for a mechanistic characterization of warfarin PK but also useful for investigating the complicated dose-exposure relationship of warfarin. Thus, the objective of this study was to develop a PBPK model for warfarin that integrates information regarding CYP2C9 genetic polymorphisms and their impact on DDIs. Generic PBPK models for both S- and R-warfarin, the two enantiomers of warfarin, were constructed in R with the mrgsolve package. As expected, a generic PBPK model structure did not adequately characterize the warfarin PK profile collected up to 15 days following the administration of a single oral dose of warfarin, especially for S-warfarin. However, following the integration of an empirical target-mediated drug disposition (TMDD) component, the PBPK-TMDD model well characterized the PK profiles collected for both S- and R-warfarin in subjects with different CYP2C9 genotypes. Following the integration of enzyme inhibition and induction effects, the PBPK-TMDD model also characterized the PK profiles of both S- and R-warfarin in various DDI settings. The developed mathematic framework may be useful in building algorithms to better inform the clinical dosing of warfarin. SIGNIFICANCE STATEMENT: The present study found that a traditional physiology-based pharmacokinetic (PBPK) model cannot sufficiently characterize the pharmacokinetic profiles of warfarin enantiomers when warfarin is administered as a single dose, but a PBPK model with a target-mediated drug disposition mechanism can. After incorporating CYP2C9 genotypes and drug-drug interaction information, the developed model is anticipated to facilitate the understanding of warfarin disposition in subjects with different CYP2C9 genotypes in the absence and presence of both cytochrome P450 inhibitors and cytochrome P450 inducers.

Population pharmacokinetics, pharmacodynamics and pharmacogenetics modelling of oxypurinol in Hmong adults with gout and/or hyperuricemia.

AIMS: The aim of this study was to quantify identifiable sources of variability, including key pharmacogenetic variants in oxypurinol pharmacokinetics and their pharmacodynamic effect on serum urate (SU). METHODS: Hmong participants (n = 34) received 100 mg allopurinol twice daily for 7 days followed by 150 mg allopurinol twice daily for 7 days. A sequential population pharmacokinetic pharmacodynamics (PKPD) analysis with non-linear mixed effects modelling was performed. Allopurinol maintenance dose to achieve target SU was simulated based on the final PKPD model. RESULTS: A one-compartment model with first-order absorption and elimination best described the oxypurinol concentration-time data. Inhibition of SU by oxypurinol was described with a direct inhibitory Emax model using steady-state oxypurinol concentrations. Fat-free body mass, estimated creatinine clearance and SLC22A12 rs505802 genotype (0.32 per T allele, 95% CI 0.13, 0.55) were found to predict differences in oxypurinol clearance. Oxypurinol concentration required to inhibit 50% of xanthine dehydrogenase activity was affected by PDZK1 rs12129861 genotype (-0.27 per A allele, 95% CI -0.38, -0.13). Most individuals with both PDZK1 rs12129861 AA and SLC22A12 rs505802 CC genotypes achieve target SU (with at least 75% success rate) with allopurinol below the maximum dose, regardless of renal function and body mass. In contrast, individuals with both PDZK1 rs12129861 GG and SLC22A12 rs505802 TT genotypes would require more than the maximum dose, thus requiring selection of alternative medications. CONCLUSIONS: The proposed allopurinol dosing guide uses individuals' fat-free mass, renal function and SLC22A12 rs505802 and PDZK1 rs12129861 genotypes to achieve target SU.

Investigating the contribution of residual unexplained variability components on bias and imprecision of parameter estimates in population pharmacokinetic mixed-effects modeling.

In a nonlinear mixed-effects modeling (NLMEM) approach of pharmacokinetic (PK) and pharmacodynamic (PD) data, two levels of random effects are generally modeled: between-subject variability (BSV) and residual unexplained variability (RUV). The goal of this simulation-estimation study was to investigate the extent to which PK and RUV model misspecification, errors in recording dosing and sampling times, and variability in drug content uniformity contribute to the estimated magnitude of RUV and PK parameter bias. A two-compartment model with first-order absorption and linear elimination was simulated as a true model. PK parameters were clearance 5.0 L/h; central volume of distribution 35 L; inter-compartmental clearance 50 L/h; peripheral volume of distribution 50 L. All parameters were assumed to have a 30% coefficient of variation (CV). One hundred in-silico subjects were administered a labeled dose of 120 mg under 4 sample collection designs. PK and RUV model misspecifications were associated with relatively larger increases in the magnitude of RUV compared to other sources for all levels of sampling design. The contribution of dose and dosing time misspecifications have negligible effects on RUV but result in higher bias in PK parameter estimates. Inaccurate sampling time data results in biased RUV and increases with the magnitude of perturbations. Combined perturbation scenarios in the studied sources will propagate the variability and accumulate in RUV magnitude and bias in PK parameter estimates. This work provides insight into the potential contributions of many factors that comprise RUV and bias in PK parameters.

Pharmacokinetic Modeling of Warfarin ІI - Model-based Analysis of Warfarin Metabolites following Warfarin Administered either Alone or Together with Fluconazole or Rifampin.

The objective of this study is to conduct a population pharmacokinetic (PK) model-based analysis on 10 warfarin metabolites (4'-, 6-, 7-, 8- and 10-hydroxylated (OH)-S- and R- warfarin), when warfarin is administered alone or together with either fluconazole or rifampin. One or two compartment PK models expanded from target mediated drug disposition (TMDD) models developed previously for warfarin enantiomers were able to sufficiently characterize the PK profiles of 10 warfarin metabolites in plasma and urine under different conditions. Model-based analysis shows CYP2C9 mediated metabolic elimination pathways are more inhibitable by fluconazole (% formation CL (CLf) of 6- and 7-OH-S-warfarin decrease: 73.2% and 74.8%) but less inducible by rifampin (% CLf of 6- and 7-OH-S-warfarin increase: 85% and 75%), compared with non-CYP2C9 mediated elimination pathways (% CLf of 10-OH-S-warfarin and CLR of S-warfarin decrease in the presence of fluconazole: 65.0% and 15.3%; % CLf of 4'- 8- and 10-OH-S-warfarin increase in the presence of rifampin: 260%, 127% and 355%), which potentially explains the CYP2C9 genotype-dependent DDIs exhibited by S-warfarin, when warfarin is administrated together with fluconazole or rifampin. Additionally, for subjects with CYP2C9 *2 and *3 variants, a model-based analysis of warfarin metabolite profiles in subjects with various CYP2C9 genotypes demonstrates CYP2C9 mediated elimination is less important and non-CYP2C9 mediated elimination is more important, compared with subjects without these variants. To our knowledge, this is so far one of the most comprehensive population-based PK analyses of warfarin metabolites in subjects with various CYP2C9 genotypes under different co-medications. Significance Statement The studies we wish to publish are potentially impactful. The need for a TMDD pharmacokinetic model and the demonstration of genotyped-dependent drug interactions may explain the extensive variability in dose-response relationships that are seen in the clinical dose adjustments of warfarin.

Pharmacokinetic Modeling of Warfarin І - Model-based Analysis of Warfarin Enantiomers with a Target Mediated Drug Disposition Model Reveals CYP2C9 Genotype-dependent Drug-drug Interactions of S-Warfarin.

The objective of this study is to characterize the impact of CYP2C9 genotype on warfarin drug-drug interactions when warfarin is taken together with fluconazole, a cytochrome P450 (CYP) inhibitor, or rifampin, a CYP inducer with a nonlinear mixed effect modeling approach. A target mediated drug disposition model with a urine compartment was necessary to characterize both S-warfarin and R-warfarin plasma and urine pharmacokinetic profiles sufficiently. Following the administration of fluconazole, our study found subjects with CYP2C9 *2 or *3 alleles experience smaller changes in S-warfarin CL compared with subjects without these alleles (69.5%, 64.8%, 59.7% and 47.8% decrease in subjects with CYP2C9 *1/*1, *1/*3, *2/*3 and *3/*3 respectively). Whereas, following the administration of rifampin, subjects with CYP2C9 *2/*3 or CYP2C9 *3/*3 experience larger changes in S-warfarin CL compared with subjects with at least one copy of CYP2C9 *1 or *1B (115%, 111%, 119%, 198% and 193% increase in subjects with CYP2C9 *1/*1, *1B/*1B, *1/*3, *2/*3 and *3/*3 respectively). The results suggest different dose adjustments are potentially required for patients with different CYP2C9 genotypes if warfarin is administered together with CYP inhibitors or inducers. Significance Statement The present study found a target mediated drug disposition model is needed to sufficiently characterize the clinical pharmacokinetic profiles of warfarin racemates under different co-treatments in subjects with various CYP2C9 genotypes, following a single dose of warfarin administration. The study also found S-warfarin, the pharmacologically more active ingredient in warfarin, exhibits CYP2C9 genotype-dependent drug-drug interactions, which indicates the dose of warfarin may need to be adjusted differently in subjects with different CYP2C9 genotypes in the presence of drug-drug interactions.

An integrated PK-PD model for cortisol and the 17-hydroxyprogesterone and androstenedione biomarkers in children with congenital adrenal hyperplasia.

AIMS: The aim of this study was to characterize the pharmacokinetic/pharmacodynamic relationships of cortisol and the adrenal biomarkers 17-hydroxyprogesterone and androstenedione in children with congenital adrenal hyperplasia (CAH). METHODS: A nonlinear mixed-effect modelling approach was used to analyse cortisol, 17-hydroxyprogesterone and androstenedione concentrations obtained over 6 hours from children with CAH (n = 50). A circadian rhythm was evident and the model leveraged literature information on circadian rhythm in untreated children with CAH. Indirect response models were applied in which cortisol inhibited the production rate of all three compounds using an Imax model. RESULTS: Cortisol was characterized by a one-compartment model with apparent clearance and volume of distribution estimated at 22.9 L/h/70 kg and 41.1 L/70 kg, respectively. The IC50 values of cortisol concentrations for cortisol, 17-hydroxyprogesterone and androstenedione were estimated to be 1.36, 0.45 and 0.75 µg/dL, respectively. The inhibitory effect was found to be more potent on 17OHP than D4A, and the IC50 values were higher in salt-wasting subjects than simple virilizers. Production rates of cortisol, 17-hydroxyprogesterone and androstenedione were higher in simple-virilizer subjects. Half-lives of cortisol, 17-hydroxyprogesterone and androstenedione were 60, 47 and 77 minutes, respectively. CONCLUSION: Rapidly changing biomarker responses to cortisol concentrations highlight that single measurements provide volatile information about a child's disease control. Our model closely captured observed cortisol, 17-hydroxyprogesterone and androstenedione concentrations. It can be used to predict concentrations over 24 hours and allows many novel exposure metrics to be calculated, e.g., AUC, AUC-above-threshold, time-within-range, etc. Our long-range goal is to uncover dose-exposure-outcome relationships that clinicians can use in adjusting hydrocortisone dose and timing.

Pharmacokinetic/pharmacodynamic data extrapolation models for improved pediatric efficacy and toxicity estimation, with application to secondary hyperparathyroidism.

In most drug development settings, the regulatory approval process is accompanied by extensive studies performed to understand the drug's pharmacokinetic (PK) and pharmacodynamic (PD) properties. In this article, we attempt to utilize the rich PK/PD data to inform the borrowing of information from adults during pediatric drug development. In pediatric settings, it is especially crucial that we are parsimonious with the patients recruited for experimentation. We illustrate our approaches in the context of clinical trials of cinacalcet for treating secondary hyperparathyroidism in pediatric and adult patients with chronic kidney disease, where we model both parathyroid hormone (efficacy endpoint) and corrected calcium levels (safety endpoint). We use population PK/PD modeling of the cinacalcet data to quantitatively assess the similarity between adults and children, and use this information in various hierarchical Bayesian adult borrowing rules whose statistical properties can then be evaluated. In particular, we simulate the bias and mean square error performance of our approaches in settings where borrowing is and is not warranted to inform guidelines for the future use of our methods.

Brain Distributional Kinetics of a Novel MDM2 Inhibitor SAR405838: Implications for Use in Brain Tumor Therapy.

Achieving an effective drug concentration in the brain is as important as targeting the right pathway when developing targeted agents for brain tumors. SAR405838 is a novel molecularly targeted agent that is in clinical trials for various solid tumors. Its application for tumors in the brain has not yet been examined, even though the target, the MDM2-p53 interaction, is attractive for tumors that could occur in the brain, including glioblastoma and brain metastases. In vitro and in vivo studies indicate that SAR405838 is a substrate of P-glycoprotein (P-gp). P-gp mediated active efflux at the blood-brain barrier plays a dominant role in limiting SAR405838 brain distribution. Even though the absence of P-gp significantly increases the drug exposure in the brain, the systemic exposure, including absorption and clearance processes, were unaffected by P-gp deletion. Model-based parameters of SAR405838 distribution across the blood-brain barrier indicate the CLout of the brain was approximately 40-fold greater than the CLin The free fraction of SAR405838 in plasma and brain were found to be low, and subsequent Kpuu values were less than unity, even in P-gp/Bcrp knockout mice. These results indicate additional efflux transporters other than P-gp and Bcrp may be limiting distribution of SAR405838 to the brain. Concomitant administration of elacridar significantly increased brain exposure, also without affecting the systemic exposure. This study characterized the brain distributional kinetics of SAR405838, a novel MDM2 inhibitor, to evaluate its potential in the treatment of primary and metastatic brain tumors. SIGNIFICANCE STATEMENT: This paper examined the brain distributional kinetics of a novel MDM2-p53 targeted agent, SAR405838, to see its possible application for brain tumors by using in vitro, in vivo, and in silico approaches. SAR405838 is found to be a substrate of P-glycoprotein (P-gp), which limits its distribution to the brain. Based on the findings in the paper, manipulation of the function of P-gp can significantly increase the brain exposure of SAR405838, which may give an insight on its potential benefit as a treatment for primary and metastatic brain cancer.

Orphan drug development: the increasing role of clinical pharmacology.

Over the last few decades there has been a paradigm shift in orphan drug research and development. The development of the regulatory framework, establishment of rare disease global networks that support drug developments, and advances in technology, has resulted in tremendous growth in orphan drug development. Nevertheless, several challenges during orphan drug development such as economic constraints; insufficient clinical information; fewer patients and thus inadequate power; etc. still exist. While the standard regulatory requirements for drug approval stays the same, applications of scientific judgment and regulatory flexibility is significantly important to help meeting some of the immense unmet medical need in rare diseases. Clinical pharmacology presents a vital role in accelerating orphan drug development and overcoming some of these challenges. This review highlights the critical contributions of clinical pharmacology in orphan drug development; for example, dose finding, optimizing clinical trial design, indication expansion, and population extrapolation. Examples of such applications are reviewed in this article.

Pharmacokinetics-pharmacodynamics of sertraline as an antifungal in HIV-infected Ugandans with cryptococcal meningitis.

The ASTRO-CM dose-finding pilot study investigated the role of adjunctive sertraline for the treatment of HIV-associated cryptococcal meningitis in HIV-infected Ugandan patients. The present study is a post hoc pharmacokinetic-pharmacodynamic analysis of the ASTRO-CM pilot study to provide insight into sertraline exposure-response-outcome relationships. We performed a population pharmacokinetic analysis using sertraline plasma concentration data and correlated various predicted PK-PD indices with the percentage change in log10 CFU/mL from baseline. Sertraline clearance was 1.95-fold higher in patients receiving antiretroviral (ART), resulting in 49% lower drug exposure. To quantify the clinical benefit of sertraline, we estimated rates of fungal clearance from cerebrospinal fluid (CSF) of ASTRO-CM patients using Poisson model and compared the clearance rates to a historical control study (COAT) in which patients received standard Cryptococcus therapy of amphotericin B (0.7-1.0 mg/kg per day) and fluconazole (800 mg/day) without sertraline. Adjunctive sertraline significantly increased CSF fungal clearance rate compared to COAT trial and sertraline effect was dose-independent with no covariate found to affect fungal clearance including ART. Study findings suggest sertraline response could be mediated by different mechanisms than directly inhibiting the initiation of protein translation as previously suggested; this is supported by the prediction of unbound sertraline concentrations is unlikely to reach MIC concentrations in the brain. Study findings also recommend against the use of higher doses of sertraline, especially those greater than the maximum FDA-approved daily dose (200 mg/day), since they unlikely provide any additional benefits and come with greater costs and risk of adverse events.

Amikacin Pharmacokinetic-Pharmacodynamic Analysis in Pediatric Cancer Patients.

We performed pharmacokinetic-pharmacodynamic (PK-PD) and simulation analyses to evaluate a standard amikacin dose of 15 mg/kg once daily in children with cancer and to determine an optimal dosing strategy. A population pharmacokinetic model was developed from clinical data collected in 34 pediatric patients and used in a simulation study to predict the population probability of various dosing regimens to achieve accepted safety (steady-state unbound trough plasma concentration [fCmin] of <10 mg/liter)- and efficacy (free, unbound plasma concentration-to-MIC ratio [fCmax/MIC] of ≥8)-linked targets. In addition, an adaptive resistance PD (ARPD) model of Pseudomonas aeruginosa was built based on literature time-kill curve data and linked to the PK model to perform PK-ARPD simulations and compare results with those of the probability approach. Using the probability approach, an amikacin dose of 60 mg/kg administered once daily is expected to achieve the target fCmax/MIC in 80% of pediatric patients weighing 8 to 70 kg with a 97.5% probability, and almost all patients were predicted to have fCmin of <10 mg/liter. However, PK-ARPD simulation predicted that 60 mg/kg/day is unlikely to suppress bacterial resistance with repeated dosing. Furthermore, PK-ARPD simulation suggested that amikacin at 90 mg/kg, given in two divided doses (45 mg/kg twice a day), is expected to hit safety and efficacy targets and is associated with a lower rate of bacterial resistance. The disagreement between the two methods is due to the inability of the probability approach to predict development of drug resistance with repeated dosing. This originates from the use of PK-PD indices based on the MIC that neglects measurement errors, ignores the time course dynamic nature of bacterial growth and killing, and incorrectly assumes the MIC to be constant during treatment.

Lamotrigine pharmacokinetics following oral and stable-labeled intravenous administration in young and elderly adult epilepsy patients: Effect of age.

OBJECTIVE: The objectives of this study were to investigate the effect of age on pharmacokinetic parameters of lamotrigine (LTG) and estimate parameter variability. METHODS: Patients (>18 years old) who were already on a steady-state dose of LTG therapy with no interacting comedications were enrolled. Patients with significant cardiac disease, severe kidney dysfunction, or moderate-to-severe liver dysfunction were excluded. Fifty milligrams of a stable-labeled intravenous LTG formulation (SL-LTG) replaced 50 mg of a patient's normal daily oral LTG dose. Thirteen blood samples were collected in each person over 96 hours. SL-LTG and unlabeled LTG concentrations were measured simultaneously by gas chromatography-mass spectrometry. Concentration-time data were analyzed by nonlinear mixed-effects modeling (NONMEM version 7.3). RESULTS: Twenty-eight patients representing 16 young (18-48 years old) and 12 elderly (63-87 years old) patients were included, yielding 382 unlabeled and 351 SL-LTG concentrations. A two-compartment model with first-order absorption and elimination adequately described the plasma concentration-time data. Bioavailability of oral LTG was approximately 74% and did not differ by age. LTG clearance was 27.2% lower in elderly than in young patients (1.80 L/h for a 70-kg patient). SIGNIFICANCE: Although LTG bioavailability was not affected by age, LTG clearance was 27.2% lower in elderly versus young patients of comparable body weight, possibly indicating lower dosages being needed in this population.

Pharmacokinetic-pharmacodynamic modelling of acute N-terminal pro B-type natriuretic peptide after doxorubicin infusion in breast cancer.

AIMS: The aim of the present study was to develop a pharmacokinetic-pharmacodynamic (PK-PD) model to characterize the relationship between plasma doxorubicin and N-terminal pro B-type natriuretic peptide (NT-proBNP) concentrations within 48 h of doxorubicin treatment. METHODS: The study enrolled 17 female patients with stages 1-3 breast cancer and receiving adjuvant doxorubicin (60 mg m(-2) ) and cyclophosphamide (600 mg m(-2) ) every 14 days for four cycles. In two consecutive cycles, plasma concentrations of doxorubicin, doxorubicinol, troponin and NT-proBNP were collected before infusion, and up to 48 h after the end of doxorubicin infusion. Nonlinear mixed-effects modelling was used to describe the PK-PD relationship of doxorubicin and NT-proBNP. RESULTS: A three-compartment parent drug with a one-compartment metabolite model best described the PK of doxorubicin and doxorubicinol. Troponin concentrations remained similar to baseline. An indirect PD model with transit compartments best described the relationship of doxorubicin exposure and acute NT-proBNP response. Estimated PD parameters were associated with large between-subject variability (total assay variability 38.8-73.9%). Patient clinical factors, including the use of enalapril, were not observed to be significantly associated with doxorubicin PK or NT-proBNP PD variability. CONCLUSION: The relationship between doxorubicin concentration and the acute NT-proBNP response was successfully described with a population PK-PD model. This model will serve as a valuable framework for future studies to identify clinical factors associated with the acute response to doxorubicin. Future studies are warranted to examine the relationship between this acute response and subsequent heart failure. Should such a relationship be established, this model could provide useful information on patients' susceptibility to doxorubicin-induced long-term cardiotoxicity.

A model-based approach to assess the exposure-response relationship of Lorenzo's oil in adrenoleukodystrophy.

AIMS: X-linked adrenoleukodystrophy (X-ALD) is a peroxisomal disorder, most commonly affecting boys, associated with increased very long chain fatty acids (C26:0) in all tissues, causing cerebral demyelination and adrenocortical insufficiency. Certain monounsaturated long chain fatty acids including oleic and erucic acids, known as Lorenzo's oil (LO), lower plasma C26:0 levels. The aims of this study were to characterize the effect of LO administration on plasma C26:0 concentrations and to determine whether there is an association between plasma concentrations of erucic acid or C26:0 and the likelihood of developing brain MRI abnormalities in asymptomatic boys. METHODS: Non-linear mixed effects modelling was performed on 2384 samples collected during an open label single arm trial. The subjects (n = 104) were administered LO daily at ~2-3 mg kg(-1) with a mean follow-up of 4.88 ± 2.76 years. The effect of erucic acid exposure on plasma C26:0 concentrations was characterized by an inhibitory fractional Emax model. A Weibull model was used to characterize the time-to-developing MRI abnormality. RESULTS: The population estimate for the fractional maximum reduction of C26:0 plasma concentrations was 0.76 (bootstrap 95% CI 0.73, 0.793). Our time-to-event analyses showed that every mg l(-1) increase in time-weighted average of erucic acid and C26:0 plasma concentrations was, respectively, associated with a 3.7% reduction and a 753% increase in the hazard of developing MRI abnormality. However, the results were not significant (P = 0.5344, 0.1509, respectively). CONCLUSIONS: LO administration significantly reduces the abnormally high plasma C26:0 concentrations in X-ALD patients. Further studies to evaluate the effect of LO on the likelihood of developing brain MRI abnormality are warranted.

A hierarchical Bayesian approach for combining pharmacokinetic/pharmacodynamic modeling and Phase IIa trial design in orphan drugs: Treating adrenoleukodystrophy with Lorenzo's oil.

X-linked adrenoleukodystrophy (X-ALD) is a rare, progressive, and typically fatal neurodegenerative disease. Lorenzo's oil (LO) is one of the few X-ALD treatments available, but little has been done to establish its clinical efficacy or indications for its use. In this article, we analyze data on 116 male asymptomatic pediatric patients who were administered LO. We offer a hierarchical Bayesian statistical approach to understand LO pharmacokinetics (PK) and pharmacodynamics (PD) resulting from an accumulation of very long-chain fatty acids. We experiment with individual- and observational-level errors and various choices of prior distributions and deal with the limitation of having just one observation per administration of the drug, as opposed to the more usual multiple observations per administration. We link LO dose to the plasma erucic acid concentrations by PK modeling, and then link this concentration to a biomarker (C26, a very long-chain fatty acid) by PD modeling. Next, we design a Bayesian Phase IIa study to estimate precisely what improvements in the biomarker can arise from various LO doses while simultaneously modeling a binary toxicity endpoint. Our Bayesian adaptive algorithm emerges as reasonably robust and efficient while still retaining good classical (frequentist) operating characteristics. Future work looks toward using the results of this trial to design a Phase III study linking LO dose to actual improvements in health status, as measured by the appearance of brain lesions observed via magnetic resonance imaging.

Pharmacokinetic-pharmacodynamic modelling of intravenous and oral topiramate and its effect on phonemic fluency in adult healthy volunteers.

AIMS: The aim was to develop a quantitative approach that characterizes the magnitude of and variability in phonemic generative fluency scores as measured by the Controlled Oral Word Association (COWA) test in healthy volunteers after administration of an oral and a novel intravenous (IV) formulation of topiramate (TPM). METHODS: Nonlinear mixed-effects modelling was used to describe the plasma TPM concentrations resulting from oral or IV administration. A pharmacokinetic-pharmacodynamic (PK-PD) model was developed sequentially to characterize the effect of TPM concentrations on COWA with different distributional assumptions. RESULTS: Topiramate was rapidly absorbed, with a median time to maximal concentration of 1 h and an oral bioavailability of ~100%. Baseline COWA score increased by an average of 12% after the third administration on drug-free sessions. An exponential model described the decline of COWA scores, which decreased by 14.5% for each 1 mg l(-1) increase in TPM concentration. The COWA scores were described equally well by both continuous normal and Poisson distributions. CONCLUSIONS: This analysis quantified the effect of TPM exposure on generative verbal fluency as measured by COWA. Repetitive administration of COWA resulted in a better performance, possibly due to a learning effect. The model predicts a 27% reduction in the COWA score at the average observed maximal plasma concentration after a 100 mg dose of TPM. The single-dose administration of relatively low TPM doses and narrow range of resultant concentrations in our study were limitations to investigating the PK-PD relationship at higher TPM exposures. Hence, the findings may not be readily generalized to the broader patient population.

Pharmacometric characterization of efavirenz developmental pharmacokinetics and pharmacogenetics in HIV-infected children.

The aim of this analysis was to create a pharmacometric model of efavirenz developmental pharmacokinetics and pharmacogenetics in HIV-infected children. The data consisted of 3,172 plasma concentrations from 96 HIV-1-infected children who participated in the Pediatric AIDS Clinical Trials Group 382 (PACTG382) study. Analyses were performed using NONMEM, and the impacts of body weight, age, race, sex, formulation, liver function, and cytochrome P450 2B6 (CYP2B6)-G516T and multidrug-resistance transporter gene (MDR1)-C3435T polymorphisms were assessed. A one-compartment model using weight-based allometry on oral clearance and apparent volume of distribution adequately described the data. A sigmoid maximum-effect (Emax) maturation model demonstrated an increase in oral clearance with age to reach 90% of its mature level by the age of 9 months. The liquid formulation bioavailability relative to the capsule was found to increase with age to reach 90% of its mature value by the age of 8 years. The CYP2B6-G516T polymorphism decreased oral clearance, while the MDR1-C3435T polymorphism demonstrated no effect.

Population pharmacokinetics of valganciclovir prophylaxis in paediatric and adult solid organ transplant recipients.

AIM: Our aims were to quantify ganciclovir pharmacokinetics in paediatric and adult kidney, liver and lung transplant patients taking a range of valganciclovir doses to prevent herpes virus infections, including a 450 mg regimen, and to identify sources of pharmacokinetic variability. METHOD: Plasma samples were collected at 2, 4, 8 and 12 weeks post-transplant and at 4, 6, 8 and 12 months post-transplant in subjects prescribed longer courses. Ganciclovir was measured by liquid chromatography/ultraviolet detection. Non-linear mixed effects modelling was used to analyze the concentration-time data and evaluate demographic and transplant-related covariates. RESULTS: A two compartment model with first order absorption best described the data. Given the range of body sizes, clearance and volume of distribution terms were scaled using standard weight-based allometric exponents. Creatinine clearance was included on apparent oral clearance. Final estimates in a standard 70 kg individual for apparent oral clearance, central volume of distribution, intercompartmental clearance and peripheral volume of distribution were 14.5 l h(-1) , 87.5 l, 4.80 l h(-1) and 42.6 l, respectively. The median terminal half-life for kidney, liver and lung transplant recipients was 9.4, 9.5 and 8.2 h, respectively. Median exposure (i.e. AUC(0,∞) in subjects taking valganciclovir 900 mg or 450 mg once daily was 57.4 and 34.3 µg ml(-1) h, respectively. CONCLUSION: Allometric scaling allowed simultaneous analysis of data from children and adults. Ganciclovir pharmacokinetics were similar among kidney, liver and lung transplant recipients. Ganciclovir exposure after valganciclovir 450 mg once daily may be suboptimal in some individuals and requires evaluation along with virologic outcomes data.