Population pharmacokinetic-pharmacodynamic modelling of mycophenolic acid in paediatric renal transplant recipients in the early post-transplant period

Body Surface Area-Based Dosing of Mycophenolate Mofetil in Pediatric Hematopoietic Stem Cell Transplant Recipients: A Prospective Population Pharmacokinetic Study

Mycophenolate mofetil (MMF) is commonly used for acute graft-versus-host disease (aGVHD) after allogeneic hematopoietic stem cell transplantation (HSCT). However, limited population pharmacokinetic (PPK) data are available for pediatric HSCT patients. This study aimed to develop a PPK model and recommend optimal oral MMF dosage in pediatric HSCT patients. This prospective study involved pediatric HSCT patients at a tertiary academic institution. Patients received oral MMF 15-20 mg/kg twice daily for aGVHD prophylaxis and treatment. The PPK analysis was conducted using a nonlinear mixed-effects modeling method. Simulation was performed considering different body surface areas (BSAs) (0.5 m2, 1.0 m2, 1.5 m2) and dosing (400 mg/m2, 600 mg/m2, 900 mg/m2 twice daily). Based on the simulation, an optimal dosage of oral MMF was suggested. A total of 20 patients and 80 samples were included in the PPK model development. A one-compartment model with first-order absorption adequately described the pharmacokinetics of mycophenolic acid (MPA). BSA was a statistically significant covariate on Vd/F. Simulation suggested the optimal dosage of oral MMF as 900 mg/m2 twice daily, respectively. A reliable PPK model was developed with good predictive performance. This model-informed optimal MMF dosage in pediatric HSCT patients can provide valuable dosing guidance in real-world clinical practice.

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.

Individualization of Mycophenolic Acid Therapy through Pharmacogenetic, Pharmacokinetic and Pharmacodynamic Testing

Mycophenolic acid (MPA) is a widely used immunosuppressive agent and exerts its effect by inhibiting inosine 5′-monophosphate dehydrogenase (IMPDH), the main regulating enzyme of purine metabolism. However, significant unexplained differences in the efficacy and tolerability of MPA therapy pose a clinical challenge. Therefore, broad pharmacogenetic, pharmacokinetic, and pharmacodynamic approaches are needed to individualize MPA therapy. In this prospective cohort study including 277 renal transplant recipients, IMPDH2 rs11706052 SNP status was assessed by genetic sequencing, and plasma MPA trough levels were determined by HPLC and IMPDH enzyme activity in peripheral blood mononuclear cells (PBMCs) by liquid chromatography–mass spectrometry. Among the 277 patients, 84 were identified with episodes of biopsy-proven rejection (BPR). No association was found between rs11706052 SNP status and graft rejection (OR 1.808, and 95% CI, 0.939 to 3.479; p = 0.076). Furthermore, there was no association between MPA plasma levels and BPR (p = 0.69). However, the patients with graft rejection had a significantly higher predose IMPDH activity in PBMCs compared to the controls without rejection at the time of biopsy (110.1 ± 50.2 vs. 95.2 ± 45.4 pmol/h; p = 0.001), and relative to the baseline IMPDH activity before transplantation (p = 0.042). Our results suggest that individualization of MPA therapy, particularly through pharmacodynamic monitoring of IMPDH activity in PBMCs, has the potential to improve the clinical outcomes of transplant patients.

Population pharmacokinetics of mycophenolate mofetil in pediatric patients early after liver transplantation

Objective: To investigate the factors influencing the pharmacokinetics of mycophenolate mofetil (MMF) in pediatric patients after liver transplantation, and to establish a population pharmacokinetics model, which can provide a reference for clinical dosage adjustment.

Methods: A prospective study in a single center was performed on pediatric patients who were administrated with mycophenolate mofetil dispersible tablets (MMFdt) for at least 4 days after liver transplantation continuously. Blood samples were collected in ethylene diamine tetraacetic acid anticoagulant tubes before dosing and 0.5, 1, 2, 4, 8, and 12 h after the morning intake of MMFdt. The concentrations of mycophenolic acid (MPA) in plasma were assayed with a validated reverse-phase high-performance liquid chromatography method. UGT1A8 518C > G, UGT1A9 -275T > A, UGT1A9 -2152C > T, UGT2B7 211G > T, SLC O 1B1 521T > C polymorphism were determined by Sanger sequencing. Nonlinear mixed effects modeling was used to establish the population pharmacokinetics (PPK) model. The predictability and stability of the model were internally evaluated by the goodness of fit plots, visual prediction check, normalized prediction errors, and bootstraps.

Results: A two-compartment model with first-order absorption and first-order elimination was established with 115 MPA concentrations from 20 pediatric patients. The final model were: CL/F (L/h) = 14.8×(WT/7.5)^0.75×(DOSE/11.16)^0.452×е^0.06, Ka (h⁻¹) = 2.02×(WT/7.5)^−0.25, Vc/F (L) = 6.01×(WT/7.5), Vp/F (L) = 269 (fixed), Q/F (L/h) = 15.4×(WT/7.5)^0.75×е^1.39. Where CL/F was the apparent clearance rate, Ka was the absorption rate constant, Vc/F was the apparent distribution volume of the central compartment, Vp/F was the apparent distribution volume of the peripheral compartment, Q/F was the atrioventricular clearance rate, WT was the body weight of the subject, and DOSE was the MMFdt administered dose. The model indicated there was large inter-individual variability in CL/F and Q/F after multiple dosing of MMFdt. Internal evaluation results showed that the final model had good stability and prediction performance.

Conclusion: A stable and predictive population pharmacokinetic model of MMFdt in pediatric patients after the early stage of liver transplantation was established. The pediatric patient’s weight and the dose of MMFdt can be a reference to adjust the MMFdt dose.

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.

Recent Advances in Therapeutic Drug Monitoring of Voriconazole, Mycophenolic Acid, and Vancomycin: A Literature Review of Pediatric Studies

The review includes studies dated 2011-2021 presenting the newest information on voriconazole (VCZ), mycophenolic acid (MPA), and vancomycin (VAN) therapeutic drug monitoring (TDM) in children. The need of TDM in pediatric patients has been emphasized by providing the information on the differences in the drugs pharmacokinetics. TDM of VCZ should be mandatory for all pediatric patients with invasive fungal infections (IFIs). Wide inter- and intrapatient variability in VCZ pharmacokinetics cause achieving and maintaining therapeutic concentration during therapy challenging in this population. Demonstrated studies showed, in most cases, VCZ plasma concentrations to be subtherapeutic, despite the updated dosages recommendations. Only repeated TDM can predict drug exposure and individualizing dosing in antifungal therapy in children. In children treated with mycophenolate mofetil (MMF), similarly as in adult patients, the role of TDM for MMF active form, MPA, has not been well established and is undergoing continued debate. Studies on the MPA TDM have been carried out in children after renal transplantation, other organ transplantation such as heart, liver, or intestine, in children after hematopoietic stem cell transplantation or cord blood transplantation, and in children with lupus, nephrotic syndrome, Henoch-Schönlein purpura, and other autoimmune diseases. MPA TDM is based on the area under the concentration-time curve; however, the proposed values differ according to the treatment indication, and other approaches such as pharmacodynamic and pharmacogenetic biomarkers have been proposed. VAN is a bactericidal agent that requires TDM to prevent an acute kidney disease. The particular group of patients is the pediatric one. For this group, the general recommendations of the dosing may not be valid due to the change of the elimination rate and volume of distribution between the subjects. The other factor is the variability among patients that concerns the free fraction of the drug. It may be caused by both the patients' population and sample preconditioning. Although VCZ, MMF, and VAN have been applied in pediatric patients for many years, there are still few issues to be solve regarding TDM of these drugs to ensure safe and effective treatment. Except for pharmacokinetic approach, pharmacodynamics and pharmacogenetics have been more often proposed for TDM.

Pharmacokinetic Model Analysis of Supralingual, Oral and Intravenous Deliveries of Mycophenolic Acid

Mycophenolic acid (MPA) is commonly used for organ rejection prophylaxis via oral administration in the clinic. Recent studies have shown that MPA also has anticancer activities. To explore new therapeutic options for oral precancerous/cancerous lesions, MPA was designed to release topically on the dorsal tongue surface via a mucoadhesive patch. The objective of this study was to establish the pharmacokinetic (PK) and tongue tissue distribution of mucoadhesive MPA patch formulation after supralingual administration in rats and also compare the PK differences between oral, intravenous, and supralingual administration of MPA. Blood samples were collected from Sprague Dawley rats before and after a single intravenous bolus injection, a single oral dose, or a mucoadhesive patch administration on the dorsal tongue surface for 4 h, all with a dose of 0.5 mg/kg of MPA. Plots of MPA plasma concentration versus time were obtained. As multiple peaks were found in all three curves, the enterohepatic recycling (EHR) model in the Phoenix software was adapted to describe their PK parameters with an individual PK analysis method. The mean half-lives of intravenous and oral administrations were 10.5 h and 7.4 h, respectively. The estimated bioavailability after oral and supralingual administration was 72.4% and 7.6%, respectively. There was a 0.5 h lag-time presented after supralingual administration. The results suggest that the systemic plasma MPA concentrations were much lower in rats receiving supralingual administration compared to those receiving doses from the other two routes, and the amount of MPA accumulated in the tongue after patch application showed a sustained drug release pattern. Studies on the dynamic of drug retention in the tongue after supralingual administration showed that ~3.8% of the dose was accumulated inside of tongue right after the patch removal, ~0.11% of the dose remained after 20 h, and ~20.6% of MPA was not released from the patches 4 h after application. The data demonstrate that supralingual application of an MPA patch can deliver a high amount of drug at the site of administration with little systemic circulation exposure, hence lowering the potential gastrointestinal side effects associated with oral administration. Thus, supralingual administration is a potential alternative route for treating oral lesions.

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

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.

PK/PD Study of Mycophenolate Mofetil in Children With Systemic Lupus Erythematosus to Inform Model-Based Precision Dosing

Objectives: To evaluate the mycophenolic acid [MPA, the active form of mycophenolate mofetil (MMF)] pharmacokinetic parameters in relation to clinical response to identify target exposure ranges in pediatric patients with systemic lupus erythematosus (SLE).

Methods: This was a retrospective study using pharmacokinetic data collected in 67 pediatric patients aged 4–18 years with SLE. Target MPA exposures for effective inhibition of SLE activity (as measured by SLE disease Activity Index (SLEDAI), active SLE was defined as a SLEDAI score of ≥6, and a controlled disease was defined as a SLEDAI score of ≤4) were assessed by receiver operating characteristic (ROC) curve and logistic regression. Exposure-response models were developed to quantitatively describe the relationship between SLEDAI score and AUC_0–12 or C_trough, respectively.

Results: The MPA AUC_0-12 in patients with active SLE was significantly lower than that in patients with inactive SLE. ROC analysis revealed that an AUC_0–12 threshold of 39 μg h/ml or a C_trough of 1.01 μg/ml was associated with the lowest risk of active SLE. Logistic regression analysis revealed that an AUC_0–12 of less than 34 μg h/ml or a C_trough of less than 1.2 μg/ml probably is associated with active SLE. The results of the exposure-response modeling also indicated that an AUC_0-12 less than 32 μg h/ml or a C_trough less than 1.1 μg/ml was associated with suboptimal clinical outcome. An AUC_0-12 above 50 μg h/ml or a C_trough above 1.7 ug/ml was associated with disease control.

Conclusion: Both AUC_0–12 and C_trough of MPA are predictive of the likelihood of active SLE in pediatric patients receiving MMF. An individualized dosing regimen of MMF, with a target AUC_0–12 or C_trough, should be considered for SLE patients.

Pharmacodynamic assessment of mycophenolic acid in resting and activated target cell population during the first year after renal transplantation

Aims

To explore the pharmacodynamics of mycophenolic acid (MPA) through inosine monophosphate dehydrogenase (IMPDH) capacity measurement and purine levels in peripheral blood mononuclear cells (PBMC) longitudinally during the first year after renal transplantation (TX).

Methods

PBMC were isolated from renal recipients 0–4 days prior to and 6–9 days, 5–7 weeks and 1 year after TX (before and 1.5 hours after dose). IMPDH capacity and purine (guanine and adenine) levels were measured in stimulated and nonstimulated PBMC.

Results

Twenty‐nine patients completed the follow‐up period, of whom 24 received MPA. In stimulated PBMC, the IMPDH capacity (pmol 10⁻⁶ cells min⁻¹) was median (interquartile range) 127 (95.8–147) before TX and thereafter 44.9 (19.2–93.2) predose and 12.1 (4.64–23.6) 1.5 hours postdose across study days after TX. The corresponding IMPDH capacity in nonstimulated PBMC was 5.71 (3.79–6.93), 3.35 (2.31–5.62) and 2.71 (1.38–4.08), respectively. Predose IMPDH capacity in nonstimulated PBMC increased with time, reaching pre‐TX values at 1 year. In stimulated PBMC, both purines were reduced before (median 39% reduction across days after TX) and after (69% reduction) dose compared to before TX. No alteration in the purine levels was observed in nonstimulated PBMC. Patients needing dose reductions during the first year had lower pre‐dose IMPDH capacity in nonstimulated PBMC (1.87 vs 3.00 pmol 10⁻⁶ cells min⁻¹, P = .049) at 6–9 days.

Conclusion

The inhibitory effect of MPA was stronger in stimulated PBMC. Nonstimulated PBMC became less sensitive to MPA during the first year after TX. Early IMPDH capacity appeared to be predictive of dose reductions.

Population Pharmacokinetics of Mycophenolic Acid: An Update

The most recent comprehensive reviews on the population pharmacokinetics of mycophenolic acid (MPA) were published in 2014. Since then, several population pharmacokinetic studies on MPA have been published. The majority of literature is still focused on the kidney transplant population, although studies have also been conducted in liver and lung transplantation, autoimmune diseases, and hematopoietic stem cell transplant. While the majority of the model building is still based on parametric non-linear mixed-effects modeling, recent studies suggest the suitability of other methodologies. Additionally, instead of just focusing on pharmacokinetic modeling, a trend toward describing the relationships between pharmacokinetic and pharmacodynamic parameters is observed. Given the importance of enterohepatic recirculation (EHR) in the pharmacokinetics of MPA, more authors have attempted to characterize this process in their models. Overall, the recent models have become more sophisticated and incorporate EHR, pharmacodynamic relationships, and metabolites while maintaining many of the population values and covariates identified previously. However, the number of MPA population pharmacokinetic models describing the enteric-coated formulation of MPA (EC-MPA) is still limited. Given the increasing use of EC-MPA, more studies are needed to fill this literature gap. In addition, few studies are yet available characterizing free MPA concentration or MPA metabolites. Given the extensive protein binding, low to intermediate extraction, and intrinsic clearance characteristics of MPA in humans, including these variables would improve the population structural models.

Pharmacokinetics and dynamics of mycophenolate mofetil after single-dose oral administration in juvenile dachshunds

Mycophenolate mofetil (MMF) is recommended as an alternative/complementary immunosuppressant. Pharmacokinetic and dynamic effects of MMF are unknown in young-aged dogs. We investigated the pharmacokinetics and pharmacodynamics of single oral dose MMF metabolite, mycophenolic acid (MPA), in healthy juvenile dogs purpose-bred for the tripeptidyl peptidase 1 gene (TPP1) mutation. The dogs were heterozygous for the mutation (nonaffected carriers). Six dogs received 13 mg/kg oral MMF and two placebo. Pharmacokinetic parameters derived from plasma MPA were evaluated. Whole-blood mitogen-stimulated T-cell proliferation was determined using a flow cytometric assay. Plasma MPA C_max (mean ± SD, 9.33 ± 7.04 μg/ml) occurred at <1 hr. The AUC_0-∞ (mean ± SD, 12.84±6.62 hr*μg/ml), MRT_inf (mean ± SD, 11.09 ± 9.63 min), T1/2 (harmonic mean ± PseudoSD 5.50 ± 3.80 min), and k/d (mean ± SD, 0.002 ± 0.001 1/min). Significant differences could not be detected between % inhibition of proliferating CD5+ T lymphocytes at any time point (p = .380). No relationship was observed between MPA concentration and % inhibition of proliferating CD5+ T lymphocytes (R = .148, p = .324). Pharmacodynamics do not support the use of MMF in juvenile dogs at the administered dose based on existing therapeutic targets.