Prediction of Free from Total Mycophenolic Acid Concentrations in Stable Renal Transplant Patients: A Population-Based Approach

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.

Population pharmacokinetics and limited sampling strategy for therapeutic drug monitoring of mycophenolate mofetil in Japanese patients with lupus nephritis

BACKGROUND

Mycophenolate mofetil (MMF), a prodrug of the immunosuppressive agent mycophenolic acid (MPA), is difficult to administer because of the pharmacokinetic complexity of MPA. Although dosage adjustment according to the 12-h area under the concentration-time curve (AUC0-12) is thought to be desirable, multiple blood samplings for AUC calculation may pose a clinical challenge. A limited sampling strategy (LSS) would provide a solution; however, little is known about MPA pharmacokinetics in lupus nephritis patients, especially in those with Asian backgrounds, or few, if any, LSSs are reported for them.

METHODS

Thirty-four adult Japanese patients receiving MMF for lupus nephritis were examined retrospectively. MPA pharmacokinetics were investigated, and a PPK model was developed using Phoenix® NLME™ software. Single and double blood sampling strategies from Bayesian estimation using the PPK model and from multiple linear regression were compared. Tolerability was also evaluated.

RESULTS

In the pharmacokinetic analysis, renal function and serum albumin had significant effects on dose-normalized AUC0-12; and serum albumin, concomitant proton pump inhibitor (PPI) and iron/magnesium oxide did on dose-normalized maximum concentration. As a PPK model, a two-compartment model was developed with a transit absorption model and first-order elimination, in which creatinine clearance and serum albumin were covariates for MPA clearance. The double sampling strategy at 1 and 4 h by multiple linear regression showed the best agreement with the observed AUC0-12 (r2 = 0.885). Of the single sampling strategies, the one at 6 h by Bayesian estimation performed best (r2 = 0.769). The tolerability evaluation showed that correlations were suggested for gastrointestinal involvement.

CONCLUSIONS

The present study developed the first PPK model of MPA for Japanese lupus nephritis patients. As for LSSs, a double sampling strategy at 1 and 4 h by multiple linear regression would work best; when only a single blood sampling is allowed, a strategy at 6 h by Bayesian estimation using the PPK model developed in this study would be best. The LSSs good enough for clinical use may facilitate safer, more effective, and individualized therapy.

A model based on machine learning for the prediction of cyclosporin A trough concentration in Chinese allo-HSCT patients

BACKGROUND

Cyclosporin A is a calcineurin inhibitor which has a narrow therapeutic window and high interindividual variability. Various population pharmacokinetic models have been reported; however, professional software and technical personnel were needed and the variables of the models were limited. Therefore, the aim of this study was to establish a model based on machine learning to predict CsA trough concentrations in Chinese allo-HSCT patients.

METHODS

A total of 7874 cases of CsA therapeutic drug monitoring data from 2069 allo-HSCT patients were retrospectively included. Sequential forward selection was used to select variable subsets, and eight different algorithms were applied to establish the prediction model.

RESULTS

XGBoost exhibited the highest prediction ability. Except for the variables that were identified by previous studies, some rarely reported variables were found, such as norethindrone, WBC, PAB, and hCRP. The prediction accuracy within ±30% of the actual trough concentration was above 0.80, and the predictive ability of the models was demonstrated to be effective in external validation.

CONCLUSION

In this study, models based on machine learning technology were established to predict CsA levels 3-4 days in advance during the early inpatient phase after HSCT. A new perspective for CsA clinical application is provided.

Utility of salivary mycophenolic acid concentration monitoring: Modeling and Monte Carlo validation approach

The results of the previous studies demonstrated an association between mycophenolic acid (MPA) exposure, serum albumin level (ALB), and adverse effects in kidney transplant patients. The aim was the identification of mathematical correlation and association between both, total and unbound MPA concentration in relation to ALB, body mass (BM), age and estimated glomerular filtration rate (eGFR) in stable kidney transplant recipients. Furthermore, investigation was conducted with the aim to clarify the role of salivary concentration (CSAL ) of MPA in adverse effect profile. In order to analyze the association between total and salivary concentration of MPA in relation to ALB, BM, age and eGFR, a least squares method for determining the correlation between these parameters was performed. In addition, derived mathematical model based on experimental data can also be performed and simulated through the Monte Carlo (MC) approach. Adverse effects were grouped according to the nature of symptoms and scored by a previously published validated system. Numerically calculated values of CSAL from the models [CSAL  = f(ALB, BM, age, eGFR, CP ) = a00  + a10 *(ALB, BM, age, eGFR) + a01 *CP ] were then compared with those from validation set of patients, where the best fitting model was for ALB [CSAL  = 54.96-1.64*ALB +13.4*CP ]. Adverse effects estimation showed the difference in esthetic score, positively correlated with CSAL in the lower ALB group (145.41 ± 219.02 vs. 354.08 ± 262.19; with statistical significance p = .014) and almost significant for gastrointestinal score (167.69 ± 174.79 vs. 347.55 ± 320.95; p = .247). The study showed that CSAL MPA may contribute to management of adverse effects, but these findings require confirmation of clinical utility.

Significant Correlations between p-Cresol Sulfate and Mycophenolic Acid Plasma Concentrations in Adult Kidney Transplant Recipients

BACKGROUND AND OBJECTIVES

Mycophenolic acid (MPA) is a commonly prescribed life-long immunosuppressant for kidney transplant recipients. The frequently observed large variations in MPA plasma exposure may lead to severe adverse outcomes; therefore, characterizations of contributing factors can potentially improve the precision dosing of MPA. Our group recently reported the potent inhibitory effects of p-cresol (a protein-bound uremic toxin that can be accumulated in kidney transplant patients) on the hepatic metabolism of MPA in human in vitro models. Based on these data, the hypothesis for this clinical investigation was that a direct correlation between p-cresol and MPA plasma exposure should be evident in adult kidney transplant recipients.

METHODS

Using a prospective and observational approach, adult kidney transplant recipients within the first year after transplant on oral mycophenolate mofetil (with tacrolimus ± prednisone) were screened for recruitment. The exclusion criteria were cold ischemia time > 30 h, malignancy, pregnancy, severe renal dysfunction (i.e., estimated glomerular filtration rate, eGFR, < 10 mL/min/1.73 m²), active graft rejection, or MPA intolerance. Patients' demographic and biochemistry data were collected. Total and free plasma concentrations of MPA, MPA glucuronide (MPAG), and total p-cresol sulfate (the predominant, quantifiable form of p-cresol in the plasma) were quantified using validated assays. Correlational and categorical analyses were performed using GraphPad Prism.

RESULTS

Forty patients (11 females) were included: donor type (living/deceased: 20/20), induction regimen (basiliximab/thymoglobulin/basiliximab followed by thymoglobulin: 35/3/2), post-transplant time (74 ± 60 days, mean ± standard deviation), age (53.7 ± 12.4 years), bodyweight (79.8 ± 18.5 kg), eGFR (51.9 ± 18.0 mL/min/1.73 m²), serum albumin (3.6 ± 0.5 g/dL), prednisone dose (18.5 ± 13.2 mg, n = 33), and tacrolimus trough concentration (9.4 ± 2.4 µg/L). Based on Spearman analysis, significant control correlations supporting the validity of our dataset were observed between total MPA trough concentration (C₀) and total MPAG C₀ (correlation coefficient [R] = 0.39), ratio of total MPAG C₀-to-total MPA C₀ and post-transplant time (R = - 0.56), total MPAG C₀ and eGFR (R = - 0.35), and p-cresol sulfate concentration and eGFR (R = - 0.70). Our primary analysis indicated the novel observation that total MPA C₀ (R = 0.39), daily dose-normalized total MPA C₀ (R = 0.32), and bodyweight-normalized total MPA C₀ (R = 0.32) were significantly correlated with plasma p-cresol sulfate concentrations. Consistently, patients categorized with elevated p-cresol sulfate concentrations (i.e., ≥ median of 3.2 µg/mL) also exhibited increased total MPA C₀ (by 57 % vs those below median), daily dose-normalized total MPA C₀ (by 89 %), and bodyweight-normalized total MPA C₀ (by 62 %). Our secondary analyses with MPA metabolites, unbound concentrations, free fractions, and MPA metabolite ratios supported additional potential interacting mechanisms.

CONCLUSION

We have identified a novel, positive association between p-cresol sulfate exposure and total MPA C₀ in adult kidney transplant recipients, which is supported by published mechanistic in vitro data. Our findings confirm a potential role of p-cresol as a significant clinical variable affecting the pharmacokinetics of MPA. These data also provide the justifications for conducting subsequent full-scale pharmacokinetic-pharmacodynamic studies to further characterize the cause-effect relationships of this interaction, which could also rule out potential confounding variables not adequately controlled in this correlational study.

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.

Estimation of Mycophenolic Acid Exposure in Heart Transplant Recipients by Population Pharmacokinetic and Limited Sampling Strategies

Purpose: The aim of this study is i) to establish a strategy to estimate the area under the curve of the dosing interval (AUC0-12h) of mycophenolic acid (MPA) in the heart transplant recipients and ii) to find the covariates that significantly affect the pharmacokinetics of MPA exposure. Methods: This single-center, prospective, open-label, observational study was conducted in 91 adult heart transplant recipients orally taking mycophenolate mofetil dispersible tablets. Samples collected intensively and sparsely were analyzed by the enzyme-multiplied immunoassay technique, and all the data were used in PPK modeling. Potential covariates were tested stepwise. The goodness-of-fit plots, the normalized prediction distribution error, and prediction-corrected visual predictive check were used for model evaluation. Optimal sampling times by ED-optimal strategy and multilinear regression (MLR) were analyzed based on the simulated data by the final PPK model. Moreover, using intensive data from 14 patients, the accuracy of AUC0-12h estimation was evaluated by Passing-Bablok regression analysis and Bland-Alman plots for both the PPK model and MLR equation. Results: A two-compartment model with first-order absorption and elimination with a lag time was chosen as the structure model. Co-medication of proton pump inhibitors (PPIs), estimated glomerular filtration rate (eGFR), and albumin (ALB) were found to significantly affect bioavailability (F), clearance of central compartment (CL/F), and the distribution volume of the central compartment (V2/F), respectively. Co-medication of PPIs decreased F by 27.6%. When eGFR decreased by 30 ml/min/1.73 m2, CL/F decreased by 23.7%. However, the impact of ALB on V2/F was limited to MPA exposure. The final model showed an adequate fitness of the data. The optimal sampling design was pre-dose and 1 and 4 h post-dose for pharmacokinetic estimation. The best-fit linear equation was finally established as follows: AUC0-12h = 3.539 × C0 + 0.288 × C0.5 + 1.349 × C1 + 6.773 × C4.5. Conclusion: A PPK model was established with three covariates in heart transplant patients. Co-medication of PPIs and eGFR had a remarkable impact on AUC0-12h of MPA. A linear equation was also concluded with four time points as an alternative way to estimate AUC0-12h for MPA.

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.

The use of mycophenolate mofetil area under the curve

PURPOSE OF REVIEW

Although mycophenolate mofetil (MMF) has been used successfully to treat a myriad of autoimmune diseases, its complex pharmacokinetics make it difficult to determine the true drug exposure for an individual patient. This review summarizes the body of literature focused on the gold standard measurement of the area under the curve (AUC) of mycophenolic acid (MPA), the active metabolite of MMF.

RECENT FINDINGS

Fixed dosing of MMF leads to highly variable drug exposure. Retrospective series have reported improved clinical outcomes when a minimum AUC value from 0 to 12 h (AUC0-12h) ≥30 mg h/l is achieved. MPA levels are affected by various drug interactions, hypoalbuminemia, and renal insufficiency and the measurement of free rather than total MPA levels is prudent in some situations. A limited number of studies employing prospective dose adjustment of MMF based on AUC0-12h measurements have yielded mixed results.

SUMMARY

Given the wide range of MPA AUC encountered in autoimmune diseases, dose adjustments of MMF based on AUC rather than fixed dosing of MMF should be considered in both clinical practice and clinical trials. Limited sampling strategies have been proposed to improve clinical feasibility of measurements, but a standard is yet to be defined.

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.

Population Pharmacokinetic Model of Plasma and Cellular Mycophenolic Acid in Kidney Transplant Patients from the CIMTRE Study

Background and Objective

Mycophenolate mofetil is widely used in kidney transplant recipients. Mycophenolate mofetil is hydrolysed by blood esterases to mycophenolic acid (MPA), the active drug. Although MPA therapeutic drug monitoring has been recommended to optimise the treatment efficacy by the area under the plasma concentration vs time curve, little is known regarding MPA concentrations in peripheral blood mononuclear cells, where MPA inhibits inosine monophosphate dehydrogenase. This study aimed to build a pharmacokinetic model using a population approach to describe MPA total and unbound concentrations in plasma and into peripheral blood mononuclear cells in 78 adult kidney transplant recipients receiving mycophenolate mofetil therapy combined with tacrolimus and prednisone.

Methods

Total and unbound plasma concentrations and peripheral blood mononuclear cell concentrations were assayed. A three-compartment model, two for plasma MPA and one for peripheral blood mononuclear cell MPA, with a zero-order absorption and a first-order elimination was used to describe the data.

Results

Mycophenolic acid average concentrations in peripheral blood mononuclear cells were well above half-maximal effective concentration for inosine monophosphate dehydrogenase and no relationship was found with the occurrence of graft rejection. Three covariates affected unbound and intracellular MPA pharmacokinetics: creatinine clearance, which has an effect on unbound MPA clearance, human serum albumin, which influences fraction unbound MPA and the ABCB1 3435 C>T (rs1045642) genetic polymorphism, which has an effect on MPA efflux transport from peripheral blood mononuclear cells.

Conclusion

This population pharmacokinetic model demonstrated the intracellular accumulation of MPA, the efflux of MPA out of the cells being dependent on P-glycoprotein transporters. Nevertheless, further studies are warranted to investigate the relevance of MPA concentrations in peripheral blood mononuclear cells to dosing regimen optimisation.

Electronic supplementary material

The online version of this article (10.1007/s40268-020-00319-y) contains supplementary material, which is available to authorized users.

Population pharmacokinetics of mycophenolic acid in paediatric patients

Mycophenolic acid (MPA) is widely used in paediatric kidney transplant patients and sometimes prescribed for additional indications. Population pharmacokinetic or pharmacodynamic modelling has been frequently used to characterize the fixed, random and covariate effects of MPA in adult patients. However, MPA population pharmacokinetic data in the paediatric population have not been systematically summarized. The objective of this narrative review was to provide an up-to-date critique of currently available paediatric MPA population pharmacokinetic models, with emphases on modelling techniques, pharmacological findings and clinical relevance. PubMed and EMBASE were searched from inception of database to May 2020, where a total of 11 studies have been identified representing kidney transplant (n = 4), liver transplant (n = 1), haematopoietic stem cell transplant (n = 1), idiopathic nephrotic syndrome (n = 2), systemic lupus erythematosus (n = 2), and a combined population consisted of kidney, liver and haematopoietic stem cell transplant patients (n = 1). Critical analyses were provided in the context of MPA absorption, distribution, metabolism, excretion and bioavailability in this paediatric database. Comparisons to adult patients were also provided. With respect to clinical utility, Bayesian estimation models (n = 6) with acceptable accuracy and precision for MPA exposure determination have also been identified and systematically evaluated. Overall, our analyses have identified unique features of MPA clinical pharmacology in the paediatric population, while recognizing several gaps that still warrant further investigations. This review can be used by pharmacologists and clinicians for improving MPA pharmacokinetic-pharmacodynamic modelling and patient care.

Mycophenolic acid interferes the transcriptional regulation and protein trafficking of maturation surface markers in dendritic cells

BACKGROUND

The ability of dendritic cells (DCs) to regulate adaptive immunity makes them interesting cells to be used as therapeutic targets modulating alloimmune responses. Mycophenolic acid (MPA) is an immunosuppressor commonly used in transplantation, and its effect on DCs has not been fully investigated.

METHODS

Monocyte-derived DCs were obtained from healthy volunteers and cultured for 7 days. Cells were treated with MPA on day 2 and matured by lipopolysaccharide (LPS) stimulation. Functionality of mature DC (mDCs) was evaluated by allogeneic mixed lymphocytes reaction. Surface expression of maturation markers (CD40, CD83, CD86, and ICAM-1) was analyzed in both immature DCs (iDCs) and mDCs by flow cytometry. To assess transcriptional regulation and protein subcellular location, RT-PCR and confocal microscopy were used, respectively.

RESULTS

MPA decreased surface expression of all maturation markers in mDCs and significantly abrogated DCs-induced allogeneic T-cell proliferation after MPA pre-treatment. In iDCs, the reduced surface protein expression after MPA paralleled with mRNA downregulation of their genes. In mDCs, the mRNA levels of ICAM-1, CD40 and CD83 were enhanced in MPA-treated mDCs with an increase in the expression of CD83 and ICAM-1 near the Golgi compared to non-treated mDCs. In contrast, mRNA levels of CD86 were diminished after MPA treatment.

CONCLUSIONS

The reduced surface markers expression in mDCs exerted by MPA produced a decline in their capacity to activate immune responses. Moreover, the inhibition of guanosine-derived nucleotide biosynthesis by MPA treatment leads to DC maturation interference by two mechanisms depending on the marker, transcriptional downregulation or disrupted intracellular protein trafficking.

Assessment of pharmacokinetic mycophenolic acid clearance models using Monte Carlo numerical analysis

Previously, we performed population pharmacokinetic analysis and indicated age, mycophenolate mofetil (MMF)/mycophenolic acid (MPA) daily dose, and presence of nifedipine in patient therapy as significant predictors of MPA apparent clearance (CL/F) variability. This study aimed to determine the reliability of previously published population pharmacokinetic models derived from similar studies. Furthermore, this study investigated correspondence between chosen population models from the literature.By means of the Monte Carlo simulation method, pharmacokinetic models from different studies are simulated and analysed in the range of standard deviations of measured system parameters as well as the range of observed model parameters taken from the comparison studies.The 1000 numerical simulations were performed for every analysed model in order to calculate the most possible MPA CL/F values according to the expected values from the performed experiment. Fitting our results with other models showed how the presence of nifedipine makes difference in MPA CL/F values.By testing the data from selected studies into our model, a similar range of expected CL/F values was obtained, which may confirm the validity of our model. The results of our population pharmacokinetic study are partially applicable in models by other researchers.

Early prognostic performance of miR155-5p monitoring for the risk of rejection: Logistic regression with a population pharmacokinetic approach in adult kidney transplant patients

Previous results from our group and others have shown that urinary pellet expression of miR155-5p and urinary CXCL-10 production could play a key role in the prognosis and diagnosis of acute rejection (AR) in kidney transplantation patients. Here, a logistic regression model was developed using NONMEM to quantify the relationships of miR155-5p urinary expression, CXCL-10 urinary concentration and tacrolimus and mycophenolic acid (MPA) exposure with the probability of AR in adult kidney transplant patients during the early post-transplant period. Owing to the contribution of therapeutic drug monitoring to achieving target exposure, neither tacrolimus nor MPA cumulative exposure was identified as a predictor of AR in the studied population. Even though CXCL-10 urinary concentration showed a trend, its effect on AR was not significant. In contrast, urinary miR155-5p expression was prognostic of clinical outcome. Monitoring miR155-5p urinary pellet expression together with immunosuppressive drug exposure could be very useful during routine clinical practice to identify patients with a potential high risk of rejection at the early stages of the post-transplant period. This early risk assessment would allow for the optimization of treatment and improved prevention of AR.

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.

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.