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

Population pharmacokinetics of mycophenolic acid and dose optimisation in adult Chinese kidney transplant recipients

1. This study aimed to establish a population pharmacokinetic (PPK) model of mycophenolic acid (MPA), quantify the effect of clinical factors and pharmacogenomics of MPA, and optimise the dosage for adult kidney transplant recipients.2. One-hundred and four adult renal transplant patients were enrolled. The PPK model was established using the Phoenix® NMLE software and the stepwise methods were filtered for significant covariates. Monte Carlo simulations were performed to optimise the dosage regimen.3. A two-compartment model with first-order absorption and elimination (including lag time) provided a more accurate description of MPA pharmacokinetics. Serum albumin (ALB) significantly affected the central apparent clearance (CL/F), whereas post-transplant time and creatinine clearance were associated with a central apparent volume of distribution (V/F). The estimated population values obtained by the final model were 17.5 L/h and 93.97 L for CL/F and V/F, respectively. Simulation results revealed that larger mycophenolate mofetil doses are required as the ALB concentration decreases. This study established a PPK model of MPA and validated it using various methods. ALB significantly affected CL/F and recommended optimal dose strategies were given based on the final model. These results provide a reference for the personalised therapy of MPA for kidney transplant patients.

Significant Effects of Renal Function on Mycophenolic Acid Total Clearance in Pediatric Kidney Transplant Recipients with Population Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVES

Mycophenolic acid (MPA) is an immunosuppressant commonly prescribed in pediatric kidney transplantation to prevent graft rejection. Large variabilities in MPA plasma exposures have been observed in this population, which could result in severe adverse effects. The majority of the MPA pharmacokinetic data have been reported in adult populations, whereas information in pediatric patients is still very limited. The objective of this study was to establish a novel, nonlinear mixed-effects model for MPA and investigate the clinical variables affecting MPA population pharmacokinetics in pediatric kidney transplant recipients.

METHODS

Data were collected retrospectively from pediatric kidney transplant patients (≤ 18 years when MPA concentrations were initially collected; on oral administration of mycophenolate mofetil) in Calgary, Alberta, Canada. Nonlinear mixed-effect modeling was conducted using stochastic approximation expectation-maximization in Monolix 2021R2 (Lixoft SAS, France) to determine population pharmacokinetic estimates, interindividual variabilities, and interoccasional variabilities. Covariate models were constructed using the Model Proposal function in Monolix in conjunction with a systematic stepwise inclusion/elimination protocol. The best model was selected based on objective function values, relative standard errors, goodness-of-fit plots, prediction-corrected visual predictive checks, and numerical predictive checks.

RESULTS

A total of 50 pediatric kidney transplant patients (25 female) with 219 MPA plasma concentration-time profiles were included. The average age (± standard deviation) and posttransplant time for the sample population were 12.8 ± 4.8 years and 762 ± 1160 days, respectively. The majority of study subjects (i.e., > 85% based on all occasions) were co-administered tacrolimus. A two-compartment, first-order absorption with lag time and linear elimination structural model with lognormal distributed proportional residual errors best described the MPA concentration-time data. The absorption rate constant (2.52 h-1 or 0.042 min-1), lag time (0.166 h or 9.96 min), volumes of distributions of the central (22.8 L) and peripheral (216 L) compartments, and intercompartment clearance (17.6 L h-1 or 0.293 L min-1) were consistent with literature values; whereas total MPA clearance (0.72 L h-1 or 0.012 L min-1) was relatively reduced, likely due to the general lack of cyclosporine interactions and the stabilized graft functions from significantly longer posttransplant time in our sample population. Of the clinical variables tested, only estimated glomerular filtration rate (eGFR) was identified a significant covariate affecting total MPA clearance with a positive, exponential relationship. The final population pharmacokinetic model was successfully evaluated/validated using a variety of complementary methods.

CONCLUSION

We have successfully constructed and validated a novel population pharmacokinetic model of MPA in pediatric kidney transplant patients. A positive, nonlinear relationship between eGFR and total MPA clearance identified in our model is likely attributed to multiple concurrent mechanisms, which warrant further systematic investigations.

Immune Monitoring of Mycophenolate Mofetil Activity in Healthy Volunteers Using Ex Vivo T Cell Function Assays

Mycophenolate mofetil (MMF) is part of the standard immunosuppressive treatment after transplantation and usually given as "one-dose-fits-all" together with a calcineurin inhibitor (CNI). Although drug concentrations are frequently monitored, there is still a group of patients who experience side effects related to excessive or insufficient immune suppression. We therefore aimed to identify biomarkers that reflect the overall immune status of the patient and might support individualized dosing. We previously studied immune biomarkers for CNIs and aimed to investigate whether these are also suitable to monitor MMF activity. Healthy volunteers received a single dose of MMF or placebo, after which IMPDH enzymatic activity, T cell proliferation, and cytokine production were measured and compared to MPA (MMF's active metabolite) concentration in three different matrices (plasma, peripheral blood mononuclear cells, and T cells). MPA concentrations in T cells exceeded those in PBMCs, but all intracellular concentrations correlated strongly with plasma concentrations. At clinically relevant MPA concentrations, IL-2 and IFN-γ production was mildly suppressed, while MPA T cell proliferation was strongly inhibited. Based on these data, it is expected that monitoring of T cell proliferation in MMF-treated transplantation patients may be a valid strategy to avoid excessive immune suppression.

Influence of Calcineurin Inhibitors and Genetic Polymorphism of Transporters on Enterohepatic Circulation and Exposure of Mycophenolic Acid in Chinese Adult Renal Allograft Recipients

There is significant enterohepatic circulation (EHC) during the disposition of mycophenolic acid (MPA). The aim of this study was to elucidate factors influencing the EHC of MPA in Chinese adult renal allograft recipients. After 2 weeks of therapy with mycophenolate mofetil or enteric-coated mycophenolate sodium, blood samples were collected from 125 patients at 0 to 12 hours post-administration and MPA concentrations were determined. The influence of calcineurin inhibitors (CNIs) and genetic polymorphisms on MPA exposure and EHC was studied. The Shapley additive explanations method was used to estimate the impact of various factors on the area under the plasma drug concentration-time curve (AUC_0-12h ) for MPA. An extreme gradient boosting (XGboost) machine learning-based model was established to predict AUC_0-12h . Results showed that the dose-normalized AUC_6-12h (dn-AUC_6-12h ) of MPA was significantly lower in patients co-administered with cyclosporine (CsA) than in patients co-administered with tacrolimus (TAC) (P < .05). For patients co-administered with TAC, patients with ABCC2 C-24T CC or SLCO1B1 T521C TT genotypes had significantly higher values of dn-AUC_6-12h (P < .05). Patients with SLCO1B3 334T/699G alleles had significantly lower dn-AUC_6-12h values than homozygotes (P < .05). By introducing body weight, age, and EHC-related factors, including co-administered CNIs and genetic polymorphism of drug transporters, as covariates in the XGboost machine learning model, the prediction performance of AUC_0-12h for MPA in Chinese adult renal allograft recipients can be improved.

Clinical Evidence on the Purported Pharmacokinetic Interactions between Corticosteroids and Mycophenolic Acid

Corticosteroids (steroids) are commonly used concurrently with mycophenolic acid (MPA) as the first-line immunosuppression therapy for the prevention of rejection in solid organ transplantations. Steroids are also commonly administered with MPA in various autoimmune disorders such as systemic lupus erythematosus and idiopathic nephrotic syndrome. Despite various review articles having suggested the presence of pharmacokinetic interactions between MPA and steroids, definitive data have not yet been demonstrated. The aim of this Current Opinion is to critically evaluate the available clinical data and propose the optimal study design for characterising the MPA-steroid pharmacokinetic interactions. The PubMed and Embase databases were searched for relevant clinical articles in English as of September 29, 2022, where a total of 8 papers have been identified as supporting and 22 as non-supporting the purported drug interaction. To objectively evaluate the data, novel assessment criteria to effectively diagnose the interaction based on known MPA pharmacology were formulated, including the availability of independent control groups, prednisolone concentrations, MPA metabolite data, unbound MPA concentrations, and the characterisations of entero-hepatic recirculation and MPA renal clearance. Overall, the majority of the identified corticosteroid data were pertaining to prednisone or prednisolone. Our assessment indicated that no conclusive mechanistic data supporting the interaction are available in the current clinical literature, and further studies are required to quantify the effects/mechanisms of steroid-tapering or withdrawal on MPA pharmacokinetics. This current opinion provides justification for further translational investigations, as this particular drug interaction has the potential to exert significant adverse outcomes in patients prescribed MPA.

Effect of Meal Timings and Meal Content on the AUC0-12h of Mycophenolic Acid: A Simulation Study

Meal timings and content related to gallbladder emptying in the enterohepatic circulation are important for explaining the high variability in mycophenolic acid exposure. The limited sampling strategy (LSS) was established to estimate the area under the plasma concentration-time curve from time 0 to 12 hours (AUC_0-12h ) of mycophenolic acid in therapeutic drug monitoring. The aim of this study is to investigate the effect of meal timings and content on the AUC_0-12h of mycophenolic acid and to assess the influence of meals on LSS. A mycophenolic acid pharmacokinetic model with a mechanism-based enterohepatic circulation process was employed to perform simulations under various assumed meal scenarios. The simulations were compared to evaluate the effect of meal timings and meal content on mycophenolic acid AUC_0-12h . Monte Carlo simulations were performed using the meal parameter with the greatest impact on mycophenolic acid AUC_0-12h as a variable. The corresponding LSS equations were established, and the predictive performance was assessed. Both the meal timings and meal content affected the mycophenolic acid AUC_0-12h , and the postdose fasting period had the greatest impact. The predictive performance of the LSS is sensitive to the postdose fasting period. Therefore, meal timings may improve the estimation of mycophenolic acid AUC_0-12h and the efficacy of therapeutic drug monitoring.

Fecal β-glucuronidase activity differs between hematopoietic cell and kidney transplantation and a possible mechanism for disparate dose requirements

The intestinal microbiota produces β-glucuronidase that plays an essential role in the metabolism of the immunosuppressant mycophenolate mofetil (MMF). This drug is commonly used in organ and hematopoietic cell transplantation (HCT), with variations in dosing across transplant types. We hypothesized that β-glucuronidase activity differs between transplant types, which may account for differences in dosing requirements. We evaluated fecal β-glucuronidase activity in patients receiving MMF post-allogeneic HCT and post-kidney transplant. Kidney transplant patients had significantly greater β-glucuronidase activity (8.48 ± 6.21 nmol/hr/g) than HCT patients (3.50 ± 3.29 nmol/hr/g; P = .001). Microbially mediated β-glucuronidase activity may be a critical determinant in the amount of mycophenolate entering the systemic circulation and an important factor to consider for precision dosing of MMF.

Important lack of difference in tacrolimus and mycophenolic acid pharmacokinetics between Aboriginal and Caucasian kidney transplant recipients

AIM

To examine whether differences in tacrolimus and mycophenolic acid (MPA) pharmacokinetics contribute to the poorer kidney transplant outcomes experienced by Aboriginal Australians.

METHODS

Concentration-time profiles for tacrolimus and MPA were prospectively collected from 43 kidney transplant recipients: 27 Aboriginal and 16 Caucasian. Apparent clearance (CL/F) and distribution volume (V/F) for each individual were derived from concentration-time profiles combined with population pharmacokinetic priors, with subsequent assessment for between-group difference in pharmacokinetics. In addition, population pharmacokinetic models were developed using the prospective dataset supplemented by previously developed structural models for tacrolimus and MPA. The change in NONMEM objective function was used to assess improvement in goodness of model fit.

RESULTS

No differences were found between Aboriginal and Caucasian groups or empirical Bayes estimates, for CL/F or V/F of MPA or tacrolimus. However, a higher prevalence of CYP3A5 expressers (26% compared with 0%) and wider between-subject variability in tacrolimus CL/F (SD = 5.00 compared with 3.25 L/h/70 kg) were observed in the Aboriginal group, though these differences failed to reach statistical significance (p = .07 and p = .08).

CONCLUSION

There were no differences in typical tacrolimus or MPA pharmacokinetics between Aboriginal and Caucasian kidney transplant recipients. This means that Bayesian dosing tools developed to optimise tacrolimus and MPA dosing in Caucasian recipients may be applied to Aboriginal recipients. In turn, this may improve drug exposure and thereby transplant outcomes in this group. Aboriginal recipients appeared to have greater between-subject variability in tacrolimus CL/F and a higher prevalence of CYP3A5 expressers, attributes that have been linked with inferior outcomes.

Sustained suppression of enterohepatic circulation of mycophenolic acid by antimicrobial-associated diarrhea in a kidney transplant recipient with Crohn's disease: A case report

Mycophenolic acid (MPA) undergoes enterohepatic circulation. A kidney transplant patient on mycophenolate mofetil was treated with tazobactam/piperacillin for pyelonephritis, and developed antimicrobial-associated diarrhea. Consequently, the MPA trough level decreased by approximately 90%. Furthermore, it took approximately a month for the MPA level to normalize even after diarrhea had resolved.

Recent lessons learned from population pharmacokinetic studies of mycophenolic acid: physiological, genomic, and drug interactions leading to the prediction of drug effects

INTRODUCTION

Mycophenolic acid (MPA) is a widely used immunosuppressant in transplantation and autoimmune disease. Highly variable pharmacokinetics have been observed with MPA, but the exact mechanisms remain largely unknown.

AREAS COVERED

The current review provided a critical, comprehensive update of recently published population pharmacokinetic/dynamic models of MPA (n=16 papers identified from PubMed and Embase, inclusive from January 2017 to August 2021), with specific emphases on the intrinsic and extrinsic factors influencing the pharmacology of MPA. The significance of the identified covariates, potential mechanisms, and comparisons to historical literature have been provided.

EXPERT OPINION

While select covariates affecting the population pharmacokinetics of MPA are consistently observed and mechanistically supported, some variables have not been regularly reported and/or lacked mechanistic explanation. Very few pharmacodynamic models were available, pointing to the need to extrapolate pharmacokinetic findings. Ideal models of MPA should consist of: i) utilizing optimal sampling points to allow the characterizations of absorption, re-absorption, and elimination phases; ii) characterizing unbound/total MPA, MPA metabolites, plasma/urinary concentrations, and genetic polymorphisms to facilitate mechanistic interpretations; and iii) incorporating actual outcomes and pharmacodynamic data to establish clinical relevance. We anticipate the field will continue to expand in the next 5 to 10 years.

On the population pharmacokinetics and the enterohepatic recirculation of hyzetimibe and its main metabolite in Chinese healthy subjects

BACKGROUNDS

Hyzetimibe (HS-25), a new drug approved for hypercholesterolemia, exhibits obvious enterohepatic recirculation (EHC) after oral administration. Up to now, little is known about the kinetics of HS-25. Therefore, we performed this population pharmacokinetic (PopPK) analysis aiming to describe the PK behavior of HS-25 and its main metabolite (M1), and to identify significant covariates contributing to the variability.

METHODS

The plasma concentration data used for modeling were obtained from an open-label, single-dose, randomized, two-period crossover bioequivalence study. PopPK modeling was performed with NONMEM 7.4.1 using nonlinear mixed effect modeling approach. Goodness of fit plots, bootstrap and visual predictive check were used for model internal validation. Data from another study was used for external validation.

RESULTS

Data from 16 male and 8 female subjects were used in the PopPK analysis. HS-25 and M1 concentrations in the modeling cohort were well described by a one-compartment model incorporating first-pass metabolism and a gallbladder compartment, accounting for the EHC process. The release kinetic of gall was mimicked by a first-order constant plus a switch on/off effect. Body weight was identified as a significant covariate effecting on the clearance and apparent distribution volume of HS-25, as well as k_mg , the transfer rate from metabolite compartment to gallbladder compartment. Internal and external validation demonstrated an acceptable predictive ability of the final model.

CONCLUSIONS

We present the first PopPK model describing HS-25 and M1 concentrations simultaneously, with the EHC process considered. The modeling and simulation results could provide reference for the clinical use of HS-25.

Discovery of a Potential Liver Fibrosis Inhibitor from a Mushroom Endophytic Fungus by Genome Mining of a Silent Biosynthetic Gene Cluster

Liver fibrosis has accounted for liver diseases and overall mortality, but no relevant drug has been developed. Filamentous fungi are important resources of natural products for pharmaceutical development. Calcarisporium arbuscula is a mushroom endophytic fungus, which primarily produces aurovertins. Here, in an aurovertin null-production mutant, one silent gene cluster (mca17) was activated by overexpression of a pathway-specific zinc finger transcriptional regulator, and a tetramic acid-type compound (1, MCA17-1) was identified. Along with detailed structural characterization, its biosynthesis was proposed to be produced from the core PKS-NRPS hybrid enzyme. Moreover, 1 suppressed the activation of LX-2 upon transforming growth factor-β (TGF-β) challenge and had stronger bioactivity than the positive control obeticholic acid (OCA) against liver fibrosis. Our work suggested that this engineered fungus could be a producer of 1 for promising pharmaceutical development, and alternatively, it would be developed as a mushroom ingredient in dietary therapy to prevent liver fibrosis.

Genetic polymorphisms in metabolic enzymes and transporters have no impact on mycophenolic acid pharmacokinetics in adult kidney transplant patients co-treated with tacrolimus: A population analysis

WHAT IS KNOWN AND OBJECTIVE

Mycophenolate mofetil, an ester prodrug of mycophenolic acid (MPA), is widely used to prevent graft rejection after kidney transplantation. The pharmacokinetic (PK) of MPA has been extensively studied, which revealed a high degree of variability. An integrated population PK (PopPK) model of MPA and its main metabolite mycophenolic acid glucuronide (MPAG) was developed using the adult patients who underwent kidney transplant and were administered oral mycophenolate mofetil combined with tacrolimus.

METHODS

In total, 917 MPA and 740 MPAG concentrations in191 adult patients were analysed via nonlinear mixed-effects modelling. The concentration-time data were adequately described using a chain compartment model, including central and peripheral compartments for MPA and a central compartment for MPAG. Stepwise forward inclusion and backward elimination procedures were used to investigate the effects of genetic polymorphisms, including in UGT1A8, UGT1A9, UGT2B7, ABCB1, ABCC2, ABCG2, SLCO1B1, SLCO1B3, and HNF1α.

RESULTS AND DISCUSSION

These genetic polymorphisms in metabolic enzymes and transporters have no obvious impact on the PK of MPA in adult patients who underwent kidney transplant and were co-treated with tacrolimus. The post-transplant time, serum albumin, and creatinine clearance were identified as significant covariates affecting the PK of MPA and MPAG, which should be considered in the clinical use of mycophenolate mofetil.

WHAT IS NEW AND CONCLUSION

We established a PopPK model of MPA and MPAG in Chinese adult patients who underwent kidney transplant and were co-treated with tacrolimus. Genetic polymorphisms in metabolic enzymes and transporters showed no obvious impact on MMF PK. A model-informed dosing strategy was proposed by the established model, and MMF dose adjustment should be based on ALB levels and the post-transplantation time.

Model-informed precision dosing to optimise immunosuppressive therapy in renal transplantation

Immunosuppressive therapy is pivotal for sustained allograft and patient survival after renal transplantation. However, optimally balanced immunosuppressive therapy is challenged by between-patient and within-patient pharmacokinetic (PK) variability. This could warrant the application of personalised dosing strategies to optimise individual patient outcomes. Pharmacometrics, the science that investigates the xenobiotic-biotic interplay using computer-aided mathematical modelling, provides options to describe and quantify this PK variability and enables identification of patient characteristics affecting immunosuppressant PK and treatment outcomes. Here, we review and critically appraise the available pharmacometric model-informed dosing solutions for the typical immunosuppressants in modern renal transplantation, to guide their initial and subsequent dosing.

Personalized Therapy for Mycophenolate: Consensus Report by the International Association of Therapeutic Drug Monitoring and Clinical Toxicology

ABSTRACT

When mycophenolic acid (MPA) was originally marketed for immunosuppressive therapy, fixed doses were recommended by the manufacturer. Awareness of the potential for a more personalized dosing has led to development of methods to estimate MPA area under the curve based on the measurement of drug concentrations in only a few samples. This approach is feasible in the clinical routine and has proven successful in terms of correlation with outcome. However, the search for superior correlates has continued, and numerous studies in search of biomarkers that could better predict the perfect dosage for the individual patient have been published. As it was considered timely for an updated and comprehensive presentation of consensus on the status for personalized treatment with MPA, this report was prepared following an initiative from members of the International Association of Therapeutic Drug Monitoring and Clinical Toxicology (IATDMCT). Topics included are the criteria for analytics, methods to estimate exposure including pharmacometrics, the potential influence of pharmacogenetics, development of biomarkers, and the practical aspects of implementation of target concentration intervention. For selected topics with sufficient evidence, such as the application of limited sampling strategies for MPA area under the curve, graded recommendations on target ranges are presented. To provide a comprehensive review, this report also includes updates on the status of potential biomarkers including those which may be promising but with a low level of evidence. In view of the fact that there are very few new immunosuppressive drugs under development for the transplant field, it is likely that MPA will continue to be prescribed on a large scale in the upcoming years. Discontinuation of therapy due to adverse effects is relatively common, increasing the risk for late rejections, which may contribute to graft loss. Therefore, the continued search for innovative methods to better personalize MPA dosage is warranted.

Disease Progression in Patients with Autosomal Dominant Retinitis Pigmentosa due to a Mutation in Inosine Monophosphate Dehydrogenase 1 (IMPDH1)

Purpose

Mutations in the inosine monophosphate dehydrogenase 1 (IMPDH1) gene are a common cause of inherited retinal degeneration (IRD). Due to species- and tissue-dependent expression of IMPDH1, there are no appropriate models of human IMPDH1 disease. Therefore, a limited understanding remains of disease expression and rates of progression for IMPDH1-related IRD.

Methods

We evaluated semiautomated kinetic and chromatic static perimetry, spectral-domain optical coherence tomography (SD-OCT), and ultra-wide field fundus images with autofluorescence in a cohort of 12 patients (ages 11–58 at first visit). Ten patients had longitudinal data for which rates of progression were estimated.

Results

Visual acuities were relatively stable over time and the photoreceptors within the central retina remained intact. Perifoveal photoreceptor loss measured over a period of years coincided with visual fields, which were constricted and progressed over time in all patients. Rod sensitivity showed a similar pattern of defect to that of the kinetic perimetry and the autofluorescence ultra-wide field imaging. Full-field electroretinograms were severely reduced and the dark-adapted rod and mixed responses were extinguished at earlier visits than the light-adapted cone responses.

Conclusions

There was variability in disease severity at the first visit, but results show that the peripheral retina is more susceptible to the deleterious consequences of an IMPDH1 mutation. Given the pattern of degeneration and the alternatively spliced isoforms of IMPDH1, potential interventions may consider targeting the periphery early in disease, modulating transcript expression, and/or preserving central vision at late stages of the disease.

Translational Relevance

These results inform clinical prognosis and offer evidence strategies toward therapeutic intervention.

Effect of Protein Binding on Exposure of Unbound and Total Mycophenolic Acid: A Population Pharmacokinetic Analysis in Chinese Adult Kidney Transplant Recipients

OBJECTIVES

The population pharmacokinetic (popPK) characteristics of total mycophenolic acid (tMPA) have been investigated in various ethnic populations. However, investigations of popPK of unbound MPA (uMPA) are few. Thus, a popPK analysis was performed to: (1) characterize the PK of uMPA and tMPA and its 7-O-mycophenolic acid glucuronide (MPAG) metabolite in kidney transplant patients cotreated with cyclosporine (CsA), and (2) identify the clinically significant covariates that explain variability in the dose-exposure relationship.

METHODS

A total of 740 uMPA, 741 tMPA, and 734 total MPAG (tMPAG) concentration-time data from 58 Chinese kidney transplant patients receiving MPA in combination with CsA were analyzed using NONMEM® software with the stochastic approximation expectation maximization (SAEM) followed by the important sampling (IMP) method. The influence of covariates was tested using a stepwise procedure.

RESULTS

The PK of uMPA and unbound MPAG (uMPAG) were characterized by a two- and one-compartment model with first-order elimination, respectively. A linear protein binding model was used to link uMPA and tMPA. Apparent clearance (CL/F) and central volume of distribution (VC/F) of uMPA (CLuMPA/F and VCuMPA/F, respectively) and protein binding rate constant (k B) were estimated to be 851 L/h [relative standard error (RSE), 7.1%], 718 L (18.5%) and 53.4/h (2.3%), respectively. For uMPAG, the population values (RSE) of CL/F (CLuMPAG) and VC/F (VCuMPAG/F) were 5.71 L/h (4.4%) and 29.9 L (7.7%), respectively. Between-subject variability (BSVs) on CLuMPA/F, VCuMPA/F, CLuMPAG/F, and VCuMPAG/F were 51.0, 80.0, 31.8 and 48.4%, respectively, whereas residual unexplained variability (RUVs) for uMPA, tMPA, and uMPAG were 47.0, 45.9, and 22.0%, respectively. Significant relationships were found between k B and serum albumin (ALB) and between CLuMPAG/F and glomerular filtration rate (GFR). Additionally, model-based simulation showed that changes in ALB concentrations substantially affected tMPA but not uMPA exposure.

CONCLUSIONS

The established model adequately described the popPK characteristics of the uMPA, tMPA, and MPAG. The estimated CLuMPA/F and unbound fraction of MPA (FUMPA) in Chinese kidney transplant recipients cotreated with CsA were comparable to those published previously in Caucasians. We recommend monitoring uMPA instead of tMPA to optimize mycophenolate mofetil (MMF) dosing for patients with lower ALB levels.

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

BACKGROUND AND OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSION

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

Systematic external evaluation of published population pharmacokinetic models of mycophenolate mofetil in adult kidney transplant recipients co-administered with tacrolimus

AIMS

Various mycophenolate mofetil (MMF) population pharmacokinetic (popPK) models have been developed to describe its PK characteristics and facilitate its optimal dosing in adult kidney transplant recipients co-administered with tacrolimus. However, the external predictive performance has been unclear. Thus, this study aimed to comprehensively evaluate the external predictability of published MMF popPK models in such populations and investigate the potential influencing factors.

METHODS

The external predictability of qualified popPK models was evaluated using an independent dataset. The evaluation included prediction- and simulation-based diagnostics, and Bayesian forecasting. In addition, factors influencing model predictability, especially the impact of structural models, were investigated.

RESULTS

Fifty full PK profiles from 45 patients were included in the evaluation dataset and 11 published popPK models were identified and evaluated. In prediction-based diagnostics, the prediction error within ±30% was less than 50% in most published models. The prediction- and variability-corrected visual predictive check and posterior predictive check showed large discrepancies between the observations and simulations in most models. Moreover, the normalized prediction distribution errors of all models did not follow a normal distribution. Bayesian forecasting demonstrated an improvement in the model predictability. Furthermore, the predictive performance of two-compartment (2CMT) models incorporating the enterohepatic circulation (EHC) process was not superior to that of conventional 2CMT models.

CONCLUSIONS

The published models showed large variability and unsatisfactory predictive performance, which indicated that therapeutic drug monitoring was necessary for MMF clinical application. Further studies incorporating potential covariates need to be conducted to investigate the key factors influencing model predictability of MMF.

A Gallbladder-Based Enterohepatic Circulation Model for Pharmacokinetic Studies

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.