Search

Purpose: The purpose of this study was to develop a population pharmacokinetic (PK) model for 3-AP, to evaluate the effect of ABCB1 polymorphisms on the pharmacokinetic profile of 3-AP, and to assess the relationship between 3AP disposition and patient covariates.Methods: A total of 40 patients with advanced cancer from two phase 1 studies were included in the population PK model building. Patients received 3-AP 25-105 mg/m(2) IV on day 1. 3-AP plasma and erythrocyte levels were sampled at 10 timepoints over a 24-h period and measured by a validated HPLC method. Data were analyzed by a nonlinear mixed-effects modeling approach using the NONMEM system.Results: 3-AP pharmacokinetics were described as a 3-compartment model with first-order elimination, with one compartment representing the plasma and another representing erythrocyte concentrations. Gender was associated with volume of distribution, in which women had a lower V2. The number of cycles administered was associated with clearance; those with decreased clearance were more likely to receive less than 2 cycles before going off study.Conclusion: This study suggests that monitoring 3-AP plasma concentrations in the first cycle and dose adjustment in those with decreased clearance may be helpful in decreasing toxicity associated with the 3-AP. [ABSTRACT FROM AUTHOR]

Purpose: Gemcitabine, a pyrimidine nucleoside, is approved for the treatment of non-small cell lung cancer, pancreatic carcinoma, and breast cancer. Chemotherapy regimens are determined experimentally with static tissue culture systems, animal models, and in Phase I clinical trials. The aim of this study was to assess for gemcitabine-induced cell death following infusion of drug under clinically-relevant conditions of infusion rate and drug exposure in an in vitro bioreactor system.Methods: To estimate an appropriate harvest time for cells from the bioreactor after drug treatment, we estimated the temporal relationship between gemcitabine treatment for 1 h and cell death at a later time point with monolayer growth assays (i.e., static culture). Afterward, 5.3 mg gemcitabine was infused over 0.5 h in the bioreactor, followed by mono-exponential decay, simulating patient concentration-time profiles (n = 4). Controls were run with drug-free media (n = 4). Cells were harvested from the bioreactor at a later time point and assessed for cell death by flow cytometry.Results: According to monolayer growth assay results, cytotoxicity became more apparent with increasing time. The E Max for cells 48 h after treatment was 50% and after 144 h, 93% (P = 0.022; t test), while flow cytometry showed complete DNA degradation by 120 h. Gemcitabine was infused in the bioreactor. The gemcitabine area under the concentration-time curve (AUC) was 56.4 microM h and the maximum concentration was 87.5 +/- 2.65 microM. Flow cytometry results were as follows: the G1 fraction decreased from 65.1 +/- 4.91 to 28.6 +/- 12% (P = 0.005) and subG1 increased from 14.1 +/- 5.28 to 42.6 +/- 9.78% (P = 0.004) relative to control. An increase in apoptotic cells was observed by TUNEL assay.Conclusions: The in vitro bioreactor system will be expanded to test additional cell lines, and will serve as a useful model system for assessing the role of drug pharmacokinetics in delivery of optimized anticancer treatment. [ABSTRACT FROM AUTHOR]

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 E max 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., (© 2023 The Authors. British Journal of Clinical Pharmacology published by John Wiley & Sons Ltd on behalf of British Pharmacological Society.)

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., (© 2023. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

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., (Copyright © 2023 by The American Society for Pharmacology and Experimental Therapeutics.)

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., (© 2022 Genentech Inc. The Journal of Clinical Pharmacology published by Wiley Periodicals LLC on behalf of American College of Clinical Pharmacology.)

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., Competing Interests: Declaration of Competing Interest The authors report no declarations of interest., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Cyclophosphamide (CY) is an alkylating agent widely used in the field of oncology and hematopoietic cell transplantation (HCT). It is recommended to use an adjusted body weight with an adjustment factor of 0.25 (ABW25) for dosing of CY in obese patients undergoing HCT. However, evidence based on the pharmacokinetics (PK) of CY to support this recommendation is lacking. We aimed to identify a dosing strategy of CY that achieves equivalent exposures among obese and nonobese patients. The present study is a secondary analysis of a previously conducted observational PK study of phosphoramide mustard (PM), the final cytotoxic metabolite of CY. Data were collected from 85 adults with hematologic malignancy who received a single infusion of CY 50 mg/kg, fludarabine, ± anti-thymocyte globulin, and a single fraction of total body irradiation as HCT conditioning therapy. A previously developed population PK model in these patients was used for simulations. Using individualized PK parameters from that analysis, simulations were performed to assess cumulative exposures of PM (i.e., area-under-the-curve [AUC]) resulting from 8 different dosing strategies according to various measures of body size: (1) "mg/kg" by total body weight (TBW); (2) "mg/kg" by ideal body weight (IBW); (3) "mg/kg" by fat free mass; (4) "mg/m 2 " by body surface area (BSA); (5) "mg/kg" by TBW combined with ABW25 (TBW-ABW25); (6) "mg/kg" by IBW combined with ABW25 (IBW-ABW25); (7) "mg/kg" by TBW combined with ABW by adjustment factor of 0.50 (TBW-ABW50); and (8) "mg" by fixed-dose. We defined equivalent exposure as the effect of obesity on PM AUC within ±20% from the PM AUC in the nonobese group, where obesity is defined based on TBW/IBW ratio (i.e., nonobese, <1.2; mildly obese, 1.2-1.5; and moderately/severely obese, >1.5). Primary and secondary outcomes were PM AUC 0-8hours and PM AUC 0-infinity , respectively. In the 85 patients, with the median age of 63 years (range 21-75), 46% were classified as mildly and 25% were moderately/severely obese based on the TBW/IBW ratio. Negative correlations (i.e., higher the extent of obesity, lower the PM AUC) were shown when dosing simulations were based on IBW, TBW-ABW25, and fixed dosing (P < .05). Positive correlations were shown when dosing was simulated by TBW (P < .05). None of the 8 dosing strategies attained equivalent PM AUC 0-8hours between patients with versus without obesity, whereas dosing by BSA and TBW-ABW50 attained equivalent PM AUC 0-infinity (P < .05). Our study predicted that the recommended ABW25 dose adjustment may result in lower exposure of CY therapy in obese patients than in nonobese. A CY dosing strategy that would result in similar PM concentrations between obese and nonobese was not identified for early exposure (i.e., PM AUC 0-8hours ). The data suggest though that CY dosing based on "mg/m 2 " by BSA or "mg/kg" by TBW-ABW50 would result in similar total exposure (i.e., PM AUC 0-infinity ) and may minimize exposure differences in obese and nonobese patients., (Copyright © 2022 The American Society for Transplantation and Cellular Therapy. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: Declaration of Competing Interest None of the authors has a financial conflict of interest or potential conflict of interest to declare., (Copyright © 2022 Elsevier B.V. All rights reserved.)

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., (Copyright © 2020 American Society for Pharmacology and Experimental Therapeutics.)

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 (CL f ) of 6- and 7-OH-S-warfarin decrease: 73.2% and 74.8%) but less inducible by rifampin (% CL f of 6- and 7-OH-S-warfarin increase: 85% and 75%), compared with non-CYP2C9 mediated elimination pathways (% CL f of 10-OH-S-warfarin and CL R of S-warfarin decrease in the presence of fluconazole: 65.0% and 15.3%; % CL f 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., (Copyright © 2020 American Society for Pharmacology and Experimental Therapeutics.)

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.

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., (© 2021, The American College of Clinical Pharmacology.)

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., (Copyright © 2021 Elsevier B.V. All rights reserved.)

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., (© 2021 The Authors. Clinical and Translational Science published by Wiley Periodicals LLC on behalf of the American Society for Clinical Pharmacology and Therapeutics.)

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., (© 2020, The American College of Clinical Pharmacology.)

A specific model for drug absorption is necessarily assumed in pharmacokinetic (PK) analyses following extravascular dosing. Unfortunately, an inappropriate absorption model may force other model parameters to be poorly estimated. An added complexity arises in population PK analyses when different individuals appear to have different absorption patterns. The aim of this study is to demonstrate that a deep neural network (DNN) can be used to prescreen data and assign an individualized absorption model consistent with either a first-order, Erlang, or split-peak process. Ten thousand profiles were simulated for each of the three aforementioned shapes and used for training the DNN algorithm with a 30% hold-out validation set. During the training phase, a 99.7% accuracy was attained, with 99.4% accuracy during in the validation process. In testing the algorithm classification performance with external patient data, a 93.7% accuracy was reached. This algorithm was developed to prescreen individual data and assign a particular absorption model prior to a population PK analysis. We envision it being used as an efficient prescreening tool in other situations that involve a model component that appears to be variable across subjects. It has the potential to reduce the time needed to perform a manual visual assignment and eliminate inter-assessor variability and bias in assigning a sub-model.

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 50 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 50 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., (© 2020 The Authors. British Journal of Clinical Pharmacology published by John Wiley & Sons Ltd on behalf of British Pharmacological Society.)

Significant in vitro and in vivo evidence supports the potential use of sertraline as an anticancer and antimicrobial agent. Yet, it is unknown whether effective sertraline concentrations are clinically achieved at therapeutic doses. The study objectives were to develop a physiologically-based pharmacokinetic (PBPK) model of sertraline and estimate the probability of achieving effective concentrations in various human tissues. A generic PBPK model consisting of perfusion-limited compartments representing the body organs linked together by blood flows and incorporated with clearance, tissue distribution, and absorption models was implemented in R using the mrgsolve package. Sertraline clearance and volume of distribution were first optimized from i.v. plasma concentration data then absorption and bioavailability were optimized from oral data. Predicted unbound sertraline concentrations at steady-state in human tissues did not reach concentrations determined in vitro, indicating therapeutic doses of sertraline are unlikely to produce concentrations required for anticancer and antimicrobial activities in humans., (© 2020 The Authors Clinical Pharmacology & Therapeutics © 2020 American Society for Clinical Pharmacology and Therapeutics.)

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.

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

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.

Background and Objectives: Strategies for modeling the enterohepatic circulation (EHC) process reported in the literature vary; however, gallbladder-based models currently provide the best physiological representation of the process. Regardless, the addition of a gallbladder to the model does not fully depict the physiology of EHC. A more physiological gallbladder-based EHC model is needed. This model should take into account a physiological representation of the bile secretion, gallbladder filling and emptying, the duration of gallbladder emptying, and irregular mealtimes. Considering all of these factors, the objectives of the present analysis were to propose a gallbladder-based EHC model and then to use that model to perform sensitivity analyses evaluating the effect of the extent of EHC on the pharmacokinetic profile and noncompartmental analysis (NCA) calculations., Methods: A gallbladder-based model that describes the EHC process was developed and used to perform determinant simulations assuming various degrees of EHC. Next, these simulations were compared to evaluate the effect of the EHC on the pharmacokinetic profiles of orally administered drugs. The influence of the EHC process on the NCA calculations was determined while assuming two sampling schemes that differed in the times at which sampling was performed in relation to meal times., Results: The presence of EHC results in nonlinearity in the system and changes the pharmacokinetic profile, affecting the maximum concentration (C max ), time to C max (T max ), and half-life estimates. Comparison of the results obtained using the two sampling schemes for a drug undergoing various degrees of EHC demonstrated a significant influence of the selected sampling times on the NCA estimations. Bias in the NCA calculations was also dependent on the sampling times used., Conclusion: Caution should be taken when designing clinical studies for drugs that undergo EHC. It may be essential to consider the timing of meals when planning pharmacokinetic studies and defining sampling times. The period over which samples are taken needs to be extended as compared to that traditionally used with other drugs. Future studies that attempt to identify the best sampling strategies in the presence of EHC are needed.

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., Competing Interests: Conflict of interest statement: J.R.R. received research support in the form of drug/placebo from Boehringer Ingelheim. C.K.F. received research support from Novo Nordisk. S.D.S. received grant funding from Astra Zeneca Pharmaceuticals. A.S.K. received research support (drug/placebo) from Astra Zeneca Pharmaceuticals and served as a consultant for Novo Nordisk, WW, and Vivus Pharmaceuticals but did not accept personal or professional income for these activities. The other authors have no disclosures.

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 K a 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., (Copyright © 2018. Published by Elsevier B.V.)

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., (Wiley Periodicals, Inc. © 2018 International League Against Epilepsy.)

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., (© 2018, The American College of Clinical Pharmacology.)

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 [ fC min ] of <10 mg/liter)- and efficacy (free, unbound plasma concentration-to-MIC ratio [ fC max /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 fC max /MIC in 80% of pediatric patients weighing 8 to 70 kg with a 97.5% probability, and almost all patients were predicted to have fC min 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., (Copyright © 2018 American Society for Microbiology.)

Multiple factors can impact warfarin therapy, including genetic variations in the drug-metabolizing enzyme cytochrome P450 2C9 (CYP2C9). Compared with individuals with the wild-type allele, CYP2C9*1, carriers of the common *3 variant have significantly impaired CYP2C9 metabolism. Genetic variations in CYP2C9, the primary enzyme governing the metabolic clearance of the more potent S-enantiomer of the racemic anticoagulant warfarin, may impact warfarin-drug interactions. To establish a baseline for such studies, plasma and urine concentrations of R- and S-warfarin and 10 warfarin metabolites were monitored for up to 360 hours following a 10-mg warfarin dose in healthy subjects with 4 different CYP2C9 genotypes: CYP2C9*1/*1 (n = 8), CYP2C9*1/*3 (n = 9), CYP2C9*2/*3 (n = 3), and CYP2C9*3/*3 (n = 4). Plasma clearance of S-warfarin, but not R-warfarin, decreased multiexponentially and in a CYP2C9 gene-dependent manner: 56%, 70%, and 75% for CYP2C9*1/*3, CYP2C9*2/*3, and CYP2C9*3/*3 genotypes, respectively, compared with CYP2C9*1/*1, resulting in pronounced differences in the S:R ratio that identified warfarin-sensitive genotypes. CYP2C9 was the primary P450 enzyme contributing to S-warfarin metabolism and a minor contributor to R-warfarin metabolism. In the presence of a defective CYP2C9 allele, switching of warfarin metabolism to other oxidative pathways and P450 enzymes for the metabolic elimination of S-warfarin was not observed. The 10-hydroxywarfarin metabolites, whose detailed pharmacokinetics are reported for the first time, exhibited a prolonged half-life with no evidence of renal excretion and displayed elimination rate-limited kinetics. Understanding the impact of CYP2C9 genetics on warfarin pharmacokinetics lays the foundation for future genotype-dependent warfarin-drug interaction studies., (© 2016, The American College of Clinical Pharmacology.)

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.

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.

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., (© 2016 The British Pharmacological Society.)

We performed an intensive prospective study designed to obtain as much data as possible on the incubation and early illness periods of primary Epstein-Barr virus (EBV) infection. Undergraduate students who lacked EBV antibody and oral EBV DNA (EBV-naive) were seen every 2 weeks during their freshman year. Clinical and behavioral data, oral washes and venous blood were obtained. EBV antibodies were quantified by enzyme immunoassay and viral loads by PCR. During a median 8 months of observation, 14/85 subjects experienced primary EBV infections (24 cases/100 person-years). The only significant risk factor for acquisition of EBV infection was deep kissing (P=0.02). Eleven subjects had infectious mononucleosis with a median duration of 21 days. Two subjects were hospitalized. Infections were initially identified in 12 subjects by finding EBV DNA in oral cells before onset of symptoms and in 2 subjects by symptom reporting. EBV DNA and viral capsid antigen (VCA) IgM and gp350 IgG antibodies were present in the blood before onset of illness. To provide a more robust evaluation of primary EBV infection in undergraduate university students, we combined data on risk factors and antibody responses from this and an earlier study that used the exact same clinical and laboratory methods. The observation that the only significant risk factor for acquisition of EBV infection was deep kissing was confirmed. Most importantly, higher amounts of gp350 antibody correlated significantly with a lower severity of infectious mononucleosis (P<0.0001), which strengthens the rationale for a gp350-based prophylactic EBV vaccine.

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., (© 2016 The British Pharmacological Society.)

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., (© 2015, The American College of Clinical Pharmacology.)

Characterizing Epstein-Barr virus (EBV) dynamics in asymptomatic immunocompetent persons provides a baseline for defining quantitative thresholds associated with EBV disease. Studying latent membrane protein (LMP)-1 sequence variation over time could establish the rates of reactivation and superinfection, and also trace transmission. Twelve asymptomatic adult subjects were evaluated prospectively nine times over 6 months. EBV serum antibodies were measured by enzyme immunoassay. EBV DNA in oral and whole-blood samples was quantitated by real-time (TaqMan) PCR and analyzed for LMP-1 sequence variability. All 11 antibody positive subjects had EBV DNA detected in their oral compartment at least once during the 6-month study. The quantities ranged from 1.70 to 4.91 log10 copies EBV per ml of oral cell pellet. One subject was continuously viremic for 79 days. Overall, EBV DNA was detected in 63 (24%) of 260 samples from 11 antibody-positive subjects and in 0/27 samples from an antibody-negative subject. The quantities in positive samples ranged from 1.7 to 4.9 log10 copies EBV per ml. EBV LMP-1 gene sequence variations in subjects were constant over time regardless of the compartment sampled. Subjects 18-30 years old had EBV DNA detected more frequently than subjects >30 years old (38/108 positive samples versus 25/152; P<0.001). In conclusion, EBV DNA shedding is common in asymptomatic adults. The younger adults shed more frequently, which may reflect a shorter time from their primary EBV infection to sampling. The LMP-1 sequence analysis method employed here could be used to trace person-to-person transmission because patterns remained almost identical over time.

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.

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., (© 2015, The American College of Clinical Pharmacology.)

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., (© 2014 The British Pharmacological Society.)

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., (Copyright © 2014 Elsevier Ireland Ltd. All rights reserved.)

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., (© 2014 The British Pharmacological Society.)

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 of LTG CL/F during pregnancy. The gestational age-associated increase in CL/F displayed a ten-fold higher rate in 77% of the women (0.118 L/h/week) compared to 23% (0.0115 L/h/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.

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(-2) min(-1) 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(-2) min(-1) 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.