Parent-metabolite Population Pharmacokinetic Model Research Articles

Impact of CYP2D6*2, CYP2D6*35, rs5758550, and related haplotypes on risperidone clearance in vivo.

The CYP2D6 gene exhibits significant polymorphism, contributing to variability in responses to drugs metabolized by CYP2D6. While CYP2D6*2 and CYP2D6*35 are presently designated as alleles encoding normal metabolism, this classification is based on moderate level evidence. Additionally, the role of the formerly called "enhancer" single nucleotide polymorphism (SNP) rs5758550 is unclear. In this study, the impacts of CYP2D6*2, CYP2D6*35 and rs5758550 on CYP2D6 activity were investigated using risperidone clearance as CYP2D6 activity marker. A joint parent-metabolite population pharmacokinetic model was used to describe 1,565 serum concentration measurements of risperidone and 9-hydroxyrisperidone in 512 subjects. Risperidone population clearance was modeled as the sum of a CYP2D6-independent clearance term and the partial clearances contributed from each individually expressed CYP2D6 allele or haplotype. In addition to the well-characterized CYP2D6 alleles (*3-*6, *9, *10 and *41), *2, *35 and two haplotypes assigned as CYP2D6*2-rs5758550G and CYP2D6*2-rs5758550A were evaluated. Each evaluated CYP2D6 allele was associated with significantly lower risperidone clearance than the reference normal function allele CYP2D6*1 (p < 0.001). Further, rs5758550 differentiated the effect of CYP2D6*2 (p = 0.005). The haplotype-specific clearances for CYP2D6*2-rs5758550A, CYP2D6*2-rs5758550G and CYP2D6*35 were estimated to 30%, 66% and 57%, respectively, relative to the clearance for CYP2D6*1. Notably, rs5758550 is in high linkage disequilibrium (R2 > 0.85) with at least 24 other SNPs and cannot be assigned as a functional SNP. CYP2D6*2 and CYP2D6*35 encode reduced risperidone clearance, and the extent of reduction for CYP2D6*2 is differentiated by rs5758550. Genotyping of these haplotypes might improve the precision of genotype-guided prediction of CYP2D6-mediated clearance.

Population Pharmacokinetic Modelling for Estimation of Remifentanil Metabolic-Ratio Using Non-steady-State Concentrations under Rapidly Adaptive Dosing.

To predict steady-state metabolite-to-drug concentration ratio (metabolic ratio) for analgesic drug remifentanil, using sparse non-steady-state data from patients with normal or impaired renal function during individualised, highly variable and rapidly adaptive intravenous infusion. A three-compartment joint parent-metabolite population pharmacokinetic model was developed using concentrations of remifentanil and its metabolite remifentanil acid from two trials. Renal function was included as an important mechanistic covariate. To address the large covariate effect and highly individualised and rapidly adaptive dosing, standardised visual predictive check was conducted on the observations and individualised visual predictive check was conducted on metabolic ratio estimates. The model was used to simulate metabolic ratio distribution in patients with various renal functions. The model, including its covariate structure, adequately described the data. The predictive checks allowed informative model evaluation. The predictedmedian (10th - 90th percentile) of remifentanil metabolic ratio was 12.5 (2.4-58.2) for patients with normal or mildly impaired renal function, or 54.3 (12.8-218.4) for patients with moderately or severely impaired renal function. The methodologies applied here allowed robust estimation of steady-state parameters using non-steady-state sparse data under highly variable adaptive dosing.

Relationship between exposure to treosulfan and its monoepoxytransformer – An insight from population pharmacokinetic study in pediatric patients before hematopoietic stem cell transplantation

Treosulfan (TREO), a structural analog of busulfan, is currently studied as a myeloablative agent in conditioning regimens before hematopoietic stem cell transplantation in pediatric patients. High exposure to TREO (>1650 mg∗h/mL) might be related to early toxicity, especially skin toxicity and mucositis. The aim of the present study was to investigate a potential relationship between exposure to TREO and its monoepoxytransformer (S,S-EBDM), as well as variability of the pharmacokinetics of these entities by means of a population pharmacokinetic approach with a non-linear mixed-effects analysis.The study included data from 14 children with malignant and non-malignant diseases treated with TREO in daily doses 10–14 g/m2. The parent-metabolite population pharmacokinetic model was developed in NONMEM 7.3 software. Upon the constructed model, an extensive simulation was performed to assess the correlation between exposure to TREO and S,S-EBDM.It was found that TREO and S,S-EBDM pharmacokinetics was best described with 2-compartmental and 1-compartmental linear models, respectively. The vast majority (>65%) of TREO was transformed to S,S-EBDM. Overall, a considerable interpatient variability of pharmacokinetic parameters was observed, especially the clearance of S,S-EBDM. A weak correlation was found between the exposure to TREO and S,S-EBDM (r = 0.1681, p < 0.0001). Also, patients with an exposure to TREO above 1650 mg∗h/mL were most likely to have also a high exposure to S,S-EBDM (35.38 μM∗h vs. 43.14 μM∗h, p < 0.0001).In summary, a parent-metabolite population pharmacokinetic model for TREO and S,S-EBDM was developed for the first time. It was shown that there is a weak correlation between exposure to TREO and S,S-EBDM. Therefore therapeutic drug monitoring of not only prodrug but also its active epoxide might be needed.

Polymorphic Expression of UGT1A9 is Associated with Variable Acetaminophen Glucuronidation in Neonates: A Population Pharmacokinetic and Pharmacogenetic Study

Acetaminophen (paracetamol, APAP) is widely used as an analgesic and antipyretic drug in children and neonates. A number of enzymes contribute to the metabolism of acetaminophen, and genetic factors might be important to explain variability in acetaminophen metabolism among individuals. The current investigation utilized a previously published parent-metabolite population pharmacokinetic model describing acetaminophen glucuronidation, sulfation, and oxidation to examine the potential role of genetic variability on the relevant metabolic pathways. Neonates were administered 30-min intravenous infusions of acetaminophen 15mg/kg every 12h (<28weeks' gestational age [GA]) or every 8h (≥28weeks GA) for 48h. A total of 18 sequence variations (SVs) in UDP-glucuronosyltransferase (UGT), sulfotransferase (SULT), and cytochrome P450 (CYP) genes from 33 neonates (aged 1-26days) were examined in a stepwise manner for an effect on the metabolic formation clearance of acetaminophen by glucuronidation (UGT), sulfation (SULT), and oxidation (CYP). The stepwise covariate modeling procedure was performed using NONMEM® version 7.3. Incorporation of genotype as a covariate for one SV located in the UGT1A9 gene promoter region (rs3832043, - 118>insT, T9>T10) significantly improved model fit (likelihood ratio test, p<0.001) and reduced between-subject variability in glucuronide formation clearance. Individuals with the UGT1A9 T10 polymorphism, indicating insertion of an additional thymidine nucleotide, had a 42% reduction in clearance to APAP-glucuronide as compared to their wild-type counterparts. This study shows a pharmacogenetic effect of an SV in the UGT1A9 promoter region on the metabolism of acetaminophen in neonates.

Population pharmacokinetics of oxcarbazepine and its monohydroxy derivative in epileptic children.

Oxcarbazepine is an antiepileptic drug with an activity mostly due to its monohydroxy derivative metabolite (MHD). A parent-metabolite population pharmacokinetic model in children was developed to evaluate the consistency between the recommended paediatric doses and the reference range for trough concentration (Ctrough ) of MHD (3-35mgl-1 ). A total of 279 plasma samples were obtained from 31 epileptic children (age 2-12 years) after a single dose of oxcarbazepine. Concentration-time data were analysed with Monolix 4.3.2. The probability to obtain Ctrough between 3-35mgl-1 was determined by Monte Carlo simulations for doses ranging from 10 to 90mgkg-1 day-1 . A parent-metabolite model with two compartments for oxcarbazepine and one compartment for MHD best described the data. Typical values for oxcarbazepine clearance, central and peripheral distribution volume and distribution clearance were 140lh-1 70kg-1 , 337l70kg-1 , 60.7l and 62.5lh-1 , respectively. Typical values for MHD clearance and distribution volume were 4.11lh-1 70kg-1 and 54.8l70kg-1 respectively. Clearances and distribution volumes of oxcarbazepine and MHD were related to body weight via empirical allometric models. Enzyme-inducing antiepileptic drugs (EIAEDs) increased MHD clearance by 29.3%. Fifty-kg children without EIAEDs may need 20-30mgkg-1 day-1 instead of the recommended target maintenance dose (30-45mgkg-1 day-1 ) to obtain Ctrough within the reference range. By contrast, 10-kg children with EIAEDs would need 90mgkg-1 day-1 instead of the maximum recommended dose of 60mgkg-1 day-1 . This population pharmacokinetic model of oxcarbazepine supports current dose recommendations, except for 10-kg children with concomitant EIAEDs and 50-kg children without EIAEDs.

Influence of genetic co-factors on the population pharmacokinetic model for clopidogrel and its active thiol metabolite

PurposeA high interindividual variability is observed in the pharmacokinetics of clopidogrel, a widely used antiplatelet drug. In the present study, a joint parent-metabolite population pharmacokinetic model was developed to adequately describe observed concentrations of clopidogrel and its active thiol metabolite (H4).MethodsThe study included 63 patients undergoing elective coronarography or percutaneous coronary intervention. The population pharmacokinetic model was developed in the NONMEM 7.3 software, and first-order conditional estimation method with interaction was applied. Also, the influence of covariates was evaluated (age, weight, body mass index (BMI), obesity defined as BMI ≥ 30 kg/m2, sex, diabetes mellitus, co-administration of PPI or statins, presence of CYP2C19*2, CYP2C19*17, CYP3A4*1G alleles, and ABCB1 3435 TT genotype).ResultsIt was found that the only significant covariate was the presence of CYP2C19*2 allele, which had an impact on lower conversion of clopidogrel to H4. As a result, predicted area under the time-concentration curve values was lower in carriers of this allele, with median 5.94 ng h/ml (interquartile range 3.92–12.51 [ng∙h/ml]) vs. 12.70 ng h/ml in non-carriers (interquartile range, 7.00–19.39 [ng∙h/ml]), respectively (p = 0.004).ConclusionsDeveloped model predicts that the only significant covariate influencing the observed concentrations and therefore the exposure to the active H4 metabolite is the presence of CYP2C19*2 allele.

Influences of CYP2D6*10 polymorphisms on the pharmacokinetics of iloperidone and its metabolites in Chinese patients with schizophrenia: a population pharmacokinetic analysis.

Iloperidone is an atypical antipsychotic drug that is mainly metabolized by CYP2D6, CYP3A4, and cytosolic enzymes. Previous studies show that extensive and poor metabolizers of CYP2D6 exhibit different plasma concentrations of iloperidone and its metabolites. The aim of this study was to develop a parent-metabolite population pharmacokinetic (PPK) model to quantify the effects of CYP2D6*10 allele on the pharmacokinetics of iloperidone and its metabolites in Chinese schizophrenia patients. Seventy Chinese schizophrenia patients were enrolled, from whom limited blood samples were collected on d 15 (0 h) and d 28 (0, 4 and 12 h after drug administration). The plasma concentrations of iloperidone and its metabolites M1 (P-88) and M2 (P-95) were simultaneously detected using a validated HPLC-MS assay. CYP2D6*10 (rs1065852) genotyping was performed. A PPK model was developed based on data from the patients using the NONMEM software (version 7.2). A one-compartment model with first-order absorption and elimination was used to describe the pharmacokinetic data related to iloperidone and its metabolites. Patients with the CYP2D6*10 T/T genotype had significantly higher concentrations of iloperidone and M1, and lower concentrations of M2 than the patients with C/C or C/T genotypes. The CYP2D6*10 genotype affected the elimination constants for transformation of iloperidone to the metabolites M1 (K23) and M2 (K24). The K23 value of the patients with T/T genotype was 1.34-fold as great as that of the patients with C/C or C/T genotype. The K24 value of the patients with C/T and T/T genotypes was 0.693- and 0.492-fold, respectively, as low as that of the patients with C/C genotype. CYP2D6*10 mutations affect the pharmacokinetics of iloperidone and its metabolites in Chinese schizophrenia patients, suggesting that the clinical doses of iloperidone for patients with CYP2D6*10 mutations need to be optimized.

Evaluation of the effect of residual error covariance for a parent-metabolite population pharmacokinetic model

Background The residual error of parent and metabolite concentrations in pharmacokinetic (PK) samples can be correlated due to issues such as mislabeled timing, sample extraction, and misspecification, although it is not typically estimated in population modeling. We investigated the impact of excluding this covariance on population estimates and clinical trial simulations. Methods Multivariate data from a single-dose oral study in Japanese and Caucasian male healthy volunteers (n=25) were described with a parent-metabolite population PK model using nonlinear mixed-effects modeling. Weight effects on clearances and volumes of distribution were included. Covariance between parent and metabolite residual errors was estimated or excluded. A multiple-dose, study was then simulated using both models. Results Both models described the observed data well. Parameter estimates for clearances and volumes were similar between models. The correlation between proportional and additive residual errors was estimated at 61 and 100%, respectively. Multiple-dose simulations of 100 trials displayed that excluding the covariance term influenced predicted covariate effect on exposure by −4 to +17%. Conclusions The exclusion of covariance in residual error may have only minor impacts on resulting population parameter estimates, however the distribution of simulated multiple-dose trial outcomes may become exaggerated, which may influence the predicted exposure in subpopulations. Clinical Pharmacology & Therapeutics (2005) 77, P91–P91; doi: 10.1016/j.clpt.2004.12.242

A population pharmacokinetic model for doxorubicin and doxorubicinol in the presence of a novel MDR modulator, zosuquidar trihydrochloride (LY335979).

To develop a population pharmacokinetic model for doxorubicin and doxorubicinol in the presence of zosuquidar.3HCl, a potent P-glycoprotein inhibitor. The population approach was used (implemented with NONMEM) to analyse doxorubicin-doxorubicinol pharmacokinetic data from 40 patients who had received zosuquidar.3HCl and doxorubicin intravenously (separately in cycle 1 and concomitantly in cycle 2 over 48 h and 0.5 h, respectively). A five-compartment pharmacokinetic model (including three compartments for doxorubicin pharmacokinetics with two pathways for doxorubicinol formation) best described the doxorubicin-doxorubicinol pharmacokinetics in the presence of zosuquidar.3HCl. Doxorubicin clearance (CL), peripheral volume of distribution (V2) and doxorubicinol apparent clearance (CLm/fm) and apparent volume of distribution (Vm/fm) were 62.3 l/h, 2360 l, 143 l/h and 3150 l, respectively, in the absence or presence of low doses of zosuquidar.3HCl (<500 mg). In the presence of high doses of zosuquidar.3HCl (>or=500 mg), these values decreased by 25%, 26%, 48% and 73%, respectively, and doxorubicinol pharmacokinetics were characterized by a delayed t(max) (24 h versus 4 h), which led to the inclusion of the parallel pathways. A decrease in the objective function ( P<0.005) was observed when the impact of zosuquidar.3HCl was accounted for. This integrated parent-metabolite population pharmacokinetic model accurately characterized the increase in doxorubicin and doxorubicinol exposure (1.33- and 2-fold, respectively) in the presence of zosuquidar.3HCl (>or=500 mg) and provided insights into the pharmacokinetic interaction, which may be useful in designing future clinical trials.