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.

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.

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.

Influence of Renal Function on Phosphoramide Mustard Exposure: A Nonlinear Mixed-Effects Analysis

Phosphoramide mustard (PM) is the final cytotoxic metabolite formed from the parent compound cyclophosphamide through a complex metabolic pathway, primarily through hepatic metabolism. Little is known about the effect of renal elimination on the disposition of PM. We evaluated the effect of renal function on PM exposure after single doses of cyclophosphamide in 85 patients undergoing allogeneic hematopoietic cell transplantation using nonlinear mixed-effects modeling. Mixed linear and nonlinear elimination pathways were required to adequately describe the disposition of PM. Creatinine clearance (CrCL) was incorporated as a covariate associated with first-order elimination, representing renal clearance (Cl_R ) of PM. For a 70-kg patient, Cl_R was 14.9 L/h, Volume of distribution was 525 L, maximum rate was 81.2 mg/h, and the concentration to achieve 50% of maximum rate was 0.51 mg/L. We conducted simulations to explore the impact of CrCL as a measure of renal function and observed that when CrCL decreases from 120 to 40 mL/min, PM area under the plasma concentration-time curve (AUC) from time 0 to 8 hours and AUC increases by 9.2% and 80.9% on average after a single dose, respectively. Our data suggest that renal function has limited influence on PM exposure during the first 8 hours after dosing but has a large impact on the total exposure. Dose adjustment of cyclophosphamide may not be necessary in hematopoietic cell transplant recipients with moderate to severe kidney dysfunction to attain targeted exposures based on AUC from time 0 to 8 hours. However, dose reduction may be necessary if demonstrated at some future time that total AUC is a better surrogate for safety or toxicity.

Findings from a pilot study of buprenorphine population pharmacokinetics: A potential effect of HIV on buprenorphine bioavailability

BACKGROUND

Buprenorphine is widely used in the treatment of opioid use disorder (OUD). There are few pharmacokinetic models of buprenorphine across diverse populations. Population pharmacokinetics (POPPK) allows for covariates to be included in pharmacokinetic studies, thereby opening the potential to evaluate the effect of comorbidities, medications, and other factors on buprenorphine pharmacokinetics. This pilot study used POPPK to explore buprenorphine pharmacokinetics in patients with and without HIV receiving buprenorphine for OUD.

METHODS

Plasma buprenorphine levels were measured in 54 patients receiving buprenorphine for OUD just prior to and 2-5 h following regular buprenorphine dosing. A linear one-compartment POPPK model with first-order estimation was used to evaluate buprenorphine clearance (CL/F) and volume of distribution (V/F). Covariates included weight and HIV status.

RESULTS

All HIV+ patients reported complete past-month adherence to taking antiretroviral therapy that included either efavirenz or nevirapine. Buprenorphine CL/F was 76% higher in HIV+ patients (n = 17) than HIV- patients (n = 37). Buprenorphine V/F was 41% higher in the HIV+ patients.

CONCLUSIONS

POPPK can be used to model buprenorphine pharmacokinetics in a real-world clinical population. While interactions between ART and buprenorphine alter buprenorphine CL/F, we also found alteration in V/F. Proportionate changes in CL/F and V/F might indicate a primary effect on bioavailability (F) rather than two separate effects. These findings indicate reduced buprenorphine bioavailability in patients with HIV.

Rapid antibody responses to Epstein-Barr virus correlate with reduced severity of primary infection

BACKGROUND

We investigated Epstein-Barr virus (EBV) antibody kinetics in university freshmen who developed laboratory-documented primary EBV infection during prospective studies and correlated these kinetics with disease severity.

METHODS

EBV-naïve participants had blood collected periodically and sera tested for EBV-specific antibodies with line blot and enzyme immunoassays. The line blot assay contained EBNA-1, p18, p23, BZLF-1, p138, and p54 antigens; the enzyme immunoassay contained viral capsid antigen and EBNA-1. Severity of illness (SOI) was graded 0 (asymptomatic) to 6 (bedridden). Participants with maximum SOI scores 0-2 were compared with those whose maximum SOI scores were 3-6. Time to first antibody response was analyzed using the semi-parametric COX model.

RESULTS

A total of 201 sera from 38 college students collected before, during, and after primary EBV infection were tested. Earlier antibody responses correlated with milder symptoms. This was most pronounced for late-developing antibodies. The median time to development of p18 IgG was significantly earlier among low-SOI participants (64 days) than high-SOI patients (119 days; P = 0.0003).). Participants with mild disease developed EBNA-1 antibodies sooner than participants with more severe disease (125 days versus >270 days; P = 0.017). Participants with mild disease also showed more rapid loss of antibodies against IgG EA p138 and p54 ≥12 weeks post-infection (P = 0.012 and P = 0.026, respectively).

CONCLUSIONS

These data suggest that rapid antibody responses to EBV correlate with reduced severity of primary EBV infection.

Pharmacogenetics-based population pharmacokinetic analysis for dose optimization of ritonavir-boosted atazanavir in Thai adult HIV-infected patients

BACKGROUND

This population pharmacokinetic-pharmacogenetic study aimed to investigate the optimal dose of RTV-boosted ATV (ATV/RTV) for Thai adult HIV-infected patients.

METHODS

A total of 1460 concentrations of ATV and RTV from 544 patients receiving an ATV/RTV-based regimen were analyzed. The CYP3A5 6986 A > G, ABCB1 3435 C > T, ABCB1 2677 G > T, SLCO1B1 521 T > C, and NR1I2 63396 C > T were genotyped. A population pharmacokinetic model was performed using a nonlinear mixed-effect model (NONMEM®). Monte Carlo simulations were conducted to compare the percentages of patients achieving the therapeutic range of ATV through concentrations (C_trough).

RESULTS

The apparent oral clearance of ATV (CL/F_ATV) without RTV was 7.69 L/h with interindividual variability (IIV) of 28.7%. Patients with CYP3A5 6986 GG had a 7.1% lower CL/F_ATV than those with AA or AG genotype. The CL/F_ATV decreased by 10.8% for females compared with males. Simulation results showed higher percentages (~70%) of patient receiving doses of 200/100 or 200/50 mg achieved the target ATV C_trough, while more patients (~40%) receiving a standard dose (300/100 mg) had ATV C_trough above this target.

CONCLUSIONS

Both CYP3A5 6986 A > G and female decreased CL/F_ATV in Thai HIV-infected patients. Simulations supported that the reduced dose of ATV/RTV was sufficient to achieve the target concentration for Thai population.

Effect of HIV, antiretrovirals, and genetics on methadone pharmacokinetics: Results from the methadone antiretroviral pharmacokinetics study

BACKGROUND

Methadone treatment of opioid use disorder in HIV-infected individuals is complicated by drug-drug interactions. Genetic and other cofactors further contribute to interindividual variability in methadone pharmacokinetics. We used population pharmacokinetics to estimate the effect of drug-drug interactions, genetics, and other cofactors on methadone pharmacokinetics in a methadone maintained population in Vietnam.

METHODS

Plasma R- and S-methadone levels were measured in 309 methadone maintained individuals just before and 2-5 h following methadone dosing. A linear one-compartment population pharmacokinetic model with first-order conditional estimation with interaction was used to evaluate methadone clearance (CL/F) and volume of distribution (V/F). The influence of covariates on parameter estimates was evaluated using stepwise covariate modeling. Covariates included HIV status, antiretroviral use (efavirenz or nevirapine), weight, BMI, age, methadone dose, and 8 single nucleotide polymorphisms in across the CYP2B6, ABCB1, and NR1I3 genes.

RESULTS

Taking either efavirenz or nevirapine increased R-methadone CL/F 220%. Nevirapine and efavirenz increased S-methadone CL/F by 404% and 273%, respectively. Variants in NR1I3 increased R- and S-methadone CL/F by approximately 20% only in patients taking efavirenz. Different alleles in ABCB1 rs2032582 either increased or decreased R-methadone CL/F by 10%. The CYP 2B6*4 variant decreased S-methadone CL/F by 18%. HIV-infection increased R- and S-methadone CL/F and V/F by 24%-39%.

CONCLUSIONS

The HIV antiretrovirals nevirapine and efavirenz significantly increase methadone clearance. Variants inNR1I3 increased the effect of efavirenz on methadone clearance. Other variants affecting methadone CL/F were also confirmed. To our knowledge, this is the first report of HIV itself affecting methadone pharmacokinetics.

Population Pharmacokinetic Analysis of N-Acetylcysteine in Pediatric Patients With Inherited Metabolic Disorders Undergoing Hematopoietic Stem Cell Transplant

N-acetylcysteine (NAC) has been used in patients with cerebral adrenoleukodystrophy as an antioxidant agent in association with hematopoietic stem cell transplant (HSCT). However, an understanding of the pharmacokinetic characteristics of intravenous NAC dosing in these patients is limited. If and how NAC pharmacokinetics change following the transplant is unknown. Toward that end, a total of 260 blood samples obtained from 18 pediatric patients with inherited metabolic disorders who underwent HSCT were included in a population pharmacokinetic analysis using nonlinear mixed-effects modeling. NAC clearance (CL) and volume of distribution (V) were explored on 3 occasions: -7, +7, and +21 days relative to transplant. Additionally, the effect of transplant procedure on NAC disposition was explored by accounting for between-occasion variability. The covariate OCC was modeled as a fixed-effect parameter on CL and/or V1. A 2-compartment model adequately described the pharmacokinetics of total NAC. Weight-based allometric scaling on pharmacokinetic parameters was assumed using standard coefficients. Estimates for CL, central (V1), and peripheral volume (V2), and intercompartment clearance were 14.7 L/h, 23.2 L, 17.1 L, 3.99 L/h, respectively, for a 70-kg person. The data only supported between-subject variability in CL (12%) and V1 (41%). Residual variability was estimated to be 16%. HSCT did not change CL and V1 significantly, and analysis across occasions did not reveal any trends. Pharmacokinetic parameter estimates were in general comparable to those reported previously in different populations. These results suggest that dosing of NAC does not need to be altered following HSCT.

A pharmacokinetic simulation study to assess the performance of a sparse blood sampling approach to quantify early drug exposure

Estimating early exposure of drugs used for the treatment of emergent conditions is challenging because blood sampling to measure concentrations is difficult. The objective of this work was to evaluate predictive performance of two early concentrations and prior pharmacokinetic (PK) information for estimating early exposure. The performance of a modeling approach was compared with a noncompartmental analysis (NCA). A simulation study was performed using literature‐based models for phenytoin (PHT), levetiracetam (LEV), and valproic acid (VPA). These models were used to simulate rich concentration‐time profiles from 0 to 2 h. Profiles without residual unexplained variability (RUV) were used to obtain the true partial area under the curve (pAUC) until 2 h after the start of drug infusion. From the profiles with the RUV, two concentrations per patient were randomly selected. These concentrations were analyzed under a population model to obtain individual population PK (PopPK) pAUCs. The NCA pAUCs were calculated using a linear trapezoidal rule. Percent prediction errors (PPEs) for the PopPK pAUCs and NCA pAUCs were calculated. A PPE within ±20% of the true value was considered a success and the number of successes was obtained for 100 simulated datasets. For PHT, LEV, and VPA, respectively, the median value of the success statistics obtained using the PopPK approach of 81%, 92%, and 88% were significantly higher than the 72%, 80%, and 67% using the NCA approach (p < 0.05; Mann–Whitney U test). This study provides a means by which early exposure can be estimated with good precision from two concentrations and a PopPK approach. It can be applied to other settings in which early exposures are of interest.

Early Exposure of Fosphenytoin, Levetiracetam, and Valproic Acid After High-Dose Intravenous Administration in Young Children With Benzodiazepine-Refractory Status Epilepticus

Fosphenytoin (FOS) and its active form, phenytoin (PHT), levetiracetam (LEV), and valproic acid (VPA) are commonly used second-line treatments of status epilepticus. However, limited information is available regarding LEV and VPA concentrations following high intravenous doses, particularly in young children. The Established Status Epilepticus Treatment Trial, a blinded, comparative effectiveness study of FOS, LEV, and VPA for benzodiazepine-refractory status epilepticus provided an opportunity to investigate early drug concentrations. Patients aged ≥2 years who continued to seizure despite receiving adequate doses of benzodiazepines were randomly assigned to FOS, LEV, or VPA infused over 10 minutes. A sparse blood-sampling approach was used, with up to 2 samples collected per patient within 2 hours following drug administration. The objective of this work was to report early drug exposure of PHT, LEV, and VPA and plasma protein binding of PHT and VPA. Twenty-seven children with median (interquartile range) age of 4 (2.5-6.5) years were enrolled. The total plasma concentrations ranged from 69 to 151.3 μg/mL for LEV, 11.3 to 26.7 μg/mL for PHT and 126 to 223 μg/mL for VPA. Free fraction ranged from 4% to 19% for PHT and 17% to 51% for VPA. This is the first report in young children of LEV concentrations with convulsive status epilepticus as well as VPA concentrations after a 40 mg/kg dose. Several challenges limited patient enrollment and blood sampling. Additional studies with a larger sample size are required to evaluate the exposure-response relationships in this emergent condition.

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 I_max 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 IC₅₀ 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 IC₅₀ 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.

Manipulation of Hydrocortisone Tablets Leads to Iatrogenic Cushing Syndrome in a 6-Year-Old Girl With CAH

Currently there are no commercially available hydrocortisone formulations for the treatment of children with congenital adrenal hyperplasia (CAH) that allow for smaller doses (0.1-1.25 mg) and incremental adjustments needed to control excess androgen production and avoid the negative effects of overtreatment. This lack of availability has led physicians to recommend dividing hydrocortisone 5-mg tablets into 4 to 6 pieces, compounding capsules or hydrocortisone suspension, or crushing 5- or 10-mg tablets in 5 or 10 mL of water. We report a case of iatrogenic Cushing syndrome in a 6-year 11-month-old girl with salt-wasting CAH treated with hydrocortisone tablets that were administered after crushing and dispersing into water to obtain the prescribed dose. She presented with poor growth, increasing body mass index (BMI), excess downy hair, round facies, and gastric ulcers. Her hydrocortisone dose was 8.1 mg/m²/day. Results for all adrenal steroid concentrations were undetectable at 8 am, 12 hours after her last dose. The year prior to presentation her parents began dissolving 10 mg of hydrocortisone in 10 mL of water and using this preparation over the course of 24 hours, which coincided with rapid increase of BMI. We switched her to a pharmacy-compounded alcohol-free hydrocortisone suspension with total daily doses ranging from 6.5 to 8.2 mg/m²/day, which resulted in resolution of her cushingoid features, a decrease in BMI, and catch-up growth. Our case highlights that manipulation of hydrocortisone tablets by parents can result in great variability in dosing and the need for commercially available pediatric formulations allowing for smaller dosing required in young children.

Population Pharmacokinetics of Sertraline in Healthy Subjects: a Model-Based Meta-analysis

Sertraline pharmacokinetics is poorly understood and highly variable due to large between-subject variability with inconsistent reports for oral bioavailability. The study objective was to characterize sertraline pharmacokinetics by developing and validating a sertraline population pharmacokinetic (PK) model in healthy subjects using published clinical PK data. We carried a systematic literature search in PubMed in October 2015 and identified 27 pharmacokinetic studies of sertraline conducted in healthy adult subjects and reported in the English language. Sixty mean plasma concentration-time profiles made of 748 plasma concentrations following IV, single, and multiple oral doses ranging from 5 to 400 mg were extracted and analyzed for dose proportionality by a log-linear model and fitted to a 2-compartment pharmacokinetic model in NONMEM using a model-based meta-analysis (MBMA) approach. After a single oral dose, sertraline C_max and AUC_∞ increased with dose proportionally between 50 and 200 mg, and bioavailability increased nonlinearly with dose from 5 to 50 mg and plateaued afterwards while T_max and t_1/2 did not change with dose. Following multiple oral doses, C_max and AUC_∞ increased proportionally with dose across the entire dose range (5-200 mg) while bioavailability, T_max, and t_1/2 remained constant with dose. Sertraline absorption was time-dependent and best described by a sigmoidal E_max function of time after dose. Study findings indicate that sertraline PK is linear in healthy adult subjects at doses ≥ 50 mg, and exposures were nonlinear only after single oral doses < 50 mg likely due to reduced bioavailability.

Comparative Evaluation of Median Versus Youden Index Dichotomization Methods: Exposure-Response Analysis of Mycophenolic Acid and Acyl-Glucuronide Metabolite

BACKGROUND AND OBJECTIVES

Dichotomization of pharmacokinetic exposure measures in exposure-response relationship studies provides results that are interpretable in clinical care. Several methods exist in the literature on how to define the cut-off values needed for the dichotomization process. Commonly, the sample median is utilized to define the dichotomizing value; however, statistical methods based on the exposure metric and its association with the outcome are argued to result in a more proper definition of the optimal cut-point. The Youden index is a recommended statistical method to define the cut-off value. The current analysis objective is to compare the dichotomization results based on the Youden index versus median methods.

METHODS

Utilizing mycophenolic acid (MPA) exposure data and its related acute rejection and leukopenia outcome variables, the current study compared the MPA exposure-response relationship outcomes when MPA exposure is dichotomized via the Youden index versus median methods. Univariate logistic models were utilized to quantify the relationships between MPA exposure, including total MPA, unbound MPA, and the acyl-glucuronide metabolite of MPA, and the probabilities of acute rejection and leukopenia.

RESULTS

The overall trend of the results of the logistic models demonstrated a general similarity in the inferred exposure-response associations when considering either the Youden index-based or the median-based dichotomization methods.

CONCLUSION

The results demonstrated in this analysis suggest that both the Youden index and the median methods provide similar conclusions when dichotomization of a continuous variable is considered. However, confirmation of these conclusions comes from future powered studies that include a larger number of subjects.

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.

Precision medicine in adult and pediatric obesity: a clinical perspective

It remains largely unknown as to why some individuals experience substantial weight loss with obesity interventions, while others receiving these same interventions do not. Person-specific characteristics likely play a significant role in this heterogeneity in treatment response. The practice of precision medicine accounts for an individual's genes, environment, and lifestyle when deciding upon treatment type and intensity in order to optimize benefit and minimize risk. In this review, we first discuss biopsychosocial determinants of obesity, as understanding the complexity of this disease is necessary for appreciating how difficult it is to develop individualized treatment plans. Next, we present literature on person-specific characteristics associated with, and predictive of, weight loss response to various obesity treatments including lifestyle modification, pharmacotherapy, metabolic and bariatric surgery, and medical devices. Finally, we discuss important gaps in our understanding of the causes of obesity in relation to the suboptimal treatment outcomes in certain patients, and offer solutions that may lead to the development of more effective and targeted obesity therapies.

Iron supplements in nursing home patients associated with reduced carbamazepine absorption

Persons in nursing homes receive a number of medications that may interfere with the pharmacokinetics of carbamazepine (CBZ). The aim of our study was to determine factors that may affect the pharmacokinetics of CBZ in elderly nursing home patients.

METHODS

CBZ concentration data collected from 60 nursing homes across the US were evaluated. Inclusion criteria included residency in a nursing home for at least 2 months, age 65 years or older, a stable dosing regimen of CBZ for at least 4 weeks (considered steady state), available CBZ concentration, and complete information regarding all co-medications. Using a nonlinear mixed-effects model, the data were adequately described by a one-compartment model with first-order absorption and elimination. Goodness-of-fit plots, plausibility of parameter estimates, visual predictive check and nonparametric bootstrap were used to evaluate the models.

MAIN FINDINGS

The final data set consisted of 345 CBZ concentrations from 99 subjects (38 males, 61 females). The population estimate of apparent clearance (CL/F) for a 70-kg person was 3.69 L/hr (RSE 6.9%). Residents were receiving either immediate (93.9%) or extended release (6.1%) formulation of CBZ and the Ka of each formulation was fixed to literature values. Age, sex, and co-medications had no effect on CL/F and apparent volume of distribution. Iron supplementation, which was taken by 16% of the residents, resulted in a 33% decrease in bioavailability (p < 0.001). No other medications were found to have an effect.

CONCLUSIONS

Results from this pharmacokinetic study indicate that use of iron supplementation is associated with a reduction in absorption of CBZ and may need to be considered when dosing CBZ in patients taking iron supplementation.

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.

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.

Mycophenolic Acid and Its Metabolites in Kidney Transplant Recipients: A Semimechanistic Enterohepatic Circulation Model to Improve Estimating Exposure

Mycophenolic acid (MPA) is an approved immunosuppressive agent widely prescribed to prevent rejection after kidney transplantation. Wide between-subject variability (BSV) in MPA exposure exists which in part may be due to variability in enterohepatic recirculation (EHC). Several modeling strategies were developed to evaluate EHC as part of MPA pharmacokinetics, however mechanistic representation of EHC is limited. These models have not provided a satisfactory representation of the physiology of EHC in their modeling assumptions. The aim of this study was i) to develop an integrated model of MPA (total and unbound) and its metabolites (MPAG and acyl-MPAG) in kidney recipients, where this model provides a more physiological representation of EHC process, and ii) to evaluate the effect of donor and recipient clinical covariates and genotypes on MPA disposition. A five-compartment model with first-order input into an unbound MPA compartment connected to the MPAG, acyl-MPAG, and gallbladder compartment best fit the data. To represent the EHC process, the model was built based on the physiological concepts related to the hepatobiliary system and the gallbladder filling and emptying processes. The effect of cyclosporine versus tacrolimus on clearance of unbound MPA was included in the base model. Covariate analysis showed creatinine clearance to be significant on oral clearance of unbound MPA. The hepatic nuclear factor 1 alpha (HNF1A) genetic single nucleotide polymorphism (SNP) (rs2393791) in the recipient significantly affected the fraction of enterohepatically-circulated drug. Oral clearance of MPAG was affected by recipient IMPDH1 SNP (rs2288553), diabetes at the time of transplant, and donor sex.

Pharmacokinetics and Saturable Absorption of Gabapentin in Nursing Home Elderly Patients

Pharmacokinetic data of gabapentin (GBP) in community-dwelling elderly patients show a significant effect of advanced age on GBP pharmacokinetics due to altered renal function. However, there are no data in elderly nursing home (NH) patients to evaluate gabapentin absorption and elimination. Our objective was to characterize the pharmacokinetics of GBP in elderly nursing home patients maintained on GBP therapy. This was a prospective pharmacokinetic study in elderly nursing home patients (≥60 years) receiving GBP for the management of chronic pain or epilepsy from seven nursing homes. Pharmacokinetic parameters were estimated by nonlinear mixed-effects modeling. A one-compartment model described the data and clearance (CL) was associated with estimated glomerular filtration rate (eGFR) (p < 0.0001). The GBP CL in elderly nursing home patients was 2.93 L/h. After adjusting for the effect of GFR, GBP CL was not affected by age, sex, body weight, or comorbidity scores. No significant effects of body size measures, age, and sex were detected on volume of distribution. Dose-dependent bioavailability of GBP was demonstrated, and the saturable absorption profile was described by a nonlinear hyperbolic function. Prediction-corrected visual predictive check (pc-VPC) suggests adequate fixed- and random-effects models that successfully simulated the mean trend and variability in gabapentin concentration-time profiles. In this analysis, the parameters of the hyperbolic nonlinearity appear to be similar between elderly and younger adults.

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 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.

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.

Pharmacokinetic-Pharmacodynamic Modeling of Intravenous and Oral Topiramate and Its Effect on the Symbol-Digit Modalities Test in Adult Healthy Volunteers

A sequential pharmacokinetic-pharmacodynamic (PK-PD) modeling approach was used to quantify the effects of a single dose of topiramate (100 or 200 mg) on working memory, attention, and psychomotor speed as measured by the Symbol-Digit Modalities Test (SDMT). Established on data pooled from 3 randomized, crossover studies in healthy subjects (19-55 years of age), using both oral and a novel stable-labeled intravenous (IV) formulation of topiramate, an inhibitory Emax model was found to characterize the topiramate concentration-SDMT score relationship well. At the EC50 of 2.85 μg/mL, this topiramate plasma concentration value was estimated to be associated with a 25.5% reduction of SDMT score relative to baseline. Age was an important determinant of the baseline SDMT score, with an estimated decrease of 1.13% in baseline SDMT score with every year of age. Moreover, this approach enabled the quantification of the practice effect observed with repeated administration of the neuropsychological test over shorter testing intervals than have previously been reported in the literature. The finding of a significant effect following a single dose of topiramate in the range widely used to treat migraine and epilepsy needs to be evaluated in a broader patient population undergoing chronic treatment, as the narrow range of resultant concentrations limits the generalizability of the findings.

Genotype-guided tacrolimus dosing in African-American kidney transplant recipients

Tacrolimus is dependent on CYP3A5 enzyme for metabolism. Expression of the CYP3A5 enzyme is controlled by several alleles including CYP3A5*1, CYP3A5*3, CYP3A5*6 and CYP3A5*7. African Americans (AAs) have on average higher tacrolimus dose requirements than Caucasians; however, some have requirements similar to Caucasians. Studies in AAs have primarily evaluated the CYP3A5*3 variant; however, there are other common nonfunctional variants in AAs (CYP3A5*6 and CYP3A5*7) that do not occur in Caucasians. These variants are associated with lower dose requirements and may explain why some AAs are metabolically similar to Caucasians. We created a tacrolimus clearance model in 354 AAs using a development and validation cohort. Time after transplant, steroid and antiviral use, age and CYP3A5*1, *3, *6 and *7 alleles were significant toward clearance. This study is the first to develop an AA-specific genotype-guided tacrolimus dosing model to personalize therapy.

Steady-state pharmacokinetics and bioavailability of immediate-release and extended-release formulations of lamotrigine in elderly epilepsy patients: Use of stable isotope methodology

A classic 2-period crossover bioavailability study was conducted to evaluate the relative and absolute bioavailability of immediate-release (IR) and extended-release (XR) lamotrigine formulations under steady-state conditions in elderly patients with epilepsy. On treatment days, each subject's morning dose (IR or XR lamotrigine) was replaced with an intravenous 50-mg dose of stable-labeled lamotrigine. Lamotrigine concentrations were measured at 13 points between 0 and 96 hours. XR and IR lamotrigine formulations were similar with respect to steady-state area under the concentration-time curve from 0 to 24 hours (AUC0-24 h ss), average concentration (Cavg, ss), and trough concentration (Cτ, ss). A 33% lower fluctuation in concentrations with XR was observed relative to IR lamotrigine. The time to peak concentration (Tmax, ss) was delayed for XR lamotrigine (3.0 vs 1.3 hours) with lower peak concentration (15% lower). The absolute bioavailability for IR and XR formulations was 73% and 92%, respectively. The formulations were bioequivalent with respect to AUC0-24 h ss, Cτ, ss, and Cavg, ss indicating that it may be possible to switch directly from IR to XR lamotrigine without changes in the total daily dose.

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.

Ethnic and genetic factors in methadone pharmacokinetics: a population pharmacokinetic study

BACKGROUND

Treatment of opiate use disorders with methadone is complicated by wide interindividual variability in pharmacokinetics. To identify potentially contributing covariates in methadone pharmacokinetics, we used population pharmacokinetic modeling to estimate clearance (CL/F) and volume of distribution (V/F) for each methadone enantiomer in an ethnically diverse methadone maintained population.

METHODS

Plasma levels of the opiate-active R-methadone and opiate-inactive S-methadone were measured in 206 methadone maintained subjects approximately two and twenty-three hours after a daily oral dose of rac-methadone. A linear one-compartment population pharmacokinetic model with first-order conditional estimation with interaction (FOCE-I) was used to evaluate methadone CL/F and V/F. The influence of covariates on parameter estimates was evaluated using stepwise covariate modeling. Covariates included ethnicity, gender, weight, BMI, age, methadone dose, and 21 single nucleotide polymorphisms in genes implicated in methadone pharmacokinetics.

RESULTS

In the final model, for each enantiomer, Hmong ethnicity reduced CL/F by approximately 30% and the rs2032582 (ABCB1 2677G>T/A) GG genotype was associated with a 20% reduction in CL/F. The presence of the rs3745274 minor allele (CYP2B6 515G>T) reduced CL/F by up to 20% for S-methadone only. A smaller effect of age was noted on CL/F for R-methadone.

CONCLUSION

This is the first report showing the influence of the rs2032582 and rs3745274 variants on methadone pharmacokinetics rather than simply dose requirements or plasma levels. Population pharmacokinetics is a valuable method for identifying the influences on methadone pharmacokinetic variability.

MODEL-BASED LAMOTRIGINE CLEARANCE CHANGES DURING PREGNANCY: CLINICAL IMPLICATION

Objective

The objective of the study was to characterize changes in the oral clearance (CL/F) of lamotrigine (LTG) over the course of pregnancy and the postpartum period through a model-based approach incorporating clinical characteristics that may influence CL/F, in support of developing clinical management guidelines.

Methods

Women receiving LTG therapy who were pregnant or planning pregnancy were enrolled. Maternal blood samples were collected at each visit. A pharmacokinetic analysis was performed using a population-based, nonlinear, mixed-effects model.

Results

A total of 600 LTG concentrations from 60 women (64 pregnancies) were included. The baseline LTG CL/F was 2.16 L/h with a between-subject variability of 40.6%. The influence of pregnancy on CL/F was described by gestational week. Two subpopulations of women emerged based on the rate of increase in LTG CL/F during pregnancy. The gestational age-associated increase in CL/F displayed a 10-fold higher rate in 77% of the women (0.118 L/h per week) compared to 23% (0.0115 L/h per week). The between-subject variability in these slopes was 43.0%. The increased CL/F at delivery declined to baseline values with a half-life of 0.55 weeks.

Interpretation

The majority of women had a substantial increase in CL/F from 2.16 to 6.88 L/h by the end of pregnancy, whereas 23% of women had a minimal increase. An increase in CL/F may correspond to decreases in LTG blood concentrations necessitating the need for more frequent dosage adjustments and closer monitoring in some pregnant women with epilepsy. Postpartum doses should be tapered to preconception dose ranges within 3 weeks of delivery.

SLC28A3 genotype and gemcitabine rate of infusion affect dFdCTP metabolite disposition in patients with solid tumours

Background:

Gemcitabine is used for the treatment of several solid tumours and exhibits high inter-individual pharmacokinetic variability. In this study, we explore possible predictive covariates on drug and metabolite disposition.

Methods:

Forty patients were enrolled. Gemcitabine and dFdU concentrations in the plasma and dFdCTP concentrations in peripheral blood mononuclear cell were measured to 72 h post infusion, and pharmacokinetic parameters were estimated by nonlinear mixed-effects modelling. Patient-specific covariates were tested in model development.

Results:

The pharmacokinetics of gemcitabine was best described by a two-compartment model with body surface area, age and NT5C2 genotype as significant covariates. The pharmacokinetics of dFdU and dFdCTP were adequately described by three-compartment models. Creatinine clearance and cytidine deaminase genotype were significant covariates for dFdU pharmacokinetics. Rate of infusion of <25 mg m⁻² min⁻¹ and the presence of homozygous major allele for SLC28A3 (CC genotype) were each associated with an almost two-fold increase in the formation clearance of dFdCTP.

Conclusion:

Prolonged dFdCTP systemic exposures (⩾72 h) were commonly observed. Infusion rate <25 mg m⁻² min⁻¹ and carriers for SLC28A3 variant were each associated with about two-fold higher dFdCTP formation clearance. The impacts of these covariates on treatment-related toxicity in more selected patient populations (that is, first-line treatment, single disease state and so on) are not yet clear.

Validation of tacrolimus equation to predict troughs using genetic and clinical factors

AIM

Tacrolimus is an immunosuppressant used in transplantation. This article reports the validation of the authors' recently developed genetics-based tacrolimus equation that predicts troughs.

METHODS

Validation was performed in an independent cohort of 795 kidney transplant recipients receiving tacrolimus. The performance of the equation to predict initial troughs was assessed by calculating the bias and precision of the equation. For all troughs in the first 6 months post-transplant, a comparison was made between the troughs predicted using the equation versus those predicted using a basic apparent clearance model with no covariates.

RESULTS

For initial troughs, the equation had a low bias (0.2 ng/ml) and high precision (1.8 ng/ml). For all troughs, the equation predicted troughs significantly better than the basic apparent clearance model.

CONCLUSION

The tacrolimus equation had good bias and precision in predicting initial troughs and performed better than a basic apparent clearance model for all the troughs.

Population pharmacokinetics of unbound and total drug concentrations following intravenously administered carbamazepine in elderly and younger adult patients with epilepsy

The objective of the study was to investigate the pharmacokinetics (PK) of unbound and total plasma carbamazepine (CBZ) concentrations following simultaneous administration of intravenous and oral formulations. We tested the hypothesis that age-related alterations in physiology and patient characteristics influence CBZ disposition and protein binding. Patients (n = 113) on maintenance therapy received a 100 mg dose of a novel, intravenous, stable-labeled (SL) CBZ formulation as partial replacement of their morning CBZ dose. A two-compartment model described unbound and total SL-CBZ data. The stable-labeled intravenous dosing methodology enabled the estimation of the CBZ clearance (CL) and volumes of distribution. The CL of CBZ was dependent on race through the model equation unbound CL (L/hour) = 11.2 × (1.30)(Race); where Race = 1 for Caucasian, 0 for African American. Total body weight explained 57% and 70% of the interindividual variability in the central and peripheral volumes of distribution, respectively. Age, sex, smoking, plasma albumin, and alpha 1-acid glycoprotein concentrations had no effect on CL, binding or volumes of distribution. The model was evaluated via bootstrap and predictive check. Results may support race specific dosing for CBZ where an average African-American individual would receive 70% of the standard dose prescribed for the Caucasian person.

Population pharmacokinetics of rifampicin in Mexican patients with tuberculosis

WHAT IS KNOWN AND OBJECTIVE

Rifampicin (RIF) shows wide variability in its pharmacokinetics. The purpose of this study was to develop and validate a population pharmacokinetic model to characterize the inter- and intra-individual variability in pharmacokinetic parameters of RIF in Mexican patients.

METHODS

Ninety-four patients receiving antituberculosis therapy participated in this prospective study. Plasma concentration-time data were described using a one-compartment model with lag time, absorption and first-order elimination. The potential influence of demographic and clinical characteristics of the patients, and the pharmaceutical formulation (A, B, C and D) on the pharmacokinetics parameters, was evaluated by non-linear mixed-effect modelling (nonmem). Seventy-seven additional patients participated in the validation of the model.

RESULTS AND DISCUSSION

The final population pharmacokinetic model obtained was as follows: apparent clearance CL/F = 8·17 L/h (1·40 as high for males), apparent distribution volume V(d)/F = 50·1 L (1·29 as high for males), absorption rate constant K(aA) = 0·391/h, K(aB,C,D) = 2·70/h, relative bioavailability F(A) = 0·468, F(B,C,D) = 1, lag time in the absorption phase T(lag) = 0·264 h. The final model improved the precision on the parameter estimates (CL/F, V(d) /F and K(a) by 31·9%, 16·7% and 92·9%, respectively). The residual variability was 27·3%.

WHAT IS NEW AND CONCLUSION

Gender was associated with changes in CL/F and V(d) /F whereas the pharmaceutical formulation was associated with changes in F and altered the K(a) . The validation data set showed that the model could be used in clinical practice for Bayesian dose adjustment of RIF in TB patients.